Showing 1–50 of 111 results for author: Davies, F

The Predicament of Absorption-dominated Reionization II: Observational Estimate of the Clumping Factor at the End of Reionization

Authors: Frederick B. Davies, Sarah E. I. Bosman, Steven R. Furlanetto

Abstract: The history of reionization reflects the cumulative injection of ionizing photons by sources and the absorption of ionizing photons by sinks. The latter process is traditionally described in terms of a "clumping factor" which encodes the average quadratic increase in the recombination rate of dense gas within the cosmic web. The recent measurement of a short mean free path of ionizing photons from… ▽ More The history of reionization reflects the cumulative injection of ionizing photons by sources and the absorption of ionizing photons by sinks. The latter process is traditionally described in terms of a "clumping factor" which encodes the average quadratic increase in the recombination rate of dense gas within the cosmic web. The recent measurement of a short mean free path of ionizing photons from stacked quasar spectra at $z\simeq6$ has placed the importance of sinks under increased scrutiny, but its connection to the recombination rate is not immediately obvious. Here we present analytic arguments to connect the clumping factor to the mean free path by invoking ionization equilibrium within the ionized phase of the intergalactic medium at the end of (and after) reionization. We find that the latest mean free path and hydrogen photoionization rate measurements at $z=5-6$ imply a global clumping factor $C\approx12$, much higher than previous determinations from radiation-hydrodynamic simulations of the reionization process. Similar values of $C$ are also derived when applying the same procedure to observations at $2<z<5$. Compared to the traditional assumption of $C=3$, high-redshift galaxies must produce roughly twice as many ionizing photons ($\approx3$ photons per baryon) to reionize the universe by $z\sim6$. This additional requirement on the ionizing photon budget may help to reconcile the reionization history with JWST observations that suggest a far greater output of ionizing photons by the most distant galaxy populations. △ Less

Submitted 11 July, 2024; v1 submitted 26 June, 2024; originally announced June 2024.

Comments: 11 pages, 5+1 figures. Typos fixed, discussion improved, appendix added

Report number: NORDITA 2024-021

Quantifying the Precision of IGM Damping Wing Measurements Towards Quasars

Authors: Timo Kist, Joseph F. Hennawi, Frederick B. Davies

Abstract: We investigate the precision with which the Lyman-$α$ damping wing signature imprinted on the spectra of high-redshift quasars (QSOs) by the foreground neutral intergalactic medium (IGM) can measure the history of cosmic reionization. We leverage a novel inference pipeline based on a generative probabilistic model for the entire spectrum (both red- and blueward of the Lyman-$α$ line), accounting f… ▽ More We investigate the precision with which the Lyman-$α$ damping wing signature imprinted on the spectra of high-redshift quasars (QSOs) by the foreground neutral intergalactic medium (IGM) can measure the history of cosmic reionization. We leverage a novel inference pipeline based on a generative probabilistic model for the entire spectrum (both red- and blueward of the Lyman-$α$ line), accounting for all relevant sources of uncertainty - the stochasticity caused by patchy reionization, the impact of the quasar's ionizing radiation on the IGM, it's unknown intrinsic spectrum, and spectral noise. Performing fast JAX-based Hamiltonian Monte-Carlo (HMC) parameter inference, we precisely measure the underlying global IGM neutral fraction as well as the lifetime of the quasar. Running a battery of tests on over a thousand mocks, we find optimal precision when running the pipeline with a six parameter PCA continuum model (five coefficients and a normalization) on $\mathrm{S}/\mathrm{N} \sim 10$ spectra, binned to a $\sim 500\,\mathrm{km}/\mathrm{s}$ velocity pixel scale, and extending at least out to the C IV $λ\,1549\,\text{Å}$ emission line. After marginalizing out nuisance parameters associated with the quasar continuum, a single spectrum constrains the IGM neutral fraction to $28.0_{-8.8}^{+8.2}\,\%$ and the quasar lifetime to $0.80_{-0.55}^{+0.22}\,\mathrm{dex}$, improving notably towards spectra with a stronger IGM damping wing imprint. Higher precision can be achieved by averaging over statistical quasar samples. We identify two primary sources of uncertainty that contribute approximately equally to the total error budget: the uncertain quasar continuum model and the stochastic distribution of neutral regions arising from both the reionization topology and the location of the quasar's ionization front. △ Less

Submitted 17 June, 2024; originally announced June 2024.

Comments: Submitted to MNRAS

Precisely Measuring the Cosmic Reionization History from IGM Damping Wings Towards Quasars

Authors: Joseph F. Hennawi, Timo Kist, Frederick B. Davies, John Tamanas

