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.


Title: The contribution of globular clusters to cosmic reionization
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 find that the averaged fesc over the lifetime of a star particle follows a similar distribution for cluster stars and non-cluster stars. Clusters tend to have low fesc in the first few Myr, presumably because they form preferentially in more extreme environments with high optical depths; the fesc increases later as feedback starts to destroy the natal cloud. On the other hand, some non-cluster stars formed between cluster complexes or in the compressed shells at the front of a superbubble can also have high fesc. We find that cluster stars on average have comparable fesc to non-cluster stars. This result is robust across several star formation models in our simulations. Our results suggest that the fraction of ionizing photons from proto-GCs to cosmic reionization is comparable to the cluster formation efficiencies in high-redshift galaxies and thus proto-GCs likely contribute an appreciable fraction of photons but are not the dominant sources for reionization.  more » « less
Award ID(s):
1752913 1910346 1715216 1715070
PAR ID:
10258058
Author(s) / Creator(s):
; ; ; ;
Date Published:
Journal Name:
Monthly Notices of the Royal Astronomical Society
Volume:
504
Issue:
3
ISSN:
0035-8711
Page Range / eLocation ID:
4062 to 4071
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. null (Ed.)
    ABSTRACT We present the escape fraction of hydrogen ionizing photons (fesc) from a sample of 34 high-resolution cosmological zoom-in simulations of galaxies at z ≥ 5 in the Feedback in Realistic Environments project, post-processed with a Monte Carlo radiative transfer code for ionizing radiation. Our sample consists of 8500 haloes in Mvir ∼ 108–$$10^{12}\, M_{\odot }$$ (M* ∼ 104–$$10^{10}\, M_{\odot }$$) at z = 5–12. We find the sample average 〈fesc〉increases with halo mass for Mvir ∼ 108–$$10^{9.5}\, M_{\odot }$$, becomes nearly constant for 109.5–$$10^{11}\, M_{\odot }$$, and decreases at $${\gtrsim}10^{11}\, M_{\odot }$$. Equivalently, 〈fesc〉 increases with stellar mass up to $$M_{\ast }\sim 10^8\, M_{\odot }$$ and decreases at higher masses. Even applying single-star stellar population synthesis models, we find a moderate 〈fesc〉 ∼ 0.2 for galaxies at $$M_{\ast }\sim 10^8\, M_{\odot }$$. Nearly half of the escaped ionizing photons come from stars 1–3 Myr old and the rest from stars 3–10 Myr old. Binaries only have a modest effect, boosting 〈fesc〉 by ∼25–35 per cent and the number of escaped photons by 60–80 per cent. Most leaked ionizing photons are from vigorously star-forming regions that usually contain a feedback-driven kpc-scale superbubble surrounded by a dense shell. The shell is forming stars while accelerated, so new stars formed earlier in the shell are already inside the shell. Young stars in the bubble and near the edge of the shell can fully ionize some low-column-density paths pre-cleared by feedback, allowing a large fraction of their ionizing photons to escape. The decrease of 〈fesc〉 at the high-mass end is due to dust attenuation, while at the low-mass end, 〈fesc〉 decreases owing to inefficient star formation and hence feedback. At fixed mass, 〈fesc〉 tends to increase with redshift. Although the absolute 〈fesc〉does not fully converge with resolution in our simulations, the mass- and redshift-dependence of 〈fesc〉 is likely robust. Our simulations produce sufficient ionizing photons for cosmic reionization. 
    more » « less
  2. 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
  3. Abstract A fundamental requirement for reionizing the Universe is that a sufficient fraction of the ionizing photons emitted by galaxies successfully escapes into the intergalactic medium. However, due to the scarcity of high-redshift observational data, the sources driving reionization remain uncertain. In this work, we calculate the ionizing escape fractions (fesc) of reionization-era galaxies from the state-of-the-art thesan simulations, which combine an accurate radiation-hydrodynamic solver (arepo-rt) with the well-tested IllustrisTNG galaxy formation model to self-consistently simulate both small-scale galaxy physics and large-scale reionization throughout a large patch of the universe ($$L_\text{box} = 95.5\, \text{cMpc}$$). This allows the formation of numerous massive haloes ($$M_\text{halo} \gtrsim 10^{10}\, {\text{M}_{\odot }}$$), which are often statistically underrepresented in previous studies but are believed to be important to achieving rapid reionization. We find that low-mass galaxies ($$M_\text{stars} \lesssim 10^7\, {\text{M}_{\odot }}$$) are the main drivers of reionization above z ≳ 7, while high-mass galaxies ($$M_\text{stars} \gtrsim 10^8\, {\text{M}_{\odot }}$$) dominate the escaped ionizing photon budget at lower redshifts. We find a strong dependence of fesc on the effective star formation rate (SFR) surface density defined as the SFR per gas mass per escape area, i.e. $$\bar{\Sigma }_\text{SFR} = \text{SFR}/M_\text{gas}/R_{200}^2$$. The variation in halo escape fractions decreases for higher mass haloes, which can be understood from the more settled galactic structure, SFR stability, and fraction of sightlines within each halo significantly contributing to the escaped flux. Dust is capable of reducing the escape fractions of massive galaxies, but the impact on the global fesc depends on the dust model. Finally, active galactic nuclei are unimportant for reionization in thesan and their escape fractions are lower than stellar ones due to being located near the centres of galaxy gravitational potential wells. 
    more » « less
  4. ABSTRACT We perform cosmological hydrodynamical simulations to study the formation of proto-globular cluster candidates in progenitors of present-day dwarf galaxies $$(M_{\rm vir} \approx 10^{10}\, {\rm M}_\odot$$ at z = 0) as part of the ‘Feedback in Realistic Environment’ (FIRE) project. Compact (r1/2 < 30 pc), relatively massive (0.5 × 105 ≲ M⋆/M⊙ ≲ 5 × 105), self-bound stellar clusters form at 11 ≳ z ≳ 5 in progenitors with $$M_{\rm vir} \approx 10^9\, \mathrm{M}_{\odot }$$. Cluster formation is triggered when at least $$10^7\, \mathrm{M}_{\odot }$$ of dense, turbulent gas reaches $$\Sigma _{\rm gas} \approx 10^4\, {\rm M}_\odot \, {\rm pc}^{-2}$$ as a result of the compressive effects of supernova feedback or from cloud–cloud collisions. The clusters can survive for $$2-3\, {\rm Gyr}$$; absent numerical effects, they could possibly survive substantially longer, perhaps to z = 0. The longest lived clusters are those that form at significant distance – several hundreds of pc – from their host galaxy. We therefore predict that globular clusters forming in progenitors of present-day dwarf galaxies will be offset from any pre-existing stars within their host dark matter haloes as opposed to deeply embedded within a well-defined galaxy. Properties of the nascent clusters are consistent with observations of some of the faintest and most compact high-redshift sources in Hubble Space Telescope lensing fields and are at the edge of what will be detectable as point sources in deep imaging of non-lensed fields with JWST. By contrast, the star clusters’ host galaxies will remain undetectable. 
    more » « less
  5. ABSTRACT The existence of globular clusters (GCs) in a few satellite galaxies, and their absence in majority of dwarf galaxies, present a challenge for models attempting to understand the origins of GCs. In addition to GC presence appearing stochastic and difficult to describe with average trends, in the smallest satellite galaxies GCs contribute a substantial fraction of total stellar mass. We investigate the stochasticity and number of GCs in dwarf galaxies using an updated version of our model that links the formation of GCs to the growth of the host galaxy mass. We find that more than 50 per cent of dwarf galaxies with stellar mass $$M_{\star }\lesssim 2\times 10^7\, \mathrm{M}_\odot$$ do not host GCs, whereas dwarfs with $$M_{\star }\sim 10^8\, \mathrm{M}_\odot$$ almost always contain some GCs, with a median number ∼10 at z  = 0. These predictions are in agreement with the observations of the Local Volume dwarfs. We also confirm the near-linear GC system mass–halo mass relation down to $$M_{\mathrm{h}}\simeq 10^8\, \mathrm{M}_\odot$$ under the assumption that GC formation and evolution in galaxies of all mass can be described by the same physical model. A detailed case study of two model dwarfs that resemble the Fornax dwarf spheroidal galaxy shows that observational samples can be notably biased by incompleteness below detection limit and at large radii. 
    more » « less