Showing 1–50 of 61 results for author: Bosman, S E I

MIDIS: MIRI uncovers Virgil, an extended source at $z\simeq 6.6$ with the photometric properties of Little Red Dots

Authors: Edoardo Iani, Pierluigi Rinaldi, Karina I. Caputi, Marianna Annunziatella, Danial Langeroodi, Jens Melinder, Pablo G. Pérez-González, Javier Álvarez-Márquez, Leindert A. Boogaard, Sarah E. I. Bosman, Luca Costantin, Thibaud Moutard, Luis Colina, Göran Östlin, Thomas R. Greve, Gillian Wright, Almudena Alonso-Herrero, Arjan Bik, Steven Gillman, Alejandro Crespo Gómez, Jens Hjorth, Alvaro Labiano, John P. Pye, Tuomo V. Tikkanen, Paul P. van der Werf

Abstract: We present Virgil, a MIRI extremely red object (MERO) detected with the F1000W filter as part of the MIRI Deep Imaging Survey (MIDIS) observations of the Hubble Ultra Deep Field (HUDF). Virgil is a Lyman-$α$ emitter (LAE) at $z_{spec} = 6.6312\pm 0.0019$ (from VLT/MUSE) with a rest-frame UV-to-optical spectral energy distribution (SED) typical of LAEs at similar redshifts. However, MIRI observatio… ▽ More We present Virgil, a MIRI extremely red object (MERO) detected with the F1000W filter as part of the MIRI Deep Imaging Survey (MIDIS) observations of the Hubble Ultra Deep Field (HUDF). Virgil is a Lyman-$α$ emitter (LAE) at $z_{spec} = 6.6312\pm 0.0019$ (from VLT/MUSE) with a rest-frame UV-to-optical spectral energy distribution (SED) typical of LAEs at similar redshifts. However, MIRI observations reveal an unexpected extremely red color at rest-frame near-infrared wavelengths, $\rm F444W - F1000W = 2.33 \pm 0.06$. Such steep rise in the near-infrared, completely missed without MIRI imaging, is poorly reproduced by models including only stellar populations and hints towards the presence of an Active Galactic Nucleus (AGN). Interestingly, the overall SED shape of Virgil resembles that of the recently discovered population of Little Red Dots (LRDs) but does not meet their compactness criterion: at rest-frame UV-optical wavelengths Virgil's morphology follows a 2D-Sérsic profile with average index $n = 0.93^{+0.85}_{-0.31}$ and $r_e = 0.43$~pkpc. Only at MIRI wavelengths Virgil is unresolved due to the coarser PSF. We also estimate a bolometric luminosity $L_{\rm bol} = (8.4-11.1)\times 10^{44}\rm~erg~s^{-1}$ and a supermassive black hole mass $M_{\rm BH} = (4-7)\times 10^7\rm ~ M_\odot$ in agreement with recently reported values for LRDs. This discovery demonstrates the crucial importance of deep MIRI surveys to find AGN amongst high-$z$ galaxies that otherwise would be completely missed and raises the question of how common Virgil-like objects could be in the early Universe. △ Less

Submitted 26 June, 2024; originally announced June 2024.

Comments: 17 pages, 10 figures, 3 tables. Submitted to ApJ

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

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

Constraints on the variation of the fine-structure constant at 3<z<10 with JWST emission-line galaxies

Authors: Linhua Jiang, Shuqi Fu, Feige Wang, Sarah E. I. Bosman, Zheng Cai, Hyunsung D. Jun, Zhiwei Pan, Fengwu Sun, Jinyi Yang, Huanian Zhang

Abstract: We present constraints on the spacetime variation of the fine-structure constant $α$ at redshifts $3<z<10$ using JWST emission-line galaxies. The galaxy sample consists of 572 high-quality spectra with strong and narrow [O III] $λλ$4959,5007 doublet emission lines from 522 galaxies, including 267 spectra at $z>5$. The [O III] doublet lines are arguably the best emission lines to probe the variatio… ▽ More We present constraints on the spacetime variation of the fine-structure constant $α$ at redshifts $3<z<10$ using JWST emission-line galaxies. The galaxy sample consists of 572 high-quality spectra with strong and narrow [O III] $λλ$4959,5007 doublet emission lines from 522 galaxies, including 267 spectra at $z>5$. The [O III] doublet lines are arguably the best emission lines to probe the variation in $α$. We divide our sample into 5 subsamples based on redshift and calculate the relative variation $Δα/α$ for the individual subsamples. The calculated $Δα/α$ values are consistent with zero within $1σ$ at all redshifts, suggesting no time variation in $α$ above a level of $(1-2) \times10^{-4}$ ($1σ$) in the past 13.2 billion years. When the whole sample is combined, the constraint is improved to be $Δα/α= (0.4\pm0.7) \times10^{-4}$. We further test the spatial variation in $α$ using four subsamples of galaxies in four different directions on the sky. The measured $Δα/α$ values are consistent with zero at a $1σ$ level of $\sim10^{-4}$. While the constraints in this work are not as stringent as those from lower-redshift quasar absorption lines in previous studies, this work uses an independent tracer and provides the first constraints on $Δα/α$ at the highest redshifts. Our analyses also indicate that the relative wavelength calibration of the JWST spectra is robust. With the growing number of emission-line galaxies from JWST, we expect to achieve stronger constraints in the future. △ Less

Submitted 14 May, 2024; originally announced May 2024.

Comments: 9 pages, 6 figures, submitted to 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

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

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

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

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

Submitted 15 March, 2024; originally announced March 2024.

Comments: 20 pages, 21 figures

A NIRCam-dark galaxy detected with the MIRI/F1000W filter in the MIDIS/JADES Hubble Ultra Deep Field

Authors: Pablo G. Pérez-González, Pierluigi Rinaldi, Karina I. Caputi, Javier Álvarez-Márquez, Marianna Annunziatella, Danial Langeroodi, Thibaud Moutard, Leindert Boogaard, Edoardo Iani, Jens Melinder, Luca Costantin, Goran Östlin, Luis Colina, Thomas R. Greve, Gillian Wright, Almudena Alonso-Herrero, Arjan Bik, Sarah E. I. Bosman, Alejandro Crespo Gómez, Daniel Dicken, Andreas Eckart, Macarena García-Marín, Steven Gillman, Manuel Güdel, Thomas Henning , et al. (10 additional authors not shown)

