Self-consistent proto-globular cluster formation in cosmological simulations of high-redshift galaxies
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 more »
Authors:
; ; ; ; ; ; ; ; ;
Award ID(s):
Publication Date:
NSF-PAR ID:
10173303
Journal Name:
Monthly Notices of the Royal Astronomical Society
Volume:
493
Issue:
3
Page Range or eLocation-ID:
4315 to 4332
ISSN:
0035-8711
1. ABSTRACT We study the escape fraction of ionizing photons (fesc) in two cosmological zoom-in simulations of galaxies in the reionization era with halo mass Mhalo ∼ 1010 and $10^{11}\, \mathrm{ M}_{\odot }$ (stellar mass M* ∼ 107 and $10^9\, \mathrm{ M}_{\odot }$) at z = 5 from the Feedback in Realistic Environments project. These simulations explicitly resolve the formation of proto-globular clusters (GCs) self-consistently, where 17–39 per cent of stars form in bound clusters during starbursts. Using post-processing Monte Carlo radiative transfer calculations of ionizing radiation, we compute fesc from cluster stars and non-cluster stars formed during a starburst over ∼100 Myr in each galaxy. We findmore »
We study a suite of extremely high-resolution cosmological Feedback in Realistic Environments simulations of dwarf galaxies ($M_{\rm halo} \lesssim 10^{10}\rm \, M_{\odot }$), run to z = 0 with $30\, \mathrm{M}_{\odot }$ resolution, sufficient (for the first time) to resolve the internal structure of individual supernovae remnants within the cooling radius. Every halo with $M_{\rm halo} \gtrsim 10^{8.6}\, \mathrm{M}_{\odot }$ is populated by a resolved stellar galaxy, suggesting very low-mass dwarfs may be ubiquitous in the field. Our ultra-faint dwarfs (UFDs; $M_{\ast }\lt 10^{5}\, \mathrm{M}_{\odot }$) have their star formation (SF) truncated early (z ≳ 2), likely by reionization, whilemore »
3. ABSTRACT We study the growth of stellar discs of Milky Way-sized galaxies using a suite of cosmological simulations. We calculate the half-mass axis lengths and axis ratios of stellar populations split by age in galaxies with stellar mass $M_{*}=10^7\!-\!10^{10}\, \mathrm{M}_{\odot }$ at redshifts z > 1.5. We find that in our simulations stars always form in relatively thin discs, and at ages below 100 Myr are contained within half-mass height z1/2 ∼ 0.1 kpc and short-to-long axial ratio z1/2/x1/2 ∼ 0.15. Disc thickness increases with the age of stellar population, reaching median z1/2 ∼ 0.8 kpc and z1/2/x1/2 ∼ 0.6 for stars oldermore »
4. ABSTRACT We present and study a large suite of high-resolution cosmological zoom-in simulations, using the FIRE-2 treatment of mechanical and radiative feedback from massive stars, together with explicit treatment of magnetic fields, anisotropic conduction and viscosity (accounting for saturation and limitation by plasma instabilities at high β), and cosmic rays (CRs) injected in supernovae shocks (including anisotropic diffusion, streaming, adiabatic, hadronic and Coulomb losses). We survey systems from ultrafaint dwarf ($M_{\ast }\sim 10^{4}\, \mathrm{M}_{\odot }$, $M_{\rm halo}\sim 10^{9}\, \mathrm{M}_{\odot }$) through Milky Way/Local Group (MW/LG) masses, systematically vary uncertain CR parameters (e.g. the diffusion coefficient κ and streaming velocity), andmore »
We use the simba cosmological galaxy formation simulation to investigate the relationship between major mergers ($\lesssim$4:1), starbursts, and galaxy quenching. Mergers are identified via sudden jumps in stellar mass M* well above that expected from in situ star formation, while quenching is defined as going from specific star formation rate (sSFR) $\gt t_{\rm H}^{-1}$ to $\lt 0.2t_{\rm H}^{-1}$, where tH is the Hubble time. At z ≈ 0–3, mergers show ∼2–3× higher SFR than a mass-matched sample of star-forming galaxies, but globally represent $\lesssim 1{{\ \rm per\ cent}}$ of the cosmic SF budget. At low masses, the increase in SFRmore »