Abstract: We introduce a new approach for analyzing the IGM damping wings imprinted on the proximity zones of quasars in the epoch of reionization (EoR). Whereas past work has typically forgone the additional constraining power afforded by the blue side continuum ($1216\,Å \lesssim λ\lesssim 1280\,Å$) and/or opted not to model the large correlated IGM transmission fluctuations in the proximity zone (… ▽ More We introduce a new approach for analyzing the IGM damping wings imprinted on the proximity zones of quasars in the epoch of reionization (EoR). Whereas past work has typically forgone the additional constraining power afforded by the blue side continuum ($1216\,Å \lesssim λ\lesssim 1280\,Å$) and/or opted not to model the large correlated IGM transmission fluctuations in the proximity zone ($λ\lesssim 1216\,Å$), we construct a generative probabilistic model for the entire spectrum accounting for all sources of error - the stochasticity induced by patchy reionization, the impact of the quasar's ionizing radiation on the IGM, the unknown intrinsic spectrum of the quasar, and spectral noise. This principled Bayesian method allows us to marginalize out nuisance parameters associated with the quasar's radiation and its unknown intrinsic spectrum to precisely measure the IGM neutral fraction, $\langle x_{\rm HI}\rangle$. A key element of our analysis is the use of dimensionality reduction (DR) to describe the intrinsic quasar spectrum via a small number of nuisance parameters. Using a large sample of $15,559$ SDSS/BOSS quasars at $z \gtrsim 2.15$ we trained and quantified the performance of six distinct DR methods, and find that a six parameter PCA model (five coefficients plus a normalization) performs best, with complex machine learning approaches providing no advantage. By conducting statistical inference on 100 realistic mock EoR quasar spectra, we demonstrate the reliability of the credibility contours that we obtain on $\langle x_{\rm HI}\rangle$ and the quasar lifetime, $t_{\rm Q}$. The new method introduced here will transform IGM damping wings into a precision probe of reionization, on the same solid methodological and statistical footing as other precision cosmological measurements. △ Less

Submitted 21 June, 2024; v1 submitted 17 June, 2024; originally announced June 2024.

Comments: Submitted to MNRAS

Optical and near-infrared spectroscopy of quasars at $z>6.5$: public data release and composite spectrum

Authors: Silvia Onorato, Joseph F. Hennawi, Jan-Torge Schindler, Jinyi Yang, Feige Wang, Aaron J. Barth, Eduardo Bañados, Anna-Christina Eilers, Sarah E. I. Bosman, Frederick B. Davies, Bram P. Venemans, Chiara Mazzucchelli, Silvia Belladitta, Fabio Vito, Emanuele Paolo Farina, Irham T. Andika, Xiaohui Fan, Fabian Walter, Roberto Decarli, Masafusa Onoue, Riccardo Nanni

Abstract: We present optical and near-infrared (NIR) spectroscopic observations for a sample of $45$ quasars at $6.50 < z \leq 7.64$ with absolute magnitudes at $1450$ Å in the range $-28.82 \leq M_{1450} \leq -24.13$ and their composite spectrum. The median redshift and $M_{1450}$ of the quasars in the sample are $z_{\rm{median}}=6.71$ and $M_{1450,\rm{median}} \simeq -26.1$, respectively. The NIR spectra… ▽ More We present optical and near-infrared (NIR) spectroscopic observations for a sample of $45$ quasars at $6.50 < z \leq 7.64$ with absolute magnitudes at $1450$ Å in the range $-28.82 \leq M_{1450} \leq -24.13$ and their composite spectrum. The median redshift and $M_{1450}$ of the quasars in the sample are $z_{\rm{median}}=6.71$ and $M_{1450,\rm{median}} \simeq -26.1$, respectively. The NIR spectra are taken with echelle spectrographs, complemented with additional data from optical long slit instruments, and then reduced consistently using the open-source Python-based spectroscopic data reduction pipeline PypeIt. The median value of the mean signal-to-noise ratios of the spectra in J, H, and K band (median $\langle \rm{SNR}_λ \rangle$) is: median $\langle \rm{SNR}_{J} \rangle=9.7$, median $\langle \rm{SNR}_{H} \rangle=10.3$, and median $\langle \rm{SNR}_{K} \rangle=11.7$; demonstrating the good data quality. This work presents the largest medium/moderate-resolution sample of quasars at $z>6.5$ from ground-based instruments. Its homogeneity and reproducibility make it ideally suited for several scientific goals, i.e., the study of the quasar proximity zones and damping wings, the Ly$α$ forest, the intergalactic medium's metal content, as well as other properties such as the distribution of SMBH masses and Eddington ratios. Our composite spectrum is compared to others at both high and low-$z$ from the literature, showing differences in the strengths of many emission lines, probably due to differences in luminosity among the samples, but a consistent continuum slope, which proves that the same spectral features are preserved in quasars at different redshift ranges. △ Less

Submitted 11 June, 2024; originally announced June 2024.

Comments: 18 pages, 10 figures, submitted to MNRAS

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 wings in the Lyman-α forest: a model-independent measurement of the neutral fraction at 5.4<z<6.1

Authors: Benedetta Spina, Sarah E. I. Bosman, Frederick B. Davies, Prakash Gaikwad, Yongda Zhu

Abstract: Recent observations have positioned the endpoint of the Epoch of Reionisation (EoR) at redshift $z \sim 5.3$. However, observations of the Lyman-$α$ forest have not yet been able to discern whether reionisation occurred slowly and late, with substantial neutral hydrogen persisting at redshift $\sim 6$, or rapidly and earlier, with the apparent late end driven by the fluctuating UV background. Gunn… ▽ More Recent observations have positioned the endpoint of the Epoch of Reionisation (EoR) at redshift $z \sim 5.3$. However, observations of the Lyman-$α$ forest have not yet been able to discern whether reionisation occurred slowly and late, with substantial neutral hydrogen persisting at redshift $\sim 6$, or rapidly and earlier, with the apparent late end driven by the fluctuating UV background. Gunn-Peterson (GP) absorption troughs are solid indicators that reionisation is not complete until $z=5.3$, but whether they contain significantly neutral gas has not yet been proven. We aim to answer this question by directly measuring, for the first time, the neutral hydrogen fraction ($x_\mathrm{HI}$) at the end of the EoR ($5 \lesssim z \lesssim 6$) in high-redshift quasars spectra. For high neutral fractions $x_\mathrm{HI}\gtrsim0.1$, GP troughs exhibit damping wing (DW) absorption extending over $1000$ km s$^{-1}$ beyond the troughs. While conclusively detected in Lyman-$α$ emission lines of quasars at $z\geq7$, DWs are challenging to observe in the general Lyman-$α$ forest due to absorption complexities and small-scale stochastic transmission features. We report the first successful identification of the stochastic DW signal adjacent to GP troughs at redshifts $z=5.6$ through careful stacking of the dark gaps in Lyman-$α$ forest. We use the signal to present a measurement of the corresponding global $x_\mathrm{HI}=0.19\pm0.07$ $(_{-0.16}^{+0.11})$ at $1σ$ $(2σ)$ at $z=5.6$ and a limit $x_\mathrm{HI}<0.44$ at $z=5.9$. The detection of this signal demonstrates the existence of substantially neutral islands near the conclusion of the EoR, unequivocally signaling a late-and-slow reionization scenario. △ Less

Submitted 20 May, 2024; originally announced May 2024.

Comments: 7 pages, 4 figures, submitted to A&A Letters

A Spatially Resolved [CII] Survey of 31 $z\sim7$ Massive Galaxies Hosting Luminous Quasars

Authors: Feige Wang, Jinyi Yang, Xiaohui Fan, Bram Venemans, Roberto Decarli, Eduardo Bañados, Fabian Walter, Aaron J. Barth, Fuyan Bian, Frederick B. Davies, Anna-Christina Eilers, Emanuele Paolo Farina, Joseph F. Hennawi, Jiang-Tao Li, Chiara Mazzucchelli, Ran Wang, Xue-Bing Wu, Minghao Yue

Abstract: The [CII] 158 $μ$m emission line and the underlying far-infrared (FIR) dust continuum are important tracers for studying star formation and kinematic properties of early galaxies. We present a survey of the [CII] emission lines and FIR continua of 31 luminous quasars at $z>6.5$ using the Atacama Large Millimeter Array (ALMA) and the NOrthern Extended Millimeter Array (NOEMA) at sub-arcsec resoluti… ▽ More The [CII] 158 $μ$m emission line and the underlying far-infrared (FIR) dust continuum are important tracers for studying star formation and kinematic properties of early galaxies. We present a survey of the [CII] emission lines and FIR continua of 31 luminous quasars at $z>6.5$ using the Atacama Large Millimeter Array (ALMA) and the NOrthern Extended Millimeter Array (NOEMA) at sub-arcsec resolution. This survey more than doubles the number of quasars with [CII] and FIR observations at these redshifts and enables statistical studies of quasar host galaxies deep into the epoch of reionization. We detect [CII] emission in 27 quasar hosts with a luminosity range of $L_{\rm [CII]}=(0.3-5.5)\times10^9~L_\odot$ and detect the FIR continuum of 28 quasar hosts with a luminosity range of $L_{\rm FIR}=(0.5-13.0)\times10^{12}~L_\odot$. Both $L_{\rm [CII]}$ and $L_{\rm FIR}$ are correlated ($ρ\simeq0.4$) with the quasar bolometric luminosity, albeit with substantial scatter. The quasar hosts detected by ALMA are clearly resolved with a median diameter of $\sim$5 kpc. About 40% of the quasar host galaxies show a velocity gradient in [CII] emission, while the rest show either dispersion-dominated or disturbed kinematics. Basic estimates of the dynamical masses of the rotation-dominated host galaxies yield $M_{\rm dyn}=(0.1-7.5)\times10^{11}~M_\odot$. Considering our findings alongside those of literature studies, we found that the ratio between $M_{\rm BH}$ and $M_{\rm dyn}$ is about ten times higher than that of local $M_{\rm BH}-M_{\rm dyn}$ relation on average but with substantial scatter (the ratio difference ranging from $\sim$0.6 to 60) and large uncertainties. △ Less

Submitted 23 April, 2024; originally announced April 2024.

Comments: accepted for publication in ApJ

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

Deep Learning the Intergalactic Medium using Lyman-alpha Forest at $ 4 \leq z \leq 5$

Authors: Fahad Nasir, Prakash Gaikwad, Frederick B. Davies, James S. Bolton, Ewald Puchwein, Sarah E. I. Bosman

Abstract: Unveiling the thermal history of the intergalactic medium (IGM) at $4 \leq z \leq 5$ holds the potential to reveal early onset HeII reionization or lingering thermal fluctuations from HI reionization. We set out to reconstruct the IGM gas properties along simulated Lyman-alpha forest data on pixel-by-pixel basis, employing deep Bayesian neural networks. Our approach leverages the Sherwood-Relics s… ▽ More Unveiling the thermal history of the intergalactic medium (IGM) at $4 \leq z \leq 5$ holds the potential to reveal early onset HeII reionization or lingering thermal fluctuations from HI reionization. We set out to reconstruct the IGM gas properties along simulated Lyman-alpha forest data on pixel-by-pixel basis, employing deep Bayesian neural networks. Our approach leverages the Sherwood-Relics simulation suite, consisting of diverse thermal histories, to generate mock spectra. Our convolutional and residual networks with likelihood metric predicts the Ly$α$ optical depth-weighted density or temperature for each pixel in the Ly$α$ forest skewer. We find that our network can successfully reproduce IGM conditions with high fidelity across range of instrumental signal-to-noise. These predictions are subsequently translated into the temperature-density plane, facilitating the derivation of reliable constraints on thermal parameters. This allows us to estimate temperature at mean cosmic density, $T_{\rm 0}$ with one sigma confidence $δT_{\rm 0} \sim 1000{\rm K}$ using only one $20$Mpc/h sightline ($Δz\simeq 0.04$) with a typical reionization history. Existing studies utilize redshift pathlength comparable to $Δz\simeq 4$ for similar constraints. We can also provide more stringent constraints on the slope ($1σ$ confidence interval $δ{\rm γ} \lesssim 0.1$) of the IGM temperature-density relation as compared to other traditional approaches. We test the reconstruction on a single high signal-to-noise observed spectrum ($20$ Mpc/h segment), and recover thermal parameters consistent with current measurements. This machine learning approach has the potential to provide accurate yet robust measurements of IGM thermal history at the redshifts in question. △ Less

Submitted 8 April, 2024; originally announced April 2024.

Comments: 17 pages, Submitted to MNRAS

Blind QSO reconstruction challenge: Exploring methods to reconstruct the Ly$α$ emission line of QSOs

Authors: Bradley Greig, Sarah E. I. Bosman, Frederick B. Davies, Dominika Ďurovčíková, Hassan Fathivavsari, Bin Liu, Romain A. Meyer, Zechang Sun, Valentina D'Odorico, Simona Gallerani, Andrei Mesinger, Yuan-Sen Ting

Abstract: Reconstructing the intrinsic Ly$α$ line flux from high-$z$ QSOs can place constraints on the neutral hydrogen content of the intergalactic medium during reionisation. There are now $\gtrsim10$ different Ly$α$ reconstruction pipelines using different methodologies to predict the Ly$α$ line flux from correlations with the spectral information redward of Ly$α$. However, there have been few attempts t… ▽ More Reconstructing the intrinsic Ly$α$ line flux from high-$z$ QSOs can place constraints on the neutral hydrogen content of the intergalactic medium during reionisation. There are now $\gtrsim10$ different Ly$α$ reconstruction pipelines using different methodologies to predict the Ly$α$ line flux from correlations with the spectral information redward of Ly$α$. However, there have been few attempts to directly compare the performance of these pipelines. Therefore, we devised a blind QSO challenge to compare these reconstruction pipelines on a uniform set of objects. Each author was provided de-identified, observed rest-frame QSO spectra with spectral information only redward of 1260Å rest-frame to ensure unbiased reconstruction. We constructed two samples of 30 QSOs, from X-Shooter and SDSS both spanning $3.5<z<4.5$. Importantly, the purpose of this comparison study was not to champion a single, best performing reconstruction pipeline but rather to explore the relative performance of these pipelines over a range of QSOs with broad observational characteristics to infer general trends. In summary, we find machine learning approaches in general provide the strongest ``best guesses" but underestimate the accompanying statistical uncertainty, although these can be recalibrated, whilst pipelines that decompose the spectral information, for example principal component or factor analysis generally perform better at predicting the Ly$α$ profile. Further, we found that reconstruction pipelines trained on SDSS QSOs performed similarly on average for both the X-Shooter and SDSS samples indicating no discernible biases owing to differences in the observational characteristics of the training set or QSO being reconstructed, although the recovered distributions of reconstructions for X-Shooter were broader likely due to an increased fraction of outliers. △ Less

Submitted 1 April, 2024; originally announced April 2024.

Comments: 30 pages, 16 figures, 4 tables. Submitted to MNRAS, comments welcome

Ultralight vector dark matter search using data from the KAGRA O3GK run

Authors: The LIGO Scientific Collaboration, the Virgo Collaboration, the KAGRA Collaboration, A. G. Abac, R. Abbott, H. Abe, I. Abouelfettouh, F. Acernese, K. Ackley, C. Adamcewicz, S. Adhicary, N. Adhikari, R. X. Adhikari, V. K. Adkins, V. B. Adya, C. Affeldt, D. Agarwal, M. Agathos, O. D. Aguiar, I. Aguilar, L. Aiello, A. Ain, P. Ajith, T. Akutsu, S. Albanesi , et al. (1778 additional authors not shown)

Abstract: Among the various candidates for dark matter (DM), ultralight vector DM can be probed by laser interferometric gravitational wave detectors through the measurement of oscillating length changes in the arm cavities. In this context, KAGRA has a unique feature due to differing compositions of its mirrors, enhancing the signal of vector DM in the length change in the auxiliary channels. Here we prese… ▽ More Among the various candidates for dark matter (DM), ultralight vector DM can be probed by laser interferometric gravitational wave detectors through the measurement of oscillating length changes in the arm cavities. In this context, KAGRA has a unique feature due to differing compositions of its mirrors, enhancing the signal of vector DM in the length change in the auxiliary channels. Here we present the result of a search for $U(1)_{B-L}$ gauge boson DM using the KAGRA data from auxiliary length channels during the first joint observation run together with GEO600. By applying our search pipeline, which takes into account the stochastic nature of ultralight DM, upper bounds on the coupling strength between the $U(1)_{B-L}$ gauge boson and ordinary matter are obtained for a range of DM masses. While our constraints are less stringent than those derived from previous experiments, this study demonstrates the applicability of our method to the lower-mass vector DM search, which is made difficult in this measurement by the short observation time compared to the auto-correlation time scale of DM. △ Less

Submitted 5 March, 2024; originally announced March 2024.

Comments: 20 pages, 5 figures

Report number: LIGO-P2300250

A Massive Protocluster Anchored by a Luminous Quasar at $z=6.63$

Authors: Feige Wang, Jinyi Yang, Joseph F. Hennawi, Xiaohui Fan, Minghao Yue, Eduardo Bañados, Shane Bechtel, Fuyan Bian, Sarah Bosman, Jaclyn B. Champagne, Frederick B. Davies, Roberto Decarli, Emanuele Paolo Farina, Chiara Mazzucchelli, Bram Venemans, Fabian Walter

Abstract: Protoclusters, the progenitors of galaxy clusters, trace large scale structures in the early Universe and are important to our understanding of structure formation and galaxy evolution. To date, only a handful of protoclusters have been identified in the Epoch of Reionization (EoR). As one of the rarest populations in the early Universe, distant quasars that host active supermassive black holes ar… ▽ More Protoclusters, the progenitors of galaxy clusters, trace large scale structures in the early Universe and are important to our understanding of structure formation and galaxy evolution. To date, only a handful of protoclusters have been identified in the Epoch of Reionization (EoR). As one of the rarest populations in the early Universe, distant quasars that host active supermassive black holes are thought to reside in the most massive dark matter halos at that cosmic epoch, and could thus potentially pinpoint some of the earliest protoclusters. In this letter, we report the discovery of a massive protocluster around a luminous quasar at $z=6.63$. This protocluster is anchored by the quasar, and includes three [CII] emitters at $z\sim6.63$, 12 spectroscopically confirmed Ly$α$ emitters (LAEs) at $6.54<z\le6.64$, and a large number of narrow-band imaging selected LAE candidates at the same redshift. This structure has an overall overdensity of $δ=3.3^{+1.1}_{-0.9}$ within $\sim35\times74$ cMpc$^2$ on the sky and an extreme overdensity of $δ>30$ in its central region (i.e., $R\lesssim2$ cMpc). We estimate that this protocluster will collapse into a galaxy cluster with a mass of $6.9^{+1.2}_{-1.4}\times10^{15}~M_\odot$ at the current epoch, more massive than the most massive clusters known in the local Universe such as Coma. In the quasar vicinity, we discover a double-peaked LAE which implies that the quasar has a UV lifetime greater than 0.8 Myrs and has already ionized its surrounding intergalactic medium. △ Less

Submitted 2 February, 2024; originally announced February 2024.

Comments: Accepted for publication in ApJL

Recognizing Blazars Using Radio Morphology from the VLA Sky Survey

Authors: Zhang-Liang Xie, Eduardo Banados, Silvia Belladitta, Chiara Mazzucchelli, Jan-Torge Schindler, Frederick B. Davies, Bram P. Venemans

Abstract: Blazars are radio-loud Active Galactic Nuclei (AGN) whose jets have a very small angle to our line of sight. Observationally, the radio emission are mostly compact or a compact-core with a 1-sided jet. With 2.5$^{\prime\prime}$ resolution at 3 GHz, the Very Large Array Sky Survey (VLASS) enables us to resolve the structure of some blazar candidates in the sky north of Decl. $-40$ deg. We introduce… ▽ More Blazars are radio-loud Active Galactic Nuclei (AGN) whose jets have a very small angle to our line of sight. Observationally, the radio emission are mostly compact or a compact-core with a 1-sided jet. With 2.5$^{\prime\prime}$ resolution at 3 GHz, the Very Large Array Sky Survey (VLASS) enables us to resolve the structure of some blazar candidates in the sky north of Decl. $-40$ deg. We introduce an algorithm to classify radio sources as either blazar-like or non-blazar-like based on their morphology in the VLASS images. We apply our algorithm to three existing catalogs, including one of known blazars (Roma-BzCAT) and two of blazar candidates identified by WISE colors and radio emission (WIBRaLS, KDEBLLACS). We show that in all three catalogs, there are objects with morphology inconsistent with being blazars. Considering all the catalogs, more than 12% of the candidates are unlikely to be blazars, based on this analysis. Notably, we show that 3% of the Roma-BzCAT "confirmed" blazars could be a misclassification based on their VLASS morphology. The resulting table with all sources and their radio morphological classification is available online. △ Less

Submitted 9 January, 2024; v1 submitted 8 January, 2024; originally announced January 2024.

Comments: 21 pages, 11 figures, accepted for publication in The Astrophysical Journal

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

Signatures of Small-scale Structure of the Pre-reionization Intergalactic Medium in $z\gtrsim7$ Quasar Proximity Zones

Authors: Frederick B. Davies, Joseph F. Hennawi

Abstract: The small-scale structure of baryons in the intergalactic medium is intimately linked to their past thermal history. Prior to the $\gtrsim10^4$ K photoheating during the epoch of reionization, cold baryons may have closely traced the clumpy cosmic web of dark matter down to scales as low as $\lesssim1$ comoving kpc, depending on the degree of heating by the X-ray background. After the passage of t… ▽ More The small-scale structure of baryons in the intergalactic medium is intimately linked to their past thermal history. Prior to the $\gtrsim10^4$ K photoheating during the epoch of reionization, cold baryons may have closely traced the clumpy cosmic web of dark matter down to scales as low as $\lesssim1$ comoving kpc, depending on the degree of heating by the X-ray background. After the passage of the ionization front, this clumpy structure can persist for $\sim10^{8}$ years. The strong Ly$α$ damping wings detected towards a few of the highest redshift quasars, in addition to their smaller-than-expected Ly$α$-transmissive proximity zones, suggest that they have ionized and heated the foreground intergalactic medium less than $10^7$ years ago. Signatures of the pre-reionization small-scale structure should thus persist in their intergalactic surroundings. Here we explore how the persistence of this clumpy structure can affect the statistics of Ly$α$ transmission inside the transparent proximity zones of $z\gtrsim7$ quasars by post-processing a suite of small-volume hydrodynamical simulations with 1D ionizing radiative transfer. We find that the Ly$α$ flux power spectrum and flux PDF statistics of ten $z=7.5$ proximity zones, with realistic observational parameters, could distinguish the gaseous structure of a $T_{\rm IGM}\sim2$ K CDM model from warm dark matter models with particle masses $m_{\rm WDM}>10$ keV and X-ray heated models with $f_{\rm X}f_{\rm abs}>0.1$ ($T_{\rm IGM}(z=7.5)\gtrsim275$ K) at the $2σ$ level. △ Less

Submitted 11 December, 2023; originally announced December 2023.

Comments: 18 pages, 17 figures, submitted to ApJ

Possible Contamination of the Intergalactic Medium Damping Wing in ULAS J1342+0928 by Proximate Damped Ly$α$ Absorption

Authors: Frederick B. Davies, Eduardo Bañados, Joseph F. Hennawi, Sarah E. I. Bosman

Abstract: The red damping wing from neutral hydrogen in the intergalactic medium is a smoking-gun signal of ongoing reionization. One potential contaminant of the intergalactic damping wing signal is dense gas associated with foreground galaxies, which can give rise to proximate damped Ly$α$ absorbers. The Ly$α$ imprint of such absorbers on background quasars is indistinguishable from the intergalactic medi… ▽ More The red damping wing from neutral hydrogen in the intergalactic medium is a smoking-gun signal of ongoing reionization. One potential contaminant of the intergalactic damping wing signal is dense gas associated with foreground galaxies, which can give rise to proximate damped Ly$α$ absorbers. The Ly$α$ imprint of such absorbers on background quasars is indistinguishable from the intergalactic medium within the uncertainty of the intrinsic quasar continuum, and their abundance at $z\gtrsim7$ is unknown. Here we show that the complex of low-ionization metal absorption systems recently discovered by deep JWST/NIRSpec observations in the foreground of the $z=7.54$ quasar ULAS~J1342$+$0928 can potentially reproduce the quasar's spectral profile close to rest-frame Ly$α$ without invoking a substantial contribution from the intergalactic medium, but only if the absorbing gas is extremely metal-poor ($[{\rm O}/{\rm H}]\sim-3.5$). Such a low oxygen abundance has never been observed in a damped Ly$α$ absorber at any redshift, but this possibility still complicates the interpretation of the spectrum. Our analysis highlights the need for deep spectroscopy of high-redshift quasars with JWST or ELT to "purify" damping wing quasar samples, an exercise which is impossible for much fainter objects like galaxies. △ Less

Submitted 11 December, 2023; originally announced December 2023.

Comments: 10 pages, 5 figures, submitted to ApJL

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

Authors: Caitlin C. Doughty, Joseph F. Hennawi, Jose Oñorbe, Frederick B. Davies, Zarija Lukić

Abstract: The epoch of hydrogen reionization is complete by $z=5$, but its progression at higher redshifts is uncertain. Measurements of Ly$α$ forest opacity show large scatter at $z<6$, suggestive of spatial fluctuations in neutral fraction ($x_\mathrm{HI}$), temperature, or ionizing background, either individually or in combination. However, these effects are degenerate, necessitating modeling these physi… ▽ More The epoch of hydrogen reionization is complete by $z=5$, but its progression at higher redshifts is uncertain. Measurements of Ly$α$ forest opacity show large scatter at $z<6$, suggestive of spatial fluctuations in neutral fraction ($x_\mathrm{HI}$), temperature, or ionizing background, either individually or in combination. However, these effects are degenerate, necessitating modeling these physics in tandem in order to properly interpret the observations. We begin this process by developing a framework for modeling the reionization history and associated temperature fluctuations, with the intention of incorporating ionizing background fluctuations at a later time. To do this, we generate several reionization histories using semi-numerical code AMBER, selecting histories with volume-weighted neutral fractions that adhere to the observed CMB optical depth and dark pixel fractions. Implementing these histories in the \texttt{Nyx} cosmological hydrodynamics code, we examine the evolution of gas within the simulation, and the associated metrics of the Ly$α$ forest opacity. We find that the pressure smoothing scale within the IGM is strongly correlated with the adiabatic index of the temperature-density relation. We find that while models with 20,000 K photoheating at reionization are better able to reproduce the shape of the observed $z=5$ 1D flux power spectrum than those with 10,000 K, they fail to match the highest wavenumbers. The simulated autocorrelation function and optical depth distributions are systematically low and narrow, respectively, compared to the observed values, but are in better agreement when the reionization history is longer in duration, more symmetric in its distribution of reionization redshifts, or if there are remaining neutral regions at $z<6$. The systematically low variance likely requires the addition of a fluctuating UVB. △ Less

Submitted 3 December, 2023; originally announced December 2023.

Comments: To be submitted to MNRAS. Comments welcome

Searching for the Imprints of AGN Feedback on the Lyman Alpha Forest Around Luminous Red Galaxies

Authors: Vikram Khaire, Teng Hu, Joseph F. Hennawi, Joseph N. Burchett, Michael Walther, Frederick Davies

Abstract: We explore the potential of using the low-redshift Lyman-$α$ (Ly$α$) forest surrounding luminous red galaxies (LRGs) as a tool to constrain active galactic nuclei (AGN) feedback models. Our analysis is based on snapshots from the Illustris and IllustrisTNG simulations at a redshift of $z=0.1$. These simulations offer an ideal platform for studying the influence of AGN feedback on the gas surroundi… ▽ More We explore the potential of using the low-redshift Lyman-$α$ (Ly$α$) forest surrounding luminous red galaxies (LRGs) as a tool to constrain active galactic nuclei (AGN) feedback models. Our analysis is based on snapshots from the Illustris and IllustrisTNG simulations at a redshift of $z=0.1$. These simulations offer an ideal platform for studying the influence of AGN feedback on the gas surrounding galaxies, as they share the same initial conditions and underlying code but incorporate different feedback prescriptions. Both simulations show significant impacts of feedback on the temperature and density of the gas around massive halos. Following our previous work, we adjusted the UV background in both simulations to align with the observed number density of Ly$α$ lines ($\rm dN/dz$) in the intergalactic medium and study the Ly$α$ forest around massive halos hosting LRGs, at impact parameters ($r_{\perp}$) ranging from 0.1 to 100 pMpc. Our findings reveal that $\rm dN/dz$, as a function of $r_{\perp}$, is approximately 1.5 to 2 times higher in IllustrisTNG compared to Illustris up to $r_{\perp}$ of $\sim 10$ pMpc. To further assess whether existing data can effectively discern these differences, we search for archival data containing spectra of background quasars probing foreground LRGs. Through a feasibility analysis based on this data, we demonstrate that ${\rm dN/dz} (r_{\perp})$ measurements can distinguish between feedback models of IllustrisTNG and Illustris with a precision exceeding 12$σ$. This underscores the potential of ${\rm dN/dz} (r_{\perp})$ measurements around LRGs as a valuable benchmark observation for discriminating between different feedback models. △ Less

Submitted 14 November, 2023; originally announced November 2023.

Comments: 21 pages (including 4 page appendix), Submitted to MNRAS

Forecasting constraints on the high-z IGM thermal state from the Lyman-$α$ forest flux auto-correlation function

Authors: Molly Wolfson, Joseph F. Hennawi, Frederick B. Davies, Zarija Lukić, Jose Oñorbe

Abstract: The auto-correlation function of the Lyman-$α$ (Ly$α$) forest flux from high-z quasars can statistically probe all scales of the intergalactic medium (IGM) just after the epoch of reionization. The thermal state of the IGM, which is determined by the physics of reionization, sets the amount of small-scale power seen in the \lya forest. To study the sensitivity of the auto-correlation function to t… ▽ More The auto-correlation function of the Lyman-$α$ (Ly$α$) forest flux from high-z quasars can statistically probe all scales of the intergalactic medium (IGM) just after the epoch of reionization. The thermal state of the IGM, which is determined by the physics of reionization, sets the amount of small-scale power seen in the \lya forest. To study the sensitivity of the auto-correlation function to the thermal state of the IGM, we compute the auto-correlation function from cosmological hydrodynamical simulations with semi-numerical models of the thermal state of the IGM. We create mock data sets of 20 quasars to forecast constraints on $T_0$ and $γ$, which characterize a tight temperature-density relation in the IGM, at $5.4 \leq z \leq 6$. At $z = 5.4$ we find that an ideal data set constrains $T_0$ to 29\% and $γ$ to 9\%. In addition, we investigate four realistic reionization scenarios that combine temperature and ultra-violet background (UVB) fluctuations at $z = 5.8$. We find that, when using mock data generated from a model that includes temperature and UVB fluctuations, we can rule out a model with no temperature or UVB fluctuations at $>1σ$ level 50.5\% of the time. △ Less

Submitted 11 September, 2023; originally announced September 2023.

Comments: 18 pages, 13 figures, comments welcome. arXiv admin note: text overlap with arXiv:2208.09013

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 Impact of the WHIM on the IGM Thermal State Determined from the Low-$z$ Lyman-$α$ Forest

Authors: Teng Hu, Vikram Khaire, Joseph F. Hennawi, Jose Onorbe, Michael Walther, Zarija Lukic, Frederick Davies

Abstract: At $z \lesssim 1$, shock heating caused by large-scale velocity flows and possibly violent feedback from galaxy formation, converts a significant fraction of the cool gas ($T\sim 10^4$ K) in the intergalactic medium (IGM) into warm-hot phase (WHIM) with $T >10^5$K, resulting in a significant deviation from the previously tight power-law IGM temperature-density relationship,… ▽ More At $z \lesssim 1$, shock heating caused by large-scale velocity flows and possibly violent feedback from galaxy formation, converts a significant fraction of the cool gas ($T\sim 10^4$ K) in the intergalactic medium (IGM) into warm-hot phase (WHIM) with $T >10^5$K, resulting in a significant deviation from the previously tight power-law IGM temperature-density relationship, $T=T_0 (ρ/ {\barρ})^{γ-1}$. This study explores the impact of the WHIM on measurements of the low-$z$ IGM thermal state, $[T_0,γ]$, based on the $b$-$N_{H I}$ distribution of the Lyman-$α$ forest. Exploiting a machine learning-enabled simulation-based inference method trained on Nyx hydrodynamical simulations, we demonstrate that [$T_0$, $γ$] can still be reliably measured from the $b$-$N_{H I}$ distribution at $z=0.1$, notwithstanding the substantial WHIM in the IGM. To investigate the effects of different feedback, we apply this inference methodology to mock spectra derived from the IllustrisTNG and Illustris simulations at $z=0.1$. The results suggest that the underlying $[T_0,γ]$ of both simulations can be recovered with biases as low as $|Δ\log(T_0/\text{K})| \lesssim 0.05$ dex, $|Δγ| \lesssim 0.1$, smaller than the precision of a typical measurement. Given the large differences in the volume-weighted WHIM fractions between the three simulations (Illustris 38\%, IllustrisTNG 10\%, Nyx 4\%) we conclude that the $b$-$N_{H I}$ distribution is not sensitive to the WHIM under realistic conditions. Finally, we investigate the physical properties of the detectable Lyman-$α$ absorbers, and discover that although their $T$ and $Δ$ distributions remain mostly unaffected by feedback, they are correlated with the photoionization rate used in the simulation. △ Less

Submitted 28 August, 2023; originally announced August 2023.

Comments: 23 pages, 24 figures

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

First measurement of the Mg II forest correlation function in the Epoch of Reionization

Authors: Suk Sien Tie, Joseph F. Hennawi, Feige Wang, Silvia Onorato, Jinyi Yang, Eduardo Bañados, Frederick B. Davies, Jose Oñorbe

Abstract: In the process of producing the roughly three ionizing photons per atom required to reionize the IGM, the same massive stars explode and eject metals into their surroundings. While the overly sensitive Lya transition makes Gunn-Peterson absorption of background quasar light an ineffective probe of reionization at z > 6, strong low-ionization transitions like the MgII doublet will give rise to a de… ▽ More In the process of producing the roughly three ionizing photons per atom required to reionize the IGM, the same massive stars explode and eject metals into their surroundings. While the overly sensitive Lya transition makes Gunn-Peterson absorption of background quasar light an ineffective probe of reionization at z > 6, strong low-ionization transitions like the MgII doublet will give rise to a detectable "metal-line forest", if metals pollute the neutral IGM. We measure the auto-correlation of the MgII forest transmission using a sample of ten ground based z >= 6.80 quasar spectra probing the redshift range 5.96 < z_MgII < 7.42 (z_MgII,median = 6.47). The correlation function exhibits strong small-scale clustering and a pronounced peak at the doublet velocity (768 km/s) arising from strong absorbers in the CGM of galaxies. After these strong absorbers are identified and masked the signal is consistent with noise. Our measurements are compared to a suite of models generated by combining a large hydrodynamical simulation with a semi-numerical reionization topology, assuming a simple uniform enrichment model. We obtain a 95% credibility upper limit of [Mg/H] < -3.73 at z_MgII,median = 6.47, assuming uninformative priors on [Mg/H] and the IGM neutral fraction x_HI. Splitting the data into low-z (5.96 < z_MgII < 6.47; z_MgII,median = 6.235) and high-z (6.47 < z_MgII < 7.42; z_MgII,median = 6.72) subsamples again yields null-detections and 95% upper limits of [Mg/H] < -3.75 and [Mg/H] < -3.45, respectively. These first measurements set the stage for an approved JWST Cycle 2 program (GO 3526) targeting a similar number of quasars that will be an order of magnitude more sensitive, making the Mgii forest an emerging powerful tool to deliver precision constraints on the reionization and enrichment history of the Universe. △ Less

Submitted 22 August, 2023; originally announced August 2023.

Comments: Submitted to MNRAS; 23 pages, 21 figures, 1 table

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

Can the Low Redshift Lyman Alpha Forest Constrain AGN Feedback Models?

Authors: Vikram Khaire, Teng Hu, Joseph F. Hennawi, Michael Walther, Frederick Davies

Abstract: We investigate the potential of low-redshift Lyman alpha (Ly$α$) forest for constraining active galactic nuclei (AGN) feedback models by analyzing the Illustris and IllustrisTNG simulation at z=0.1. These simulations are ideal for studying the impact of AGN feedback on the intergalactic medium (IGM) as they share initial conditions with significant differences in the feedback prescriptions. Both s… ▽ More We investigate the potential of low-redshift Lyman alpha (Ly$α$) forest for constraining active galactic nuclei (AGN) feedback models by analyzing the Illustris and IllustrisTNG simulation at z=0.1. These simulations are ideal for studying the impact of AGN feedback on the intergalactic medium (IGM) as they share initial conditions with significant differences in the feedback prescriptions. Both simulations reveal that the IGM is significantly impacted by AGN feedback. Specifically, feedback is stronger in Illustris and results in reducing cool baryon fraction to 23% relative to 39% in IllustrisTNG. However, when comparing various statistics of Ly$α$ forest such as 2D and marginalized distributions of Doppler widths and H I column density, line density, and flux power spectrum with real data, we find that most of these statistics are largely insensitive to the differences in feedback models. This lack of sensitivity arises because of the fundamental degeneracy between the fraction of cool baryons and the H I photoionization rate ($Γ_{\rm HI}$) as their product determines the optical depth of the Ly$α$ forest. Since the $Γ_{\rm HI}$ cannot be precisely predicted from first principles, it needs to be treated as a nuisance parameter adjusted to match the observed Ly$α$ line density. After adjusting $Γ_{\rm HI}$, the distinctions in the considered statistics essentially fade away. Only the Ly$α$ flux power spectrum at small spatial scales exhibits potentially observable differences, although this may be specific to the relatively extreme feedback model employed in Illustris. Without independent constraints on either $Γ_{\rm HI}$ or cool baryon fraction, constraining AGN feedback with low-redshift Ly$α$ forest will be very challenging. △ Less

Submitted 21 November, 2023; v1 submitted 8 June, 2023; originally announced June 2023.

Comments: 20 Pages, 11 Figures, Accepted for publication in MNRAS

Convergence of small scale Ly$α$ structure at high-$z$ under different reionization scenarios

Authors: Caitlin C. Doughty, Joseph F. Hennawi, Frederick B. Davies, Zarija Lukić, Jose Oñorbe

Abstract: The Ly$α$ forest (LAF) at $z>5$ probes the thermal and reionization history of the intergalactic medium (IGM) and the nature of dark matter, but its interpretation requires comparison to cosmological hydrodynamical simulations. At high-$z$, convergence of these simulations is more exacting since transmission is dominated by underdense voids that are challenging to resolve. With evidence mounting f… ▽ More The Ly$α$ forest (LAF) at $z>5$ probes the thermal and reionization history of the intergalactic medium (IGM) and the nature of dark matter, but its interpretation requires comparison to cosmological hydrodynamical simulations. At high-$z$, convergence of these simulations is more exacting since transmission is dominated by underdense voids that are challenging to resolve. With evidence mounting for a late end to reionization, small structures down to the sub-kpc level may survive to later times than conventionally thought due to the reduced time for pressure smoothing to impact the gas, further tightening simulation resolution requirements. We perform a suite of simulations using the Eulerian cosmological hydrodynamics code Nyx, spanning domain sizes of 1.25-10 $h^{-1}$ Mpc and 5-80 $h^{-1}$ kpc cells, and explore the interaction of these variables with the timing of reionization on the properties of the matter distribution and the simulated LAF at $z=5.5$. In observable Ly$α$ power, convergence within 10% is achieved for $k< 0.1$ s/km, but larger $k$ shows deviation of up to 20 percent. While a later reionization retains more small structure in the density field, because of the greater thermal broadening there is little difference in the convergence of LAF power between early ($z=9$) and later ($z=6$) reionizations. We conclude that at $z\sim5.5$, resolutions of 10 kpc are necessary for convergence of LAF power at $k<0.1$ s/km, while higher-$k$ modes require higher resolution, and that the timing of reionization does not significantly impact convergence given realistic photoheating. △ Less

Submitted 25 May, 2023; originally announced May 2023.

Comments: Submitted to MNRAS. Comments welcome

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)

