ABSTRACT We report the formation of bound star clusters in a sample of high-resolution cosmological zoom-in simulations of z ≥ 5 galaxies from the Feedback In Realistic Environments project. We find that bound clusters preferentially form in high-pressure clouds with gas surface densities over $10^4\, \mathrm{ M}_{\odot }\, {\rm pc}^{-2}$, where the cloud-scale star formation efficiency is near unity and young stars born in these regions are gravitationally bound at birth. These high-pressure clouds are compressed by feedback-driven winds and/or collisions of smaller clouds/gas streams in highly gas-rich, turbulent environments. The newly formed clusters follow a power-law mass function of dN/dM ∼ M−2. The cluster formation efficiency is similar across galaxies with stellar masses of ∼107–$10^{10}\, \mathrm{ M}_{\odot }$ at z ≥ 5. The age spread of cluster stars is typically a few Myr and increases with cluster mass. The metallicity dispersion of cluster members is ∼0.08 dex in $\rm [Z/H]$ and does not depend on cluster mass significantly. Our findings support the scenario that present-day old globular clusters (GCs) were formed during relatively normal star formation in high-redshift galaxies. Simulations with a stricter/looser star formation model form a factor of a few more/fewer bound clusters per stellar mass formed, while the shape of the mass function is unchanged. Simulations with a lower local star formation efficiency form more stars in bound clusters. The simulated clusters are larger than observed GCs due to finite resolution. Our simulations are among the first cosmological simulations that form bound clusters self-consistently in a wide range of high-redshift galaxies.
more »
« less
Stellar populations in STARFORGE: the origin and evolution of star clusters and associations
Most stars form in highly clustered environments within molecular clouds, but eventually disperse into the distributed stellar field population. Exactly how the stellar distribution evolves from the embedded stage into gas-free associations and (bound) clusters is poorly understood. We investigate the long-term evolution of stars formed in the starforge simulation suite – a set of radiation-magnetohydrodynamic simulations of star-forming turbulent clouds that include all key stellar feedback processes inherent to star formation. We use nbody6++gpu to follow the evolution of the young stellar systems after gas removal. We use HDBSCAN to define stellar groups and analyse the stellar kinematics to identify the true bound star clusters. The conditions modeled by the simulations, i.e. global cloud surface densities below 0.15 g cm−2, star formation efficiencies below 15 per cent, and gas expulsion time-scales shorter than a free fall time, primarily produce expanding stellar associations and small clusters. The largest star clusters, which have ∼1000 bound members, form in the densest and lowest velocity dispersion clouds, representing ∼32 and 39 per cent of the stars in the simulations, respectively. The cloud’s early dynamical state plays a significant role in setting the classical star formation efficiency versus bound fraction relation. All stellar groups follow a narrow mass-velocity dispersion power-law relation at 10 Myr with a power-law index of 0.21. This correlation result in a distinct mass–size relationship for bound clusters. We also provide valuable constraints on the gas dispersal time-scale during the star formation process and analyse the implications for the formation of bound systems.
more » « less- NSF-PAR ID:
- 10478392
- Publisher / Repository:
- Oxford University Press
- Date Published:
- Journal Name:
- Monthly Notices of the Royal Astronomical Society
- Volume:
- 527
- Issue:
- 3
- ISSN:
- 0035-8711
- Format(s):
- Medium: X Size: p. 6732-6751
- Size(s):
- ["p. 6732-6751"]
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
more » « less
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.
-
ABSTRACT We study the spatially resolved (sub-kpc) gas velocity dispersion (σ)–star formation rate (SFR) relation in the FIRE-2 (Feedback in Realistic Environments) cosmological simulations. We specifically focus on Milky Way-mass disc galaxies at late times (z ≈ 0). In agreement with observations, we find a relatively flat relationship, with σ ≈ 15–30 km s−1 in neutral gas across 3 dex in SFRs. We show that higher dense gas fractions (ratios of dense gas to neutral gas) and SFRs are correlated at constant σ. Similarly, lower gas fractions (ratios of gas to stellar mass) are correlated with higher σ at constant SFR. The limits of the σ–ΣSFR relation correspond to the onset of strong outflows. We see evidence of ‘on-off’ cycles of star formation in the simulations, corresponding to feedback injection time-scales of 10–100 Myr, where SFRs oscillate about equilibrium SFR predictions. Finally, SFRs and velocity dispersions in the simulations agree well with feedback-regulated and marginally stable gas disc (Toomre’s Q = 1) model predictions, and the simulation data effectively rule out models assuming that gas turns into stars at (low) constant efficiency (i.e. 1 per cent per free-fall time). And although the simulation data do not entirely exclude gas accretion/gravitationally powered turbulence as a driver of σ, it appears to be subdominant to stellar feedback in the simulated galaxy discs at z ≈ 0.more » « less
-
null (Ed.)ABSTRACT We study the role of group infall in the assembly and dynamics of galaxy clusters in ΛCDM. We select 10 clusters with virial mass M200 ∼ 1014 $\rm M_\odot$ from the cosmological hydrodynamical simulation Illustris and follow their galaxies with stellar mass M⋆ ≥ 1.5 × 108 $\rm M_\odot$. A median of ${\sim}38{{\ \rm per\ cent}}$ of surviving galaxies at z = 0 is accreted as part of groups and did not infall directly from the field, albeit with significant cluster-to-cluster scatter. The evolution of these galaxy associations is quick, with observational signatures of their common origin eroding rapidly in 1–3 Gyr after infall. Substructure plays a dominant role in fostering the conditions for galaxy mergers to happen, even within the cluster environment. Integrated over time, we identify (per cluster) an average of 17 ± 9 mergers that occur in infalling galaxy associations, of which 7 ± 3 occur well within the virial radius of their cluster hosts. The number of mergers shows large dispersion from cluster to cluster, with our most massive system having 42 mergers above our mass cut-off. These mergers, which are typically gas rich for dwarfs and a combination of gas rich and gas poor for M⋆ ∼ 1011 $\rm M_\odot$, may contribute significantly within ΛCDM to the formation of specific morphologies, such as lenticulars (S0) and blue compact dwarfs in groups and clusters.more » « less
-
more » « less
ABSTRACT We study the formation and early evolution of star clusters that have a wide range of masses and background cloud mass surface densities, Σcloud, which help set the initial sizes, densities, and velocity dispersions of the natal gas clumps. Initial clump masses of 300, 3000, and 30 000 M⊙ are considered, from which star clusters are born with an assumed 50 per cent overall star formation efficiency and with 50 per cent primordial binarity. This formation is gradual, i.e. with a range of star formation efficiencies per free-fall time from 1 to 100 per cent, so that the formation time can range from 0.7 Myr for low-mass, high-Σcloud clumps to ∼30 Myr for high-mass, low-Σcloud clumps. Within this framework of the turbulent clump model, for a given Σcloud, clumps of higher mass are of lower initial volume density, but their dynamical evolution leads to higher bound fractions and causes them to form much higher density cluster cores and maintain these densities for longer periods. This results in systematic differences in the evolution of binary properties, degrees of mass segregation, and rates of creation of dynamically ejected runaways. We discuss the implications of these results for observed star clusters and stellar populations.