Abstract: We report the discovery of Cerberus, an extremely red object detected with the MIRI Deep Imaging Survey (MIDIS) observations in the F1000W filter of the Hubble Ultra Deep Field. The object is detected at $S/N\sim6$, with $\mathrm{F1000W}\sim27$ mag, and undetected in the NIRCam data gathered by the JWST Advanced Deep Extragalactic Survey, JADES, fainter than the 30.0-30.5 mag $5σ$ detection limits… ▽ More We report the discovery of Cerberus, an extremely red object detected with the MIRI Deep Imaging Survey (MIDIS) observations in the F1000W filter of the Hubble Ultra Deep Field. The object is detected at $S/N\sim6$, with $\mathrm{F1000W}\sim27$ mag, and undetected in the NIRCam data gathered by the JWST Advanced Deep Extragalactic Survey, JADES, fainter than the 30.0-30.5 mag $5σ$ detection limits in individual bands, as well as in the MIDIS F560W ultra-deep data ($\sim$29 mag, $5σ$). Analyzing the spectral energy distribution built with low-$S/N$ ($<5$) measurements in individual optical-to-mid-infrared filters and higher $S/N$ ($\gtrsim5$) in stacked NIRCam data, we discuss the possible nature of this red NIRCam-dark source using a battery of codes. We discard the possibility of Cerberus being a Solar System body based on the $<0.016$" proper motion in the 1-year apart JADES and MIDIS observations. A sub-stellar Galactic nature is deemed unlikely, given that the Cerberus' relatively flat NIRCam-to-NIRCam and very red NIRCam-to-MIRI flux ratios are not consistent with any brown dwarf model. The extragalactic nature of Cerberus offers 3 possibilities: (1) A $z\sim0.4$ galaxy with strong emission from polycyclic aromatic hydrocarbons; the very low inferred stellar mass, $\mathrm{M}_\star=10^{5-6}$ M$_\odot$, makes this possibility highly improbable. (2) A dusty galaxy at $z\sim4$ with an inferred stellar mass $\mathrm{M}_\star\sim10^{8}$ M$_\odot$. (3) A galaxy with observational properties similar to those of the reddest little red dots discovered around $z\sim7$, but Cerberus lying at $z\sim15$, with the rest-frame optical dominated by emission from a dusty torus or a dusty starburst. △ Less

Submitted 26 May, 2024; v1 submitted 26 February, 2024; originally announced February 2024.

Comments: Accepted for publication in ApJL, results unchanged, manuscript shortened to publish as a letter

A SPectroscopic survey of biased halos In the Reionization Era (ASPIRE): Impact of Galaxies on the CGM Metal Enrichment at z > 6 Using the JWST and VLT

Authors: Siwei Zou, Zheng Cai, Feige Wang, Xiaohui Fan, Jaclyn B. Champagne, Joseph F. Hennawi, Jan-Torge Schindler, Emanuele P. Farina, Jinyi Yang, Kohei Inayoshi, Eduardo Banados, Sarah E. I. Bosman, Zihao Li, Xiaojing Lin, Yunjing Wu, Fengwu Sun, Zi-Yi Guo, Girish Kulkarni, Melanie Habouzit, Stephane Charlot, Jacopo Chevallard, Thomas Connor, Anna-Christina Eilers, Linhua Jiang, Xiangyu Jin , et al. (5 additional authors not shown)

Abstract: We characterize the multiphase circumgalactic medium and galaxy properties at z = 6.0-6.5 in four quasar fields from the James Webb Space Telescope A SPectroscopic survey of biased halos In the Reionization Era (ASPIRE) program. We use the Very Large Telescope/X-shooter spectra of quasar J0305-3150 to identify one new metal absorber at z = 6.2713 with multiple transitions (OI, MgI, FeII and CII).… ▽ More We characterize the multiphase circumgalactic medium and galaxy properties at z = 6.0-6.5 in four quasar fields from the James Webb Space Telescope A SPectroscopic survey of biased halos In the Reionization Era (ASPIRE) program. We use the Very Large Telescope/X-shooter spectra of quasar J0305-3150 to identify one new metal absorber at z = 6.2713 with multiple transitions (OI, MgI, FeII and CII). They are combined with the published absorbing systems in Davies et al. (2023a) at the same redshift range to form of a sample of nine metal absorbers at z = 6.03 to 6.49. We identify eight galaxies within 1000 km s$^{-1}$ and 350 kpc around the absorbing gas from the ASPIRE spectroscopic data, with their redshifts secured by [OIII]($λλ$4959, 5007) doublets and H$β$ emission lines. Our spectral energy distribution fitting indicates that the absorbing galaxies have stellar mass ranging from 10$^{7.2}$ to 10$^{8.8}M_{\odot}$ and metallicity between 0.02 and 0.4 solar. Notably, the z = 6.2713 system in the J0305-3150 field resides in a galaxy overdensity region, which contains two (tentatively) merging galaxies within 350 kpc and seven galaxies within 1 Mpc. We measure the relative abundances of $α$ elements to iron ([$α$/Fe]) and find that the CGM gas in the most overdense region exhibits a lower [$α$/Fe] ratio. Our modeling of the galaxy's chemical abundance favors a top-heavy stellar initial mass function, and hints that we may be witnessing the contribution of the first generation Population III stars to the CGM at the end of reionization epoch. △ Less

Submitted 31 January, 2024; originally announced February 2024.

Comments: 21 pages, 4 figures in the main text. Accepted for publication in ApJL

Widespread AGN feedback in a forming brightest cluster galaxy at $z=4.1$ unveiled by JWST

Authors: Aayush Saxena, Roderik A. Overzier, Montserrat Villar-Martín, Tim Heckman, Namrata Roy, Kenneth J. Duncan, Huub Röttgering, George Miley, Catarina Aydar, Philip Best, Sarah E. I. Bosman, Alex J. Cameron, Krisztina Éva Gabányi, Andrew Humphrey, Sandy Morais, Masafusa Onoue, Laura Pentericci, Victoria Reynaldi, Bram Venemans

Abstract: We present rest-frame optical spectroscopy using JWST/NIRSpec IFU for the radio galaxy TN J1338-1942 at z=4.1, one of the most luminous galaxies in the early Universe with powerful extended radio jets. Previous observations showed evidence for strong, large-scale outflows on the basis of its large (~150 kpc) halo detected in Ly-alpha, and high velocity [O II] emission features detected in ground-b… ▽ More We present rest-frame optical spectroscopy using JWST/NIRSpec IFU for the radio galaxy TN J1338-1942 at z=4.1, one of the most luminous galaxies in the early Universe with powerful extended radio jets. Previous observations showed evidence for strong, large-scale outflows on the basis of its large (~150 kpc) halo detected in Ly-alpha, and high velocity [O II] emission features detected in ground-based IFU data. Our NIRSpec/IFU observations spatially resolve the emission line properties across the host galaxy in great detail. We find at least five concentrations of line emission, coinciding with discrete continuum features previously detected in imaging from HST and JWST, over an extent of ~2'' (~15 kpc). The spectral diagnostics enabled by NIRSpec unambiguously trace the activity of the obscured AGN plus interaction between the interstellar medium and the radio jet as the dominant mechanisms for the ionization state and kinematics of the gas in the system. A secondary region of very high ionization lies at roughly 5 kpc distance from the nucleus, and within the context of an expanding cocoon enveloping the radio lobe, this may be explained by strong shock-ionization of the entrained gas. However, it could also signal the presence of a second obscured AGN, which may also offer an explanation for an intriguing outflow feature seen perpendicular to the radio axis. The presence of a dual SMBH system in this galaxy would support that large galaxies in the early Universe quickly accumulated their mass through the merging of smaller units (each with their own SMBH), at the centers of large overdensities. The inferred black hole mass to stellar mass ratio of 0.01-0.1 for TNJ1338 points to a more rapid assembly of black holes compared to the stellar mass of galaxies at high redshifts, consistent with other recent observations. △ Less

Submitted 22 January, 2024; originally announced January 2024.

Comments: 17 pages, 11 figures, submitted to MNRAS, comments welcome!

JWST Reveals Powerful Feedback from Radio Jets in a Massive Galaxy at z = 4.1

Authors: Namrata Roy, Timothy Heckman, Roderik Overzier, Aayush Saxena, Kenneth Duncan, George Miley, Montserrat Villar Martín, Krisztina Éva Gabányi, Catarina Aydar, Sarah E. I. Bosman, Huub Rottgering, Laura Pentericci, Masafusa Onoue, Victoria Reynaldi