A SPectroscopic survey of biased halos In the Reionization Era (ASPIRE): JWST Reveals a Filamentary Structure around a z=6.61 Quasar

Authors: Feige Wang, Jinyi Yang, Joseph F. Hennawi, Xiaohui Fan, Fengwu Sun, Jaclyn B. Champagne, Tiago Costa, Melanie Habouzit, Ryan Endsley, Zihao Li, Xiaojing Lin, Romain A. Meyer, Jan-Torge Schindler, Yunjing Wu, Eduardo Bañados, Aaron J. Barth, Aklant K. Bhowmick, Rebekka Bieri, Laura Blecha, Sarah Bosman, Zheng Cai, Luis Colina, Thomas Connor, Frederick B. Davies, Roberto Decarli , et al. (34 additional authors not shown)

Abstract: We present the first results from the JWST ASPIRE program (A SPectroscopic survey of biased halos In the Reionization Era). This program represents an imaging and spectroscopic survey of 25 reionization-era quasars and their environments by utilizing the unprecedented capabilities of NIRCam Wide Field Slitless Spectroscopy (WFSS) mode. ASPIRE will deliver the largest ($\sim280~{\rm arcmin}^2$) gal… ▽ More We present the first results from the JWST ASPIRE program (A SPectroscopic survey of biased halos In the Reionization Era). This program represents an imaging and spectroscopic survey of 25 reionization-era quasars and their environments by utilizing the unprecedented capabilities of NIRCam Wide Field Slitless Spectroscopy (WFSS) mode. ASPIRE will deliver the largest ($\sim280~{\rm arcmin}^2$) galaxy redshift survey at 3-4 $μ$m among JWST Cycle-1 programs and provide extensive legacy values for studying the formation of the earliest supermassive black holes (SMBHs), the assembly of galaxies, early metal enrichment, and cosmic reionization. In this first ASPIRE paper, we report the discovery of a filamentary structure traced by the luminous quasar J0305-3150 and ten [OIII] emitters at $z=6.6$. This structure has a 3D galaxy overdensity of $δ_{\rm gal}=12.6$ over 637 cMpc$^3$, one of the most overdense structures known in the early universe, and could eventually evolve into a massive galaxy cluster. Together with existing VLT/MUSE and ALMA observations of this field, our JWST observations reveal that J0305-3150 traces a complex environment where both UV-bright and dusty galaxies are present, and indicate that the early evolution of galaxies around the quasar is not simultaneous. In addition, we discovered 31 [OIII] emitters in this field at other redshifts, $5.3<z<6.7$, with half of them situated at $z\sim5.4$ and $z\sim6.2$. This indicates that star-forming galaxies, such as [OIII] emitters, are generally clustered at high redshifts. These discoveries demonstrate the unparalleled redshift survey capabilities of NIRCam WFSS and the potential of the full ASPIRE survey dataset. △ Less

