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

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

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

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

Submitted 20 June, 2024; originally announced June 2024.

Comments: 24 pages, 10 figures

Quantifying Lyman-$α$ emissions from reionization fronts

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

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

Submitted 20 June, 2024; originally announced June 2024.

Comments: 21 pages, 6 figures

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

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

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

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

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

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

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

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

Submitted 14 May, 2024; originally announced May 2024.

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

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

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

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

Submitted 18 April, 2024; originally announced April 2024.

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

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

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

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

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

Comments: Accepted by A&A

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

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

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

Submitted 15 March, 2024; originally announced March 2024.

Comments: 20 pages, 21 figures

Chemical Diagnostics to Unveil Environments Enriched by First Stars

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

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

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

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

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

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

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

Submitted 13 December, 2023; originally announced December 2023.

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

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

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

Abstract: Recent measurements of the ionizing photon mean free path (MFP) based on composite quasar spectra may point to reionization ending at $z<6$. These measurements are challenging because they rely on assumptions about the proximity zones of the quasars. For example, some quasars might have been close to neutral patches where reionization was still ongoing ("neutral islands"), and it is unclear how th… ▽ More Recent measurements of the ionizing photon mean free path (MFP) based on composite quasar spectra may point to reionization ending at $z<6$. These measurements are challenging because they rely on assumptions about the proximity zones of the quasars. For example, some quasars might have been close to neutral patches where reionization was still ongoing ("neutral islands"), and it is unclear how they would affect the measurements. We address this question with mock MFP measurements from radiative transfer simulations. We find that, even in the presence of neutral islands, our mock MFP measurements agree to within $30~\%$ with the true spatially averaged MFP in our simulations, which includes opacity from both the ionized medium and the islands. The inferred MFP is sensitive at the $<~50\%$ level to assumptions about quasar environments and lifetimes for realistic models. We demonstrate that future analyses with improved data may require explicitly modeling the effects of neutral islands on the composite spectra, and we outline a method for doing this. Lastly, we quantify the effects of neutral islands on Lyman-series transmission, which has been modeled with optically thin simulations in previous MFP analyses. Neutral islands can suppress transmission at $λ_{\rm rest}<912$ Angstroms significantly, up to a factor of 2 for $z_{\rm qso}=6$ in a plausible reionization scenario, owing to absorption by many closely spaced lines as quasar light redshifts into resonance. However, the suppression is almost entirely degenerate with the spectrum normalization and thus does not significantly bias the inferred MFP. △ Less

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

Comments: Accepted to MNRAS, minor changes

Extreme mass-ratio inspiral and waveforms for a spinning body into a Kerr black hole via osculating geodesics and near-identity transformations

Authors: Lisa V. Drummond, Philip Lynch, Alexandra G. Hanselman, Devin R. Becker, Scott A. Hughes

Abstract: Understanding the orbits of spinning bodies in curved spacetime is important for modeling binary black hole systems with small mass ratios. At zeroth order in mass ratio, the smaller body moves on a geodesic. Post-geodesic effects are needed to model the system accurately. One very important post-geodesic effect is the gravitational self-force, which describes the small body's interaction with its… ▽ More Understanding the orbits of spinning bodies in curved spacetime is important for modeling binary black hole systems with small mass ratios. At zeroth order in mass ratio, the smaller body moves on a geodesic. Post-geodesic effects are needed to model the system accurately. One very important post-geodesic effect is the gravitational self-force, which describes the small body's interaction with its own contribution to a binary's spacetime. Another post-geodesic effect, the spin-curvature force, is due to the smaller body's spin coupling to spacetime curvature. In this paper, we combine the leading orbit-averaged backreaction of point-particle gravitational-wave emission with the spin-curvature force to construct the worldline and gravitational waveform for a spinning body spiraling into a Kerr black hole. We use an osculating geodesic integrator, which treats the worldline as evolution through a sequence of geodesic orbits, as well as near-identity transformations, which eliminate dependence on orbital phases, allowing for fast computation of inspirals. The resulting inspirals and waveforms include all critical dynamical effects which govern such systems (orbit and precession frequencies, inspiral, strong-field gravitational-wave amplitudes), and as such form an effective first model for the inspiral of spinning bodies into Kerr black holes. We emphasize that our present calculation is not self consistent, since we neglect effects which enter at the same order as effects we include. However, our analysis demonstrates that the impact of spin-curvature forces can be incorporated into EMRI waveform tools with relative ease. The calculation is sufficiently modular that it should not be difficult to include neglected post-geodesic effects as efficient tools for computing them become available. (Abridged) △ Less

Submitted 12 October, 2023; originally announced October 2023.

Comments: 38 pages, including 17 pages of Appendices, 13 figures. Submitted to Physics Review D. This posting and submission supersedes arXiv:2305.08919, in response to helpful critical referee feedback

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

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

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

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

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

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

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

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

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

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

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

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

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

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

Submitted 6 September, 2023; originally announced September 2023.

Comments: 17 pages, 13 figures, submitted to MNRAS

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

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

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

Submitted 24 August, 2023; originally announced August 2023.

Comments: 20 pages, 17 figures. Accepted to ApJ

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

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

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

Submitted 8 August, 2023; originally announced August 2023.

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

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

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

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

Submitted 28 June, 2023; originally announced June 2023.

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

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

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

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

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

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

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

Extreme mass-ratio inspiral of a spinning body into a Kerr black hole I: Evolution along generic trajectories

Authors: Lisa V. Drummond, Alexandra G. Hanselman, Devin R. Becker, Scott A. Hughes

Abstract: The study of spinning bodies moving in curved spacetime has relevance to binary black hole systems with large mass ratios, as well as being of formal interest. At zeroth order in a binary's mass ratio, the smaller body moves on a geodesic of the larger body's spacetime. Post-geodesic corrections describing forces driving the small body's worldline away from geodesics must be incorporated to model… ▽ More The study of spinning bodies moving in curved spacetime has relevance to binary black hole systems with large mass ratios, as well as being of formal interest. At zeroth order in a binary's mass ratio, the smaller body moves on a geodesic of the larger body's spacetime. Post-geodesic corrections describing forces driving the small body's worldline away from geodesics must be incorporated to model the system accurately. An important post-geodesic effect is the gravitational self-force, which describes the small body's interaction with its own spacetime curvature. This effect includes the backreaction due to gravitational-wave emission that leads to the inspiral of the small body into the black hole. When a spinning body orbits a black hole, its spin couples to spacetime curvature. This introduces another post-geodesic correction known as the spin-curvature force. An osculating geodesic integrator that includes both the backreaction due to gravitational waves and spin-curvature forces can be used to generate a spinning-body inspiral. In this paper, we use an osculating geodesic integrator to combine the leading backreaction of gravitational waves with the spin-curvature force. Our analysis only includes the leading orbit-averaged dissipative backreaction, and examines the spin-curvature force to leading order in the small body's spin. This is sufficient to build generic inspirals of spinning bodies, and serves as a foundation for further work examining how to include secondary spin in large-mass-ratio waveform models. △ Less

