Studying the Lyman α optical depth fluctuations at z ∼ 5.5 using fast semi-numerical methods

We present a computationally efficient and fast semi-numerical technique for simulating the Lyman α (Ly α) absorption optical depth in presence of neutral hydrogen ‘islands’ left over from reionization at redshifts 5 ≲ z ≲ 6. The main inputs to the analysis are (i) a semi-numerical photon-conserving model of ionized regions during reionization (named script) along with a prescription for simulating the shadowing by neutral islands and (ii) the fluctuating Gunn–Peterson approximation to model the Ly α absorption. Our model is then used for simulating the large-scale fluctuations in the effective optical depth as observed along sightlines towards high-z quasars. Our model is fully described by five parameters. By setting two of them to default values and varying the other three, we obtain the constraints on reionization history at 5 ≲ z ≲ 6 as allowed by the data. We confirm that reionization is not complete before z ∼ 5.6 at ≳2σ confidence, with the exact confidence limits depending on how the non-detections of the flux in the data are treated. We also confirm that the completion of reionization can be as late as z ∼ 5.2. With further improvements in the model and with more sightlines at z ∼ 6, we can take advantage of the computational efficiency of our analysis to obtain more stringent constraints on the ionization fraction at the tail end of reionization

A comparison of quasar emission reconstruction techniques for z ≥ 5.0 Lyman α and Lyman β transmission

Reconstruction techniques for intrinsic quasar continua are crucial for the precision study of Lyman-$\alpha$ (Ly-$\alpha$) and Lyman-$\beta$ (Ly-$\beta$) transmission at $z>5.0$, where the $\lambda5.7$ with IR X-Shooter spectroscopy, obtaining well-characterised measurements for the mean flux transmission at $4.7<z<6.3$. Our results demonstrate the importance of testing and, when relevant, training, continuum reconstruction techniques in a systematic way

A tale of seven narrow spikes and a long trough: Constraining the timing of the percolation of HII bubbles at the tail end of reionization with ULAS J1120+0641

High-signal to noise observations of the Ly$\alpha$ forest transmissivity in the z = 7.085 QSO ULAS J1120+0641 show seven narrow transmission spikes followed by a long 240 cMpc/h trough. Here we use radiative transfer simulations of cosmic reionization previously calibrated to match a wider range of Ly$\alpha$ forest data to show that the occurrence of seven transmission spikes in the narrow redshift range z = 5.85 - 6.1 is very sensitive to the exact timing of reionization. Occurrence of the spikes requires the most under dense regions of the IGM to be already fully ionised. The rapid onset of a long trough at z = 6.12 requires a strong decrease of the photo-ionisation rate at z$\sim$6.1 in this line-of-sight, consistent with the end of percolation at this redshift. The narrow range of reionisation histories that we previously found to be consistent with a wider range of Ly$\alpha$ forest data have a reasonable probability of showing seven spikes and the mock absorption spectra provide an excellent match to the spikes and the trough in the observed spectrum of ULAS J1120+0641. Despite the large overall opacity of Ly$\alpha$ at z > 5.8, larger samples of high signal-to-noise observations of rare transmission spikes should therefore provide important further insights into the exact timing of the percolation of HII bubbles at the tail-end of reionizatio

Re-examining the case for neutral gas near the redshift 7 quasar ULAS J1120+0641

Signs of damping-wing absorption attenuating the Lyman α emission line of the first known z ∼ 7 quasar, ULAS J1120+0641, recently provided exciting evidence of a significantly neutral intergalactic medium (IGM). This long-awaited signature of reionization was inferred, in part, from a deficit of flux in the quasar's Lyman α emission line based on predictions from a composite of lower redshift quasars. The composite sample was chosen based on its C IV emission line properties; however, as the original study by Mortlock et al. noted, the composite contained a slight velocity offset in C IV compared to ULAS J1120+0641. Here we test whether this offset may be related to the predicted strength of the Lyman α emission line. We confirm the significant (∼10 per cent at rms) scatter in Lyman α flux for quasars of a given C IV velocity and equivalent width found by Mortlock et al. We further find that among lower redshift objects chosen to more closely match the C IV properties of ULAS J1120+0641, its Lyman α emission falls within the observed distribution of fluxes. Among lower redshift quasars chosen to more closely match in C IV velocity and equivalent width, we find that ULAS J1120+0641 falls within the observed distribution of Lyman α emission line strengths. This suggests that damping-wing absorption may not be present, potentially weakening the case for neutral gas around this object. Larger samples of z > 7 quasars may therefore be needed to establish a clearer picture of the IGM neutral fraction at these redshifts

Three Lyα Emitting Galaxies within a Quasar Proximity Zone at z ~ 5.8

Quasar proximity zones at z > 5.5 correspond to overdense and overionized environments. Galaxies found inside proximity zones can therefore display features that would otherwise be masked by absorption in the intergalactic medium. We demonstrate the utility of this quasar-galaxy synergy by reporting the discovery of the first three “proximate Lyα emitters” (LAEs) within the proximity zone of quasar J0836+0054 at z = 5.795 (Aerith A, B, and C). Aerith A, located behind the quasar with an impact parameter D^{\perp} = 278 \pm 8 pkpc, provides the first detection of an Lyα transverse proximity effect. We model the transmission and show that it constrains the onset of J0836ʼs quasar phase to 0.2Myr < 28Myr < t in the past. The second object, Aerith B at a distance D < 912 pkpc from the quasar, displays a bright and broad double-peaked Lyα emission line. The peak separation implies a low ionizing f_{esc} \leqslant 1%. We fit the Lyα line with an outflowing shell model, finding a typical central density log N_{HI}/cm^{-2} = 19.3_{-0.2}^{+0.8}, outflow velocity v_{out} = 16_{-11}^{+4}km s^{-1}, and gas temperature log T/K = 3.8_{-0.7}^{+0.8} compared to 2 < z < 3 analog LAEs. We detect object Aerith C via an Lyα emission line at z = 5.726. This corresponds with the edge of the quasar’s proximity zone (Dz < 0.02), suggesting that the proximity zone is truncated by a density fluctuation. Via the analyses conducted here, we illustrate how proximate LAEs offer unique insight into the ionizing properties of both quasars and galaxies during hydrogen reionization