Submitted 19 April, 2023; originally announced April 2023.

Comments: accepted for publication in ApJL

A SPectroscopic survey of biased halos In the Reionization Era (ASPIRE): A First Look at the Rest-frame Optical Spectra of $z > 6.5$ Quasars Using JWST

Authors: Jinyi Yang, Feige Wang, Xiaohui Fan, Joseph F. Hennawi, Aaron J. Barth, Eduardo Bañados, Fengwu Sun, Weizhe Liu, Zheng Cai, Linhua Jiang, Zihao Li, Masafusa Onoue, Jan-Torge Schindler, Yue Shen, Yunjing Wu, Aklant K. Bhowmick, Rebekka Bieri, Laura Blecha, Sarah Bosman, Jaclyn B. Champagne, Luis Colina, Thomas Connor, Tiago Costa, Frederick B. Davies, Roberto Decarli , et al. (31 additional authors not shown)

Abstract: Studies of rest-frame optical emission in quasars at $z>6$ have historically been limited by the wavelengths accessible by ground-based telescopes. The James Webb Space Telescope (JWST) now offers the opportunity to probe this emission deep into the reionization epoch. We report the observations of eight quasars at $z>6.5$ using the JWST/NIRCam Wide Field Slitless Spectroscopy, as a part of the ''… ▽ More Studies of rest-frame optical emission in quasars at $z>6$ have historically been limited by the wavelengths accessible by ground-based telescopes. The James Webb Space Telescope (JWST) now offers the opportunity to probe this emission deep into the reionization epoch. We report the observations of eight quasars at $z>6.5$ using the JWST/NIRCam Wide Field Slitless Spectroscopy, as a part of the ''A SPectroscopic survey of biased halos In the Reionization Era (ASPIRE)" program. Our JWST spectra cover the quasars' emission between rest frame $\sim$ 4100 and 5100 Å. The profiles of these quasars' broad H$β$ emission lines span a FWHM from 3000 to 6000 $\rm{km~s^{-1}}$. The H$β$-based virial black hole (BH) masses, ranging from 0.6 to 2.1 billion solar masses, are generally consistent with their MgII-based BH masses. The new measurements based on the more reliable H$β$ tracer thus confirm the existence of billion solar-mass BHs in the reionization epoch. In the observed [OIII] $λλ$4960,5008 doublets of these luminous quasars, broad components are more common than narrow core components ($\le~1200~\rm{km~s^{-1}}$), and only one quasar shows stronger narrow components than broad. Two quasars exhibit significantly broad and blueshifted [OIII] emission, thought to trace galactic-scale outflows, with median velocities of $-610~\rm{km~s^{-1}}$ and $-1430~\rm{km~s^{-1}}$ relative to the [CII] $158\,μ$m line. All eight quasars show strong optical FeII emission, and follow the Eigenvector 1 relations defined by low-redshift quasars. The entire ASPIRE program will eventually cover 25 quasars and provide a statistical sample for the studies of the BHs and quasar spectral properties. △ Less

