Showing 1–50 of 255 results for author: Hennawi, J F

Quasar Impostors: Two Extremely UV-Bright ($M_{\rm UV}\approx-23.5$) Reionisation-Epoch Galaxies Powered by Very Massive Stars

Authors: Daming Yang, Joseph F. Hennawi, Sarah E. I. Bosman, Frederick B. Davies, Rychard Bouwens, Timo Kist, Eduardo Bañados, Jiamu Huang, Alice E. Shapley, Marianne Vestergaard, Hiddo S. B. Algera, Silvia Belladitta, Anna-Christina Eilers, Xiaohui Fan, Francesco Guarneri, Xiangyu Jin, Romain A. Meyer, Elia Pizzati, Huub Röttgering, Jan-Torge Schindler, Mauro Stefanon, Feige Wang, Jinyi Yang

Abstract: The extreme bright end of the galaxy UV luminosity function during reionisation remains poorly constrained, particularly where the galaxy and quasar luminosity functions overlap and source classification becomes ambiguous. We present JWST/NIRSpec and ALMA Band-6 observations of J1450-0144 ($z=6.627$) and J1429-0104 ($z=6.796$), two $M_{\rm UV}\simeq-23.5$ sources originally classified as faint qua… ▽ More The extreme bright end of the galaxy UV luminosity function during reionisation remains poorly constrained, particularly where the galaxy and quasar luminosity functions overlap and source classification becomes ambiguous. We present JWST/NIRSpec and ALMA Band-6 observations of J1450-0144 ($z=6.627$) and J1429-0104 ($z=6.796$), two $M_{\rm UV}\simeq-23.5$ sources originally classified as faint quasars by SHELLQs. NIRSpec reveals blue UV continua, strong P Cygni profiles in N V, Si IV, and C IV, broad He II $\lambda1640$ emission with rest-frame equivalent widths of $8.8\pm1.2$ and $3.7\pm1.1$ Å, respectively, and narrow nebular lines, reclassifying both as extremely UV-luminous galaxies. Standard population-synthesis models cannot simultaneously reproduce the strong He II and wind features, whereas models incorporating very massive stars (VMS; $M\gtrsim100\,M_\odot$) with dedicated wind prescriptions can. These models favor a star-formation duration of 2-4 Myr for J1450-0144, with a broader allowed range for J1429-0104, stellar masses of $\log(M_\star/M_\odot)\approx9.2$-$9.9$, and star-formation rates of $\simeq300$-$540\,M_\odot\,{\rm yr}^{-1}$. Under the same VMS wind models, equivalent-width diagnostics imply $M_{\rm up}\gtrsim225\,M_\odot$ for J1429-0104, while J1450-0144 lies beyond even the $M_{\rm up}=475\,M_\odot$ grid. ALMA detects luminous [C II] 158 $μ$m emission in both systems, with $L_{\rm [CII]}\approx0.8$ and $4.1\times10^{9}\,L_\odot$, respectively. J1429-0104 additionally shows bright dust continuum, with both [C II] and dust offset by $\sim5.4$ kpc from its UV emission. These sources demonstrate that VMS can power some of the most UV-luminous galaxies at cosmic dawn and show that source classifications, and hence the inferred demographics of both galaxies and quasars in the crossover regime, need revisiting. △ Less

Submitted 18 August, 2026; originally announced August 2026.

Comments: 18 figures, 4 tables. Submitted

Euclid: Discovery of 31 new quasars at $6.6 < z < 7.8$

Authors: D. Yang, J. F. Hennawi, F. Guarneri, J. Wolf, S. Belladitta, J. -T. Schindler, A. C. N. Hughes, E. Bañados, D. J. Mortlock, J. Yang, F. Wang, X. Fan, K. Jahnke, D. Stern, C. J. Willott, A. J. Barth, H. J. A. Rottgering, R. G. Varadaraj, R. Decarli, A. -C. Eilers, M. Ezziati, Y. Fu, J. Huang, X. Jin, Y. Kang , et al. (189 additional authors not shown)

Abstract: We report the discovery of 31 new high-$z$ quasars in the redshift range $6.6 < z < 7.8$. These quasars were selected from approximately 3000 deg$^2$ of sky covered during the first 1.5 years of the Euclid Wide Survey, representing the initial results of the Euclid high-$z$ quasar search. Our candidate selection employed multiple machine-learning and probabilistic techniques applied to the Euclid… ▽ More We report the discovery of 31 new high-$z$ quasars in the redshift range $6.6 < z < 7.8$. These quasars were selected from approximately 3000 deg$^2$ of sky covered during the first 1.5 years of the Euclid Wide Survey, representing the initial results of the Euclid high-$z$ quasar search. Our candidate selection employed multiple machine-learning and probabilistic techniques applied to the Euclid $I_E$, $Y_E$, $J_E$, and $H_E$ images, supplemented by ancillary $z$-band data when available. Spectroscopic follow-up observations were carried out with Keck, Magellan, and the Large Binocular Telescope (LBT). Among the new discoveries, there are 12 quasars at $z \geq 7$, more than doubling the number of previously known quasars at $z \geq 7$. The newly discovered quasars exhibit $21.2 < J_E < 23.2$ ($-25.5 < M_{1450} < -23.6$), extending quasar studies to the faint end of the quasar luminosity function (QLF) at $z \gtrsim 7$. The quasar with the highest-$z$, EUCL J172902.75+641018.1 at $z \approx 7.77$, sets the new redshift record for the most distant quasar ever reported. These discoveries demonstrate Euclid's transformative role in high-$z$ quasar discovery and set the stage for future follow-up studies of the early galaxies hosting quasars, supermassive black hole growth, and the intergalactic medium in the epoch of reionisation. △ Less

Submitted 3 July, 2026; originally announced July 2026.

Comments: Accepted for publication in A&A

Journal ref: A&A, 711, A104 (2026)

Euclid: A UV-faint quasar in a highly luminous star-forming host galaxy at $z \approx 7.7$

Authors: S. Belladitta, R. Decarli, E. Bañados, F. Walter, D. Yang, F. Guarneri, K. Jahnke, S. Bisogni, S. E. I. Bosman, X. Fan, Y. Fu, J. F. Hennawi, Y. Matsuoka, D. J. Mortlock, M. Onoue, J. -T. Schindler, L. Spinoglio, D. Stern, F. Wang, G. Vietri, C. J. Willott, J. Wolf, J. Yang, R. A. A. Bowler, K. I. Caputi , et al. (158 additional authors not shown)

Abstract: Constraining the co-evolution of supermassive black holes and their host galaxies in the first billion years after the Big Bang is essential for understanding the formation of the earliest cosmic structures. Here, we present IRAM/NOrthern Extended Millimeter Array (NOEMA) observations of the $z \approx 7.7$ quasar EUCL\,J125308.55+705432.3, recently discovered in the first data release of the Eucl… ▽ More Constraining the co-evolution of supermassive black holes and their host galaxies in the first billion years after the Big Bang is essential for understanding the formation of the earliest cosmic structures. Here, we present IRAM/NOrthern Extended Millimeter Array (NOEMA) observations of the $z \approx 7.7$ quasar EUCL\,J125308.55+705432.3, recently discovered in the first data release of the Euclid Wide Survey. We report the most distant detections of [CII] 158$μ\mathrm{m}$ and cold dust emission in a quasar host to date. The [CII] emission line sets the systemic redshift at $z=7.6980\pm0.0004$. The source exhibits luminosities of $L_{\rm FIR}=3.6\times10^{12}\,L_{\odot}$ and $L_{[CII]}=2\times10^9\,L_{\odot}$, respectively, a dust mass of 1.4$\times 10^{8}\,M_{\odot}$, and a dynamical mass in the range 0.33-1.3$\times10^{10}\,M_{\odot}$. Remarkably, despite being nearly two magnitudes fainter in the rest-frame UV ($M_{1450}=-24.06$) than previously known $z\approx7.5$ quasars (<M$_{1450}$> $\sim-$26.5), EUCL\,J125308.55+705432.3 exhibits the brightest [CII] emission among them. This indicates that the host galaxy is actively star-forming, with a star-formation rate $>250\,M_{\odot}\,\mathrm{yr}^{-1}$, consistent with recent findings that UV-faint quasars at $z>6$ preferentially reside in [CII]-luminous galaxies. The UV-faintness likely reflects dust obscuration or sub-Eddington accretion, rather than lower host mass, suggesting these systems are at a different stage in their evolution compared to UV-bright quasars. These IRAM/NOEMA observations highlight the power of combining Euclid's wide-area quasar discovery potential with submillimetre follow-up observations to characterise the host galaxies of early supermassive black holes across a broader redshift and luminosity range than previously accessible. △ Less

Submitted 3 July, 2026; originally announced July 2026.

Comments: 12 pages, 6 figures, 1 table, Accepted for Publication in A&A

A Measurement of the Thermal and Ionization State of the IGM at $z < 0.5$

Authors: Teng Hu, Vikram Khaire, Joseph F. Hennawi, Todd M. Tripp, Jose Onorbe, Michael Walther, Zarija Lukic

Abstract: We apply a machine-learning-based inference method that exploits the joint Doppler parameter-column density (b-NHI) distribution from Lya forest decomposition to measure the thermal and ionization state of the intergalactic medium (IGM) in four redshift bins spanning z = 0.06 to 0.48, using 82 archival quasar spectra from the Cosmic Origin Spectrograph (COS) on board Hubble Space Telescope (HST).… ▽ More We apply a machine-learning-based inference method that exploits the joint Doppler parameter-column density (b-NHI) distribution from Lya forest decomposition to measure the thermal and ionization state of the intergalactic medium (IGM) in four redshift bins spanning z = 0.06 to 0.48, using 82 archival quasar spectra from the Cosmic Origin Spectrograph (COS) on board Hubble Space Telescope (HST). Our results show that the low-z IGM (z < 0.5) is extremely hot and nearly isothermal, with log(T0/K) = 4.45 (+0.08 / -0.12) [T0 = 28183 (+5700 / -6804) K] and gamma = 1.06 (+0.13 / -0.09) at z = 0.1. This temperature lies approx 7sigma (and 7 times) above the canonical prediction (log T0 approx 3.60, i.e. T0 ~ 4000 K, with gamma ~ 1.6 at z = 0), where the IGM is expected to have cooled long after He II reionization. We also measure the hydrogen photoionization rate to be log (GammaHI/s^-1) = -13.70 (+0.10 / -0.08) at z = 0.1, which is about approx 4sigma below the range predicted by current UV-background synthesis models (approx -13.3). To investigate the discrepancy between these high temperatures and theoretical models, we assess the impact of small-scale turbulence. By exploring a parameter grid in turbulent velocity (vtur) and GammaHI, we find that a standard IGM thermal and ionization state combined with unresolved turbulence of vtur simeq 15 km s^-1 can successfully reproduce the observed line widths at z = 0.1. Comparisons with high-resolution Space Telescope Imaging Spectrograph (STIS) expanded data indicate that the observed line widths are unlikely to be caused by instrumental resolution effects. Our findings suggest that either new heating mechanisms or unresolved turbulence are required to explain the unexpectedly broad Lya lines observed in the low-z IGM. △ Less

Submitted 3 June, 2026; originally announced June 2026.

Comments: 20 pages, 16 figures

The Impact of Cosmic Variance and Satellites on JWST Clustering Measurements at Redshift around 6

Authors: Jiamu Huang, Elia Pizzati, Joseph F. Hennawi, Joop Schaye, Matthieu Schaller, Benjamin Snyder, Yi Kang

Abstract: We present a framework for inferring the dark matter halo masses of quasars and [O III]-emitting galaxies from JWST/NIRCam Wide Field Slitless Spectroscopy (WFSS) clustering measurements at z approximately 6. Using the FLAMINGO-10k N-body simulation, we construct mock realizations of quasar and galaxy catalogs that incorporate realistic selection functions, spatial coverage, and sensitivity limits… ▽ More We present a framework for inferring the dark matter halo masses of quasars and [O III]-emitting galaxies from JWST/NIRCam Wide Field Slitless Spectroscopy (WFSS) clustering measurements at z approximately 6. Using the FLAMINGO-10k N-body simulation, we construct mock realizations of quasar and galaxy catalogs that incorporate realistic selection functions, spatial coverage, and sensitivity limits matched to the ASPIRE survey. These mocks enable accurate measurements of the quasar-galaxy cross-correlation and galaxy auto-correlation functions, with covariance matrices derived from 1000 realizations that capture both cosmic variance and bin-to-bin correlations. We employ Bayesian inference to fit the correlation functions and infer the minimum halo masses for quasars and galaxies. Our results demonstrate that Poisson pair-count uncertainties, commonly adopted in high-redshift clustering studies, significantly underestimate the true measurement errors. The dominant missing component is cosmic variance: even the diagonal of the full covariance matrix exceeds the Poisson expectation, with off-diagonal bin-to-bin correlations contributing a smaller additional correction. In particular, 1) the commonly used Poisson error on the correlation functions underestimates the true uncertainty by a factor of approximately 3; 2) the uncertainties on the inferred minimum halo masses are underestimated by a factor of approximately 1.5-3 when adopting Poisson errors instead of the full covariance matrix; 3) the inferred QSO halo mass is robust to whether central and satellite [O III]-emitters share a common mass threshold. Our framework provides a more complete error budget for JWST/WFSS clustering analyses, enabling robust constraints on the host halo masses and duty cycles of high-redshift quasars and emission-line galaxies. △ Less

Submitted 11 May, 2026; originally announced May 2026.

Comments: 18 pages, 14 figures

Detection of an Extended Ly$α$ Halo around a $\textit{z}=6.64$ Broad Absorption Line Quasar with the Keck Cosmic Web Imager

