Search

The stochastic nature of star formation and photon propagation in high-redshift galaxies can result in sizable galaxy-to-galaxy scatter in their properties. Ignoring this scatter by assuming mean quantities can bias estimates of their emissivity and corresponding observables. We construct a flexible, semi-empirical model, sampling scatter around the following mean relations: (i) the conditional halo mass function (CHMF); (ii) the stellar-to-halo mass relation (SHMR); (iii) galaxy star formation main sequence (SFMS); (iv) fundamental metallicity relation (FMR); (v) conditional intrinsic luminosity; and (vi) photon escape fraction. In our fiducial model, ignoring scatter in these galaxy properties overestimates the duration of the EoR, delaying its completion by up to $\Delta z$ ~ 2. We quantify the relative importance of each of the above sources of scatter in determining the ionizing, soft-band X-ray and Lyman Werner (LW) emissivities as a function of scale and redshift. We find that scatter around the SFMS is important for all bands, especially at the highest redshifts where the emissivity is dominated by the faintest, most "bursty" galaxies. Ignoring this scatter would underestimate the mean emissivity and its standard deviation computed over 5 cMpc regions by factors of up to $\sim$2-10 at $5< z < 15$. Scatter around the X-ray luminosity to star formation rate relation is important for determining X-ray emissivity, accounting for roughly half of its mean and standard deviation. The importance of scatter in the ionizing escape fraction depends on its functional form, while scatter around the SHMR contributes at the level of ~10-20%. Although scatter does flatten the UV luminosity functions, shifting the bright end by 1-2 magnitudes, the level of scatter in our fiducial model is insufficient to fully explain recent estimates from JWST photometry (consistent with previous studies)., Comment: 17 pages, 13 figures, submitted to A&A

The attenuation of Ly$\alpha$ 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$\alpha$ 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$\alpha$ emission profiles (obtained via covariance matrix reconstruction of the Ly$\alpha$ 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$\alpha$ forest) along with other observational probes considered in the literature., Comment: 21 pages, 9 figures and 2 tables. Accepted for publication in MNRAS

Reconstructing the intrinsic Ly$\alpha$ line flux from high-$z$ QSOs can place constraints on the neutral hydrogen content of the intergalactic medium during reionisation. There are now $\gtrsim10$ different Ly$\alpha$ reconstruction pipelines using different methodologies to predict the Ly$\alpha$ line flux from correlations with the spectral information redward of Ly$\alpha$. However, there have been few attempts to directly compare the performance of these pipelines. Therefore, we devised a blind QSO challenge to compare these reconstruction pipelines on a uniform set of objects. Each author was provided de-identified, observed rest-frame QSO spectra with spectral information only redward of 1260\AA\ rest-frame to ensure unbiased reconstruction. We constructed two samples of 30 QSOs, from X-Shooter and SDSS both spanning $3.5

The detection of the 21-cm signal at $z\gtrsim6$ will reveal insights into the properties of the first galaxies responsible for driving reionisation. To extract this information, we perform parameter inference which requires embedding 3D simulations of the 21-cm signal within a Bayesian inference pipeline. Presently, when performing inference we must choose which sources of uncertainty to sample and which to hold fixed. Since the astrophysics of galaxies are much more uncertain than those of the underlying halo-mass function (HMF), we usually parameterise and model the former while fixing the latter. However, in doing so we may bias our inference of the properties of these first galaxies. In this work, we explore the consequences of assuming an incorrect choice of HMF and quantify the relative biases in our inferred astrophysical model parameters when considering the wrong HMF. We then relax this assumption by constructing a generalised five parameter model for the HMF and simultaneously recover these parameters along with our underlying astrophysical model. For this analysis, we use 21cmFAST and perform Simulation-Based Inference by applying marginal neural ratio estimation to learn the likelihood-to-evidence ratio using Swyft. Using a mock 1000 hour observation of the 21-cm power spectrum from the forthcoming Square Kilometre Array, conservatively assuming foreground wedge avoidance, we find assuming the incorrect HMF can bias the recovered astrophysical parameters by up to $\sim3-4\sigma$ even when including independent information from observed luminosity functions. When considering our generalised HMF model, we recover constraints on our astrophysical parameters with a factor of $\sim2-4$ larger marginalised uncertainties. Importantly, these constraints are unbiased, agnostic to the underlying HMF and therefore more conservative., Comment: 27 pages, 14 figures, 3 tables and 3 appendices. Submitted to MNRAS, comments welcome

Enlightening our understanding of the first galaxies responsible for driving reionisation requires detecting the 21-cm signal from neutral hydrogen. Interpreting the wealth of information embedded in this signal requires Bayesian inference. Parameter inference from the 21-cm signal is primarily restricted to the spherically averaged power spectrum (1D PS) owing to its relatively straightforward derivation of an analytic likelihood function enabling traditional Monte-Carlo Markov-Chain (MCMC) approaches. However, in recent years, simulation-based inference (SBI) has become feasible which removes the necessity of having an analytic likelihood, enabling more complex summary statistics of the 21-cm signal to be used for Bayesian inference. In this work, we use SBI, specifically marginal neural ratio estimation to learn the likelihood-to-evidence ratio with Swyft, to explore parameter inference using the cylindrically averaged 2D PS. Since the 21-cm signal is anisotropic, the 2D PS should yield more constraining information compared to the 1D PS which isotropically averages the signal. For this, we consider a mock 1000 hr observation of the 21-cm signal using the SKA and compare the performance of the 2D PS relative to the 1D PS. Additionally, we explore two separate foreground mitigation strategies, perfect foreground removal and wedge avoidance. We find the 2D PS outperforms the 1D PS by improving the marginalised uncertainties on individual astrophysical parameters by up to $\sim30-40$ per cent irrespective of the foreground mitigation strategy. Primarily, these improvements stem from how the 2D PS distinguishes between the transverse, $k_{\perp}$, and redshift dependent, $k_{\parallel}$ information which enables greater sensitivity to the complex reionisation morphology., Comment: 16 pages, 6 figures and 1 table. Submitted to MNRAS, comments welcome