Submitted 19 April, 2023; originally announced April 2023.

Comments: 13 pages, 4 figures, accepted for publication in ApJL

Detecting and Characterizing Young Quasars. III. The Impact of Gravitational Lensing Magnification

Authors: Minghao Yue, Anna-Christina Eilers, Robert A. Simcoe, Sirio Belli, Frederick B. Davies, David DePalma, Joseph F. Hennawi, Charlotte A. Mason, Julian B. Muñoz, Erica J. Nelson, Sandro Tacchella

Abstract: We test the impact of gravitational lensing on the lifetime estimates of seven high-redshift quasars at redshift $z\gtrsim6$. The targeted quasars are identified by their small observed proximity zone sizes, which indicate extremely short quasar lifetimes $(t_Q\lesssim10^5 \text{ yrs})$. However, these estimates of quasar lifetimes rely on the assumption that the observed luminosities of the quasa… ▽ More We test the impact of gravitational lensing on the lifetime estimates of seven high-redshift quasars at redshift $z\gtrsim6$. The targeted quasars are identified by their small observed proximity zone sizes, which indicate extremely short quasar lifetimes $(t_Q\lesssim10^5 \text{ yrs})$. However, these estimates of quasar lifetimes rely on the assumption that the observed luminosities of the quasars are intrinsic and not magnified by gravitational lensing, which would bias the lifetime estimates towards younger ages. In order to test possible effects of gravitational lensing, we obtain high-resolution images of the seven quasars with the {\em Hubble Space Telescope (HST)} and look for signs of strong lensing. We do not find any evidence of strong lensing, i.e., all quasars are well-described by point sources, and no foreground lensing galaxy is detected. We estimate that the strong lensing probabilities for these quasars are extremely small $(\sim1.4\times10^{-5})$, and show that weak lensing changes the estimated quasar lifetimes by only $\lesssim0.2$ dex. We thus confirm that the short lifetimes of these quasars are intrinsic. The existence of young quasars indicates a high obscured fraction, radiatively inefficient accretion, and/or flickering light curves for high-redshift quasars. We further discuss the impact of lensing magnification on measurements of black hole masses and Eddington ratios of quasars. △ Less

Submitted 18 April, 2023; originally announced April 2023.

Comments: 14 pages, 4 figures, 3 tables. Accepted by the Astrophysical Journal

On the recent discovery claim of a new $z>7$ quasar

Authors: Sarah E. I. Bosman, Frederick B. Davies, Eduardo Bañados