Abstract: We report observations of a powerful ionized gas outflow in a z = 4.1 luminous ($ L_{1.4GHz} \sim 10^{28.3} \ W \ Hz^{-1}$) radio galaxy TNJ1338-1942 hosting an obscured quasar using the Near Infrared Spectrograph (NIRSpec) on board JWST. We spatially resolve a large-scale (~15 kpc) outflow and measure resolved outflow rates. The outflowing gas shows velocities exceeding 900 $ km \ s^{-1}$ and bro… ▽ More We report observations of a powerful ionized gas outflow in a z = 4.1 luminous ($ L_{1.4GHz} \sim 10^{28.3} \ W \ Hz^{-1}$) radio galaxy TNJ1338-1942 hosting an obscured quasar using the Near Infrared Spectrograph (NIRSpec) on board JWST. We spatially resolve a large-scale (~15 kpc) outflow and measure resolved outflow rates. The outflowing gas shows velocities exceeding 900 $ km \ s^{-1}$ and broad line profiles with line widths exceeding 1200 $ km \ s^{-1}$ located at ~10 kpc projected distance from the central nucleus. The outflowing nebula spatially overlaps with the brightest radio lobe, indicating that the powerful radio jets are responsible for the extraordinary kinematics exhibited by the ionized gas. The ionized gas is possibly ionized by the central obscured quasar with a contribution from shocks. The spatially resolved mass outflow rate shows that the region with the broadest line profiles exhibits the strongest outflow rates, with an integrated mass outflow rate of ~500 $ M_{\odot} \ yr^{-1}$. Our hypothesis is that an over-pressured shocked jet fluid expands laterally to create an expanding ellipsoidal "cocoon" that causes the surrounding gas to accelerate outwards. The total kinetic energy injected by the radio jet is about 3 orders of magnitude larger than the total kinetic energy measured in the outflowing ionized gas. This implies that kinetic energy must be transferred inefficiently from the jets to the gas. The bulk of the deposited energy possibly lies in the form of hot (~$ 10^7$ K) X-ray-emitting gas. △ Less

Submitted 21 January, 2024; originally announced January 2024.

Comments: 24 pages, 15 figures, 2 tables. Under review in ApJ

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

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

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

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

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

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

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

MIDIS: Unveiling the Role of Strong Ha-emitters during the Epoch of Reionization with JWST

Authors: P. Rinaldi, K. I. Caputi, E. Iani, L. Costantin, S. Gillman, P. G. Perez-Gonzalez, G. Ostlin, L. Colina, T. R. Greve, H. U. Noorgard-Nielsen, G. S. Wright, J. Alvarez-Marquez, A. Eckart, M. Garcia-Marin, J. Hjorth, O. Ilbert, S. Kendrew, A. Labiano, O. Le Fevre, J. Pye, T. Tikkanen, F. Walter, P. van der Werf, M. Ward, M. Annunziatella , et al. (18 additional authors not shown)

Abstract: By using the ultra-deep \textit{JWST}/MIRI image at 5.6 $μm$ in the Hubble eXtreme Deep Field, we constrain the role of strong H$α$-emitters (HAEs) during Cosmic Reionization at $z\simeq7-8$. Our sample of HAEs is comprised of young ($<35\;\rm Myr$) galaxies, except for one single galaxy ($\approx 300\;\rm Myr$), with low stellar masses ($\lesssim 10^{9}\;\rm M_{\odot}$). These HAEs show a wide ra… ▽ More By using the ultra-deep \textit{JWST}/MIRI image at 5.6 $μm$ in the Hubble eXtreme Deep Field, we constrain the role of strong H$α$-emitters (HAEs) during Cosmic Reionization at $z\simeq7-8$. Our sample of HAEs is comprised of young ($<35\;\rm Myr$) galaxies, except for one single galaxy ($\approx 300\;\rm Myr$), with low stellar masses ($\lesssim 10^{9}\;\rm M_{\odot}$). These HAEs show a wide range of UV-$β$ slopes, with a median value of $β= -2.15\pm0.21$ which broadly correlates with stellar mass. We estimate the ionizing photon production efficiency ($ξ_{ion,0}$) of these sources (assuming $f_{esc,LyC} = 0\%$), which yields a median value $\rm log_{10}(ξ_{ion,0}/(Hz\;erg^{-1})) = 25.50^{+0.10}_{-0.12}$. We show that $ξ_{ion,0}$ positively correlates with EW$_{0}$(H$α$) and specific star formation rate (sSFR). Instead $ξ_{ion,0}$ weakly anti-correlates with stellar mass and $β$. Based on the $β$ values, we predict $f_{esc, LyC}=4\%^{+3}_{-2}$, which results in $\rm log_{10}(ξ_{ion}/(Hz\;erg^{-1})) = 25.55^{+0.11}_{-0.13}$. Considering this and related findings from the literature, we find a mild evolution of $ξ_{ion}$with redshift. Additionally, our results suggest that these HAEs require only modest escape fractions ($f_{esc, rel}$) of 6$-$15\% to reionize their surrounding intergalactic medium. By only considering the contribution of these HAEs, we estimated their total ionizing emissivity ($\dot{N}_{ion}$) as $\dot{N}_{ion} = 10^{50.53 \pm 0.45}; \text{s}^{-1}\text{Mpc}^{-3}$. When comparing their $\dot{N}_{ion}$ with "non-H$α$ emitter" galaxies across the same redshift range, we find that that strong, young, and low-mass emitters may have played an important role during Cosmic Reionization. △ Less

Submitted 13 June, 2024; v1 submitted 27 September, 2023; originally announced September 2023.

Comments: 21 pages, 11 Figures, 1 table. Accepted for publication in ApJ. Updated version 13/06/2024 -- found a mistake in the reported eq. 2 (now corrected)

MIDIS: JWST NIRCam and MIRI unveil the stellar population properties of Ly$α$-emitters and Lyman-Break galaxies at z ~ 3-7

Authors: Edoardo Iani, Karina I. Caputi, Pierluigi Rinaldi, Marianna Annunziatella, Leindert A. Boogaard, Göran Östlin, Luca Costantin, Steven Gillman, Pablo G. Pérez-González, Luis Colina, Gillian Wright, Almudena Alonso-Herrero, Javier Álvarez-Márquez, Arjan Bik, Sarah E. I. Bosman, Alejandro Crespo-Gómez, Andreas Eckart, Thomas R. Greve, Thomas K. Henning, Jens Hjorth, Iris Jermann, Alvaro Labiano, Danial Langeroodi, Jens Melinder, Thibaud Moutard , et al. (7 additional authors not shown)

