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.

Quantifying Lyman-$α$ emissions from reionization fronts

Imaging reionization's last phases with I-front Lyman-$α$ emissions

The Impact of Neutral Intergalactic Gas on Lyman-$α$ Intensity Mapping During Reionization

Scattering of Lyα Photons through the Reionizing Intergalactic Medium: I. Spectral Energy Distribution

First Light And Reionisation Epoch Simulations (FLARES) XIII: The Lyman-continuum emission of high-redshift galaxies

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

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

Dark-ages Reionization and Galaxy Formation Simulation XX. The Ly$α$ IGM transmission properties and environment of bright galaxies during the Epoch of Reionization

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

Extended red wings and the visibility of reionization-epoch Lyman-$α$ emitters

Reionization morphology and intrinsic velocity offsets allow transmission of Lyman-α emission from JADES-GS-z13-1-LA

The hydrodynamic response of small-scale structure to reionization drives large IGM temperature fluctuations that persist to z = 4

The Evolution of the Lyman-alpha Luminosity Function during Reionization

Can 21-Cm Observations Discriminate Between High-Mass and Low-Mass Galaxies As Reionization Sources?

Quasar Proximity Zones and Patchy Reionization

The reionizing bubble size distribution around galaxies

Inferring the IGM thermal history during reionisation with the Lyman-$α$ forest power spectrum at redshift $z \simeq 5$

Deciphering spatially resolved Lyman-alpha profiles in reionization analogs: the Sunburst Arc at cosmic noon

Non-equilibrium temperature evolution of ionization fronts during the Epoch of Reionization

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

Deciphering Lyman-α emission deep into the epoch of reionization