Radio interferometers targeting the 21cm brightness temperature fluctuations at high redshift are subject to systematic effects that operate over a range of different timescales. These can be isolated by designing appropriate Fourier filters that operate in fringe-rate (FR) space, the Fourier pair of local sidereal time (LST). Applications of FR filtering include separating effects that are correlated with the rotating sky vs. those relative to the ground, down-weighting emission in the primary beam sidelobes, and suppressing noise. FR filtering causes the noise contributions to the visibility data to become correlated in time however, making interpretation of subsequent averaging and error estimation steps more subtle. In this paper, we describe fringe rate filters that are implemented using discrete prolate spheroidal sequences, and designed for two different purposes -- beam sidelobe/horizon suppression (the `mainlobe' filter), and ground-locked systematics removal (the `notch' filter). We apply these to simulated data, and study how their properties affect visibilities and power spectra generated from the simulations. Included is an introduction to fringe-rate filtering and a demonstration of fringe-rate filters applied to simple situations to aid understanding., Comment: 21 pages, 18 figures, submitted to Monthly Notices of the Royal Astronomical Society

The cosmic 21-cm signal will bring data-driven advances to studies of the Cosmic Dawn (CD) and Epoch of Reionization (EoR). Radio telescopes such as the SKA will eventually map the HI fluctuations over the first billion years - the majority of our observable Universe. With such large data volumes, it becomes increasingly important to develop "optimal" summary statistics, allowing us to learn as much as possible about the CD and EoR. In this work we compare the constraining power of several 21-cm summary statistics, using the determinant of the Fisher information matrix, $\det F$. Since we do not have an established "fiducial" model for the astrophysics of the first galaxies, we compute the distribution of $\det F$ across the prior volume. Using a large database of cosmic 21-cm lightcones that include realizations of telescope noise, we compare the following summaries: (i) the spherically-averaged power spectrum (1DPS), (ii) the cylindrically-averaged power spectrum (2DPS), (iii) the 2D Wavelet scattering transform (WST), (iv) a recurrent neural network (RNN), (v) an information-maximizing neural network (IMNN), and (vi) the combination of 2DPS and IMNN. Our best performing individual summary is the 2DPS, having relatively high Fisher information throughout parameter space. Although capable of achieving the highest Fisher information for some parameter choices, the IMNN does not generalize well, resulting in a broad distribution. Our best results are achieved with the concatenation of the 2DPS and IMNN. The combination of only these two complimentary summaries reduces the recovered parameter variances on average by factors of $\sim$6.5 - 9.5, compared with using each summary independently. Finally, we point out that that the common assumption of a constant covariance matrix when doing Fisher forecasts using 21-cm summaries can significantly underestimate parameter constraints., Comment: 11 pages, 8 figures, submitted to A&A

This paper presents the design and deployment of the Hydrogen Epoch of Reionization Array (HERA) phase II system. HERA is designed as a staged experiment targeting 21 cm emission measurements of the Epoch of Reionization. First results from the phase I array are published as of early 2022, and deployment of the phase II system is nearing completion. We describe the design of the phase II system and discuss progress on commissioning and future upgrades. As HERA is a designated Square Kilometer Array (SKA) pathfinder instrument, we also show a number of "case studies" that investigate systematics seen while commissioning the phase II system, which may be of use in the design and operation of future arrays. Common pathologies are likely to manifest in similar ways across instruments, and many of these sources of contamination can be mitigated once the source is identified.

Detection of the faint 21 cm line emission from the Cosmic Dawn and Epoch of Reionisation will require not only exquisite control over instrumental calibration and systematics to achieve the necessary dynamic range of observations but also validation of analysis techniques to demonstrate their statistical properties and signal loss characteristics. A key ingredient in achieving this is the ability to perform high-fidelity simulations of the kinds of data that are produced by the large, many-element, radio interferometric arrays that have been purpose-built for these studies. The large scale of these arrays presents a computational challenge, as one must simulate a detailed sky and instrumental model across many hundreds of frequency channels, thousands of time samples, and tens of thousands of baselines for arrays with hundreds of antennas. In this paper, we present a fast matrix-based method for simulating radio interferometric measurements (visibilities) at the necessary scale. We achieve this through judicious use of primary beam interpolation, fast approximations for coordinate transforms, and a vectorised outer product to expand per-antenna quantities to per-baseline visibilities, coupled with standard parallelisation techniques. We validate the results of this method, implemented in the publicly-available matvis code, against a high-precision reference simulator, and explore its computational scaling on a variety of problems., Comment: 25 pages, 20 figures, submitted to RAS Techniques and Instruments, matvis is publicly available at https://github.com/HERA-Team/matvis

Observations with radio arrays that target the 21-cm signal originating from the early Universe suffer from a variety of systematic effects. An important class of these are reflections and spurious couplings between antennas. We apply a Hamiltonian Monte Carlo sampler to the modelling and mitigation of these systematics in simulated Hydrogen Epoch of Reionisation Array (HERA) data. This method allows us to form statistical uncertainty estimates for both our models and the recovered visibilities, which is an important ingredient in establishing robust upper limits on the Epoch of Reionisation (EoR) power spectrum. In cases where the noise is large compared to the EoR signal, this approach can constrain the systematics well enough to mitigate them down to the noise level for both systematics studied. Where the noise is smaller than the EoR, our modelling can mitigate the majority of the reflections with there being only a minor level of residual systematics, while cross-coupling sees essentially complete mitigation. Our approach performs similarly to existing filtering/fitting techniques used in the HERA pipeline, but with the added benefit of rigorously propagating uncertainties. In all cases it does not significantly attenuate the underlying signal., Comment: 19 pages, 14 figures, submitted to MNRAS

The key to detecting neutral hydrogen during the epoch of reionization (EoR) is to separate the cosmological signal from the dominating foreground radiation. We developed direct optimal mapping (DOM) to map interferometric visibilities; it contains only linear operations, with full knowledge of point spread functions from visibilities to images. Here, we demonstrate a fast Fourier transform-based image power spectrum and its window functions computed from the DOM images. We use noiseless simulation, based on the Hydrogen Epoch of Reionization Array Phase I configuration, to study the image power spectrum properties. The window functions show $<10^{-11}$ of the integrated power leaks from the foreground-dominated region into the EoR window; the 2D and 1D power spectra also verify the separation between the foregrounds and the EoR., Comment: Published in ApJ

Recent years have witnessed rapid progress in observations of the Epoch of Reionization (EoR). These have enabled high-dimensional inference of galaxy and intergalactic medium (IGM) properties during the first billion years of our Universe. However, even using efficient, semi-numerical simulations, traditional inference approaches that compute 3D lightcones on-the-fly can take $10^5$ core hours. Here we present 21cmEMU: an emulator of several summary observables from the popular 21cmFAST simulation code. 21cmEMU takes as input nine parameters characterizing EoR galaxies, and outputs the following summary statistics: (i) the IGM mean neutral fraction; (ii) the 21-cm power spectrum; (iii) the mean 21-cm spin temperature; (iv) the sky-averaged (global) 21-cm signal; (v) the ultraviolet (UV) luminosity functions (LFs); and (vi) the Thomson scattering optical depth to the cosmic microwave background (CMB). All observables are predicted with sub-percent median accuracy, with a reduction of the computational cost by a factor of over 10$^4$. After validating inference results, we showcase a few applications, including: (i) quantifying the relative constraining power of different observational datasets; (ii) seeing how recent claims of a late EoR impact previous inferences; and (iii) forecasting upcoming constraints from the sixth observing season of the Hydrogen Epoch of Reionization Array (HERA) telescope. 21cmEMU is publicly-available, and is included as an alternative simulator in the public 21CMMC sampler., Comment: 21 pages, 13 figures, submitted to MNRAS

Recently, the Lyman-$\alpha$ (Ly$\alpha$) forest flux auto-correlation function has been shown to be sensitive to the mean free path of hydrogen-ionizing photons, $\lambda_{\text{mfp}}$, for simulations at $z \geq 5.4$. Measuring $\lambda_{\text{mfp}}$ at these redshifts will give vital information on the ending of reionization. Here we present the first observational measurements of the Ly$\alpha$ 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$\alpha$ 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 $\lambda_{\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 $\lambda_{\text{mfp}}$ from this data., Comment: 17 pages, 13 figures, submitted to MNRAS

The redshifted 21cm signal from the Cosmic Dawn is expected to provide unprecedented insights into early Universe astrophysics and cosmology. Here we explore how dark matter can heat the intergalactic medium before the first galaxies, leaving a distinctive imprint in the 21cm power spectrum. We provide the first dedicated Fisher matrix forecasts on the sensitivity of the Hydrogen Epoch of Reionization Array (HERA) telescope to dark matter decays. We show that with 1000 hours of observation, HERA has the potential to improve current cosmological constraints on the dark matter decay lifetime by up to three orders of magnitude. Even in extreme scenarios with strong X-ray emission from early-forming, metal-free galaxies, the bounds on the decay lifetime would be improved by up to two orders of magnitude. Overall, HERA shall improve on existing limits for dark matter masses below $2$ GeV$/c^2$ for decays into $e^+e^-$ and below few MeV$/c^2$ for decays into photons., Comment: 33 pages, 12 figures

The mean free path of ionizing photons, $\lambda_{\rm mfp}$, is a critical parameter for modeling the intergalactic medium (IGM) both during and after reionization. We present direct measurements of $\lambda_{\rm mfp}$ from QSO spectra over the redshift range $5

The patchy kinetic Sunyaev-Zel'dovich (kSZ) signal is an integral probe of the timing and morphology of the epoch of reionization (EoR). Recent observations have claimed a low signal-to-noise (S/N) measurement, with a dramatic increase in S/N expected in the near future. In this work, we quantify what we can learn about the EoR from the kSZ signal. We perform Bayesian inference by sampling galaxy properties and using forward-models of the kSZ as well as other EoR and galaxy observations in the likelihood. Including the recent kSZ measurement obtained by the South Pole Telescope ($\mathcal{D}_{3000}^{\rm{pkSZ}} = 1.1_{-0.7}^{+1.1} \mu$K$^2$) shifts the posterior distribution in favor of faster and later reionization models, resulting in lower values of the optical depth to the CMB: $\tau_e = 0.052_{-0.008}^{+0.009}$ with a 68$\%$ confidence interval (C.I.). The combined EoR and UV luminosity function observations also imply a typical ionizing escape fraction of $0.04_{-0.03}^{+0.05}$ (95$\%$ C.I.), without a strong dependence on halo mass. We show how the patchy kSZ power from our posterior depends on the commonly-used parameters of reionization. For a given midpoint and duration, the EoR morphology only has a few percent impact on the patchy kSZ power in our posterior. However, a physical model is needed to obtain tight constraints from the current low S/N patchy kSZ measurement, as it allows us to take advantage of complimentary high-$z$ observations. Future high S/N detections of the patchy kSZ should decrease the current uncertainties on the timing of the EoR by factors of $\sim$2 - 3., Comment: 14 pages, 6 figures, submitted to MNRAS

We explore how the characteristics of the cross-correlation functions between the 21cm emission from the spin-flip transition of neutral hydrogen (HI) and early Lyman-$\alpha$ (Ly$\alpha$) radiation emitting galaxies (Ly$\alpha$ emitters, LAEs) depend on the reionisation history and topology and the simulated volume. For this purpose, we develop an analytic expression for the 21cm-LAE cross-correlation function and compare it to results derived from different Astraeus and 21cmFAST reionisation simulations covering a physically plausible range of scenarios where either low-mass ($<10^{9.5}M_\odot$) or massive ($>10^{9.5}M_\odot$) galaxies drive reionisation. Our key findings are: (i) the negative small-scale ($<2$ cMpc) cross-correlation amplitude scales with the intergalactic medium's (IGM) average HI fraction ($\langle\chi_\mathrm{HI}\rangle$) and spin-temperature weighted overdensity in neutral regions ($\langle1+\delta\rangle_\mathrm{HI}$); (ii) the inversion point of the cross-correlation function traces the peak of the size distribution of ionised regions around LAEs; (iii) the cross-correlation amplitude at small scales is sensitive to the reionisation topology, with its anti-correlation or correlation decreasing the stronger the ionising emissivity of the underlying galaxy population is correlated to the cosmic web gas distribution (i.e. the more low-mass galaxies drive reionisation); (iv) the required simulation volume to not underpredict the 21cm-LAE anti-correlation amplitude when the cross-correlation is derived via the cross-power spectrum rises as the size of ionised regions and their variance increases. Our analytic expression can serve two purposes: to test whether simulation volumes are sufficiently large, and to act as a fitting function when cross-correlating future 21cm signal Square Kilometre Array and LAE galaxy observations., Comment: 13 pages, 5 figures; accepted for publication in MNRAS

We investigate the impact of field-to-field variation, deriving from cosmic variance, in measured Lyman-$\alpha$ emitter (LAE) luminosity functions (LFs) and this variation's impact on inferences of the neutral fraction of the intergalactic medium (IGM) during reionization. We post-process a z=7 IGM simulation to populate the dark matter halos with LAEs. These LAEs have realistic UV magnitudes, Ly$\alpha$ fluxes, and Ly$\alpha$ line profiles. We calculate the attenuation of Ly$\alpha$ emission in universes with varying IGM neutral fraction, $\bar{\rm{x}}_{\rm{HI}}$. In a $\bar{\rm{x}}_{\rm{HI}}=0.3$ simulation, we perform 100 realizations of a mock 2 square degree survey with a redshift window $\Delta z = 0.5$ and flux limit $\rm{f}_{Ly\alpha}>1\times10^{-17}\:\rm{ergs}\:\: \rm{s}^{-1} \: \rm{cm}^{-2}$; such a survey is typical in depth and volume of the largest LAE surveys conducted today. For each realization, we compute the LAE LF and use it to recover the input $\bar{\rm{x}}_{\rm{HI}}$. Comparing the inferred values of $\bar{\rm{x}}_{\rm{HI}}$ across the ensemble of the surveys, we find that cosmic variance, deriving from large-scale structure and variation in the neutral gas along the sightline, imposes a floor in the uncertainty of $\Delta \bar{\rm{x}}_{\rm{HI}} \sim 0.2$ when $\bar{\rm{x}}_{\rm{HI}}$ $=0.3$. We explore mitigation strategies to decrease this uncertainty, such as increasing the volume, decreasing the flux limit, or probing the volume with many independent fields. Increasing the area and/or depth of the survey does not mitigate the uncertainty, but composing a survey with many independent fields is effective. This finding highlights the best strategy for LAE surveys aiming at constraining $\bar{\rm{x}}_{\rm{HI}}$ of the universe during reionization., Comment: 17 pages, 13 figures

Observations of the cosmic 21-cm power spectrum (PS) are starting to enable precision Bayesian inference of galaxy properties and physical cosmology, during the first billion years of our Universe. Here we investigate the impact of common approximations about the likelihood used in such inferences, including: (i) assuming a Gaussian functional form; (ii) estimating the mean from a single realization; and (iii) estimating the (co)variance at a single point in parameter space. We compare "classical" inference that uses an explicit likelihood with simulation based inference (SBI) that estimates the likelihood from a training set. Our forward-models include: (i) realizations of the cosmic 21-cm signal computed with 21cmFAST by varying UV and X-ray galaxy parameters together with the initial conditions; (ii) realizations of the telescope noise corresponding to a 1000 h integration with SKA1-Low; (iii) the excision of Fourier modes corresponding to a foreground-dominated, horizon "wedge". We find that the 1D PS likelihood is well described by a Gaussian accounting for covariances between wavemodes and redshift bins (higher order correlations are small). However, common approaches of estimating the forward-modeled mean and (co)variance from a random realization or at a single point in parameter space result in biased and over-constrained posteriors. Our best results come from using SBI to fit a non-Gaussian likelihood with a Gaussian mixture neural density estimator. Such SBI can be performed with up to an order of magnitude fewer simulations than classical, explicit likelihood inference. Thus SBI provides accurate posteriors at a comparably low computational cost., Comment: 16 pages, 12 figures, updated to match the version accepted for publication in MNRAS

Constraining when and how reionisation began is pivotal for understanding when the first galaxies formed. Lyman-alpha (Ly$\alpha$) emission from galaxies is currently our most promising probe of these early stages. At z>7 the majority of galaxies detected with Ly$\alpha$ are in candidate overdensities. Here we quantify the probability of these galaxies residing in large ionised bubbles. We create (1.6 Gpc)$^3$ reionising intergalactic medium (IGM) simulations, providing sufficient volume to robustly measure bubble size distributions around UV-bright galaxies and rare overdensities. We find $M_{\rm UV} \lesssim -16$ galaxies and overdensities are $\gtrsim$10-1000x more likely to trace ionised bubbles compared to randomly selected positions. The brightest galaxies and strongest overdensities have bubble size distributions with highest characteristic size and least scatter. We compare two models: gradual reionisation driven by numerous UV-faint galaxies versus more rapid reionisation by rarer brighter galaxies, producing larger bubbles at fixed neutral fraction. We demonstrate that recently observed z~7 overdensities are highly likely to trace large ionised bubbles, corroborated by their high Ly$\alpha$ detection rates. However, the z~8.7 association of Ly$\alpha$ emitters in EGS and GN-z11, with Ly$\alpha$ at z=10.6, are unlikely to trace large bubbles in our fiducial model -- 11% and 7% probability of >1 proper Mpc bubbles, respectively. Ly$\alpha$ detections at such high redshifts could be explained by: a less neutral IGM than previously expected; larger ionised regions at fixed neutral fraction; or if intrinsic Ly$\alpha$ flux is unusually strong in these galaxies. We discuss how to test these scenarios with JWST and the prospects for using upcoming wide-area surveys to distinguish between reionisation models., Comment: 15 pages (+ 3 pages in Appendix), 13 figures, submitted to MNRAS

The upcoming Square Kilometre Array Observatory (SKAO) will produce images of neutral hydrogen distribution during the epoch of reionization by observing the corresponding 21-cm signal. However, the 21-cm signal will be subject to instrumental limitations such as noise and galactic foreground contamination which pose a challenge for accurate detection. In this study, we present the SegU-Net v2 framework, an enhanced version of our convolutional neural network, built to identify neutral and ionized regions in the 21-cm signal contaminated with foreground emission. We trained our neural network on 21-cm image data processed by a foreground removal method based on Principal Component Analysis achieving an average classification accuracy of 71 per cent between redshift $z=7$ to $11$. We tested SegU-Net v2 against various foreground removal methods, including Gaussian Process Regression, Polynomial Fitting, and Foreground-Wedge Removal. Results show comparable performance, highlighting SegU-Net v2's independence on these pre-processing methods. Statistical analysis shows that a perfect classification score with $AUC=95\%$ is possible for $8(10\, {\rm cMpc})^3$ at $z>9$, for follow-up studies with infrared/optical telescopes to detect these sources., Comment: 19 pages, 10 figures, 3 tables, 2 appendixes

Radio interferometers aiming to measure the power spectrum of the redshifted 21 cm line during the Epoch of Reionisation (EoR) need to achieve an unprecedented dynamic range to separate the weak signal from overwhelming foreground emissions. Calibration inaccuracies can compromise the sensitivity of these measurements to the effect that a detection of the EoR is precluded. An alternative to standard analysis techniques makes use of the closure phase, which allows one to bypass antenna-based direction-independent calibration. Similarly to standard approaches, we use a delay spectrum technique to search for the EoR signal. Using 94 nights of data observed with Phase I of the Hydrogen Epoch of Reionization Array (HERA), we place approximate constraints on the 21 cm power spectrum at $z=7.7$. We find at 95% confidence that the 21 cm EoR brightness temperature is $\le$(372)$^2$ "pseudo" mK$^2$ at 1.14 "pseudo" $h$ Mpc$^{-1}$, where the "pseudo" emphasises that these limits are to be interpreted as approximations to the actual distance scales and brightness temperatures. Using a fiducial EoR model, we demonstrate the feasibility of detecting the EoR with the full array. Compared to standard methods, the closure phase processing is relatively simple, thereby providing an important independent check on results derived using visibility intensities, or related., Comment: 16 pages, 14 figures, accepted for publication by MNRAS

The 21-cm signal from neutral hydrogen in the early universe will provide unprecedented information about the first stars and galaxies. Extracting this information, however, requires accounting for many unknown astrophysical processes. Semi-numerical simulations are key for exploring the vast parameter space of said processes. These simulations use approximate techniques such as excursion-set and perturbation theory to model the 3D evolution of the intergalactic medium, at a fraction of the computational cost of hydrodynamic and/or radiative transfer simulations. However, exploring the enormous parameter space of the first galaxies can still be computationally expensive. Here we introduce 21cmfish, a Fisher-matrix wrapper for the semi-numerical simulation 21cmFAST. 21cmfish facilitates efficient parameter forecasts, scaling to significantly higher dimensionalities than MCMC approaches, assuming a multi-variate Gaussian posterior. Our method produces comparable parameter uncertainty forecasts to previous MCMC analyses but requires ~10$^4$x fewer simulations. This enables a rapid way to prototype analyses adding new physics and/or additional parameters. We carry out a forecast for HERA using the largest astrophysical parameter space to-date, with 10 free parameters, spanning both population II and III star formation. We find X-ray parameters for the first galaxies could be measured to sub-percent precision, and, though they are highly degenerate, the stellar-to-halo mass relation and ionizing photon escape fraction for population II and III galaxies can be constrained to ~10% precision (logarithmic quantities). Using a principal component analysis we find HERA is most sensitive to the product of the ionizing escape fraction and the stellar-to-halo mass fraction for population II galaxies., Comment: 18 pages, accepted for publication in MNRAS. Fisher-matrix wrapper for 21cmFAST available here: https://github.com/charlottenosam/21cmfish

Radio Frequency Interference (RFI) is one of the systematic challenges preventing 21cm interferometric instruments from detecting the Epoch of Reionization. To mitigate the effects of RFI on data analysis pipelines, numerous inpaint techniques have been developed to restore RFI corrupted data. We examine the qualitative and quantitative errors introduced into the visibilities and power spectrum due to inpainting. We perform our analysis on simulated data as well as real data from the Hydrogen Epoch of Reionization Array (HERA) Phase 1 upper limits. We also introduce a convolutional neural network that capable of inpainting RFI corrupted data in interferometric instruments. We train our network on simulated data and show that our network is capable at inpainting real data without requiring to be retrained. We find that techniques that incorporate high wavenumbers in delay space in their modeling are best suited for inpainting over narrowband RFI. We also show that with our fiducial parameters Discrete Prolate Spheroidal Sequences (DPSS) and CLEAN provide the best performance for intermittent ``narrowband'' RFI while Gaussian Progress Regression (GPR) and Least Squares Spectral Analysis (LSSA) provide the best performance for larger RFI gaps. However we caution that these qualitative conclusions are sensitive to the chosen hyperparameters of each inpainting technique. We find these results to be consistent in both simulated and real visibilities. We show that all inpainting techniques reliably reproduce foreground dominated modes in the power spectrum. Since the inpainting techniques should not be capable of reproducing noise realizations, we find that the largest errors occur in the noise dominated delay modes. We show that in the future, as the noise level of the data comes down, CLEAN and DPSS are most capable of reproducing the fine frequency structure in the visibilities of HERA data., Comment: 21 pages, 13 figures

We report the most sensitive upper limits to date on the 21 cm epoch of reionization power spectrum using 94 nights of observing with Phase I of the Hydrogen Epoch of Reionization Array (HERA). Using similar analysis techniques as in previously reported limits (HERA Collaboration 2022a), we find at 95% confidence that $\Delta^2(k = 0.34$ $h$ Mpc$^{-1}$) $\leq 457$ mK$^2$ at $z = 7.9$ and that $\Delta^2 (k = 0.36$ $h$ Mpc$^{-1}) \leq 3,496$ mK$^2$ at $z = 10.4$, an improvement by a factor of 2.1 and 2.6 respectively. These limits are mostly consistent with thermal noise over a wide range of $k$ after our data quality cuts, despite performing a relatively conservative analysis designed to minimize signal loss. Our results are validated with both statistical tests on the data and end-to-end pipeline simulations. We also report updated constraints on the astrophysics of reionization and the cosmic dawn. Using multiple independent modeling and inference techniques previously employed by HERA Collaboration (2022b), we find that the intergalactic medium must have been heated above the adiabatic cooling limit at least as early as $z = 10.4$, ruling out a broad set of so-called "cold reionization" scenarios. If this heating is due to high-mass X-ray binaries during the cosmic dawn, as is generally believed, our result's 99% credible interval excludes the local relationship between soft X-ray luminosity and star formation and thus requires heating driven by evolved low-metallicity stars., Comment: 57 pages, 37 figures. Updated to match the accepted ApJ version. Corresponding author: Joshua S. Dillon

Combining the visibilities measured by an interferometer to form a cosmological power spectrum is a complicated process. In a delay-based analysis, the mapping between instrumental and cosmological space is not a one-to-one relation. Instead, neighbouring modes contribute to the power measured at one point, with their respective contributions encoded in the window functions. To better understand the power measured by an interferometer, we assess the impact of instrument characteristics and analysis choices on these window functions. Focusing on the Hydrogen Epoch of Reionization Array (HERA) as a case study, we find that long-baseline observations correspond to enhanced low-k tails of the window functions, which facilitate foreground leakage, whilst an informed choice of bandwidth and frequency taper can reduce said tails. With simple test cases and realistic simulations, we show that, apart from tracing mode mixing, the window functions help accurately reconstruct the power spectrum estimator of simulated visibilities. The window functions depend strongly on the beam chromaticity, and less on its spatial structure - a Gaussian approximation, ignoring side lobes, is sufficient. Finally, we investigate the potential of asymmetric window functions, down-weighting the contribution of low-k power to avoid foreground leakage. The window functions presented here correspond to the latest HERA upper limits for the full Phase I data. They allow an accurate reconstruction of the power spectrum measured by the instrument and will be used in future analyses to confront theoretical models and data directly in cylindrical space., Comment: 18 pages, 19 figures, accepted for publication in MNRAS

Motivated by the desire for wide-field images with well-defined statistical properties for 21cm cosmology, we implement an optimal mapping pipeline that computes a maximum likelihood estimator for the sky using the interferometric measurement equation. We demonstrate this direct optimal mapping with data from the Hydrogen Epoch of Reionization (HERA) Phase I observations. After validating the pipeline with simulated data, we develop a maximum likelihood figure-of-merit for comparing four sky models at 166MHz with a bandwidth of 100kHz. The HERA data agree with the GLEAM catalogs to <10%. After subtracting the GLEAM point sources, the HERA data discriminate between the different continuum sky models, providing most support for the model of Byrne et al. 2021. We report the computation cost for mapping the HERA Phase I data and project the computation for the HERA 320-antenna data; both are feasible with a modern server. The algorithm is broadly applicable to other interferometers and is valid for wide-field and non-coplanar arrays., Comment: 16 pages, 10 figures, 2 tables, published on ApJ

High-energy emission associated with star formation has been proposed as a significant source of interstellar medium (ISM) ionization in low-metallicity starbursts and an important contributor to the heating of the intergalactic medium (IGM) in the high-redshift ($z > 8$) Universe. Using Chandra observations of a sample of 30 galaxies at $D \approx$~200--450 Mpc that have high specific star-formation rates of 3--9 Gyr$^{-1}$ and metallicities near $Z \approx 0.3 Z_\odot$, we provide new measurements of the average 0.5--8 keV spectral shape and normalization per unit star-formation rate (SFR). We model the sample-combined X-ray spectrum as a combination of hot gas and high-mass X-ray binary (HMXB) populations and constrain their relative contributions. We derive scaling relations of $\log L_{\rm 0.5-8 keV}^{\rm HMXB}$/SFR $= 40.19 \pm 0.06$ and $\log L_{\rm 0.5-2 keV}^{\rm gas}$/SFR $= 39.58^{+0.17}_{-0.28}$; significantly elevated compared to local relations. The HMXB scaling is also somewhat higher than $L_{\rm 0.5-8 keV}^{\rm HMXB}$-SFR-$Z$ relations presented in the literature, potentially due to our galaxies having relatively low HMXB obscuration and young and X-ray luminous stellar populations. The elevation of the hot gas scaling relation is at the level expected for diminished attenuation due to a reduction of metals; however, we cannot conclude that an $L_{\rm 0.5-2 keV}^{\rm gas}$-SFR-$Z$ relation is driven solely by changes in ISM metal content. Finally, we present SFR-scaled spectral models (both emergent and intrinsic) that span the X-ray--to--IR band, providing new benchmarks for studies of the impact of ISM ionization and IGM heating in the early Universe., Comment: Accepted for publication in ApJ (25 pages, 9 figures, and 4 tables)

X-rays from High-Mass X-ray Binaries (HMXBs) are likely the main source of heating of the intergalactic medium (IGM) during Cosmic Dawn (CD), before the completion of reionization. This Epoch of Heating (EoH; $z\sim 10-15$) should soon be detected via the redshifted 21-cm line from neutral hydrogen, allowing us to indirectly study the properties of HMXBs in the unseen, first galaxies. Low-redshift observations, as well as theoretical models, imply that the integrated X-ray luminosity to star formation rate of HMXBs ($L_{\rm X}/{\rm SFR}$) should increase in metal-poor environments, typical of early galaxies. Here we study the impact of the metallicity ($Z$) dependence of $L_{\rm X}/{\rm SFR}$ during the EoH. For our fiducial models, galaxies with star formation rates of order $10^{-3} - 10^{-1}$ $M_\odot$ yr$^{-1}$ and metallicities of order $10^{-3} - 10^{-2}$ $Z_\odot$ are the dominant contributors to the X-ray background (XRB) during this period. Different $L_{\rm X}/{\rm SFR}$-$Z$ relations result in factors of $\sim$ 3 differences in these ranges, as well as in the mean IGM temperature and the large-scale 21-cm power, at a given redshift. We compute mock 21-cm observations adopting as a baseline a 1000h integration with the upcoming Square Kilometer Array (SKA), for two different $L_{\rm X}/{\rm SFR}$-$Z$ relations. We perform inference on these mock observations using the common simplification of a constant $L_{\rm X}/{\rm SFR}$, finding that constant $L_{\rm X}/{\rm SFR}$ models can recover the IGM evolution of the more complicated $L_{\rm X}/{\rm SFR}$-$Z$ simulations only during the EoH. At $z<10$, where the typical galaxies are more polluted, constant $L_{\rm X}/{\rm SFR}$ models over-predict the XRB and its relative contribution to the early stages of the reionization.

A limiting systematic effect in 21-cm interferometric experiments is the chromaticity due to the coupling between the sky and the instrument. This coupling is sourced by the instrument primary beam; therefore it is important to know the beam to extremely high precision. Here we demonstrate how known beam uncertainties can be characterized using databases of beam models. In this introductory work, we focus on beam errors arising from physically offset and/or broken antennas within a station. We use the public code OSKAR to generate an "ideal" SKA beam formed from 256 antennas regularly-spaced in a 35-m circle, as well as a large database of "perturbed" beams sampling distributions of broken/offset antennas. We decompose the beam errors ("ideal" minus "perturbed") using Principal Component Analysis (PCA) and Kernel PCA (KPCA). Using 20 components, we find that PCA/KPCA can reduce the residual of the beam in our datasets by 60-90% compared with the assumption of an ideal beam. Using a simulated observation of the cosmic signal plus foregrounds, we find that assuming the ideal beam can result in 1% error in the EoR window and 10% in the wedge of the 2D power spectrum. When PCA/KPCA is used to characterize the beam uncertainties, the error in the power spectrum shrinks to below 0.01% in the EoR window and <1% in the wedge. Our framework can be used to characterize and then marginalize over uncertainties in the beam for robust next-generation 21-cm parameter estimation., Comment: 12 pages, 9 main figures

Upper limits from the current generation of interferometers targeting the 21-cm signal from high redshifts have recently begun to rule out physically realistic, though still extreme, models of the Epoch of Reionization (EoR). While inferring the detailed properties of the first galaxies is one of the most important motivations for measuring the high-$z$ 21-cm signal, they can also provide useful constraints on the properties of the intergalactic medium (IGM). Motivated by this, we build a simple, phenomenological model for the 21-cm power spectrum that works directly in terms of IGM properties, which bypasses the computationally expensive 3-D semi-numerical modeling generally employed in inference pipelines and avoids explicit assumptions about galaxy properties. The key simplifying assumptions are that (i) the ionization field is binary, and composed of spherical bubbles with an abundance described well by a parametric bubble size distribution, and (ii) that the spin temperature of the ``bulk'' IGM outside bubbles is uniform. Despite the simplicity of the model, the mean ionized fraction and spin temperature of the IGM recovered from mock 21-cm power spectra generated with \textsc{21cmfast} are generally in good agreement with the true input values. This suggests that it is possible to obtain comparable constraints on the IGM using models with very different assumptions, parameters, and priors. Our approach will thus be complementary to semi-numerical models as upper limits continue to improve in the coming years., Comment: 23 pages, 14 figures, accepted for publication in MNRAS, code publicly available at https://github.com/mirochaj/micro21cm

The excursion set reionization framework is widely used, due to its speed and accuracy in reproducing the 3D topology of reionization. However, it is known that it does not conserve photon number. Here, we introduce an efficient, on-the-fly recipe to approximately account for photon conservation. Using a flexible galaxy model shown to reproduce current high-$z$ observables, we quantify the bias in the inferred reionization history and galaxy properties resulting from the non-conservation of ionizing photons. Using a mock 21-cm observation, we perform inference with and without correcting for ionizing photon conservation. We find that ignoring photon conservation results in very modest biases in the inferred galaxy properties, for our fiducial model. The notable exception is in the power-law scaling of the ionizing escape fraction with halo mass, which can be biased from the true value by $\sim2.4\sigma$ (corresponding to $\sim0.2$ in the power-law index). Our scheme is implemented in the public code ${\tt 21cmFAST}$., Comment: 9 pages, 4 figures, 1 table, accepted for publication in MNRAS

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

The formation of the first galaxies during cosmic dawn and reionization (at redshifts $z=5-30$), triggered the last major phase transition of our universe, as hydrogen evolved from cold and neutral to hot and ionized. The 21-cm line of neutral hydrogen will soon allow us to map these cosmic milestones and study the galaxies that drove them. To aid in interpreting these observations, we upgrade the publicly available code {\tt 21cmFAST}. We introduce a new, flexible parametrization of the additive feedback from: an inhomogeneous, $H_2$-dissociating (Lyman-Werner; LW) background; and dark matter -- baryon relative velocities; which recovers results from recent, small-scale hydrodynamical simulations with both effects. We perform a large, "best-guess" simulation as the 2021 installment of the Evolution of 21-cm Structure (EOS) project. This improves the previous release with a galaxy model that reproduces the observed UV luminosity functions (UVLFs), and by including a population of molecular-cooling galaxies. The resulting 21-cm global signal and power spectrum are significantly weaker, primarily due to a more rapid evolution of the star-formation rate density required to match the UVLFs. Nevertheless, we forecast high signal-to-noise detections for both HERA and the SKA. We demonstrate how the stellar-to-halo mass relation of the unseen, first galaxies can be inferred from the 21-cm evolution. Finally, we show that the spatial modulation of X-ray heating due to relative velocities provides a unique acoustic signature that is detectable at $z \approx 10-15$ in our fiducial model. Ours are the first public simulations with joint inhomogeneous LW and relative-velocity feedback across the entire cosmic dawn and reionization., Comment: 23+3 pages, 23+3 figures. Matches accepted version. Simulation output available at https://www.dropbox.com/sh/dqh9r6wb0s68jfo/AACc9ZCqsN0SQ_JJN7GRVuqDa?dl=0

Measuring the density of the intergalactic medium using quasar sightlines in the epoch of reionization is challenging due to the saturation of Lyman-$\alpha$ absorption. Near a luminous quasar, however, the enhanced radiation creates a proximity zone observable in the quasar spectra where the Lyman-$\alpha$ 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 $\sigma_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 $\sigma_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}$., Comment: 16 pages, 11 figures, submitted to ApJ, comments are welcome

We present a framework for identifying and flagging malfunctioning antennas in large radio interferometers. We outline two distinct categories of metrics designed to detect outliers along known failure modes of large arrays: cross-correlation metrics, based on all antenna pairs, and auto-correlation metrics, based solely on individual antennas. We define and motivate the statistical framework for all metrics used, and present tailored visualizations that aid us in clearly identifying new and existing systematics. We implement these techniques using data from 105 antennas in the Hydrogen Epoch of Reionization Array (HERA) as a case study. Finally, we provide a detailed algorithm for implementing these metrics as flagging tools on real data sets., Comment: 31 pages, 17 figures

Recently, the Hydrogen Epoch of Reionization Array (HERA) collaboration has produced the experiment's first upper limits on the power spectrum of 21-cm fluctuations at z~8 and 10. Here, we use several independent theoretical models to infer constraints on the intergalactic medium (IGM) and galaxies during the epoch of reionization (EoR) from these limits. We find that the IGM must have been heated above the adiabatic cooling threshold by z~8, independent of uncertainties about the IGM ionization state and the nature of the radio background. Combining HERA limits with galaxy and EoR observations constrains the spin temperature of the z~8 neutral IGM to 27 K < T_S < 630 K (2.3 K < T_S < 640 K) at 68% (95%) confidence. They therefore also place a lower bound on X-ray heating, a previously unconstrained aspects of early galaxies. For example, if the CMB dominates the z~8 radio background, the new HERA limits imply that the first galaxies produced X-rays more efficiently than local ones (with soft band X-ray luminosities per star formation rate constrained to L_X/SFR = { 10^40.2, 10^41.9 } erg/s/(M_sun/yr) at 68% confidence), consistent with expectations of X-ray binaries in low-metallicity environments. The z~10 limits require even earlier heating if dark-matter interactions (e.g., through millicharges) cool down the hydrogen gas. Using a model in which an extra radio background is produced by galaxies, we rule out (at 95% confidence) the combination of high radio and low X-ray luminosities of L_{r,\nu}/SFR > 3.9 x 10^24 W/Hz/(M_sun/yr) and L_X/SFR<10^40 erg/s/(M_sun/yr). The new HERA upper limits neither support nor disfavor a cosmological interpretation of the recent EDGES detection. The analysis framework described here provides a foundation for the interpretation of future HERA results., Comment: 40 pages, 19 figures, accepted to ApJ

The presence of excess scatter in the Ly-$\alpha$ 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-$\alpha$ 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-$\alpha$ systems. We present improved measurements of the mean Ly-$\alpha$ transmission over $4.93.5 \sigma$). 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., Comment: Accepted for publication in MNRAS. Full measurement datasets are available on the first author's website

