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.

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

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

Reionization relics in the cross-correlation between the Ly$α$ forest and 21 cm intensity mapping in the post-reionization era

Reionization Relics in the Cross-Correlation Between the Lyα Forest and 21 Cm Intensity Mapping in the Post-Reionization Era

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

Possible evidence for a large-scale enhancement in the Lyman-$α$ forest power spectrum at redshift $\mathbf{\textit{z}\geq 4}$

Constraining cosmology and ionization history with combined 21 cm power spectrum and global signal measurements

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

Separating the Memory of Reionisation from Cosmology in the Lyα Forest Power Spectrum at the Post-Reionisation Era

One-dimensional power spectrum from first DESI Lyman-α forest

Cosmological Constraints from the eBOSS Lyman-$α$ Forest using the PRIYA Simulations

Separating the memory of reionization from cosmology in the Ly$\alpha$ forest power spectrum at the post-reionization era

Fisher matrix forecasts on the astrophysics of galaxies during the epoch of reionisation from the 21-cm power spectra

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

Forecasting constraints on the high-z IGM thermal state from the Lyman-$α$ forest flux auto-correlation function

First Constraints on the Epoch of Reionization Using the non-Gaussianity of the Kinematic Sunyaev-Zel{'}dovich Effect from the South Pole Telescope and {\it Herschel}-SPIRE Observations

Light-Cone Effect on the Reionization 21-Cm Power Spectrum

Towards unveiling the Cosmic Reionization: the ionizing photon production efficiency ($ξ_{ion}$) of Low-mass H$α$ emitters at $z \sim 2.3$

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

Improved Constraints on the 21 cm EoR Power Spectrum and the X-Ray Heating of the IGM with HERA Phase I Observations

First Constraints on the Epoch of Reionization Using the Non-Gaussianity of the Kinematic Sunyaev-Zel'dovich Effect from the South Pole Telescope and Herschel-SPIRE Observations