Large Ly α opacity fluctuations and low CMB τ in models of late reionization with large islands of neutral hydrogen extending to z &lt; 5.5

High-redshift QSO spectra show large spatial fluctuations in the Ly-alpha opacity of the intergalactic medium on surprisingly large scales at z>~5.5. We present a radiative transfer simulation of cosmic reionization driven by galaxies that reproduces this large scatter and the rapid evolution of the Ly-alpha opacity distribution at 5<z<6. The simulation also reproduces the low Thomson scattering optical depth reported by the latest CMB measurement and is consistent with the observed short near-zones and strong red damping wings in the highest-redshift QSOs. It also matches the rapid disappearance of observed Ly-alpha emission by galaxies at z>~6. Reionization is complete at z=5.3 in our model, and 50% of the volume of the Universe is ionized at z=7. Agreement with the Ly-alpha forest data in such a late reionization model requires a rapid evolution of the ionizing emissivity of galaxies that peaks at z~6.8. The late end of reionization results in a large scatter in the photoionisation rate and the neutral hydrogen fraction at redshifts as low as z<~5.5 with large residual neutral 'islands' that can produce very long Gunn-Peterson troughs resembling those seen in the data

The role of galaxies and AGNs in reionizing the IGM – II. Metal-tracing the faint sources of reionization at 5 ≲ z ≲ 6

We present a new method to study the contribution of faint sources to the ultraviolet background using the 1D correlation of metal absorbers with the intergalactic medium transmission in a quasi-stellar object (QSO) sightline. We take advantage of a sample of 25 high signal-to-noise ratio QSO spectra to retrieve 150 triply-ionized carbon (C IV) absorbers at 4.5 ≲ z ≲ 6.2, of which 37 systems whose expected H I absorption lie in the Lyman-α forest. We derive improved constraints on the cosmic density of C IV at 4.3 < z < 6.2 and infer from abundance matching that C IV absorbers trace MUV ≲ −16 galaxies. Correlation with the Lyman-α forest of the QSOs indicates that these objects are surrounded by a highly opaque region at r ≲ 5 cMpc h−1 followed by an excess of transmission at r ≳ 10 cMpc h−1 detected at 2.7σ. This is in contrast to equivalent measurements at lower redshifts where only the opaque trough is detected. We interpret this excess as a statistical enhancement of the local photoionization rate due to clustered faint galaxies around the C IV absorbers. Using the analytical framework described in Paper I of this series, we derive a constraint on the average product of the escape fraction and the Lyman continuum photon production efficiency of the galaxy population clustered around the C IV absorbers, log⟨fescξion⟩/[erg−1Hz]=25.01+0.30−0.19⁠. This implies that faint galaxies beyond the reach of current facilities may have harder radiation fields and/or larger escape fractions than currently detected objects at the end of the reionization epoch

Fluorescent C II* 1335A emission spectroscopically resolved in a galaxy at z = 5.754

We report the discovery of the first spectroscopically resolved C II /C II* 1334, 1335A doublet in the Lyman-break galaxy J0215-0555 at z = 5.754. The separation of the resonant and fluorescent emission channels was possible thanks to the large redshift of the source and long integration time, as well as the small velocity width of the feature, 0.6 +- 0.2A. We model this emission and find that at least two components are required to reproduce the combination of morphologies of C II* emission, C II absorption and emission, and Lyman-alpha emission from the object. We suggest that the close alignment between the fluorescence and Lyman-alpha emission could indicate an ionisation escape channel within the object. While the faintness of such a C II /C II* doublet makes it prohibitively difficult to pursue for similar systems with current facilities, we suggest it can become a valuable porosity diagnostic in the era of JWST and the upcoming generations of ELTs

The role of galaxies and AGN in reionising the IGM -- III : IGM-galaxy cross-correlations at z~6 from 8 quasar fields with DEIMOS and MUSE

We present improved results of the measurement of the correlation between galaxies and the intergalactic medium transmission at the end of reionization. We have gathered a sample of 13 spectroscopically confirmed Lyman-break galaxies (LBGs) and 21 Lyman-α emitters (LAEs) at angular separations 20 arcsec ≲ θ ≲ 10 arcmin (∼0.1–4 pMpc at z ∼ 6) from the sightlines to eight background z ≳ 6 quasars. We report for the first time the detection of an excess of Lyman-α transmission spikes at ∼10–60 cMpc from LAEs (3.2σ) and LBGs (1.9σ). We interpret the data with an improved model of the galaxy–Lyman-α transmission and two-point cross-correlations, which includes the enhanced photoionization due to clustered faint sources, enhanced gas densities around the central bright objects and spatial variations of the mean free path. The observed LAE(LBG)–Lyman-α transmission spike two-point cross-correlation function (2PCCF) constrains the luminosity-averaged escape fraction of all galaxies contributing to reionization to ⟨fesc⟩MUV−20(2σ)⁠) is necessary to reproduce the observed 2PCCF and that reionization might be driven by different sub-populations around LBGs and LAEs at z ∼ 6