Authors: Raymond P. Remigio, Aaron J. Barth, Feige Wang, Jinyi Yang, Joseph F. Hennawi, Ryan J. Cooke, Eduardo Banados, Xiaohui Fan, Emanuele Paolo Farina

Abstract: We present the first results from a program searching for extended Ly$α$ halos around high redshift ($ z \gtrsim 6.5$) quasars using the red channel of the Keck Cosmic Web Imager (KCWI). Our observations reveal a Ly$α$ halo extending to $\simeq11$ pkpc around the $z=6.64$ broad absorption line quasar J0910$-$0414. The Ly$α$ velocity field displays a rotation-like gradient, and the gas velocity dis… ▽ More We present the first results from a program searching for extended Ly$α$ halos around high redshift ($ z \gtrsim 6.5$) quasars using the red channel of the Keck Cosmic Web Imager (KCWI). Our observations reveal a Ly$α$ halo extending to $\simeq11$ pkpc around the $z=6.64$ broad absorption line quasar J0910$-$0414. The Ly$α$ velocity field displays a rotation-like gradient, and the gas velocity dispersion is consistent with gravitationally dominated motion ($σ_{\mathrm{Lyα}}<300$ km s$^{-1}$). Comparison with the $[\mathrm{C\;II}]$ kinematics of the host galaxy core from ALMA observations shows that the Ly$α$-emitting gas extends over a much larger region, shows distinct kinematics, and has a smaller velocity dispersion ($σ_{\mathrm{Lyα}} \simeq 0.6σ_{\mathrm{[C\;II]}}$). The Ly$α$ spectral region of the quasar is largely obscured by a deep $\mathrm{N\;V}$ absorption trough, and as a result, roughly $55\%$ of the total Ly$α$ flux is from the extended halo. These observations demonstrate the potential of KCWI for probing the cool gas reservoir that fuels the growth of quasars and their hosts in the epoch of reionization. △ Less

Submitted 7 May, 2026; originally announced May 2026.

Comments: 10 pages, 5 figures, 1 table. Accepted for publication into The Astrophysical Journal Letters

Three Hundred Quasars from the Couch: A first look at high-redshift quasar discovery with SPHEREx

Authors: Frederick B. Davies, Sarah E. I. Bosman, Arpita Ganguly, Eduardo Bañados, Silvia Belladitta, Daniel Stern, Javier A. Acevedo Barroso, Daming Yang, Joseph F. Hennawi, Feige Wang, Jinyi Yang, Xiaohui Fan

Abstract: Photometric selection of luminous high-redshift ($z\gtrsim4$) quasars is plagued by contamination from numerous low-mass Galactic stars, reddened lower-redshift quasars, as well as compact luminous red galaxies. Confirmation of these rare objects thus requires extensive spectroscopic campaigns on 4 and 8-meter-class telescopes with relatively low success rates. Here we demonstrate the utility of S… ▽ More Photometric selection of luminous high-redshift ($z\gtrsim4$) quasars is plagued by contamination from numerous low-mass Galactic stars, reddened lower-redshift quasars, as well as compact luminous red galaxies. Confirmation of these rare objects thus requires extensive spectroscopic campaigns on 4 and 8-meter-class telescopes with relatively low success rates. Here we demonstrate the utility of SPHEREx spectrophotometric survey data for quasar confirmation with no ground-based follow-up required, "from the couch," applied to candidates from a purposefully simplistic photometric and astrometric Gaia+WISE selection down to low Galactic latitudes ($|b|\geq8^\circ$). Primarily from the detection of their strong broad H$α$ emission lines, we discover 87 new luminous $4.0 < z < 5.7$ quasars with median $M_\text{1450} = -27.5$, including 19 quasars at $z>5$, and recover 219 previously published quasars at $z>4$. We validate our SPHEREx selection with a 100% confirmation rate in ground-based spectroscopic follow-up of 29 of our new $z>4$ quasars, including 11 unpublished archival spectra. We also discover 203 additional lower-redshift quasars at $0.3 < z < 4$, consisting primarily of relatively rare highly-reddened and strong broad-absorption-line objects that are likely missed by traditional quasar surveys. Finally, we show that the Ly$α$ absorption breaks and H$α$ lines of luminous quasars are already detectable at redshifts $5.7\lesssim z\lesssim6.5$ after the completion of only the first of four all-sky surveys to be performed by SPHEREx during its planned two-year mission. △ Less

Submitted 10 March, 2026; originally announced March 2026.

Comments: 12+4 pages, 11+2 figures, submitted to A&A

ASPIRE: The Environments and Dark Matter Halos of Luminous Quasars in the Epoch of Reionization

Authors: Feige Wang, Jaclyn B. Champagne, Jiamu Huang, Jinyi Yang, Joseph F. Hennawi, Xiaohui Fan, Haowen Zhang, Tiago Costa, Roberto Decarli, Melanie Habouzit, Fengwu Sun, Eduardo Banados, Xiangyu Jin, Koki Kakiichi, Romain A. Meyer, Yunjing Wu, Silvia Belladitta, Laura Blecha, Sarah E. I. Bosman, Zheng Cai, Thomas Connor, Frederick B. Davies, Anna-Christina Eilers, Zoltan Haiman, Hyunsung D. Jun , et al. (20 additional authors not shown)

Abstract: We present a systematic study of the environments of 25 luminous quasars at $z > 6.5$ from the ASPIRE program. Using JWST/NIRCam WFSS data, we identified 487 galaxies at $5.3 \lesssim z \lesssim 7.0$ exhibiting [OIII] emission. Among these, 122 [OIII] emitters lie within $|Δv_{\rm los}| < 1000~{\rm km~s^{-1}}$ of the quasars, corresponding to a $\sim9.4$-fold enhancement relative to the average ga… ▽ More We present a systematic study of the environments of 25 luminous quasars at $z > 6.5$ from the ASPIRE program. Using JWST/NIRCam WFSS data, we identified 487 galaxies at $5.3 \lesssim z \lesssim 7.0$ exhibiting [OIII] emission. Among these, 122 [OIII] emitters lie within $|Δv_{\rm los}| < 1000~{\rm km~s^{-1}}$ of the quasars, corresponding to a $\sim9.4$-fold enhancement relative to the average galaxy density at other redshifts. Furthermore, we identified 16 [CII]-emitting galaxies at the quasar redshifts from ALMA mosaic observations. A cross-correlation function (CCF) analysis between quasars and [OIII]+[CII] emitters yields a cross-correlation length of $r_0^{\rm QG} = 8.68^{+0.51}_{-0.55}~h^{-1}~\mathrm{cMpc}$ and a auto-correlation of $r_0^{\rm{QQ}}=15.76^{+2.48}_{-2.70}~h^{-1}~{\rm cMpc}$, indicating that $z \sim 7$ quasars reside in dark matter halos with $M_{\rm halo} = 10^{12.27^{+0.21}_{-0.26}}~M_\odot$. Notably, the number of [OIII]-emitting galaxies at quasar redshifts varies significantly from field to field, ranging from zero to twenty, highlighting a diverse quasar environment. Remarkably, seven quasars trace significant galaxy overdensities (i.e., protoclusters), with $δ_{\rm gal} > 5$ within a volume of $V \sim 500~{\rm cMpc^3}$. We also find that $|Δv_{\rm los}|$ increases rapidly toward smaller galaxy-quasar separations in protocluster fields, consistent with galaxy kinematics around extremely massive halos in cosmological simulations. By combining JWST and ALMA data, we reveal the complex and diverse environments of these early quasars, providing robust evidence that the earliest luminous quasars are effective tracers of galaxy overdensities, albeit with substantial field-to-field variation. △ Less

Submitted 30 June, 2026; v1 submitted 4 February, 2026; originally announced February 2026.

Comments: Accepted for publication in ApJ

A close look at the black hole masses and hot dusty toruses of the first quasars with MIRI-MRS

Authors: Sarah E. I. Bosman, Javier Álvarez-Márquez, Frederick B. Davies, Klaudia Protušová, Joseph F. Hennawi, Jinyi Yang, Benedetta Spina, Luis Colina, Xiaohui Fan, Göran Östlin, Fabian Walter, Feige Wang, Martin Ward, Almudena Alonso Herrero, Aaron J. Barth, Silvia Belladitta, Leindert Boogaard, Karina I. Caputi, Thomas Connor, Dominika Ďurovčíková, Anna-Christina Eilers, Alejandro Crespo Gómez, Jens Hjorth, Hyunsung D. Jun, Danial Langeroodi , et al. (7 additional authors not shown)

Abstract: The presence of supermassive black holes (SMBHs, $M_\text{BH}\sim10^9 M_\odot$) at $z>7$ remains a puzzle. While their existence appears to require exotic formation or growth processes, it is possible that BH mass estimates are incorrect due to differences from the low-$z$ quasars where BH mass scaling relations are calibrated. In this work, we employ JWST MIRI-MRS spectroscopy to measure the rest… ▽ More The presence of supermassive black holes (SMBHs, $M_\text{BH}\sim10^9 M_\odot$) at $z>7$ remains a puzzle. While their existence appears to require exotic formation or growth processes, it is possible that BH mass estimates are incorrect due to differences from the low-$z$ quasars where BH mass scaling relations are calibrated. In this work, we employ JWST MIRI-MRS spectroscopy to measure the rest-frame optical/IR properties of the four highest-redshift known luminous type-1 quasars at $7.08\leq z<7.64$. We use three new broad lines to measure updated BH masses, H$α$, Pa$α$ and Pa$β$, finding them to be in the range $(4-15)\cdot10^8 M_\odot$. Our black hole mass estimates from all tracers agree with each other and with previous, less accurate, ground-based measurements based on MgII. The flux ratios of the H lines deviate from expectations for case A and B recombination in the same way as in $z<3$ quasars, indicating similar physical conditions in the Broad Line Region. Rest-frame near-IR continuum emission from a hot dusty torus surrounding the accretion disc is unambiguously detected in all four objects. We model the emission with SKIRTOR and constrain the inclination (face-on) and the opening angle ($θ=40-60^\circ$) of the tori. These constraints are consistent for the four objects and with expectations from luminous quasars. We estimate a total dust mass $(1-4)\cdot10^6 M_\odot$ in the tori, corresponding to $(0.2-7)\%$ of the total dust in the quasar host galaxies. Given observed accretion rates, these SMBHs will deplete their tori in only $\sim5$ Myr. Overall, we confirm that $z>7$ SMBHs in quasars could not have grown from stellar-remnant BHs if the radiative efficiency of accretion is $10\%$. We also find no evidence that inferred BH masses and accretion processes in $z>7$ quasars differ significantly from their near-identical counterparts at $z<3$. △ Less

Submitted 4 November, 2025; originally announced November 2025.

Comments: Submitted to ApJ; 26 pages, 12 figures

Reliable Parameter Inference for the Epoch of Reionization using Balanced Neural Ratio Estimation

Authors: Diego González-Hernández, Molly Wolfson, Joseph F. Hennawi

Abstract: We present an application of the Balanced Neural Ratio Estimation (BNRE) algorithm to improve the statistical validity of parameter estimates used to characterize the Epoch of Reionization, where the common assumption of a multivariate Gaussian likelihood leads to overconfident and biased posterior distributions. Using a two-parameter model of the Ly$α$ forest autocorrelation function, we show tha… ▽ More We present an application of the Balanced Neural Ratio Estimation (BNRE) algorithm to improve the statistical validity of parameter estimates used to characterize the Epoch of Reionization, where the common assumption of a multivariate Gaussian likelihood leads to overconfident and biased posterior distributions. Using a two-parameter model of the Ly$α$ forest autocorrelation function, we show that BNRE yields posterior distributions that are significantly better calibrated than those obtained under the Gaussian likelihood assumption, as verified through the Test of Accuracy with Random Points (TARP) and Simulation-Based Calibration (SBC) diagnostics. These results demonstrate the potential of Simulation-Based Inference (SBI) methods, and in particular BNRE, to provide statistically robust parameter constraints within existing astrophysical modeling frameworks. △ Less

Submitted 11 November, 2025; v1 submitted 4 November, 2025; originally announced November 2025.

Comments: Machine Learning and the Physical Sciences Workshop, NeurIPS 2025

The luminosity function and clustering of bright quasars in the FLAMINGO cosmological simulations

Authors: Boyi Ding, Elia Pizzati, Joop Schaye, Joseph F. Hennawi, William McDonald, Matthieu Schaller

