An official website of the United States government
ABSTRACT Recent quasar absorption line observations suggest that reionization may end as late as $$z \approx 5.3$$. As a means to search for large neutral hydrogen islands at $$z\ \lt\ 6$$, we revisit long dark gaps in the Ly $$\beta$$ forest in Very Large Telescope/X-Shooter and Keck/Echellette Spectrograph and Imager quasar spectra. We stack the Ly $$\alpha$$ forest corresponding to both edges of these Ly $$\beta$$ dark gaps and identify a damping wing-like extended absorption profile. The average redshift of the stacked forest is $z=5.8$. By comparing these observations with reionization simulations, we infer that such a damping wing-like feature can be naturally explained if these gaps are at least partially created by neutral islands. Conversely, simulated dark gaps lacking neutral hydrogen struggle to replicate the observed damping wing features. Furthermore, this damping wing-like profile implies that the volume-averaged neutral hydrogen fraction must be $$\langle x_{\rm H\,{\small {I}}} \rangle \ge 6.1 \pm 3.9~{{\ \rm per\ cent}}$$ at $z = 5.8$. Our results offer robust evidence that reionization extends below $z=6$.
more »
« less
The effect of reionization on direct measurements of the mean free path
ABSTRACT Recent measurements of the ionizing photon mean free path (MFP) based on composite quasar spectra may point to reionization ending at z < 6. These measurements are challenging because they rely on assumptions about the proximity zones of the quasars. For example, some quasars might have been close to neutral patches where reionization was still ongoing (‘neutral islands’), and it is unclear how they would affect the measurements. We address this question with mock MFP measurements from radiative transfer simulations. We find that, even in the presence of neutral islands, our mock MFP measurements agree to within $$30~{{\ \rm per\ cent}}$$ with the true spatially averaged MFP in our simulations, which includes opacity from both the ionized medium and the islands. The inferred MFP is sensitive at the $$\lt ~50~{{\ \rm per\ cent}}$$ level to assumptions about quasar environments and lifetimes for realistic models. We demonstrate that future analyses with improved data may require explicitly modelling the effects of neutral islands on the composite spectra, and we outline a method for doing this. Lastly, we quantify the effects of neutral islands on Lyman-series transmission, which has been modelled with optically thin simulations in previous MFP analyses. Neutral islands can suppress transmission at λrest < 912 Å significantly, up to a factor of 2 for zqso = 6 in a plausible reionization scenario, owing to absorption by many closely spaced lines as quasar light redshifts into resonance. However, the suppression is almost entirely degenerate with the spectrum normalization and thus does not significantly bias the inferred MFP.
more »
« less
- PAR ID:
- 10506396
- Publisher / Repository:
- Oxford University Press
- Date Published:
- Journal Name:
- Monthly Notices of the Royal Astronomical Society
- Volume:
- 530
- Issue:
- 4
- ISSN:
- 0035-8711
- Format(s):
- Medium: X Size: p. 5209-5219
- Size(s):
- p. 5209-5219
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
ABSTRACT We introduce a probabilistic approach to select 6 ≤ $$z$$ ≤ 8 quasar candidates for spectroscopic follow-up, which is based on density estimation in the high-dimensional space inhabited by the optical and near-infrared photometry. Densities are modelled as Gaussian mixtures with principled accounting of errors using the extreme deconvolution (XD) technique, generalizing an approach successfully used to select lower redshift ($$z$$ ≤ 3) quasars. We train the probability density of contaminants on 1902 071 7-d flux measurements from the 1076 deg2 overlapping area from the Dark Energy Camera Legacy Survey (DECaLS) ($$z$$), VIKING (YJHKs), and unWISE (W1W2) imaging surveys, after requiring they dropout of DECaLS g and r, whereas the distribution of high-$$z$$ quasars are trained on synthetic model photometry. Extensive simulations based on these density distributions and current estimates of the quasar luminosity function indicate that this method achieves a completeness of $$\ge 56{{\ \rm per\ cent}}$$ and an efficiency of $$\ge 5{{\ \rm per\ cent}}$$ for selecting quasars at 6 < $$z$$ < 8 with JAB < 21.5. Among the classified sources are 8 known 6 < $$z$$ < 7 quasars, of which 2/8 are selected suggesting a completeness $$\simeq 25{{\ \rm per\ cent}}$$, whereas classifying the 6 known (JAB < 21.5) quasars at $$z$$ > 7 from the entire sky, we select 5/6 or a completeness of $$\simeq 80{{\ \rm per\ cent}}$$. The failure to select the majority of 6 < $$z$$ < 7 quasars arises because our quasar density model is based on an empirical quasar spectral energy distribution model that underestimates the scatter in the distribution of fluxes. This new approach to quasar selection paves the way for efficient spectroscopic follow-up of Euclid quasar candidates with ground-based telescopes and James Webb Space Telescope.more » « less
-
ABSTRACT Fluctuations in Lyman-α (Ly α) forest transmission towards high-z quasars are partially sourced from spatial fluctuations in the ultraviolet background, the level of which are set by the mean free path of ionizing photons (λmfp). The autocorrelation function of Ly α forest flux characterizes the strength and scale of transmission fluctuations and, as we show, is thus sensitive to λmfp. Recent measurements at z ∼ 6 suggest a rapid evolution of λmfp at z > 5.0 which would leave a signature in the evolution of the autocorrelation function. For this forecast, we model mock Ly α forest data with properties similar to the XQR-30 extended data set at 5.4 ≤ z ≤ 6.0. At each z, we investigate 100 mock data sets and an ideal case where mock data matches model values of the autocorrelation function. For ideal data with λmfp = 9.0 cMpc at z = 6.0, we recover $$\lambda _{\text{mfp}}=12^{+6}_{-3}$$ cMpc. This precision is comparable to direct measurements of λmfp from the stacking of quasar spectra beyond the Lyman limit. Hypothetical high-resolution data leads to a $$\sim 40~{{\ \rm per\ cent}}$$ reduction in the error bars over all z. The distribution of mock values of the autocorrelation function in this work is highly non-Gaussian for high-z, which should caution work with other statistics of the high-z Ly α forest against making this assumption. We use a rigorous statistical method to pass an inference test, however future work on non-Gaussian methods will enable higher precision measurements.more » « less
-
ABSTRACT 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 disfavours 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.more » « less
-
ABSTRACT Recent observations from the EIGER JWST program have measured for the first time the quasar–galaxy cross-correlation function at $$z\approx 6$$. The autocorrelation function of faint $$z\approx 6$$ quasars was also recently estimated. These measurements provide key insights into the properties of quasars and galaxies at high redshift and their relation with the host dark matter haloes. In this work, we interpret these data building upon an empirical quasar population model that has been applied successfully to quasar clustering and demographic measurements at $$z\approx 2\!-\!4$$. We use a new, large-volume N-body simulation with more than a trillion particles, FLAMINGO-10k, to model quasars and galaxies simultaneously. We successfully reproduce observations of $$z\approx 6$$ quasars and galaxies (i.e. their clustering properties and luminosity functions), and infer key quantities such as their luminosity–halo mass relation, the mass function of their host haloes, and their duty cycle/occupation fraction. Our key findings are (i) quasars reside on average in $$\approx 10^{12.5}\, {\rm M}_{\odot }$$ haloes (corresponding to $$\approx 5\sigma$$ fluctuations in the initial conditions of the linear density field), but the distribution of host halo masses is quite broad; (ii) the duty cycle of (UV-bright) quasar activity is relatively low ($$\approx 1~{{\ \rm per\ cent}}$$); (iii) galaxies (that are bright in [O iii]) live in much smaller haloes ($$\approx 10^{10.9}\, {\rm M}_{\odot }$$) and have a larger duty cycle (occupation fraction) of $$\approx 13~{{\ \rm per\ cent}}$$. Finally, we focus on the inferred properties of quasars and present a homogeneous analysis of their evolution with redshift. The picture that emerges reveals a strong evolution of the host halo mass and duty cycle of quasars at $$z\approx 2\!-\!6$$, and calls for new investigations of the role of quasar activity across cosmic time.more » « less