Submitted 19 October, 2023; v1 submitted 15 May, 2023; originally announced May 2023.

Comments: 23 pages, 8 figures. Note that this posting is superseded by arXiv:2310.08438, in response to helpful critical referee feedback

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

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

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

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

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

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

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

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

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

Submitted 3 May, 2023; originally announced May 2023.

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

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

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

Submitted 1 May, 2023; originally announced May 2023.

Comments: 16 pages, 9 figures, Accepted in MNRAS

Report number: TIFR/TH/23-5

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

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

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

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

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

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

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

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

Submitted 5 March, 2023; originally announced March 2023.

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

ASCENT - A balloon-borne hard X-ray imaging spectroscopy telescope using transition edge sensor microcalorimeter detectors

Authors: Fabian Kislat, Daniel Becker, Douglas Bennett, Adrika Dasgupta, Joseph Fowler, Christopher L. Fryer, Johnathon Gard, Ephraim Gau, Danielle Gurgew, Keon Harmon, Takayuki Hayashi, Scott Heatwole, Md Arman Hossen, Henric Krawczynski, R. James Lanzi, Jason Legere, John A. B. Mates, Mark McConnell, Johanna Nagy, Takashi Okajima, Toshiki Sato, Daniel Schmidt, Sean Spooner, Daniel Swetz, Keisuke Tamura , et al. (4 additional authors not shown)

Abstract: Core collapse supernovae are thought to be one of the main sources in the galaxy of elements heavier than iron. Understanding the origin of the elements is thus tightly linked to our understanding of the explosion mechanism of supernovae and supernova nucleosynthesis. X-ray and gamma-ray observations of young supernova remnants, combined with improved theoretical modeling, have resulted in enormou… ▽ More Core collapse supernovae are thought to be one of the main sources in the galaxy of elements heavier than iron. Understanding the origin of the elements is thus tightly linked to our understanding of the explosion mechanism of supernovae and supernova nucleosynthesis. X-ray and gamma-ray observations of young supernova remnants, combined with improved theoretical modeling, have resulted in enormous improvements in our knowledge of these events. The isotope ${}^{44}$Ti is one of the most sensitive probes of the innermost regions of the core collapse engine, and its spatial and velocity distribution are key observables. Hard X-ray imaging spectroscopy with the Nuclear Spectroscopic Telescope Array (NuSTAR) has provided new insights into the structure of the supernova remnant Cassiopeia A (Cas A), establishing the convective nature of the supernova engine. However, many questions about the details of this engine remain. We present here the concept for a balloon-borne follow-up mission called ASCENT (A SuperConducting ENergetic x-ray Telescope). ASCENT uses transition edge sensor gamma-ray microcalorimeter detectors with a demonstrated 55 eV Full Width Half Maximum (FWHM) energy resolution at 97 keV. This 8--16-fold improvement in energy resolution over NuSTAR will allow high resolution imaging and spectroscopy of the ${}^{44}$Ti emission. This will allow a detailed reconstruction of gamma-ray line redshifts, widths, and shapes, allowing us to address questions such as: What is the source of the neutron star "kicks"? What is the dominant production pathway for ${}^{44}$Ti? Is the engine of Cas A unique? △ Less

Submitted 4 January, 2023; originally announced January 2023.

Comments: 23 pages, 11 figures

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

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

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

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

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

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

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

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

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

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

The 511-CAM Mission: A Pointed 511 keV Gamma-Ray Telescope with a Focal Plane Detector Made of Stacked Transition Edge Sensor Microcalorimeter Arrays

Authors: Farzane Shirazi, Md. Arman Hossen, Daniel Becker, Daniel Schmidt, Daniel Swetz, Douglas Bennett, Dana Braun, Ephraim Gau, Fabian Kislat, Johnathon Gard, John Mates, Joel Weber, Nicole Rodriguez Cavero, Sohee Chun, Lindsey Lisalda, Andrew West, Bhupal Dev, Francesc Ferrer, Richard Bose, Joel Ullom, Henric Krawczynski