Abstract: Koptelova et al. 2022 (K22) recently claimed a new quasar discovery at $z=7.46$. After careful consideration of the publicly-available data underlying K22's claim, we find that the observations were contaminated by a moving Solar System object, likely a main-belt asteroid. In the absence of the contaminated photometry, there is no evidence for the nearby, persistent WISE source being a high-redshi… ▽ More Koptelova et al. 2022 (K22) recently claimed a new quasar discovery at $z=7.46$. After careful consideration of the publicly-available data underlying K22's claim, we find that the observations were contaminated by a moving Solar System object, likely a main-belt asteroid. In the absence of the contaminated photometry, there is no evidence for the nearby, persistent WISE source being a high-redshift object; in fact, a detection of the source in DELS $z$-band rules out a redshift $z>7.3$. We present our findings as a cautionary tale of the dangers of passing asteroids for photometric selections. △ Less

Submitted 12 April, 2023; v1 submitted 11 April, 2023; originally announced April 2023.

Comments: RNAAS; 2 pages, 1 figure. The asteroid has been conclusively identified as (43838) 1993 FW49

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

When Spectral Modeling Meets Convolutional Networks: A Method for Discovering Reionization-era Lensed Quasars in Multi-band Imaging Data

Authors: Irham Taufik Andika, Knud Jahnke, Arjen van der Wel, Eduardo Bañados, Sarah E. I. Bosman, Frederick B. Davies, Anna-Christina Eilers, Anton Timur Jaelani, Chiara Mazzucchelli, Masafusa Onoue, Jan-Torge Schindler

Abstract: Over the last two decades, around 300 quasars have been discovered at $z\gtrsim6$, yet only one has identified as being strongly gravitationally lensed. We explore a new approach -- enlarging the permitted spectral parameter space, while introducing a new spatial geometry veto criterion -- which is implemented via image-based deep learning. We first apply this approach to a systematic search for r… ▽ More Over the last two decades, around 300 quasars have been discovered at $z\gtrsim6$, yet only one has identified as being strongly gravitationally lensed. We explore a new approach -- enlarging the permitted spectral parameter space, while introducing a new spatial geometry veto criterion -- which is implemented via image-based deep learning. We first apply this approach to a systematic search for reionization-era lensed quasars, using data from the Dark Energy Survey, the Visible and Infrared Survey Telescope for Astronomy Hemisphere Survey, and the Wide-field Infrared Survey Explorer.Our search method consists of two main parts: (i) the preselection of the candidates based on their spectral energy distributions (SEDs) using catalog-level photometry and (ii) relative probabilities calculation of the candidates being a lens or some contaminant, utilizing a convolutional neural network (CNN) classification. The training data sets are constructed by painting deflected point-source lights over actual galaxy images, to generate realistic galaxy-quasar lens models, optimized to find systems with small image separations, i.e., Einstein radii of $θ_\mathrm{E} \leq 1$ arcsec. Visual inspection is then performed for sources with CNN scores of $P_\mathrm{lens} > 0.1$, which leads us to obtain 36 newly selected lens candidates, which are awaiting spectroscopic confirmation. These findings show that automated SED modeling and deep learning pipelines, supported by modest human input, are a promising route for detecting strong lenses from large catalogs that can overcome the veto limitations of primarily dropout-based SED selection approaches. △ Less

Submitted 6 January, 2023; v1 submitted 26 November, 2022; originally announced November 2022.

Comments: 24 pages, 17 figures, and 2 tables. Accepted for publication in The Astrophysical Journal. We welcome comments from the reader

(Nearly) Model-Independent Constraints on the Neutral Hydrogen Fraction in the Intergalactic Medium at $z\sim 5-7$ Using Dark Pixel Fractions in Ly$α$ and Ly$β$ Forests

Authors: Xiangyu Jin, Jinyi Yang, Xiaohui Fan, Feige Wang, Eduardo Bañados, Fuyan Bian, Frederick B. Davies, Anna-Christina Eilers, Emanuele Paolo Farina, Joseph F. Hennawi, Fabio Pacucci, Bram Venemans, Fabian Walter

Abstract: Cosmic reionization was the last major phase transition of hydrogen from neutral to highly ionized in the intergalactic medium (IGM). Current observations show that the IGM is significantly neutral at $z>7$, and largely ionized by $z\sim5.5$. However, most methods to measure the IGM neutral fraction are highly model-dependent, and are limited to when the volume-averaged neutral fraction of the IGM… ▽ More Cosmic reionization was the last major phase transition of hydrogen from neutral to highly ionized in the intergalactic medium (IGM). Current observations show that the IGM is significantly neutral at $z>7$, and largely ionized by $z\sim5.5$. However, most methods to measure the IGM neutral fraction are highly model-dependent, and are limited to when the volume-averaged neutral fraction of the IGM is either relatively low ($\bar{x}_{\rm HI} \lesssim 10^{-3}$) or close to unity ($\bar{x}_{\rm HI}\sim 1$). In particular, the neutral fraction evolution of the IGM at the critical redshift range of $z=6-7$ is poorly constrained. We present new constraints on $\bar{x}_{\rm HI}$ at $z\sim5.1-6.8$, by analyzing deep optical spectra of $53$ quasars at $5.73<z<7.09$. We derive model-independent upper limits on the neutral hydrogen fraction based on the fraction of "dark" pixels identified in the Lyman $α$ (Ly$α$) and Lyman $β$ (Ly$β$) forests, without any assumptions on the IGM model or the intrinsic shape of the quasar continuum. They are the first model-independent constraints on the IGM neutral hydrogen fraction at $z\sim6.2-6.8$ using quasar absorption measurements. Our results give upper limits of $\bar{x}_{\rm HI}(z=6.3) < 0.79\pm0.04$ (1$σ$), $\bar{x}_{\rm HI} (z=6.5) < 0.87\pm0.03$ (1$σ$), and $\bar{x}_{\rm HI} (z=6.7) < 0.94^{+0.06}_{-0.09}$ (1$σ$). The dark pixel fractions at $z>6.1$ are consistent with the redshift evolution of the neutral fraction of the IGM derived from the Planck 2018. △ Less

Submitted 22 November, 2022; originally announced November 2022.

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

Lyman-alpha Emitters in Ionized Bubbles: Constraining the Environment and Ionized Fraction

Authors: A. C. Trapp, Steven R. Furlanetto, Frederick B. Davies

Abstract: Lyman-alpha emitters (LAEs) are excellent probes of the reionization process, as they must be surrounded by large ionized bubbles in order to be visible during the reionization era. Large ionized regions are thought to correspond to over-dense regions and may be protoclusters, making them interesting test-beds for early massive structures. Close associations containing several LAEs are often assum… ▽ More Lyman-alpha emitters (LAEs) are excellent probes of the reionization process, as they must be surrounded by large ionized bubbles in order to be visible during the reionization era. Large ionized regions are thought to correspond to over-dense regions and may be protoclusters, making them interesting test-beds for early massive structures. Close associations containing several LAEs are often assumed to mark over-dense, ionized bubbles. Here, we develop the first framework to quantify the ionization and density fields of high-z galaxy associations. We explore the interplay between (i) the large-scale density of a survey field, (ii) Poisson noise due to the small number density of bright sources at high redshifts (z~7), and (iii) the effects of the ionized fraction on the observation of LAEs. We use Bayesian statistics, a simple model of reionization, and a Monte-Carlo simulation to construct a more comprehensive method for calculating the large-scale density of LAE regions than previous works. We find that Poisson noise has a strong effect on the inferred density of a region and show how the ionized fraction can be inferred. We then apply our framework to the strongest association yet identified: Hu et al. (2021) found 14 LAEs in a volume of ~50,000 cMpc^3 inside the COSMOS field at z~7. We show that this is most likely a 2.5-sigma over-density inside of an ionized or nearly ionized bubble. We also show that this LAE association implies that the global ionized fraction is Q = 0.60 (+0.08,-0.09), within the context of a simple reionization model. △ Less

Submitted 12 October, 2022; originally announced October 2022.

Comments: 14 pages, 14 figures, submitted to MNRAS

A Generative Model for Quasar Spectra

Authors: Anna-Christina Eilers, David W. Hogg, Bernhard Schölkopf, Daniel Foreman-Mackey, Frederick B. Davies, Jan-Torge Schindler

Abstract: We build a multi-output generative model for quasar spectra and the properties of their black hole engines, based on a Gaussian process latent-variable model. This model treats every quasar as a vector of latent properties such that the spectrum and all physical properties of the quasar are associated with non-linear functions of those latent parameters; the Gaussian process kernel functions defin… ▽ More We build a multi-output generative model for quasar spectra and the properties of their black hole engines, based on a Gaussian process latent-variable model. This model treats every quasar as a vector of latent properties such that the spectrum and all physical properties of the quasar are associated with non-linear functions of those latent parameters; the Gaussian process kernel functions define priors on the function space. Our generative model is trained with a justifiable likelihood function that allows us to treat heteroscedastic noise and missing data correctly, which is crucial for all astrophysical applications. It can predict simultaneously unobserved spectral regions, as well as the physical properties of quasars in held-out test data. We apply the model to rest-frame ultraviolet and optical quasar spectra for which precise black hole masses (based on reverberation mapping measurements) are available. Unlike reverberation-mapping studies, which require multi-epoch data, our model predicts black hole masses from single-epoch spectra, even with limited spectral coverage. We demonstrate the capabilities of the model by predicting black hole masses and unobserved spectral regions. We find that we predict black hole masses at close to the best possible accuracy. △ Less

Submitted 6 September, 2022; originally announced September 2022.

Comments: accepted for publication in ApJ

Forecasting constraints on the mean free path of ionizing photons at $z \geq 5.4$ from the Lyman-$α$ forest flux auto-correlation function

Authors: Molly Wolfson, Joseph F. Hennawi, Frederick B. Davies, Jose Oñorbe