Abstract: We study the stellar population properties of 182 spectroscopically-confirmed (MUSE/VLT) Lyman-$α$ emitters (LAEs) and 450 photometrically-selected Lyman-Break galaxies (LBGs) at z = 2.8 - 6.7 in the Hubble eXtreme Deep Field (XDF). Leveraging the combined power of HST and JWST NIRCam and MIRI observations, we analyse their rest-frame UV-through-near-IR spectral energy distributions (SEDs) with MI… ▽ More We study the stellar population properties of 182 spectroscopically-confirmed (MUSE/VLT) Lyman-$α$ emitters (LAEs) and 450 photometrically-selected Lyman-Break galaxies (LBGs) at z = 2.8 - 6.7 in the Hubble eXtreme Deep Field (XDF). Leveraging the combined power of HST and JWST NIRCam and MIRI observations, we analyse their rest-frame UV-through-near-IR spectral energy distributions (SEDs) with MIRI playing a crucial role in robustly assessing the LAE's stellar mass and ages. Our LAEs are low-mass objects (log$_{10}$(M$_\star$[M$_\odot$]) ~ 7.5), with little or no dust extinction (E(B - V) ~ 0.1) and a blue UV continuum slope ($β$ ~ -2.2). While 75% of our LAEs are young (< 100 Myr), the remaining 25% have significantly older stellar populations (> 100 Myr). These old LAEs are statistically more massive, less extinct and have lower specific star formation rate (sSFR) compared to young LAEs. Besides, they populate the M$_\star$ - SFR plane along the main-sequence (MS) of star-forming galaxies, while young LAEs populate the starburst region. The comparison between the LAEs properties to those of a stellar-mass matched sample of LBGs shows no statistical difference between these objects, except for the LBGs redder UV continuum slope and marginally larger E(B - V) values. Interestingly, 48% of the LBGs have ages < 10 Myr and are classified as starbursts, but lack detectable Ly$α$ emission. This is likely due to HI resonant scattering and/or selective dust extinction. Overall, we find that JWST observations are crucial in determining the properties of LAEs and shedding light on the properties and similarities between LAEs and LBGs. △ Less

Submitted 25 January, 2024; v1 submitted 15 September, 2023; originally announced September 2023.

Comments: 22 pages, 15 figures, 5 tables. Accepted for publication in APJ

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

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

First rest-frame infrared spectrum of a z>7 quasar: JWST/MRS observations of J1120+0641

Authors: Sarah E. I. Bosman, Javier Álvarez-Márquez, Luis Colina, Fabian Walter, Almudena Alonso-Herrero, Martin J. Ward, Göran Östlin, Thomas R. Greve, Gillian Wright, Arjan Bik, Leindert Boogaard, Karina I. Caputi, Luca Costantin, Andreas Eckart, Macarena García-Marín, Steven Gillman, Manuel Güdel, Thomas Henning, Jens Hjorth, Edoardo Iani, Olivier Ilbert, Iris Jermann, Alvaro Labiano, Pierre-Olivier Lagage, Danial Langeroodi , et al. (7 additional authors not shown)

Abstract: We present a JWST/MRS spectrum of the quasar J1120+0641 at z=7.0848, the first spectroscopic observation of a reionisation-era quasar in the rest-frame infrared ($0.6<λ<3.4μ$m). In the context of the mysterious fast assembly of the first supermassive black holes at z>7, our observations enable for the first time the detection of hot torus dust, the H$α$ emission line, and the Paschen-series broad… ▽ More We present a JWST/MRS spectrum of the quasar J1120+0641 at z=7.0848, the first spectroscopic observation of a reionisation-era quasar in the rest-frame infrared ($0.6<λ<3.4μ$m). In the context of the mysterious fast assembly of the first supermassive black holes at z>7, our observations enable for the first time the detection of hot torus dust, the H$α$ emission line, and the Paschen-series broad emission lines in a quasar at z>7. Hot torus dust is clearly detected as an upturn in the continuum emission at $λ_{\text{rest}}\simeq1.3μ$m, leading to a black-body temperature of $T=1413.5^{+5.7}_{-7.4}$K. Compared to similarly-luminous quasars at 0<z<6, the hot dust in J1120+0641 is somewhat elevated in temperature (top 1%). The temperature is more typical among 6<z<6.5 quasars (top 25%), leading us to postulate a weak evolution in the hot dust temperature at z>6 ($2σ$ significance). We measure the black hole mass of J1120+0641 based on the H$α$ Balmer line, $M_{\text{BH}}=1.52\pm0.17\cdot 10^9 M_\odot$, which is in good agreement with the previous rest-UV MgII black hole mass measurement. The black hole mass based on the Paschen-series lines is also consistent, indicating no significant extinction in the rest-frame UV measurement. The broad H$α$, Pa-$α$ and Pa-$β$ emission lines are consistent with an origin in a common broad-line region (BLR) with density log$N_H/\text{cm}^{-3}\geq 12$, ionisation parameter $-7<$log$U<-4$, and extinction E(B-V)$\lesssim 0.1$mag. These BLR parameters are consistent with similarly-bright quasars at 0<z<4. Overall, we find that both J1120+0641's hot dust torus and hydrogen BLR properties show no significant peculiarity when compared to luminous quasars down to z=0. The quasar accretion structures must have therefore assembled very quickly, as they appear fully "mature" less than 760 million years after the Big Bang. △ Less

Submitted 26 July, 2023; originally announced July 2023.

Comments: 19 pages, 10 figures; submitted

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)

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

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

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

Submitted 1 May, 2023; originally announced May 2023.

Comments: 16 pages, 9 figures, Accepted in MNRAS

Report number: TIFR/TH/23-5

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

Life beyond 30: probing the -20<M_UV<-17 luminosity function at 8<z<13 with the NIRCam parallel field of the MIRI Deep Survey

Authors: Pablo G. Pérez-González, Luca Costantin, Danial Langeroodi, Pierluigi Rinaldi, Marianna Annunziatella, Olivier Ilbert, Luis Colina, Hans Ulrik Noorgaard-Nielsen, Thomas Greve, Göran Ostlin, Gillian Wright, Almudena Alonso-Herrero, Javier Álvarez-Márquez, Karina I. Caputi, Andreas Eckart, Olivier Le Fèvre, Álvaro Labiano, Macarena García-Marín, Jens Hjorth, Sarah Kendrew, John P. Pye, Tuomo Tikkanen, Paul van der Werf, Fabian Walter, Martin Ward , et al. (19 additional authors not shown)

Abstract: We present the ultraviolet luminosity function and an estimate of the cosmic star formation rate density at $8<z<13$ derived from deep NIRCam observations taken in parallel with the MIRI Deep Survey (MDS) of the Hubble Ultra Deep Field (HUDF), NIRCam covering the parallel field 2 (HUDF-P2). Our deep (40 hours) NIRCam observations reach a F277W magnitude of 30.8 ($5σ$), more than 2 magnitudes deepe… ▽ More We present the ultraviolet luminosity function and an estimate of the cosmic star formation rate density at $8<z<13$ derived from deep NIRCam observations taken in parallel with the MIRI Deep Survey (MDS) of the Hubble Ultra Deep Field (HUDF), NIRCam covering the parallel field 2 (HUDF-P2). Our deep (40 hours) NIRCam observations reach a F277W magnitude of 30.8 ($5σ$), more than 2 magnitudes deeper than JWST public datasets already analyzed to find high redshift galaxies. We select a sample of 44 $z>8$ galaxy candidates based on their dropout nature in the F115W and/or F150W filters, a high probability for their photometric redshifts, estimated with three different codes, being at $z>8$, good fits based on $χ^2$ calculations, and predominant solutions compared to $z<8$ alternatives. We find mild evolution in the luminosity function from $z\sim13$ to $z\sim8$, i.e., only a small increase in the average number density of $\sim$0.2 dex, while the faint-end slope and absolute magnitude of the knee remain approximately constant, with values $α=-2.2\pm0.1$ and $M^*=-20.8\pm0.2$ mag. Comparing our results with the predictions of state-of-the-art galaxy evolution models, we find two main results: (1) a slower increase with time in the cosmic star formation rate density compared to a steeper rise predicted by models; (2) nearly a factor of 10 higher star formation activity concentrated in scales around 2 kpc in galaxies with stellar masses $\sim10^8$ M$_\odot$ during the first 350 Myr of the Universe, $z\sim12$, with models matching better the luminosity density observational estimations $\sim$150 Myr later, by $z\sim9$. △ Less

