Showing 1–50 of 98 results for author: Becker, G D

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

Authors: Bayu Wilson, Anson D'Aloisio, George D. Becker, Christopher Cain, Eli Visbal

Abstract: Long troughs observed in the $z > 5.5$ Ly$α$ and Ly$β$ forests are thought to be caused by the last remaining neutral patches during the end phases of reionization -- termed neutral islands. If this is true, then the longest troughs mark locations where we are most likely to observe the reionizing intergalactic medium (IGM). A key feature of the neutral islands is that they are bounded by ionizati… ▽ More Long troughs observed in the $z > 5.5$ Ly$α$ and Ly$β$ forests are thought to be caused by the last remaining neutral patches during the end phases of reionization -- termed neutral islands. If this is true, then the longest troughs mark locations where we are most likely to observe the reionizing intergalactic medium (IGM). A key feature of the neutral islands is that they are bounded by ionization fronts (I-fronts) which emit Lyman series lines. In this paper, we explore the possibility of directly imaging the outline of neutral islands with a narrowband survey targeting Ly$α$. In a companion paper, we quantified the intensity of I-front Ly$α$ emissions during reionization and its dependence on the spectrum of incident ionizing radiation and I-front speed. Here we apply those results to reionization simulations to model the emissions from neutral islands. We find that neutral islands would appear as diffuse structures that are tens of comoving Mpc across, with surface brightnesses in the range $\approx 1 - 5\times 10^{-21}$ erg s$^{-1}$ cm$^{-2}$ arcsec$^{-2}$. The islands are brighter if the spectrum of ionizing radiation driving the I-fronts is harder, and/or if the I-fronts are moving faster. We develop mock observations for current and futuristic observatories and find that, while extremely challenging, detecting neutral islands is potentially within reach of an ambitious observing program with wide-field narrowband imaging. Our results demonstrate the potentially high impact of low-surface brightness observations for studying reionization. △ Less

Submitted 20 June, 2024; originally announced June 2024.

Comments: 24 pages, 10 figures

Quantifying Lyman-$α$ emissions from reionization fronts

Authors: Bayu Wilson, Anson D'Aloisio, George D. Becker, Christopher Cain, Eli Visbal

Abstract: During reionization, intergalactic ionization fronts (I-fronts) are sources of Ly$α$ line radiation produced by collisional excitation of hydrogen atoms within the fronts. In principle, detecting this emission could provide direct evidence for a reionizing intergalactic medium (IGM). In this paper, we use a suite of high-resolution one-dimensional radiative transfer simulations run on cosmological… ▽ More During reionization, intergalactic ionization fronts (I-fronts) are sources of Ly$α$ line radiation produced by collisional excitation of hydrogen atoms within the fronts. In principle, detecting this emission could provide direct evidence for a reionizing intergalactic medium (IGM). In this paper, we use a suite of high-resolution one-dimensional radiative transfer simulations run on cosmological density fields to quantify the parameter space of I-front Ly$α$ emission. We find that the Ly$α$ production efficiency -- the ratio of emitted Ly$α$ flux to incident ionizing flux driving the front -- depends mainly on the I-front speed and the spectral index of the ionizing radiation. IGM density fluctuations on scales smaller than the typical I-front width produce scatter in the efficiency, but they do not significantly boost its mean value. The Ly$α$ flux emitted by an I-front is largest if 3 conditions are met simultaneously: (1) the incident ionizing flux is large; (2) the incident spectrum is hard, consisting of more energetic photons; (3) the I-front is traveling through a cosmological over-density, which causes it to propagate more slowly. We present a convenient parameterization of the efficiency in terms of I-front speed and incident spectral index. We make these results publicly available as an interpolation table and we provide a simple fitting function for a representative ionizing background spectrum. Our results can be applied as a sub-grid model for I-front Ly$α$ emissions in reionization simulations with spatial and/or temporal resolutions too coarse to resolve I-front structure. In a companion paper, we use our results to explore the possibility of directly imaging Ly$α$ emission around neutral islands during the last phases of reionization. △ Less

Submitted 20 June, 2024; originally announced June 2024.

Comments: 21 pages, 6 figures

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

Authors: Yongda Zhu, George D. Becker, Sarah E. I. Bosman, Christopher Cain, Laura C. Keating, Fahad Nasir, Valentina D'Odorico, Eduardo Bañados, Fuyan Bian, Manuela Bischetti, James S. Bolton, Huanqing Chen, Anson D'Aloisio, Frederick B. Davies, Rebecca L. Davies, Anna-Christina Eilers, Xiaohui Fan, Prakash Gaikwad, Bradley Greig, Martin G. Haehnelt, Girish Kulkarni, Samuel Lai, Ewald Puchwein, Yuxiang Qin, Emma V. Ryan-Weber , et al. (6 additional authors not shown)

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<6$, we revisit long dark gaps in the Ly$β$ forest in VLT/X-Shooter and Keck/ESI quasar spectra. We stack the Ly$α$ forest corresponding to both edges of these Ly$β$ dark gaps and identify a damping wing-like extended absorption profi… ▽ More 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<6$, we revisit long dark gaps in the Ly$β$ forest in VLT/X-Shooter and Keck/ESI quasar spectra. We stack the Ly$α$ forest corresponding to both edges of these Ly$β$ 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 HI} \rangle \geq 6.1 \pm 3.9\%$ at $z = 5.8$. Our results offer robust evidence that reionization extends below $z=6$. △ Less

Submitted 28 June, 2024; v1 submitted 20 May, 2024; originally announced May 2024.

Comments: 8 pages, 5 figures, 1 table; accepted for publication in MNRAS Letters

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

Authors: George D. Becker, James S. Bolton, Yongda Zhu, Seyedazim Hashemi

Abstract: Multiple observations now suggest that the hydrogen reionization may have ended well below redshift six. While there has previously been no conclusive proof of extended neutral islands in the $z < 6$ intergalactic medium, it is possible that such islands give rise to the giant Ly$α$ absorption troughs seen in the spectra of high-redshift quasars. Here we present evidence that the deepest and longe… ▽ More Multiple observations now suggest that the hydrogen reionization may have ended well below redshift six. While there has previously been no conclusive proof of extended neutral islands in the $z < 6$ intergalactic medium, it is possible that such islands give rise to the giant Ly$α$ absorption troughs seen in the spectra of high-redshift quasars. Here we present evidence that the deepest and longest-known Ly$α$ trough at $z < 6$, towards ULAS J0148+0600 (J0148), is associated with damping wing absorption. The evidence comes from a window of strong Ly$α$ transmission at the edge of the J0148 proximity zone. We show that the relatively smooth profile of this transmission window is highly unlikely to arise from resonant absorption alone, but is consistent with the presence of a damping wing. We further argue that the damping wing is unlikely to arise from a compact source due to the lack of associated metal lines, and is more likely to arise from an extended neutral island associated with the giant Ly$α$ trough. We investigate the physical conditions that may give rise to the strong transmission window, and speculate that it may signal an usually deep void, nearby ionizing sources, and/or the recent passage of an ionization front. △ Less

Submitted 14 May, 2024; originally announced May 2024.

Comments: 17 pages, 18 figures; resubmitted to MNRAS after addressing referee comments

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

Authors: Bradley Greig, Andrei Mesinger, Eduardo Bañados, George D. Becker, Sarah E. I. Bosman, Huanqing Chen, Frederick B. Davies, Valentina D'Odorico, Anna-Christina Eilers, Simona Gallerani, Martin G. Haehnelt, Laura Keating, Samuel Lai, Yuxiang Qin, Emma Ryan-Weber, Sindhu Satyavolu, Feige Wang, Jinyi Yang, Yongda Zhu

Abstract: The attenuation of Ly$α$ photons by neutral hydrogen in the intergalactic medium (IGM) at $z\gtrsim5$ continues to be a powerful probe for studying the epoch of reionisation. Given a framework to estimate the intrinsic (true) Ly$α$ emission of high-$z$ sources, one can infer the ionisation state of the IGM during reionisation. In this work, we use the enlarged XQR-30 sample of 42 high-resolution a… ▽ More The attenuation of Ly$α$ photons by neutral hydrogen in the intergalactic medium (IGM) at $z\gtrsim5$ continues to be a powerful probe for studying the epoch of reionisation. Given a framework to estimate the intrinsic (true) Ly$α$ emission of high-$z$ sources, one can infer the ionisation state of the IGM during reionisation. In this work, we use the enlarged XQR-30 sample of 42 high-resolution and high-SNR QSO spectra between $5.8\lesssim\,z\lesssim\,6.6$ obtained with VLT/X-Shooter to place constraints on the IGM neutral fraction. This is achieved using our existing Bayesian QSO reconstruction framework which accounts for uncertainties such as the: (i) posterior distribution of predicted intrinsic Ly$α$ emission profiles (obtained via covariance matrix reconstruction of the Ly$α$ and N V emission lines from unattenuated high-ionisation emission line profiles; C IV, Si IV + O IV] and C III]) and (ii) distribution of ionised regions within the IGM using synthetic damping wing profiles drawn from a $1.6^3$ Gpc$^3$ reionisation simulation. Following careful quality control, we used 23 of the 42 available QSOs to obtain constraints/limits on the IGM neutral fraction during the tail-end of reionisation. Our median and 68th percentile constraints on the IGM neutral fraction are: $0.20\substack{+0.14\\-0.12}$ and $0.29\substack{+0.14\\-0.13}$ at $z = 6.15$~and 6.35. Further, we also report 68th percentile upper-limits of $\bar{x}_{\mathrm{H\,{\scriptscriptstyle I}}} < 0.21$, 0.20, 0.21 and 0.18 at $z = 5.8, 5.95, 6.05$~and 6.55. These results imply reionisation is still ongoing at $5.8\lesssim\,z\lesssim\,6.55$, consistent with previous results from XQR-30 (dark fraction and Ly$α$ forest) along with other observational probes considered in the literature. △ Less

Submitted 18 April, 2024; originally announced April 2024.

Comments: 21 pages, 9 figures and 2 tables. Accepted for publication in MNRAS

Evidence of Pop~III stars' chemical signature in neutral gas at z~6. A study based on the E-XQR-30 spectroscopic sample

Authors: Alessio Sodini, Valentina D'Odorico, Stefania Salvadori, Irene Vanni, Manuela Bischetti, Guido Cupani, Rebecca Davies, George D. Becker, Eduardo Bañados, Sarah Bosman, Frederick Davies, Emanuele Paolo Farina, Andrea Ferrara, Laura Keating, Girish Kulkarni, Samuel Lai, Emma Ryan-Weber, Alma Maria Sebastian, Fabian Walter

Abstract: This study explores the metal enrichment signatures attributed to the first generation of stars (PopIII) in the Universe, focusing on the E-XQR-30 sample. We aim to identify traces of Pop III metal enrichment by analyzing neutral gas in the interstellar medium of primordial galaxies and their satellite clumps, detected in absorption. To chase the chemical signature of PopIII stars, we studied meta… ▽ More This study explores the metal enrichment signatures attributed to the first generation of stars (PopIII) in the Universe, focusing on the E-XQR-30 sample. We aim to identify traces of Pop III metal enrichment by analyzing neutral gas in the interstellar medium of primordial galaxies and their satellite clumps, detected in absorption. To chase the chemical signature of PopIII stars, we studied metal absorption systems in the E-XQR-30 sample, selected through the detection of the OI absorption line at 1302A. The OI line is a reliable tracer of HI and allowed us to overcome the challenges posed by the Lyman-$α$ forest's increasing saturation at redshifts above $\sim5$ to identify Damped Lyman-$α$ systems (DLA). We detected and analyzed 29 OI systems at $z\geq5.4$, differentiating between proximate DLAs (PDLA) and intervening DLAs. Voigt function fits were applied to obtain ionic column densities, and relative chemical abundances were determined for 28 systems. These were then compared with the predictions of theoretical models. Our findings expand the study of OI systems at $z\geq5.4$ fourfold. No systematic differences were observed in the average chemical abundances between PDLAs and intervening DLAs. The chemical abundances in our sample align with literature systems at $z>4.5$, suggesting a similar enrichment pattern for this class of absorption systems. A comparison between these DLA-analogues at $4.5<z<6.5$ with a sample of very metal-poor DLAs at $2<z<4.5$ shows in general similar average values for the relative abundances, with the exception of [C/O], [Si/Fe] and [Si/O] which are significantly larger for the high-$z$ sample. Furthermore, the dispersion of the measurements significantly increases in the high-redshift bin. This increase is predicted by the theoretical models and indicates a potential retention of PopIII signatures in the probed gas. (Abridged) △ Less

Submitted 4 June, 2024; v1 submitted 16 April, 2024; originally announced April 2024.

Comments: Accepted by A&A

E-XQR-30: The evolution of MgII, CII and OI across 2<z<6

Authors: Alma Maria Sebastian, Emma Ryan-Weber, Rebecca L. Davies, George D. Becker, Laura C. Keating, Valentina D'Odorico, Romain A. Meyer, Sarah E. I. Bosman, Guido Cupani, Girish Kulkarni, Martin G. Haehnelt, Samuel Lai, Anna-Christina Eilers, Manuela Bischetti, Simona Gallerani

