skip to main content


Title: Simulating the effect of photoheating feedback during reionization
Abstract We present self-consistent radiation hydrodynamic simulations of hydrogen reionization performed with arepo-rt complemented by a state-of-the-art galaxy formation model. We examine how photoheating feedback, due to reionization, shapes the galaxies properties. Our fiducial model completes reionization by z ≈ 6 and matches observations of the Ly α forest, the cosmic microwave background electron scattering optical depth, the high-redshift ultraviolet (UV) luminosity function, and stellar mass function. Contrary to previous works, photoheating suppresses star formation rates by more than $50{{\ \rm per\ cent}}$ only in haloes less massive than ∼108.4 M⊙ (∼108.8 M⊙) at z = 6 (z = 5), suggesting inefficient photoheating feedback from photons within galaxies. The use of a uniform UV background that heats up the gas at z ≈ 10.7 generates an earlier onset of suppression of star formation compared to our fiducial model. This discrepancy can be mitigated by adopting a UV background model with a more realistic reionization history. In the absence of stellar feedback, photoheating alone is only able to quench haloes less massive than ∼109 M⊙ at z ≳ 5, implying that photoheating feedback is sub-dominant in regulating star formation. In addition, stellar feedback, implemented as a non-local galactic wind scheme in the simulations, weakens the strength of photoheating feedback by reducing the amount of stellar sources. Most importantly, photoheating does not leave observable imprints in the UV luminosity function, stellar mass function, or the cosmic star formation rate density. The feasibility of using these observables to detect imprints of reionization therefore requires further investigation.  more » « less
Award ID(s):
1814053 1814259
PAR ID:
10104826
Author(s) / Creator(s):
; ; ; ;
Date Published:
Journal Name:
Monthly Notices of the Royal Astronomical Society
Volume:
488
Issue:
1
ISSN:
0035-8711
Page Range / eLocation ID:
419 to 437
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. null (Ed.)
    ABSTRACT The James Webb Space Telescope (JWST) is expected to observe galaxies at z > 10 that are presently inaccessible. Here, we use a self-consistent empirical model, the universemachine, to generate mock galaxy catalogues and light-cones over the redshift range z = 0−15. These data include realistic galaxy properties (stellar masses, star formation rates, and UV luminosities), galaxy–halo relationships, and galaxy–galaxy clustering. Mock observables are also provided for different model parameters spanning observational uncertainties at z < 10. We predict that Cycle 1 JWST surveys will very likely detect galaxies with M* > 107 M⊙ and/or M1500 < −17 out to at least z ∼ 13.5. Number density uncertainties at z > 12 expand dramatically, so efforts to detect z > 12 galaxies will provide the most valuable constraints on galaxy formation models. The faint-end slopes of the stellar mass/luminosity functions at a given mass/luminosity threshold steepen as redshift increases. This is because observable galaxies are hosted by haloes in the exponentially falling regime of the halo mass function at high redshifts. Hence, these faint-end slopes are robustly predicted to become shallower below current observable limits (M* < 107 M⊙ or M1500 > −17). For reionization models, extrapolating luminosity functions with a constant faint-end slope from M1500 = −17 down to M1500 = −12 gives the most reasonable upper limit for the total UV luminosity and cosmic star formation rate up to z ∼ 12. We compare to three other empirical models and one semi-analytic model, showing that the range of predicted observables from our approach encompasses predictions from other techniques. Public catalogues and light-cones for common fields are available online. 
    more » « less
  2. ABSTRACT

    We measure escape fractions, fesc, of ionizing radiation from galaxies in the sphinx suite of cosmological radiation-hydrodynamical simulations of reionization, resolving haloes with $M_{\rm vir}\gtrsim 7.5 \times 10^7 \ {\rm {M}_{\odot }}$ with a minimum cell width of ≈10 pc. Our new and largest 20 co-moving Mpc wide volume contains tens of thousands of star-forming galaxies with halo masses up to a few times 1011 M⊙. The simulated galaxies agree well with observational constraints of the ultraviolet (UV) luminosity function in the Epoch of Reionization. The escape fraction fluctuates strongly in individual galaxies over time-scales of a few Myr, due to its regulation by supernova and radiation feedback, and at any given time a tiny fraction of star-forming galaxies emits a large fraction of the ionizing radiation escaping into the intergalactic medium. Statistically, fesc peaks in intermediate-mass, intermediate-brightness, and low-metallicity galaxies (M* ≈ 107 M⊙, M1500 ≈ −17, Z ≲ 5 × 10−3 Z⊙), dropping strongly for lower and higher masses, brighter and dimmer galaxies, and more metal-rich galaxies. The escape fraction correlates positively with both the short-term and long-term specific star formation rate. According to sphinx, galaxies too dim to be yet observed, with ${M_{1500}}\gtrsim -17$, provide about 55 per cent of the photons contributing to reionization. The global averaged fesc naturally decreases with decreasing redshift, as predicted by UV background models and low-redshift observations. This evolution is driven by decreasing specific star formation rates over cosmic time.

     
    more » « less
  3. ABSTRACT We present an updated model of the cosmic ionizing background from the UV to the X-rays. Relative to our previous model, the new model provides a better match to a large number of up-to-date empirical constraints, including: (1) new galaxy and AGN luminosity functions; (2) stellar spectra including binary stars; (3) obscured and unobscured AGN; (4) a measurement of the non-ionizing UV background; (5) measurements of the intergalactic H i and He ii photoionization rates at z ∼ 0−6; (6) the local X-ray background; and (7) improved measurements of the intergalactic opacity. In this model, AGN dominate the H i ionizing background at z ≲ 3 and star-forming galaxies dominate it at higher redshifts. Combined with the steeply declining AGN luminosity function beyond z ∼ 2, the slow evolution of the H i ionization rate inferred from the high-redshift H i Ly α forest requires an escape fraction from star-forming galaxies that increases with redshift (a population-averaged escape fraction of $\approx 1{{\ \rm per\ cent}}$ suffices to ionize the intergalactic medium at z = 3 when including the contribution from AGN). We provide effective photoionization and photoheating rates calibrated to match the Planck 2018 reionization optical depth and recent constraints from the He ii Ly α forest in hydrodynamic simulations. 
    more » « less
  4. ABSTRACT

    We introduce the Astrid  simulation, a large-scale cosmological hydrodynamic simulation in a $250 \, h^{-1}\mathrm{Mpc}$ box with 2 × 55003 particles. Astrid contains a large number of high redshift galaxies, which can be compared to future survey data, and resolves galaxies in haloes more massive than $2\times 10^9 \, \mathrm{M}_{\odot }$. Astrid  has been run from z = 99 to 3. As a particular focus is modelling the high redshift Universe, it contains models for inhomogeneous hydrogen and helium reionization, baryon relative velocities and massive neutrinos, as well as supernova and AGN feedback. The black hole model includes mergers driven by dynamical friction rather than repositioning. We briefly summarize the implemented models, and the technical choices we took when developing the simulation code. We validate the model, showing good agreement with observed ultraviolet luminosity functions, galaxy stellar mass functions and specific star formation rates (SFRs). We show that the redshift at which a given galaxy underwent hydrogen reionization has a large effect on the halo gas fraction. Finally, at z = 6, haloes with $M \sim 2\times 10^9 \, \mathrm{M}_{\odot }$ which have been reionized have an SFR 1.5 times greater than those which have not yet been reionized.

     
    more » « less
  5. ABSTRACT

    The feedback loop between the galaxies producing the background radiation field for reionization and their growth is crucial, particularly for low-mass haloes. Despite this, the vast majority of galaxy formation studies employ a spatially uniform, time-varying reionizing background, with the majority of reionization studies employing galaxy formation models only required to work at high redshift. This paper uses the well-studied TNG galaxy formation model, calibrated at low redshift, coupled to the arepo-rt code, to self-consistently solve the coupled problems of galaxy evolution and reionization, evaluating the impact of patchy (and slow) reionization on early galaxies. thesan-hr is an extension of the thesan project to higher resolution (a factor of 50 increase, with a baryonic mass of mb ≈ 104 M⊙), to additionally enable the study of ‘mini-haloes’ with virial temperatures Tvir < 104 K. Comparing the self-consistent model to a uniform UV background, we show that galaxies in thesan-hr are predicted to be larger in physical extent (by a factor ∼2), less metal enriched (by ∼0.2 dex), and less abundant (by a factor ∼10 at M1500 =   − 10) by z = 5. We show that differences in star formation and enrichment patterns lead to significantly different predictions for star formation in low mass haloes, low-metallicity star formation, and even the occupation fraction of haloes. We posit that cosmological galaxy formation simulations aiming to study early galaxy formation (z ≳ 3) must employ a spatially inhomogeneous UV background to accurately reproduce galaxy properties.

     
    more » « less