The mean free path of ionizing photons, λmfp, is a key factor in the photoionization of the intergalactic medium (IGM). At z ≳ 5, however, λmfp may be short enough that measurements towards QSOs are biased by the QSO proximity effect. We present new direct measurements of λmfp that address this bias and extend up to z ∼ 6 for the first time. Our measurements at z ∼ 5 are based on data from the Giant Gemini GMOS survey and new Keck LRIS observations of low-luminosity QSOs. At z ∼ 6 we use QSO spectra from Keck ESI and VLT X-Shooter. We measure $\lambda _{\rm mfp} = 9.09^{+1.62}_{-1.28}$ proper Mpc and $0.75^{+0.65}_{-0.45}$ proper Mpc (68 per cent confidence) at z = 5.1 and 6.0, respectively. The results at z = 5.1 are consistent with existing measurements, suggesting that bias from the proximity effect is minor at this redshift. At z = 6.0, however, we find that neglecting the proximity effect biases the result high by a factor of two or more. Our measurement at z = 6.0 falls well below extrapolations from lower redshifts, indicating rapid evolution in λmfp over 5 < z < 6. This evolution disfavors models in which reionization ended early enough that the IGM had time to fully relax hydrodynamically by z = 6, but is qualitatively consistent with models wherein reionization completed at z = 6 or even significantly later. Our mean free path results are most consistent with late reionization models wherein the IGM is still 20 per cent neutral at z = 6, although our measurement at z = 6.0 is even lower than these models prefer.

We present a new investigation of the intergalactic medium (IGM) near the end of reionization using “dark gaps” in the Lyα forest. Using spectra of 55 QSOs at z em > 5.5, including new data from the XQR-30 VLT Large Programme, we identify gaps in the Lyα forest where the transmission averaged over 1 comoving h −1 Mpc bins falls below 5%. Nine ultralong (L > 80 h −1 Mpc) dark gaps are identified at z < 6. In addition, we quantify the fraction of QSO spectra exhibiting gaps longer than 30 h −1 Mpc, F 30, as a function of redshift. We measure F 30 ≃ 0.9, 0.6, and 0.15 at z = 6.0, 5.8, and 5.6, respectively, with the last of these long dark gaps persisting down to z ≃5.3. Comparing our results with predictions from hydrodynamical simulations, we find that the data are consistent with models wherein reionization extends significantly below redshift six. Models wherein the IGM is essentially fully reionized that retain large-scale fluctuations in the ionizing UV background at z ≲6 are also potentially consistent with the data. Overall, our results suggest that signatures of reionization in the form of islands of neutral hydrogen and/or large-scale fluctuations in the ionizing background remain present in the IGM until at least z ≃ 5.3.

We present a new investigation of the intergalactic medium near reionization using dark gaps in the Lyβ forest. With its lower optical depth, Lyβ offers a potentially more sensitive probe to any remaining neutral gas compared to the commonly used Lyα line. We identify dark gaps in the Lyβ forest using spectra of 42 QSOs at z em > 5.5, including new data from the XQR-30 VLT Large Programme. Approximately 40% of these QSO spectra exhibit dark gaps longer than 10 h −1 Mpc at z ≃ 5.8. By comparing the results to predictions from simulations, we find that the data are broadly consistent both with models where fluctuations in the Lyα forest are caused solely by ionizing ultraviolet background fluctuations and with models that include large neutral hydrogen patches at z < 6 due to a late end to reionization. Of particular interest is a very long (L = 28 h −1 Mpc) and dark (τ eff ≳ 6) gap persisting down to z ≃ 5.5 in the Lyβ forest of the z = 5.85 QSO PSO J025−11. This gap may support late reionization models with a volume-weighted average neutral hydrogen fraction of 〈x H I〉 ≳ 5% by z = 5.6. Finally, we infer constraints on 〈x H I〉 over 5.5 ≲ z ≲ 6.0 based on the observed Lyβ dark gap length distribution and a conservative relationship between gap length and neutral fraction derived from simulations. We find 〈x H I〉 ≤ 0.05, 0.17, and 0.29 at z ≃ 5.55, 5.75, and 5.95, respectively. These constraints are consistent with models where reionization ends significantly later than z = 6.

The observed large-scale scatter in Lyα opacity of the intergalactic medium at z < 6 implies large fluctuations in the neutral hydrogen fraction that are unexpected long after reionization has ended. A number of models have emerged to explain these fluctuations that make testable predictions for the relationship between Lyα opacity and density. We present selections of z = 5.7 Lyα-emitting galaxies (LAEs) in the fields surrounding two highly opaque quasar sightlines with long Lyα troughs. The fields lie toward the z = 6.0 quasar ULAS J0148+0600, for which we reanalyze previously published results using improved photometric selection, and toward the z = 6.15 quasar SDSS J1250+3130, for which results are presented here for the first time. In both fields, we report a deficit of LAEs within 20 h −1 Mpc of the quasar. The association of highly opaque sightlines with galaxy underdensities in these two fields is consistent with models in which the scatter in Lyα opacity is driven by large-scale fluctuations in the ionizing UV background or by an ultra-late reionization that has not yet concluded at z = 5.7.

We present Keck/OSIRIS infrared IFU observations of the z = 3.153 sub-DLA DLA2233+131, previously detected in absorption to a background quasar and studied with single slit spectroscopy and PMAS integral field spectroscopy (IFU). We used the Laser Guide Star Adaptive Optics (LGSAO) and OSIRIS IFU to reduce the point-spread function of the background quasar to FWHM∼0.15 and marginally resolve extended, foreground DLA emission. We detect [O III]λ5007 emission with a flux F [O III]λ5007 = (2.4 ± 0.5) × 10 −17 ergs s −1 cm −2 , as well as unresolved [O III]λ4959 and Hβλ4861 emission. Using a composite spectrum over the emission region, we measure dynamical mass ∼ 3.1×10 9 M . We make several estimates of star formation rate using [O III]λ5007 and Hβλ4861 emission, and measure a star formation rate of ∼ 7.1 − 13.6 M yr −1 . We map [O III]λ5007 and Hβλ4861 emission and the corresponding velocity fields to search for signs of kinematic structure. These maps allow for a more detailed kinematic analysis than previously possible for this galaxy. While some regions show slightly red and blue-shifted emission indicative of potential edge-on disk rotation, the data are insufficient to support this interpretation.