Abstract: The 511 keV gamma-ray emission from the galactic center region may fully or partially originate from the annihilation of positrons from dark matter particles with electrons from the interstellar medium. Alternatively, the positrons could be created by astrophysical sources, involving exclusively standard model physics. We describe here a new concept for a 511 keV mission called 511-CAM (511 keV ga… ▽ More The 511 keV gamma-ray emission from the galactic center region may fully or partially originate from the annihilation of positrons from dark matter particles with electrons from the interstellar medium. Alternatively, the positrons could be created by astrophysical sources, involving exclusively standard model physics. We describe here a new concept for a 511 keV mission called 511-CAM (511 keV gamma-ray CAmera using Micro-calorimeters) that combines focusing gamma-ray optics with a stack of Transition Edge Sensor (TES) microcalorimeter arrays in the focal plane. The 511-CAM detector assembly has a projected 511 keV energy resolution of 390 eV Full Width Half Maximum (FWHM) or better, and improves by a factor of at least 11 on the performance of state-of-the-art Ge-based Compton telescopes. Combining this unprecedented energy resolution with sub-arcmin angular resolutions afforded by Laue lens or channeling optics could make substantial contributions to identifying the origin of the 511 keV emission by discovering and characterizing point sources and measuring line-of-sight velocities of the emitting plasmas. △ Less

Submitted 12 September, 2023; v1 submitted 29 June, 2022; originally announced June 2022.

Comments: 21 pages, 11 figures, 2 tables; revised manuscript that addressed the referee's comments

Journal ref: Journal of Astronomical Telescopes, Instruments, and Systems, Volume 9, id. 024006 (2023)

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

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

Abstract: We present a new investigation of the intergalactic medium (IGM) near reionization using dark gaps in the Lyman-$β$ (Ly$β$) forest. With its lower optical depth, Ly$β$ offers a potentially more sensitive probe to any remaining neutral gas compared to commonly used Ly$α$ line. We identify dark gaps in the Ly$β$ forest using spectra of 42 QSOs at $z_{\rm em}>5.5$, including new data from the XQR-30… ▽ More We present a new investigation of the intergalactic medium (IGM) near reionization using dark gaps in the Lyman-$β$ (Ly$β$) forest. With its lower optical depth, Ly$β$ offers a potentially more sensitive probe to any remaining neutral gas compared to commonly used Ly$α$ line. We identify dark gaps in the Ly$β$ forest using spectra of 42 QSOs at $z_{\rm em}>5.5$, including new data from the XQR-30 VLT Large Programme. Approximately $40\%$ of these QSO spectra exhibit dark gaps longer than $10h^{-1}{\rm Mpc}$ at $z\simeq5.8$. By comparing the results to predictions from simulations, we find that the data are broadly consistent both with models where fluctuations in the Ly$α$ forest are caused solely by ionizing ultraviolet background (UVB) fluctuations and with models that include large neutral hydrogen patches at $z<6$ due to a late end to reionization. Of particular interest is a very long ($L=28h^{-1}{\rm Mpc}$) and dark ($τ_{\rm eff} \gtrsim 6$) gap persisting down to $z\simeq 5.5$ in the Ly$β$ forest of the $z_{\rm}=5.85$ QSO PSO J025$-$11. This gap may support late reionization models with a volume-weighted average neutral hydrogen fraction of $ \langle x_{\rm HI}\rangle \gtrsim 5\%$ by $z=5.6$. Finally, we infer constraints on $\langle x_{\rm HI}\rangle$ over $5.5 \lesssim z \lesssim 6.0$ based on the observed Ly$β$ dark gap length distribution and a conservative relationship between gap length and neutral fraction derived from simulations. We find $\langle x_{\rm HI}\rangle \leq 0.05$, 0.17, and 0.29 at $z\simeq 5.55$, 5.75, and 5.95, respectively. These constraints are consistent with models where reionization ends significantly later than $z = 6$. △ Less

Submitted 9 May, 2022; originally announced May 2022.

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

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

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

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

Submitted 26 October, 2021; originally announced October 2021.

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

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

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

Abstract: The observed large-scale scatter in Lyman $α$ opacity of the intergalactic medium at $z<6$ implies large fluctuations in the neutral hydrogen fraction that are unexpected long after reionization has ended. A number of models have emerged to explain these fluctuations that make testable predictions for the relationship between Ly$α$ opacity and density. We present selections of $z=5.7$ Ly$α$-emitti… ▽ More The observed large-scale scatter in Lyman $α$ opacity of the intergalactic medium at $z<6$ implies large fluctuations in the neutral hydrogen fraction that are unexpected long after reionization has ended. A number of models have emerged to explain these fluctuations that make testable predictions for the relationship between Ly$α$ opacity and density. We present selections of $z=5.7$ Ly$α$-emitting galaxies (LAEs) in the fields surrounding two highly opaque quasar sightlines with long Ly$α$ troughs. The fields lie towards the $z=6.0$ quasar ULAS J0148+0600, for which we re-analyze previously published results using improved photometric selection, and towards the $z=6.15$ quasar SDSS J1250+3130, for which results are presented here for the first time. In both fields, we report a deficit of LAEs within 20 $h^{-1}$ Mpc of the quasar. The association of highly opaque sightlines with galaxy underdensities in these two fields is consistent with models in which the scatter in Ly$α$ opacity is driven by large-scale fluctuations in the ionizing UV background, or by an ultra-late reionization that has not yet concluded at $z=5.7$. △ Less

Submitted 27 September, 2021; originally announced September 2021.

Comments: Accepted to ApJ. 17 pages, 13 figures

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

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

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

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

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

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

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

Abstract: We present a new investigation of the intergalactic medium (IGM) near the end of reionization using "dark gaps" in the Lyman-alpha (Ly$α$) forest. Using spectra of 55 QSOs at $z_{\rm em}>5.5$, including new data from the XQR-30 VLT Large Programme, we identify gaps in the Ly$α$ forest where the transmission averaged over 1 comoving $h^{-1}\,{\rm Mpc}$ bins falls below 5%. Nine ultra-long (… ▽ More We present a new investigation of the intergalactic medium (IGM) near the end of reionization using "dark gaps" in the Lyman-alpha (Ly$α$) forest. Using spectra of 55 QSOs at $z_{\rm em}>5.5$, including new data from the XQR-30 VLT Large Programme, we identify gaps in the Ly$α$ forest where the transmission averaged over 1 comoving $h^{-1}\,{\rm Mpc}$ bins falls below 5%. Nine ultra-long ($L > 80~h^{-1}\,{\rm Mpc}$) dark gaps are identified at $z<6$. In addition, we quantify the fraction of QSO spectra exhibiting gaps longer than $30~h^{-1}\,{\rm Mpc}$, $F_{30}$, as a function of redshift. We measure $F_{30} \simeq 0.9$, 0.6, and 0.15 at $z = 6.0$, 5.8, and 5.6, respectively, with the last of these long dark gaps persisting down to $z \simeq 5.3$. Comparing our results with predictions from hydrodynamical simulations, we find that the data are consistent with models wherein reionization extends significantly below redshift six. Models wherein the IGM is essentially fully reionized that retain large-scale fluctuations in the ionizing UV background at $z \lesssim 6$ are also potentially consistent with the data. Overall, our results suggest that signature of reionization in the form of islands of neutral hydrogen and/or large-scale fluctuations in the ionizing background remain present in the IGM until at least $z \simeq 5.3$. △ Less

Submitted 13 September, 2021; originally announced September 2021.

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

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

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

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

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

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

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

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

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

Submitted 21 May, 2021; originally announced May 2021.

Comments: 8 pages, 5 figures, submitted to ApJL

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

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

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

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

Comments: 8 pages, 3 figures

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

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

Abstract: The mean free path of ionizing photons, $λ_{\rm mfp}$, is a key factor in the photoionization of the intergalactic medium (IGM). At $z \gtrsim 5$, however, $λ_{\rm mfp}$ may be short enough that measurements towards QSOs are biased by the QSO proximity effect. We present new direct measurements of $λ_{\rm mfp}$ that address this bias and extend up to $z \sim 6$ for the first time. Our measurements… ▽ More The mean free path of ionizing photons, $λ_{\rm mfp}$, is a key factor in the photoionization of the intergalactic medium (IGM). At $z \gtrsim 5$, however, $λ_{\rm mfp}$ may be short enough that measurements towards QSOs are biased by the QSO proximity effect. We present new direct measurements of $λ_{\rm mfp}$ that address this bias and extend up to $z \sim 6$ for the first time. Our measurements at $z \sim 5$ are based on data from the Giant Gemini GMOS survey and new Keck LRIS observations of low-luminosity QSOs. At $z \sim 6$ we use QSO spectra from Keck ESI and VLT X-Shooter. We measure $λ_{\rm mfp} = 9.09^{+1.62}_{-1.28}$ proper Mpc and $0.75^{+0.65}_{-0.45}$ proper Mpc (68% confidence) at $z = 5.1$ and 6.0, respectively. The results at $z = 5.1$ are consistent with existing measurements, suggesting that bias from the proximity effect is minor at this redshift. At $z = 6.0$, however, we find that neglecting the proximity effect biases the result high by a factor of two or more. Our measurement at $z = 6.0$ falls well below extrapolations from lower redshifts, indicating rapid evolution in $λ_{\rm mfp}$ over $5 < z < 6$. This evolution disfavors models in which reionization ended early enough that the IGM had time to fully relax hydrodynamically by $z = 6$, but is qualitatively consistent with models wherein reionization completed at $z = 6$ or even significantly later. Our mean free path results are most consistent with late reionization models wherein the IGM is still 20% neutral at $z=6$, although our measurement at $z = 6.0$ is even lower than these models prefer. △ Less

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

Comments: Accepted for publication in MNRAS

The Atacama Cosmology Telescope: Summary of DR4 and DR5 Data Products and Data Access

Authors: Maya Mallaby-Kay, Zachary Atkins, Simone Aiola, Stefania Amodeo, Jason E. Austermann, James A. Beall, Daniel T. Becker, J. Richard Bond, Erminia Calabrese, Grace E. Chesmore, Steve K. Choi, Kevin T. Crowley, Omar Darwish, Edwawd V. Denison, Mark J. Devlin, Shannon M. Duff, Adriaan J. Duivenvoorden, Jo Dunkley, Simone Ferraro, Kyra Fichman, Patricio A. Gallardo, Joseph E. Golec, Yilun Guan, Dongwon Han, Matthew Hasselfield , et al. (35 additional authors not shown)

Abstract: Two recent large data releases for the Atacama Cosmology Telescope (ACT), called DR4 and DR5, are available for public access. These data include temperature and polarization maps that cover nearly half the sky at arcminute resolution in three frequency bands; lensing maps and component-separated maps covering ~ 2,100 deg^2 of sky; derived power spectra and cosmological likelihoods; a catalog of o… ▽ More Two recent large data releases for the Atacama Cosmology Telescope (ACT), called DR4 and DR5, are available for public access. These data include temperature and polarization maps that cover nearly half the sky at arcminute resolution in three frequency bands; lensing maps and component-separated maps covering ~ 2,100 deg^2 of sky; derived power spectra and cosmological likelihoods; a catalog of over 4,000 galaxy clusters; and supporting ancillary products including beam functions and masks. The data and products are described in a suite of ACT papers; here we provide a summary. In order to facilitate ease of access to these data we present a set of Jupyter IPython notebooks developed to introduce users to DR4, DR5, and the tools needed to analyze these data. The data products (excluding simulations) and the set of notebooks are publicly available on the NASA Legacy Archive for Microwave Background Data Analysis (LAMBDA); simulation products are available on the National Energy Research Scientific Computing Center (NERSC). △ Less

Submitted 29 April, 2021; v1 submitted 4 March, 2021; originally announced March 2021.

Comments: Accepted to ApJS. 21 pages, 8 figures. Data and notebooks available on LAMBDA https://lambda.gsfc.nasa.gov/product/act/

The design of the Ali CMB Polarization Telescope receiver

Authors: Maria Salatino, Jason E. Austermann, Keith L. Thompson, Peter A. R. Ade, Xiran Bai, James A. Beall, Dan T. Becker, Yifu Cai, Zhi Chang, Ding Chen, Pisin Chen, Jake Connors, Jacques Delabrouille, Bradley Dober, Shannon M. Duff, Guanhua Gao, Shamik Ghosh, Richard C. Givhan, Gene C. Hilton, Bin Hu, Johannes Hubmayr, Ethan D. Karpel, Chao-Lin Kuo, Hong Li, Mingzhe Li , et al. (50 additional authors not shown)

Abstract: Ali CMB Polarization Telescope (AliCPT-1) is the first CMB degree-scale polarimeter to be deployed on the Tibetan plateau at 5,250m above sea level. AliCPT-1 is a 90/150 GHz 72 cm aperture, two-lens refracting telescope cooled down to 4 K. Alumina lenses, 800mm in diameter, image the CMB in a 33.4° field of view on a 636mm wide focal plane. The modularized focal plane consists of dichroic polariza… ▽ More Ali CMB Polarization Telescope (AliCPT-1) is the first CMB degree-scale polarimeter to be deployed on the Tibetan plateau at 5,250m above sea level. AliCPT-1 is a 90/150 GHz 72 cm aperture, two-lens refracting telescope cooled down to 4 K. Alumina lenses, 800mm in diameter, image the CMB in a 33.4° field of view on a 636mm wide focal plane. The modularized focal plane consists of dichroic polarization-sensitive Transition-Edge Sensors (TESes). Each module includes 1,704 optically active TESes fabricated on a 150mm diameter silicon wafer. Each TES array is read out with a microwave multiplexing readout system capable of a multiplexing factor up to 2,048. Such a large multiplexing factor has allowed the practical deployment of tens of thousands of detectors, enabling the design of a receiver that can operate up to 19 TES arrays for a total of 32,376 TESes. AliCPT-1 leverages the technological advancements in the detector design from multiple generations of previously successful feedhorn-coupled polarimeters, and in the instrument design from BICEP-3, but applied on a larger scale. The cryostat receiver is currently under integration and testing. During the first deployment year, the focal plane will be populated with up to 4 TES arrays. Further TES arrays will be deployed in the following years, fully populating the focal plane with 19 arrays on the fourth deployment year. Here we present the AliCPT-1 receiver design, and how the design has been optimized to meet the experimental requirements. △ Less

Submitted 23 January, 2021; originally announced January 2021.

Comments: Proc. SPIE, 11453, 114532A (2020)

Journal ref: Proceedings of SPIE, 'Millimeter, Submillimeter, and Far-Infrared Detectors and Instrumentation for Astronomy X', Volume 11453, 114532A (2020)

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

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

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

Submitted 19 January, 2021; originally announced January 2021.

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

Design and pre-flight performance of SPIDER 280 GHz receivers

Authors: E. C. Shaw, P. A. R. Ade, S. Akers, M. Amiri, J. Austermann, J. Beall, D. T. Becker, S. J. Benton, A. S. Bergman, J. J. Bock, J. R. Bond, S. A. Bryan, H. C. Chiang, C. R. Contaldi, R. S. Domagalski, O. Doré, S. M. Duff, A. J. Duivenvoorden, H. K. Eriksen, M. Farhang, J. P. Filippini, L. M. Fissel, A. A. Fraisse, K. Freese, M. Galloway , et al. (57 additional authors not shown)

Abstract: In this work we describe upgrades to the Spider balloon-borne telescope in preparation for its second flight, currently planned for December 2021. The Spider instrument is optimized to search for a primordial B-mode polarization signature in the cosmic microwave background at degree angular scales. During its first flight in 2015, Spider mapped ~10% of the sky at 95 and 150 GHz. The payload for th… ▽ More In this work we describe upgrades to the Spider balloon-borne telescope in preparation for its second flight, currently planned for December 2021. The Spider instrument is optimized to search for a primordial B-mode polarization signature in the cosmic microwave background at degree angular scales. During its first flight in 2015, Spider mapped ~10% of the sky at 95 and 150 GHz. The payload for the second Antarctic flight will incorporate three new 280 GHz receivers alongside three refurbished 95- and 150 GHz receivers from Spider's first flight. In this work we discuss the design and characterization of these new receivers, which employ over 1500 feedhorn-coupled transition-edge sensors. We describe pre-flight laboratory measurements of detector properties, and the optical performance of completed receivers. These receivers will map a wide area of the sky at 280 GHz, providing new information on polarized Galactic dust emission that will help to separate it from the cosmological signal. △ Less

Submitted 22 December, 2020; originally announced December 2020.

Comments: 13 pages, 8 figures; as published in the conference proceedings for SPIE Millimeter, Submillimeter, and Far-Infrared Detectors and Instrumentation for Astronomy X (2020)

Journal ref: "Design and pre-flight performance of SPIDER 280 GHz receivers," Proc. SPIE 11453, Millimeter, Submillimeter, and Far-Infrared Detectors and Instrumentation for Astronomy X, 114532F (13 December 2020)

A Microwave SQUID Multiplexer Optimized for Bolometric Applications

Authors: B. Dober, Z. Ahmed, K. Arnold, D. T. Becker, D. A. Bennett, J. A. Connors, A. Cukierman, J. M. D'Ewart, S. M. Duff, J. E. Dusatko, J. C. Frisch, J. D. Gard, S. W. Henderson, R. Herbst, G. C. Hilton, J. Hubmayr, Y. Li, J. A. B. Mates, H. McCarrick, C. D Reintsema, M. Silva-Feaver, L. Ruckman, J. N. Ullom, L. R. Vale, D. D. Van Winkle , et al. (5 additional authors not shown)

Abstract: A microwave SQUID multiplexer ($μ$MUX) has been optimized for coupling to large arrays of superconducting transition-edge sensor (TES) bolometers. We present the scalable cryogenic multiplexer chip design in a 1820-channel multiplexer configuration for the 4-8 GHz rf band. The key metrics of yield, sensitivity, and crosstalk are determined through measurements of 455 readout channels, which span 4… ▽ More A microwave SQUID multiplexer ($μ$MUX) has been optimized for coupling to large arrays of superconducting transition-edge sensor (TES) bolometers. We present the scalable cryogenic multiplexer chip design in a 1820-channel multiplexer configuration for the 4-8 GHz rf band. The key metrics of yield, sensitivity, and crosstalk are determined through measurements of 455 readout channels, which span 4-5 GHz. The median white-noise level is 45 pA/$\sqrt{\textrm{Hz}}$, evaluated at 2 Hz, with a 1/f knee $\leq$ 20 mHz after common-mode subtraction. The white-noise level decreases the sensitivity of a TES bolometer optimized for detection of the cosmic microwave background at 150 GHz by only 3%. The measured crosstalk between any channel pair is $\leq$ 0.3%. △ Less

Submitted 19 January, 2021; v1 submitted 15 October, 2020; originally announced October 2020.

Comments: 6 pages, 5 figures

The Atacama Cosmology Telescope: A Catalog of > 4000 Sunyaev-Zel'dovich Galaxy Clusters

Authors: M. Hilton, C. Sifón, S. Naess, M. Madhavacheril, M. Oguri, E. Rozo, E. Rykoff, T. M. C. Abbott, S. Adhikari, M. Aguena, S. Aiola, S. Allam, S. Amodeo, A. Amon, J. Annis, B. Ansarinejad, C. Aros-Bunster, J. E. Austermann, S. Avila, D. Bacon, N. Battaglia, J. A. Beall, D. T. Becker, G. M. Bernstein, E. Bertin , et al. (124 additional authors not shown)

Abstract: We present a catalog of 4195 optically confirmed Sunyaev-Zel'dovich (SZ) selected galaxy clusters detected with signal-to-noise > 4 in 13,211 deg$^2$ of sky surveyed by the Atacama Cosmology Telescope (ACT). Cluster candidates were selected by applying a multi-frequency matched filter to 98 and 150 GHz maps constructed from ACT observations obtained from 2008-2018, and confirmed using deep, wide-a… ▽ More We present a catalog of 4195 optically confirmed Sunyaev-Zel'dovich (SZ) selected galaxy clusters detected with signal-to-noise > 4 in 13,211 deg$^2$ of sky surveyed by the Atacama Cosmology Telescope (ACT). Cluster candidates were selected by applying a multi-frequency matched filter to 98 and 150 GHz maps constructed from ACT observations obtained from 2008-2018, and confirmed using deep, wide-area optical surveys. The clusters span the redshift range 0.04 < z < 1.91 (median z = 0.52). The catalog contains 222 z > 1 clusters, and a total of 868 systems are new discoveries. Assuming an SZ-signal vs. mass scaling relation calibrated from X-ray observations, the sample has a 90% completeness mass limit of M500c > 3.8 x 10$^{14}$ MSun, evaluated at z = 0.5, for clusters detected at signal-to-noise ratio > 5 in maps filtered at an angular scale of 2.4'. The survey has a large overlap with deep optical weak-lensing surveys that are being used to calibrate the SZ-signal mass-scaling relation, such as the Dark Energy Survey (4566 deg$^2$), the Hyper Suprime-Cam Subaru Strategic Program (469 deg$^2$), and the Kilo Degree Survey (825 deg$^2$). We highlight some noteworthy objects in the sample, including potentially projected systems; clusters with strong lensing features; clusters with active central galaxies or star formation; and systems of multiple clusters that may be physically associated. The cluster catalog will be a useful resource for future cosmological analyses, and studying the evolution of the intracluster medium and galaxies in massive clusters over the past 10 Gyr. △ Less

Submitted 2 December, 2020; v1 submitted 23 September, 2020; originally announced September 2020.

Comments: 35 pages, 27 figures, accepted for publication in ApJS; v1.0 catalogs will be available from LAMBDA https://lambda.gsfc.nasa.gov/product/act/actpol_prod_table.cfm; v1.0 catalogs available from https://astro.ukzn.ac.za/~mjh/ACTDR5/v1.0/ until then

The Atacama Cosmology Telescope: Weighing distant clusters with the most ancient light

Authors: Mathew S. Madhavacheril, Cristóbal Sifón, Nicholas Battaglia, Simone Aiola, Stefania Amodeo, Jason E. Austermann, James A. Beall, Daniel T. Becker, J. Richard Bond, Erminia Calabrese, Steve K. Choi, Edward V. Denison, Mark J. Devlin, Simon R. Dicker, Shannon M. Duff, Adriaan J. Duivenvoorden, Jo Dunkley, Rolando Dünner, Simone Ferraro, Patricio A. Gallardo, Yilun Guan, Dongwon Han, J. Colin Hill, Gene C. Hilton, Matt Hilton , et al. (36 additional authors not shown)

Abstract: We use gravitational lensing of the cosmic microwave background (CMB) to measure the mass of the most distant blindly-selected sample of galaxy clusters on which a lensing measurement has been performed to date. In CMB data from the the Atacama Cosmology Telescope (ACT) and the Planck satellite, we detect the stacked lensing effect from 677 near-infrared-selected galaxy clusters from the Massive a… ▽ More We use gravitational lensing of the cosmic microwave background (CMB) to measure the mass of the most distant blindly-selected sample of galaxy clusters on which a lensing measurement has been performed to date. In CMB data from the the Atacama Cosmology Telescope (ACT) and the Planck satellite, we detect the stacked lensing effect from 677 near-infrared-selected galaxy clusters from the Massive and Distant Clusters of WISE Survey (MaDCoWS), which have a mean redshift of $ \langle z \rangle = 1.08$. There are no current optical weak lensing measurements of clusters that match the distance and average mass of this sample. We detect the lensing signal with a significance of $4.2 σ$. We model the signal with a halo model framework to find the mean mass of the population from which these clusters are drawn. Assuming that the clusters follow Navarro-Frenk-White density profiles, we infer a mean mass of $\langle M_{500c}\rangle = \left(1.7 \pm 0.4 \right)\times10^{14}\,\mathrm{M}_\odot$. We consider systematic uncertainties from cluster redshift errors, centering errors, and the shape of the NFW profile. These are all smaller than 30% of our reported uncertainty. This work highlights the potential of CMB lensing to enable cosmological constraints from the abundance of distant clusters populating ever larger volumes of the observable Universe, beyond the capabilities of optical weak lensing measurements. △ Less

Submitted 1 November, 2020; v1 submitted 16 September, 2020; originally announced September 2020.

Comments: 14 pages, 3 figures, matches version accepted in ApJL, code available at https://github.com/ACTCollaboration/madcows_lensing/

The Atacama Cosmology Telescope: Modeling the Gas Thermodynamics in BOSS CMASS galaxies from Kinematic and Thermal Sunyaev-Zel'dovich Measurements

Authors: Stefania Amodeo, Nicholas Battaglia, Emmanuel Schaan, Simone Ferraro, Emily Moser, Simone Aiola, Jason E. Austermann, James A. Beall, Rachel Bean, Daniel T. Becker, Richard J. Bond, Erminia Calabrese, Victoria Calafut, Steve K. Choi, Edward V. Denison, Mark Devlin, Shannon M. Duff, Adriaan J. Duivenvoorden, Jo Dunkley, Rolando Dünner, Patricio A. Gallardo, Kirsten R. Hall, Dongwon Han, J. Colin Hill, Gene C. Hilton , et al. (30 additional authors not shown)

Abstract: The thermal and kinematic Sunyaev-Zel'dovich effects (tSZ, kSZ) probe the thermodynamic properties of the circumgalactic and intracluster medium (CGM and ICM) of galaxies, groups, and clusters, since they are proportional, respectively, to the integrated electron pressure and momentum along the line-of-sight. We present constraints on the gas thermodynamics of CMASS galaxies in the Baryon Oscillat… ▽ More The thermal and kinematic Sunyaev-Zel'dovich effects (tSZ, kSZ) probe the thermodynamic properties of the circumgalactic and intracluster medium (CGM and ICM) of galaxies, groups, and clusters, since they are proportional, respectively, to the integrated electron pressure and momentum along the line-of-sight. We present constraints on the gas thermodynamics of CMASS galaxies in the Baryon Oscillation Spectroscopic Survey (BOSS) using new measurements of the kSZ and tSZ signals obtained in a companion paper. Combining kSZ and tSZ measurements, we measure within our model the amplitude of energy injection $εM_\star c^2$, where $M_\star$ is the stellar mass, to be $ε=(40\pm9)\times10^{-6}$, and the amplitude of the non-thermal pressure profile to be $α_{\rm Nth}<0.2$ (2$σ$), indicating that less than 20% of the total pressure within the virial radius is due to a non-thermal component. We estimate the effects of including baryons in the modeling of weak-lensing galaxy cross-correlation measurements using the best-fit density profile from the kSZ measurement. Our estimate reduces the difference between the original theoretical model and the weak-lensing galaxy cross-correlation measurements in arXiv:1611.08606 by half but does not fully reconcile it. Comparing the tSZ measurements to cosmological simulations, we find that simulations underestimate the CGM pressure at large radii while they fare better in comparison with the kSZ measurements. This suggests that the energy injected via feedback models in the simulations that we compared against does not sufficiently heat the gas at these radii. We do not find significant disagreement at smaller radii. These measurements provide novel tests of current and future simulations. This work demonstrates the power of joint, high signal-to-noise kSZ and tSZ observations, upon which future cross-correlation studies will improve. △ Less

Submitted 9 February, 2023; v1 submitted 11 September, 2020; originally announced September 2020.

Comments: Corrected error in the algorithm that calculates the kSZ temperature profile for a given GNFW density model. The value of $\log_{\rm10} ρ_0$ changed by 0.75$σ$ in Tab.II and Fig.2, and affected the results in the left panels of Fig. 6. Conclusions are unchanged. Erratum published at https://link.aps.org/doi/10.1103/PhysRevD.107.049903

Journal ref: Phys. Rev. D 103, 063514 (2021)

The Atacama Cosmology Telescope: Combined kinematic and thermal Sunyaev-Zel'dovich measurements from BOSS CMASS and LOWZ halos

Authors: Emmanuel Schaan, Simone Ferraro, Stefania Amodeo, Nick Battaglia, Simone Aiola, Jason E. Austermann, James A. Beall, Rachel Bean, Daniel T. Becker, Richard J. Bond, Erminia Calabrese, Victoria Calafut, Steve K. Choi, Edward V. Denison, Mark J. Devlin, Shannon M. Duff, Adriaan J. Duivenvoorden, Jo Dunkley, Rolando Dünner, Patricio A. Gallardo, Yilun Guan, Dongwon Han, J. Colin Hill, Gene C. Hilton, Matt Hilton , et al. (33 additional authors not shown)

Abstract: The scattering of cosmic microwave background (CMB) photons off the free-electron gas in galaxies and clusters leaves detectable imprints on high resolution CMB maps: the thermal and kinematic Sunyaev-Zel'dovich effects (tSZ and kSZ respectively). We use combined microwave maps from the Atacama Cosmology Telescope (ACT) DR5 and Planck in combination with the CMASS and LOWZ galaxy catalogs from the… ▽ More The scattering of cosmic microwave background (CMB) photons off the free-electron gas in galaxies and clusters leaves detectable imprints on high resolution CMB maps: the thermal and kinematic Sunyaev-Zel'dovich effects (tSZ and kSZ respectively). We use combined microwave maps from the Atacama Cosmology Telescope (ACT) DR5 and Planck in combination with the CMASS and LOWZ galaxy catalogs from the Baryon Oscillation Spectroscopic Survey (BOSS DR10 and DR12), to study the gas associated with these galaxy groups. Using individual reconstructed velocities, we perform a stacking analysis and reject the no-kSZ hypothesis at 6.5$σ$, the highest significance to date. This directly translates into a measurement of the electron number density profile, and thus of the gas density profile. Despite the limited signal to noise, the measurement shows at high significance that the gas density profile is more extended than the dark matter density profile, for any reasonable baryon abundance (formally $>90σ$ for the cosmic baryon abundance). We simultaneously measure the tSZ signal, i.e. the electron thermal pressure profile of the same CMASS objects, and reject the no-tSZ hypothesis at 10$σ$. We combine tSZ and kSZ measurements to estimate the electron temperature to 20% precision in several aperture bins, and find it comparable to the virial temperature. In a companion paper, we analyze these measurements to constrain the gas thermodynamics and the properties of feedback inside galaxy groups. We present the corresponding LOWZ measurements in this paper, ruling out a null kSZ (tSZ) signal at 2.9 (13.9)$σ$, and leave their interpretation to future work. Our stacking software ThumbStack is publicly available at https://github.com/EmmanuelSchaan/ThumbStack and directly applicable to future Simons Observatory and CMB-S4 data. △ Less

Submitted 16 February, 2021; v1 submitted 11 September, 2020; originally announced September 2020.

Comments: Accepted in Physical Review D, Editors' Suggestion

Journal ref: Phys. Rev. D 103, 063513 (2021)

Probing galaxy evolution in massive clusters using ACT and DES: splashback as a cosmic clock

Authors: Susmita Adhikari, Tae-hyeon Shin, Bhuvnesh Jain, Matt Hilton, Eric Baxter, Chihway Chang, Risa H. Wechsler, Nick Battaglia, J. Richard Bond, Sebastian Bocquet, Joseph DeRose, Steve K. Choi, Mark Devlin, Jo Dunkley, August E. Evrard, Simone Ferraro, J. Colin Hill, John P. Hughes, Patricio A. Gallardo, Martine Lokken, Amanda MacInnis, Jeffrey McMahon, Mathew S. Madhavacheril, Frederico Nati, Laura B. Newburgh , et al. (91 additional authors not shown)

Abstract: We measure the projected number density profiles of galaxies and the splashback feature in clusters selected by the Sunyaev--Zeldovich (SZ) effect from the Advanced Atacama Cosmology Telescope (AdvACT) survey using galaxies observed by the Dark Energy Survey (DES). The splashback radius for the complete galaxy sample is consistent with theoretical measurements from CDM-only simulations, and is loc… ▽ More We measure the projected number density profiles of galaxies and the splashback feature in clusters selected by the Sunyaev--Zeldovich (SZ) effect from the Advanced Atacama Cosmology Telescope (AdvACT) survey using galaxies observed by the Dark Energy Survey (DES). The splashback radius for the complete galaxy sample is consistent with theoretical measurements from CDM-only simulations, and is located at $2.4^{+0.3}_{-0.4}$ Mpc $h^{-1}$. We split the sample based on galaxy color and find significant differences in the profile shapes. Red galaxies and those in the green valley show a splashback-like minimum in their slope profile consistent with theoretical predictions, while the bluest galaxies show a weak feature that appears at a smaller radius. We develop a mapping of galaxies to subhalos in $N$-body simulations by splitting subhalos based on infall time onto the cluster halos. We find that the location of the steepest slope and differences in the shapes of the profiles can be mapped to differences in the average time of infall of galaxies of different colors. The minima of the slope in the galaxy profiles trace a discontinuity in the phase space of dark matter halos. By relating spatial profiles to infall time for galaxies of different colours, we can use splashback as a clock to understand galaxy quenching. We find that red galaxies have on average been in their clusters for over $3.2 ~\rm Gyrs$, green galaxies about $2.2 ~\rm Gyrs$, while blue galaxies have been accreted most recently and have not reached apocenter. Using the information from the complete radial profiles, we fit a simple quenching model and find that the onset of galaxy quenching in clusters occurs after a delay of about a gigayear, and that galaxies quench rapidly thereafter with an exponential timescale of $0.6$ Gyr. △ Less

Submitted 26 August, 2020; originally announced August 2020.

Comments: 23 pages, 10 figures, to be submitted to ApJ

The Atacama Cosmology Telescope: Delensed Power Spectra and Parameters

Authors: Dongwon Han, Neelima Sehgal, Amanda MacInnis, Alexander van Engelen, Blake D. Sherwin, Mathew S. Madhavacheril, Simone Aiola, Nicholas Battaglia, James A. Beall, Daniel T. Becker, Erminia Calabrese, Steve K. Choi, Omar Darwish, Edward V. Denison, Mark J. Devlin, Jo Dunkley, Simone Ferraro, Anna E. Fox, Matthew Hasselfield, J. Colin Hill, Gene C. Hilton, Matt Hilton, Renée Hložek, Johannes Hubmayr, John P. Hughes , et al. (17 additional authors not shown)

Abstract: We present LCDM cosmological parameter constraints obtained from delensed microwave background power spectra. Lensing maps from a subset of DR4 data from the Atacama Cosmology Telescope (ACT) are used to undo the lensing effect in ACT spectra observed at 150 and 98 GHz. At 150 GHz, we remove the lensing distortion with an effective efficiency of 30% (TT), 30% (EE), 26% (TE) and 20% (BB); this resu… ▽ More We present LCDM cosmological parameter constraints obtained from delensed microwave background power spectra. Lensing maps from a subset of DR4 data from the Atacama Cosmology Telescope (ACT) are used to undo the lensing effect in ACT spectra observed at 150 and 98 GHz. At 150 GHz, we remove the lensing distortion with an effective efficiency of 30% (TT), 30% (EE), 26% (TE) and 20% (BB); this results in detections of the delensing effect at 8.7 sigma (TT), 5.1 sigma (EE), 2.6 sigma (TE), and 2.4 sigma (BB) significance. The combination of 150 and 98 GHz TT, EE, and TE delensed spectra is well fit by a standard LCDM model. We also measure the shift in best-fit parameters when fitting delensed versus lensed spectra; while this shift does not inform our ability to measure cosmological parameters, it does provide a three-way consistency check among the lensing inferred from the best-fit parameters, the lensing in the CMB power spectrum, and the reconstructed lensing map. This shift is predicted to be zero when fitting with the correct model since both lensed and delensed spectra originate from the same region of sky. Fitting with a LCDM model and marginalizing over foregrounds, we find that the shift in cosmological parameters is consistent with zero. Our results show that gravitational lensing of the microwave background is internally consistent within the framework of the standard cosmological model. △ Less

Submitted 13 November, 2020; v1 submitted 28 July, 2020; originally announced July 2020.

Comments: 29 pages, 17 figures, version matches that accepted by JCAP

Journal ref: JCAP, Issue 01, article id. 031 (2021)

The Atacama Cosmology Telescope: DR5 maps of 18,000 square degrees of the microwave sky from ACT 2008-2018 data

Authors: Sigurd Naess, Simone Aiola, Jason E. Austermann, Nick Battaglia, James A. Beall, Daniel T. Becker, Richard J. Bond, Erminia Calabrese, Steve K. Choi, Nicholas F. Cothard, Kevin T. Crowley, Omar Darwish, Rahul Datta, Edward V. Denison, Mark Devlin, Cody J. Duell, Shannon M. Duff, Adriaan J. Duivenvoorden, Jo Dunkley, Rolando Dünner, Anna E. Fox, Patricio A. Gallardo, Mark Halpern, Dongwon Han, Matthew Hasselfield , et al. (37 additional authors not shown)

Abstract: This paper presents a maximum-likelihood algorithm for combining sky maps with disparate sky coverage, angular resolution and spatially varying anisotropic noise into a single map of the sky. We use this to merge hundreds of individual maps covering the 2008-2018 ACT observing seasons, resulting in by far the deepest ACT maps released so far. We also combine the maps with the full Planck maps, res… ▽ More This paper presents a maximum-likelihood algorithm for combining sky maps with disparate sky coverage, angular resolution and spatially varying anisotropic noise into a single map of the sky. We use this to merge hundreds of individual maps covering the 2008-2018 ACT observing seasons, resulting in by far the deepest ACT maps released so far. We also combine the maps with the full Planck maps, resulting in maps that have the best features of both Planck and ACT: Planck's nearly white noise on intermediate and large angular scales and ACT's high-resolution and sensitivity on small angular scales. The maps cover over 18,000 square degrees, nearly half the full sky, at 100, 150 and 220 GHz. They reveal 4,000 optically-confirmed clusters through the Sunyaev Zel'dovich effect (SZ) and 18,500 point source candidates at $> 5σ$, the largest single collection of SZ clusters and millimeter wave sources to date. The multi-frequency maps provide millimeter images of nearby galaxies and individual Milky Way nebulae, and even clear detections of several nearby stars. Other anticipated uses of these maps include, for example, thermal SZ and kinematic SZ cluster stacking, CMB cluster lensing and galactic dust science. The method itself has negligible bias. However, due to the preliminary nature of some of the component data sets, we caution that these maps should not be used for precision cosmological analysis. The maps are part of ACT DR5, and are available on LAMBDA at https://lambda.gsfc.nasa.gov/product/act/actpol_prod_table.cfm. There is also a web atlas at https://phy-act1.princeton.edu/public/snaess/actpol/dr5/atlas. △ Less

Submitted 17 February, 2021; v1 submitted 14 July, 2020; originally announced July 2020.

Comments: 38 pages, 29 figures, data release on lambda. Published in JCAP

The Atacama Cosmology Telescope: A Measurement of the Cosmic Microwave Background Power Spectra at 98 and 150 GHz

Authors: Steve K. Choi, Matthew Hasselfield, Shuay-Pwu Patty Ho, Brian Koopman, Marius Lungu, Maximilian H. Abitbol, Graeme E. Addison, Peter A. R. Ade, Simone Aiola, David Alonso, Mandana Amiri, Stefania Amodeo, Elio Angile, Jason E. Austermann, Taylor Baildon, Nick Battaglia, James A. Beall, Rachel Bean, Daniel T. Becker, J Richard Bond, Sarah Marie Bruno, Erminia Calabrese, Victoria Calafut, Luis E. Campusano, Felipe Carrero , et al. (114 additional authors not shown)

Abstract: We present the temperature and polarization angular power spectra of the CMB measured by the Atacama Cosmology Telescope (ACT) from 5400 deg$^2$ of the 2013-2016 survey, which covers $>$15000 deg$^2$ at 98 and 150 GHz. For this analysis we adopt a blinding strategy to help avoid confirmation bias and, related to this, show numerous checks for systematic error done before unblinding. Using the like… ▽ More We present the temperature and polarization angular power spectra of the CMB measured by the Atacama Cosmology Telescope (ACT) from 5400 deg$^2$ of the 2013-2016 survey, which covers $>$15000 deg$^2$ at 98 and 150 GHz. For this analysis we adopt a blinding strategy to help avoid confirmation bias and, related to this, show numerous checks for systematic error done before unblinding. Using the likelihood for the cosmological analysis we constrain secondary sources of anisotropy and foreground emission, and derive a "CMB-only" spectrum that extends to $\ell=4000$. At large angular scales, foreground emission at 150 GHz is $\sim$1% of TT and EE within our selected regions and consistent with that found by Planck. Using the same likelihood, we obtain the cosmological parameters for $Λ$CDM for the ACT data alone with a prior on the optical depth of $τ=0.065\pm0.015$. $Λ$CDM is a good fit. The best-fit model has a reduced $χ^2$ of 1.07 (PTE=0.07) with $H_0=67.9\pm1.5$ km/s/Mpc. We show that the lensing BB signal is consistent with $Λ$CDM and limit the celestial EB polarization angle to $ψ_P =-0.07^{\circ}\pm0.09^{\circ}$. We directly cross correlate ACT with Planck and observe generally good agreement but with some discrepancies in TE. All data on which this analysis is based will be publicly released. △ Less

Submitted 23 November, 2020; v1 submitted 14 July, 2020; originally announced July 2020.

Comments: 44 pages, 27 figures, products available on the NASA LAMBDA website, version accepted for publication in JCAP