skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.

Explore Research Products in the PAR
It may take a few hours for recently added research products to appear in PAR search results.


Title: A model for the formation of stellar associations and clusters from giant molecular clouds
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 ratio of the strength of gravity to stellar feedback is the key parameter setting the masses of star clusters, and of the various feedback channels direct stellar radiation (photon momentum and photoionization) is the most important on GMC scales.  more » « less
Award ID(s):
1652522 1910346 1752913
PAR ID:
10278870
Author(s) / Creator(s):
; ; ; ; ; ;
Date Published:
Journal Name:
Monthly Notices of the Royal Astronomical Society
ISSN:
0035-8711
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. ; ; ; ; (, Monthly Notices of the Royal Astronomical Society)
    null (Ed.)
    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 and SFE similar to star-forming GMCs in our Galaxy, we predict M50 to be $$\gt 20 \, \mathrm{M}_{\odot }$$, an order of magnitude larger than observed ($$\sim 2 \, \mathrm{M}_\odot$$), together with an excess of brown dwarfs. Moreover, M50 is sensitive to initial cloud properties and evolves strongly in time within a given cloud, predicting much larger IMF variations than are observationally allowed. We conclude that physics beyond MHD turbulence and gravity are necessary ingredients for the IMF. 
    more » « less
  2. ; (, Monthly Notices of the Royal Astronomical Society)
    Abstract Most simulations of galaxies and massive giant molecular clouds (GMCs) cannot explicitly resolve the formation (or predict the main-sequence masses) of individual stars. So they must use some prescription for the amount of feedback from an assumed population of massive stars (e.g. sampling the initial mass function, IMF). We perform a methods study of simulations of a star-forming GMC with stellar feedback from UV radiation, varying only the prescription for determining the luminosity of each stellar mass element formed (according to different IMF sampling schemes). We show that different prescriptions can lead to widely varying (factor of ∼3) star formation efficiencies (on GMC scales) even though the average mass-to-light ratios agree. Discreteness of sources is important: radiative feedback from fewer, more-luminous sources has a greater effect for a given total luminosity. These differences can dominate over other, more widely recognized differences between similar literature GMC-scale studies (e.g. numerical methods, cloud initial conditions, presence of magnetic fields). Moreover the differences in these methods are not purely numerical: some make different implicit assumptions about the nature of massive star formation, and this remains deeply uncertain in star formation theory. 
    more » « less
  3. ; ; ; ; ; ; ; (, Monthly Notices of the Royal Astronomical Society)
    ABSTRACT Giant molecular clouds (GMCs) are well studied in the local Universe, however, exactly how their properties vary during galaxy evolution is poorly understood due to challenging resolution requirements, both observational and computational. We present the first time-dependent analysis of GMCs in a Milky Way-like galaxy and an Large Magellanic Cloud (LMC)-like dwarf galaxy of the FIRE-2 (Feedback In Realistic Environments) simulation suite, which have sufficient resolution to predict the bulk properties of GMCs in cosmological galaxy formation self-consistently. We show explicitly that the majority of star formation outside the galactic centre occurs within self-gravitating gas structures that have properties consistent with observed bound GMCs. We find that the typical cloud bulk properties such as mass and surface density do not vary more than a factor of 2 in any systematic way after the first Gyr of cosmic evolution within a given galaxy from its progenitor. While the median properties are constant, the tails of the distributions can briefly undergo drastic changes, which can produce very massive and dense self-gravitating gas clouds. Once the galaxy forms, we identify only two systematic trends in bulk properties over cosmic time: a steady increase in metallicity produced by previous stellar populations and a weak decrease in bulk cloud temperatures. With the exception of metallicity, we find no significant differences in cloud properties between the Milky Way-like and dwarf galaxies. These results have important implications for cosmological star and star cluster formation and put especially strong constraints on theories relating the stellar initial mass function to cloud properties. 
    more » « less
  4. ; ; ; ; ; (, Monthly Notices of the Royal Astronomical Society)
    Abstract Measurements of the star formation efficiency (SFE) of giant molecular clouds (GMCs) in the Milky Way generally show a large scatter, which could be intrinsic or observational. We use magnetohydrodynamic simulations of GMCs (including feedback) to forward-model the relationship between the true GMC SFE and observational proxies. We show that individual GMCs trace broad ranges of observed SFE throughout collapse, star formation, and disruption. Low measured SFEs ($${\ll} 1\hbox{ per cent}$$) are ‘real’ but correspond to early stages; the true ‘per-freefall’ SFE where most stars actually form can be much larger. Very high ($${\gg} 10\hbox{ per cent}$$) values are often artificially enhanced by rapid gas dispersal. Simulations including stellar feedback reproduce observed GMC-scale SFEs, but simulations without feedback produce 20× larger SFEs. Radiative feedback dominates among mechanisms simulated. An anticorrelation of SFE with cloud mass is shown to be an observational artefact. We also explore individual dense ‘clumps’ within GMCs and show that (with feedback) their bulk properties agree well with observations. Predicted SFEs within the dense clumps are ∼2× larger than observed, possibly indicating physics other than feedback from massive (main-sequence) stars is needed to regulate their collapse. 
    more » « less
  5. ; ; ; ; ; ; ; (, Monthly Notices of the Royal Astronomical Society)
    ABSTRACT The properties of young star clusters formed within a galaxy are thought to vary in different interstellar medium conditions, but the details of this mapping from galactic to cluster scales are poorly understood due to the large dynamic range involved in galaxy and star cluster formation. We introduce a new method for modelling cluster formation in galaxy simulations: mapping giant molecular clouds (GMCs) formed self-consistently in a FIRE-2 magnetohydrodynamic galaxy simulation on to a cluster population according to a GMC-scale cluster formation model calibrated to higher resolution simulations, obtaining detailed properties of the galaxy’s star clusters in mass, metallicity, space, and time. We find $$\sim 10{{\ \rm per\ cent}}$$ of all stars formed in the galaxy originate in gravitationally bound clusters overall, and this fraction increases in regions with elevated Σgas and ΣSFR, because such regions host denser GMCs with higher star formation efficiency. These quantities vary systematically over the history of the galaxy, driving variations in cluster formation. The mass function of bound clusters varies – no single Schechter-like or power-law distribution applies at all times. In the most extreme episodes, clusters as massive as 7 × 106 M⊙ form in massive, dense clouds with high star formation efficiency. The initial mass–radius relation of young star clusters is consistent with an environmentally dependent 3D density that increases with Σgas and ΣSFR. The model does not reproduce the age and metallicity statistics of old ($$\gt 11\rm Gyr$$) globular clusters found in the Milky Way, possibly because it forms stars more slowly at z > 3. 
    more » « less
  • Citation Formats
  • MLA
  • APA
  • Chicago
  • BibTeX
  • Save / Share this Record
  • Send to Email