Submitted 22 May, 2023; v1 submitted 5 February, 2023; originally announced February 2023.

Comments: Accepted for publication in ApJL, in memoriam Hans Ulrik Noorgaard-Nielsen and Olivier Le Fèvre

The fraction and kinematics of broad absorption line quasars across cosmic time

Authors: Manuela Bischetti, Fabrizio Fiore, Chiara Feruglio, Valentina D'Odorico, Nahum Arav, Tiago Costa, Kastytis Zubovas, George Becker, Sarah E. I. Bosman, Guido Cupani, Rebecca Davies, Anna-Christina Eilers, Emanuele Paolo Farina, Andrea Ferrara, Massimo Gaspari, Chiara Mazzucchelli, Masafusa Onoue, Enrico Piconcelli, Maria-Vittoria Zanchettin, Yongda Zhu

Abstract: Luminous quasars are powerful targets to investigate the role of feedback from supermassive black-holes (BHs) in regulating the growth phases of BHs themselves and of their host galaxies, up to the highest redshifts. Here we investigate the cosmic evolution of the occurrence and kinematics of BH-driven outflows, as traced by broad absorption line (BAL) features, due to the C IV ionic transition. W… ▽ More Luminous quasars are powerful targets to investigate the role of feedback from supermassive black-holes (BHs) in regulating the growth phases of BHs themselves and of their host galaxies, up to the highest redshifts. Here we investigate the cosmic evolution of the occurrence and kinematics of BH-driven outflows, as traced by broad absorption line (BAL) features, due to the C IV ionic transition. We exploit a sample of 1935 quasars quasars at $z=2.1-6.6$ with bolometric luminosity log($L_{\rm bol}/$erg s$^{-1})\gtrsim46.5$, drawn from the Sloan Digital Sky Survey and from the X-shooter legacy survey of Quasars at Reionisation (XQR-30). We consider rest-frame optical bright quasars to minimise observational biases due to quasar selection criteria. We apply a homogeneous BAL identification analysis, based on employing composite template spectra to estimate the quasar intrinsic emission. We find a BAL quasar fraction close to 20\% at $z\sim2-4$, while it increases to almost 50\% at $z\sim6$. The velocity and width of the BAL features also increase at $z\gtrsim4.5$. We exclude that the redshift evolution of the BAL properties is due to differences in terms of quasar luminosity and accretion rate. These results suggest significant BH feedback occurring in the 1 Gyr old Universe, likely affecting the growth of BHs and, possibly, of their host galaxies, as supported by models of early BH and galaxy evolution. △ Less

Submitted 18 April, 2023; v1 submitted 23 January, 2023; originally announced January 2023.

Comments: Accepted for publication in ApJ

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

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

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

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

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

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

Detection of stellar light from quasar host galaxies at redshifts above 6

Authors: Xuheng Ding, Masafusa Onoue, John D. Silverman, Yoshiki Matsuoka, Takuma Izumi, Michael A. Strauss, Knud Jahnke, Camryn L. Phillips, Junyao Li, Marta Volonteri, Zoltan Haiman, Irham Taufik Andika, Kentaro Aoki, Shunsuke Baba, Rebekka Bieri, Sarah E. I. Bosman, Connor Bottrell, Anna-Christina Eilers, Seiji Fujimoto, Melanie Habouzit, Masatoshi Imanishi, Kohei Inayoshi, Kazushi Iwasawa, Nobunari Kashikawa, Toshihiro Kawaguchi , et al. (19 additional authors not shown)

Abstract: The detection of starlight from the host galaxies of quasars during the reionization epoch ($z>6$) has been elusive, even with deep HST observations. The current highest redshift quasar host detected, at $z=4.5$, required the magnifying effect of a foreground lensing galaxy. Low-luminosity quasars from the Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP) mitigate the challenge of detecting the… ▽ More The detection of starlight from the host galaxies of quasars during the reionization epoch ($z>6$) has been elusive, even with deep HST observations. The current highest redshift quasar host detected, at $z=4.5$, required the magnifying effect of a foreground lensing galaxy. Low-luminosity quasars from the Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP) mitigate the challenge of detecting their underlying, previously-undetected host galaxies. Here we report rest-frame optical images and spectroscopy of two HSC-SSP quasars at $z>6$ with JWST. Using NIRCam imaging at 3.6$μ$m and 1.5$μ$m and subtracting the light from the unresolved quasars, we find that the host galaxies are massive (stellar masses of $13\times$ and $3.4\times$ $10^{10}$ M$_{\odot}$, respectively), compact, and disk-like. NIRSpec medium-resolution spectroscopy shows stellar absorption lines in the more massive quasar, confirming the detection of the host. Velocity-broadened gas in the vicinity of these quasars enables measurements of their black hole masses ($1.4\times 10^9$ and $2.0\times$ $10^{8}$ M$_{\odot}$, respectively). Their location in the black hole mass - stellar mass plane is consistent with the distribution at low redshift, suggesting that the relation between black holes and their host galaxies was already in place less than a billion years after the Big Bang. △ Less

Submitted 23 June, 2023; v1 submitted 25 November, 2022; originally announced November 2022.

Comments: Matched to the published Nature version of the article. 27 pages, 4 main figures, 1 table, 6 supplementary figures, 2 supplementary table

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

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

Widespread, strong outflows in XQR-30 quasars at the Reionisation epoch

Authors: M. Bischetti, C. Feruglio, V. D'Odorico, N. Arav, E. Bañados, G. Becker, S. E. I. Bosman, S. Carniani, S. Cristiani, G. Cupani, R. Davies, A. C. Eilers, E. P. Farina, A. Ferrara, R. Maiolino, C. Mazzucchelli, A. Mesinger, R. Meyer, M. Onoue, E. Piconcelli, E. Ryan-Weber, J-T. Schindler, F. Wang, J. Yang, Y. Zhu , et al. (1 additional authors not shown)

Abstract: Luminous quasars powered by accretion onto billion solar mass black holes already exist at the epoch of Reionisation, when the Universe was 0.5-1 Gyr old. These objects likely reside in over-dense regions of the Universe, and will grow to form today's giant galaxies. How their huge black holes formed in such short times is debated, particularly as they lie above the local black hole mass-galaxy dy… ▽ More Luminous quasars powered by accretion onto billion solar mass black holes already exist at the epoch of Reionisation, when the Universe was 0.5-1 Gyr old. These objects likely reside in over-dense regions of the Universe, and will grow to form today's giant galaxies. How their huge black holes formed in such short times is debated, particularly as they lie above the local black hole mass-galaxy dynamical mass correlation, thus following the black hole-dominance growth path. It is unknown what slowed down the black hole growth, leading towards the symbiotic growth observed in the local Universe, and when this process started, although black hole feedback is a likely driver. This deadlock is due to the lack of large, homogeneous samples of high-redshift quasars with high-quality, broad-band spectroscopic information. Here we report results from a VLT/X-shooter survey of 30 quasars at redshift 5.8$\le$z$\le$6.6 (XQR-30). About 50% of their spectra reveal broad blue-shifted absorption line (BAL) throughs, tracing powerful ionised winds. The BAL fraction in z$\gtrsim$6 quasars is 2-3 times higher than in quasars at z~2-4.5. XQR-30 BAL quasars exhibit extreme outflow velocities, up to 17% of the light speed, rarely observed at lower redshift. These outflows inject large amounts of energy into the galaxy interstellar medium, which can contrast nuclear gas accretion, slowing down the black-hole growth. The star-formation rate in high-z quasar hosts is generally $>$100 M$_\odot$/yr, so these galaxies are growing at a fast rate. The BAL phase may then mark the beginning of significant feedback, acting first on black hole growth and possibly later on galaxy growth. The red optical colors of BAL quasars at z$\gtrsim$6 indeed suggest that these systems are dusty and may be caught during an initial quenching phase of obscured accretion. △ Less

