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
Testing feedback from star clusters in simulations of the Milky Way formation
We present a suite of galaxy formation simulations that directly model star cluster formation and disruption. Starting from a model previously developed by our group, here we introduce several improvements to the prescriptions for cluster formation and feedback, then test these updates using a large suite of cosmological simulations of Milky Way mass galaxies. We perform a differential analysis with the goal of understanding how each of the updates affects star cluster populations. Two key parameters are the momentum boost of supernova feedback fboost and star formation efficiency per free-fall time ϵff. We find that fboost has a strong influence on the galactic star formation rate, with higher values leading to less star formation. The efficiency ϵff does not have a significant impact on the global star formation rate, but dramatically changes cluster properties, with increasing ϵff leading to a higher maximum cluster mass, shorter age spread of stars within clusters, and higher integrated star formation efficiencies. We also explore the redshift evolution of the observable cluster mass function, finding that most massive clusters have formed at high redshift z > 4. Extrapolation of cluster disruption to z = 0 produces good agreement with both the Galactic globular cluster mass function and age–metallicity relation. Our results emphasize the importance of using small-scale properties of galaxies to calibrate subgrid models of star cluster formation and feedback.
more » « less- Award ID(s):
- 1909063
- PAR ID:
- 10409859
- Publisher / Repository:
- Oxford University Press
- Date Published:
- Journal Name:
- Monthly Notices of the Royal Astronomical Society
- Volume:
- 514
- Issue:
- 1
- ISSN:
- 0035-8711
- Page Range / eLocation ID:
- p. 280-301
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
more » « less
ABSTRACT Recent observations indicate that galactic outflows are ubiquitous in high-redshift (high-z) galaxies, including normal star-forming galaxies, quasar hosts, and dusty star-forming galaxies (DSFGs). However, the impact of outflows on the evolution of their hosts is still an open question. Here, we analyse the star-formation histories and galactic outflow properties of galaxies in massive haloes ($10^{12}\, {\rm M}_{\odot }\ \lt\ M_{\rm vir}\ \lt\ 5\times 10^{12}\, {\rm M}_{\odot }$) at z ≳ 5.5 in three zoom-in cosmological simulations from the MassiveFIRE suite, as part of the Feedback In Realistic Environments (FIRE) project. The simulations were run with the FIRE-2 model, which does not include feedback from active galactic nuclei. The simulated galaxies resemble z > 4 DSFGs, with star-formation rates of $\sim\!{1000}\ {\rm M}_{\odot }\, \rm yr^{-1}$ and molecular gas masses of Mmol ∼ 1010 M⊙. However, the simulated galaxies are characterized by higher circular velocities than those observed in high-z DSFGs. The mass loading factors from stellar feedback are of the order of ∼0.1, implying that stellar feedback is inefficient in driving galactic outflows and gas is consumed by star formation on much shorter time-scales than it is expelled from the interstellar medium. We also find that stellar feedback is highly inefficient in self-regulating star formation in this regime, with an average integrated star formation efficiency (SFE) per dynamical time of 30 per cent. Finally, compared with FIRE-2 galaxies hosted in similarly massive haloes at lower redshift, we find lower mass loading factors and higher SFEs in the high-z sample. We argue that both effects originate from the higher total and gas surface densities that characterize high-z massive systems.
-
Abstract In a series of three papers, we introduced a novel cluster formation model that describes the formation, growth, and disruption of star clusters in high-resolution cosmological simulations. We tested this model on a Milky Way-sized galaxy and found that various properties of young massive clusters, such as the mass function and formation efficiency, are consistent with observations in the local universe. Interestingly, most massive clusters – globular cluster candidates – are preferentially formed during major merger events. We follow the dynamical evolution of clusters in the galactic tidal field. Due to tidal disruption, the cluster mass function evolves from initial power law to a peaked shape. The surviving clusters at z = 0 show a broad range of metallicity [Fe/H] from -3 to -0.5. A robust prediction of the model is the age–metallicity relation, in which metal-rich clusters are systematically younger than metal-poor clusters by up to 3 Gyr.more » « less
-
Abstract Dust temperature is an important property of the interstellar medium (ISM) of galaxies. It is required when converting (sub)millimetre broad-band flux to total infrared luminosity (LIR), and hence star formation rate, in high-redshift galaxies. However, different definitions of dust temperatures have been used in the literature, leading to different physical interpretations of how ISM conditions change with, e.g. redshift and star formation rate. In this paper, we analyse the dust temperatures of massive ($M_{\rm star} \gt 10^{10}\, \mathrm{M}_{\odot }$) $z$ = 2–6 galaxies with the help of high-resolution cosmological simulations from the Feedback in Realistic Environments (fire) project. At $z$ ∼ 2, our simulations successfully predict dust temperatures in good agreement with observations. We find that dust temperatures based on the peak emission wavelength increase with redshift, in line with the higher star formation activity at higher redshift, and are strongly correlated with the specific star formation rate. In contrast, the mass-weighted dust temperature, which is required to accurately estimate the total dust mass, does not strongly evolve with redshift over $z$ = 2–6 at fixed IR luminosity but is tightly correlated with LIR at fixed $z$. We also analyse an ‘equivalent’ dust temperature for converting (sub)millimetre flux density to total IR luminosity, and provide a fitting formula as a function of redshift and dust-to-metal ratio. We find that galaxies of higher equivalent (or higher peak) dust temperature (‘warmer dust’) do not necessarily have higher mass-weighted temperatures. A ‘two-phase’ picture for interstellar dust can explain the different scaling relations of the various dust temperatures.more » « less
-
Abstract Using stellar population synthesis models to infer star formation histories (SFHs), we analyze photometry and spectroscopy of a large sample of quiescent galaxies that are members of Sunyaev–Zel’dovich (SZ)-selected galaxy clusters across a wide range of redshifts. We calculate stellar masses and mass-weighted ages for 837 quiescent cluster members at 0.3 < z < 1.4 using rest-frame optical spectra and the Python-based Prospector framework, from 61 clusters in the SPT-GMOS Spectroscopic Survey (0.3 < z < 0.9) and three clusters in the SPT Hi-z cluster sample (1.25 < z < 1.4). We analyze spectra of subpopulations divided into bins of redshift, stellar mass, cluster mass, and velocity-radius phase-space location, as well as by creating composite spectra of quiescent member galaxies. We find that quiescent galaxies in our data set sample a diversity of SFHs, with a median formation redshift (corresponding to the lookback time from the redshift of observation to when a galaxy forms 50% of its mass, t 50 ) of z = 2.8 ± 0.5, which is similar to or marginally higher than that of massive quiescent field and cluster galaxy studies. We also report median age–stellar mass relations for the full sample (age of the universe at t 50 (Gyr) = 2.52 (±0.04)–1.66 (±0.12) log 10 ( M /10 11 M ⊙ )) and recover downsizing trends across stellar mass; we find that massive galaxies in our cluster sample form on aggregate ∼0.75 Gyr earlier than lower-mass galaxies. We also find marginally steeper age–mass relations at high redshifts, and report a bigger difference in formation redshifts across stellar mass for fixed environment, relative to formation redshifts across environment for fixed stellar mass.more » « less