Abstract: Intervening metal absorbers in quasar spectra at $z > 6$ can be used as probes to study the chemical enrichment of the Universe during the Epoch of Reionization (EoR). This work presents the comoving line densities ($dn/dX$) of low ionisation absorbers, namely, Mg II (2796Å), C II (1334Å) and O I (1302Å) across $2 <z < 6$ using the E-XQR-30 metal absorber catalog prepared from 42 XSHOOTER quasar s… ▽ More Intervening metal absorbers in quasar spectra at $z > 6$ can be used as probes to study the chemical enrichment of the Universe during the Epoch of Reionization (EoR). This work presents the comoving line densities ($dn/dX$) of low ionisation absorbers, namely, Mg II (2796Å), C II (1334Å) and O I (1302Å) across $2 <z < 6$ using the E-XQR-30 metal absorber catalog prepared from 42 XSHOOTER quasar spectra at $5.8 < z < 6.6$. Here, we analyse 280 Mg II ($1.9 < z < 6.4$), 22 C II ($5.2 < z < 6.4$) and 10 O I ($5.3 < z < 6.4$) intervening absorbers, thereby building up on previous studies with improved sensitivity of 50% completeness at an equivalent width of $W > 0.03$Å. For the first time, we present the comoving line densities of 131 weak ($W < 0.3$Å) intervening Mg II absorbers at $1.9 < z < 6.4$ which exhibit constant evolution with redshift similar to medium ($0.3 < W < 1.0$Å) absorbers. However, the cosmic mass density of Mg II - dominated by strong Mg II systems - traces the evolution of global star formation history from redshift 1.9 to 5.5. E-XQR-30 also increases the absorption path length by a factor of 50% for C II and O I whose line densities show a rising trend towards $z > 5$, in agreement with previous works. In the context of a decline in metal enrichment of the Universe at $z > 5$, the overall evolution in the incidence rates of absorption systems can be explained by a weak - possibly soft fluctuating - UV background. Our results, thereby, provide evidence for a late reionization continuing to occur in metal-enriched and therefore, biased regions in the Universe. △ Less

Submitted 15 March, 2024; originally announced March 2024.

Comments: 20 pages, 21 figures

Chemical Diagnostics to Unveil Environments Enriched by First Stars

Authors: Irene Vanni, Stefania Salvadori, Valentina D'Odorico, George D. Becker, Guido Cupani

Abstract: Unveiling the chemical fingerprints of the first (Pop III) stars is crucial for indirectly studying their properties and probing their massive nature. In particular, very massive Pop III stars explode as energetic Pair-Instability Supernovae (PISNe), allowing their chemical products to escape in the diffuse medium around galaxies, opening the possibility to observe their fingerprints in distant ga… ▽ More Unveiling the chemical fingerprints of the first (Pop III) stars is crucial for indirectly studying their properties and probing their massive nature. In particular, very massive Pop III stars explode as energetic Pair-Instability Supernovae (PISNe), allowing their chemical products to escape in the diffuse medium around galaxies, opening the possibility to observe their fingerprints in distant gas clouds. Recently, three z > 6.3 absorbers with abundances consistent with an enrichment from PISNe have been observed with JWST. In this Letter, we present novel chemical diagnostics to uncover environments mainly imprinted by PISNe. Furthermore, we revise the JWST low-resolution measurements by analyzing the publicly available high-resolution X-Shooter spectra for two of these systems. Our results reconcile the chemical abundances of these absorbers with those from literature, which are found to be consistent with an enrichment dominated (> 50% metals) by normal Pop II SNe. We show the power of our novel diagnostics in isolating environments uniquely enriched by PISNe from those mainly polluted by other Pop III and Pop II SNe. When the subsequent enrichment from Pop II SNe is included, however, we find that the abundances of PISN-dominated environments partially overlap with those predominantly enriched by other Pop III and Pop II SNe. We dub these areas confusion regions. Yet, the odd-even abundance ratios [Mg,Si/Al] are extremely effective in pinpointing PISN- dominated environments and allowed us to uncover, for the first time, an absorber consistent with a combined enrichment by a PISN and another Pop III SN for all the six measured elements. △ Less

Submitted 16 May, 2024; v1 submitted 28 February, 2024; originally announced February 2024.

Comments: 16 pages, 10 figures, accepted for publication on ApJL

Constraints on the Evolution of the Ionizing Background and Ionizing Photon Mean Free Path at the End of Reionization

Authors: Frederick B. Davies, Sarah E. I. Bosman, Prakash Gaikwad, Fahad Nasir, Joseph F. Hennawi, George D. Becker, Martin G. Haehnelt, Valentina D'Odorico, Manuela Bischetti, Anna-Christina Eilers, Laura C. Keating, Girish Kulkarni, Samuel Lai, Chiara Mazzucchelli, Yuxiang Qin, Sindhu Satyavolu, Feige Wang, Jinyi Yang, Yongda Zhu

Abstract: The variations in Ly$α$ forest opacity observed at $z>5.3$ between lines of sight to different background quasars are too strong to be caused by fluctuations in the density field alone. The leading hypothesis for the cause of this excess variance is a late, ongoing reionization process at redshifts below six. Another model proposes strong ionizing background fluctuations coupled to a short, spatia… ▽ More The variations in Ly$α$ forest opacity observed at $z>5.3$ between lines of sight to different background quasars are too strong to be caused by fluctuations in the density field alone. The leading hypothesis for the cause of this excess variance is a late, ongoing reionization process at redshifts below six. Another model proposes strong ionizing background fluctuations coupled to a short, spatially varying mean free path of ionizing photons, without explicitly invoking incomplete reionization. With recent observations suggesting a short mean free path at $z\sim6$, and a dramatic improvement in $z>5$ Ly$α$ forest data quality, we revisit this latter possibility. Here we apply the likelihood-free inference technique of approximate Bayesian computation to jointly constrain the hydrogen photoionization rate $Γ_{\rm HI}$ and the mean free path of ionizing photons $λ_{\rm mfp}$ from the effective optical depth distributions at $z=5.0$-$6.1$ from XQR-30. We find that the observations are well-described by fluctuating mean free path models with average mean free paths that are consistent with the steep trend implied by independent measurements at $z\sim5$-$6$, with a concomitant rapid evolution of the photoionization rate. △ Less

Submitted 13 December, 2023; originally announced December 2023.

Comments: 23 pages, 13 figures, resubmitted to ApJ after referee's comments

The effect of reionization on direct measurements of the mean free path

Authors: Joshua T. Roth, Anson D'Aloisio, Christopher Cain, Bayu Wilson, Yongda Zhu, George D. Becker

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 th… ▽ More 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~\%$ 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 $<~50\%$ level to assumptions about quasar environments and lifetimes for realistic models. We demonstrate that future analyses with improved data may require explicitly modeling 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 modeled with optically thin simulations in previous MFP analyses. Neutral islands can suppress transmission at $λ_{\rm rest}<912$ Angstroms significantly, up to a factor of 2 for $z_{\rm qso}=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. △ Less

Submitted 23 April, 2024; v1 submitted 10 November, 2023; originally announced November 2023.

Comments: Accepted to MNRAS, minor changes

A SPectroscopic survey of biased halos In the Reionization Era (ASPIRE): JWST Discovers an Overdensity around a Metal Absorption-selected Galaxy at $z\sim5.5$

Authors: Yunjing Wu, Feige Wang, Zheng Cai, Xiaohui Fan, Kristian Finlator, Jinyi Yang, Joseph F. Hennawi, Fengwu Sun, Jaclyn B. Champagne, Xiaojing Lin, Zihao Li, Zuyi Chen, Eduardo Bañados, George D. Becker, Sarah E. I. Bosman, Gstavo Bruzual, Stephane Charlot, Hsiao-Wen Chen, Jacopo Chevallard, Anna-Christina Eilers, Emanuele Paolo Farina, Xiangyu Jin, Hyunsung D. Jun, Koki Kakiichi, Mingyu Li , et al. (5 additional authors not shown)

Abstract: The launch of ${\it JWST}$ opens a new window for studying the connection between metal-line absorbers and galaxies at the end of the Epoch of Reionization (EoR). Previous studies have detected absorber-galaxy pairs in limited quantities through ground-based observations. To enhance our understanding of the relationship between absorbers and their host galaxies at $z>5$, we utilized the NIRCam Wid… ▽ More The launch of ${\it JWST}$ opens a new window for studying the connection between metal-line absorbers and galaxies at the end of the Epoch of Reionization (EoR). Previous studies have detected absorber-galaxy pairs in limited quantities through ground-based observations. To enhance our understanding of the relationship between absorbers and their host galaxies at $z>5$, we utilized the NIRCam Wide Field Slitless Spectroscopy (WFSS) to search for absorber-associated galaxies by detecting their rest-frame optical emission lines (e.g., [OIII] + H$β$). We report the discovery of a MgII-associated galaxy at $z=5.428$ using data from the ${\it JWST}$ ASPIRE program. The MgII absorber is detected on the spectrum of quasar J0305--3150 with a rest-frame equivalent width of 0.74$\mathring{A}$. The associated galaxy has an [OIII] luminosity of $10^{42.5}\ {\rm erg\ s^{-1}}$ with an impact parameter of 24.9 proper kiloparsecs (pkpc). The joint ${\it HST}$-${\it JWST}$ spectral energy distribution (SED) implies a stellar mass and star-formation rate of ${\rm M_* \approx 10^{8.8}}$ ${\rm M_{\odot}}$, ${\rm SFR}\approx 10\ {\rm M_{\odot}\ yr^{-1}}$. Its [OIII] equivalent width and stellar mass are typical of [OIII] emitters at this redshift. Furthermore, connecting the outflow starting time to the SED-derived stellar age, the outflow velocity of this galaxy is $\sim300\ {\rm km\ s^{-1}}$, consistent with theoretical expectations. We identified six additional [OIII] emitters with impact parameters of up to $\sim300$ pkpc at similar redshifts ($|dv|<1000\ {\rm km\ s^{-1}}$). The observed number is consistent with that in cosmological simulations. This pilot study suggests that systematically investigating the absorber-galaxy connection within the ASPIRE program will provide insights into the metal-enrichment history in the early universe. △ Less

Submitted 8 November, 2023; v1 submitted 28 September, 2023; originally announced September 2023.

Comments: Accepted for publication in ApJL. Main text 8 pages, 4 figures. For more information of the JWST ASPIRE program please check https://aspire-quasar.github.io/index.html

Unveiling Dark Matter free-streaming at the smallest scales with high redshift Lyman-alpha forest

Authors: Vid Iršič, Matteo Viel, Martin G. Haehnelt, James S. Bolton, Margherita Molaro, Ewald Puchwein, Elisa Boera, George D. Becker, Prakash Gaikwad, Laura C. Keating, Girish Kulkarni

Abstract: This study introduces novel constraints on the free-streaming of thermal relic warm dark matter (WDM) from Lyman-$α$ forest flux power spectra. Our analysis utilises a high-resolution, high-redshift sample of quasar spectra observed using the HIRES and UVES spectrographs ($z=4.2-5.0$). We employ a Bayesian inference framework and a simulation-based likelihood that encompasses various parameters in… ▽ More This study introduces novel constraints on the free-streaming of thermal relic warm dark matter (WDM) from Lyman-$α$ forest flux power spectra. Our analysis utilises a high-resolution, high-redshift sample of quasar spectra observed using the HIRES and UVES spectrographs ($z=4.2-5.0$). We employ a Bayesian inference framework and a simulation-based likelihood that encompasses various parameters including the free-streaming of dark matter, cosmological parameters, the thermal history of the intergalactic medium, and inhomogeneous reionization, to establish lower limits on the mass of a thermal relic WDM particle of $5.7\;\mathrm{keV}$ (at 95\% C.L.). This result surpasses previous limits from the Lyman-$α$ forest through reduction of the measured uncertainties due to a larger statistical sample and by measuring clustering to smaller scales ($k_{\rm max}=0.2\;\mathrm{km^{-1}\,s}$). The approximately two-fold improvement due to the expanded statistical sample suggests that the effectiveness of Lyman-$α$ forest constraints on WDM models at high redshifts are limited by the availability of high-quality quasar spectra. Restricting the analysis to comparable scales and thermal history priors as in prior studies ($k_{\rm max}<0.1\;\mathrm{km^{-1}\,s}$) lowers the bound on the WDM mass to $4.1\;\mathrm{keV}$. As the precision of the measurements increases, it becomes crucial to examine the instrumental and modelling systematics. On the modelling front, we argue that the impact of the thermal history uncertainty on the WDM particle mass constraint has diminished due to improved independent observations. At the smallest scales, the primary source of modeling systematic arises from the structure in the peculiar velocity of the intergalactic medium and inhomogeneous reionization. △ Less

Submitted 15 February, 2024; v1 submitted 8 September, 2023; originally announced September 2023.

Comments: 22 pages, 14 figures, 2 tables; published in PRD on Feb 8 2024

Journal ref: Phys. Rev. D, 109, 043511 (2024)

Measurements of the $z > 5$ Lyman-$α$ forest flux auto-correlation functions from the extended XQR-30 data set

Authors: Molly Wolfson, Joseph F. Hennawi, Sarah E. I. Bosman, Frederick B. Davies, Zarija Lukić, George D. Becker, Huanqing Chen, Guido Cupani, Valentina D'Odorico, Anna-Christina Eilers, Martin G. Haehnelt, Laura C. Keating, Girish Kulkarni, Samuel Lai, Andrei Mesinger, Fabian Walter, Yongda Zhu

Abstract: Recently, the Lyman-$α$ (Ly$α$) forest flux auto-correlation function has been shown to be sensitive to the mean free path of hydrogen-ionizing photons, $λ_{\text{mfp}}$, for simulations at $z \geq 5.4$. Measuring $λ_{\text{mfp}}$ at these redshifts will give vital information on the ending of reionization. Here we present the first observational measurements of the Ly$α$ forest flux auto-correlat… ▽ More Recently, the Lyman-$α$ (Ly$α$) forest flux auto-correlation function has been shown to be sensitive to the mean free path of hydrogen-ionizing photons, $λ_{\text{mfp}}$, for simulations at $z \geq 5.4$. Measuring $λ_{\text{mfp}}$ at these redshifts will give vital information on the ending of reionization. Here we present the first observational measurements of the Ly$α$ forest flux auto-correlation functions in ten redshift bins from $5.1 \leq z \leq 6.0$. We use a sample of 35 quasar sightlines at $z > 5.7$ from the extended XQR-30 data set, this data has signal-to-noise ratios of $> 20$ per spectral pixel. We carefully account for systematic errors in continuum reconstruction, instrumentation, and contamination by damped Ly$α$ systems. With these measurements, we introduce software tools to generate auto-correlation function measurements from any simulation. For an initial comparison, we show our auto-correlation measurements with simulation models for recently measured $λ_{\text{mfp}}$ values and find good agreements. Further work in modeling and understanding the covariance matrices of the data is necessary to get robust measurements of $λ_{\text{mfp}}$ from this data. △ Less