Submitted 17 May, 2022; v1 submitted 29 April, 2022; originally announced May 2022.

Comments: This preprint has not undergone peer review or any post-submission improvements or corrections. The Version of Record of this article is published in Nature, and is available online at https://doi.org/10.1038/s41586-022-04608-1

Journal ref: Nature, 605, 244-247 (2022)

Chemical Abundance of z~6 Quasar Broad-Line Regions in the XQR-30 Sample

Authors: Samuel Lai, Fuyan Bian, Christopher A. Onken, Christian Wolf, Chiara Mazzucchelli, Eduardo Banados, Manuela Bischetti, Sarah E. I. Bosman, George Becker, Guido Cupani, Valentina D'Odorico, Anna-Christina Eilers, Xiaohui Fan, Emanuele Paolo Farina, Masafusa Onoue, Jan-Torge Schindler, Fabian Walter, Feige Wang, Jinyi Yang, Yongda Zhu

Abstract: The elemental abundances in the broad-line regions of high-redshift quasars trace the chemical evolution in the nuclear regions of massive galaxies in the early universe. In this work, we study metallicity-sensitive broad emission-line flux ratios in rest-frame UV spectra of 25 high-redshift (5.8 < z < 7.5) quasars observed with the VLT/X-shooter and Gemini/GNIRS instruments, ranging over… ▽ More The elemental abundances in the broad-line regions of high-redshift quasars trace the chemical evolution in the nuclear regions of massive galaxies in the early universe. In this work, we study metallicity-sensitive broad emission-line flux ratios in rest-frame UV spectra of 25 high-redshift (5.8 < z < 7.5) quasars observed with the VLT/X-shooter and Gemini/GNIRS instruments, ranging over $\log(M_{\rm{BH}}/M_{\odot})= 8.4-9.8$ in black hole mass and $\log(L_{\rm{bol}}/\rm{erg\, s}^{-1})= 46.7-47.7$ in bolometric luminosity. We fit individual spectra and composites generated by binning across quasar properties: bolometric luminosity, black hole mass, and blueshift of the \civ\, line, finding no redshift evolution in the emission-line ratios by comparing our high-redshift quasars to lower-redshift (2.0 < z < 5.0) results presented in the literature. Using Cloudy-based locally optimally-emitting cloud photoionisation model relations between metallicity and emission-line flux ratios, we find the observable properties of the broad emission lines to be consistent with emission from gas clouds with metallicity that are at least 2-4 times solar. Our high-redshift measurements also confirm that the blueshift of the CIV emission line is correlated with its equivalent width, which influences line ratios normalised against CIV. When accounting for the CIV blueshift, we find that the rest-frame UV emission-line flux ratios do not correlate appreciably with the black hole mass or bolometric luminosity. △ Less

Submitted 7 April, 2022; originally announced April 2022.

Comments: MNRAS accepted, 24 pages, 17 figures, 4 tables

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

Constraining IGM enrichment and metallicity with the C IV forest correlation function

Authors: Suk Sien Tie, Joseph F. Hennawi, Koki Kakiichi, Sarah E. I. Bosman

Abstract: The production and distribution of metals in the diffuse intergalactic medium (IGM) have implications for galaxy formation models and the baryon (re)cycling process. The relative abundance of metals in high versus low-ionization states has also been argued to be sensitive to the Universe's reionization history. However, measurements of the background metallicity of the IGM at z~4 are sparse and in… ▽ More The production and distribution of metals in the diffuse intergalactic medium (IGM) have implications for galaxy formation models and the baryon (re)cycling process. The relative abundance of metals in high versus low-ionization states has also been argued to be sensitive to the Universe's reionization history. However, measurements of the background metallicity of the IGM at z~4 are sparse and in poor agreement with one another, and reduced sensitivity in the near-IR renders detecting individual metal absorbers nearly impossible. We present a new clustering-based technique that enables the detection of these weak IGM absorbers by statistically averaging over all spectral pixels, here applied to the C IV forest. We simulate the z=4.5 IGM with different models of inhomogeneous metal distributions and investigate its two-point correlation function (2PCF) using mock skewers of the C IV forest. The 2PCF demonstrates a clear peak at the doublet separation of the C IV line. The peak amplitude scales quadratically with metallicity, while enrichment morphology affects both the shape and amplitude of the 2PCF. The effect of enrichment topology can also be framed in terms of the metal mass- and volume-filling factors, and we show their trends as a function of the enrichment topology. For models consistent with the distribution of metals at z~3, we find that we can constrain [C/H] to within 0.2 dex, log$\,M_{\rm{min}}$ to within 0.4 dex, and $R$ to within 15%. We show that strong absorbers arising from the circumgalactic medium of galaxies can be easily identified and masked, allowing one to recover the underlying IGM signal. The auto-correlation of the metal-line forest presents a new and compelling avenue to simultaneously constrain IGM metallicity and enrichment topology with high precision at z>4, thereby pushing such measurements into the Epoch of Reionization. △ Less

Submitted 25 January, 2022; originally announced January 2022.

Comments: 19 pages, 17 figures, 1 appendix. Submitted to MNRAS

ALMA 200 pc imaging of a z~7 quasar reveals a compact, disk-like host galaxy

Authors: Fabian Walter, Marcel Neeleman, Roberto Decarli, Bram Venemans, Romain Meyer, Axel Weiss, Eduardo Banados, Sarah E. I. Bosman, Chris Carilli, Xiaohui Fan, Dominik Riechers, Hans-Walter Rix, Todd A. Thompson

Abstract: We present 0".035 resolution (~200 pc) imaging of the 158 um [CII] line and the underlying dust continuum of the z=6.9 quasar J234833.34-305410.0. The 18 h ALMA observations reveal extremely compact emission (diameter ~1 kpc) that is consistent with a simple, almost face-on, rotation-supported disk with a significant velocity dispersion of ~160 km/s. The gas mass in just the central 200 pc is ~4x1… ▽ More We present 0".035 resolution (~200 pc) imaging of the 158 um [CII] line and the underlying dust continuum of the z=6.9 quasar J234833.34-305410.0. The 18 h ALMA observations reveal extremely compact emission (diameter ~1 kpc) that is consistent with a simple, almost face-on, rotation-supported disk with a significant velocity dispersion of ~160 km/s. The gas mass in just the central 200 pc is ~4x10^9 M_sun, about a factor two higher than that of the central supermassive black hole. Consequently we do not resolve the black hole's sphere of influence, and find no kinematic signature of the central supermassive black hole. Kinematic modeling of the [CII] line shows that the dynamical mass at large radii is consistent with the gas mass, leaving little room for a significant mass contribution by stars and/or dark matter. The Toomre-Q parameter is less than unity throughout the disk, and thus is conducive to star formation, consistent with the high infrared luminosity of the system. The dust in the central region is optically thick, at a temperature >132 K. Using standard scaling relations of dust heating by star formation, this implies an unprecedented high star formation rate density of >10^4 M_sun / yr / kpc^2. Such a high number can still be explained with the Eddington limit for star formation under certain assumptions, but could also imply that the central supermassive black hole contributes to the heating of the dust in the central 110 pc. △ Less

