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: On the rise and fall of galactic ionizing output at the end of reionization
ABSTRACT Quasar absorption spectra measurements suggest that reionization proceeded rapidly, ended late at z ∼ 5.5, and was followed by a flat ionizing background evolution. Simulations that reproduce this behaviour often rely on a fine-tuned galaxy ionizing emissivity, which peaks at z ∼ 6–7 and drops a factor of 1.5–2.5 by z ∼ 5. This is puzzling since the abundance of galaxies is observed to grow monotonically during this period. Explanations for this include effects such as dust obscuration of ionizing photon escape and feedback from photoheating of the IGM. We explore the possibility that this drop in emissivity is instead an artefact of one or more modelling deficiencies in reionization simulations. These include possibly incorrect assumptions about the ionizing spectrum and/or inaccurate modelling of IGM clumping. Our results suggest that the need for a drop could be alleviated if simulations are underestimating the IGM opacity from massive, star-forming haloes. Other potential modelling issues either have a small effect or require a steeper drop when remedied. We construct an illustrative model in which the emissivity is nearly flat at reionization’s end, evolving only ∼0.05 dex at 5 < z < 7. More realistic scenarios, however, require a ∼0.1–0.3 dex drop. We also study the evolution of the Ly α effective optical depth distribution and compare to recent measurements. We find that models that feature a hard ionizing spectrum and/or are driven by faint, low-bias sources most easily reproduce the mean transmission and optical depth distribution of the forest simultaneously.  more » « less
Award ID(s):
2045600
PAR ID:
10517081
Author(s) / Creator(s):
; ; ; ;
Publisher / Repository:
MNRAS
Date Published:
Journal Name:
Monthly Notices of the Royal Astronomical Society
Volume:
531
Issue:
1
ISSN:
0035-8711
Page Range / eLocation ID:
1951 to 1970
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. ABSTRACT The high-redshift intergalactic medium (IGM) and the primeval galaxy population are rapidly becoming the new frontier of extragalactic astronomy. We investigate the IGM properties and their connection to galaxies at z ≥ 5.5 under different assumptions for the ionizing photon escape and the nature of dark matter, employing our novel thesan radiation-hydrodynamical simulation suite, designed to provide a comprehensive picture of the emergence of galaxies in a full reionization context. Our simulations have realistic ‘late’ reionization histories, match available constraints on global IGM properties, and reproduce the recently observed rapid evolution of the mean free path of ionizing photons. We additionally examine high-z Lyman-α transmission. The optical depth evolution is consistent with data, and its distribution suggests an even-later reionization than simulated, although with a strong sensitivity to the source model. We show that the effects of these two unknowns can be disentangled by characterizing the spectral shape and separation of Lyman-α transmission regions, opening up the possibility to observationally constrain both. For the first time in simulations, thesan reproduces the modulation of the Lyman-α flux as a function of galaxy distance, demonstrating the power of coupling a realistic galaxy formation model with proper radiation hydrodynamics. We find this feature to be extremely sensitive on the timing of reionization, while being relatively insensitive to the source model. Overall, thesan produces a realistic IGM and galaxy population, providing a robust framework for future analysis of the high-z Universe. 
    more » « less
  2. Abstract It has been claimed that traditional models struggle to explain the tentative detection of the 21 cm absorption trough centered atz∼ 17 measured by the EDGES collaboration. On the other hand, it has been shown that the EDGES results are consistent with an extrapolation of a declining UV luminosity density, following a simple power law of deep Hubble Space Telescope observations of 4 <z< 9 galaxies. We here explore the conditions by which the EDGES detection is consistent with current reionization and post-reionization observations, including the neutral hydrogen fraction atz∼ 6–8, Thomson-scattering optical depth, and ionizing emissivity atz∼ 5. By coupling a physically motivated source model derived from radiative transfer hydrodynamic simulations of reionization to a Markov Chain Monte Carlo sampler, we find that it is entirely possible to reconcile existing high-redshift (cosmic dawn) and low-redshift (reionization) constraints. In particular, we find that high contributions from low-mass halos along with high photon escape fractions are required to simultaneously reproduce cosmic dawn and reionization constraints. Our analysis further confirms that low-mass galaxies produce a flatter emissivity evolution, which leads to an earlier onset of reionization with a gradual and longer duration, resulting in a higher optical depth. While the models dominated by faint galaxies successfully reproduce the measured globally averaged quantities over the first one billion years, they underestimate the late redshift-instantaneous measurements in efficiently star-forming and massive systems. We show that our (simple) physically motivated semianalytical prescription produces results that are consistent with the (sophisticated) state-of-the-artTHESANradiation-magnetohydrodynamic simulation of the reionization. 
    more » « less
  3. ABSTRACT We present the Sherwood–Relics simulations, a new suite of large cosmological hydrodynamical simulations aimed at modelling the intergalactic medium (IGM) during and after the cosmic reionization of hydrogen. The suite consists of over 200 simulations that cover a wide range of astrophysical and cosmological parameters. It also includes simulations that use a new lightweight hybrid scheme for treating radiative transfer effects. This scheme follows the spatial variations in the ionizing radiation field, as well as the associated fluctuations in IGM temperature and pressure smoothing. It is computationally much cheaper than full radiation hydrodynamics simulations, and circumvents the difficult task of calibrating a galaxy formation model to observational constraints on cosmic reionization. Using this hybrid technique, we study the spatial fluctuations in IGM properties that are seeded by patchy cosmic reionization. We investigate the relevant physical processes and assess their impact on the z > 4 Lyman-α forest. Our main findings are: (i) consistent with previous studies patchy reionization causes large-scale temperature fluctuations that persist well after the end of reionization, (ii) these increase the Lyman-α forest flux power spectrum on large scales, and (iii) result in a spatially varying pressure smoothing that correlates well with the local reionization redshift. (iv) Structures evaporated or puffed up by photoheating cause notable features in the Lyman-α forest, such as flat-bottom or double-dip absorption profiles. 
    more » « less
  4. ABSTRACT The mean free path of ionizing photons, λmfp, is a key factor in the photoionization of the intergalactic medium (IGM). At z ≳ 5, however, λmfp may be short enough that measurements towards QSOs are biased by the QSO proximity effect. We present new direct measurements of λmfp that address this bias and extend up to z ∼ 6 for the first time. Our measurements at z ∼ 5 are based on data from the Giant Gemini GMOS survey and new Keck LRIS observations of low-luminosity QSOs. At z ∼ 6 we use QSO spectra from Keck ESI and VLT X-Shooter. We measure $$\lambda _{\rm mfp} = 9.09^{+1.62}_{-1.28}$$ proper Mpc and $$0.75^{+0.65}_{-0.45}$$ proper Mpc (68 per cent confidence) at z = 5.1 and 6.0, respectively. The results at z = 5.1 are consistent with existing measurements, suggesting that bias from the proximity effect is minor at this redshift. At z = 6.0, however, we find that neglecting the proximity effect biases the result high by a factor of two or more. Our measurement at z = 6.0 falls well below extrapolations from lower redshifts, indicating rapid evolution in λmfp over 5 < z < 6. This evolution disfavours models in which reionization ended early enough that the IGM had time to fully relax hydrodynamically by z = 6, but is qualitatively consistent with models wherein reionization completed at z = 6 or even significantly later. Our mean free path results are most consistent with late reionization models wherein the IGM is still 20 per cent neutral at z = 6, although our measurement at z = 6.0 is even lower than these models prefer. 
    more » « less
  5. ABSTRACT We compare a sample of five high-resolution, high S/N  Ly α forest spectra of bright 6 < z < ∼6.5 QSOs aimed at spectrally resolving the last remaining transmission spikes at z > 5 with those obtained from mock absorption spectra from the Sherwoodand Sherwood–Relics simulation suites of hydrodynamical simulations of the intergalactic medium (IGM). We use a profile-fitting procedure for the inverted transmitted flux, 1 − F, similar to the widely used Voigt profile fitting of the transmitted flux F at lower redshifts, to characterize the transmission spikes that probe predominately underdense regions of the IGM. We are able to reproduce the width and height distributions of the transmission spikes, both with optically thin simulations of the post-reionization Universe using a homogeneous UV background and full radiative transfer simulations of a late reionization model. We find that the width of the fitted components of the simulated transmission spikes is very sensitive to the instantaneous temperature of the reionized IGM. The internal structures of the spikes are more prominent in low temperature models of the IGM. The width distribution of the observed transmission spikes, which require high spectral resolution (≤ 8  km s−1) to be resolved, is reproduced for optically thin simulations with a temperature at mean density of T0 = (11 000 ± 1600, 10 500 ± 2100, 12 000 ± 2200) K at z = (5.4, 5.6, 5.8). This is weakly dependent on the slope of the temperature-density relation, which is favoured to be moderately steeper than isothermal. In the inhomogeneous, late reionization, full radiative transfer simulations where islands of neutral hydrogen persist to z ∼ 5.3, the width distribution of the observed transmission spikes is consistent with the range of T0 caused by spatial fluctuations in the temperature–density relation. 
    more » « less