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.

Damping wings in the Lyman-α forest: a model-independent measurement of the neutral fraction at 5.4<z<6.1

Damping Wing-Like Features in the Stacked Ly$α$ Forest: Potential Neutral Hydrogen Islands at $z<6$

IGM damping wing constraints on the tail end of reionisation from the enlarged XQR-30 sample

Long Dark Gaps in the Ly$β$ Forest at $z<6$: Evidence of Ultra Late Reionization from XQR-30 Spectra

Damping wing absorption associated with a giant Ly α trough at z < 6: direct evidence for late-ending reionization

Damping wing absorption associated with a giant Ly$α$ trough at $z < 6$: direct evidence for late-ending reionization

The Universe is at Most 88% Neutral at z=10.6

Chronicling the reionization history at $6\lesssim z \lesssim 7$ with emergent quasar damping wings

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

Reionization of the Intergalactic Medium and the Damping Wing of the Gunn-Peterson Trough

JWST Measurements of Neutral Hydrogen Fractions and Ionized Bubble Sizes at $z=7-12$ Obtained with Ly$α$ Damping Wing Absorptions in 27 Bright Continuum Galaxies

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

The Characteristic Shape of Damping Wings During Reionization

Ly$α$ Emission at $z=7-13$: Clear Ly$α$ Equivalent Width Evolution Indicating the Late Cosmic Reionization History

Possible Contamination of the Intergalactic Medium Damping Wing in ULAS J1342+0928 by Proximate Damped Ly$α$ Absorption

A Significantly Neutral Intergalactic Medium Around the Luminous Z = 7 Quasar J0252–0503

Neutral Fraction of Hydrogen in the Intergalactic Medium Surrounding High-Redshift Gamma-Ray Burst 210905A

Measuring the photo-ionization rate, neutral fraction and mean free path of HI ionizing photons at $4.9 \leq z \leq 6.0$ from a large sample of XShooter and ESI spectra

A massive, neutral gas reservoir permeating a galaxy proto-cluster after the reionization era

Constraining Reionization with Lyα Damping-Wing Absorption in Galaxy Spectra: A Machine Learning Model Based on Reionization Simulations