Submitted 17 January, 2022; originally announced January 2022.

Comments: Accepted for publication in The Astrophysical Journal

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

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

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

Submitted 26 October, 2021; originally announced October 2021.

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

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

Constraints on the mean free path of ionising photons at $z\sim6$ using limits on individual free paths

Authors: Sarah E. I. Bosman

Abstract: The recent measurement of an ionising mean free path $λ_{\text{mfp}}<1$ pMpc at $z=6$ challenges our understanding of the small-scale structure of the intergalactic medium (IGM) at the end of reionisation. We introduce a new method to constrain \mfp at $z=6$ by using lower limits on the individual free paths of ionisation around quasars. Lyman-limit absorbers with a density sufficient to halt ioni… ▽ More The recent measurement of an ionising mean free path $λ_{\text{mfp}}<1$ pMpc at $z=6$ challenges our understanding of the small-scale structure of the intergalactic medium (IGM) at the end of reionisation. We introduce a new method to constrain \mfp at $z=6$ by using lower limits on the individual free paths of ionisation around quasars. Lyman-limit absorbers with a density sufficient to halt ionising photons produce strong absorption in the 6 lowest-energy Lyman transitions, in the absence of which a robust lower limit can be placed on the individual free path. Applying this method to a set of $26$ quasars at $5.5<z<6.5$, we find that $80\%$ of bright quasars ($M_{1450}<-26.5$) require individual free paths larger than $2$ pMpc. We model the relation between opacity $κ$ and photo-ionisation rate $Γ$ via the parameter $ξ$ such that $κ\proptoΓ^{-ξ}$, and pose joint limits on \mfp and $ξ$. For the nominal value of $ξ=2/3$, we constrain $λ_{\text{mfp}}>0.31 \ (0.18)$ pMpc at $2σ\ (3σ)$: a much tighter lower bound than obtained through traditional stacking methods. Our constraints get significantly stronger for lower values of $ξ$. New constraints on \mfp and $ξ$ are crucial to our understanding of the reionisation-era IGM. △ Less

Submitted 31 August, 2021; v1 submitted 27 August, 2021; originally announced August 2021.

Comments: Submitted to ApJ; 15 pages, 5 figures

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

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

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

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

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

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

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

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

Submitted 21 May, 2021; originally announced May 2021.

Comments: 8 pages, 5 figures, submitted to ApJL

Reionization and galaxy inference from the high-redshift Lyα forest

Authors: Yuxiang Qin, Andrei Mesinger, Sarah E. I. Bosman, Matteo Viel

Abstract: The transmission of Lyman-α (Lyα) in the spectra of distant quasars depends on the density, temperature, and ionization state of the intergalactic medium (IGM). Therefore, high-redshift (z > 5) Lyα forests could be invaluable in studying the late stages of the epoch of reionization (EoR), as well as properties of the sources that drive it. Indeed, high-quality quasar spectra have now firmly establ… ▽ More The transmission of Lyman-α (Lyα) in the spectra of distant quasars depends on the density, temperature, and ionization state of the intergalactic medium (IGM). Therefore, high-redshift (z > 5) Lyα forests could be invaluable in studying the late stages of the epoch of reionization (EoR), as well as properties of the sources that drive it. Indeed, high-quality quasar spectra have now firmly established the existence of large-scale opacity fluctuations at z > 5, whose physical origins are still debated. Here we introduce a Bayesian framework capable of constraining the EoR and galaxy properties by forward-modelling the high-z Lyα forest. Using priors from galaxy and CMB observations, we demonstrate that the final overlap stages of the EoR (when >95% of the volume was ionized) should occur at z < 5.6, in order to reproduce the large-scale opacity fluctuations seen in forest spectra. However, it is the combination of patchy reionization and the inhomogeneous UV background that produces the longest Gunn-Peterson troughs. Lyα forest observations tighten existing constraints on the characteristic ionizing escape fraction of galaxies, with the combined observations suggesting f_{\rm esc} \approx 7^4_3%, and disfavoring a strong evolution with the galaxy's halo (or stellar) mass. △ Less

Submitted 22 January, 2021; originally announced January 2021.

Comments: 19 pages, 11 figures, submitted to MNRAS, comments are welcome

A comparison of quasar emission reconstruction techniques for $z\geq5.0$ Lyman-$α$ and Lyman-$β$ transmission

Authors: S. E. I. Bosman, D. Ďurovčíková, F. B. Davies, A. C. Eilers

Abstract: Reconstruction techniques for intrinsic quasar continua are crucial for the precision study of Lyman-$α$ (Ly-$α$) and Lyman-$β$ (Ly-$β$) transmission at $z>5.0$, where the $λ<1215 A$ emission of quasars is nearly completely absorbed. While the number and quality of spectroscopic observations has become theoretically sufficient to quantify Ly-$α$ transmission at $5.0<z<6.0$ to better than $1\%$, th… ▽ More Reconstruction techniques for intrinsic quasar continua are crucial for the precision study of Lyman-$α$ (Ly-$α$) and Lyman-$β$ (Ly-$β$) transmission at $z>5.0$, where the $λ<1215 A$ emission of quasars is nearly completely absorbed. While the number and quality of spectroscopic observations has become theoretically sufficient to quantify Ly-$α$ transmission at $5.0<z<6.0$ to better than $1\%$, the biases and uncertainties arising from predicting the unabsorbed continuum are not known to the same level. In this paper, we systematically evaluate eight reconstruction techniques on a unified testing sample of $2.7<z<3.5$ quasars drawn from eBOSS. The methods include power-law extrapolation, stacking of neighbours, and six variants of Principal Component Analysis (PCA) using direct projection, fitting of components, or neural networks to perform weight mapping. We find that power-law reconstructions and the PCA with fewest components and smallest training sample display the largest biases in the Ly-$α$ forest ($-9.58\%/+8.22\%$ respectively). Power-law extrapolations have larger scatters than previously assumed of $+13.1\%/-13.2\%$ over Ly-$α$ and $+19.9\%/-20.1\%$ over Ly-$β$. We present two new PCAs which achieve the best current accuracies of $9\%$ for Ly-$α$ and $17\%$ for Ly-$β$. We apply the eight techniques after accounting for wavelength-dependent biases and scatter to a sample $19$ quasars at $z>5.7$ with IR X-Shooter spectroscopy, obtaining well-characterised measurements for the mean flux transmission at $4.7<z<6.3$. Our results demonstrate the importance of testing and, when relevant, training, continuum reconstruction techniques in a systematic way. △ Less

Submitted 17 March, 2021; v1 submitted 18 June, 2020; originally announced June 2020.

Comments: Accepted for publication in MNRAS; new Appendix added, conclusions unchanged. 13 pages plus Appendix

