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 morphology of reionization in a dynamically clumpy universe
ABSTRACT A recent measurement of the Lyman-limit mean free path at z = 6 suggests it may have been very short, motivating a better understanding of the role that ionizing photon sinks played in reionization. Accurately modelling the sinks in reionization simulations is challenging because of the large dynamic range required if ∼104−108M⊙ gas structures contributed significant opacity. Thus, there is no consensus on how important the sinks were in shaping reionization’s morphology. We address this question with a recently developed radiative transfer code that includes a dynamical sub-grid model for the sinks based on radiative hydrodynamics simulations. Compared to assuming a fully pressure-smoothed intergalactic medium, our dynamical treatment reduces ionized bubble sizes by $$10-20~{{\ \rm per\ cent}}$$ under typical assumptions about reionization’s sources. Near reionization’s midpoint, the 21 cm power at k ∼ 0.1 hMpc−1 is similarly reduced. These effects are more modest than the $$30-60~{{\ \rm per\ cent}}$$ suppression resulting from the higher recombination rate if pressure smoothing is neglected entirely. Whether the sinks played a significant role in reionization’s morphology depends on the nature of its sources. For example, if reionization was driven by bright (MUV < −17) galaxies, the sinks reduce the large-scale 21 cm power by at most 20  per cent, even if pressure smoothing is neglected. Conveniently, when bright sources contribute significantly, the morphology in our dynamical treatment can be reproduced accurately with a uniform sub-grid clumping factor that yields the same ionizing photon budget. By contrast, if MUV ∼ −13 galaxies drove reionization, the uniform clumping model can err by up to 40  per cent.  more » « less
Award ID(s):
2045600
PAR ID:
10429294
Author(s) / Creator(s):
; ; ;
Date Published:
Journal Name:
Monthly Notices of the Royal Astronomical Society
Volume:
522
Issue:
2
ISSN:
0035-8711
Page Range / eLocation ID:
2047 to 2064
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract Reionization leads to large spatial fluctuations in the intergalactic temperature that can persist well after its completion. We study the imprints of such fluctuations on the $$z$$ ∼ 5 Ly α forest flux power spectrum using a set of radiation-hydrodynamic simulations that model different reionization scenarios. We find that large-scale coherent temperature fluctuations bring $${\sim}20\text{--}60{{\ \rm per\ cent}}$$ extra power at k ∼ 0.002 s km−1, with the largest enhancements in the models where reionization is extended or ends the latest. On smaller scales (k ≳ 0.1 s km−1), we find that temperature fluctuations suppress power by $${\lesssim}10{{\ \rm per\ cent}}$$. We find that the shape of the power spectrum is mostly sensitive to the reionization mid-point rather than temperature fluctuations from reionization’s patchiness. However, for all of our models with reionization mid-points of $$z$$ ≤ 8 ($$z$$ ≤ 12), the shape differences are $${\lesssim}20{{\ \rm per\ cent}}$$ ($${\lesssim}40{{\ \rm per\ cent}}$$) because of a surprisingly well-matched cancellation between thermal broadening and pressure smoothing that occurs for realistic thermal histories. We also consider fluctuations in the ultraviolet background, finding their impact on the power spectrum to be much smaller than temperature fluctuations at k ≳ 0.01 s km−1. Furthermore, we compare our models to power spectrum measurements, finding that none of our models with reionization mid-points of $$z$$ < 8 is strongly preferred over another and that all of our models with mid-points of $$z$$ ≥ 8 are excluded at 2.5σ. Future measurements may be able to distinguish between viable reionization models if they can be performed at lower k or, alternatively, if the error bars on the high-k power can be reduced by a factor of 1.5. 
    more » « less
  2. ABSTRACT We present a new analysis of the rest-frame ultraviolet (UV) and optical spectra of a sample of three z > 8 galaxies discovered behind the gravitational lensing cluster RX J2129.4+0009. We combine these observations with z > 7.5 galaxies from the literature, for which similar measurements are available. As already pointed out in other studies, the high [O iii]λ5007/[O ii]λ3727 ratios (O32) and steep UV continuum slopes (β) are consistent with the values observed for low-redshift Lyman continuum emitters, suggesting that such galaxies contribute to the ionizing budget of the intergalactic medium. We construct a logistic regression model to estimate the probability of a galaxy being a Lyman continuum emitter based on the measured MUV, β, and O32. Using this probability and the UV luminosity function, we construct an empirical model that estimates the contribution of high-redshift galaxies to reionization. The preferred scenario in our analysis shows that at z ∼ 8, the average escape fraction of the galaxy population [i.e. including both LyC emitters (LCEs) and non-emitters] varies with MUV, with intermediate UV luminosity (−19 < MUV < −16) galaxies having larger escape fraction. Galaxies with faint UV luminosity (−16 < MUV < −13.5) contribute most of the ionizing photons. The relative contribution of faint versus bright galaxies depends on redshift, with the intermediate UV galaxies becoming more important over time. UV bright galaxies, although more likely to be LCEs at a given log(O32) and β, contribute the least of the total ionizing photon budget. 
    more » « less
  3. 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
  4. ABSTRACT New JWST observations are revealing the first galaxies to be prolific producers of ionizing photons, which we argue gives rise to a tension between different probes of reionization. Over the last two decades, a consensus has emerged where star-forming galaxies are able to generate enough photons to drive reionization, given reasonable values for their number densities, ionizing efficiencies $$\xi _{\rm ion}$$ (per unit ultraviolet luminosity), and escape fractions $$f_{\rm esc}$$. However, some new JWST observations infer high values of $$\xi _{\rm ion}$$ during reionization and an enhanced abundance of earlier ($$z\gtrsim 9$$) galaxies, dramatically increasing the number of ionizing photons produced at high z. Simultaneously, recent low-z studies predict significant escape fractions for faint reionization-era galaxies. Put together, we show that the galaxies we have directly observed ($$M_{\rm UV} < -15$$) not only can drive reionization, but would end it too early. That is, our current galaxy observations, taken at face value, imply an excess of ionizing photons and thus a process of reionization in tension with the cosmic microwave background and Lyman-$$\alpha$$ forest. Considering galaxies down to $$M_{\rm UV}\approx -11$$, below current observational limits, only worsens this tension. We discuss possible avenues to resolve this photon budget crisis, including systematics in either theory or observations. 
    more » « less
  5. ABSTRACT The nebular recombination line H α is widely used as a star formation rate (SFR) indicator in the local and high-redshift Universe. We present a detailed H α radiative transfer study of high-resolution isolated Milky-Way and Large Magellanic Cloud simulations that include radiative transfer, non-equilibrium thermochemistry, and dust evolution. We focus on the spatial morphology and temporal variability of the H α emission, and its connection to the underlying gas and star formation properties. The H α and H β radial and vertical surface brightness profiles are in excellent agreement with observations of nearby galaxies. We find that the fraction of H α emission from collisional excitation amounts to fcol ∼ 5–$$10{{\ \rm per\ cent}}$$, only weakly dependent on radius and vertical height, and that scattering boosts the H α luminosity by $$\sim 40{{\ \rm per\ cent}}$$. The dust correction via the Balmer decrement works well (intrinsic H α emission recoverable within 25 per cent), though the dust attenuation law depends on the amount of attenuation itself both on spatially resolved and integrated scales. Important for the understanding of the H α–SFR connection is the dust and helium absorption of ionizing radiation (Lyman continuum [LyC] photons), which are about $$f_{\rm abs}\approx 28{{\ \rm per\ cent}}$$ and $$f_{\rm He}\approx 9{{\ \rm per\ cent}}$$, respectively. Together with an escape fraction of $$f_{\rm esc}\approx 6{{\ \rm per\ cent}}$$, this reduces the available budget for hydrogen line emission by nearly half ($$f_{\rm H}\approx 57{{\ \rm per\ cent}}$$). We discuss the impact of the diffuse ionized gas, showing – among other things – that the extraplanar H α emission is powered by LyC photons escaping the disc. Future applications of this framework to cosmological (zoom-in) simulations will assist in the interpretation of spectroscopy of high-redshift galaxies with the upcoming James Webb Space Telescope. 
    more » « less