We report upper-limits on the Epoch of Reionization (EoR) 21 cm power spectrum at redshifts 7.9 and 10.4 with 18 nights of data ($\sim36$ hours of integration) from Phase I of the Hydrogen Epoch of Reionization Array (HERA). The Phase I data show evidence for systematics that can be largely suppressed with systematic models down to a dynamic range of $\sim10^9$ with respect to the peak foreground power. This yields a 95% confidence upper limit on the 21 cm power spectrum of $\Delta^2_{21} \le (30.76)^2\ {\rm mK}^2$ at $k=0.192\ h\ {\rm Mpc}^{-1}$ at $z=7.9$, and also $\Delta^2_{21} \le (95.74)^2\ {\rm mK}^2$ at $k=0.256\ h\ {\rm Mpc}^{-1}$ at $z=10.4$. At $z=7.9$, these limits are the most sensitive to-date by over an order of magnitude. While we find evidence for residual systematics at low line-of-sight Fourier $k_\parallel$ modes, at high $k_\parallel$ modes we find our data to be largely consistent with thermal noise, an indicator that the system could benefit from deeper integrations. The observed systematics could be due to radio frequency interference, cable sub-reflections, or residual instrumental cross-coupling, and warrant further study. This analysis emphasizes algorithms that have minimal inherent signal loss, although we do perform a careful accounting in a companion paper of the small forms of loss or bias associated with the pipeline. Overall, these results are a promising first step in the development of a tuned, instrument-specific analysis pipeline for HERA, particularly as Phase II construction is completed en route to reaching the full sensitivity of the experiment., Comment: Accepted to ApJ. https://reionization.org/science/public-data-release-1/

