Showing 1–50 of 91 results for author: Eilers, A

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

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

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

Submitted 11 June, 2024; originally announced June 2024.

Comments: 18 pages, 10 figures, submitted to MNRAS

A quasar-galaxy merger at $z\sim 6.2$: rapid host growth via accretion of two massive satellite galaxies

Authors: Roberto Decarli, Federica Loiacono, Emanuele Paolo Farina, Massimo Dotti, Alessandro Lupi, Romain A. Meyer, Marco Mignoli, Antonio Pensabene, Michael A. Strauss, Bram Venemans, Jinyi Yang, Fabian Walter, Julien Wolf, Eduardo Bañados, Laura Blecha, Sarah Bosman, Chris L. Carilli, Andrea Comastri, Thomas Connor, Tiago Costa, Anna-Christina Eilers, Xiaohui Fan, Roberto Gilli, Hyunsung D. Jun, Weizhe Liu , et al. (16 additional authors not shown)

Abstract: We present JWST/NIRSpec Integral Field Spectroscopy in the rest-frame optical bands of the system PJ308-21, a quasar at $z=6.2342$ caught as its host galaxy interacts with companion galaxies. We detect spatially extended emission of several emission lines (H$α$, H$β$, [OIII], [NII], [SII], HeII), which we use to study the properties of the ionized phase of the interstellar medium: the source and h… ▽ More We present JWST/NIRSpec Integral Field Spectroscopy in the rest-frame optical bands of the system PJ308-21, a quasar at $z=6.2342$ caught as its host galaxy interacts with companion galaxies. We detect spatially extended emission of several emission lines (H$α$, H$β$, [OIII], [NII], [SII], HeII), which we use to study the properties of the ionized phase of the interstellar medium: the source and hardness of the photoionizing radiation field, metallicity, dust reddening, electron density and temperature, and star formation. We also marginally detect continuum starlight emission associated with the companion sources. We find that at least two independent satellite galaxies are part of the system. While the quasar host appears highly enriched and obscured, with AGN-like photoionization conditions, the western companion shows minimal dust extinction, low metallicity ($Z\sim0.4$ Z$_\odot$), and star-formation driven photoionization. The eastern companion shows higher extinction and metallicity ($Z\sim0.8$ Z$_\odot$) compared to the western companion, and it is at least partially photoionized by the nearby quasar. We do not find any indication of AGN in the companion sources. Our study shows that while the quasar host galaxy is already very massive ($M_{\rm dyn}>10^{11}$ M$_\odot$), it is still rapidly building up by accreting two relatively massive ($M_{\rm star}\sim 10^{10}$ M$_\odot$) companion sources. This dataset showcases the power of JWST in exposing the build-up of massive galaxies in the first Gyr of the Universe. △ Less

Submitted 10 June, 2024; originally announced June 2024.

Comments: 15 pages, 16 figures. Accepted for publication in A&A

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

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

Abstract: Recent quasar absorption line observations suggest that reionization may end as late as $z \approx 5.3$. As a means to search for large neutral hydrogen islands at $z<6$, we revisit long dark gaps in the Ly$β$ forest in VLT/X-Shooter and Keck/ESI quasar spectra. We stack the Ly$α$ forest corresponding to both edges of these Ly$β$ dark gaps and identify a damping wing-like extended absorption profi… ▽ More Recent quasar absorption line observations suggest that reionization may end as late as $z \approx 5.3$. As a means to search for large neutral hydrogen islands at $z<6$, we revisit long dark gaps in the Ly$β$ forest in VLT/X-Shooter and Keck/ESI quasar spectra. We stack the Ly$α$ forest corresponding to both edges of these Ly$β$ dark gaps and identify a damping wing-like extended absorption profile. The average redshift of the stacked forest is $z=5.8$. By comparing these observations with reionization simulations, we infer that such a damping wing-like feature can be naturally explained if these gaps are at least partially created by neutral islands. Conversely, simulated dark gaps lacking neutral hydrogen struggle to replicate the observed damping wing features. Furthermore, this damping wing-like profile implies that the volume-averaged neutral hydrogen fraction must be $\langle x_{\rm HI} \rangle \geq 6.1 \pm 3.9\%$ at $z = 5.8$. Our results offer robust evidence that reionization extends below $z=6$. △ Less

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

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

3D structure of the Milky Way out to 10 kpc from the Sun. Catalogue of large molecular clouds in the Galactic Plane

Authors: Sara Rezaei Kh., Henrik Beuther, Robert A. Benjamin, Anna-Christina Eilers, Thomas Henning, Maria J. Jiménez-Donaire, Marc-Antoine Miville-Deschênes

Abstract: Understanding the 3D structure of the Milky Way is a crucial step in deriving properties of the star-forming regions, as well as the Galaxy as a whole. We present a novel 3D map of the Milky Way plane that extends to 10 kpc distance from the Sun. We leverage the wealth of information in the near-IR APOGEE dataset and combine that with our state-of-the-art 3D mapping technique using Bayesian statis… ▽ More Understanding the 3D structure of the Milky Way is a crucial step in deriving properties of the star-forming regions, as well as the Galaxy as a whole. We present a novel 3D map of the Milky Way plane that extends to 10 kpc distance from the Sun. We leverage the wealth of information in the near-IR APOGEE dataset and combine that with our state-of-the-art 3D mapping technique using Bayesian statistics and the Gaussian process to provide a large-scale 3D map of the dust in the Milky Way. Our map stretches across 10 kpc along both the X and Y axes, and 750 pc in the Z direction, perpendicular to the Galactic plane. Our results reveal multi-scale over-densities as well as large cavities in the Galactic plane and shed new light on the Galactic structure and spiral arms. We also provide a catalogue of large molecular clouds identified by our map with accurate distance and volume density estimates. Utilising volume densities derived from this map, we explore mass distribution across various Galactocentric radii. A general decline towards the outer Galaxy is observed, followed by local peaks, some aligning with established features like the Molecular Ring and segments of the spiral arms. Moreover, this work explores extragalactic observational effects on derived properties of molecular clouds by demonstrating the potential biases arising from column density measurements in inferring properties of these regions, and opens exciting avenues for further exploration and analysis, offering a deeper perspective on the complex processes that shape our galaxy and beyond. △ Less

Submitted 15 May, 2024; originally announced May 2024.

Comments: Revised version after first referee report; to appear in A&A. Includes SRKh's goodbye letter to Astronomy!

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

Stacking X-ray Observations of "Little Red Dots": Implications for their AGN Properties

Authors: Minghao Yue, Anna-Christina Eilers, Tonima Tasnim Ananna, Christos Panagiotou, Erin Kara, Takamitsu Miyaji

Abstract: Recent {\em James Webb Space Telescope (JWST)} observations have revealed a population of compact extragalactic objects at $z\gtrsim4$ with red near-infrared colors, which have been dubbed as ``Little Red Dots" (LRDs). The spectroscopically-selected LRDs exhibit broad H$α$ emission lines, which has led to the conclusion that type-I active galactic nuclei (AGN) are harbored in the galaxies' dust-re… ▽ More Recent {\em James Webb Space Telescope (JWST)} observations have revealed a population of compact extragalactic objects at $z\gtrsim4$ with red near-infrared colors, which have been dubbed as ``Little Red Dots" (LRDs). The spectroscopically-selected LRDs exhibit broad H$α$ emission lines, which has led to the conclusion that type-I active galactic nuclei (AGN) are harbored in the galaxies' dust-reddened cores. However, other mechanisms, like strong outflowing winds, could also produce broad H$α$ emission lines, and thus, the nature of LRDs is still under debate. In this work, we test the AGN hypothesis for LRDs by stacking the archival {\em Chandra} observations of 19 spectroscopically-selected LRDs. We obtain non-detections in both the soft $(0.5-2\text{ keV})$ and the hard $(2-8\text{ keV})$ X-ray band after stacking, and put upper limits on the X-ray luminosities of these LRDs. We find that the soft band upper limit is $\sim1$ dex lower than the expected level from the $L_\text{X}-L_{\text{H}α}$ relation, while the hard band upper limit is consistent with this relation. This result suggests that LRDs and typical type-I AGNs have significantly different properties. We consider it unlikely that absorption is the main reason for the observed low $L_\text{X}/L_{\text{H}α}$ ratios of LRDs. Alternatively, we discuss two plausible hypotheses: (1) LRDs have intrinsically weak X-ray emissions, (2) fast, galactic-scale outflows have a major contribution to the observed broad H$α$ lines. Our findings indicate that empirical relations (e.g., for black hole mass measurements) established for typical type-I AGNs should be used with caution when analyzing the properties of LRDs. △ Less

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

