Abstract:The Lyman alpha (Lya) forest in the spectra of z>5 quasars provides a powerful probe of the late stages of the Epoch of Reionization (EoR). With the recent advent of exquisite datasets such as XQR-30, many models have struggled to reproduce the observed large-scale fluctuations in the Lya opacity. Here we introduce a Bayesian analysis framework that forward-models large-scale lightcones of IGM properties, and accounts for unresolved sub-structure in the Lya opacity by calibrating to higher-resolution hydrodynamic simulations. Our models directly connect physically-intuitive galaxy properties with the corresponding IGM evolution, without having to tune "effective" parameters or calibrate out the mean transmission. The forest data, in combination with UV luminosity functions and the CMB optical depth, are able to constrain global IGM properties at percent level precision in our fiducial model. Unlike many other works, we recover the forest observations without evoking a rapid drop in the ionizing emissivity from z~7 to 5.5, which we attribute to our sub-grid model for recombinations. In this fiducial model, reionization ends at $z=5.44\pm0.02$ and the EoR mid-point is at $z=7.7\pm0.1$. The ionizing escape fraction increases towards faint galaxies, showing a mild redshift evolution at fixed UV magnitude, Muv. Half of the ionizing photons are provided by galaxies fainter than Muv~-12, well below direct detection limits of optical/NIR instruments including JWST. We also show results from an alternative galaxy model that does not allow for a redshift evolution in the ionizing escape fraction. Despite being decisively disfavored by the Bayesian evidence, the posterior of this model is in qualitative agreement with that from our fiducial model. We caution however that our conclusions regarding the early stages of the EoR and which sources reionized the Universe are more model-dependent.

The forest at EndEoR: The effect of Lyman Limit Systems on the End of Reionisation

The Lyman-limit photon mean free path at the end of late reionisation in the Sherwood-Relics simulations

The 21-cm signal during the end stages of reionization

The impact of varying inhomogeneous reionization histories on metrics of Ly$α$ opacity

The Lyman-limit photon mean free path at the end of late reionization in the Sherwood–Relics simulations

Percent-level timing of reionization: self-consistent, implicit-likelihood inference from XQR-30+ Ly$α$ forest data

Constraints on the Evolution of the Ionizing Background and Ionizing Photon Mean Free Path at the End of Reionization

The Impact of Cosmic Variance on Inferences of Global Neutral Fraction Derived from Ly$α$ Luminosity Functions During Reionization

Convergence of small scale Ly$α$ structure at high-$z$ under different reionization scenarios

Efficient Reionization in a Large Hydrodynamic Galaxy Formation Simulation

Lyα Forest Power Spectrum As an Emerging Window into the Epoch of Reionization and Cosmic Dawn

The effect of helium reionization on the Ly α forest hydrogen flux statistics

Extracting the Astrophysics of Reionization from the Lyα Forest Power Spectrum: a First Forecast

Large scale opacity fluctuations in the Lyman alpha forest: evidence for QSOs dominating the ionising UV background at z ~ 5.5-6 ?

The Global History of Reionisation

The effect of helium reionization on the Lyman-alpha forest hydrogen flux statistics

Cosmic Reionization on Computers: The Evolution of Ionizing Background and Mean Free Path

Imprints of temperature fluctuations on the $z\sim5$ Lyman-$α$ forest: a view from radiation-hydrodynamic simulations of reionization

On the correlation between Ly$α$ forest opacity and galaxy density in late reionization models