Imaging the cosmic 21 cm signal will map out the first billion years of our Universe. The resulting 3D lightcone (LC) will encode the properties of the unseen first galaxies and physical cosmology. Here, we build on previous work using neural networks (NNs) to infer astrophysical parameters directly from 21 cm LC images. We introduce recurrent neural networks (RNNs), capable of efficiently characterizing the evolution along the redshift axis of 21 cm LC images. Using a large database of simulated cosmic 21 cm LCs, we compare the relative performance in parameter estimation of different network architectures. These including two types of RNNs, which differ in their complexity, as well as a more traditional convolutional neural network (CNN). For the ideal case of no instrumental effects, our simplest and easiest to train RNN performs the best, with a mean squared parameter estimation error (MSE) that is lower by a factor of $\ge 2$ compared with the other architectures studied here, and a factor of $\ge 8$ lower than the previously-studied CNN. We also corrupt the cosmic signal by adding noise expected from a 1000 h integration with the Square Kilometre Array, as well as excising a foreground-contaminated 'horizon wedge'. Parameter prediction errors increase when the NNs are trained on these contaminated LC images, though recovery is still good even in the most pessimistic case (with $R^2 \ge 0.5-0.95$). However, we find no notable differences in performance between network architectures on the contaminated images. We argue this is due to the size of our data set, highlighting the need for larger data sets and/or better data augmentation in order to maximize the potential of NNs in 21 cm parameter estimation., Comment: 15 pages, 11 figures, updated to match the version published in MNRAS, minor changes