Submitted 6 September, 2023; originally announced September 2023.

Comments: 17 pages, 13 figures, submitted to MNRAS

The relationship between IGM Lyman-alpha opacity and galaxy density near the end of reionization

Authors: Holly M. Christenson, George D. Becker, Anson D'Aloisio, Frederick B. Davies, Yongda Zhu, Elisa Boera, Fahad Nasir, Steven R. Furlanetto, Matthew A. Malkan

Abstract: Observed scatter in the Lyman-alpha opacity of quasar sightlines at $z<6$ has motivated measurements of the correlation between Ly$α$ opacity and galaxy density, as models that predict this scatter make strong and sometimes opposite predictions for how they should be related. Our previous work associated two highly opaque Ly$α$ troughs at $z\sim5.7$ with a deficit of Lyman-$α$ emitting galaxies (L… ▽ More Observed scatter in the Lyman-alpha opacity of quasar sightlines at $z<6$ has motivated measurements of the correlation between Ly$α$ opacity and galaxy density, as models that predict this scatter make strong and sometimes opposite predictions for how they should be related. Our previous work associated two highly opaque Ly$α$ troughs at $z\sim5.7$ with a deficit of Lyman-$α$ emitting galaxies (LAEs). In this work, we survey two of the most highly transmissive lines of sight at this redshift, towards the $z=6.02$ quasar SDSS J1306+0356 and the $z=6.17$ quasar PSO J359-06. We find that both fields are underdense in LAEs within 10 $h^{-1}$ Mpc of the quasar sightline, somewhat less extensive than underdensities associated with Ly$α$ troughs. We combine our observations with three additional fields from the literature, and find that while fields with extreme opacities are generally underdense, moderate opacities span a wider density range. The results at high opacities are consistent with models that invoke UV background fluctuations and/or late reionization to explain the observed scatter in IGM Ly$α$ opacities. There is tension at low opacities, however, as the models tend to associate lower IGM Ly$α$ opacities with higher densities. Although the number of fields surveyed is still small, the low-opacity results may support a scenario in which the ionizing background in low-density regions increases more rapidly than some models suggest after becoming ionized. Elevated gas temperatures from recent reionization may also be making these regions more transparent. △ Less

Submitted 24 August, 2023; originally announced August 2023.

Comments: 20 pages, 17 figures. Accepted to ApJ

Probing Ultra-late Reionization: Direct Measurements of the Mean Free Path over $5<z<6$

Authors: Yongda Zhu, George D. Becker, Holly M. Christenson, Anson D'Aloisio, Sarah E. I. Bosman, Tom Bakx, Valentina D'Odorico, Manuela Bischetti, Christopher Cain, Frederick B. Davies, Rebecca L. Davies, Anna-Christina Eilers, Xiaohui Fan, Prakash Gaikwad, Martin G. Haehnelt, Laura C. Keating, Girish Kulkarni, Samuel Lai, Hai-Xia Ma, Andrei Mesinger, Yuxiang Qin, Sindhu Satyavolu, Tsutomu T. Takeuchi, Hideki Umehata, Jinyi Yang

Abstract: The mean free path of ionizing photons, $λ_{\rm mfp}$, is a critical parameter for modeling the intergalactic medium (IGM) both during and after reionization. We present direct measurements of $λ_{\rm mfp}$ from QSO spectra over the redshift range $5<z<6$, including the first measurements at $z\simeq5.3$ and 5.6. Our sample includes data from the XQR-30 VLT large program, as well as new Keck/ESI o… ▽ More The mean free path of ionizing photons, $λ_{\rm mfp}$, is a critical parameter for modeling the intergalactic medium (IGM) both during and after reionization. We present direct measurements of $λ_{\rm mfp}$ from QSO spectra over the redshift range $5<z<6$, including the first measurements at $z\simeq5.3$ and 5.6. Our sample includes data from the XQR-30 VLT large program, as well as new Keck/ESI observations of QSOs near $z \sim 5.5$, for which we also acquire new [C II] 158$μ$m redshifts with ALMA. By measuring the Lyman continuum transmission profile in stacked QSO spectra, we find $λ_{\rm mfp} = 9.33_{-1.80}^{+2.06}$, $5.40_{-1.40}^{+1.47}$, $3.31_{-1.34}^{+2.74}$, and $0.81_{-0.48}^{+0.73}$ pMpc at $z=5.08$, 5.31, 5.65, and 5.93, respectively. Our results demonstrate that $λ_{\rm mfp}$ increases steadily and rapidly with time over $5<z<6$. Notably, we find that $λ_{\rm mfp}$ deviates significantly from predictions based on a fully ionized and relaxed IGM as late as $z=5.3$. By comparing our results to model predictions and indirect $λ_{\rm mfp}$ constraints based on IGM Ly$α$ opacity, we find that the $λ_{\rm mfp}$ evolution is consistent with scenarios wherein the IGM is still undergoing reionization and/or retains large fluctuations in the ionizing UV background well below redshift six. △ Less

Submitted 8 August, 2023; originally announced August 2023.

Comments: 17 pages, 6 figures, 2 tables; accepted for publication in ApJ

XQR-30: Black Hole Masses and Accretion Rates of 42 z>6 Quasars

Authors: C. Mazzucchelli, M. Bischetti, V. D'Odorico, C. Feruglio, J. -T. Schindler, M. Onoue, E. Bañados, G. D. Becker, F. Bian, S. Carniani, R. Decarli, A. -C. Eilers, E. P. Farina, S. Gallerani, S. Lai, R. A. Meyer, S. Rojas-Ruiz, S. Satyavolu, B. P. Venemans, F. Wang, J. Yang, Y. Zhu

Abstract: We present bolometric luminosities, black hole masses and Eddington ratios for 42 luminous quasars at z>6 using high signal-to-noise ratio VLT/X-Shooter spectra, acquired in the enlarged ESO Large Programme XQR-30. In particular, we derive bolometric luminosities from the rest-frame 3000 A, luminosities using a bolometric correction from the literature, and the black hole masses by modelling the s… ▽ More We present bolometric luminosities, black hole masses and Eddington ratios for 42 luminous quasars at z>6 using high signal-to-noise ratio VLT/X-Shooter spectra, acquired in the enlarged ESO Large Programme XQR-30. In particular, we derive bolometric luminosities from the rest-frame 3000 A, luminosities using a bolometric correction from the literature, and the black hole masses by modelling the spectral regions around the CIV 1549A and the MgII 2798A emission lines, with scaling relations calibrated in the local universe. We find that the black hole masses derived from both emission lines are in the same range, and the scatter of the measurements agrees with expectations from the scaling relations. The MgII-derived masses are between ~(0.8-12) x 10^9 Msun, and the derived Eddington ratios are within ~0.13-1.73, with a mean (median) of 0.84 (0.72). By comparing the total sample of quasars at z>5.8, from this work and from the literature, to a bolometric luminosity distribution-matched sample at z~1.5, we find that quasars at high redshift host slightly less massive black holes which accrete slightly more rapidly than at lower-z, with a difference in the mean Eddington ratios of the two samples of ~0.27, in agreement with recent literature work. △ Less

Submitted 28 June, 2023; originally announced June 2023.

Comments: 9 pages; 5 figures; accepted for publication in A&A

Journal ref: A&A 676, A71 (2023)

A comprehensive metallicity analysis of J0332-3557: establishing a z ~ 4 anchor for direct gas metallicity and C/O abundance investigations

Authors: Annalisa Citro, Danielle A. Berg, Dawn K. Erb, Matthew W. Auger, George D. Becker, Bethan L. James, Evan D. Skillman

Abstract: We provide one of the most comprehensive metallicity studies at z $\sim$ 4 by analyzing the UV/optical HST photometry, and rest-frame VLT-FORS2 ultraviolet and VLT-XSHOOTER optical spectra of J0332-3557, a gravitationally lensed galaxy magnified by a factor of 20. With a 5$σ$ detection of the auroral O III]1666 line, we are able to derive a direct gas metallicity estimate for our target. We find Z… ▽ More We provide one of the most comprehensive metallicity studies at z $\sim$ 4 by analyzing the UV/optical HST photometry, and rest-frame VLT-FORS2 ultraviolet and VLT-XSHOOTER optical spectra of J0332-3557, a gravitationally lensed galaxy magnified by a factor of 20. With a 5$σ$ detection of the auroral O III]1666 line, we are able to derive a direct gas metallicity estimate for our target. We find Z_gas = 12 + log(O/H) = 8.26 $\pm$ 0.06, which is compatible with an increasing of both the gas fraction and the outflow metal loading factor from z $\sim$ 0 to z $\sim$ 4. J0332-3557 is the most metal-rich individual galaxy at z $\sim$ 4 for which the C/O ratio has been measured. We derive a low log(C/O) = - 1.02 $\pm$ 0.2, which suggests that J0332-3557 is in the early stages of ISM carbon enrichment driven mostly by massive stars. The low C/O abundance also indicates that J0332-3557 is characterized by a low star formation efficiency, higher yields of oxygen, and longer burst duration. We find that EW(C III])1906,9 is as low as $\sim$ 3 Å. The main drivers of the low EW(C III])1906,9 are the higher gas metallicity and the low C/O abundance. J0332-3557 is characterized by one diffuse and two more compact regions $\sim$ 1 kpc in size. We find that the carbon emission mostly originates in the compact knots. Our study on J0332-3557 serves as an anchor for studies investigating the evolution of metallicity and C/O abundance across different redshifts. △ Less

Submitted 30 April, 2024; v1 submitted 23 May, 2023; originally announced May 2023.

Comments: 27 pages, 17 figures. Accepted for publication in ApJ

XQR-30: the ultimate XSHOOTER quasar sample at the reionization epoch

Authors: Valentina D'Odorico, E. Banados, G. D. Becker, M. Bischetti, S. E. I. Bosman, G. Cupani, R. Davies, E. P. Farina, A. Ferrara, C. Feruglio, C. Mazzucchelli, E. Ryan-Weber, J. -T. Schindler, A. Sodini, B. P. Venemans, F. Walter, H. Chen, S. Lai, Y. Zhu, F. Bian, S. Campo, S. Carniani, S. Cristiani, F. Davies, R. Decarli , et al. (24 additional authors not shown)

Abstract: The final phase of the reionization process can be probed by rest-frame UV absorption spectra of quasars at z>6, shedding light on the properties of the diffuse intergalactic medium within the first Gyr of the Universe. The ESO Large Programme "XQR-30: the ultimate XSHOOTER legacy survey of quasars at z~5.8-6.6" dedicated ~250 hours of observations at the VLT to create a homogeneous and high-quali… ▽ More The final phase of the reionization process can be probed by rest-frame UV absorption spectra of quasars at z>6, shedding light on the properties of the diffuse intergalactic medium within the first Gyr of the Universe. The ESO Large Programme "XQR-30: the ultimate XSHOOTER legacy survey of quasars at z~5.8-6.6" dedicated ~250 hours of observations at the VLT to create a homogeneous and high-quality sample of spectra of 30 luminous quasars at z~6, covering the rest wavelength range from the Lyman limit to beyond the MgII emission. Twelve quasar spectra of similar quality from the XSHOOTER archive were added to form the enlarged XQR-30 sample, corresponding to a total of ~350 hours of on-source exposure time. The median effective resolving power of the 42 spectra is R~11400 and 9800 in the VIS and NIR arm, respectively. The signal-to-noise ratio per 10 km/s pixel ranges from ~11 to 114 at $λ\simeq 1285$ Årest frame, with a median value of ~29. We describe the observations, data reduction and analysis of the spectra, together with some first results based on the E-XQR-30 sample. New photometry in the H and K bands are provided for the XQR-30 quasars, together with composite spectra whose characteristics reflect the large absolute magnitudes of the sample. The composite and the reduced spectra are released to the community through a public repository, and will enable a range of studies addressing outstanding questions regarding the first Gyr of the Universe. △ Less

Submitted 19 June, 2023; v1 submitted 8 May, 2023; originally announced May 2023.

Comments: 21 pages, 10 figures. Final version accepted by MNRAS

Journal ref: MNRAS, Volume 523, Issue 1, pp.1399-1420 (2023)

Evidence of First Stars-enriched Gas in High-redshift Absorbers

Authors: A. Saccardi, S. Salvadori, V. D'Odorico, G. Cupani, M. Fumagalli, T. A. M. Berg, G. D. Becker, S. Ellison, S. Lopez