Abstract: Cosmological hydrodynamical simulations are essential tools for studying the formation and evolution of galaxies and their central supermassive black holes. While they reproduce many key observed properties of galaxies, their limited volumes have hindered comprehensive studies of the AGN and quasar populations. In this work, we leverage the FLAMINGO simulation suite, focusing on its large $(2.8$… ▽ More Cosmological hydrodynamical simulations are essential tools for studying the formation and evolution of galaxies and their central supermassive black holes. While they reproduce many key observed properties of galaxies, their limited volumes have hindered comprehensive studies of the AGN and quasar populations. In this work, we leverage the FLAMINGO simulation suite, focusing on its large $(2.8$ $\mathrm{Gpc})^3$ volume, to investigate two key observables of quasar activity: the quasar luminosity function (QLF) and quasar clustering. FLAMINGO reproduces the observed QLF at low redshift ($z \lesssim 1$) and for faint quasars ($L_\mathrm{bol} \lesssim 10^{45}$ $\mathrm{erg s^{-1}}$), but significantly underpredicts the abundance of bright quasars at $z \approx 1$-$3$. Introducing a 0.75 dex log-normal luminosity scatter to represent unresolved small-scale variability boosts the number of bright quasars by upscattering lower-luminosity systems, thereby improving agreement with observations at the bright end. A decomposition of the QLF by black hole mass reveals that this boost is primarily driven by low-mass black holes radiating above the Eddington limit. Nevertheless, limitations remain in fully reproducing the rise and decline of the bright quasar population over cosmic time and in matching the black hole masses inferred from quasar spectra. Thanks to FLAMINGO's large volume, we can robustly sample rare, luminous quasars and measure their spatial clustering for $\log_{10} L_\mathrm{bol}/\mathrm{erg s^{-1}} \gtrsim 45.5$. The simulation reproduces the observed clustering across $0 \lesssim z \lesssim 3$, and the reduced luminosity dependence introduced by scatter aligns with observational trends. However, it underpredicts the clustering strength at $z \approx 4$, consistent with other models and possibly reflecting high-redshift observational uncertainties. △ Less

Submitted 12 February, 2026; v1 submitted 28 October, 2025; originally announced October 2025.

Comments: 22 pages, 13 figures. Accepted for publication in MNRAS; in press

QUEST (Quasar Unsupervised Encoder and Synthesis Tool): A machine learning framework to generate quasar spectra

Authors: F. Guarneri, J. T. Schindler, R. A. Meyer, D. Yang, J. F. Hennawi, L. Lucie-Smith, S. E. I. Bosman, F. B. Davies

Abstract: Quasars at the redshift frontier (z > 7.0) are fundamental probes of black hole (BH) growth and evolution but notoriously difficult to identify. At these redshifts, machine learning-based selection methods have proven to be efficient, but require appropriate training sets to express their full potential. Here, we present QUEST, a Variational Auto-Encoder capable of generating realistic quasar spec… ▽ More Quasars at the redshift frontier (z > 7.0) are fundamental probes of black hole (BH) growth and evolution but notoriously difficult to identify. At these redshifts, machine learning-based selection methods have proven to be efficient, but require appropriate training sets to express their full potential. Here, we present QUEST, a Variational Auto-Encoder capable of generating realistic quasar spectra that can be post-processed for generating synthetic photometry and for spectral imputation. We start from the SDSS DR16Q catalogue, pre-process the spectra, and vet the sample to obtain a clean data set. After training the model, we investigate the properties of its latent space to understand whether it has learnt relevant physics. We provide a pipeline to generate photometry from the sampled spectra, compare it with actual quasar photometry, and showcase the capabilities of the model in reconstructing and extending quasar spectra. The trained network faithfully reproduces the input spectrum, both in terms of sample median and variance. By examining the latent space, we find correlations with continuum and bolometric luminosity, BH mass, redshift, continuum slope, and emission line properties. When used to generate photometry, we find results in excellent agreement with the control sample. The model provides satisfactory results in reconstructing emission lines: estimates of the BH mass from the reconstructed spectra are in good agreement with those from the original spectra. Furthermore, when spectra with broad absorption line features are reconstructed, the model successfully interpolates over the absorption systems. Compared with previous work, we find excellent agreement between the spectra sampled from our model and the output of their results. However, QUEST does not require any ad-hoc tuning, and is capable of reproducing the full variety of spectra available in the training set. △ Less

Submitted 27 October, 2025; originally announced October 2025.

Comments: 20 pages (15 main text + 1 references & aknowledgments + 4 appendix). Abstract abridged for Arxiv. Submitted to A&A, comments welcome! We release the tool as a python package, available at: https://github.com/cosmic-dawn-group/QUEST

Journal ref: A&A 709, A241 (2026)

BEES: Quasar lifetime measurements from extended rest-optical emission line nebulae at $z\sim6$

Authors: Dominika Ďurovčíková, Anna-Christina Eilers, Yuzo Ishikawa, Minghao Yue, Marianne Vestergaard, Frederick B. Davies, Jan-Torge Schindler, Xiaohui Fan, Fabrizio Arrigoni Battaia, Marta Volonteri, Robert A. Simcoe, Joseph F. Hennawi, Laura Blecha, Irham T. Andika, Sarah E. I. Bosman, Rebekka Bieri

Abstract: Measurements of quasar lifetimes at high redshift indicate that the earliest billion-solar-mass supermassive black holes (SMBHs) have only been active as luminous quasars for less than a million years. Recently, extended Ly$α$ nebulae around $z\sim6$ quasars have revealed that these short observed lifetimes are unlikely a sightline-dependent effect. However, the interpretation of Ly$α$ emission is… ▽ More Measurements of quasar lifetimes at high redshift indicate that the earliest billion-solar-mass supermassive black holes (SMBHs) have only been active as luminous quasars for less than a million years. Recently, extended Ly$α$ nebulae around $z\sim6$ quasars have revealed that these short observed lifetimes are unlikely a sightline-dependent effect. However, the interpretation of Ly$α$ emission is not straightforward due to its resonant nature. In this work, we use rest-frame optical emission lines, which more directly trace photoionization by the quasar, to unambiguously validate the short line-of-sight quasar lifetimes observed at early cosmic epochs. We use deep James Webb Space Telescope/NIRSpec IFU observations of five $z\sim 6$ quasars with small proximity zones to search for their extended emission line nebulae in H$α$ and [O III]$5007$, and detect extended emission in both emission lines around four quasars in our sample. We then use the light-crossing time of these nebulae to measure quasar lifetimes along transverse sightlines. Using their H$α$ nebulae, we also confirm that recombination is likely the dominant emission mechanism behind their previously detected Ly$α$ nebulae. Our results confirm the existence of high-redshift quasars with extremely short lifetimes, $t_{\rm Q} \lesssim 10^{5}\ {\rm yr}$, hosting billion-solar-mass black holes, indicating that rapid accretion is likely responsible for the assembly of SMBHs in the early Universe. △ Less

Submitted 10 October, 2025; originally announced October 2025.

Comments: 22 pages, 10 figures, 4 tables

A first look at quasar-galaxy clustering at $z\simeq7.3$

Authors: Jan-Torge Schindler, Joseph F. Hennawi, Frederick B. Davies, Sarah E. I. Bosman, Feige Wang, Jinyi Yang, Anna-Christina Eilers, Xiaohui Fan, Koki Kakiichi, Elia Pizzati, Riccardo Nanni

Abstract: We present JWST observations of the environments surrounding two high-redshift quasars, J0252$-$0503 at $z = 7.0$ and J1007$+$2115 at $z = 7.5$, which enable the first constraints on quasar-galaxy clustering at $z \sim 7.3$. Galaxies in the vicinity of the quasars are selected through ground-based and JWST/NIRCam imaging and then spectroscopically confirmed with JWST/NIRSpec using the multi-shutte… ▽ More We present JWST observations of the environments surrounding two high-redshift quasars, J0252$-$0503 at $z = 7.0$ and J1007$+$2115 at $z = 7.5$, which enable the first constraints on quasar-galaxy clustering at $z \sim 7.3$. Galaxies in the vicinity of the quasars are selected through ground-based and JWST/NIRCam imaging and then spectroscopically confirmed with JWST/NIRSpec using the multi-shutter assembly (MSA). Over both fields, we identified 51 $z>5$ galaxies, of which eight are found within a $Δv_{\textrm{LOS}}=\pm1500 \rm{km} \rm{s}^{-1}$ line-of-sight velocity window from the quasars and another eight in the background. The galaxy J0252\_8713, located just $7\,\rm{pkpc}$ and $Δv_{\textrm{LOS}} \approx 360\,\rm{km}\,\rm{s}^{-1}$ from quasar J0252$-$0503, emerges as a compelling candidate for one of the most distant quasar-galaxy mergers. Combining the galaxy discoveries over the two fields, we measure the quasar-galaxy cross-correlation and obtain a correlation length of $r_0^{\rm{QG}}\approx7.6_{-1.6}^{+1.7}\,h^{-1}\,\rm{cMpc}$, based on a power-law model with a fixed slope of $γ_{\rm{QG}} = 2.0$. Under the assumption that quasars and galaxies trace the same underlying dark matter density fluctuations, we infer a minimum dark matter halo mass for $z\simeq7.3$ quasars of $\log_{10}(M_{\textrm{halo, min}}/\textrm{M}_{\odot})= 11.6_{-0.7}^{+0.6}$ in a halo model framework. Compared to measurements from EIGER at $\langle z \rangle = 6.25$ and ASPIRE at $\langle z \rangle = 6.7$ (where $\log_{10}(M_{\textrm{halo, min}}/\textrm{M}_{\odot}) \gtrsim 12.3$), our clustering results provide tentative evidence for a non-monotonic redshift evolution of quasar clustering properties. We further estimate a quasar duty cycle of $f_{\rm{duty}}\approx0.05\%$, consistent with constraints from quasar proximity zones and IGM damping wings. (abridged) △ Less

Submitted 20 March, 2026; v1 submitted 9 October, 2025; originally announced October 2025.

Comments: 16 pages, 7+2 figures; accepted by A&A

Journal ref: A&A 708, A160 (2026)

Using Neural Emulators and Hamiltonian Monte Carlo to constrain the Epoch of Reionization's History with the Ly$α$ Forest Power Spectrum

Authors: Diego González-Hernández, Caitlin Doughty, Molly Wolfson, Joseph F. Hennawi, Zhenyu Jin

Abstract: The Lyman-alpha (Ly$α$) forest at $z \sim 5$ offers a primary probe to constrain the history of the Epoch of Reionization (EoR), retaining thermal and ionization signatures imprinted by the reionization process. In this work, we present a new inference framework based on JAX that combines forward-modeled Ly$α$ forest observables with differentiable neural emulators and Hamiltonian Monte Carlo (HMC… ▽ More The Lyman-alpha (Ly$α$) forest at $z \sim 5$ offers a primary probe to constrain the history of the Epoch of Reionization (EoR), retaining thermal and ionization signatures imprinted by the reionization process. In this work, we present a new inference framework based on JAX that combines forward-modeled Ly$α$ forest observables with differentiable neural emulators and Hamiltonian Monte Carlo (HMC). We construct a dataset of 501 low-resolution simulations generated with user-defined reionization histories and compute a set of 1D Ly$α$ power spectra and model-dependent covariance matrices. We then train two independent neural emulators that achieve sub-percent errors across relevant scales and combine them with HMC to efficiently perform parameter estimation. We validate this framework by applying it to a suite of mock observations, demonstrating that the true parameters are reliably recovered. While this work is limited by the low resolution of the simulations used, our results highlight the potential of this method for inferring the reionization history from high-redshift Ly$α$ forest measurements. Future improvements in our reionization models will further enhance its ability to extract constraints from observational datasets. △ Less

Submitted 16 September, 2025; originally announced September 2025.

Comments: Submitted to MNRAS

Journal ref: Mon Not R Astron Soc (2026)

The Light Echo of a High-Redshift Quasar mapped with Lyman-$α$ Tomography

Authors: Anna-Christina Eilers, Minghao Yue, Jorryt Matthee, Joseph F. Hennawi, Frederick B. Davies, Robert A. Simcoe, Richard Teague, Rongmon Bordoloi, Gabriel Brammer, Yi Kang, Daichi Kashino, Ruari Mackenzie, Rohan P. Naidu, Benjamin Navarrete

Abstract: Ultra-violet (UV) radiation from accreting black holes ionizes the intergalactic gas around early quasars, carving out highly ionized bubbles in their surroundings. Any changes in a quasar's luminosity are therefore predicted to produce outward-propagating ionization gradients, affecting the Lyman-$α$ (Ly$α$) absorption opacity near the quasar's systemic redshift. This "proximity effect" is well-d… ▽ More Ultra-violet (UV) radiation from accreting black holes ionizes the intergalactic gas around early quasars, carving out highly ionized bubbles in their surroundings. Any changes in a quasar's luminosity are therefore predicted to produce outward-propagating ionization gradients, affecting the Lyman-$α$ (Ly$α$) absorption opacity near the quasar's systemic redshift. This "proximity effect" is well-documented in rest-UV quasar spectra but only provides a one-dimensional probe along our line-of-sight. Here we present deep spectroscopic observations with the James Webb Space Telescope (JWST) of galaxies in the background of a super-luminous quasar at $z_{\rm QSO}\approx6.3$, which reveal the quasar's "light echo" with Ly$α$ tomography in the transverse direction. This transverse proximity effect is detected for the first time towards multiple galaxy sightlines, allowing us to map the extent and geometry of the quasar's ionization cone. We obtain constraints on the orientation and inclination of the cone, as well as an upper limit on the obscured solid angle fraction of $f_{\rm obsc}<91\%$. Additionally, we find a timescale of the quasar's UV radiation of $t_{\rm QSO}=10^{5.6^{+0.1}_{-0.3}}$ years, which is significantly shorter than would be required to build up the central supermassive black hole (SMBH) with conventional growth models, but is consistent with independent measurements of the quasars' duty cycle. Our inferred obscured fraction disfavors a scenario where short quasar lifetimes can be explained exclusively by geometric obscuration, and instead supports the idea that radiatively inefficient accretion or growth in initially heavily enshrouded cocoons plays a pivotal role in early SMBH growth. Our results pave the way for novel studies of quasars' ionizing geometries and radiative histories at early cosmic times. △ Less

Submitted 5 September, 2025; originally announced September 2025.

Comments: Submitted to AAS Journals

First constraints on the local ionization topology in front of two quasars at z ~ 7.5

Authors: Timo Kist, Joseph F. Hennawi, Frederick B. Davies, Eduardo Bañados, Sarah E. I. Bosman, Zheng Cai, Anna-Christina Eilers, Xiaohui Fan, Zoltán Haiman, Hyunsung D. Jun, Yichen Liu, Jinyi Yang, Feige Wang

Abstract: Thus far, Lyman-$α$ damping wings towards quasars have been used to probe the \textit{global} ionization state of the foreground intergalactic medium (IGM). A new parameterization has demonstrated that the damping wing signature also carries \textit{local} information about the distribution of neutral hydrogen (HI) in front of the quasar before it started shining. Leveraging a recently introduced… ▽ More Thus far, Lyman-$α$ damping wings towards quasars have been used to probe the \textit{global} ionization state of the foreground intergalactic medium (IGM). A new parameterization has demonstrated that the damping wing signature also carries \textit{local} information about the distribution of neutral hydrogen (HI) in front of the quasar before it started shining. Leveraging a recently introduced Bayesian \texttt{JAX}-based Hamiltonian Monte Carlo (HMC) inference framework, we derive constraints on the Lorentzian-weighted HI column density $N_\mathrm{HI}^\mathrm{DW}$, the quasar's distance $r_\mathrm{patch}$ to the first neutral patch and its lifetime $t_\mathrm{Q}$ based on JWST/NIRSpec spectra of the two $z \sim 7.5$ quasars J1007+2115 and J1342+0928. After folding in model-dependent topology information, we find that J1007+2115 (and J1342+0928) is most likely to reside in a $\langle x_\mathrm{HI} \rangle = 0.32_{-0.20}^{+0.22}$ ($0.58_{-0.23}^{+0.23}$) neutral IGM while shining for a remarkably short lifetime of $\log_{10} t_\mathrm{Q} /\mathrm{yr} = 4.14_{-0.18}^{+0.74}$ (an intermediate lifetime of $5.64_{-0.43}^{+0.25}$) along a sightline with $\log_{10} N_\mathrm{HI}^\mathrm{DW} /\mathrm{cm}^{-2} = 19.70_{-0.86}^{+0.35}$ ($20.24_{-0.22}^{+0.25}$) and $r_\mathrm{patch} = 28.9_{-14.4}^{+54.0} \,\mathrm{cMpc}$ ($10.9_{-5.9}^{+5.6} \,\mathrm{cMpc}$). In light of the potential presence of local absorbers in the foreground of J1342+0928 as has been recently suggested, we also demonstrate how the Lorentzian-weighted column density $N_\mathrm{HI}^\mathrm{DW}$ provides a natural means for quantifying their contribution to the observed damping wing signal. △ Less

Submitted 29 August, 2025; originally announced August 2025.

Comments: submitted to MNRAS

Journal ref: Mon Not R Astron Soc (2025)

Inferring local quasar IGM damping wing constraints

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

Abstract: Lyman-$α$ damping wings towards quasars are a highly sensitive probe of the neutral hydrogen (HI) content in the foreground intergalactic medium (IGM), not only constraining the global timing of reionization but also the \textit{local} ionization topology near the quasar. Near-optimal extraction of this information is possible with the help of two recently introduced reionization model-independent… ▽ More Lyman-$α$ damping wings towards quasars are a highly sensitive probe of the neutral hydrogen (HI) content in the foreground intergalactic medium (IGM), not only constraining the global timing of reionization but also the \textit{local} ionization topology near the quasar. Near-optimal extraction of this information is possible with the help of two recently introduced reionization model-independent summary statistics of the HI distribution in the IGM \textit{before} the quasar started shining, complemented with the quasar's lifetime encoding the effect of its ionizing radiation as a third parameter. We introduce a fully Bayesian JAX-based Hamiltonian Monte Carlo (HMC) inference framework that allows us to jointly reconstruct the quasar's unknown continuum and constrain these local damping wing statistics. We put forward a probabilistic framework that allows us to tie these local constraints to any specific reionization model and obtain model-dependent constraints on the global timing of reionization. We demonstrate that we are able to constrain the (Lorentzian-weighted) HI column density in front of the quasar to a precision of $0.69_{-0.30}^{+0.06}\,\mathrm{dex}$ and its original distance to the first neutral patch before the quasar started shining to $31.4_{-28.1}^{+10.7}\,\mathrm{cMpc}$ (if a noticeable damping wing is present in the spectrum), extracting hitherto unused local information from the IGM damping wing imprint. Once tied to a specific reionization model, we find that the statistical fidelity of our constraints on the global IGM neutral fraction and the lifetime of the quasar improves, while retaining the same precision as achieved by pipelines that infer these parameters directly. △ Less

Submitted 29 August, 2025; originally announced August 2025.

Comments: submitted to MNRAS

Journal ref: Mon Not R Astron Soc (2026)

When Stars Mimic Monsters: Spectral Evidence for an $η$ Carinae-like Giant Eruption in SBS 0335$-$052 E

Authors: Zixuan Peng, Crystal L. Martin, Jiamu Huang, Nikolaus Z. Prusinski, Chenliang Huang, Zhuyun Zhuang, Yuan Li, Tin Long Sunny Wong, Jiayang Yang, Joseph F. Hennawi

Abstract: SBS 0335$-$052 E is an extremely low-metallicity ($Z\sim0.04\,Z_{\odot}$) blue compact dwarf galaxy. An active galactic nucleus has been proposed to explain the broad H$α$ emission and near-infrared (NIR) time variability in super star clusters 1 and 2 (SSCs 1&2). However, Peng et al. discovered broad wings in the forbidden [O III] $\lambda5007$ emission (up to $\sim5\,000\,\rm{km\,s^{-1}}$), chal… ▽ More SBS 0335$-$052 E is an extremely low-metallicity ($Z\sim0.04\,Z_{\odot}$) blue compact dwarf galaxy. An active galactic nucleus has been proposed to explain the broad H$α$ emission and near-infrared (NIR) time variability in super star clusters 1 and 2 (SSCs 1&2). However, Peng et al. discovered broad wings in the forbidden [O III] $\lambda5007$ emission (up to $\sim5\,000\,\rm{km\,s^{-1}}$), challenging the broad-line region interpretation. We present new KCWI/KCRM integral-field spectroscopy to directly compare spectra across multiple SSCs. The nebula surrounding SSCs 1&2 shows unique features. The Ly$β$-pumped O I $\lambda8446$ emission constrains $τ_{\rm\,Lyα}\sim10^8$. Multiple ionization states of iron are detected from Fe$^{+}$ to Fe$^{+4}$. Stellar photoionization models can reproduce the [Fe III]/[Fe II] and [Fe IV]/[Fe III] line ratios at high density ($n_e\sim10^6\,\rm{cm^{-3}}$), but they fail to account for most of the [Fe V] emission. The broad H$α$ wings exhibit an exponential profile; the asymmetric wings extend from $\sim-5\,000\,\rm{km\,s^{-1}}$ to $\sim10\,000\,\rm{km\,s^{-1}}$. Thomson scattering in a radially expanding medium provides a good fit with $v_w\sim200\,\rm{km\,s^{-1}}$, optical depth $τ_e\sim10$, and an outer to inner radius of 10. Enhanced N/O and potentially depleted Fe/O ratios are consistent with CNO-cycled ejecta from massive stars and with dust formation, respectively. We propose that mass loss from a massive star interacting with its circumstellar medium drives a shock that powers the NIR variability, the luminous X-ray point source, and the [Fe V] emission. If confirmed, the proposed stellar eruption would be a distant example of an $η$ Carinae-like giant eruption, and the first in an ultra-low metallicity environment. △ Less

Submitted 27 March, 2026; v1 submitted 5 August, 2025; originally announced August 2025.

Comments: 28 pages, 13 figures, 3 tables. Accepted for publication in ApJ

Journal ref: ApJ 1001, 218 (2026)

Luminous Mid-IR Selected Obscured Quasars at Cosmic Noon in SDSS Stripe82 II: Spectroscopic Diversity and Broad H$α$ Emissions

Authors: Ben Wang, Yuzo Ishikawa, Joseph F. Hennawi, Zheng Cai, Gordon T. Richards, Nadia L. Zakamska, Daming Yang, Jan-Torge Schindler

Abstract: We present a multiwavelength spectroscopic survey of 23 luminous mid-infrared-selected Type-2 quasars at redshifts z = 0.88 to 3.49. The targets were selected in the SDSS Stripe 82 field based on their bright WISE W4 detections (flux > 5 mJy) and extremely faint or red optical counterparts (e.g., r > 23 or r - W4 > 8.4), designed to identify heavily obscured quasars. Deep near-infrared (Gemini/GNI… ▽ More We present a multiwavelength spectroscopic survey of 23 luminous mid-infrared-selected Type-2 quasars at redshifts z = 0.88 to 3.49. The targets were selected in the SDSS Stripe 82 field based on their bright WISE W4 detections (flux > 5 mJy) and extremely faint or red optical counterparts (e.g., r > 23 or r - W4 > 8.4), designed to identify heavily obscured quasars. Deep near-infrared (Gemini/GNIRS) and optical (Keck/LRIS and KCWI) spectroscopy confirm 23 out of 24 candidates as Type-2 quasars in this redshift range, including 12 objects at z > 2. The spectra exhibit strong rest-frame UV and optical emission lines (Ly-alpha, C IV, [O III], H-alpha) with a wide range of line widths, indicating significant spectral diversity. Approximately one-third of the sample (8 of 23) shows broad H-alpha emission (FWHM > 2000 km/s) despite their Type-2 classification, while the rest have only narrow lines (FWHM < 2000 km/s) characteristic of classical obscured quasars. Notably, these broad-line Type-2 quasars share similar spectral energy distributions with the JWST-discovered "little red dot" (LRD) AGNs, suggesting that our sample could be lower-redshift analogues of the heavily obscured broad-line AGNs uncovered by JWST. We also find that the [O III] 5007 angstrom emission is relatively weak for their high bolometric luminosities, deviating from trends seen in lower-redshift Type-2 QSOs. A new composite spectrum for Type-2 QSOs is built using our sample. Overall, our results demonstrate that mid-IR selection efficiently uncovers a diverse population of obscured quasars and that spectroscopic follow-up is crucial for revealing their true nature. This study provides new insights into heavily obscured SMBH growth at cosmic noon and bridges the gap to the obscured AGN populations now being revealed by JWST. △ Less

Submitted 4 August, 2025; originally announced August 2025.

Escape fractions from unattenuated Ly$α$ emitters around luminous $z>6$ quasars

Authors: Minghao Yue, Anna-Christina Eilers, Jorryt Matthee, Rohan P. Naidu, Rongmon Bordoloi, Frederick B. Davies, Joseph F. Hennawi, Daichi Kashino, Ruari Mackenzie, Robert A. Simcoe

Abstract: Ionized proximity zones around luminous quasars provide a unique laboratory to characterize the Ly$α$ emission lines from $z>6$ galaxies without significant attenuation from the intergalactic medium (IGM). However, Ly$α$ line measurements for galaxies within high-redshift quasars' proximity zones have been rare so far. Here we present deep spectroscopic observations obtained with the NIRSpec/MSA i… ▽ More Ionized proximity zones around luminous quasars provide a unique laboratory to characterize the Ly$α$ emission lines from $z>6$ galaxies without significant attenuation from the intergalactic medium (IGM). However, Ly$α$ line measurements for galaxies within high-redshift quasars' proximity zones have been rare so far. Here we present deep spectroscopic observations obtained with the NIRSpec/MSA instrument on the James Webb Space Telescope (JWST) of galaxies in two $z>6$ quasar fields. We measure the Ly$α$ line fluxes for 50 galaxies at $6<z<7$ with UV absolute magnitude $M_\text{UV}<-19$ (median $M_\text{UV}=-19.97$), among which 15 are located near the luminous quasars, i.e. within $Δv<2500\rm\,km\,s^{-1}$. We find that galaxies near the quasars show significant flux bluewards of the systemic Ly$α$ wavelength, and have higher Ly$α$ equivalent width compared to galaxies at similar redshifts that are not located within the quasars' environment. Our result indicates little or no redshift evolution for the Ly$α$-emitter fraction from $z\sim6.4$ to $z\sim5$. Leveraging the low IGM opacity in the quasars' vicinity, we evaluate the Ly$α$ escape fraction ($f_\text{esc}^{\text{Ly}α}$) of high-redshift galaxies. Our analysis suggests that galaxies at $\langle z\rangle\approx6.4$ have an average $f_\text{esc}^{\text{Ly}α}=0.14\pm0.04$. This value is consistent with reionization models where the Lyman continuum escape fraction is low $(f_\text{esc}^\text{LyC}\lesssim0.1)$ for luminous galaxies, and where the most luminous galaxies have only a minor contribution to the total ionizing photon budget. △ Less

Submitted 22 September, 2025; v1 submitted 7 July, 2025; originally announced July 2025.

Comments: Accepted by ApJL

Homogeneous measurements of proximity zone sizes for 59 quasars in the Epoch of Reionization

Authors: Silvia Onorato, Joseph F. Hennawi, Elia Pizzati, Bram P. Venemans, Anna-Christina Eilers

Abstract: The overionized regions surrounding high-redshift quasars, known as proximity zones, provide a window into the interaction between supermassive black holes (SMBHs) and the intergalactic medium (IGM) during the epoch of reionization (EoR). We present new homogeneous measurements of proximity zone sizes ($R_{\mathrm{p}}$) for a sample of $59$ quasars spanning redshifts $5.77 \leq z \leq 7.54$ (media… ▽ More The overionized regions surrounding high-redshift quasars, known as proximity zones, provide a window into the interaction between supermassive black holes (SMBHs) and the intergalactic medium (IGM) during the epoch of reionization (EoR). We present new homogeneous measurements of proximity zone sizes ($R_{\mathrm{p}}$) for a sample of $59$ quasars spanning redshifts $5.77 \leq z \leq 7.54$ (median $z = 6.59$). For $15$ of these sources, we measure $R_{\mathrm{p}}$ for the first time. The quasars in our catalog have absolute magnitudes at rest-frame $1450$ Å in the range $-29.13 \leq M_{1450} \leq -25.20$ (median $M_{1450} \simeq -26.49$), providing one of the most extensive data sets for exploring $R_{\mathrm{p}}$ at these epochs. The distribution of $R_{\mathrm{p}}$ values shows a large scatter at fixed redshift and luminosity, likely reflecting variations in quasar lifetimes ($t_{\mathrm{Q}}$), IGM density fluctuations, and IGM neutral fraction. We fit a bivariate power-law model to a large sample of $100$ objects to study the dependence of $R_{\mathrm{p}}$ with both $M_{1450}$ and $z$: we find that the evolution of $R_{\mathrm{p}}$ with luminosity is in agreement with the models ($R_{\mathrm{p}} \propto 10^{-0.4 M_{1450}/2.87}$), while the evolution of $R_{\mathrm{p}}$ with $z$ is steeper than previous works ($R_{\mathrm{p}} \propto (1+z)^{-2.44}$). We identify $13$ quasars with small proximity zone size, defined using the residuals of our fit. In all cases, except for J2211$-$6320, we rule out the presence of associated dense absorbers that prematurely truncate $R_{\mathrm{p}}$, and suggest a short $t_{\mathrm{Q}}$ ($\lesssim 10^4$ yr) as a possible explanation for their small proximity zone sizes. △ Less

Submitted 3 April, 2026; v1 submitted 14 May, 2025; originally announced May 2025.

Comments: 33 pages, 15 figures, published on MNRAS

Quasar lifetime measurements from extended Ly$α$ nebulae at $z\sim 6$

Authors: Dominika Ďurovčíková, Anna-Christina Eilers, Romain A. Meyer, Emanuele Paolo Farina, Eduardo Bañados, Frederick B. Davies, Joseph F. Hennawi, Chiara Mazzucchelli, Robert A. Simcoe, Fabian Walter

Abstract: The existence of billion-solar-mass black holes hosted in luminous quasars within the first gigayear of cosmic history poses a challenge to our understanding of supermassive black hole (SMBH) growth. The problem is further exacerbated by the very short quasar lifetimes of $t_{\rm Q}\lesssim 10^6$ years, as derived from the extent of their proximity zone (PZ) sizes observed in the quasars' rest-UV… ▽ More The existence of billion-solar-mass black holes hosted in luminous quasars within the first gigayear of cosmic history poses a challenge to our understanding of supermassive black hole (SMBH) growth. The problem is further exacerbated by the very short quasar lifetimes of $t_{\rm Q}\lesssim 10^6$ years, as derived from the extent of their proximity zone (PZ) sizes observed in the quasars' rest-UV spectra. However, the quasar lifetime estimates based on the extents of the proximity zones may be underestimated, as time-variable obscuration effects might have limited the quasars' emission along our sightline in the past. In this work, we present independent quasar lifetime measurements for six quasars at $z \sim 6$ leveraging the extended nebular emission perpendicular to our line-of-sight. We use observations from the Very Large Telescope/Multi-Unit Spectroscopic Explorer (MUSE) to search for extended Ly$α$ emission in the circumgalactic medium around quasars with small proximity zones and estimate their lifetimes as the light travel time between the SMBH and the outer edge of the nebula. We find agreement between the independent lifetime estimates. For one object we find a proximate absorption system prematurely truncating the extent of the quasar's proximity zone, which thus results in an expected discrepancy between the lifetime estimates. Our results provide further evidence that the quasars' current accretion episode has only recently begun, challenging our models of SMBH growth. △ Less

Submitted 8 September, 2025; v1 submitted 30 April, 2025; originally announced May 2025.

Comments: 15 pages, 8 figures

Journal ref: ApJ 990 174 (2025)

A local, topology-independent parameterization of quasar IGM damping wings

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

Abstract: Lyman-$α$ damping wings towards quasars provide a unique probe of reionization because their strength correlates strongly with the global volume-averaged neutral hydrogen (HI) fraction of the intergalactic medium (IGM). Cosmic variance in the IGM, however, is a major source of stochasticity since the local neutral environment around a quasar varies significantly even at fixed global neutral fracti… ▽ More Lyman-$α$ damping wings towards quasars provide a unique probe of reionization because their strength correlates strongly with the global volume-averaged neutral hydrogen (HI) fraction of the intergalactic medium (IGM). Cosmic variance in the IGM, however, is a major source of stochasticity since the local neutral environment around a quasar varies significantly even at fixed global neutral fraction. We show that the IGM damping wing carries additional information about this local ionization topology, unexploited by current analysis frameworks. We introduce a set of two new physically motivated summary statistics encoding the local information about the HI distribution in the IGM before it is altered by ionization radiation from the quasar, encompassing 1) the HI column density, weighted by a Lorentzian profile mimicking the frequency dependence of the Lyman-$α$ cross section, and 2) the distance from the quasar to the first neutral patch. This description, when combined with the quasar's lifetime as a third parameter, reduces the IGM transmission scatter in the damping wing region of the spectrum to $\lesssim 1\,\%$ across the full range of physical parameter space. We introduce a simple procedure for generating synthetic HI sightlines around quasars and demonstrate that the resulting damping wing profiles are statistically indistinguishable from a realistic reionization topology. This opens the door for optimally extracting the salient local information encoded in the imprint in a model-independent fashion. In the context of a specific reionization model, measurements of these local parameters can be translated into constraints on the global timing of reionization, but in addition, they provide information about the reionization topology, hitherto unused. A marginally modified version of our framework can also be employed in the context of damping wings towards galaxies. △ Less

Submitted 20 April, 2025; originally announced April 2025.

Comments: Submitted to MNRAS

Journal ref: Mon Not R Astron Soc (2025) 2316-2339

JWST ASPIRE: How Did Galaxies Complete Reionization? Evidence for Excess IGM Transmission around ${\rm [O\,{\scriptstyle III}]}$ Emitters during Reionization

Authors: Koki Kakiichi, Xiangyu Jin, Feige Wang, Romain A. Meyer, Enrico Garaldi, Sarah E. I. Bosman, Frederick B. Davies, Xiaohui Fan, Maxime Trebitsch, Jinyi Yang, Eduardo Bañados, Jaclyn B. Champagne, Anna-Christina Eilers, Joseph F. Hennawi, Fengwu Sun, Yunjing Wu, Siwei Zou, Rahul Kannan, Aaron Smith, George D. Becker, Valentina D'Odorico, Thomas Connor, Weizhe Liu, Klaudia Protušová, Fabian Walter , et al. (1 additional authors not shown)

Abstract: The spatial correlation between galaxies and the Ly$α$ forest provides insights into how galaxies reionized the Universe. Here, we present initial results on the spatial cross-correlation between [OIII] emitters and Ly$α$ forest at 5.4<z<6.5 from the JWST ASPIRE NIRCam/F356W Grism Spectroscopic Survey in z>6.5 QSO fields. Using data from five QSO fields, we find $2σ$ evidence for excess Ly$α$ fore… ▽ More The spatial correlation between galaxies and the Ly$α$ forest provides insights into how galaxies reionized the Universe. Here, we present initial results on the spatial cross-correlation between [OIII] emitters and Ly$α$ forest at 5.4<z<6.5 from the JWST ASPIRE NIRCam/F356W Grism Spectroscopic Survey in z>6.5 QSO fields. Using data from five QSO fields, we find $2σ$ evidence for excess Ly$α$ forest transmission at ~20-40 cMpc around [OIII] emitters at z=5.86, indicating that [OIII] emitters reside within a highly ionized IGM. At smaller scales, the Ly$α$ forest is preferentially absorbed, suggesting gas overdensities around [OIII] emitters. Comparing with models including THESAN simulations, we interpret the observed cross-correlation as evidence for significant large-scale fluctuations of the IGM and the late end of reionization at z<6, characterized by ionized bubbles over 50 cMpc around [OIII] emitters. The required UV background necessitates an unseen population of faint galaxies around the [OIII] emitters. Furthermore, we find that the number of observed [OIII] emitters near individual transmission spikes is insufficient to sustain reionization in their surroundings, even assuming all [OIII] emitters harbour AGN with 100 % LyC escape fractions. Despite broad agreement, a careful analysis of ASPIRE and THESAN, using the observed host halo mass from the clustering of [OIII] emitters, suggests that the simulations underpredict the observed excess IGM transmission around [OIII] emitters, challenging our model of reionization. Potential solutions include larger ionized bubbles at z<6, more enhanced large-scale UV background or temperature fluctuations of the IGM, and possibly a patchy early onset of reionization at z>10. Current observational errors are dominated by cosmic variance, meaning future analyses of more QSO fields from JWST will improve the results. △ Less

Submitted 10 March, 2025; originally announced March 2025.

Comments: 31 pages, 22 figures, submitted to the Open Journal of Astrophysics

Luminous Mid-IR Selected Type-2 Quasars at Cosmic Noon in SDSS Stripe82 I: Selection, Composite Photometry, and Spectral Energy Distributions

Authors: Ben Wang, Joseph F. Hennawi, Zheng Cai, Gordon T. Richards, Jan-Torge Schindler, Nadia L. Zakamska, Yuzo Ishikawa, Hollis B. Akins, Zechang Sun

Abstract: We analyze 23 spectroscopically confirmed Type-2 quasars (QSOs) selected from the WISE 22$\rm μ$m band in the SDSS Stripe 82 region, focusing on their multi-band photometry and spectral energy distributions (SEDs). These objects were selected to be IR-luminous ($\rm flux_{W4} > 5mJy$, i.e., $12.62 < W4 < 14.62 \rm\ AB \, magnitude$), optically faint ($r > 23$) or with red color ($r - W4 >8.38$). G… ▽ More We analyze 23 spectroscopically confirmed Type-2 quasars (QSOs) selected from the WISE 22$\rm μ$m band in the SDSS Stripe 82 region, focusing on their multi-band photometry and spectral energy distributions (SEDs). These objects were selected to be IR-luminous ($\rm flux_{W4} > 5mJy$, i.e., $12.62 < W4 < 14.62 \rm\ AB \, magnitude$), optically faint ($r > 23$) or with red color ($r - W4 >8.38$). Gemini/GNIRS observations were conducted for all 24 candidates, and 18/24 were also observed with Keck/LRIS. The observations confirm 23 to be real Type-2 QSOs in the redshift range $0.88 - 2.99$ (12 are at $z>2$). We collect multi-band photometry and conduct SED fitting. The composite photometry probes the wavelength from 0.1$\rm μ$m to 10$\rm μ$m at the rest frame. The IR emission is dominated by dust torus implying an average torus luminosity for the sample of $L_{\rm torus} 10^{46.84} \rm erg/s$. The origin of the rest-UV/optical light is not definitive, but we present three possible scenarios: scattered light, stellar emission, and the reddened accretion disk. Assuming an obscured:unobscured ratio of approximately 1:1, our targets have $L_{\rm bol} = 10^{46.28} \rm erg \,s^{-1} - 10^{47.49} \rm erg \,s^{-1}$ and around SMBH masses $\rm 10^{8.18} M_{\odot} - 10^{9.39} M_{\odot}$, assuming they accreate at the Eddington limit. Compared to previous Type-2 AGN SEDs, our targets have a brighter dust torus and redder optical-IR color. By comparing the SED to the results from JWST `little red dots' (LRDs), we find that these IR-selected Type-2 QSOs have similar SED shapes to the LRDs. This pilot Type-2 QSO survey demonstrates that mid-IR selection is an efficient way to find luminous Type-2 QSOs at $z>2$. Finally, the composite photometry and Type-2 QSOs SED model generated by this sample provide a guide for finding more Type-2 QSOs at higher redshift. △ Less

Submitted 22 May, 2025; v1 submitted 23 January, 2025; originally announced January 2025.

A unique window into the Epoch of Reionisation: A double-peaked Lyman-$α$ emitter in the proximity zone of a quasar at $z\sim 6.6$

Authors: Klaudia Protušová, Sarah E. I. Bosman, Feige Wang, Romain A. Meyer, Jaclyn B. Champagne, Frederick B. Davies, Anna-Christina Eilers, Xiaohui Fan, Joseph F. Hennawi, Xiangyu Jin, Hyunsung D. Jun, Koki Kakiichi, Zihao Li, Weizhe Liu, Jinyi Yang

Abstract: We present a study of a double-peaked Ly$α$ emitter, named LAE-11, found in the proximity zone of QSO J0910-0414 at $z\sim6.6$. We use a combination of deep photometric data from Subaru Telescope, HST, and JWST with spectroscopic data from Keck/DEIMOS, NIRCam WFSS and NIRSpec MSA to characterise the ionising and general properties of the galaxy, as well as the quasar environment surrounding it. We… ▽ More We present a study of a double-peaked Ly$α$ emitter, named LAE-11, found in the proximity zone of QSO J0910-0414 at $z\sim6.6$. We use a combination of deep photometric data from Subaru Telescope, HST, and JWST with spectroscopic data from Keck/DEIMOS, NIRCam WFSS and NIRSpec MSA to characterise the ionising and general properties of the galaxy, as well as the quasar environment surrounding it. We detect Ly$α$, H$β$, [OIII] doublet, and H$α$ emission lines in the various spectral datasets. The presence of a double-peaked Ly$α$ line in the spectrum allows us to characterise the opening angle and lifetime of the QSO as $θ_Q>49.62$° and $t_Q > 3.8\times10^5$ years. LAE-11 is a fairly bright (M$_\rm{UV} = -19.84^{+0.14}_{-0.16}$), blue galaxy with a UV slope of $β= -2.61^{+0.06}_{-0.08}$. Since the galaxy is located in a quasar-ionised region, we have a unique opportunity to measure the escape fraction of Lyman Continuum photons using the un-attenuated double-peaked Ly$α$ emission profile and its equivalent width at such high redshift. We employ diagnostics which do not rely on the detection of Ly$α$ for comparison, and find that all tracers of ionising photon leakage agree within 1$σ$ uncertainty. We measure a moderate escape of Lyman Continuum photons from LAE-11 of $f_\rm{esc}^\rm{LyC} = (8 - 33)\%$. Detections of both H$α$ and H$β$ emission lines allow for separate measurements of the ionising photon production efficiency, resulting with $\log(ξ_\rm{ion}/\rm{Hz~erg^{-1}}) = 25.59\pm0.08$ and $25.65\pm0.09$, for H$α$ and H$β$, respectively. The total ionising output of LAE-11, $\log(f_\rm{esc}^\rm{LyC}ξ_\rm{ion, Hα}/\rm{Hz~erg^{-1}}) = 24.85^{+0.24}_{-0.34}$, is higher than the value of $24.3 - 24.8$ which is traditionally assumed as needed to drive Reionisation forward. △ Less

Submitted 19 December, 2024; v1 submitted 16 December, 2024; originally announced December 2024.

Comments: Submitted to A&A. Abstract abridged for arxiv

Journal ref: A&A 700, A218 (2025)

A SPectroscopic survey of biased halos In the Reionization Era (ASPIRE): Spectroscopically Complete Census of Obscured Cosmic Star Formation Rate Density at $z=4-6$

Authors: Fengwu Sun, Feige Wang, Jinyi Yang, Jaclyn B. Champagne, Roberto Decarli, Xiaohui Fan, Eduardo Bañados, Zheng Cai, Luis Colina, Eiichi Egami, Joseph F. Hennawi, Xiangyu Jin, Hyunsung D. Jun, Yana Khusanova, Mingyu Li, Zihao Li, Xiaojing Lin, Weizhe Liu, Romain A. Meyer, Maria A. Pudoka, George H. Rieke, Yue Shen, Wei Leong Tee, Bram Venemans, Fabian Walter , et al. (3 additional authors not shown)

Abstract: We present a stringent measurement of the dust-obscured star-formation rate density (SFRD) at $z=4-6$ from the ASPIRE JWST Cycle-1 medium and ALMA Cycle-9 large program. We obtained JWST/NIRCam grism spectroscopy and ALMA 1.2-mm continuum map along 25 independent quasar sightlines, covering a total survey area of $\sim$35 arcmin$^2$ where we search for dusty star-forming galaxies (DSFGs) at… ▽ More We present a stringent measurement of the dust-obscured star-formation rate density (SFRD) at $z=4-6$ from the ASPIRE JWST Cycle-1 medium and ALMA Cycle-9 large program. We obtained JWST/NIRCam grism spectroscopy and ALMA 1.2-mm continuum map along 25 independent quasar sightlines, covering a total survey area of $\sim$35 arcmin$^2$ where we search for dusty star-forming galaxies (DSFGs) at $z = 0 - 7$. We identify eight DSFGs in seven fields at $z=4-6$ through the detection of H$α$ or [O III] $λ$5008 lines, including fainter lines such as H$β$, [O III] $λ$4960, [N II] $λ$6585, [S II] $λλ$6718,6733 for six sources. With this spectroscopically complete DSFG sample at $z=4-6$ and negligible impact from cosmic variance (shot noise), we measure the infrared luminosity function (IRLF) down to $L_\mathrm{IR} \sim 2\times10^{11}$ $L_\odot$. We find flattening of IRLF at $z=4-6$ towards the faint end (power-law slope $α= 0.59_{-0.45}^{+0.39}$). We determine the dust-obscured cosmic SFRD at this epoch as $\log[ρ_\mathrm{SFR,IR} / (\mathrm{M}_\odot\,\mathrm{yr}^{-1}\,\mathrm{Mpc}^{-3})] = -1.52_{-0.13}^{+0.14}$. This is significantly higher than previous determination using ALMA data in the Hubble Ultra Deep Field, which is void of DSFGs at $z=4-6$ because of strong cosmic variance (shot noise). We conclude that the majority ($66\pm7$%) of cosmic star formation at $z \sim 5$ is still obscured by dust. We also discuss the uncertainty of SFRD propagated from far-IR spectral energy distribution and IRLF at the bright end, which will need to be resolved with future ALMA and JWST observations. △ Less

Submitted 9 December, 2024; originally announced December 2024.

Comments: 27 pages, 13 figures, 3 tables. Accepted for publication in ApJ

Extreme Deconvolution Reimagined: Conditional Densities via Neural Networks and an Application in Quasar Classification

Authors: Yi Kang, Joseph F. Hennawi, Jan-Torge Schindler, John Tamanas, Riccardo Nanni

Abstract: Density estimation is a fundamental problem that arises in many areas of astronomy, with applications ranging from selecting quasars using color distributions to characterizing stellar abundances. Astronomical observations are inevitably noisy; however, the density of a noise-free feature is often the desired outcome. The extreme-deconvolution (XD) method can be used to deconvolve the noise and ob… ▽ More Density estimation is a fundamental problem that arises in many areas of astronomy, with applications ranging from selecting quasars using color distributions to characterizing stellar abundances. Astronomical observations are inevitably noisy; however, the density of a noise-free feature is often the desired outcome. The extreme-deconvolution (XD) method can be used to deconvolve the noise and obtain noise-free density estimates by fitting a mixture of Gaussians to data where each sample has non-identical (heteroscedastic) Gaussian noise. However, XD does not generalize to cases where some feature dimensions have highly non-Gaussian distribution, and no established method exists to overcome this limitation. We introduce a possible solution using neural networks to perform Gaussian mixture modeling of the Gaussian-like dimensions conditioned on those non-Gaussian features. The result is the CondXD algorithm, a generalization of XD that performs noise-free conditional density estimation. We apply CondXD to a toy model and find that it is more accurate than other approaches. We further test our method on a real-world high redshift quasar versus contaminant classification problem. Specifically, we estimate noise-free densities in flux-ratio (i.e., color) space for contaminants, conditioned on their magnitude. Our results are comparable to the existing method, which divides the samples into magnitude bins and applies XD separately in each bin, and our method is approximately ten times faster. Overall, our method has the potential to significantly improve estimating conditional densities and enable new discoveries in astronomy. △ Less

Submitted 3 December, 2024; originally announced December 2024.

A Little Red Dot at $\mathbf{z=7.3}$ within a Large Galaxy Overdensity

Authors: Jan-Torge Schindler, Joseph F. Hennawi, Frederick B. Davies, Sarah E. I. Bosman, Ryan Endsley, Feige Wang, Jinyi Yang, Aaron J. Barth, Anna-Christina Eilers, Xiaohui Fan, Koki Kakiichi, Michael Maseda, Elia Pizzati, Riccardo Nanni

Abstract: The nature of "Little Red Dots" and their relation to other forms of accreting supermassive black holes remain an open question. Here we report the discovery of a Little Red Dot at $z=7.3$. It is attenuated by moderate amounts of dust, $A_V = {2.79}\,\textrm{mag}$, with an intrinsic bolometric luminosity of $10^{46.6}\,\textrm{erg}\,\textrm{s}^{-1}$ and a SMBH mass of… ▽ More The nature of "Little Red Dots" and their relation to other forms of accreting supermassive black holes remain an open question. Here we report the discovery of a Little Red Dot at $z=7.3$. It is attenuated by moderate amounts of dust, $A_V = {2.79}\,\textrm{mag}$, with an intrinsic bolometric luminosity of $10^{46.6}\,\textrm{erg}\,\textrm{s}^{-1}$ and a SMBH mass of $5\times10^8\,\textrm{M}_\odot$. Most notably, this object is embedded in an overdensity of eight nearby galaxies, allowing us to calculate a spectroscopic estimate of the clustering of galaxies around Little Red Dots. We find a Little Red Dot-galaxy cross-correlation length of $r_0\!=\!8\pm2\,\textrm{h}^{-1}\,\textrm{cMpc}$, comparable to that of $z\!\sim\!6$ UV-luminous quasars. The resulting estimate of their minimum dark matter halo mass of $\log_{10}(M_{\textrm{halo, min}}/\textrm{M}_{\odot})= 12.0_{-1.0}^{+0.8}$ indicates that nearly all halos above this mass must host actively accreting SMBHs at $z\approx7$, in strong contrast with the far smaller duty cycle of luminous quasars ($<1\%$). Our results, taken at face value, motivate a picture in which SMBHs in Little Red Dot phases could serve as the obscured precursors of UV-luminous quasars, which provides a natural explanation for the short UV-luminous lifetimes inferred from both quasar clustering and quasar proximity zones. △ Less

Submitted 19 September, 2025; v1 submitted 18 November, 2024; originally announced November 2024.

Comments: The article is published in Nature Astronomy as Schindler, JT., Hennawi, J.F., Davies, F.B. et al. A little red dot at z = 7.3 within a large galaxy overdensity. Nat Astron (2025). https://doi.org/10.1038/s41550-025-02660-1

Neural network emulator to constrain the high-$z$ IGM thermal state from Lyman-$α$ forest flux auto-correlation function

Authors: Zhenyu Jin, Molly Wolfson, Joseph F. Hennawi, Diego González-Hernández

Abstract: We present a neural network emulator to constrain the thermal parameters of the intergalactic medium (IGM) at $\displaystyle{5.4}\le{z}\le{6.0}$ using the Lyman-$\displaystyleα$ (Ly$\displaystyleα$) forest flux auto-correlation function. Our auto-differentiable JAX-based framework accelerates the surrogate model generation process using approximately 100 sparsely sampled Nyx hydrodynamical simulat… ▽ More We present a neural network emulator to constrain the thermal parameters of the intergalactic medium (IGM) at $\displaystyle{5.4}\le{z}\le{6.0}$ using the Lyman-$\displaystyleα$ (Ly$\displaystyleα$) forest flux auto-correlation function. Our auto-differentiable JAX-based framework accelerates the surrogate model generation process using approximately 100 sparsely sampled Nyx hydrodynamical simulations with varying combinations of thermal parameters, i.e., the temperature at mean density $\displaystyle{T}_{0}$, the slope of the temperature$\displaystyle-$density relation $\displaystyleγ$, and the mean transmission flux $\displaystyle{\langle}{F}{\rangle}$. We show that this emulator has a typical accuracy of 1.0% across the specified redshift range. Bayesian inference of the IGM thermal parameters, incorporating emulator uncertainty propagation, is further expedited using NumPyro Hamiltonian Monte Carlo. We compare both the inference results and computational cost of our framework with the traditional nearest-neighbor interpolation approach applied to the same set of mock Ly$α$ flux. By examining the credibility contours of the marginalized posteriors for $\displaystyle{T}_{0},γ,\text{and}{\langle}{F}{\rangle}$ obtained using the emulator, the statistical reliability of measurements is established through inference on 100 realistic mock data sets of the auto-correlation function. △ Less

Submitted 20 December, 2024; v1 submitted 8 October, 2024; originally announced October 2024.

A Quasar-Anchored Protocluster at z=6.6 in the ASPIRE Survey: II. An Environmental Analysis of Galaxy Properties in an Overdense Structure

Authors: Jaclyn B. Champagne, Feige Wang, Jinyi Yang, Xiaohui Fan, Joseph F. Hennawi, Fengwu Sun, Eduardo Bañados, Sarah E. I. Bosman, Tiago Costa, Melanie Habouzit, Xiangyu Jin, Hyunsung D. Jun, Mingyu Li, Weizhe Liu, Federica Loiacono, Alessandro Lupi, Chiara Mazzucchelli, Maria Pudoka, Sofia Rojas-Ruiz, Wei Leong Tee, Maxime Trebitsch, Haowen Zhang, Ming-Yang Zhuang, Siwei Zou

Abstract: We present paper II comprising a 35 arcmin$^2$ \textit{JWST}/NIRCam imaging and wide-field slitless spectroscopy mosaic centered on J0305$-$3150, a luminous quasar at $z=6.61$. The F356W grism data reveals 124 [OIII]+H$β$ emitters at $5.3<z<7$, 53 of which constitute a protocluster spanning (10 cMpc)$^2$ across $6.5<z<6.8$. We find no evidence of any broad-line AGN in individual galaxies or stacki… ▽ More We present paper II comprising a 35 arcmin$^2$ \textit{JWST}/NIRCam imaging and wide-field slitless spectroscopy mosaic centered on J0305$-$3150, a luminous quasar at $z=6.61$. The F356W grism data reveals 124 [OIII]+H$β$ emitters at $5.3<z<7$, 53 of which constitute a protocluster spanning (10 cMpc)$^2$ across $6.5<z<6.8$. We find no evidence of any broad-line AGN in individual galaxies or stacking, reporting a median H$β$ FWHM of 585 $\pm$ 152 km s$^{-1}$; however, the mass-excitation diagram and ``little red dot" color and compactness criteria suggest that there are a few AGN candidates on the outskirts of the protocluster. We fit the spectral energy distributions (SEDs) of the [OIII] emitters with \texttt{Prospector} and \texttt{Bagpipes}, and find that none of the SED-derived properties (stellar mass, age, or star formation rate) correlates with proximity to the quasar. While there is no correlation between galaxy age and local galaxy density, we find modest correlations between local galaxy density with increasing stellar mass, decreasing 10-to-100 Myr star formation rate ratios and decreasing nebular line equivalent widths. We further find that the protocluster galaxies are consistent with being more massive, older, and hosting higher star formation rates than the field sample at the 3$σ$ level, distributed in a filamentary structure which supports inside-out formation of the protocluster. There is modest evidence that galaxy evolution proceeds differently as a function of the density of local environment within protoclusters during the epoch of reionization, and the central quasar has little effect on the galaxy properties of the surrounding structure. △ Less

Submitted 20 February, 2025; v1 submitted 4 October, 2024; originally announced October 2024.

Comments: accepted to ApJ. 24 pages, 13 figures. Paper I: arXiv:2410.03826

A Quasar-Anchored Protocluster at z=6.6 in the ASPIRE Survey: I. Properties of [OIII] Emitters in a 10 Mpc Overdensity Structure

Authors: Jaclyn B. Champagne, Feige Wang, Haowen Zhang, Jinyi Yang, Xiaohui Fan, Joseph F. Hennawi, Fengwu Sun, Eduardo Bañados, Sarah E. I. Bosman, Tiago Costa, Anna-Christina Eilers, Ryan Endsley, Xiangyu Jin, Hyunsung D. Jun, Mingyu Li, Xiaojing Lin, Weizhe Liu, Federica Loiacono, Alessandro Lupi, Chiara Mazzucchelli, Maria Pudoka, Klaudia Protušovà, Sofia Rojas-Ruiz, Wei Leong Tee, Maxime Trebitsch , et al. (3 additional authors not shown)

Abstract: ASPIRE (A SPectroscopic survey of bIased halos in the Reionization Era) is a quasar legacy survey primarily using \textit{JWST} to target a sample of 25 $z>6$ quasars with NIRCam slitless spectroscopy and imaging. The first study in this series found evidence of a strong overdensity of galaxies around J0305$-$3150, a luminous quasar at $z=6.61$, within a single NIRCam pointing obtained in JWST Cyc… ▽ More ASPIRE (A SPectroscopic survey of bIased halos in the Reionization Era) is a quasar legacy survey primarily using \textit{JWST} to target a sample of 25 $z>6$ quasars with NIRCam slitless spectroscopy and imaging. The first study in this series found evidence of a strong overdensity of galaxies around J0305$-$3150, a luminous quasar at $z=6.61$, within a single NIRCam pointing obtained in JWST Cycle 1. Here, we present the first results of a JWST Cycle 2 mosaic that covers 35 arcmin$^2$ with NIRCam imaging/WFSS of the same field to investigate the spatial extent of the putative protocluster. The F356W grism data targets [OIII]+H$β$ at $5.3<z<7$ and reveals a population of 124 line emitters down to a flux limit of 1.2$\times$10$^{-18}$ erg s$^{-1}$ cm$^{-2}$. Fifty-three of these galaxies lie at $6.5<z<6.8$ spanning 10 cMpc on the sky, corresponding to an overdensity within a 2500 cMpc$^3$ volume of 12.5 $\pm$ 2.6, anchored by the quasar. Comparing to the [OIII] luminosity function from the Emission line galaxies and Intergalactic Gas in the Epoch of Reionization (EIGER) project, we find a dearth of faint [OIII] emitters at log(L/erg\,s$^{-1}$) $<$ 42.3, which we suggest is consistent with either bursty star formation causing galaxies to scatter around the grism detection limit or modest suppression from quasar feedback. While we find a strong filamentary overdensity of [OIII] emitters consistent with a protocluster, we suggest that we could be insensitive to a population of older, more massive Lyman-break galaxies with weak nebular emission on scales $>10$ cMpc. △ Less

Submitted 20 February, 2025; v1 submitted 4 October, 2024; originally announced October 2024.

Comments: accepted for publication in ApJ. 18 pages, 7 figures. see paper II: arXiv:2410.03827

A SPectroscopic survey of biased halos In the Reionization Era (ASPIRE): JWST Supports Earlier Reionization around [OIII] Emitters

Authors: Xiangyu Jin, Jinyi Yang, Xiaohui Fan, Feige Wang, Koki Kakiichi, Romain A. Meyer, George D. Becker, Siwei Zou, Eduardo Bañados, Jaclyn B. Champagne, Valentina D'Odorico, Minghao Yue, Sarah E. I. Bosman, Zheng Cai, Anna-Christina Eilers, Joseph F. Hennawi, Hyunsung D. Jun, Mingyu Li, Zihao Li, Weizhe Liu, Maria Pudoka, Sindhu Satyavolu, Fengwu Sun, Wei Leong Tee, Yunjing Wu

Abstract: Understanding when and how reionization happened is crucial for studying the early structure formation and the properties of first galaxies in the Universe. At $z>5.5$, the observed IGM optical depth shows a significant scatter, indicating an inhomogeneous reionization process. However, the nature of the inhomogeneous reionization remains debated. ASPIRE is a JWST Cycle 1 program that has spectros… ▽ More Understanding when and how reionization happened is crucial for studying the early structure formation and the properties of first galaxies in the Universe. At $z>5.5$, the observed IGM optical depth shows a significant scatter, indicating an inhomogeneous reionization process. However, the nature of the inhomogeneous reionization remains debated. ASPIRE is a JWST Cycle 1 program that has spectroscopically identified $>400$ [OIII] emitters in 25 quasar fields at $z>6.5$. Combined with deep ground-based optical spectroscopy of ASPIRE quasars, ASPIRE program provides the current largest sample for IGM-galaxy connection studies during cosmic reionization. We present the first results of IGM effective optical depth measurements around [OIII] emitters using 14 ASPIRE quasar fields. We find the IGM transmission is tightly related with reionization-era galaxies to the extent that significant excess of Ly$α$ transmission exists around [OIII] emitters. We measure the stacked IGM effective optical depth of IGM patches associated with [OIII] emitters and find they reach the same IGM effective optical depth at least dz~0.1 ahead of those IGM patches where no [OIII] emitters are detected, supporting earlier reionization around [OIII] emitters. Our results indicate an enhancement in IGM Ly$α$ transmission around [OIII] emitters at scales beyond 25 $h^{-1}$ cMpc, consistent with the predicted topology of reionization from fluctuating UV background (UVB) models. △ Less

Submitted 2 October, 2024; originally announced October 2024.

Comments: 27 pages, 19 figures, 1 table, accepted for publication in ApJ

"Little red dots" cannot reside in the same dark matter halos as comparably luminous unobscured quasars

Authors: Elia Pizzati, Joseph F. Hennawi, Joop Schaye, Anna-Christina Eilers, Jiamu Huang, Jan-Torge Schindler, Feige Wang

Abstract: The James Webb Space Telescope (JWST) has uncovered a new population of candidate broad-line AGN emerging in the early Universe, named ''little red dots'' (LRDs) because of their compactness and red colors at optical wavelengths. LRDs appear to be surprisingly abundant ($\approx 10^{-5}\,\mathrm{cMpc}^{-3}$) given that their inferred bolometric luminosities largely overlap with the ones of the UV-… ▽ More The James Webb Space Telescope (JWST) has uncovered a new population of candidate broad-line AGN emerging in the early Universe, named ''little red dots'' (LRDs) because of their compactness and red colors at optical wavelengths. LRDs appear to be surprisingly abundant ($\approx 10^{-5}\,\mathrm{cMpc}^{-3}$) given that their inferred bolometric luminosities largely overlap with the ones of the UV-luminous quasars identified at high $z$ in wide-field spectroscopic surveys. In this work, we investigate how the population of LRDs and/or other UV-obscured AGN relates to the one of unobscured, UV-selected quasars. By comparing their number densities, we infer an extremely large and rapidly evolving obscured:unobscured ratio, ranging from $\approx20:1$ at $z\approx4$ to $\approx2300:1$ at $z\approx7$, and possibly extending out to very high ($\approx10^{47}\,\mathrm{erg}\,\mathrm{s}^{-1}$) bolometric luminosities. This large obscured:unobscured ratio is incompatible with the UV-luminous duty cycle measured for unobscured quasars at $z\approx4-6$, suggesting that LRDs are too abundant to be hosted by the same halos as unobscured quasars. This implies that either (a) the bolometric luminosities of LRDs are strongly overestimated (possibly because LRDs are dominated by stellar galaxy light) or (b) LRDs follow different scaling relations than the ones of UV-selected quasars, representing a new population of accreting SMBHs emerging in the early Universe. A direct comparison between the clustering of LRDs and the one of faint UV-selected quasars will ultimately confirm these findings, and shed light on key properties of LRDs such as their host mass distribution and duty cycle. We provide a mock analysis for the clustering of LRDs and show that it is feasible with current and upcoming JWST surveys. △ Less

Submitted 2 May, 2025; v1 submitted 26 September, 2024; originally announced September 2024.

Comments: 15 pages. Accepted for publication in MNRAS

A SPectroscopic survey of biased halos In the Reionization Era (ASPIRE): Broad-line AGN at $z=4-5$ revealed by JWST/NIRCam WFSS

Authors: Xiaojing Lin, Feige Wang, Xiaohui Fan, Zheng Cai, Jaclyn B. Champagne, Fengwu Sun, Marta Volonteri, Jinyi Yang, Joseph F. Hennawi, Eduardo Bañados, Aaron Barth, Anna-Christina Eilers, Emanuele Paolo Farina, Weizhe Liu, Xiangyu Jin, Hyunsung D. Jun, Alessandro Lupi, Koki Kakiichi, Chiara Mazzucchelli, Masafusa Onoue, Zhiwei Pan, Elia Pizzati, Sofía Rojas-Ruiz, Jan-Torge Schindler, Benny Trakhtenbrot , et al. (11 additional authors not shown)

Abstract: Low-luminosity AGNs with low-mass black holes (BHs) in the early universe are fundamental to understanding the BH growth and their co-evolution with the host galaxies. Utilizing JWST NIRCam Wide Field Slitless Spectroscopy (WFSS), we perform a systematic search for broad-line ${\rm Hα}$ emitters (BHAEs) at $z\approx 4-5$ in 25 fields of the ASPIRE (A SPectroscopic survey of biased halos In the Rei… ▽ More Low-luminosity AGNs with low-mass black holes (BHs) in the early universe are fundamental to understanding the BH growth and their co-evolution with the host galaxies. Utilizing JWST NIRCam Wide Field Slitless Spectroscopy (WFSS), we perform a systematic search for broad-line ${\rm Hα}$ emitters (BHAEs) at $z\approx 4-5$ in 25 fields of the ASPIRE (A SPectroscopic survey of biased halos In the Reionization Era) project, covering a total area of 275 arcmin$^2$. We identify 16 BHAEs with FWHM of the broad components spanning from $\sim$ 1000 km s$^{-1}$ to 3000 km s$^{-1}$. Assuming the broad linewidths arise due to Doppler broadening around BHs, the implied BH masses range from $10^7$ to $10^{8}~M_\odot$, with broad ${\rm Hα}$-converted bolometric luminosity of $10^{44.5}-10^{45.5}$ erg s$^{-1}$ and Eddington ratios of $0.07-0.47$. The spatially extended structure of the F200W stacked image may trace the stellar light from the host galaxies. The ${\rm Hα}$ luminosity function indicates an increasing AGN fraction towards the higher ${\rm Hα}$ luminosities. We find possible evidence for clustering of BHAEs: two sources are at the same redshift with a projected separation of 519 kpc; one BHAE appears as a composite system residing in an overdense region with three close companion ${\rm Hα}$ emitters. Three BHAEs exhibit blueshifted absorption troughs indicative of the presence of high-column-density gas. We find the broad-line and photometrically selected BHAE samples exhibit different distributions in the optical continuum slopes, which can be attributed to their different selection methods. The ASPIRE broad-line ${\rm Hα}$ sample provides a good database for future studies of faint AGN populations at high redshift. △ Less

Submitted 24 July, 2024; originally announced July 2024.

Comments: 19 pages, 13 figures, 4 tables. Accepted by the ApJ

A blazar in the epoch of reionization

Authors: Eduardo Banados, Emmanuel Momjian, Thomas Connor, Silvia Belladitta, Roberto Decarli, Chiara Mazzucchelli, Bram P. Venemans, Fabian Walter, Feige Wang, Zhang-Liang Xie, Aaron J. Barth, Anna-Christina Eilers, Xiaohui Fan, Yana Khusanova, Jan-Torge Schindler, Daniel Stern, Jinyi Yang, Irham Taufik Andika, Chris Carilli, Emanuele P. Farina, Andrew Fabian, Joseph F. Hennawi, Antonio Pensabene, Sofia Rojas-Ruiz

Abstract: Relativistic jets are thought to play a crucial role in the formation and evolution of massive galaxies and supermassive black holes. Blazars, which are quasars with jets aligned along our line of sight, provide insights into the jetted population and have been observed up to redshifts of z=6.1. Here, we report the discovery and multi-wavelength characterization of the blazar VLASS J041009.05-0139… ▽ More Relativistic jets are thought to play a crucial role in the formation and evolution of massive galaxies and supermassive black holes. Blazars, which are quasars with jets aligned along our line of sight, provide insights into the jetted population and have been observed up to redshifts of z=6.1. Here, we report the discovery and multi-wavelength characterization of the blazar VLASS J041009.05-013919.88 at z=7 (Universe's age ~750 Myr), powered by a ~7x10^8 Msun black hole. The presence of this high-redshift blazar implies a large population of similar but unaligned jetted sources in the early Universe. Our findings suggest two possible scenarios: in one, the jet in J0410-0139 is intrinsically low-power but appears highly luminous due to relativistic beaming, suggesting that most UV-bright quasars at this redshift host jets. Alternatively, if J0410-0139 represents an intrinsically powerful radio source, there should be hundreds to thousands of radio-quiet quasars at z~7 with properties similar to J0410-0139, a prediction in tension with observed quasar densities based on their UV luminosity function. These results support the hypothesis that rapid black hole growth in the early Universe may be driven by jet-enhanced or obscured super-Eddington accretion, potentially playing a key role in forming massive black holes during the epoch of reionization. △ Less

Submitted 16 December, 2024; v1 submitted 9 July, 2024; originally announced July 2024.

Comments: Updated to match accepted/final version in Nature Astronomy

Quantifying the Precision of IGM Damping Wing Measurements Towards Quasars

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

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

Submitted 17 June, 2024; originally announced June 2024.

Comments: Submitted to MNRAS

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

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

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

Submitted 15 April, 2025; v1 submitted 17 June, 2024; originally announced June 2024.

Comments: Accepted by MNRAS

[OIII] emission in z=2 quasars with and without Broad Absorption Lines

Authors: Matthew J. Temple, Amy L. Rankine, Manda Banerji, Joseph F. Hennawi, Paul C. Hewett, James H. Matthews, Riccardo Nanni, Claudio Ricci, Gordon T. Richards

Abstract: Understanding the links between different phases of outflows from active galactic nuclei is a key goal in extragalactic astrophysics. Here we compare [OIII] $λλ$4960,5008 outflow signatures in quasars with and without Broad Absorption Lines (BALs), aiming to test how the broad absorption troughs seen in the rest-frame ultraviolet are linked to the narrow line region outflows seen in the rest-frame… ▽ More Understanding the links between different phases of outflows from active galactic nuclei is a key goal in extragalactic astrophysics. Here we compare [OIII] $λλ$4960,5008 outflow signatures in quasars with and without Broad Absorption Lines (BALs), aiming to test how the broad absorption troughs seen in the rest-frame ultraviolet are linked to the narrow line region outflows seen in the rest-frame optical. We present new near-infrared spectra from Magellan/FIRE which cover [OIII] in 12 quasars with 2.1 < z < 2.3, selected to have strong outflow signatures in CIV $λ$1550. Combining with data from the literature, we build a sample of 73 BAL, 115 miniBAL and 125 non-BAL QSOs with 1.5 < z < 2.6. The strength and velocity width of [OIII] correlate strongly with the CIV emission properties, but no significant difference is seen in the [OIII] emission-line properties between the BALs, non-BALs and miniBALs once the dependence on CIV emission is taken into account. A weak correlation is observed between the velocities of CIV BALs and [OIII] emission, which is accounted for by the fact that both outflow signatures correlate with the underlying CIV emission properties. Our results add to the growing evidence that BALs and non-BALs are drawn from the same parent population and are consistent with a scenario wherein BAL troughs are intermittent tracers of persistent quasar outflows, with a part of such outflow becoming optically thick along our line-of-sight for sporadic periods of time within which BALs are observed. △ Less

Submitted 17 June, 2024; originally announced June 2024.

Comments: 14 pages, 6 figures. Accepted for publication in MNRAS

Journal ref: MNRAS 532, 424-437 (2024)

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 of the mean signal-to-noise ratios per $110$ km s$^{-1}$ pixel in the 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. Despite the diversity in instrumental setups and spectral quality, the data set is uniformly processed and well-characterized, making it ideally suited for several scientific goals, including 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 29 May, 2025; v1 submitted 11 June, 2024; originally announced June 2024.

Comments: 21 pages, 11 figures, published in MNRAS

MAMMOTH-Subaru. II. Diverse Populations of Circumgalactic Ly$α$ Nebulae at Cosmic Noon

Authors: Mingyu Li, Haibin Zhang, Zheng Cai, Yongming Liang, Nobunari Kashikawa, Ke Ma, Xiaohui Fan, J. Xavier Prochaska, Bjorn H. C. Emonts, Xin Wang, Yunjing Wu, Shiwu Zhang, Qiong Li, Sean D. Johnson, Minghao Yue, Fabrizio Arrigoni Battaia, Sebastiano Cantalupo, Joseph F. Hennawi, Satoshi Kikuta, Yuanhang Ning, Masami Ouchi, Rhythm Shimakawa, Ben Wang, Weichen Wang, Zheng Zheng , et al. (1 additional authors not shown)

Abstract: Circumgalactic Lyman-alpha (Ly$α$) nebulae are gaseous halos around galaxies exhibiting luminous extended Ly$α$ emission. This work investigates Ly$α$ nebulae from deep imaging of $\sim12~\mathrm{deg}^2$ sky, targeted by the MAMMOTH-Subaru survey. Utilizing the wide-field capability of Hyper Suprime-Cam (HSC), we present one of the largest blind Ly$α$ nebula selections, including QSO nebulae, Ly… ▽ More Circumgalactic Lyman-alpha (Ly$α$) nebulae are gaseous halos around galaxies exhibiting luminous extended Ly$α$ emission. This work investigates Ly$α$ nebulae from deep imaging of $\sim12~\mathrm{deg}^2$ sky, targeted by the MAMMOTH-Subaru survey. Utilizing the wide-field capability of Hyper Suprime-Cam (HSC), we present one of the largest blind Ly$α$ nebula selections, including QSO nebulae, Ly$α$ blobs, and radio galaxy nebulae down to typical $2σ$ Ly$α$ surface brightness of $(5-10)\times10^{-18}\mathrm{~erg~s^{-1}~cm^{-2}~arcsec^{-2}}$. The sample contains 117 nebulae with Ly$α$ sizes of 40 - 400 kpc, and the most gigantic one spans about 365 kpc, referred to as the Ivory Nebula. Combining multiwavelength data, we investigate diverse nebula populations and associated galaxies. We find a small fraction of Ly$α$ nebulae have QSOs ($\sim7\%$), luminous infrared galaxies ($\sim1\%$), and radio galaxies ($\sim 2\%$). Remarkably, among the 28 enormous Ly$α$ nebulae (ELANe) exceeding 100 kpc, about 80\% are associated with UV-faint galaxies ($M_\mathrm{UV} > -22$), categorized as Type II ELANe. We underscore that Type II ELANe constitute the majority but remain largely hidden in current galaxy and QSO surveys. Dusty starburst and obscured AGN activity are proposed to explain the nature of Type II ELANe. The SED of stacking all Ly$α$ nebulae also reveals signs of massive dusty star-forming galaxies with obscured AGNs. We propose a model to explain the dusty nature where the diverse populations of Ly$α$ nebulae capture massive galaxies at different evolutionary stages undergoing violent assembling. Ly$α$ nebulae provide critical insights into the formation and evolution of today's massive cluster galaxies at cosmic noon. △ Less

Submitted 26 September, 2024; v1 submitted 21 May, 2024; originally announced May 2024.

Comments: Accepted by ApJS after minor revision; 26 pages, 10 figures, 3 tables,

AGNfitter-rx: Modelling the radio-to-X-ray SEDs of AGNs

Authors: L. N. Martínez-Ramírez, G. Calistro Rivera, Elisabeta Lusso, F. E. Bauer, Emanuele Nardini, Johannes Buchner, Michael J. I. Brown, Juan C. B. Pineda, Matthew J. Temple, Manda Banerji, M. Stalevski, Joseph F. Hennawi

Abstract: We present new frontiers in the modelling of the spectral energy distributions (SED) of active galaxies by introducing the radio-to-X-ray fitting capabilities of the publicly available Bayesian code AGNfitter. The new code release, called AGNfitter-rx, models the broad-band photometry covering the radio, infrared (IR), optical, ultraviolet (UV) and X-ray bands consistently, using a combination of… ▽ More We present new frontiers in the modelling of the spectral energy distributions (SED) of active galaxies by introducing the radio-to-X-ray fitting capabilities of the publicly available Bayesian code AGNfitter. The new code release, called AGNfitter-rx, models the broad-band photometry covering the radio, infrared (IR), optical, ultraviolet (UV) and X-ray bands consistently, using a combination of theoretical and semi-empirical models of the AGN and host galaxy emission. This framework enables the detailed characterization of four physical components of the active nuclei: the accretion disk, the hot dusty torus, the relativistic jets/core radio emission, and the hot corona; alongside modeling three components within the host galaxy: stellar populations, cold dust, and the radio emission from the star-forming regions. Applying AGNfitter-rx to a diverse sample of 36 AGN SEDs at z<0.7 from the AGN SED ATLAS, we investigate and compare the performance of state-of-the-art torus and accretion disk emission models on fit quality and inferred physical parameters. We find that clumpy torus models that include polar winds and semi-empirical accretion disk templates including emission line features significantly increase the fit quality in 67% of the sources, by effectively reducing by $2σ$ fit residuals in the $1.5-5 μ\rm m$ and $0.7 μ\rm m$ regimes.We demonstrate that, by applying AGNfitter-rx on photometric data, we are able to estimate inclination and opening angles of the torus, consistent with spectroscopic classifications within the AGN unified model, as well as black hole mass estimates in agreement with virial estimates based on H$α$. The wavelength coverage and the flexibility for the inclusion of state-of-the-art theoretical models make AGNfitter-rx a unique tool for the further development of SED modelling for AGNs in present and future radio-to-X-ray galaxy surveys. △ Less

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

Comments: 20 pages, 10 figures, Accepted for publication by A&A

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

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

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

Submitted 23 April, 2024; originally announced April 2024.

Comments: accepted for publication in ApJ

A unified model for the clustering of quasars and galaxies at $z\approx6$

Authors: Elia Pizzati, Joseph F. Hennawi, Joop Schaye, Matthieu Schaller, Anna-Christina Eilers, Feige Wang, Carlos S. Frenk, Willem Elbers, John C. Helly, Ruari Mackenzie, Jorryt Matthee, Rongmon Bordoloi, Daichi Kashino, Rohan P. Naidu, Minghao Yue

Abstract: Recent observations from the EIGER JWST program have measured for the first time the quasar-galaxy cross-correlation function at $z\approx6$. The auto-correlation function of faint $z\approx6$ quasars was also recently estimated. These measurements provide key insights into the properties of quasars and galaxies at high redshift and their relation with the host dark matter halos. In this work, we… ▽ More Recent observations from the EIGER JWST program have measured for the first time the quasar-galaxy cross-correlation function at $z\approx6$. The auto-correlation function of faint $z\approx6$ quasars was also recently estimated. These measurements provide key insights into the properties of quasars and galaxies at high redshift and their relation with the host dark matter halos. In this work, we interpret these data building upon an empirical quasar population model that has been applied successfully to quasar clustering and demographic measurements at $z\approx2-4$. We make use of a new, large-volume N-body simulation with more than a trillion particles, FLAMINGO-10k, to model quasars and galaxies simultaneously. We successfully reproduce observations of $z\approx6$ quasars and galaxies (i.e., their clustering properties and luminosity functions), and infer key quantities such as their luminosity-halo mass relation, the mass function of their host halos, and their duty cycle/occupation fraction. Our key findings are: (i) quasars reside on average in $\approx10^{12.5}\,{\rm M}_\odot$ halos (corresponding to $\approx5σ$ fluctuations in the initial conditions of the linear density field), but the distribution of host halo masses is quite broad; (ii) the duty cycle of (UV-bright) quasar activity is relatively low ($\approx1\%$); (iii) galaxies (that are bright in [OIII]) live in much smaller halos ($\approx10^{10.9}\,{\rm M}_\odot$) and have a larger duty cycle (occupation fraction) of $\approx13\%$. Finally, we focus on the inferred properties of quasars and present a homogeneous analysis of their evolution with redshift. The picture that emerges reveals a strong evolution of the host halo mass and duty cycle of quasars at $z\approx2-6$, and calls for new investigations of the role of quasar activity across cosmic time. △ Less

Submitted 5 October, 2024; v1 submitted 18 March, 2024; originally announced March 2024.

Comments: Accepted for publication in MNRAS; 15 pages, 4 figures

EIGER VI. The Correlation Function, Host Halo Mass and Duty Cycle of Luminous Quasars at $z\gtrsim6$

Authors: Anna-Christina Eilers, Ruari Mackenzie, Elia Pizzati, Jorryt Matthee, Joseph F. Hennawi, Haowen Zhang, Rongmon Bordoloi, Daichi Kashino, Simon J. Lilly, Rohan P. Naidu, Robert A. Simcoe, Minghao Yue, Carlos S. Frenk, John C. Helly, Matthieu Schaller, Joop Schaye

Abstract: We expect luminous ($M_{1450}\lesssim-26.5$) high-redshift quasars to trace the highest density peaks in the early universe. Here, we present observations of four $z\gtrsim6$ quasar fields using JWST/NIRCam in imaging and widefield slitless spectroscopy mode and report a wide range in the number of detected [OIII]-emitting galaxies in the quasars' environments, ranging between a density enhancemen… ▽ More We expect luminous ($M_{1450}\lesssim-26.5$) high-redshift quasars to trace the highest density peaks in the early universe. Here, we present observations of four $z\gtrsim6$ quasar fields using JWST/NIRCam in imaging and widefield slitless spectroscopy mode and report a wide range in the number of detected [OIII]-emitting galaxies in the quasars' environments, ranging between a density enhancement of $δ\approx65$ within a $2$ cMpc radius - one of the largest proto-clusters during the Epoch of Reionization discovered to date - to a density contrast consistent with zero, indicating the presence of a UV-luminous quasar in a region comparable to the average density of the universe. By measuring the two-point cross-correlation function of quasars and their surrounding galaxies, as well as the galaxy auto-correlation function, we infer a correlation length of quasars at $\langle z\rangle=6.25$ of $r_0^{\rm QQ}=22.0^{+3.0}_{-2.9}~{\rm cMpc}\,h^{-1}$, while we obtain a correlation length of the [OIII]-emitting galaxies of $r_0^{\rm GG}=4.1\pm0.3~{\rm cMpc}\,h^{-1}$. By comparing the correlation functions to dark-matter-only simulations we estimate the minimum mass of the quasars' host dark matter halos to be $\log_{10}(M_{\rm halo, min}/M_\odot)=12.43^{+0.13}_{-0.15}$ (and $\log_{10}(M_{\rm halo, min}^{\rm [OIII]}/M_\odot) = 10.56^{+0.05}_{-0.03}$ for the [OIII]-emitters), indicating that (a) luminous quasars do not necessarily reside within the most overdense regions in the early universe, and that (b) the UV-luminous duty cycle of quasar activity at these redshifts is $f_{\rm duty}\ll1$. Such short quasar activity timescales challenge our understanding of early supermassive black hole growth and provide evidence for highly dust-obscured growth phases or episodic, radiatively inefficient accretion rates. △ Less

Submitted 4 September, 2024; v1 submitted 12 March, 2024; originally announced March 2024.

Comments: ApJ accepted

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

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

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

Submitted 2 February, 2024; originally announced February 2024.

Comments: Accepted for publication in ApJL

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

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

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

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

Submitted 13 December, 2023; originally announced December 2023.

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

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

Authors: Frederick B. Davies, Joseph F. Hennawi

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

Submitted 11 December, 2023; originally announced December 2023.

Comments: 18 pages, 17 figures, submitted to ApJ