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.

Cosmic Reionization On Computers: Biases and Uncertainties in the Measured Mean Free Path at the End Stage of Reionization

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

Robustness of direct measurements of the mean free path of ionizing photons in the epoch of reionization

The effect of reionization on direct measurements of the mean free path

Probing Ultra-late Reionization: Direct Measurements of the Mean Free Path over $5<z<6$

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

Quasar Proximity Zones and Patchy Reionization

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

The Predicament of Absorption-dominated Reionization II: Observational Estimate of the Clumping Factor at the End of Reionization

Cosmic Reionization On Computers: Statistics, Physical Properties and Environment of Lyman Limit Systems at $z\sim6$

Cosmic Reionization on Computers: Statistics, Physical Properties, and Environments of Lyman Limit Systems at z ∼ 6

Implications of $z \sim 6$ Quasar Proximity Zones for the Epoch of Reionization and Quasar Lifetimes

Measuring the Mean Free Path of HI Ionizing Photons at 3.2≤ Z≤4.6 with DESI Y1 Quasars

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

Measuring the Mean Free Path of HI Ionizing Photons at $3.2\leq z\leq4.6$ with DESI Y1 Quasars

Precisely Measuring the Cosmic Reionization History from IGM Damping Wings Towards Quasars

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

Modeling Quasar Proximity Zones in a Realistic Cosmological Environment with a Self-consistent Light Curve

On the rise and fall of galactic ionizing output at the end of reionization