Abstract: The first stars were born from chemically pristine gas. They were likely massive, and thus they rapidly exploded as supernovae, enriching the surrounding gas with the first heavy elements. In the Local Group, the chemical signatures of the first stellar population were identified among low-mass, long-lived, very metal-poor ([Fe/H]<-2) stars, characterized by high abundances of carbon over iron ([C… ▽ More The first stars were born from chemically pristine gas. They were likely massive, and thus they rapidly exploded as supernovae, enriching the surrounding gas with the first heavy elements. In the Local Group, the chemical signatures of the first stellar population were identified among low-mass, long-lived, very metal-poor ([Fe/H]<-2) stars, characterized by high abundances of carbon over iron ([C/Fe]>+0.7): the so-called carbon-enhanced metal-poor stars. Conversely, a similar carbon excess caused by first-star pollution was not found in dense neutral gas traced by absorption systems at different cosmic time. Here we present the detection of 14 very metal-poor, optically thick absorbers at redshift z~3-4. Among these, 3 are carbon-enhanced and reveal an overabundance with respect to Fe of all the analyzed chemical elements (O, Mg, Al, and Si). Their relative abundances show a distribution with respect to [Fe/H] that is in very good agreement with those observed in nearby very metal-poor stars. All the tests we performed support the idea that these C-rich absorbers preserve the chemical yields of the first stars. Our new findings suggest that the first-star signatures can survive in optically thick but relatively diffuse absorbers, which are not sufficiently dense to sustain star formation and hence are not dominated by the chemical products of normal stars. △ Less

Submitted 3 May, 2023; originally announced May 2023.

New quasar proximity zone size measurements at $z\sim 6$ using the enlarged XQR-30 sample

Authors: Sindhu Satyavolu, Anna-Christina Eilers, Girish Kulkarni, Emma Ryan-Weber, Rebecca L. Davies, George D. Becker, Sarah E. I. Bosman, Bradley Greig, Chiara Mazzucchelli, Eduardo Bañados, Manuela Bischetti, Valentina D'Odorico, Xiaohui Fan, Emanuele Paolo Farina, Martin G. Haehnelt, Laura C. Keating, Samuel Lai, Fabian Walter

Abstract: Proximity zones of high-redshift quasars are unique probes of their central supermassive black holes as well as the intergalactic medium in the last stages of reionization. We present 22 new measurements of proximity zones of quasars with redshifts between 5.8 and 6.6, using the enlarged XQR-30 sample of high-resolution, high-SNR quasar spectra. The quasars in our sample have UV magnitudes of… ▽ More Proximity zones of high-redshift quasars are unique probes of their central supermassive black holes as well as the intergalactic medium in the last stages of reionization. We present 22 new measurements of proximity zones of quasars with redshifts between 5.8 and 6.6, using the enlarged XQR-30 sample of high-resolution, high-SNR quasar spectra. The quasars in our sample have UV magnitudes of $M_{1450}\sim -27$ and black hole masses of $10^9$$\unicode{x2013}$$10^{10}$ M$_\odot$. Our inferred proximity zone sizes are 2$\unicode{x2013}$7 physical Mpc, with a typical uncertainty of less than 0.5 physical Mpc, which, for the first time, also includes uncertainty in the quasar continuum. We find that the correlation between proximity zone sizes and the quasar redshift, luminosity, or black hole mass, indicates a large diversity of quasar lifetimes. Two of our proximity zone sizes are exceptionally small. The spectrum of one of these quasars, with $z=6.02$, displays, unusually for this redshift, damping wing absorption without any detectable metal lines, which could potentially originate from the IGM. The other quasar has a high-ionization absorber $\sim$0.5 pMpc from the edge of the proximity zone. This work increases the number of proximity zone measurements available in the last stages of cosmic reionization to 87. This data will lead to better constraints on quasar lifetimes and obscuration fractions at high redshift, which in turn will help probe the seed mass and formation redshift of supermassive black holes. △ Less

Submitted 1 May, 2023; originally announced May 2023.

Comments: 16 pages, 9 figures, Accepted in MNRAS

Report number: TIFR/TH/23-5

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

Authors: Prakash Gaikwad, Martin G. Haehnelt, Frederick B. Davies, Sarah E. I. Bosman, Margherita Molaro, Girish Kulkarni, Valentina D'Odorico, George D. Becker, Rebecca L. Davies, Fahad Nasir, James S. Bolton, Laura C. Keating, Vid Iršič, Ewald Puchwein, Yongda Zhu, Shikhar Asthana, Jinyi Yang, Samuel Lai, Anna-Christina Eilers

Abstract: We measure the mean free path ($λ_{\rm mfp,HI}$), photo-ionization rate ($\langle Γ_{\rm HI} \rangle$) and neutral fraction ($\langle f_{\rm HI} \rangle$) of hydrogen in 12 redshift bins at $4.85<z<6.05$ from a large sample of moderate resolution XShooter and ESI QSO absorption spectra. The fluctuations in ionizing radiation field are modeled by post-processing simulations from the Sherwood suite… ▽ More We measure the mean free path ($λ_{\rm mfp,HI}$), photo-ionization rate ($\langle Γ_{\rm HI} \rangle$) and neutral fraction ($\langle f_{\rm HI} \rangle$) of hydrogen in 12 redshift bins at $4.85<z<6.05$ from a large sample of moderate resolution XShooter and ESI QSO absorption spectra. The fluctuations in ionizing radiation field are modeled by post-processing simulations from the Sherwood suite using our new code "EXtended reionization based on the Code for Ionization and Temperature Evolution" (EX-CITE). EX-CITE uses efficient Octree summation for computing intergalactic medium attenuation and can generate large number of high resolution $Γ_{\rm HI}$ fluctuation models. Our simulation with EX-CITE shows remarkable agreement with simulations performed with the radiative transfer code Aton and can recover the simulated parameters within $1σ$ uncertainty. We measure the three parameters by forward-modeling the Ly$α$ forest and comparing the effective optical depth ($τ_{\rm eff, HI}$) distribution in simulations and observations. The final uncertainties in our measured parameters account for the uncertainties due to thermal parameters, modeling parameters, observational systematics and cosmic variance. Our best fit parameters show significant evolution with redshift such that $λ_{\rm mfp,HI}$ and $\langle f_{\rm HI} \rangle$ decreases and increases by a factor $\sim 6$ and $\sim 10^{4}$, respectively from $z \sim 5$ to $z \sim 6$. By comparing our $λ_{\rm mfp,HI}$, $\langle Γ_{\rm HI} \rangle$ and $\langle f_{\rm HI} \rangle$ evolution with that in state-of-the-art Aton radiative transfer simulations and the Thesan and CoDa-III simulations, we find that our best fit parameter evolution is consistent with a model in which reionization completes by $z \sim 5.2$. △ Less

Submitted 23 August, 2023; v1 submitted 4 April, 2023; originally announced April 2023.

Comments: 30 pages (+14 pages appendices), 14 figures (+18 figures appendices); Accepted for publication in MNRAS; Main results are summarized in Fig. 10, Fig. 11 and Table 3

Examining the Decline in the C IV Content of the Universe over 4.3 < z < 6.3 using the E-XQR-30 Sample

Authors: Rebecca L. Davies, Emma Ryan-Weber, Valentina D'Odorico, Sarah E. I. Bosman, Romain A. Meyer, George D. Becker, Guido Cupani, Laura C. Keating, Manuela Bischetti, Frederick B. Davies, Anna-Christina Eilers, Emanuele Paolo Farina, Martin G. Haehnelt, Andrea Pallottini, Yongda Zhu

Abstract: Intervening CIV absorbers are key tracers of metal-enriched gas in galaxy halos over cosmic time. Previous studies suggest that the CIV cosmic mass density ($Ω_{\rm CIV}$) decreases slowly over 1.5 $\lesssim z\lesssim$ 5 before declining rapidly at $z\gtrsim$ 5, but the cause of this downturn is poorly understood. We characterize the $Ω_{\rm CIV}$ evolution over 4.3 $\lesssim z\lesssim$ 6.3 using… ▽ More Intervening CIV absorbers are key tracers of metal-enriched gas in galaxy halos over cosmic time. Previous studies suggest that the CIV cosmic mass density ($Ω_{\rm CIV}$) decreases slowly over 1.5 $\lesssim z\lesssim$ 5 before declining rapidly at $z\gtrsim$ 5, but the cause of this downturn is poorly understood. We characterize the $Ω_{\rm CIV}$ evolution over 4.3 $\lesssim z\lesssim$ 6.3 using 260 absorbers found in 42 XSHOOTER spectra of $z\sim$ 6 quasars, of which 30 come from the ESO Large Program XQR-30. The large sample enables us to robustly constrain the rate and timing of the downturn. We find that $Ω_{\rm CIV}$ decreases by a factor of 4.8 $\pm$ 2.0 over the ~300 Myr interval between $z\sim$ 4.7 and $z\sim$ 5.8. The slope of the column density (log N) distribution function does not change, suggesting that CIV absorption is suppressed approximately uniformly across 13.2 $\leq$ log N/cm$^{-2}$ < 15.0. Assuming that the carbon content of galaxy halos evolves as the integral of the cosmic star formation rate density (with some delay due to stellar lifetimes and outflow travel times), we show that chemical evolution alone could plausibly explain the fast decline in $Ω_{\rm CIV}$ over 4.3 $\lesssim z\lesssim$ 6.3. However, the CIV/CII ratio decreases at the highest redshifts, so the accelerated decline in $Ω_{\rm CIV}$ at $z\gtrsim$ 5 may be more naturally explained by rapid changes in the gas ionization state driven by evolution of the UV background towards the end of hydrogen reionization. △ Less

Submitted 5 March, 2023; originally announced March 2023.

Comments: 20 pages, 8 figures, 4 tables. Published in MNRAS

The XQR-30 Metal Absorber Catalog: 778 Absorption Systems Spanning 2 < z < 6.5

Authors: Rebecca L. Davies, Emma Ryan-Weber, Valentina D'Odorico, Sarah E. I. Bosman, Romain A. Meyer, George D. Becker, Guido Cupani, Manuela Bischetti, Alma M. Sebastian, Anna-Christina Eilers, Emanuele Paola Farina, Feige Wang, Jinyi Yang, Yongda Zhu

Abstract: Intervening metal absorption lines in the spectra of z > 6 quasars are fundamental probes of the ionization state and chemical composition of circumgalactic and intergalactic gas near the end of the reionization epoch. Large absorber samples are required to robustly measure typical absorber properties and to refine models of the synthesis, transport, and ionization of metals in the early Universe.… ▽ More Intervening metal absorption lines in the spectra of z > 6 quasars are fundamental probes of the ionization state and chemical composition of circumgalactic and intergalactic gas near the end of the reionization epoch. Large absorber samples are required to robustly measure typical absorber properties and to refine models of the synthesis, transport, and ionization of metals in the early Universe. The "Ultimate XSHOOTER legacy survey of quasars at z~5.8-6.6" (XQR-30) has obtained high signal-to-noise spectra of 30 luminous quasars, nearly quadrupling the existing sample of 12 high quality z~6 quasar spectra. We use this unprecedented sample to construct a catalog of 778 systems showing absorption in one or more of MgII (360 systems), FeII (184), CII (46), CIV (479), SiIV (127), and NV (13) which span 2 < z < 6.5. This catalog significantly expands on existing samples of z > 5 absorbers, especially for CIV and SiIV which are important probes of the ionizing photon background at high redshift. The sample is 50% (90%) complete for rest-frame equivalent widths W > 0.03AA (0.09AA). We publicly release the absorber catalog along with completeness statistics and a Python script to compute the absorption search path for different ions and redshift ranges. This dataset is a key legacy resource for studies of enriched gas from the era of galaxy assembly to cosmic noon, and paves the way for even higher redshift studies with the James Webb Space Telescope and 30m-class telescopes. △ Less

Submitted 19 January, 2023; v1 submitted 28 November, 2022; originally announced November 2022.

Comments: 29 pages and 14 figures including 2 appendices. Accepted for publication in MNRAS. The metal absorber catalogue can be downloaded from https://github.com/XQR-30/Metal-catalogue. Contents unchanged from v1 apart from change to author list

The Sherwood-Relics simulations: overview and impact of patchy reionization and pressure smoothing on the intergalactic medium

Authors: Ewald Puchwein, James S. Bolton, Laura C. Keating, Margherita Molaro, Prakash Gaikwad, Girish Kulkarni, Martin G. Haehnelt, Vid Iršič, Tomáš Šoltinský, Matteo Viel, Dominique Aubert, George D. Becker, Avery Meiksin

Abstract: We present the Sherwood-Relics simulations, a new suite of large cosmological hydrodynamical simulations aimed at modelling the intergalactic medium (IGM) during and after the cosmic reionization of hydrogen. The suite consists of over 200 simulations that cover a wide range of astrophysical and cosmological parameters. It also includes simulations that use a new lightweight hybrid scheme for trea… ▽ More We present the Sherwood-Relics simulations, a new suite of large cosmological hydrodynamical simulations aimed at modelling the intergalactic medium (IGM) during and after the cosmic reionization of hydrogen. The suite consists of over 200 simulations that cover a wide range of astrophysical and cosmological parameters. It also includes simulations that use a new lightweight hybrid scheme for treating radiative transfer effects. This scheme follows the spatial variations in the ionizing radiation field, as well as the associated fluctuations in IGM temperature and pressure smoothing. It is computationally much cheaper than full radiation hydrodynamics simulations and circumvents the difficult task of calibrating a galaxy formation model to observational constraints on cosmic reionization. Using this hybrid technique, we study the spatial fluctuations in IGM properties that are seeded by patchy cosmic reionization. We investigate the relevant physical processes and assess their impact on the z > 4 Lyman-alpha forest. Our main findings are: (i) Consistent with previous studies patchy reionization causes large scale temperature fluctuations that persist well after the end of reionization, (ii) these increase the Lyman-alpha forest flux power spectrum on large scales, and (iii) result in a spatially varying pressure smoothing that correlates well with the local reionization redshift. (iv) Structures evaporated or puffed up by photoheating cause notable features in the Lyman-alpha forest, such as flat-bottom or double-dip absorption profiles. △ Less

Submitted 28 February, 2023; v1 submitted 26 July, 2022; originally announced July 2022.

Comments: Version accepted for publication in MNRAS, some discussions added, data added to Fig. 8. 24 pages, 15 figures. Images and movies may be viewed at https://www.nottingham.ac.uk/astronomy/sherwood-relics/

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

Authors: Yongda Zhu, George D. Becker, Sarah E. I. Bosman, Laura C. Keating, Valentina D'Odorico, Rebecca L. Davies, Holly M. Christenson, Eduardo Bañados, Fuyan Bian, Manuela Bischetti, Huanqing Chen, Frederick B. Davies, Anna-Christina Eilers, Xiaohui Fan, Prakash Gaikwad, Bradley Greig, Martin G. Haehnelt, Girish Kulkarni, Samuel Lai, Andrea Pallottini, Yuxiang Qin, Emma Ryan-Weber, Fabian Walter, Feige Wang, Jinyi Yang

Abstract: We present a new investigation of the intergalactic medium (IGM) near reionization using dark gaps in the Lyman-$β$ (Ly$β$) forest. With its lower optical depth, Ly$β$ offers a potentially more sensitive probe to any remaining neutral gas compared to commonly used Ly$α$ line. We identify dark gaps in the Ly$β$ forest using spectra of 42 QSOs at $z_{\rm em}>5.5$, including new data from the XQR-30… ▽ More We present a new investigation of the intergalactic medium (IGM) near reionization using dark gaps in the Lyman-$β$ (Ly$β$) forest. With its lower optical depth, Ly$β$ offers a potentially more sensitive probe to any remaining neutral gas compared to commonly used Ly$α$ line. We identify dark gaps in the Ly$β$ forest using spectra of 42 QSOs at $z_{\rm em}>5.5$, including new data from the XQR-30 VLT Large Programme. Approximately $40\%$ of these QSO spectra exhibit dark gaps longer than $10h^{-1}{\rm Mpc}$ at $z\simeq5.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 (UVB) 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=28h^{-1}{\rm Mpc}$) and dark ($τ_{\rm eff} \gtrsim 6$) gap persisting down to $z\simeq 5.5$ in the Ly$β$ forest of the $z_{\rm}=5.85$ QSO PSO J025$-$11. This gap may support late reionization models with a volume-weighted average neutral hydrogen fraction of $ \langle x_{\rm HI}\rangle \gtrsim 5\%$ by $z=5.6$. Finally, we infer constraints on $\langle x_{\rm HI}\rangle$ over $5.5 \lesssim z \lesssim 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 $\langle x_{\rm HI}\rangle \leq 0.05$, 0.17, and 0.29 at $z\simeq 5.55$, 5.75, and 5.95, respectively. These constraints are consistent with models where reionization ends significantly later than $z = 6$. △ Less

Submitted 9 May, 2022; originally announced May 2022.

Comments: 17 pages, 13 figures, accepted for publication in ApJ. The Figure set and machine readable table are available at https://ydzhuastro.github.io/lyb.html

Measuring the Density Fields around Bright Quasars at $z\sim 6$ with XQR-30 Spectra

Authors: Huanqing Chen, Anna-Christina Eilers, Sarah E. I. Bosman, Nickolay Y. Gnedin, Xiaohui Fan, Feige Wang, Jinyi Yang, Valentina D'Odorico, George D. Becker, Manuela Bischetti, Chiara Mazzucchelli, Andrei Mesinger, Andrea Pallottini

Abstract: Measuring the density of the intergalactic medium using quasar sightlines in the epoch of reionization is challenging due to the saturation of Lyman-$α$ absorption. Near a luminous quasar, however, the enhanced radiation creates a proximity zone observable in the quasar spectra where the Lyman-$α$ absorption is not saturated. In this study, we use $10$ high-resolution ($R\gtrsim 10,000$)… ▽ More Measuring the density of the intergalactic medium using quasar sightlines in the epoch of reionization is challenging due to the saturation of Lyman-$α$ absorption. Near a luminous quasar, however, the enhanced radiation creates a proximity zone observable in the quasar spectra where the Lyman-$α$ absorption is not saturated. In this study, we use $10$ high-resolution ($R\gtrsim 10,000$) $z\sim 6$ quasar spectra from the extended XQR-30 sample to measure the density field in the quasar proximity zones. We find a variety of environments within $3$ pMpc distance from the quasars. We compare the observed density cumulative distribution function (CDF) with models from the $\textit{Cosmic Reionization on Computers}$ simulation, and find a good agreement between $1.5$ to $3$pMpc from the quasar. This region is far away from the quasar hosts and hence approaching the mean density of the universe, which allows us to use the CDF to set constraints on the cosmological parameter $σ_8=0.6\pm0.3$. The uncertainty is mainly due to the limited number of high-quality quasar sightlines currently available. Utilizing the more than $>200$ known quasars at $z\gtrsim 6$, this method will allow us in the future to tighten the constraint on $σ_8$ to the percent level. In the region closer to the quasar within $1.5$ pMpc, we find the density is higher than predicted in the simulation by $1.23 \pm 0.17$, suggesting the typical host dark matter halo mass of a bright quasar ($M_{\rm 1450}<-26.5$) at $z\sim 6$ is $\log_{\rm 10} (M/M_\odot)=12.5^{+0.4}_{-0.7}$. △ Less

Submitted 26 October, 2021; originally announced October 2021.

Comments: 16 pages, 11 figures, submitted to ApJ, comments are welcome

Constraints on the End of Reionization from the Density Fields Surrounding Two Highly Opaque Quasar Sightlines

Authors: Holly M. Christenson, George D. Becker, Steven R. Furlanetto, Frederick B. Davies, Matthew A. Malkan, Yongda Zhu, Elisa Boera, Adam Trapp

Abstract: The observed large-scale scatter in Lyman $α$ 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$α$-emitti… ▽ More The observed large-scale scatter in Lyman $α$ 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 towards the $z=6.0$ quasar ULAS J0148+0600, for which we re-analyze previously published results using improved photometric selection, and towards 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$. △ Less

Submitted 27 September, 2021; originally announced September 2021.

Comments: Accepted to ApJ. 17 pages, 13 figures

SDSS J1059+4251, a highly magnified z ~ 2.8 star-forming galaxy: ESI observations of the rest-frame UV spectrum

Authors: Annalisa Citro, Dawn K. Erb, Max Pettini, Matthew W. Auger, George D. Becker, Bethan L. James

Abstract: Detailed analyses of high-redshift galaxies are challenging due to their faintness, but this difficulty can be overcome with gravitational lensing, in which the magnification of the flux enables high signal-to-noise ratio (S/N) spectroscopy. We present the rest-frame ultraviolet (UV) Keck Echellette Spectrograph and Imager (ESI) spectrum of the newly discovered z = 2.79 lensed galaxy SDSS J1059+42… ▽ More Detailed analyses of high-redshift galaxies are challenging due to their faintness, but this difficulty can be overcome with gravitational lensing, in which the magnification of the flux enables high signal-to-noise ratio (S/N) spectroscopy. We present the rest-frame ultraviolet (UV) Keck Echellette Spectrograph and Imager (ESI) spectrum of the newly discovered z = 2.79 lensed galaxy SDSS J1059+4251. With an observed magnitude F814W = 18.8 and a magnification factor μ= 31 \pm 3, J1059+4251 is both highly magnified and intrinsically luminous, about two magnitudes brighter than M* at z ~ 2-3. With stellar mass M* = (3.22 \pm 0.20) \times 10^10 M_sun, UV star formation rate SFR=50 \pm 7 M_sun yr^-1, and stellar metallicity Z*~ 0.15-0.5 Z_sun, J1059+4251 is typical of bright star-forming galaxies at similar redshifts. Thanks to the high S/N and the spectral resolution of the ESI spectrum, we are able to separate the interstellar and stellar features and derive properties that would be inaccessible without the aid of the lensing. We find evidence of a gas outflow with speeds up to -1000 km s^-1, and of an inflow that is probably due to accreting material seen along a favorable line of sight. We measure relative elemental abundances from the interstellar absorption lines and find that alpha-capture elements are overabundant compared to iron-peak elements, suggestive of rapid star formation. However, this trend may also be affected by dust depletion. Thanks to the high data quality, our results represent a reliable step forward in the characterization of typical galaxies at early cosmic epochs. △ Less

Submitted 22 September, 2021; v1 submitted 20 September, 2021; originally announced September 2021.

Comments: 23 pages, 13 figures, Accepted for publication in ApJ

Chasing the Tail of Cosmic Reionization with Dark Gap Statistics in the Ly$α$ Forest over $5 < z < 6$

Authors: Yongda Zhu, George D. Becker, Sarah E. I. Bosman, Laura C. Keating, Holly M. Christenson, Eduardo Bañados, Fuyan Bian, Frederick B. Davies, Valentina D'Odorico, Anna-Christina Eilers, Xiaohui Fan, Martin G. Haehnelt, Girish Kulkarni, Andrea Pallottini, Yuxiang Qin, Feige Wang, Jinyi Yang

Abstract: We present a new investigation of the intergalactic medium (IGM) near the end of reionization using "dark gaps" in the Lyman-alpha (Ly$α$) forest. Using spectra of 55 QSOs at $z_{\rm 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}\,{\rm Mpc}$ bins falls below 5%. Nine ultra-long (… ▽ More We present a new investigation of the intergalactic medium (IGM) near the end of reionization using "dark gaps" in the Lyman-alpha (Ly$α$) forest. Using spectra of 55 QSOs at $z_{\rm 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}\,{\rm Mpc}$ bins falls below 5%. Nine ultra-long ($L > 80~h^{-1}\,{\rm Mpc}$) dark gaps are identified at $z<6$. In addition, we quantify the fraction of QSO spectra exhibiting gaps longer than $30~h^{-1}\,{\rm Mpc}$, $F_{30}$, as a function of redshift. We measure $F_{30} \simeq 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 \simeq 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 \lesssim 6$ are also potentially consistent with the data. Overall, our results suggest that signature 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 \simeq 5.3$. △ Less

Submitted 13 September, 2021; originally announced September 2021.

Comments: 30 pages, 19 figures, 3 tables, accepted for publication in ApJ

Hydrogen reionisation ends by $z=5.3$: Lyman-$α$ optical depth measured by the XQR-30 sample

Authors: Sarah E. I. Bosman, Frederick B. Davies, George D. Becker, Laura C. Keating, Rebecca L. Davies, Yongda Zhu, Anna-Christina Eilers, Valentina D'Odorico, Fuyan Bian, Manuela Bischetti, Stefano V. Cristiani, Xiaohui Fan, Emanuele P. Farina, Martin G. Haehnelt, Joseph F. Hennawi, Girish Kulkarni, Andrei Mesinger, Romain A. Meyer, Masafusa Onoue, Andrea Pallottini, Yuxiang Qin, Emma Ryan-Weber, Jan-Torge Schindler, Fabian Walter, Feige Wang , et al. (1 additional authors not shown)

Abstract: The presence of excess scatter in the Ly-$α$ forest at $z\sim 5.5$, together with the existence of sporadic extended opaque Gunn-Peterson troughs, has started to provide robust evidence for a late end of hydrogen reionisation. However, low data quality and systematic uncertainties complicate the use of Ly-$α$ transmission as a precision probe of reionisation's end stages. In this paper, we assembl… ▽ More The presence of excess scatter in the Ly-$α$ forest at $z\sim 5.5$, together with the existence of sporadic extended opaque Gunn-Peterson troughs, has started to provide robust evidence for a late end of hydrogen reionisation. However, low data quality and systematic uncertainties complicate the use of Ly-$α$ transmission as a precision probe of reionisation's end stages. In this paper, we assemble a sample of 67 quasar sightlines at $z>5.5$ with high signal-to-noise ratios of $>10$ per $\leq 15$ km s$^{-1}$ spectral pixel, relying largely on the new XQR-30 quasar sample. XQR-30 is a large program on VLT/X-Shooter which obtained deep (SNR $>20$ per pixel) spectra of 30 quasars at $z>5.7$. We carefully account for systematics in continuum reconstruction, instrumentation, and contamination by damped Ly-$α$ systems. We present improved measurements of the mean Ly-$α$ transmission over $4.9<z<6.1$. Using all known systematics in a forward modelling analysis, we find excellent agreement between the observed Ly-$α$ transmission distributions and the homogeneous-UVB simulations Sherwood and Nyx up to $z\leq5.2$ ($<1 σ$), and mild tension ($\sim 2.5 σ$) at $z=5.3$. Homogeneous UVB models are ruled out by excess Ly-$α$ transmission scatter at $z\geq5.4$ with high confidence ($>3.5 σ$). Our results indicate that reionisation-related fluctuations, whether in the UVB, residual neutral hydrogen fraction, and/or IGM temperature, persist in the intergalactic medium until at least $z=5.3$ ($t=1.1$ Gyr after the Big Bang). This is further evidence for a late end to reionisation. △ Less

Submitted 11 April, 2022; v1 submitted 8 August, 2021; originally announced August 2021.

Comments: Accepted for publication in MNRAS. Full measurement datasets are available on the first author's website

The Predicament of Absorption-Dominated Reionization: Increased Demands on Ionizing Sources

Authors: Frederick B. Davies, Sarah E. I. Bosman, Steven R. Furlanetto, George D. Becker, Anson D'Aloisio

Abstract: The reionization epoch concludes when ionizing photons reach every corner of the Universe. Reionization has generally been assumed to be limited primarily by the rate at which galaxies produce ionizing photons, but the recent measurement of a surprisingly short ionizing photon mean free path of $0.75^{+0.65}_{-0.45}$ proper Mpc at $z = 6$ by Becker et al. (2021) suggests that absorption by residua… ▽ More The reionization epoch concludes when ionizing photons reach every corner of the Universe. Reionization has generally been assumed to be limited primarily by the rate at which galaxies produce ionizing photons, but the recent measurement of a surprisingly short ionizing photon mean free path of $0.75^{+0.65}_{-0.45}$ proper Mpc at $z = 6$ by Becker et al. (2021) suggests that absorption by residual neutral hydrogen in the otherwise ionized intergalactic medium may play a much larger role than previously expected. Here we show that consistency between this short mean free path and the coeval dark pixel fraction in the Ly$α$ forest requires a cumulative output of $6.1^{+11}_{-2.4}$ ionizing photons per baryon by reionization's end. This represents a dramatic increase in the ionizing photon budget over previous estimates, greatly exacerbating the tension with measurements of the ionizing output from galaxies at later times. Translating this constraint into the instantaneous ionizing production from galaxies in our model, we find $\log_{10}f_{\rm esc}ξ_{\rm ion}/\text{(erg/Hz)}^{-1} =25.02_{-0.21}^{+0.45}$ at $z\sim6$. Even with optimistic assumptions about the ionizing production efficiency of early stellar populations, and assuming the galaxy luminosity function extends to extremely faint sources ($M_{\text{UV}}\leq-11$), complete reionization requires the escape fraction of ionizing photons to exceed $20\%$ across the galaxy population. This is far larger than observed in any galaxy population at lower redshifts, requiring rapid evolution in galaxy properties after the first billion years of cosmic time. This tension cannot be completely relieved within existing observational constraints on the hydrogen neutral fraction and mean free path. △ Less

Submitted 21 May, 2021; originally announced May 2021.

Comments: 8 pages, 5 figures, submitted to ApJL

A short mean free path at $z=6$ favors late and rapid reionization by faint galaxies

Authors: Christopher Cain, Anson D'Aloisio, Nakul Gangolli, George D. Becker

Abstract: Recent measurements of the ionizing photon mean free path ($λ_{912}^{\rm mfp}$) at $5 < z < 6$ suggest that the IGM was rapidly evolving at those times. We use radiative transfer simulations to explore the implications for reionization, with a focus on the short value of $λ_{912}^{\rm mfp} = 3.57^{+3.09}_{-2.14}$ cMpc/$h$ at $z=6$. We introduce a model for sub-resolution ionizing photon sinks base… ▽ More Recent measurements of the ionizing photon mean free path ($λ_{912}^{\rm mfp}$) at $5 < z < 6$ suggest that the IGM was rapidly evolving at those times. We use radiative transfer simulations to explore the implications for reionization, with a focus on the short value of $λ_{912}^{\rm mfp} = 3.57^{+3.09}_{-2.14}$ cMpc/$h$ at $z=6$. We introduce a model for sub-resolution ionizing photon sinks based on radiative hydrodynamics simulations of small-scale IGM clumping. We argue that the rapid evolution in $λ_{912}^{\rm mfp}$ at $z=5-6$, together with constraints on the metagalactic ionizing background, favors a late reionization process in which the neutral fraction evolved rapidly in the latter half. We also argue that the short $λ_{912}^{\rm mfp}(z=6)$ points to faint galaxies as the primary drivers of reionizaton. Our preferred model, with $λ_{912}^{\rm mfp}(z=6) = 6.5$ Mpc/$h$, has a midpoint of $z= 7.1$ and ends at $z= 5.1$. It requires 3 ionizing photons per H atom to complete reionization and a LyC photon production efficiency of $\log(f^{\rm eff}_{\rm esc} ξ_{\rm ion}/[\mathrm{erg}^{-1} \mathrm{Hz}]) = 24.8$ at $z>6$. Recovering $λ_{912}^{\rm mfp}(z=6)$ as low as the measured central value may require an increase in IGM clumpiness beyond predictions from simulations, with a commensurate increase in the photon budget. △ Less

Submitted 17 September, 2021; v1 submitted 21 May, 2021; originally announced May 2021.

Comments: 8 pages, 3 figures

The mean free path of ionizing photons at 5 < z < 6: evidence for rapid evolution near reionization

Authors: George D. Becker, Anson D'Aloisio, Holly M. Christenson, Yongda Zhu, Gábor Worseck, James S. Bolton

Abstract: The mean free path of ionizing photons, $λ_{\rm mfp}$, is a key factor in the photoionization of the intergalactic medium (IGM). At $z \gtrsim 5$, however, $λ_{\rm mfp}$ may be short enough that measurements towards QSOs are biased by the QSO proximity effect. We present new direct measurements of $λ_{\rm mfp}$ that address this bias and extend up to $z \sim 6$ for the first time. Our measurements… ▽ More The mean free path of ionizing photons, $λ_{\rm mfp}$, is a key factor in the photoionization of the intergalactic medium (IGM). At $z \gtrsim 5$, however, $λ_{\rm mfp}$ may be short enough that measurements towards QSOs are biased by the QSO proximity effect. We present new direct measurements of $λ_{\rm mfp}$ that address this bias and extend up to $z \sim 6$ for the first time. Our measurements at $z \sim 5$ are based on data from the Giant Gemini GMOS survey and new Keck LRIS observations of low-luminosity QSOs. At $z \sim 6$ we use QSO spectra from Keck ESI and VLT X-Shooter. We measure $λ_{\rm mfp} = 9.09^{+1.62}_{-1.28}$ proper Mpc and $0.75^{+0.65}_{-0.45}$ proper Mpc (68% 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 $λ_{\rm 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% neutral at $z=6$, although our measurement at $z = 6.0$ is even lower than these models prefer. △ Less

Submitted 15 September, 2021; v1 submitted 30 March, 2021; originally announced March 2021.

Comments: Accepted for publication in MNRAS

Sub-damped Lyman alpha systems in the XQ-100 survey II -- Chemical evolution at 2.4<z<4.3

Authors: Trystyn A. M. Berg, Michele Fumagalli, Valentina D'Odorico, Sara L. Ellison, Sebastian Lopez, George D. Becker, Lise Christensen, Guido Cupani, Kelly D. Denney, Ruben Sanchez-Ramirez, Gabor Worseck

Abstract: We present the measured gas-phase metal column densities in 155 sub-damped Lyman alpha systems (subDLAs) with the aim to investigate the contribution of subDLAs to the chemical evolution of the Universe. The sample was identified within the absorber-blind XQ-100 quasar spectroscopic survey over the redshift range 2.4<=z<=4.3. Using all available column densities of the ionic species investigated (… ▽ More We present the measured gas-phase metal column densities in 155 sub-damped Lyman alpha systems (subDLAs) with the aim to investigate the contribution of subDLAs to the chemical evolution of the Universe. The sample was identified within the absorber-blind XQ-100 quasar spectroscopic survey over the redshift range 2.4<=z<=4.3. Using all available column densities of the ionic species investigated (mainly CIV, SiII, MgII, SiIV, AlII, FeII, CII, and OI; in order of decreasing detection frequency), we estimate the ionization-corrected gas-phase metallicity of each system using Markov Chain Monte Carlo techniques to explore a large grid of Cloudy ionization models. Without accounting for ionization and dust depletion effects, we find that the HI-weighted gas-phase metallicity evolution of subDLAs are consistent with damped Lyman alpha systems (DLAs). When ionization corrections are included, subDLAs are systematically more metal-poor than DLAs (between ~0.5 sigma and ~3 sigma significance) by up to ~1.0 dex over the redshift range 3<=z<=4.3. The correlation of gas-phase [Si/Fe] with metallicity in subDLAs appears to be consistent with that of DLAs, suggesting that the two classes of absorbers have a similar relative dust depletion pattern. As previously seen for Lyman limit systems, the gas-phase [C/O] in subDLAs remains constantly solar for all metallicities indicating that both subDLAs and Lyman limit systems could trace carbon-rich ejecta, potentially in circumgalactic environments. △ Less

Submitted 19 January, 2021; originally announced January 2021.

Comments: Accepted for publication in MNRAS. 64 pages (20 pages of main text, 44 pages of Figures in appendix). Machine-readable versions of Tables 2 and 3 are available in the source files or available online on MNRAS

Probing the thermal state of the intergalactic medium at $z>5$ with the transmission spikes in high-resolution Ly$α$ forest spectra

Authors: Prakash Gaikwad, Michael Rauch, Martin G. Haehnelt, Ewald Puchwein, James S. Bolton, Laura C. Keating, Girish Kulkarni, Vid Iršič, Eduardo Bañados, George D. Becker, Elisa Boera, Fakhri S. Zahedy, Hsiao-Wen Chen, Robert F. Carswell, Jonathan Chardin, Alberto Rorai

Abstract: We compare a sample of five high-resolution, high S/N Ly$α$ forest spectra of bright $6<z \lesssim 6.5$ QSOs aimed at spectrally resolving the last remaining transmission spikes at $z>5$ with those obtained from mock absorption spectra from the Sherwood and Sherwood-Relics suites of hydrodynamical simulations of the intergalactic medium (IGM). We use a profile fitting procedure for the inverted tr… ▽ More We compare a sample of five high-resolution, high S/N Ly$α$ forest spectra of bright $6<z \lesssim 6.5$ QSOs aimed at spectrally resolving the last remaining transmission spikes at $z>5$ with those obtained from mock absorption spectra from the Sherwood and Sherwood-Relics suites of hydrodynamical simulations of the intergalactic medium (IGM). We use a profile fitting procedure for the inverted transmitted flux, $1-F$, similar to the widely used Voigt profile fitting of the transmitted flux $F$ at lower redshifts, to characterise the transmission spikes that probe predominately underdense regions of the IGM. We are able to reproduce the width and height distributions of the transmission spikes, both with optically thin simulations of the post-reionization Universe using a homogeneous UV background and full radiative transfer simulations of a late reionization model. We find that the width of the fitted components of the simulated transmission spikes is very sensitive to the instantaneous temperature of the reionized IGM. The internal structures of the spikes are more prominant in low temeperature models of the IGM. The width distribution of the observed transmission spikes, which require high spectral resolution ($\leq $ 8 km/s) to be resolved, is reproduced for optically thin simulations with a temperature at mean density of $T_0= (11000 \pm 1600,10500\pm 2100,12000 \pm 2200)$ K at $z= (5.4,5.6,5.8)$. This is weakly dependent on the slope of the temperature-density relation, which is favoured to be moderately steeper than isothermal. In the inhomogeneous, late reionization, full radiative transfer simulations where islands of neutral hydrogen persist to $z\sim5.3$, the width distribution of the observed transmission spikes is consistent with the range of $T_0$ caused by spatial fluctuations in the temperature-density relation. △ Less

Submitted 31 March, 2020; v1 submitted 27 January, 2020; originally announced January 2020.

Comments: 20 pages (+10 pages appendices), 12 figures (+14 figures appendices); Accepted in MNRAS; Main results are summarized in Fig. 10, Fig. 12 and Table 3

Sub-damped Lyman alpha systems in the XQ-100 survey I -- Identification and contribution to the cosmological HI budget

Authors: Trystyn A. M. Berg, Sara L. Ellison, Ruben Sanchez-Ramirez, Sebastian Lopez, Valentina D'Odorico, George D. Becker, Lise Christensen, Guido Cupani, Kelly D. Denney, Gabor Worseck

Abstract: Sub-damped Lyman alpha systems (subDLAs; HI column densities of 19.0<=logN(HI)<20.3) are rarely included in the cosmic HI census performed at redshifts z>=1.5, yet are expected to contribute significantly to the overall HI mass budget of the Universe. In this paper, we present a blindly selected sample of 155 subDLAs found along 100 quasar sightlines (with a redshift path length X=475) in the XQ-1… ▽ More Sub-damped Lyman alpha systems (subDLAs; HI column densities of 19.0<=logN(HI)<20.3) are rarely included in the cosmic HI census performed at redshifts z>=1.5, yet are expected to contribute significantly to the overall HI mass budget of the Universe. In this paper, we present a blindly selected sample of 155 subDLAs found along 100 quasar sightlines (with a redshift path length X=475) in the XQ-100 survey to investigate the contribution of subDLAs to the HI mass density of the Universe. The impact of X-Shooter's spectral resolution on sub-DLA identification is evaluated, and found to be sufficient for reliably finding absorbers with logN(HI)>=18.9. We compared the implications of searching for subDLAs solely using HI absorption versus the use of additional metal lines to confirm the identification, and found that metal-selection techniques would have missed 75 subDLAs. Using a bootstrap-Monte Carlo simulation, we computed the column density distribution function (f(N,X)) and the cosmological HI mass density of subDLAs and compared with our previous work based on the XQ-100 damped Lyman alpha systems. We do not find any significant redshift evolution in f(N,X) or cosmological HI mass density for subDLAs. However, subDLAs contribute 10-20 per cent of the total cosmological HI mass density measured at redshifts 2<z<5 (agreeing with previous measurements), and thus have a small but significant contribution to the HI budget of the Universe. △ Less

Submitted 17 July, 2019; originally announced July 2019.

Comments: Accepted for publication in MNRAS. 55 pages (41 pages of appendix containing subDLA catalogue and evolution curves)

The Evolution of O I over 3.2 < z < 6.5: Reionization of the Circumgalactic Medium

Authors: George D. Becker, Max Pettini, Marc Rafelski, Valentina D'Odorico, Elisa Boera, Lise Christensen, Guido Cupani, Sara L. Ellison, Emanuele P. Farina, Michele Fumagalli, Sebastian López, Marcel Neeleman, Emma V. Ryan-Weber, Gábor Worseck

Abstract: We present a survey for metal absorption systems traced by neutral oxygen over $3.2 < z < 6.5$. Our survey uses Keck/ESI and VLT/X-Shooter spectra of 199 QSOs with redshifts up to 6.6. In total we detect 74 O I absorbers, of which 57 are separated from the background QSO by more than 5000 km/s. We use a maximum likelihood approach to fit the distribution of O I $λ$1302 equivalent widths in bins of… ▽ More We present a survey for metal absorption systems traced by neutral oxygen over $3.2 < z < 6.5$. Our survey uses Keck/ESI and VLT/X-Shooter spectra of 199 QSOs with redshifts up to 6.6. In total we detect 74 O I absorbers, of which 57 are separated from the background QSO by more than 5000 km/s. We use a maximum likelihood approach to fit the distribution of O I $λ$1302 equivalent widths in bins of redshift, and from this determine the evolution in number density of absorbers with $W_{1302} > 0.05$ angstroms. We find that the number density does not monotonically increase with decreasing redshift, as would naively be expected from the buildup of metal-enriched circumgalactic gas with time. The number density over $4.9 < z < 5.7$ is a factor of 1.7--4.1 lower (68% confidence) than over $5.7 < z < 6.5$, with a lower value at $z < 5.7$ favored with 99% confidence. This decrease suggests that the fraction of metals in a low-ionization phase is larger at $z \sim 6$ than at lower redshifts. Absorption from highly ionized metals traced by C IV is also weaker in higher-redshift O I systems, supporting this picture. The evolution of O I absorbers implies that metal-enriched circumgalactic gas at $z \sim 6$ is undergoing an ionization transition driven by a strengthening ultraviolet background. This in turn suggests that the reionization of the diffuse intergalactic medium may still be ongoing at or only recently ended by this epoch. △ Less

Submitted 5 July, 2019; originally announced July 2019.

Comments: 43 pages, 85 figures, submitted to ApJ

Studying the Reionization Epoch with QSO Absorption Lines

Authors: George D. Becker, Anson D'Aloisio, Frederick B. Davies, Joseph F. Hennawi, Robert A. Simcoe

Abstract: Absorption signatures in the spectra of QSOs are one of our most powerful tools for studying galactic and intergalactic environments at high redshifts. With the discovery of QSOs out to z > 7, QSO absorption lines are now tracing the end stages of reionization on multiple fronts using the hydrogen Lyman-$α$ forest and heavy element absorbers. Next-generation QSO absorption line studies with large… ▽ More Absorption signatures in the spectra of QSOs are one of our most powerful tools for studying galactic and intergalactic environments at high redshifts. With the discovery of QSOs out to z > 7, QSO absorption lines are now tracing the end stages of reionization on multiple fronts using the hydrogen Lyman-$α$ forest and heavy element absorbers. Next-generation QSO absorption line studies with large optical/IR telescopes will reveal in detail how the first galaxies emerged form the cosmic web, transformed their circum- and inter-galactic environments, and completed the last major phase transition of the Universe. These efforts will complement other upcoming studies of reionization, such as those with JWST, ALMA, and redshifted 21cm experiments. △ Less

Submitted 12 March, 2019; originally announced March 2019.

Comments: Science whitepaper submitted to the Astro2020 Decadal Survey

Three new VHS-DES Quasars at 6.7 < z < 6.9 and Emission Line Properties at z > 6.5

Authors: S. L. Reed, M. Banerji, G. D. Becker, P. C. Hewett, P. Martini, R. G. McMahon, E. Pons, M. Rauch, T. M. C. Abbott, S. Allam, J. Annis, S. Avila, E. Bertin, D. Brooks, E. Buckley-Geer, A. Carnero Rosell, M. Carrasco Kind, J. Carretero, F. J. Castander, C. E. Cunha, C. B. D'Andrea, L. N. da Costa, J. De Vicente, S. Desai, H. T. Diehl , et al. (36 additional authors not shown)

Abstract: We report the results from a search for z > 6.5 quasars using the Dark Energy Survey (DES) Year 3 dataset combined with the VISTA Hemisphere Survey (VHS) and WISE All-Sky Survey. Our photometric selection method is shown to be highly efficient in identifying clean samples of high-redshift quasars leading to spectroscopic confirmation of three new quasars - VDESJ 0244-5008 (z=6.724), VDESJ 0020-365… ▽ More We report the results from a search for z > 6.5 quasars using the Dark Energy Survey (DES) Year 3 dataset combined with the VISTA Hemisphere Survey (VHS) and WISE All-Sky Survey. Our photometric selection method is shown to be highly efficient in identifying clean samples of high-redshift quasars leading to spectroscopic confirmation of three new quasars - VDESJ 0244-5008 (z=6.724), VDESJ 0020-3653 (z=6.834) and VDESJ 0246-5219 (z=6.90) - which were selected as the highest priority candidates in the survey data without any need for additional follow-up observations. The new quasars span the full range in luminosity covered by other z>6.5 quasar samples (J AB = 20.2 to 21.3; M1450 = -25.6 to -26.6). We have obtained spectroscopic observations in the near infrared for VDESJ 0244-5008 and VDESJ 0020-3653 as well as our previously identified quasar, VDESJ 0224-4711 at z=6.50 from Reed et al. (2017). We use the near infrared spectra to derive virial black-hole masses from the full-width-half-maximum of the MgII line. These black-hole masses are ~ 1 - 2 x 10$^9$M$_\odot$. Combining with the bolometric luminosities of these quasars of L$_{\rm{bol}}\simeq$ 1 - 3 x 10$^{47}$implies that the Eddington ratios are high - $\simeq$0.6-1.1. We consider the C\textrm{\textsc{IV}} emission line properties of the sample and demonstrate that our high-redshift quasars do not have unusual C\textrm{\textsc{IV}} line properties when compared to carefully matched low-redshift samples. Our new DES+VHS $z>6.5$ quasars now add to the growing census of luminous, rapidly accreting supermassive black-holes seen well into the epoch of reionisation. △ Less

Submitted 11 January, 2019; originally announced January 2019.

Comments: 12 pages, 9 figures

Revealing reionization with the thermal history of the intergalactic medium: new constraints from the Lyman-$α$ flux power spectrum

Authors: Elisa Boera, George D. Becker, James S. Bolton, Fahad Nasir

Abstract: We present a new investigation of the thermal history of the intergalactic medium (IGM) during and after reionization using the Lyman-$α$ forest flux power spectrum at $4.0\lesssim z\lesssim5.2$. Using a sample of 15 high-resolution spectra, we measure the flux power down to the smallest scales ever probed at these redshifts ($-1\lesssim \log(k/$km$^{-1}$s)$\lesssim -0.7$). These scales are highly… ▽ More We present a new investigation of the thermal history of the intergalactic medium (IGM) during and after reionization using the Lyman-$α$ forest flux power spectrum at $4.0\lesssim z\lesssim5.2$. Using a sample of 15 high-resolution spectra, we measure the flux power down to the smallest scales ever probed at these redshifts ($-1\lesssim \log(k/$km$^{-1}$s)$\lesssim -0.7$). These scales are highly sensitive to both the instantaneous temperature of the IGM and the total energy injected per unit mass during and after reionization. We measure temperatures at the mean density of $T_{0}\sim7000$-8000 K, consistent with no significant temperature evolution for redshifts $4.2\lesssim z\lesssim5.0$. We also present the first observational constraints on the integrated IGM thermal history, finding that the total energy input per unit mass increases from $u_{0}\sim4.6$ ${\rm eV}$ $m_{\rm p}^{-1}$ to 7.3 eV $m_{\rm p}^{-1}$ from $z\sim 6$ to $4.2$ assuming a $Λ$-CDM cosmology. We show how these results can be used simultaneously to obtain information on the timing and the sources of the reionization process. Our first proof of concept using simplistic models of instantaneous reionization produces results comparable to and consistent with the recent Planck constraints, favoring models with $z_{\rm rei}\sim 8.5^{+1.1}_{-0.8}$. △ Less

Submitted 30 January, 2019; v1 submitted 18 September, 2018; originally announced September 2018.

Comments: 36 pages, 35 figures, replaced with the accepted version (ApJ); 2 additional sections in the Appendix, updated references

Evidence for Large-Scale Fluctuations in the Metagalactic Ionizing Background Near Redshift Six

Authors: George D. Becker, Frederick B. Davies, Steven R. Furlanetto, Matthew A. Malkan, Elisa Boera, Craig Douglass

Abstract: The observed scatter in intergalactic Lyman-$α$ opacity at $z \lesssim 6$ requires large-scale fluctuations in the neutral fraction of the intergalactic medium (IGM) after the expected end of reionization. Post-reionization models that explain this scatter invoke fluctuations in either the ionizing ultraviolet background (UVB) or IGM temperature. These models make very different predictions, howev… ▽ More The observed scatter in intergalactic Lyman-$α$ opacity at $z \lesssim 6$ requires large-scale fluctuations in the neutral fraction of the intergalactic medium (IGM) after the expected end of reionization. Post-reionization models that explain this scatter invoke fluctuations in either the ionizing ultraviolet background (UVB) or IGM temperature. These models make very different predictions, however, for the relationship between Lyman-$α$ opacity and local density. Here we test these models using Lyman-$α$ emitting galaxies (LAEs) to trace the density field surrounding the longest and most opaque known Lyman-$α$ trough at $z < 6$. Using deep Subaru Hyper Suprime-Cam narrow-band imaging, we find a highly significant deficit of $z \simeq 5.7$ LAEs within 20 Mpc/$h$ of the trough. The results are consistent with a model in which the scatter in Lyman-$α$ opacity near $z \sim 6$ is driven by large-scale UVB fluctuations, and disfavor a scenario in which the scatter is primarily driven by variations in IGM temperature. UVB fluctuations at this epoch present a boundary condition for reionization models, and may help shed light on the nature of the ionizing sources. △ Less

Submitted 23 March, 2018; originally announced March 2018.

Comments: 13 pages, 9 figure, submitted to ApJ, comments welcome

New constraints on Lyman-α opacity with a sample of 62 quasars at z > 5.7

Authors: Sarah E. I. Bosman, Xiaohui Fan, Linhua Jiang, Sophie L. Reed, Yoshiki Matsuoka, George D. Becker, Martin G. Haehnelt

Abstract: We present measurements of the mean and scatter of the IGM Lyman-α opacity at 4.9 < z < 6.1 along the lines of sight of 62 quasars at z > 5.7, the largest sample assembled at these redshifts to date by a factor of two. The sample size enables us to sample cosmic variance at these redshifts more robustly than ever before. The spectra used here were obtained by the SDSS, DES-VHS and SHELLQs collabor… ▽ More We present measurements of the mean and scatter of the IGM Lyman-α opacity at 4.9 < z < 6.1 along the lines of sight of 62 quasars at z > 5.7, the largest sample assembled at these redshifts to date by a factor of two. The sample size enables us to sample cosmic variance at these redshifts more robustly than ever before. The spectra used here were obtained by the SDSS, DES-VHS and SHELLQs collaborations, drawn from the ESI and X-Shooter archives, reused from previous studies or observed specifically for this work. We measure the effective optical depth of Lyman-α in bins of 10, 30, 50 and 70 cMpc h-1, construct cumulative distribution functions under two treatments of upper limits on flux and explore an empirical analytic fit to residual Lyman-α transmission. We verify the consistency of our results with those of previous studies via bootstrap re-sampling and confirm the existence of tails towards high values in the opacity distributions, which may persist down to z = 5.2. Comparing our results with predictions from cosmological simulations, we find further strong evidence against models that include a spatially uniform ionizing background and temperature-density relation. We also compare to IGM models that include either a fluctuating UVB dominated by rare quasars or temperature fluctuations due to patchy reionization. Although both models produce better agreement with the observations, neither fully captures the observed scatter in IGM opacity. Our sample of 62 z > 5.7 quasar spectra opens many avenues for future study of the reionisation epoch. △ Less

Submitted 24 May, 2018; v1 submitted 22 February, 2018; originally announced February 2018.

Comments: 23 pages, updated to match MNRAS accepted version

A new Measurement of the Intergalactic Temperature at $z\sim 2.55-2.95$

Authors: Alberto Rorai, Robert F. Carswell, Martin G. Haehnelt, George D. Becker, James S. Bolton, Michael T. Murphy

Abstract: We present two measurements of the temperature-density relationship (TDR) of the intergalactic medium (IGM) in the redshift range $2.55 < z < 2.95$ using a sample of 13 high-quality quasar spectra and high resolution numerical simulations of the IGM. Our approach is based on fitting the neutral hydrogen column density $N_{HI}$ and the Doppler parameter $b$ of the absorption lines in the \mlya\ for… ▽ More We present two measurements of the temperature-density relationship (TDR) of the intergalactic medium (IGM) in the redshift range $2.55 < z < 2.95$ using a sample of 13 high-quality quasar spectra and high resolution numerical simulations of the IGM. Our approach is based on fitting the neutral hydrogen column density $N_{HI}$ and the Doppler parameter $b$ of the absorption lines in the \mlya\ forest. The first measurement is obtained using a novel Bayesian scheme which takes into account the statistical correlations between the parameters characterising the lower cut-off of the $b-N_{HI}$ distribution and the power-law parameters $T_0$ and $γ$ describing the TDR. This approach yields $T_0/ 10^3\, {\rm K}=15.6 \pm 4.4 $ and $γ=1.45 \pm 0.17$ independent of the assumed pressure smoothing of the small scale density field. In order to explore the information contained in the overall $b-N_{HI}$ distribution rather than only the lower cut-off, we obtain a second measurement based on a similar Bayesian analysis of the median Doppler parameter for separate column-density ranges of the absorbers. In this case we obtain $T_0/ 10^3\, {\rm K}=14.6 \pm 3.7$ and $γ=1.37 \pm 0.17$ in good agreement with the first measurement. Our Bayesian analysis reveals strong anti-correlations between the inferred $T_0$ and $γ$ for both methods as well as an anti-correlation of the inferred $T_0$ and the pressure smoothing length for the second method, suggesting that the measurement accuracy can in the latter case be substantially increased if independent constraints on the smoothing are obtained. Our results are in good agreement with other recent measurements of the thermal state of the IGM probing similar (over-)density ranges. △ Less

Submitted 2 November, 2017; originally announced November 2017.

Comments: Accepted for publication in MNRAS

Determining the Nature of Late Gunn-Peterson Troughs with Galaxy Surveys

Authors: Frederick B. Davies, George D. Becker, Steven R. Furlanetto

Abstract: Recent observations have discovered long (up to ~110 Mpc/h), opaque Gunn-Peterson troughs in the z ~ 5.5 Lyman-alpha forest, which are challenging to explain with conventional models of the post-reionization intergalactic medium. Here we demonstrate that observations of the galaxy populations in the vicinity of the deepest troughs can distinguish two competing models for these features: deep voids… ▽ More Recent observations have discovered long (up to ~110 Mpc/h), opaque Gunn-Peterson troughs in the z ~ 5.5 Lyman-alpha forest, which are challenging to explain with conventional models of the post-reionization intergalactic medium. Here we demonstrate that observations of the galaxy populations in the vicinity of the deepest troughs can distinguish two competing models for these features: deep voids where the ionizing background is weak due to fluctuations in the mean free path of ionizing photons would show a deficit of galaxies, while residual temperature variations from extended, inhomogeneous reionization would show an overdensity of galaxies. We use large (~550 Mpc/h) semi-numerical simulations of these competing explanations to predict the galaxy populations in the largest of the known troughs at z ~ 5.7. We quantify the strong correlation of Lyman-alpha effective optical depth and galaxy surface density in both models and estimate the degree to which realistic surveys can measure such a correlation. While a spectroscopic galaxy survey is ideal, we also show that a relatively inexpensive narrowband survey of Lyman-alpha-emitting galaxies is ~90% likely to distinguish between the competing models. △ Less

Submitted 29 August, 2017; originally announced August 2017.

Comments: 12 pages, 16 figures. Submitted to ApJ

A deep search for metals near redshift 7: the line-of-sight towards ULAS J1120+0641

Authors: Sarah E. I. Bosman, George D. Becker, Martin G. Haehnelt, Paul C. Hewett, Richard G. McMahon, Daniel J. Mortlock, Chris Simpson, Bram P. Venemans

Abstract: We present a search for metal absorption line systems at the highest redshifts to date using a deep (30h) VLT/X-Shooter spectrum of the z = 7.084 quasi-stellar object (QSO) ULAS J1120+0641. We detect seven intervening systems at z > 5.5, with the highest-redshift system being a C IV absorber at z = 6.51. We find tentative evidence that the mass density of C IV remains flat or declines with redshif… ▽ More We present a search for metal absorption line systems at the highest redshifts to date using a deep (30h) VLT/X-Shooter spectrum of the z = 7.084 quasi-stellar object (QSO) ULAS J1120+0641. We detect seven intervening systems at z > 5.5, with the highest-redshift system being a C IV absorber at z = 6.51. We find tentative evidence that the mass density of C IV remains flat or declines with redshift at z < 6, while the number density of C II systems remains relatively flat over 5 < z < 7. These trends are broadly consistent with models of chemical enrichment by star formation-driven winds that include a softening of the ultraviolet background towards higher redshifts. We find a larger number of weak ( W_rest < 0.3A ) Mg II systems over 5.9 < z < 7.0 than predicted by a power-law fit to the number density of stronger systems. This is consistent with trends in the number density of weak Mg II systems at z = 2.5, and suggests that the mechanisms that create these absorbers are already in place at z = 7. Finally, we investigate the associated narrow Si IV, C IV, and N V absorbers located near the QSO redshift, and find that at least one component shows evidence of partial covering of the continuum source. △ Less

Submitted 24 May, 2017; originally announced May 2017.

Comments: 18 pages, 23 figures, accepted for publication in MNRAS

Journal ref: MNRAS 470 (2): 1919-1934 (2017)

First constraints on fuzzy dark matter from Lyman-$α$ forest data and hydrodynamical simulations

Authors: Vid Iršič, Matteo Viel, Martin G. Haehnelt, James S. Bolton, George D. Becker

Abstract: We present constraints on the masses of extremely light bosons dubbed fuzzy dark matter from Lyman-$α$ forest data. Extremely light bosons with a De Broglie wavelength of $\sim 1$ kpc have been suggested as dark matter candidates that may resolve some of the current small scale problems of the cold dark matter model. For the first time we use hydrodynamical simulations to model the Lyman-$α$ flux… ▽ More We present constraints on the masses of extremely light bosons dubbed fuzzy dark matter from Lyman-$α$ forest data. Extremely light bosons with a De Broglie wavelength of $\sim 1$ kpc have been suggested as dark matter candidates that may resolve some of the current small scale problems of the cold dark matter model. For the first time we use hydrodynamical simulations to model the Lyman-$α$ flux power spectrum in these models and compare with the observed flux power spectrum from two different data sets: the XQ-100 and HIRES/MIKE quasar spectra samples. After marginalization over nuisance and physical parameters and with conservative assumptions for the thermal history of the IGM that allow for jumps in the temperature of up to $5000\rm\,K$, XQ-100 provides a lower limit of 7.1$\times 10^{-22}$ eV, HIRES/MIKE returns a stronger limit of 14.3$\times 10^{-22}$ eV, while the combination of both data sets results in a limit of 20 $\times 10^{-22}$ eV (2$σ$ C.L.). The limits for the analysis of the combined data sets increases to 37.5$\times 10^{-22}$ eV (2$σ$ C.L.) when a smoother thermal history is assumed where the temperature of the IGM evolves as a power-law in redshift. Light boson masses in the range $1-10 \times10^{-22}$ eV are ruled out at high significance by our analysis, casting strong doubts that FDM helps solve the "small scale crisis" of the cold dark matter models. △ Less

Submitted 21 July, 2017; v1 submitted 14 March, 2017; originally announced March 2017.

Comments: 7 pages, 3 figures; supplemental material: 2 pages, 3 figures

Journal ref: Phys. Rev. Lett. 119, 031302 (2017)

Observations of the Lyman series forest towards the redshift 7.1 quasar ULAS J1120+0641

Authors: R. Barnett, S. J. Warren, G. D. Becker, D. J. Mortlock, P. C. Hewett, R. G. McMahon, C. Simpson, B. P. Venemans

Abstract: We present a 30h integration Very Large Telescope X-shooter spectrum of the Lyman series forest towards the $z = 7.084$ quasar ULAS J1120+0641. The only detected transmission at ${\rm S/N}>5$ is confined to seven narrow spikes in the Ly$α$ forest, over the redshift range $5.858<z<6.122$, just longward of the wavelength of the onset of the Ly$β$ forest. There is also a possible detection of one fur… ▽ More We present a 30h integration Very Large Telescope X-shooter spectrum of the Lyman series forest towards the $z = 7.084$ quasar ULAS J1120+0641. The only detected transmission at ${\rm S/N}>5$ is confined to seven narrow spikes in the Ly$α$ forest, over the redshift range $5.858<z<6.122$, just longward of the wavelength of the onset of the Ly$β$ forest. There is also a possible detection of one further unresolved spike in the Ly$β$ forest at $z=6.854$, with ${\rm S/N}=4.5$. We also present revised Hubble Space Telescope F814W photometry of the source. The summed flux from the transmission spikes is in agreement with the F814W photometry, so all the transmission in the Lyman series forest may have been detected. There is a Gunn-Peterson (GP) trough in the Ly$α$ forest from $z=6.122$ all the way to the quasar near zone at $z=7.04$. The trough, of comoving length $240\,h^{-1}$Mpc, is over twice as long as the next longest known GP trough. We combine the spectroscopic and photometric results to constrain the evolution of the Ly$α$ effective optical depth with redshift, $τ_{\rm GP}^{\rm eff}$ extending a similar analysis by Simpson et al. We find $τ_{\rm GP}^{\rm eff} \propto (1+z)^ξ$ where $ξ= 11.2^{+0.4}_{-0.6}$, for $z > 5.5$. The data nevertheless provide only a weak limit on the volume-weighted hydrogen intergalactic (IGM) neutral fraction at $z\sim 6.5$, $x_{\rm HI} > 10^{-4}$, similar to limits at redshift $z\sim6$ from less distant quasars. The new observations cannot extend measurements of the neutral fraction of the IGM to higher values because absorption in the Ly$α$ forest is already saturated near $z\sim6$. For higher neutral fractions, other methods such as measuring the red damping wing of the IGM will be required. △ Less

Submitted 13 February, 2017; originally announced February 2017.

Comments: Accepted for publication in A&A. 11 pages, 10 figures

New Constraints on the free-streaming of warm dark matter from intermediate and small scale Lyman-$α$ forest data

Authors: Vid Iršič, Matteo Viel, Martin G. Haehnelt, James S. Bolton, Stefano Cristiani, George D. Becker, Valentina D'Odorico, Guido Cupani, Tae-Sun Kim, Trystyn A. M. Berg, Sebastian López, Sara Ellison, Lise Christensen, Kelly D. Denny, Gábor Worseck

Abstract: We present new measurements of the free-streaming of warm dark matter (WDM) from Lyman-$α$ flux-power spectra. We use data from the medium resolution, intermediate redshift XQ-100 sample observed with the X-shooter spectrograph ($z=3 - 4.2$) and the high-resolution, high-redshift sample used in Viel et al. (2013) obtained with the HIRES/MIKE spectrographs ($z=4.2 - 5.4$). Based on further improved… ▽ More We present new measurements of the free-streaming of warm dark matter (WDM) from Lyman-$α$ flux-power spectra. We use data from the medium resolution, intermediate redshift XQ-100 sample observed with the X-shooter spectrograph ($z=3 - 4.2$) and the high-resolution, high-redshift sample used in Viel et al. (2013) obtained with the HIRES/MIKE spectrographs ($z=4.2 - 5.4$). Based on further improved modelling of the dependence of the Lyman-$α$ flux-power spectrum on the free-streaming of dark matter, cosmological parameters, as well as the thermal history of the intergalactic medium (IGM) with hydrodynamical simulations, we obtain the following limits, expressed as the equivalent mass of thermal relic WDM particles. The XQ-100 flux power spectrum alone gives a lower limit of 1.4 keV, the re-analysis of the HIRES/MIKE sample gives 4.1 keV while the combined analysis gives our best and significantly strengthened lower limit of 5.3 keV (all 2$σ$ C.L.). The further improvement in the joint analysis is partly due to the fact that the two data sets have different degeneracies between astrophysical and cosmological parameters that are broken when the data sets are combined, and more importantly on chosen priors on the thermal evolution. These results all assume that the temperature evolution of the IGM can be modelled as a power law in redshift. Allowing for a non-smooth evolution of the temperature of the IGM with sudden temperature changes of up to 5000K reduces the lower limit for the combined analysis to 3.5 keV. A WDM with smaller thermal relic masses would require, however, a sudden temperature jump of $5000\,K$ or more in the narrow redshift interval $z=4.6-4.8$, in disagreement with observations of the thermal history based on high-resolution resolution Lyman-$α$ forest data and expectations for photo-heating and cooling in the low density IGM at these redshifts. △ Less

Submitted 8 June, 2018; v1 submitted 6 February, 2017; originally announced February 2017.

Comments: 15 pages, 7 figures; version 3 - corrected a typo in Table II for re-analysis of HIRES/MIKE data with Thermal priors

Journal ref: Phys. Rev. D 96, 023522 (2017)

The Lyman-alpha forest power spectrum from the XQ-100 Legacy Survey

Authors: Vid Iršič, Matteo Viel, Trystyn A. M. Berg, Valentina D'Odorico, Martin G. Haehnelt, Stefano Cristiani, Guido Cupani, Tae-Sun Kim, Sebastian López, Sara Ellison, George D. Becker, Lise Christensen, Kelly D. Denny, Gábor Worseck, James S. Bolton

Abstract: We present the Lyman-$α$ flux power spectrum measurements of the XQ-100 sample of quasar spectra obtained in the context of the European Southern Observatory Large Programme "Quasars and their absorption lines: a legacy survey of the high redshift universe with VLT/XSHOOTER". Using $100$ quasar spectra with medium resolution and signal-to-noise ratio we measure the power spectrum over a range of r… ▽ More We present the Lyman-$α$ flux power spectrum measurements of the XQ-100 sample of quasar spectra obtained in the context of the European Southern Observatory Large Programme "Quasars and their absorption lines: a legacy survey of the high redshift universe with VLT/XSHOOTER". Using $100$ quasar spectra with medium resolution and signal-to-noise ratio we measure the power spectrum over a range of redshifts $z = 3 - 4.2$ and over a range of scales $k = 0.003 - 0.06\,\mathrm{s\,km^{-1}}$. The results agree well with the measurements of the one-dimensional power spectrum found in the literature. The data analysis used in this paper is based on the Fourier transform and has been tested on synthetic data. Systematic and statistical uncertainties of our measurements are estimated, with a total error (statistical and systematic) comparable to the one of the BOSS data in the overlapping range of scales, and smaller by more than $50\%$ for higher redshift bins ($z>3.6$) and small scales ($k > 0.01\,\mathrm{s\,km^{-1}}$). The XQ-100 data set has the unique feature of having signal-to-noise ratios and resolution intermediate between the two data sets that are typically used to perform cosmological studies, i.e. BOSS and high-resolution spectra (e.g. UVES/VLT or HIRES). More importantly, the measured flux power spectra span the high redshift regime which is usually more constraining for structure formation models. △ Less

Submitted 6 February, 2017; originally announced February 2017.

Comments: 17 pages, 12 figures. Accepted in MNRAS

Exploring the Thermal State of the Low-Density Intergalactic Medium at z=3 with an Ultra-High Signal-to-Noise QSO Spectrum

Authors: A. Rorai, G. D. Becker, M. G. Haehnelt, R. F. Carswell, J. S. Bolton, S. Cristiani, V. D'Odorico, G. Cupani, P. Barai, F. Calura, T. -S. Kim, E. Pomante, E. Tescari, M. Viel

Abstract: At low densities the standard ionisation history of the intergalactic medium (IGM) predicts a decreasing temperature of the IGM with decreasing density once hydrogen (and helium) reionisation is complete. Heating the high-redshift, low-density IGM above the temperature expected from photo-heating is difficult, and previous claims of high/rising temperatures in low density regions of the Universe b… ▽ More At low densities the standard ionisation history of the intergalactic medium (IGM) predicts a decreasing temperature of the IGM with decreasing density once hydrogen (and helium) reionisation is complete. Heating the high-redshift, low-density IGM above the temperature expected from photo-heating is difficult, and previous claims of high/rising temperatures in low density regions of the Universe based on the probability density function (PDF) of the opacity in Lyman-$α$ forest data at $2<z<4$ have been met with considerable scepticism, particularly since they appear to be in tension with other constraints on the temperature-density relation (TDR). We utilize here an ultra-high signal-to-noise spectrum of the QSO HE0940-1050 and a novel technique to study the low opacity part of the PDF. We show that there is indeed evidence (at 90% confidence level) that a significant volume fraction of the under-dense regions at $z \sim 3$ has temperatures as high or higher than those at densities comparable to the mean and above. We further demonstrate that this conclusion is nevertheless consistent with measurements of a slope of the TDR in over-dense regions that imply a decreasing temperature with decreasing density, as expected if photo-heating of ionised hydrogen is the dominant heating process. We briefly discuss implications of our findings for the need to invoke either spatial temperature fluctuations, as expected during helium reionization, or additional processes that heat a significant volume fraction of the low-density IGM. △ Less

Submitted 11 November, 2016; originally announced November 2016.

Comments: Accepted for publication on MNRAS