Abstract: Fluctuations in Lyman-$α$ (Ly$α$) forest transmission towards high-$z$ quasars are partially sourced from spatial fluctuations in the ultraviolet background (UVB), the level of which are set by the mean free path of ionizing photons ($λ_{\text{mfp}}$). The auto-correlation function of Ly$α$ forest flux characterizes the strength and scale of transmission fluctuations and, as we show, is thus sensi… ▽ More Fluctuations in Lyman-$α$ (Ly$α$) forest transmission towards high-$z$ quasars are partially sourced from spatial fluctuations in the ultraviolet background (UVB), the level of which are set by the mean free path of ionizing photons ($λ_{\text{mfp}}$). The auto-correlation function of Ly$α$ forest flux characterizes the strength and scale of transmission fluctuations and, as we show, is thus sensitive to $λ_{\text{mfp}}$. Recent measurements at $z \sim 6$ suggest a rapid evolution of $λ_{\text{mfp}}$ at $z>5.0$ which would leave a signature in the evolution of the auto-correlation function. For this forecast, we model mock Ly$α$ forest data with properties similar to the XQR-30 extended data set at $5.4 \leq z \leq 6.0$. At each $z$ we investigate 100 mock data sets and an ideal case where mock data matches model values of the auto-correlation function. For ideal data with $λ_{\text{mfp}}=9.0$ cMpc at $z=6.0$, we recover $λ_{\text{mfp}}=12^{+6}_{-3}$ cMpc. This precision is comparable to direct measurements of $λ_{\text{mfp}}$ from the stacking of quasar spectra beyond the Lyman limit. Hypothetical high-resolution data leads to a $\sim40\%$ reduction in the error bars over all $z$. The distribution of mock values of the auto-correlation function in this work is highly non-Gaussian for high-$z$, which should caution work with other statistics of the high-$z$ Ly$α$ forest against making this assumption. We use a rigorous statistical method to pass an inference test, however future work on non-Gaussian methods will enable higher precision measurements. △ Less

Submitted 3 March, 2023; v1 submitted 18 August, 2022; originally announced August 2022.

Comments: Accepted for publication in MNRAS, 16 pages, 11 figures

Measuring the thermal and ionization state of the low-$z$ IGM using likelihood free inference

Authors: Teng Hu, Vikram Khaire, Joseph F. Hennawi, Michael Walther, Hector Hiss, Justin Alsing, Jose Oñorbe, Zarija Lukic, Frederick Davies

Abstract: We present a new approach to measure the power-law temperature density relationship $T=T_0 (ρ/ \barρ)^{γ-1}$ and the UV background photoionization rate $Γ_{\rm HI}$ of the IGM based on the Voigt profile decomposition of the Ly$α$ forest into a set of discrete absorption lines with Doppler parameter $b$ and the neutral hydrogen column density $N_{\rm HI}$. Previous work demonstrated that the shape… ▽ More We present a new approach to measure the power-law temperature density relationship $T=T_0 (ρ/ \barρ)^{γ-1}$ and the UV background photoionization rate $Γ_{\rm HI}$ of the IGM based on the Voigt profile decomposition of the Ly$α$ forest into a set of discrete absorption lines with Doppler parameter $b$ and the neutral hydrogen column density $N_{\rm HI}$. Previous work demonstrated that the shape of the $b$-$N_{\rm HI}$ distribution is sensitive to the IGM thermal parameters $T_0$ and $γ$, whereas our new inference algorithm also takes into account the normalization of the distribution, i.e. the line-density d$N$/d$z$, and we demonstrate that precise constraints can also be obtained on $Γ_{\rm HI}$. We use density-estimation likelihood-free inference (DELFI) to emulate the dependence of the $b$-$N_{\rm HI}$ distribution on IGM parameters trained on an ensemble of 624 Nyx hydrodynamical simulations at $z = 0.1$, which we combine with a Gaussian process emulator of the normalization. To demonstrate the efficacy of this approach, we generate hundreds of realizations of realistic mock HST/COS datasets, each comprising 34 quasar sightlines, and forward model the noise and resolution to match the real data. We use this large ensemble of mocks to extensively test our inference and empirically demonstrate that our posterior distributions are robust. Our analysis shows that by applying our new approach to existing Ly$α$ forest spectra at $z\simeq 0.1$, one can measure the thermal and ionization state of the IGM with very high precision ($σ_{\log T_0} \sim 0.08$ dex, $σ_γ\sim 0.06$, and $σ_{\log Γ_{\rm HI}} \sim 0.07$ dex). △ Less

Submitted 14 July, 2022; originally announced July 2022.

Comments: 20 pages, 17 figures, accepted for publication in MNRAS

The X-shooter/ALMA Sample of Quasars in the Epoch of Reionization. II. Black Hole Masses, Eddington Ratios, and the Formation of the First Quasars

Authors: Emanuele Paolo Farina, Jan-Torge Schindler, Fabian Walter, Eduardo Bañados, Frederick B. Davies, Roberto Decarli, Anna-Christina Eilers, Xiaohui Fan, Joseph F. Hennawi, Chiara Mazzucchelli, Romain A. Meyer, Benny Trakhtenbrot, Marta Volonteri, Feige Wang, Gabor Worseck, Jinyi Yang, Thales A. Gutcke, Bram P. Venemans, Sarah E. I. Bosman, Tiago Costa, Gisella De Rosa, Alyssa B. Drake, Masafusa Onoue

Abstract: We present measurements of black hole masses and Eddington ratios for a sample of 38 bright (M$_{1450}$ < -24.4 mag) quasars at 5.8 < z < 7.5, derived from VLT/X-shooter near-IR spectroscopy of their broad CIV and MgII emission lines. The black hole masses (on average M$_{BH}$ ~ 4.6 x 10$^9$ M$_{\odot}$) and accretion rates (with Eddington ratios ranging between 0.1 and 1.0) are broadly consistent… ▽ More We present measurements of black hole masses and Eddington ratios for a sample of 38 bright (M$_{1450}$ < -24.4 mag) quasars at 5.8 < z < 7.5, derived from VLT/X-shooter near-IR spectroscopy of their broad CIV and MgII emission lines. The black hole masses (on average M$_{BH}$ ~ 4.6 x 10$^9$ M$_{\odot}$) and accretion rates (with Eddington ratios ranging between 0.1 and 1.0) are broadly consistent with that of similarly luminous 0.3 < z < 2.3 quasars, but there is evidence for a mild increase in the median Eddington ratio going towards z > 6. Combined with deep ALMA observations of the [CII] 158 $μ$m line from the quasar host galaxies and VLT/MUSE investigations of the extended Ly$α$ halos, this study provides fundamental clues to models of the formation and growth of the first massive galaxies and black holes. Compared to local scaling relations, z > 5.7 black holes appear to be over-massive with respect to their host galaxies, and their accretion properties do not change with host galaxy morphology. Under the assumption that the kinematics of the T ~ 10$^4$ K gas, traced by the extended Ly$α$ halos, are dominated by the gravitational potential of the dark matter halo, we report a similar relation between the black hole mass and circular velocity to the one reported for z ~ 0 galaxies. These results paint a picture where the first supermassive black holes reside in massive halos at z > 6 and lead the first stages of galaxy formation by rapidly growing in mass with a duty cycle of order unity. However, this duty cycle needs to drastically drop towards lower redshifts, while the host galaxies continue forming stars at a rate of hundreds of solar masses per year, sustained by the large reservoirs of cool gas surrounding them. △ Less

Submitted 19 September, 2022; v1 submitted 11 July, 2022; originally announced July 2022.

Comments: Accepted for publication on AAS journals. 24 pages + appendices, 15 figures, 5 tables (including full list of z>5.7 quasars with MgII based black hole mass estimates). For access to the data and codes used in this work, please contact the authors

One pot chemical vapor deposition of high optical quality large area monolayer Janus transition metal dichalcogenides

Authors: Ziyang Gan, Ioannis Paradisanos, Ana Estrada-Real, Julian Picker, Emad Najafidehaghani, Francis Davies, Christof Neumann, Cedric Robert, Peter Wiecha, Kenji Watanabe, Takashi Taniguchi, Xavier Marie, Arkady V. Krasheninnikov, Bernhard Urbaszek, Antony George, Andrey Turchanin

Abstract: We report one-pot chemical vapor deposition (CVD) growth of large-area Janus SeMoS monolayers, with the asymmetric top (Se) and bottom (S) chalcogen atomic planes with respect to the central transition metal (Mo) atoms. The formation of these two-dimensional semiconductor monolayers takes place upon the thermodynamic equilibrium-driven exchange of the bottom Se atoms of the initially grown MoSe2 s… ▽ More We report one-pot chemical vapor deposition (CVD) growth of large-area Janus SeMoS monolayers, with the asymmetric top (Se) and bottom (S) chalcogen atomic planes with respect to the central transition metal (Mo) atoms. The formation of these two-dimensional semiconductor monolayers takes place upon the thermodynamic equilibrium-driven exchange of the bottom Se atoms of the initially grown MoSe2 single crystals on gold foils with S atoms. The growth process is characterized by complementary experimental techniques including Raman and X-ray photoelectron spectroscopy and the growth mechanisms are rationalized by first principle calculations. The remarkably high optical quality of the synthesized Janus monolayers is demonstrated by optical and magneto-optical measurements which reveal the strong exciton-phonon coupling and enable to obtain the exciton g-factor of -3.3. △ Less

Submitted 10 May, 2022; originally announced May 2022.

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

Staring at the Shadows of Archaic Galaxies: Damped Ly$α$ and Metal Absorbers toward a Young $z \sim 6$ Weak-line Quasar

Authors: Irham Taufik Andika, Knud Jahnke, Eduardo Bañados, Sarah E. I. Bosman, Frederick B. Davies, Anna-Christina Eilers, Emanuele Paolo Farina, Masafusa Onoue, Arjen van der Wel

Abstract: We characterize the Ly$α$ halo and absorption systems toward PSO J083+11, a unique $z=6.3401$ weak-line quasar, using Gemini/GNIRS, Magellan/FIRE, and VLT/MUSE data. Strong absorptions by hydrogen and several metal lines (e.g., CII, MgII, and OI) are discovered in the spectrum, which indicates the presence of: (i) a proximate sub-damped Ly$α$ (sub-DLA) system at $z=6.314$ and (ii) a MgII absorber… ▽ More We characterize the Ly$α$ halo and absorption systems toward PSO J083+11, a unique $z=6.3401$ weak-line quasar, using Gemini/GNIRS, Magellan/FIRE, and VLT/MUSE data. Strong absorptions by hydrogen and several metal lines (e.g., CII, MgII, and OI) are discovered in the spectrum, which indicates the presence of: (i) a proximate sub-damped Ly$α$ (sub-DLA) system at $z=6.314$ and (ii) a MgII absorber at $z=2.2305$. To describe the observed damping wing signal, we model the Ly$α$ absorption with a combination of a sub-DLA with the neutral hydrogen column density of $\log N_\mathrm{HI} = 20.03 \pm 0.30$ cm$^{-2}$ and absorption from the intergalactic medium with a neutral fraction of around 10 percent. The sub-DLA toward PSO J083+11 has an abundance ratio of [C/O] $=-0.04 \pm 0.33$ and metallicity of [O/H] $=-2.19 \pm 0.44$, similar to those of low-redshift metal-poor DLAs. These measurements suggest that the sub-DLA might truncate PSO J083+11's proximity zone size and complicate the quasar lifetime measurement. However, this quasar shows no sign of a Ly$α$ halo in the MUSE data cube, where the estimated $1σ$ limit of surface brightness is $2.76 \times 10^{-18}$ erg s$^{-1}$ cm$^{-2}$ arcsec$^{-2}$ at aperture size of 1 arcsecond, or equivalent to a Ly$α$ luminosity of $\leq 43.46$ erg s$^{-1}$. This non-detection, while being only weak independent evidence on its own, is at least consistent with a young quasar scenario, as expected for a quasar with a short accretion timescale. △ Less