Precision calibration poses challenges to experiments probing the redshifted 21-cm signal of neutral hydrogen from the Cosmic Dawn and Epoch of Reionization (z~30-6). In both interferometric and global signal experiments, systematic calibration is the leading source of error. Though many aspects of calibration have been studied, the overlap between the two types of instruments has received less attention. We investigate the sky based calibration of total power measurements with a HERA dish and an EDGES style antenna to understand the role of auto-correlations in the calibration of an interferometer and the role of sky in calibrating a total power instrument. Using simulations we study various scenarios such as time variable gain, incomplete sky calibration model, and primary beam model. We find that temporal gain drifts, sky model incompleteness, and beam inaccuracies cause biases in the receiver gain amplitude and the receiver temperature estimates. In some cases, these biases mix spectral structure between beam and sky resulting in spectrally variable gain errors. Applying the calibration method to the HERA and EDGES data, we find good agreement with calibration via the more standard methods. Although instrumental gains are consistent with beam and sky errors similar in scale to those simulated, the receiver temperatures show significant deviations from expected values. While we show that it is possible to partially mitigate biases due to model inaccuracies by incorporating a time-dependent gain model in calibration, the resulting errors on calibration products are larger and more correlated. Completely addressing these biases will require more accurate sky and primary beam models., Comment: 16 pages, 13 figures, 1 table; accepted for publication in MNRAS main journal

