ABSTRACT Understanding the evolution of self-gravitating, isothermal, magnetized gas is crucial for star formation, as these physical processes have been postulated to set the initial mass function (IMF). We present a suite of isothermal magnetohydrodynamic (MHD) simulations using the gizmo code that follow the formation of individual stars in giant molecular clouds (GMCs), spanning a range of Mach numbers found in observed GMCs ($\mathcal {M} \sim 10\!-\!50$). As in past works, the mean and median stellar masses are sensitive to numerical resolution, because they are sensitive to low-mass stars that contribute a vanishing fraction of the overall stellar mass. The mass-weighted median stellar mass M50 becomes insensitive to resolution once turbulent fragmentation is well resolved. Without imposing Larson-like scaling laws, our simulations find $M_\mathrm{50} \,\, \buildrel\propto \over \sim \,\,M_\mathrm{0} \mathcal {M}^{-3} \alpha _\mathrm{turb}\, \mathrm{SFE}^{1/3}$ for GMC mass M0, sonic Mach number $\mathcal {M}$, virial parameter αturb, and star formation efficiency SFE = M⋆/M0. This fit agrees well with previous IMF results from the ramses, orion2, and sphng codes. Although M50 has no significant dependence on the magnetic field strength at the cloud scale, MHD is necessary to prevent a fragmentation cascade that results in non-convergent stellar masses. For initial conditions andmore »
Less wrong: a more realistic initial condition for simulations of turbulent molecular clouds
Simulations of isolated giant molecular clouds (GMCs) are an important tool for studying the dynamics of star formation, but their turbulent initial conditions (ICs) are uncertain. Most simulations have either initialized a velocity field with a prescribed power spectrum on a smooth density field (failing to model the full structure of turbulence) or ‘stirred’ turbulence with periodic boundary conditions (which may not model real GMC boundary conditions). We develop and test a new GMC simulation setup (called turbsphere) that combines advantages of both approaches: we continuously stir an isolated cloud to model the energy cascade from larger scales, and use a static potential to confine the gas. The resulting cloud and surrounding envelope achieve a quasi-equilibrium state with the desired hallmarks of supersonic ISM turbulence (e.g. density PDF and a ∼k−2 velocity power spectrum), whose bulk properties can be tuned as desired. We use the final stirred state as initial conditions for star formation simulations with self-gravity, both with and without continued driving and protostellar jet feedback, respectively. We then disentangle the respective effects of the turbulent cascade, simulation geometry, external driving, and gravity/MHD boundary conditions on the resulting star formation. Without external driving, the new setup obtains results more »
- Publication Date:
- NSF-PAR ID:
- 10361719
- Journal Name:
- Monthly Notices of the Royal Astronomical Society
- Volume:
- 510
- Issue:
- 4
- Page Range or eLocation-ID:
- p. 4767-4778
- ISSN:
- 0035-8711
- Publisher:
- Oxford University Press
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
ABSTRACT Most observed stars are part of a multiple star system, but the formation of such systems and the role of environment and various physical processes is still poorly understood. We present a suite of radiation-magnetohydrodynamic simulations of star-forming molecular clouds from the STARFORGE project that include stellar feedback with varied initial surface density, magnetic fields, level of turbulence, metallicity, interstellar radiation field, simulation geometry and turbulent driving. In our fiducial cloud, the raw simulation data reproduces the observed multiplicity fractions for Solar-type and higher mass stars, similar to previous works. However, after correcting for observational incompleteness the simulation underpredicts these values. The discrepancy is likely due to the lack of disc fragmentation, as the simulation only resolves multiples that form either through capture or core fragmentation. The raw mass distribution of companions is consistent with randomly drawing from the initial mass function for the companions of $\gt 1\, \mathrm{M}_{\rm \odot }$ stars. However, accounting for observational incompleteness produces a flatter distribution similar to observations. We show that stellar multiplicity changes as the cloud evolves and anticorrelates with stellar density. This relationship also explains most multiplicity variations between runs, i.e. variations in the initial conditions that increase stellar density (increasedmore »
-
null (Ed.)Abstract We present a large suite of MHD simulations of turbulent, star-forming giant molecular clouds (GMCs) with stellar feedback, extending previous work by simulating 10 different random realizations for each point in the parameter space of cloud mass and size. It is found that once the clouds disperse due to stellar feedback, both self-gravitating star clusters and unbound stars generally remain, which arise from the same underlying continuum of substructured stellar density, ie. the hierarchical cluster formation scenario. The fraction of stars that are born within gravitationally-bound star clusters is related to the overall cloud star formation efficiency set by stellar feedback, but has significant scatter due to stochastic variations in the small-scale details of the star-forming gas flow. We use our numerical results to calibrate a model for mapping the bulk properties (mass, size, and metallicity) of self-gravitating GMCs onto the star cluster populations they form, expressed statistically in terms of cloud-level distributions. Synthesizing cluster catalogues from an observed GMC catalogue in M83, we find that this model predicts initial star cluster masses and sizes that are in good agreement with observations, using only standard IMF and stellar evolution models as inputs for feedback. Within our model, the ratiomore »
-
Abstract Molecular clouds are supported by turbulence and magnetic fields, but quantifying their influence on cloud life cycle and star formation efficiency (SFE) remains an open question. We perform radiation magnetohydrodynamic simulations of star-forming giant molecular clouds (GMCs) with UV radiation feedback, in which the propagation of UV radiation via ray tracing is coupled to hydrogen photochemistry. We consider 10 GMC models that vary in either initial virial parameter (1 ≤ α vir,0 ≤ 5) or dimensionless mass-to-magnetic flux ratio (0.5 ≤ μ Φ,0 ≤ 8 and ∞ ); the initial mass 10 5 M ⊙ and radius 20 pc are fixed. Each model is run with five different initial turbulence realizations. In most models, the duration of star formation and the timescale for molecular gas removal (primarily by photoevaporation) are 4–8 Myr. Both the final SFE ( ε * ) and time-averaged SFE per freefall time ( ε ff ) are reduced by strong turbulence and magnetic fields. The median ε * ranges between 2.1% and 9.5%. The median ε ff ranges between 1.0% and 8.0%, and anticorrelates with α vir,0 , in qualitative agreement with previous analytic theory and simulations. However, the time-dependent α vir ( t )more »
-
Dimensional analysis suggests that the dissipation length scale ( $\ell _{{\it\epsilon}}=u_{\star }^{3}/{\it\epsilon}$ ) is the appropriate scale for the shear-production range of the second-order streamwise structure function in neutrally stratified turbulent shear flows near solid boundaries, including smooth- and rough-wall boundary layers and shear layers above canopies (e.g. crops, forests and cities). These flows have two major characteristics in common: (i) a single velocity scale, i.e. the friction velocity ( $u_{\star }$ ) and (ii) the presence of large eddies that scale with an external length scale much larger than the local integral length scale. No assumptions are made about the local integral scale, which is shown to be proportional to $\ell _{{\it\epsilon}}$ for the scaling analysis to be consistent with Kolmogorov’s result for the inertial subrange. Here ${\it\epsilon}$ is the rate of dissipation of turbulent kinetic energy (TKE) that represents the rate of energy cascade in the inertial subrange. The scaling yields a log-law dependence of the second-order streamwise structure function on ( $r/\ell _{{\it\epsilon}}$ ), where $r$ is the streamwise spatial separation. This scaling law is confirmed by large-eddy simulation (LES) results in the roughness sublayer above a model canopy, where the imbalance between local production and dissipation of TKEmore »