Submitted 25 April, 2022; v1 submitted 4 April, 2022; originally announced April 2022.

Comments: 22 pages, 12 figures, and 1 table. Replaced to match version accepted for publication in The Astronomical Journal. We welcome the comments from the reader. Related paper: arXiv:2009.07784

Journal ref: The Astronomical Journal, 163:251 (17pp), 2022 June

HeII Ly$α$ Transmission Spikes and Absorption Troughs in Eight High-resolution Spectra Probing the End of HeII Reionization

Authors: Kirill Makan, Gábor Worseck, Frederick B. Davies, Joseph F. Hennawi, J. Xavier Prochaska, Philipp Richter

Abstract: We present statistics of HeII Lya transmission spikes and large-scale absorption troughs using archival high-resolution ($R=λ/Δλ\simeq 12,500$-$18,000$) far-UV spectra of eight HeII-transparent quasars obtained with the Cosmic Origins Spectrograph on the Hubble Space Telescope. The sample covers the redshift range 2.5<z<3.8, thereby probing the rapidly evolving HeII absorption at the end of the He… ▽ More We present statistics of HeII Lya transmission spikes and large-scale absorption troughs using archival high-resolution ($R=λ/Δλ\simeq 12,500$-$18,000$) far-UV spectra of eight HeII-transparent quasars obtained with the Cosmic Origins Spectrograph on the Hubble Space Telescope. The sample covers the redshift range 2.5<z<3.8, thereby probing the rapidly evolving HeII absorption at the end of the HeII reionization epoch. The measured lengths of the troughs decrease dramatically from L>100cMpc at z>3 to L~30cMpc at z~2.7, signaling a significant progression of HeII reionization at these redshifts. Furthermore, unexpectedly long L~65cMpc troughs detected at z~2.9 suggest that the UV background fluctuates at larger scales than predicted by current models. By comparing the measured incidence of transmission spikes to predictions from forward-modeled mock spectra created from the outputs of a (146cMpc)^3 optically thin Nyx hydrodynamical simulation employing different UV background models, we infer the redshift evolution of the HeII photoionization rate $Γ_\mathrm{He\,II}(z)$. The photoionization rate decreases with increasing redshift from $\simeq 4.6\times 10^{-15}\mathrm{\,s^{-1}}$ at z~2.6 to $\simeq 1.2 \times 10^{-15}\mathrm{\,s^{-1}}$ at z~3.2, in agreement with previous inferences from the HeII effective optical depth, and following expected trends of current models of a fluctuating HeII-ionizing background. △ Less

Submitted 9 February, 2022; v1 submitted 14 December, 2021; originally announced December 2021.

Comments: Accepted for publication in ApJ following minor revision (added 2 new figures); 16 pages, 14 figures

IGM damping wing constraints on reionisation from covariance reconstruction of two $z\gtrsim7$ QSOs

Authors: Bradley Greig, Andrei Mesinger, Frederick B. Davies, Feige Wang, Jinyi Yang, Joseph F. Hennawi

Abstract: Bright, high redshift ($z>6$) QSOs are powerful probes of the ionisation state of the intervening intergalactic medium (IGM). The detection of Ly$α$ damping wing absorption imprinted in the spectrum of high-z QSOs can provide strong constraints on the epoch of reionisation (EoR). In this work, we perform an independent Ly$α$ damping wing analysis of two known $z>7$ QSOs; DESJ0252-0503 at $z=7.00$… ▽ More Bright, high redshift ($z>6$) QSOs are powerful probes of the ionisation state of the intervening intergalactic medium (IGM). The detection of Ly$α$ damping wing absorption imprinted in the spectrum of high-z QSOs can provide strong constraints on the epoch of reionisation (EoR). In this work, we perform an independent Ly$α$ damping wing analysis of two known $z>7$ QSOs; DESJ0252-0503 at $z=7.00$ (Wang et al.) and J1007+2115 at $z=7.51$ (Yang et al.). For this, we utilise our existing Bayesian framework which simultaneously accounts for uncertainties in: (i) the intrinsic Ly$α$ emission profile (reconstructed from a covariance matrix of measured emission lines; extended in this work to include NV) and (ii) the distribution of ionised (H\,{\scriptsize II}) regions within the IGM using a $1.6^3$ Gpc$^3$ reionisation simulation. This approach is complementary to that used in the aforementioned works as it focuses solely redward of Ly$α$ ($1218 < λ< 1230$Å) making it more robust to modelling uncertainties while also using a different methodology for (i) and (ii). We find, for a fiducial EoR morphology, $\bar{x}_{\rm HI} = 0.64\substack{+0.19 \\ -0.23}$ (68 per cent) at $z=7$ and $\bar{x}_{\rm HI} = 0.27\substack{+0.21 \\ -0.17}$ at $z=7.51$ consistent within $1σ$ to the previous works above, though both are slightly lower in amplitude. Following the inclusion of NV into our reconstruction pipeline, we perform a reanalysis of ULASJ1120+0641 at $z=7.09$ (Mortlock et al.) and ULASJ1342+0928 at $z=7.54$ (Bañados et al.) finding $\bar{x}_{\rm HI} = 0.44\substack{+0.23 \\ -0.24}$ at $z=7.09$ and $\bar{x}_{\rm HI} = 0.31\substack{+0.18 \\ -0.19}$ at $z=7.54$. Finally, we combine the QSO damping wing constraints for all four $z\gtrsim7$ QSOs to obtain a single, unified constraint of $\bar{x}_{\rm HI} = 0.49\substack{+0.11 \\ -0.11}$ at $z=7.29$. △ Less

Submitted 7 December, 2021; originally announced December 2021.

Comments: 15 pages, 9 figures. Submitted to MNRAS, comments welcome

Improving IGM temperature constraints using wavelet analysis on high-redshift quasars

Authors: Molly Wolfson, Joseph F. Hennawi, Frederick B. Davies, Jose Oñorbe, Hector Hiss, Zarija Lukić

Abstract: The thermal state of the intergalactic medium (IGM) contains vital information about the epoch of reionization, one of the most transformative yet poorly understood periods in the young universe. This thermal state is encoded in the small-scale structure of Lyman-$α$ (Ly$α$) absorption in quasar spectra. The 1D flux power spectrum measures the average small-scale structure along quasar sightlines.… ▽ More The thermal state of the intergalactic medium (IGM) contains vital information about the epoch of reionization, one of the most transformative yet poorly understood periods in the young universe. This thermal state is encoded in the small-scale structure of Lyman-$α$ (Ly$α$) absorption in quasar spectra. The 1D flux power spectrum measures the average small-scale structure along quasar sightlines. At high redshifts, where the opacity is large, averaging mixes high signal-to-noise ratio transmission spikes with noisy absorption troughs. Wavelet amplitudes are an alternate statistic that maintains spatial information while quantifying fluctuations at the same spatial frequencies as the power spectrum, giving them the potential to more sensitively measure the small-scale structure. Previous Ly$α$ forest studies using wavelet amplitude probability density functions (PDFs) used limited spatial frequencies and neglected strong correlations between PDF bins and across wavelets scales, resulting in sub-optimal and unreliable parameter inference. Here we present a novel method for performing statistical inference using wavelet amplitude PDFs that spans the full range of spatial frequencies probed by the power spectrum and that fully accounts for these correlations. We applied this procedure to realistic mock data drawn from a simple thermal model parameterized by the temperature at mean density, $T_0$, and find that wavelets deliver 1$σ$ constraints on $T_0$ that are on average 7% more sensitive at $z=5$ (12% at $z=6$) than those from the power spectrum. We consider the possibility of combing wavelet PDFs with the power, but find that this does not lead to improved sensitivity. △ Less

Submitted 6 October, 2021; originally announced October 2021.

Comments: Accepted for publication in MNRAS, 22 Pages, 16 Figures

Probing Early Super-massive Black Hole Growth and Quasar Evolution with Near-infrared Spectroscopy of 37 Reionization-era Quasars at 6.3 < z <= 7.64

Authors: Jinyi Yang, Feige Wang, Xiaohui Fan, Aaron J. Barth, Joseph F. Hennawi, Riccardo Nanni, Fuyan Bian, Frederick B. Davies, Emanuele P. Farina, Jan-Torge Schindler, Eduardo Banados, Roberto Decarli, Anna-Christina Eilers, Richard Green, Hengxiao Guo, Linhua Jiang, Jiang-Tao Li, Bram Venemans, Fabian Walter, Xue-Bing Wu, Minghao Yue

Abstract: We report the results of near-infrared spectroscopic observations of 37 quasars in the redshift range $6.3< z\le7.64$, including 32 quasars at $z>6.5$, forming the largest quasar near-infrared spectral sample at this redshift. The spectra, taken with Keck, Gemini, VLT, and Magellan, allow investigations of central black hole mass and quasar rest-frame ultraviolet spectral properties. The black hol… ▽ More We report the results of near-infrared spectroscopic observations of 37 quasars in the redshift range $6.3< z\le7.64$, including 32 quasars at $z>6.5$, forming the largest quasar near-infrared spectral sample at this redshift. The spectra, taken with Keck, Gemini, VLT, and Magellan, allow investigations of central black hole mass and quasar rest-frame ultraviolet spectral properties. The black hole masses derived from the MgII emission lines are in the range $(0.3-3.6)\times10^{9}\,M_{\odot}$, which requires massive seed black holes with masses $\gtrsim10^{3-4}\,M_{\odot}$, assuming Eddington accretion since $z=30$. The Eddington ratio distribution peaks at $λ_{\rm Edd}\sim0.8$ and has a mean of 1.08, suggesting high accretion rates for these quasars. The CIV - MgII emission line velocity differences in our sample show an increase of CIV blueshift towards higher redshift, but the evolutionary trend observed from this sample is weaker than the previous results from smaller samples at similar redshift. The FeII/MgII flux ratios derived for these quasars up to $z=7.6$, compared with previous measurements at different redshifts, do not show any evidence of strong redshift evolution, suggesting metal-enriched environments in these quasars. Using this quasar sample, we create a quasar composite spectrum for $z>6.5$ quasars and find no significant redshift evolution of quasar broad emission lines and continuum slope, except for a blueshift of the CIV line. Our sample yields a strong broad absorption line quasar fraction of $\sim$24%, higher than the fractions in lower redshift quasar samples, although this could be affected by small sample statistics and selection effects. △ Less

Submitted 28 September, 2021; originally announced September 2021.

Comments: Accepted for publication in ApJ. 26 pages, 9 figures

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

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