Lyman-alpha transmission properties of the intergalactic medium in the CoDaII simulation

Authors: Max Gronke, Pierre Ocvirk, Charlotte Mason, Jorryt Matthee, Sarah E. I. Bosman, Jenny G. Sorce, Joseph Lewis, Kyungjin Ahn, Dominique Aubert, Taha Dawoodbhoy, Ilian T. Iliev, Paul R. Shapiro, Gustavo Yepes

Abstract: The decline in abundance of Lyman-$α$ (Ly$α$) emitting galaxies at $z \gtrsim 6$ is a powerful and commonly used probe to constrain the progress of cosmic reionization. We use the CoDaII simulation, which is a radiation hydrodynamic simulation featuring a box of $\sim 94$ comoving Mpc side length, to compute the Ly$α$ transmission properties of the intergalactic medium (IGM) at $z\sim 5.8$ to $7$.… ▽ More The decline in abundance of Lyman-$α$ (Ly$α$) emitting galaxies at $z \gtrsim 6$ is a powerful and commonly used probe to constrain the progress of cosmic reionization. We use the CoDaII simulation, which is a radiation hydrodynamic simulation featuring a box of $\sim 94$ comoving Mpc side length, to compute the Ly$α$ transmission properties of the intergalactic medium (IGM) at $z\sim 5.8$ to $7$. Our results mainly confirm previous studies, i.e., we find a declining Ly$α$ transmission with redshift and a large sightline-to-sightline variation. However, motivated by the recent discovery of blue Ly$α$ peaks at high redshift, we also analyze the IGM transmission on the blue side, which shows a rapid decline at $z\gtrsim 6$ of the blue transmission. This low transmission can be attributed not only to the presence of neutral regions but also to the residual neutral hydrogen within ionized regions, for which a density even as low as $n_{\rm HI}\sim 10^{-9}\,\mathrm{cm}^{-3}$ (sometimes combined with kinematic effects) leads to a significantly reduced visibility. Still, we find that $\sim 1\%$ of sightlines towards $M_{\mathrm{1600AB}}\sim -21$ galaxies at $z\sim 7$ are transparent enough to allow a transmission of a blue Ly$α$ peak. We discuss our results in the context of the interpretation of observations. △ Less

Submitted 21 October, 2021; v1 submitted 29 April, 2020; originally announced April 2020.

Comments: 13 pages, 12 figures. Published version

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 508, Issue 3, December 2021, Pages 3697-3709

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

Authors: T. Roy Choudhury, Aseem Paranjape, Sarah E. I. Bosman

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

Submitted 17 December, 2020; v1 submitted 19 March, 2020; originally announced March 2020.

Comments: Accepted in MNRAS. Matches the accepted version

Four direct measurements of the fine-structure constant 13 billion years ago

Authors: Michael R. Wilczynska, John K. Webb, Matthew Bainbridge, Sarah E. I. Bosman, John D. Barrow, Robert F. Carswell, Mariusz P. Dabrowski, Vincent Dumont, Ana Catarina Leite, Chung-Chi Lee, Katarzyna Leszczynska, Jochen Liske, Konrad Marosek, Carlos J. A. P. Martins, Dinko Milakovic, Paolo Molaro, Luca Pasquini

Abstract: Observations of the redshift z=7.085 quasar J1120+0641 have been used to search for variations of the fine structure constant, alpha, over the redshift range 5.5 to 7.1. Observations at z=7.1 probe the physics of the universe when it was only 0.8 billion years old. These are the most distant direct measurements of alpha to date and the first measurements made with a near-IR spectrograph. A new AI… ▽ More Observations of the redshift z=7.085 quasar J1120+0641 have been used to search for variations of the fine structure constant, alpha, over the redshift range 5.5 to 7.1. Observations at z=7.1 probe the physics of the universe when it was only 0.8 billion years old. These are the most distant direct measurements of alpha to date and the first measurements made with a near-IR spectrograph. A new AI analysis method has been employed. Four measurements from the X-SHOOTER spectrograph on the European Southern Observatory's Very Large Telescope (VLT) directly constrain any changes in alpha relative to the value measured on Earth (alpha_0). The weighted mean strength of the electromagnetic force over this redshift range in this location in the universe is da/a = (alpha_z - alpha_0)/alpha_0 = (-2.18 +/- 7.27) X 10^{-5}, i.e. we find no evidence for a temporal change from the 4 new very high redshift measurements. When the 4 new measurements are combined with a large existing sample of lower redshift measurements, a new limit on possible spatial variation of da/a is marginally preferred over a no-variation model at the 3.7 sigma level. △ Less

Submitted 17 March, 2020; originally announced March 2020.

Comments: Accepted for publication in Science Advances

Three Lyman-alpha emitting galaxies within a quasar proximity zone at z=5.8

Authors: Sarah E. I. Bosman, Koki Kakiichi, Romain A. Meyer, Max Gronke, Nicolas Laporte, Richard S. Ellis

Abstract: Quasar proximity zones at $z>5.5$ correspond to over-dense and over-ionized environments. Galaxies found inside proximity zones can therefore display features which would otherwise be masked by absorption in the IGM. We demonstrate the utility of this quasar-galaxy synergy by reporting the discovery of the first three `proximate Lyman-$α$ emitters' (LAEs) within the proximity zone of quasar J0836… ▽ More Quasar proximity zones at $z>5.5$ correspond to over-dense and over-ionized environments. Galaxies found inside proximity zones can therefore display features which would otherwise be masked by absorption in the IGM. We demonstrate the utility of this quasar-galaxy synergy by reporting the discovery of the first three `proximate Lyman-$α$ emitters' (LAEs) within the proximity zone of quasar J0836 at $z=5.802$ (\textit{Aerith A, B} and \textit{C}). \textit{Aerith A}, located behind the quasar with an impact parameter $D_\perp = 278$ pkpc, provides the first detection of a Lyman-$α$ transverse proximity effect. We model the transmission and show it constrains the onset of J0836's quasar phase to $0.2 \text{Myr}<t<20\text{Myr}$ in the past. The second object, \textit{Aerith B} at a distance $D=750$ pkpc from the quasar, displays a bright, broad double-peaked \lal emission line. Based on relations calibrated at $z\leq3$, the peak separation implies a low ionizing $f_{\text{esc}} \lesssim 1\%$, the most direct such constraint on a reionization-era galaxy. We fit the Ly-$α$ line with an outflowing shell model, finding a completely typical central density $\text{log N}_{\text{HI}}/\text{cm}^{-2} = 19.3_{-0.2}^{+0.8}$, outflow velocity $v=16_{-11}^{+4}$ km s$^{-1}$, and gas temperature $\text{log} T/\text{K} = 3.8_{-0.7}^{+0.8}$ compared to $2<z<3$ analogue LAEs. Finally, we detect an emission line at $λ=8177$ Å in object \textit{Aerith C} which, if it is \lal at $z=5.726$, would correspond closely with the end of the quasar's proximity zone ($Δz<0.02$ from the boundary) and suggests the quasar influences the IGM up to $\sim85$ cMpc away, making it the largest quasar proximity zone. Via the analyses conducted here, we illustrate how proximate LAEs offer unique insight into the ionizing properties of both quasars and galaxies during the epoch of reionization. △ Less

Submitted 24 December, 2019; originally announced December 2019.

Comments: Submitted to ApJ. Comments welcome