We describe the validation of the HERA Phase I software pipeline by a series of modular tests, building up to an end-to-end simulation. The philosophy of this approach is to validate the software and algorithms used in the Phase I upper limit analysis on wholly synthetic data satisfying the assumptions of that analysis, not addressing whether the actual data meet these assumptions. We discuss the organization of this validation approach, the specific modular tests performed, and the construction of the end-to-end simulations. We explicitly discuss the limitations in scope of the current simulation effort. With mock visibility data generated from a known analytic power spectrum and a wide range of realistic instrumental effects and foregrounds, we demonstrate that the current pipeline produces power spectrum estimates that are consistent with known analytic inputs to within thermal noise levels (at the 2 sigma level) for k > 0.2 h/Mpc for both bands and fields considered. Our input spectrum is intentionally amplified to enable a strong `detection' at k ~0.2 h/Mpc -- at the level of ~25 sigma -- with foregrounds dominating on larger scales, and thermal noise dominating at smaller scales. Our pipeline is able to detect this amplified input signal after suppressing foregrounds with a dynamic range (foreground to noise ratio) of > 10^7. Our validation test suite uncovered several sources of scale-independent signal loss throughout the pipeline, whose amplitude is well-characterized and accounted for in the final estimates. We conclude with a discussion of the steps required for the next round of data analysis., Comment: 32 pages, 20 figures. Submitted to the Astrophysical Journal

The number density of small dark matter (DM) halos hosting faint high-redshift galaxies is sensitive to the DM free-streaming properties. However, constraining these DM properties is complicated by degeneracies with the uncertain baryonic physics governing star formation. In this work, we use a flexible astrophysical model and a Bayesian inference framework to analyse ultra-violet (UV) luminosity functions (LFs) at z=6-8. We vary the complexity of the galaxy model (single vs double power law for the stellar -- halo mass relation) as well as the matter power spectrum (cold DM vs thermal relic warm DM), comparing their Bayesian evidences. Adopting a conservatively wide prior range for the WDM particle mass, we show that the UV LFs at z=6-8 only weakly favour CDM over WDM. We find that particle masses of $\lesssim$2 keV are rejected at a 95% credible level in all models that have a WDM-like power spectrum cutoff. This bound should increase to ~2.5 keV with the James Webb Space Telescope (JWST)., Comment: prepared for submission to MNRAS

As current- and next-generation astronomical instruments come online, they will generate an unprecedented deluge of data. Analyzing these data in real time presents unique conceptual and computational challenges, and their long-term storage and archiving is scientifically essential for generating reliable, reproducible results. We present here the real-time processing (RTP) system for the Hydrogen Epoch of Reionization Array (HERA), a radio interferometer endeavoring to provide the first detection of the highly redshifted 21 cm signal from Cosmic Dawn and the Epoch of Reionization by an interferometer. The RTP system consists of analysis routines run on raw data shortly after they are acquired, such as calibration and detection of radio-frequency interference (RFI) events. RTP works closely with the Librarian, the HERA data storage and transfer manager which automatically ingests data and transfers copies to other clusters for post-processing analysis. Both the RTP system and the Librarian are public and open source software, which allows for them to be modified for use in other scientific collaborations. When fully constructed, HERA is projected to generate over 50 terabytes (TB) of data each night, and the RTP system enables the successful scientific analysis of these data., Comment: 15 pages, 3 figures, published in Astronomy and Computing

Precise measurements of the 21 cm power spectrum are crucial for understanding the physical processes of hydrogen reionization. Currently, this probe is being pursued by low-frequency radio interferometer arrays. As these experiments come closer to making a first detection of the signal, error estimation will play an increasingly important role in setting robust measurements. Using the delay power spectrum approach, we have produced a critical examination of different ways that one can estimate error bars on the power spectrum. We do this through a synthesis of analytic work, simulations of toy models, and tests on small amounts of real data. We find that, although computed independently, the different error bar methodologies are in good agreement with each other in the noise-dominated regime of the power spectrum. For our preferred methodology, the predicted probability distribution function is consistent with the empirical noise power distributions from both simulated and real data. This diagnosis is mainly in support of the forthcoming HERA upper limit, and also is expected to be more generally applicable., Comment: 35 Pages, 9 Figures, 4 Tables. Replaced with accepted ApJ version; some clarifying text added in response to referee comments with no changes to results

We present the first application of the isosceles bispectrum to MCMC parameter inference from the cosmic 21-cm signal. We extend the MCMC sampler 21cmMC to use the fast bispectrum code, BiFFT, when computing the likelihood. We create mock 1000h observations with SKA1-low, using PyObs21 to account for uv-sampling and thermal noise. Assuming the spin temperature is much higher than that of the CMB, we consider two different reionization histories for our mock observations: fiducial and late-reionization. For both models we find that bias on the inferred parameter means and 1-$\sigma$ credible intervals can be substantially reduced by using the isosceles bispectrum (calculated for a wide range of scales and triangle shapes) together with the power spectrum (as opposed to just using one of the statistics). We find that making the simplifying assumption of a Gaussian likelihood with a diagonal covariance matrix does not notably bias parameter constraints for the three-parameter reionization model and basic instrumental effects considered here. This is true even if we use extreme (unlikely) initial conditions which would be expected to amplify biases. We also find that using the cosmic variance error calculated with Monte-Carlo simulations using the fiducial model parameters whilst assuming the late-reionization model for the simulated data also does not strongly bias the inference. This implies we may be able to sparsely sample and interpolate the cosmic variance error over the parameter space, substantially reducing computational costs. All codes used in this work are publicly-available., Comment: 12 pages, 11 figures (submitted to MNRAS)

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

This brief code paper presents a new Python-wrapped version of the popular 21cm cosmology simulator, 21cmFAST. The new version, v3+, maintains the same core functionality of previous versions of 21cmFAST, but features a simple and intuitive interface, and a great deal more flexibility. This evolution represents the work of a formalized collaboration, and the new version, available publicly on GitHub, provides a single point-of-reference for all future upgrades and community-added features. In this paper, we describe simple usage of 21cmFAST, some of its new features, and provide a simple performance benchmark., Comment: 6 pages, 3 figures, published in JOSS

The first generation of galaxies is expected to form in minihalos, accreting gas through ${\rm H}_2$ cooling, and possessing unique properties. Although unlikely to be directly detected in UV/infrared surveys, the radiation from these molecular-cooling galaxies (MCGs) could leave an imprint in the 21-cm signal from the Cosmic Dawn. Here we quantify their detectability with upcoming radio interferometers. We generate mock 21-cm power spectra using a model for both MCGs as well as more massive, atomic-cooling galaxies (AGCs), allowing both populations to have different properties and scaling relations. The galaxy parameters are chosen so as to be consistent with: (i) high-redshift UV luminosity functions; (ii) the upper limit on the neutral fraction from QSO spectra; (iii) the Thomson scattering optical depth to the CMB; and (iv) the timing of the recent putative EDGES detection. The latter implies a significant contribution of MCGs to the Cosmic Dawn, if confirmed to be cosmological. We then perform Bayesian inference on two models including and ignoring MCG contributions. Comparing their Bayesian evidences, we find a strong preference for the model including MCGs, despite the fact that it has more free parameters. This suggests that if MCGs indeed play a significant role in the Cosmic Dawn, it should be possible to infer their properties from upcoming 21-cm power spectra. Our study illustrates how these observations can discriminate among uncertain galaxy formation models with varying complexities, by maximizing the Bayesian evidence., Comment: 12 pages, 3 figures, submitted to MNRAS, comments are welcome

The Epoch of Reionization (EoR) depends on the complex astrophysics governing the birth and evolution of the first galaxies and structures in the intergalactic medium. EoR models rely on cosmic microwave background (CMB) observations, and in particular the large-scale E-mode polarization power spectra (EE PS), to help constrain their highly uncertain parameters. However, rather than directly forward-modelling the EE PS, most EoR models are constrained using a summary statistic -- the Thompson scattering optical depth, $\tau_e$. Compressing CMB observations to $\tau_e$ requires adopting a basis set for the EoR history. The common choice is the unphysical, redshift-symmetric hyperbolic tangent (Tanh) function, which differs in shape from physical EoR models based on hierarchical structure formation. Combining public EoR and CMB codes, 21cmFAST and CLASS, here we quantify how inference using the $\tau_e$ summary statistic impacts the resulting constraints on galaxy properties and EoR histories. Using the last Planck 2018 data release, we show that the marginalized constraints on the EoR history are more sensitive to the choice of the basis set (Tanh vs physical model) than to the CMB likelihood statistic ($\tau_e$ vs PS). For example, EoR histories implied by the growth of structure show a small tail of partial reionization extending to higher redshifts. However, biases in inference using $\tau_e$ are negligible for the Planck 2018 data. Using EoR constraints from high-redshift observations including the quasar dark fraction, galaxy UV luminosity functions and CMB EE PS, our physical model recovers $\tau_e=0.0569^{+0.0081}_{-0.0066}$., Comment: 11 pages, 5 figures, Updated to match the published version. Added Gamma evolution; Only minor changes