Comments: Corrected for errors in author information

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

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

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

Submitted 18 April, 2024; originally announced April 2024.

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

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.4}\,{\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 ($\lesssim1\%$); (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 $\approx20\%$. 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 18 March, 2024; originally announced March 2024.

Comments: Submitted to MNRAS; 15 pages, 4 figures

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

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

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

Submitted 15 March, 2024; originally announced March 2024.

Comments: 20 pages, 21 figures

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 our universe, and therefore to reside in proto-clusters encompassing an abundance of galaxies in close vicinity. 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 detecte… ▽ More We expect luminous ($M_{1450}\lesssim-26.5$) high-redshift quasars to trace the highest density peaks in our universe, and therefore to reside in proto-clusters encompassing an abundance of galaxies in close vicinity. 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 $δ>100$ 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}=21.3^{+2.7}_{-2.6}~{\rm cMpc}\,h^{-1}$, while we obtain a correlation length of the [OIII]-emitting galaxies of $r_0^{\rm GG}=4.2\pm0.1~{\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.30\pm0.14$ (and $\log_{10}(M_{\rm halo, min}^{\rm [OIII]}/M_\odot) = 10.72\pm0.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 12 March, 2024; originally announced March 2024.

Comments: Submitted to AAS Journals

A quasar-galaxy merger at $z\sim6.2$: black hole mass and quasar properties from the NIRSpec spectrum

Authors: Federica Loiacono, Roberto Decarli, Marco Mignoli, Emanuele Paolo Farina, Eduardo Bañados, Sarah Bosman, Anna-Christina Eilers, Jan-Torge Schindler, Michael A. Strauss, Marianne Vestergaard, Feige Wang, Laura Blecha, Chris L. Carilli, Andrea Comastri, Thomas Connor, Tiago Costa, Massimo Dotti, Xiaohui Fan, Roberto Gilli, Hyunsung D. Jun, Weizhe Liu, Alessandro Lupi, Madeline A. Marshall, Chiara Mazzucchelli, Romain A. Meyer , et al. (9 additional authors not shown)

Abstract: We present JWST/NIRSpec integral field data of the quasar PJ308-21 at $z=6.2342$. As shown by previous ALMA and HST imaging, the quasar has two companion sources, interacting with the quasar host galaxy. The high-resolution G395H/290LP NIRSpec spectrum covers the $2.87-5.27\ \rm μm$ wavelength range and shows the rest-frame optical emission of the quasar with exquisite quality ($S/N\sim 100-400$ p… ▽ More We present JWST/NIRSpec integral field data of the quasar PJ308-21 at $z=6.2342$. As shown by previous ALMA and HST imaging, the quasar has two companion sources, interacting with the quasar host galaxy. The high-resolution G395H/290LP NIRSpec spectrum covers the $2.87-5.27\ \rm μm$ wavelength range and shows the rest-frame optical emission of the quasar with exquisite quality ($S/N\sim 100-400$ per spectral element). Based on the H$β$ line from the broad line region, we obtain an estimate of the black hole mass $M_{\rm BH,Hβ}\sim 2.7\times 10^{9}\ \rm M_{\odot}$. This value is within a factor $\lesssim 1.5$ of the H$α$-based black hole mass from the same spectrum ($M_{\rm BH, Hα}\sim 1.93\times 10^{9}\ \rm M_{\odot}$) and is consistent with a previous estimate relying on the MgII $λ2799$ ($M_{\rm BH, MgII}\sim 2.65\times 10^{9}\ \rm M_{\odot}$). All these $M_{\rm BH}$ are within the $\sim 0.5$ dex intrinsic scatter of the adopted mass calibrations. The high Eddington ratio of PJ308-21 $λ_{\rm Edd,Hβ}\sim 0.67$ ($λ_{\rm Edd,Hα}\sim 0.96$) is in line with the overall quasar population at $z \gtrsim 6$. The relative strengths of the [OIII], FeII and H$β$ lines are consistent with the empirical "Eigenvector 1" correlations as observed for low redshift quasars. We find evidence for blueshifted [OIII] $λ5007$ emission with a velocity offset $Δv_{\rm [OIII]}=-1922\pm 39$ km s$^{-1}$ from the systemic velocity and a $\rm FWHM([OIII])=2776^{+75}_{-74}$ km s$^{-1}$. This may be the signature of an outflow from the nuclear region, despite the true values of $Δv_{\rm [OIII]}$ and $\rm FWHM([OIII])$ are likely more uncertain due to the blending with H$β$ and FeII lines. Our study demonstrates the unique capabilities of NIRSpec in capturing quasar spectra at cosmic dawn and studying their properties in unprecedented detail. △ Less

Submitted 20 February, 2024; originally announced February 2024.

Comments: 17 pages, 9 figures, accepted for publication in A&A

Journal ref: A&A 685, A121 (2024)

Evidence for Episodic Black Hole Growth of Reionization-Era Quasars observed with Magellan/FIRE

Authors: Leah Bigwood, Anna-Christina Eilers, Robert A. Simcoe

Abstract: Observations of high-redshift quasars hosting billion solar mass black holes at $z\gtrsim6$ challenge our understanding of early supermassive black hole (SMBH) growth. In this work, we conduct a near-infrared spectroscopic study of $19$ quasars at $6.2\lesssim z\lesssim 7.5$, using the Folded-port InfraRed Echellette (FIRE) instrument on the $6.5$-meter Magellan/Baade Telescope. We estimate the si… ▽ More Observations of high-redshift quasars hosting billion solar mass black holes at $z\gtrsim6$ challenge our understanding of early supermassive black hole (SMBH) growth. In this work, we conduct a near-infrared spectroscopic study of $19$ quasars at $6.2\lesssim z\lesssim 7.5$, using the Folded-port InfraRed Echellette (FIRE) instrument on the $6.5$-meter Magellan/Baade Telescope. We estimate the single-epoch masses of the quasars' SMBHs by means of the MgII emission line and find black hole masses of $M_{\text{BH}} \approx(0.2-4.8)\,\times\,10^9\,M_\odot$. Furthermore, we measure the sizes of the quasars' proximity zones, which are regions of enhanced transmitted flux bluewards of the Ly$α\,$ emission line, ionized by the quasars' radiation itself. While it has been shown that the proximity zone sizes correlate with the quasars' lifetimes due to the finite response time of the intergalactic medium to the quasars' radiation, we do not find any correlation between the proximity zone sizes and the black hole mass, which suggests that quasar activity and the concomitant black hole growth are intermittent and episodic. △ Less

Submitted 18 February, 2024; originally announced February 2024.

The THESAN project: galaxy sizes during the epoch of reionization

Authors: Xuejian Shen, Mark Vogelsberger, Josh Borrow, Yongao Hu, Evan Erickson, Rahul Kannan, Aaron Smith, Enrico Garaldi, Lars Hernquist, Takahiro Morishita, Sandro Tacchella, Oliver Zier, Guochao Sun, Anna-Christina Eilers, Hui Wang

Abstract: We investigate galaxy sizes at redshift $z\gtrsim 6$ with the cosmological radiation-magneto-hydrodynamic simulation suite THESAN(-HR). These simulations simultaneously capture the reionization of the large-scale intergalactic medium and resolved galaxy properties. The intrinsic size ($r^{\ast}_{1/2}$) of simulated galaxies increases moderately with stellar mass at… ▽ More We investigate galaxy sizes at redshift $z\gtrsim 6$ with the cosmological radiation-magneto-hydrodynamic simulation suite THESAN(-HR). These simulations simultaneously capture the reionization of the large-scale intergalactic medium and resolved galaxy properties. The intrinsic size ($r^{\ast}_{1/2}$) of simulated galaxies increases moderately with stellar mass at $M_{\ast} \lesssim 10^{8}\,{\rm M}_{\odot}$ and decreases fast at larger masses, resulting in a hump feature at $M_{\ast}\sim 10^{8}\,{\rm M}_{\odot}$ that is insensitive to redshift. Low-mass galaxies are in the initial phase of size growth and are better described by a spherical shell model with feedback-driven gas outflows competing with the cold inflows. In contrast, massive galaxies fit better with the disk formation model. They generally experience a phase of rapid compaction and gas depletion, likely driven by internal disk instability rather than external processes. We identify four compact quenched galaxies in the $(95.5\,{\rm cMpc})^{3}$ volume of THESAN-1 at $z\simeq 6$, and their quenching follows reaching a characteristic stellar surface density akin to the massive compact galaxies at cosmic noon. Compared to observations, we find that the median UV effective radius ($R^{\rm UV}_{\rm eff}$) of simulated galaxies is at least three times larger than the observed ones at $M_{\ast}\lesssim 10^{9}\,{\rm M}_{\odot}$ or $M_{\rm UV}\gtrsim -20$ at $6 \lesssim z \lesssim 10$. This inconsistency, related to the hump feature of the intrinsic size--mass relation, persists across many other cosmological simulations with different galaxy formation models and demonstrates the potential of using galaxy morphology to constrain the physics of galaxy formation at high redshifts. △ Less

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

Comments: 21 pages, 17 figures. Version submitted to MNRAS

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

Chronicling the reionization history at $6\lesssim z \lesssim 7$ with emergent quasar damping wings

Authors: Dominika Ďurovčíková, Anna-Christina Eilers, Huanqing Chen, Sindhu Satyavolu, Girish Kulkarni, Robert A. Simcoe, Laura C. Keating, Martin G. Haehnelt, Eduardo Bañados

Abstract: The spectra of high-redshift ($z\gtrsim 6$) quasars contain valuable information on the progression of the Epoch of Reionization (EoR). At redshifts $z<6$, the observed Lyman-series forest shows that the intergalactic medium (IGM) is nearly ionized, while at $z>7$ the observed quasar damping wings indicate high neutral gas fractions. However, there remains a gap in neutral gas fraction constraints… ▽ More The spectra of high-redshift ($z\gtrsim 6$) quasars contain valuable information on the progression of the Epoch of Reionization (EoR). At redshifts $z<6$, the observed Lyman-series forest shows that the intergalactic medium (IGM) is nearly ionized, while at $z>7$ the observed quasar damping wings indicate high neutral gas fractions. However, there remains a gap in neutral gas fraction constraints at $6\lesssim z \lesssim 7$ where the Lyman series forest becomes saturated but damping wings have yet to fully emerge. In this work, we use a sample of 18 quasar spectra at redshifts $6.0<z<7.1$ to close this gap. We apply neural networks to reconstruct the quasars' continuum emission around the partially absorbed Lyman $α$ line to normalize their spectra, and stack these continuum-normalized spectra in three redshift bins. To increase the robustness of our results, we compare the stacks to a grid of models from two hydrodynamical simulations, ATON and CROC, and we measure the volume-averaged neutral gas fraction, $\bar{x}_{\rm HI}$, while jointly fitting for the mean quasar lifetime, $t_{\rm Q}$, for each stacked spectrum. We chronicle the evolution of neutral gas fraction using the ATON (CROC) models as follows: $\bar{x}_{\rm HI} = 0.21_{-0.07}^{+0.17}$ ($\bar{x}_{\rm HI} = 0.10_{<10^{-4}}^{+0.73}$) at $\langle z \rangle =6.10$, $\bar{x}_{\rm HI} = 0.21_{-0.07}^{+0.33}$ ($\bar{x}_{\rm HI} =0.57_{-0.47}^{+0.26}$) at $\langle z \rangle =6.46$, and $\bar{x}_{\rm HI} = 0.37_{-0.17}^{+0.17}$ ($\bar{x}_{\rm HI} =0.57_{-0.21}^{+0.26}$) at $\langle z \rangle =6.87$. At the same time we constrain the average quasar lifetime to be $t_{\rm Q} \lesssim 7\ {\rm Myr}$ across all redshift bins, in good agreement with previous studies. △ Less

Submitted 18 January, 2024; originally announced January 2024.

Comments: 28 pages, 14 figures

A New Population of Mid-Infrared-Selected Tidal Disruption Events: Implications for Tidal Disruption Event Rates and Host Galaxy Properties

Authors: Megan Masterson, Kishalay De, Christos Panagiotou, Erin Kara, Iair Arcavi, Anna-Christina Eilers, Danielle Frostig, Suvi Gezari, Iuliia Grotova, Zhu Liu, Adam Malyali, Aaron M. Meisner, Andrea Merloni, Megan Newsome, Arne Rau, Robert A. Simcoe, Sjoert van Velzen

Abstract: Most tidal disruption events (TDEs) are currently found in time-domain optical and soft X-ray surveys, both of which are prone to significant obscuration. The infrared (IR), however, is a powerful probe of dust-enshrouded environments, and hence, we recently performed a systematic search of NEOWISE mid-IR data for nearby, obscured TDEs within roughly 200 Mpc. We identified 18 TDE candidates in gal… ▽ More Most tidal disruption events (TDEs) are currently found in time-domain optical and soft X-ray surveys, both of which are prone to significant obscuration. The infrared (IR), however, is a powerful probe of dust-enshrouded environments, and hence, we recently performed a systematic search of NEOWISE mid-IR data for nearby, obscured TDEs within roughly 200 Mpc. We identified 18 TDE candidates in galactic nuclei, using difference imaging to uncover nuclear variability amongst significant host galaxy emission. These candidates were selected based on the following IR light curve properties: (1) $L_\mathrm{W2}\gtrsim10^{42}$ erg s$^{-1}$ at peak, (2) fast rise, followed by a slow, monotonic decline, (3) no significant prior variability, and (4) no evidence for AGN activity in WISE colors. The majority of these sources showed no variable optical counterpart, suggesting that optical surveys indeed miss numerous obscured TDEs. Using narrow line ionization levels and variability arguments, we identified 6 sources as possible underlying AGN, yielding a total of 12 TDEs in our gold sample. This gold sample yields a lower limit on the IR-selected TDE rate of $2.0\pm0.3\times10^{-5}$ galaxy$^{-1}$ year$^{-1}$ ($1.3\pm0.2\times10^{-7}$ Mpc$^{-3}$ year$^{-1}$), which is comparable to optical and X-ray TDE rates. The IR-selected TDE host galaxies do not show a green valley overdensity nor a preference for quiescent, Balmer strong galaxies, which are both overrepresented in optical and X-ray TDE samples. This IR-selected sample represents a new population of dusty TDEs that have historically been missed by optical and X-ray surveys and helps alleviate tensions between observed and theoretical TDE rates and the so-called missing energy problem. △ Less

Submitted 2 January, 2024; originally announced January 2024.

Comments: 28 pages, 13 figures (+ 7 pages, 4 figures in appendix), accepted for publication in ApJ

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

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

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

Submitted 13 December, 2023; originally announced December 2023.

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

ALMA 400pc Imaging of a z=6.5 Massive Warped Disk Galaxy

Authors: Marcel Neeleman, Fabian Walter, Roberto Decarli, Alyssa B. Drake, Anna-Christina Eilers, Romain A. Meyer, Bram P. Venemans

Abstract: We present 0.075 (~400 pc) resolution ALMA observations of the [CII] and dust continuum emission from the host galaxy of the z=6.5406 quasar, P036+03. We find that the emission arises from a thin, rotating disk with an effective radius of 0.21" (1.1 kpc). The velocity dispersion of the disk is consistent with a constant value of 66.4+-1.0 km/s, yielding a scale height of 80+-30 pc. The [CII] veloc… ▽ More We present 0.075 (~400 pc) resolution ALMA observations of the [CII] and dust continuum emission from the host galaxy of the z=6.5406 quasar, P036+03. We find that the emission arises from a thin, rotating disk with an effective radius of 0.21" (1.1 kpc). The velocity dispersion of the disk is consistent with a constant value of 66.4+-1.0 km/s, yielding a scale height of 80+-30 pc. The [CII] velocity field reveals a distortion that we attribute to a warp in the disk. Modeling this warped disk yields an inclination estimate of 40.4+-1.3 degrees and a rotational velocity of 116+-3 km/s. The resulting dynamical mass estimate of (1.96+-0.10) x 10^10 Msun is lower than previous estimates, which strengthens the conclusion that the host galaxy is less massive than expected based on local scaling relations between the black hole mass and the host galaxy mass. Using archival MUSE Ly-alpha observations, we argue that counter-rotating halo gas could provide the torque needed to warp the disk. We further detect a region with excess (15-sigma) dust continuum emission, which is located 1.3 kpc northwest of the galaxy's center and is gravitationally unstable (Toomre-Q < 0.04). We posit this is a star-forming region whose formation was triggered by the warp, because the region is located within a part of the warped disk where gas can efficiently lose angular momentum. The combined ALMA and MUSE imaging provides a unique view of how gas interactions within the disk-halo interface can influence the growth of massive galaxies within the first billion years of the universe. △ Less

Submitted 1 November, 2023; originally announced November 2023.

Comments: Accepted for publication in the Astrophysical Journal

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

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

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

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

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

EIGER V. Characterizing the Host Galaxies of Luminous Quasars at $z\gtrsim6$

Authors: Minghao Yue, Anna-Christina Eilers, Robert A. Simcoe, Ruari Mackenzie, Jorryt Matthee, Daichi Kashino, Rongmon Bordoloi, Simon J. Lilly, Rohan P. Naidu

Abstract: We report {\em JWST}/NIRCam measurements of quasar host galaxy emissions and supermassive black hole (SMBH) masses for six quasars at $5.9<z<7.1$ in the \textit{Emission-line galaxies and Intergalactic Gas in the Epoch of Reionization} (EIGER) project. We obtain deep NIRCam imaging in the F115W, F200W, and F356W bands, as well as F356W grism spectroscopy of the quasars. We use bright unsaturated s… ▽ More We report {\em JWST}/NIRCam measurements of quasar host galaxy emissions and supermassive black hole (SMBH) masses for six quasars at $5.9<z<7.1$ in the \textit{Emission-line galaxies and Intergalactic Gas in the Epoch of Reionization} (EIGER) project. We obtain deep NIRCam imaging in the F115W, F200W, and F356W bands, as well as F356W grism spectroscopy of the quasars. We use bright unsaturated stars to construct models of the point spread function (PSF) and estimate the errors of these PSFs. We then measure or constrain the fluxes and morphology of the quasar host galaxies by fitting the quasar images as a point source plus an exponential disk. We successfully detect the host galaxy of three quasars, which have host-to-quasar flux ratios of $\sim1\%-5\%$. Spectral Energy Distribution (SED) fitting suggests that these quasar host galaxies have stellar masses of $M_*\gtrsim10^{10}M_\odot$. For quasars with host galaxy non-detections, we estimate the upper limits of their stellar masses. We use the grism spectra to measure the {\hb} line profile and the continuum luminosity, then estimate the SMBH masses for the quasars. Our results indicate that the positive relation between SMBH masses and host galaxy stellar masses already exists at redshift $z\gtrsim6$. The quasars in our sample show a high black hole to stellar mass ratio of $M_\text{BH}/M_*\sim0.15$, which is about $\sim2$ dex higher than local relations. We find that selection effects only contribute partially to the high $M_\text{BH}/M_*$ ratios of high-redshift quasars. This result hints at a possible redshift evolution of the $M_\text{BH}-M_*$ relation. △ Less

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

Comments: Accepted by ApJ

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

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

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

Submitted 6 September, 2023; originally announced September 2023.

Comments: 17 pages, 13 figures, submitted to MNRAS

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

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

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

Submitted 8 August, 2023; originally announced August 2023.

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

EIGER IV: The cool 10$^4$K circumgalactic environment of high-$z$ galaxies reveals remarkably efficient IGM enrichment

Authors: Rongmon Bordoloi, Robert A. Simcoe, Jorryt Matthee, Daichi Kashino, Ruari Mackenzie, Simon J. Lilly, Anna-Christina Eilers, Bin Liu, David DePalma, Minghao Yue, Rohan P. Naidu

Abstract: We report new observations of the cool diffuse gas around 29, $2.3<z<6.3$ galaxies, using deep JWST/NIRCam slitless grism spectroscopy around the sightline to the quasar J0100+2802. The galaxies span a stellar mass range of $7.1 \leq \log M_{*}/M_{sun} \leq 10.7$, and star-formation rates of $-0.1 < \log \; SFR/M_{sun}yr^{-1} \; <2.3$. We find galaxies for seven MgII absorption systems within 300… ▽ More We report new observations of the cool diffuse gas around 29, $2.3<z<6.3$ galaxies, using deep JWST/NIRCam slitless grism spectroscopy around the sightline to the quasar J0100+2802. The galaxies span a stellar mass range of $7.1 \leq \log M_{*}/M_{sun} \leq 10.7$, and star-formation rates of $-0.1 < \log \; SFR/M_{sun}yr^{-1} \; <2.3$. We find galaxies for seven MgII absorption systems within 300 kpc of the quasar sightline. The MgII radial absorption profile falls off sharply with radii, with most of the absorption extending out to 2-3$R_{200}$ of the host galaxies. Six out of seven MgII absorption systems are detected around galaxies with $\log M_{*}/M_{sun} >$9. MgII absorption kinematics are shifted from the systemic redshift of host galaxies with a median absolute velocity of 135 km/s and standard deviation of 85 km/s. The high kinematic offset and large radial separation ($R> 1.3 R_{200}$), suggest that five out of the seven MgII absorption systems are gravitationally not bound to the galaxies. In contrast, most cool circumgalactic media at $z<1$ are gravitationally bound. The high incidence of unbound MgII gas in this work suggests that towards the end of reionization, galaxy halos are in a state of remarkable disequilibrium, and are highly efficient in enriching the intergalactic medium. Two strongest MgII absorption systems are detected at $z\sim$ 4.22 and 4.5, the former associated with a merging galaxy system and the latter associated with three kinematically close galaxies. Both these galaxies reside in local galaxy over-densities, indicating the presence of cool MgII absorption in two "proto-groups" at $z>4$. △ Less

Submitted 8 January, 2024; v1 submitted 3 July, 2023; originally announced July 2023.

Comments: 23 pages, 12 figures, 4 tables, Accepted for publication in ApJ

Streamlined Lensed Quasar Identification in Multiband Images via Ensemble Networks

Authors: Irham Taufik Andika, Sherry H. Suyu, Raoul Cañameras, Alejandra Melo, Stefan Schuldt, Yiping Shu, Anna-Christina Eilers, Anton Timur Jaelani, Minghao Yue

Abstract: Quasars experiencing strong lensing offer unique viewpoints on subjects related to the cosmic expansion rate, the dark matter profile within the foreground deflectors, and the quasar host galaxies. Unfortunately, identifying them in astronomical images is challenging since they are overwhelmed by the abundance of non-lenses. To address this, we have developed a novel approach by ensembling cutting… ▽ More Quasars experiencing strong lensing offer unique viewpoints on subjects related to the cosmic expansion rate, the dark matter profile within the foreground deflectors, and the quasar host galaxies. Unfortunately, identifying them in astronomical images is challenging since they are overwhelmed by the abundance of non-lenses. To address this, we have developed a novel approach by ensembling cutting-edge convolutional networks (CNNs) -- for instance, ResNet, Inception, NASNet, MobileNet, EfficientNet, and RegNet -- along with vision transformers (ViTs) trained on realistic galaxy-quasar lens simulations based on the Hyper Suprime-Cam (HSC) multiband images. While the individual model exhibits remarkable performance when evaluated against the test dataset, achieving an area under the receiver operating characteristic curve of $>$97.3% and a median false positive rate of 3.6%, it struggles to generalize in real data, indicated by numerous spurious sources picked by each classifier. A significant improvement is achieved by averaging these CNNs and ViTs, resulting in the impurities being downsized by factors up to 50. Subsequently, combining the HSC images with the UKIRT, VISTA, and unWISE data, we retrieve approximately 60 million sources as parent samples and reduce this to 892,609 after employing a photometry preselection to discover $z>1.5$ lensed quasars with Einstein radii of $θ_\mathrm{E}<5$ arcsec. Afterward, the ensemble classifier indicates 3080 sources with a high probability of being lenses, for which we visually inspect, yielding 210 prevailing candidates awaiting spectroscopic confirmation. These outcomes suggest that automated deep learning pipelines hold great potential in effectively detecting strong lenses in vast datasets with minimal manual visual inspection involved. △ Less

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

Comments: Accepted for publication in the Astronomy & Astrophysics journal. 28 pages, 11 figures, and 3 tables. We welcome comments from the reader

Journal ref: A&A 678, A103 (2023)

Quaia, the Gaia-unWISE Quasar Catalog: An All-Sky Spectroscopic Quasar Sample

Authors: Kate Storey-Fisher, David W. Hogg, Hans-Walter Rix, Anna-Christina Eilers, Giulio Fabbian, Michael Blanton, David Alonso

Abstract: We present a new, all-sky quasar catalog, Quaia, that samples the largest comoving volume of any existing spectroscopic quasar sample. The catalog draws on the 6,649,162 quasar candidates identified by the Gaia mission that have redshift estimates from the space observatory's low-resolution BP/RP spectra. This initial sample is highly homogeneous and complete, but has low purity, and 18% of even t… ▽ More We present a new, all-sky quasar catalog, Quaia, that samples the largest comoving volume of any existing spectroscopic quasar sample. The catalog draws on the 6,649,162 quasar candidates identified by the Gaia mission that have redshift estimates from the space observatory's low-resolution BP/RP spectra. This initial sample is highly homogeneous and complete, but has low purity, and 18% of even the bright ($G<20.0$) confirmed quasars have discrepant redshift estimates ($|Δz/(1+z)|>0.2$) compared to those from the Sloan Digital Sky Survey (SDSS). In this work, we combine the Gaia candidates with unWISE infrared data (based on the Wide-field Infrared Survey Explorer survey) to construct a catalog useful for cosmological and astrophysical quasar studies. We apply cuts based on proper motions and Gaia and unWISE colors, reducing the number of contaminants by $\sim$4$\times$. We improve the redshifts by training a $k$-nearest neighbors model on SDSS redshifts, and achieve estimates on the $G<20.0$ sample with only 6% (10%) catastrophic errors with $|Δz/(1+z)|>0.2$ ($0.1$), a reduction of $\sim$3$\times$ ($\sim$2$\times$) compared to the Gaia redshifts. The final catalog has 1,295,502 quasars with $G<20.5$, and 755,850 candidates in an even cleaner $G<20.0$ sample, with accompanying rigorous selection function models. We compare Quaia to existing quasar catalogs, showing that its large effective volume makes it a highly competitive sample for cosmological large-scale structure analyses. The catalog is publicly available at https://zenodo.org/records/10403370. △ Less

Submitted 18 March, 2024; v1 submitted 30 June, 2023; originally announced June 2023.

Comments: Published in the Astrophysical Journal; updated to match published version. Catalog available at https://doi.org/10.5281/zenodo.10403370. Animation available at https://cosmo.nyu.edu/ksf/quaia.mp4

Journal ref: Storey-Fisher (2024) 964, 69

Constraining cosmology with the Gaia-unWISE Quasar Catalog and CMB lensing: structure growth

Authors: David Alonso, Giulio Fabbian, Kate Storey-Fisher, Anna-Christina Eilers, Carlos García-García, David W. Hogg, Hans-Walter Rix

Abstract: We study the angular clustering of Quaia, a Gaia- and unWISE-based catalog of over a million quasars with an exceptionally well-defined selection function. With it, we derive cosmology constraints from the amplitude and growth of structure across cosmic time. We divide the sample into two redshift bins, centered at $z=1.0$ and $z=2.1$, and measure both overdensity auto-correlations and cross-corre… ▽ More We study the angular clustering of Quaia, a Gaia- and unWISE-based catalog of over a million quasars with an exceptionally well-defined selection function. With it, we derive cosmology constraints from the amplitude and growth of structure across cosmic time. We divide the sample into two redshift bins, centered at $z=1.0$ and $z=2.1$, and measure both overdensity auto-correlations and cross-correlations with maps of the Cosmic Microwave Background convergence measured by Planck. From these data, and including a prior from measurements of the baryon acoustic oscillations scale, we place constraints on the amplitude of the matter power spectrum $σ_8=0.766\pm 0.034$, and on the matter density parameter $Ω_m=0.343^{+0.017}_{-0.019}$. These measurements are in reasonable agreement with \planck at the $\sim$ 1.4$σ$ level, and are found to be robust with respect to observational and theoretical uncertainties. We find that our slightly lower value of $σ_8$ is driven by the higher-redshift sample, which favours a low amplitude of matter fluctuations. We present plausible arguments showing that this could be driven by contamination of the CMB lensing map by high-redshift extragalactic foregrounds, which should also affect other cross-correlations with tracers of large-scale structure beyond $z\sim1.5$. Our constraints are competitive with those from state-of-the-art 3$\times$2-point analyses, but arise from a range of scales and redshifts that is highly complementary to those covered by cosmic shear data and most galaxy clustering samples. This, coupled with the unprecedented combination of volume and redshift precision achieved by Quaia allows us to break the usual degeneracy between $Ω_m$ and $σ_8$. △ Less

Submitted 3 July, 2023; v1 submitted 30 June, 2023; originally announced June 2023.

Comments: Submitted to JCAP, corrects links to public data release. The per-bin selection functions used in this analysis, and the complete set of power spectra, covariance, and associated metadata used in our fiducial analysis are publicly available at https://doi.org/10.5281/zenodo.8098635 . The public Quaia catalog can be found at https://doi.org/10.5281/zenodo.8060755

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

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

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

Submitted 28 June, 2023; originally announced June 2023.

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

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

Little Red Dots: an abundant population of faint AGN at z~5 revealed by the EIGER and FRESCO JWST surveys

Authors: Jorryt Matthee, Rohan P. Naidu, Gabriel Brammer, John Chisholm, Anna-Christina Eilers, Andy Goulding, Jenny Greene, Daichi Kashino, Ivo Labbe, Simon J. Lilly, Ruari Mackenzie, Pascal A. Oesch, Andrea Weibel, Stijn Wuyts, Mengyuan Xiao, Rongmon Bordoloi, Rychard Bouwens, Pieter van Dokkum, Garth Illingworth, Ivan Kramarenko, Michael V. Maseda, Charlotte Mason, Romain A. Meyer, Erica J. Nelson, Naveen A. Reddy , et al. (3 additional authors not shown)

Abstract: Characterising the prevalence and properties of faint active galactic nuclei (AGN) in the early Universe is key for understanding the formation of supermassive black holes (SMBHs) and determining their role in cosmic reionization. We perform a spectroscopic search for broad H$α$ emitters at $z\approx4-6$ using deep JWST/NIRCam imaging and wide field slitless spectroscopy from the EIGER and FRESCO… ▽ More Characterising the prevalence and properties of faint active galactic nuclei (AGN) in the early Universe is key for understanding the formation of supermassive black holes (SMBHs) and determining their role in cosmic reionization. We perform a spectroscopic search for broad H$α$ emitters at $z\approx4-6$ using deep JWST/NIRCam imaging and wide field slitless spectroscopy from the EIGER and FRESCO surveys. We identify 20 H$α$ lines at $z=4.2-5.5$ that have broad components with line widths from $\sim1200-3700$ km s$^{-1}$, contributing $\sim30-90$ % of the total line flux. We interpret these broad components as being powered by accretion onto SMBHs with implied masses $\sim10^{7-8}$ M$_{\odot}$. In the UV luminosity range M$_{\rm UV}=-21$ to $-18$, we measure number densities of $\approx10^{-5}$ cMpc$^{-3}$. This is an order of magnitude higher than expected from extrapolating quasar UV luminosity functions. Yet, such AGN are found in only $<1$ % of star-forming galaxies at $z\sim5$. The SMBH mass function agrees with large cosmological simulations. In two objects we detect narrow red- and blue-shifted H$α$ absorption indicative, respectively, of dense gas fueling SMBH growth and outflows. We may be witnessing early AGN feedback that will clear dust-free pathways through which more massive blue quasars are seen. We uncover a strong correlation between reddening and the fraction of total galaxy luminosity arising from faint AGN. This implies that early SMBH growth is highly obscured and that faint AGN are only minor contributors to cosmic reionization. △ Less

Submitted 1 February, 2024; v1 submitted 8 June, 2023; originally announced June 2023.

Comments: 26 pages, 20 figures. Accepted for publication in ApJ. Various improvements and updated citations. Main additions: colours (Fig 11), broad Halpha LF (Fig 18, Sect 5.2), comparisons to simulations (Sect 5.3)

Compact [C II] emitters around a C IV absorption complex at redshift 5.7

Authors: Daichi Kashino, Simon J. Lilly, Robert A. Simcoe, Rongmon Bordoloi, Ruari Mackenzie, Jorryt Matthee, Anna-Christina Eilers

Abstract: The physical conditions of the circumgalactic medium are probed by intervening absorption-line systems in the spectrum of background quasi-stellar objects out to the epoch of cosmic reionization. A correlation between the ionization state of the absorbing gas and the nature of the nearby galaxies has been suggested by the sources detected either in Lyalpha or [C ii] 158 m near to respectively high… ▽ More The physical conditions of the circumgalactic medium are probed by intervening absorption-line systems in the spectrum of background quasi-stellar objects out to the epoch of cosmic reionization. A correlation between the ionization state of the absorbing gas and the nature of the nearby galaxies has been suggested by the sources detected either in Lyalpha or [C ii] 158 m near to respectively highly-ionized and neutral absorbers. This is also likely linked to the global changes in the incidence of absorption systems of different types and the process of cosmic reionization. Here we report the detection of two [C ii]-emitting galaxies at redshift $z \sim 5.7$ that are associated with a complex high-ionization C iv absorption system. These objects are part of an overdensity of galaxies and have compact sizes (< 2.4 kpc) and narrow line widths (FWHM $\sim$ 62--64 km s-1). Hydrodynamic simulations predict that similar narrow [C ii] emission may arise from the heating of small ($\lesssim$ 3 kpc) clumps of cold neutral medium or a compact photodissociation region. The lack of counterparts in the rest-frame ultraviolet indicates severe obscuration of the sources that are exciting the [C ii] emission. These results may suggest a connection between the properties of the [C ii] emission, the rare overdensity of galaxies and the unusual high ionization state of the gas in this region. △ Less

Submitted 22 May, 2023; originally announced May 2023.

Comments: Published in Nature on 10 May 2023; authors' version; link to the paper: https://www.nature.com/articles/s41586-023-05901-3

Journal ref: Nature 617, 261-264 (2023)

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

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

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

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

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

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

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

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

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

Submitted 1 May, 2023; originally announced May 2023.

Comments: 16 pages, 9 figures, Accepted in MNRAS

Report number: TIFR/TH/23-5

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

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

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

Submitted 19 April, 2023; originally announced April 2023.

Comments: accepted for publication in ApJL

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

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

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

Submitted 19 April, 2023; originally announced April 2023.

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

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

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

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

Submitted 18 April, 2023; originally announced April 2023.

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

Measuring the photo-ionization rate, neutral fraction and mean free path of HI ionizing photons at $4.9 \leq z \leq 6.0$ from a large sample of XShooter and ESI spectra

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

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

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

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

The dark matter profile of the Milky Way inferred from its circular velocity curve

Authors: Xiaowei Ou, Anna-Christina Eilers, Lina Necib, Anna Frebel

Abstract: In this paper, we construct the circular velocity curve of the Milky Way out to $\sim 30$ kpc, providing an updated model of the dark matter density profile. We derive precise parallaxes for 120,309 stars with a data-driven model, using APOGEE DR17 spectra combined with Gaia DR3, 2MASS, and WISE photometry. At outer galactic radii up to 30 kpc, we find a significantly faster decline in the circula… ▽ More In this paper, we construct the circular velocity curve of the Milky Way out to $\sim 30$ kpc, providing an updated model of the dark matter density profile. We derive precise parallaxes for 120,309 stars with a data-driven model, using APOGEE DR17 spectra combined with Gaia DR3, 2MASS, and WISE photometry. At outer galactic radii up to 30 kpc, we find a significantly faster decline in the circular velocity curve compared to the inner parts. This decline is better fit with a cored Einasto profile with a slope parameter $0.91^{+0.04}_{-0.05}$ than a generalized Navarro-Frenk-White (NFW) profile. The virial mass of the best-fit dark matter halo profile is only $1.81^{+0.06}_{-0.05}\times10^{11}$ $M_{\odot}$, significantly lower than what a generalized NFW profile delivers. We present a study of the potential systematics, affecting mainly large radii. Such a low mass for the Galaxy is driven by the functional forms tested, given that it probes beyond our measurements. It is found to be in tension with mass measurements from globular clusters, dwarf satellites, and streams. Our best-fit profile also lowers the expected dark matter annihilation signal flux from the galactic centre by more than an order of magnitude, compared to an NFW profile-fit. In future work, we will explore profiles with more flexible functional forms to more fully leverage the circular velocity curve and observationally constrain the properties of the Milky Way's dark matter halo. △ Less

Submitted 26 May, 2023; v1 submitted 22 March, 2023; originally announced March 2023.

Comments: 17 pages, 11 figures. Fixed typos, updated values, and included additional analyses

SPT-CL J2215-3537: A Massive Starburst at the Center of the Most Distant Relaxed Galaxy Cluster

Authors: Michael S. Calzadilla, Lindsey E. Bleem, Michael McDonald, Michael D. Gladders, Adam B. Mantz, Steven W. Allen, Matthew B. Bayliss, Anna-Christina Eilers, Benjamin Floyd, Julie Hlavacek-Larrondo, Gourav Khullar, Keunho J. Kim, Guillaume Mahler, Keren Sharon, Taweewat Somboonpanyakul, Brian Stalder, Antony A. Stark

Abstract: We present the discovery of the most distant, dynamically relaxed cool core cluster, SPT-CL J2215-3537 (SPT2215) and its central brightest cluster galaxy (BCG) at z=1.16. Using new X-ray observations, we demonstrate that SPT2215 harbors a strong cool core, with a central cooling time of 200 Myr (at 10 kpc) and a maximal intracluster medium cooling rate of 1900+/-400 Msun/yr. This prodigious coolin… ▽ More We present the discovery of the most distant, dynamically relaxed cool core cluster, SPT-CL J2215-3537 (SPT2215) and its central brightest cluster galaxy (BCG) at z=1.16. Using new X-ray observations, we demonstrate that SPT2215 harbors a strong cool core, with a central cooling time of 200 Myr (at 10 kpc) and a maximal intracluster medium cooling rate of 1900+/-400 Msun/yr. This prodigious cooling may be responsible for fueling extended, star-forming filaments observed in Hubble Space Telescope imaging. Based on new spectrophotometric data, we detect bright [O II] emission in the BCG, implying an unobscured star formation rate (SFR) of 320^{+230}_{-140} Msun/yr. The detection of a weak radio source (2.0+/-0.8 mJy at 0.8 GHz) suggests ongoing feedback from an active galactic nucleus (AGN), though the implied jet power is less than half the cooling luminosity of the hot gas, consistent with cooling overpowering heating. The extreme cooling and SFR of SPT2215 is rare among known cool core clusters, and it is even more remarkable that we observe these at such high redshift, when most clusters are still dynamically disturbed. The high mass of this cluster, coupled with the fact that it is dynamically relaxed with a highly-isolated BCG, suggests that it is an exceptionally rare system that must have formed very rapidly in the early Universe. Combined with the high SFR, SPT2215 may be a high-z analog of the Phoenix cluster, potentially providing insight into the limits of AGN feedback and star formation in the most massive galaxies. △ Less

Submitted 17 March, 2023; originally announced March 2023.

Comments: 14 pages, 5 figures, 1 table. Accepted in ApJ

Rapidly evolving Galactic plane outbursts in NEOWISE: Revisiting the Galactic nova rate with the first all-sky search in the mid-infrared

Authors: Liam Zuckerman, Kishalay De, Anna-Christina Eilers, Aaron M. Meisner, Christos Panagiotou

Abstract: The Galactic nova rate is intimately linked to our understanding of its chemical enrichment and progenitor channels of Type Ia supernovae. Yet past estimates have varied by more than an order of magnitude ($\approx10-300$ yr$^{-1}$) owing to limitations in both discovery methods as well as assumptions regarding the Galactic dust distribution and extragalactic stellar populations. Recent estimates… ▽ More The Galactic nova rate is intimately linked to our understanding of its chemical enrichment and progenitor channels of Type Ia supernovae. Yet past estimates have varied by more than an order of magnitude ($\approx10-300$ yr$^{-1}$) owing to limitations in both discovery methods as well as assumptions regarding the Galactic dust distribution and extragalactic stellar populations. Recent estimates utilizing synoptic near-infrared surveys have begun to provide a glimpse of a consensus ($\approx25-50$ yr$^{-1}$); however, a consistent estimate remains lacking. Here, we present the first all-sky search for Galactic novae using 8 years of data from the NEOWISE mid-infrared (MIR) survey. Operating at $3.4$ and $4.6$ $μ$m where interstellar extinction is negligible, the 6-month cadence NEOWISE dataset offers unique sensitivity to discover slowly evolving novae across the entire Galaxy. Using a novel image subtraction pipeline together with systematic selection criteria, we identify a sample of 49 rapidly evolving MIR outbursts as candidate Galactic novae. While 27 of these sources are known novae, the remaining are previously missed nova candidates discovered in this work. The unknown events are spatially clustered along the densest and most heavily obscured regions of the Galaxy where previous novae are severely underrepresented. We use simulations of the NEOWISE survey strategy, the pipeline detection efficiency, and our criteria to derive a Galactic nova rate of $47.9^{+3.1}_{-8.3}$ yr$^{-1}$. The discovery of these exceptionally bright (yet overlooked) nova candidates confirm emerging suggestions that optical surveys have been highly incomplete in searches for Galactic novae, highlighting the potential for MIR searches in revealing the demographics of Galactic stellar outbursts. △ Less

Submitted 18 March, 2023; v1 submitted 15 March, 2023; originally announced March 2023.

Comments: 13 pages, 4 figures, 1 table, 3 figures in Appendix. Submitted to MNRAS. Comments welcome

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

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

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

Submitted 5 March, 2023; originally announced March 2023.

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

A luminous dust-obscured Tidal Disruption Event candidate in a star forming galaxy at 42 Mpc

Authors: Christos Panagiotou, Kishalay De, Megan Masterson, Erin Kara, Michael Calzadilla, Anna-Christina Eilers, Danielle Frostig, Viraj Karambelkar, Mansi Kasliwal, Nathan Lourie, Aaron M. Meisner, Robert A. Simcoe, Robert Stein, Jeffry Zolkower

Abstract: While the vast majority of Tidal Disruption Events (TDEs) has been identified by wide-field sky surveys in the optical and X-ray bands, recent studies indicate that a considerable fraction of TDEs may be dust obscured, and thus preferentially detected in the infrared (IR) wavebands. In this Letter, we present the discovery of a luminous mid-IR nuclear flare (termed WTP 14adbjsh) identified in a sy… ▽ More While the vast majority of Tidal Disruption Events (TDEs) has been identified by wide-field sky surveys in the optical and X-ray bands, recent studies indicate that a considerable fraction of TDEs may be dust obscured, and thus preferentially detected in the infrared (IR) wavebands. In this Letter, we present the discovery of a luminous mid-IR nuclear flare (termed WTP 14adbjsh) identified in a systematic transient search of archival images from the NEOWISE mid-IR survey. The source reached a peak luminosity of $L \simeq 10^{43} \text{erg s}^{-1}$ at 4.6 $μ$m in 2015, before fading in the IR with a TDE-like $F \propto t^{-5/3}$ decline, radiating a total of more than $ 3\times 10^{51}$ erg in the last 7 years. The transient event took place in the nearby galaxy NGC 7392, at a distance of around 42 Mpc; yet, no optical or X-ray flare is detected. We interpret the transient as the nearest TDE candidate detected in the last decade, which was missed at other wavelengths due to dust obscuration, hinting at the existence of TDEs that have been historically overlooked. Unlike most previously detected TDEs, the transient was discovered in a star forming galaxy, corroborating earlier suggestions that dust obscuration suppresses significantly the detection of TDEs in these environments. Our results demonstrate that the study of IR-detected TDEs is critical in order to obtain a complete understanding of the physics of TDEs, and to conclude whether TDEs occur preferentially in a particular class of galaxies. △ Less

Submitted 5 March, 2023; originally announced March 2023.

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

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

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

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

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

Comments: Accepted for publication in ApJ

EIGER III. JWST/NIRCam observations of the ultra-luminous high-redshift quasar J0100+2802

Authors: Anna-Christina Eilers, Robert A. Simcoe, Minghao Yue, Ruari Mackenzie, Jorryt Matthee, Dominika Durovcikova, Daichi Kashino, Rongmon Bordoloi, Simon J. Lilly

Abstract: We present the first rest-frame optical spectrum of a high-redshift quasar observed with JWST/NIRCam in Wide Field Slitless (WFSS) mode. The observed quasar, J0100+2802, is the most luminous quasar known at $z>6$. We measure the mass of the central supermassive black hole (SMBH) by means of the rest-frame optical H$β$ emission line, and find consistent mass measurements of the quasar's SMBH of… ▽ More We present the first rest-frame optical spectrum of a high-redshift quasar observed with JWST/NIRCam in Wide Field Slitless (WFSS) mode. The observed quasar, J0100+2802, is the most luminous quasar known at $z>6$. We measure the mass of the central supermassive black hole (SMBH) by means of the rest-frame optical H$β$ emission line, and find consistent mass measurements of the quasar's SMBH of $M_\bullet\approx10^{10}\,M_\odot$ when compared to the estimates based on the properties of rest-frame UV emission lines CIV and MgII, which are accessible from ground-based observatories. To this end, we also present a newly reduced rest-frame UV spectrum of the quasar observed with X-Shooter/VLT and FIRE/Magellan for a total of 16.8 hours. We readdress the question whether this ultra-luminous quasar could be effected by strong gravitational lensing making use of the diffraction limited NIRCam images in three different wide band filters (F115W, F200W, F356W), which improves the achieved spatial resolution compared to previous images taken with the Hubble Space Telescope by a factor of two. We do not find any evidence for a foreground deflecting galaxy, nor for multiple images of the quasar, and determine the probability for magnification due to strong gravitational lensing with image separations below the diffraction limit of $Δθ\lesssim 0.05''$ to be $\lesssim 2.2\times 10^{-3}$. Our observations therefore confirm that this quasar hosts a ten billion solar mass black hole less than $1$ Gyr after the Big Bang, which is challenging to explain with current black hole formation models. △ Less

Submitted 18 May, 2023; v1 submitted 29 November, 2022; originally announced November 2022.

Comments: Accepted for publication in ApJ

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

Submitted 22 November, 2022; originally announced November 2022.

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

EIGER II. first spectroscopic characterisation of the young stars and ionised gas associated with strong H$β$ and [OIII] line-emission in galaxies at z=5-7 with JWST

Authors: Jorryt Matthee, Ruari Mackenzie, Robert A. Simcoe, Daichi Kashino, Simon J. Lilly, Rongmon Bordoloi, Anna-Christina Eilers

Abstract: We present emission-line measurements and physical interpretations for a sample of 117 [OIII] emitting galaxies at $z=5.33-6.93$, using the first deep JWST/NIRCam wide field slitless spectroscopic observations. Our 9.7-hour integration is centered upon the $z=6.3$ quasar J0100+2802 -- the first of six fields targeted by the EIGER survey -- and covers $λ=3-4$ microns. We detect 133 [OIII] doublets,… ▽ More We present emission-line measurements and physical interpretations for a sample of 117 [OIII] emitting galaxies at $z=5.33-6.93$, using the first deep JWST/NIRCam wide field slitless spectroscopic observations. Our 9.7-hour integration is centered upon the $z=6.3$ quasar J0100+2802 -- the first of six fields targeted by the EIGER survey -- and covers $λ=3-4$ microns. We detect 133 [OIII] doublets, but merge pairs within $\approx$10 kpc and 600 km s$^{-1}$, motivated by their small scale clustering excess. We detect H$β$ in 68 and H$γ$ emission in two galaxies. The galaxies are characterised by a UV luminosity M$_{\rm UV}\sim-19.6$ ($-17.7$ to $-22.3$), stellar mass ~$10^8$ $(10^{6.8-10.1})$ M$_{\odot}$, H$β$ and [OIII] EWs $\approx$ 850 Angstrom (up to 3000 Angstrom), young ages (~100 Myr), a highly excited interstellar medium ([OIII]/H$β\approx6$) and low dust attenuations. These high EWs are very rare in the local Universe, but we show they are ubiquitous at $z\sim6$ based on the measured number densities. The stacked spectrum reveals H$γ$ and [OIII]$_{4364}$ which shows that the galaxies are typically dust and metal poor (E(B-V)=0.1, 12+log(O/H)=7.4) with a high electron temperature ($2\times10^4$ K) and a production efficiency of ionising photons ($ξ_{\rm ion}=10^{25.3}$ Hz erg$^{-1}$). We further show the existence of a strong mass-metallicity relation. The young highly ionising stellar populations, moderately low metallicities, low dust attenuations and high ionisation state in z~6 galaxies conspire to maximise the [OIII] output from galaxies, yielding an [OIII] luminosity density at z~6 that is significantly higher than at z~2, despite the order of magnitude decline in cosmic star formation. Thus, [OIII] emission-line surveys with JWST prove a highly efficient method to trace the galaxy density in the epoch of reionization. △ Less

Submitted 17 March, 2023; v1 submitted 15 November, 2022; originally announced November 2022.

Comments: Accepted for publication in ApJ. Main text 22 pages, 20 figures. Main results in Figs 14 (Xi_ion), 15 (MEx diagram),17 (MZR), 19 ([OIII] luminosity density)

EIGER I. a large sample of [OIII]-emitting galaxies at $5.3 < z < 6.9$ and direct evidence for local reionization by galaxies

Authors: Daichi Kashino, Simon J. Lilly, Jorryt Matthee, Anna-Christina Eilers, Ruari Mackenzie, Rongmon Bordoloi, Robert A. Simcoe

Abstract: We present a first sample of 117 [OIII]$λλ$4960,5008-selected star-forming galaxies at $5.33 < z < 6.93$ detected in JWST/NIRCam 3.5$μ$m slitless spectroscopy of a $6.5 \times 3.4$ arcmin$^2$ field centered on the hyperluminous quasar SDSS J0100+2802, obtained as part of the EIGER (Emission-line galaxies and Intergalactic Gas in the Epoch of Reionization) survey. Three prominent galaxy overdensiti… ▽ More We present a first sample of 117 [OIII]$λλ$4960,5008-selected star-forming galaxies at $5.33 < z < 6.93$ detected in JWST/NIRCam 3.5$μ$m slitless spectroscopy of a $6.5 \times 3.4$ arcmin$^2$ field centered on the hyperluminous quasar SDSS J0100+2802, obtained as part of the EIGER (Emission-line galaxies and Intergalactic Gas in the Epoch of Reionization) survey. Three prominent galaxy overdensities are observed, one of them at the redshift of the quasar. Galaxies are found within 200 pkpc and 105 km s$^{-1}$ of four known metal absorption-line systems in this redshift range. We focus on the role of the galaxies in ionizing the surrounding intergalactic medium (IGM) during the later stages of cosmic reionization and construct the mean Ly$α$ and Ly$β$ transmission as a function of distance from the galaxies. At the lowest redshifts in our study, $5.3 < z < 5.7$, the IGM transmission rises monotonically with distance from the galaxies. This is as expected when galaxies reside at peaks in the overdensity field of an IGM that is ionized by more or less uniform ionizing background, and has been seen at lower redshifts. In contrast, at $5.7 < z < 6.14$, the transmission of both Ly$α$ and Ly$β$ first increases with distance, but then peaks at a distance of 5 cMpc before declining. This peak in transmission is qualitatively similar to that seen (albeit at smaller distances and higher redshifts) in the THESAN simulations. Finally, in the region $6.15 < z < 6.26$ where the additional ionizing radiation from the quasar dominates, the monotonic increase in transmission with distance is re-established. This result is interpreted to represent evidence that the transmission of the IGM at $z \sim 5.9$ towards J0100+2802 results from the ``local'' ionizing radiation of galaxies that dominates over the much reduced cosmic background. △ Less

Submitted 15 November, 2022; originally announced November 2022.

Comments: 24 pages, 14 figures, submitted to ApJ, comments welcome

A Generative Model for Quasar Spectra

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

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

Submitted 6 September, 2022; originally announced September 2022.

Comments: accepted for publication in ApJ

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

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

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

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

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