Modeling the Radio Background from the First Black Holes at Cosmic Dawn: Implications for the 21 cm Absorption Amplitude

We estimate the 21 cm Radio Background from accretion onto the first intermediate-mass Black Holes between $z\approx 30$ and $z\approx 16$. Combining potentially optimistic, but plausible, scenarios for black hole formation and growth with empirical correlations between luminosity and radio-emission observed in low-redshift active galactic nuclei, we find that a model of black holes forming in molecular cooling halos is able to produce a 21 cm background that exceeds the Cosmic Microwave Background (CMB) at $z \approx 17$ though models involving larger halo masses are not entirely excluded. Such a background could explain the surprisingly large amplitude of the 21 cm absorption feature recently reported by the EDGES collaboration. Such black holes would also produce significant X-ray emission and contribute to the $0.5-2$ keV soft X-ray background at the level of $\approx 10^{-13}-10^{-12}$ erg sec$^{-1}$ cm$^{-2}$ deg$^{-2}$, consistent with existing constraints. In order to avoid heating the IGM over the EDGES trough, these black holes would need to be obscured by Hydrogen column depths of $N_\text{H} \sim 5 \times 10^{23} \text{cm}^{-2}$. Such black holes would avoid violating contraints on the CMB optical depth from Planck if their UV photon escape fractions were below $f_{\text{esc}} \lesssim 0.1$, which would be a natural result of $N_\text{H} \sim 5 \times 10^{23} \text{cm}^{-2}$ imposed by an unheated IGM.Comment: 10 pages, 7 figures, accepted to ApJ, replacement to match submitted versio

Spectral Redundancy for Calibrating Interferometers and Suppressing the Foreground Wedge in 21\,cm Cosmology

Observations of 21\,cm line from neutral hydrogen promise to be an exciting new probe of astrophysics and cosmology during the Cosmic Dawn and through the Epoch of Reionization (EoR) to when dark energy accelerates the expansion of the Universe. At each of these epochs, separating bright foregrounds from the cosmological signal is a primary challenge that requires exquisite calibration. In this paper, we present a new calibration method called \textsc{nucal} that extends redundant-baseline calibration, allowing spectral variation in antenna responses to be solved for by using correlations between visibilities measuring the same angular Fourier modes at different frequencies. By modeling the chromaticity of the beam-weighted sky with a tunable set of discrete prolate spheroidal sequences (DPSS), we develop a calibration loop that optimizes for spectrally smooth calibrated visibilities. Crucially, this technique does not require explicit models of the sky or the primary beam. With simulations that incorporate realistic source and beam chromaticity, we show that this method solves for unsmooth bandpass features, exposes narrowband interference systematics, and suppresses smooth-spectrum foregrounds below the level of 21\,cm reionization models, even within much of the so-called "wedge" region where current foreground mitigation techniques struggle. We show that this foreground subtraction can be performed with minimal cosmological signal loss for certain well-sampled angular Fourier modes, making spectral-redundant calibration a promising technique for current and next-generation 21\,cm intensity mapping experiments.Comment: 22 pages, 10 figures, Submitted to MNRA

Modeling the Radio Background from the First Black Holes at Cosmic Dawn: Implications for the 21 cm Absorption Amplitude

We estimate the 21 cm radio background from accretion onto the first intermediate-mass black holes between z ≈ 30 and z ≈ 16. Combining potentially optimistic, but plausible, scenarios for black hole formation and growth with empirical correlations between luminosity and radio emission observed in low-redshift active galactic nuclei, we find that a model of black holes forming in molecular cooling halos is able to produce a 21 cm background that exceeds the cosmic microwave background (CMB) at z ≈ 17, though models involving larger halo masses are not entirely excluded. Such a background could explain the surprisingly large amplitude of the 21 cm absorption feature recently reported by the EDGES collaboration. Such black holes would also produce significant X-ray emission and contribute to the 0.5–2 keV soft X-ray background at the level of ≈10^(−13)–10^(−12) erg s^(−1)cm^(−2) deg^(−2), consistent with existing constraints. In order to avoid heating the intergalactic medium (IGM) over the EDGES trough, these black holes would need to be obscured by hydrogen column depths of N_H ~ 5 × 10^(23) cm^(−2). Such black holes would avoid violating constraints on the CMB optical depth from Planck if their UV photon escape fractions were below f_(esc) ≾ 0.1, which would be a natural result of N_H ~ 5 × 10^(23) cm^(−2) being imposed by an unheated IGM

DAYENU: a simple filter of smooth foregrounds for intensity mapping power spectra

We introduce DPSS Approximate lazY filtEriNg of foregroUnds (DAYENU), a linear, spectral filter for H I intensity mapping that achieves the desirable foreground mitigation and error minimization properties of inverse co-variance weighting with minimal modelling of the underlying data. Beyond 21-cm power-spectrum estimation, our filter is suitable for any analysis where high dynamic-range removal of spectrally smooth foregrounds in irregularly (or regularly) sampled data is required, something required by many other intensity mapping techniques. Our filtering matrix is diagonalized by Discrete Prolate Spheroidal Sequences which are an optimal basis to model band-limited foregrounds in 21-cm intensity mapping experiments in the sense that they maximally concentrate power within a finite region of Fourier space. We show that DAYENU enables the access of large-scale line-of-sight modes that are inaccessible to tapered discrete Fourier transform estimators. Since these modes have the largest SNRs,DAYENU significantly increases the sensitivity of 21-cm analyses over tapered Fourier transforms. Slight modifications allow us to use DAYENU as a linear replacement for iterative delay CLEAN ing (DAYENUREST). We refer readers to the Code section at the end of this paper for links to examples and code

The Hydrogen Epoch of Reionization Array Dish II: Characterization of Spectral Structure with Electromagnetic Simulations and its science Implications

We use time-domain electromagnetic simulations to determine the spectral characteristics of the Hydrogen Epoch of Reionization Arrays (HERA) antenna. These simulations are part of a multi-faceted campaign to determine the effectiveness of the dish's design for obtaining a detection of redshifted 21 cm emission from the epoch of reionization. Our simulations show the existence of reflections between HERA's suspended feed and its parabolic dish reflector that fall below -40 dB at 150 ns and, for reasonable impedance matches, have a negligible impact on HERA's ability to constrain EoR parameters. It follows that despite the reflections they introduce, dishes are effective for increasing the sensitivity of EoR experiments at relatively low cost. We find that electromagnetic resonances in the HERA feed's cylindrical skirt, which is intended to reduce cross coupling and beam ellipticity, introduces significant power at large delays ($-40$ dB at 200 ns) which can lead to some loss of measurable Fourier modes and a modest reduction in sensitivity. Even in the presence of this structure, we find that the spectral response of the antenna is sufficiently smooth for delay filtering to contain foreground emission at line-of-sight wave numbers below $k_\parallel \lesssim 0.2$ $h$Mpc$^{-1}$, in the region where the current PAPER experiment operates. Incorporating these results into a Fisher Matrix analysis, we find that the spectral structure observed in our simulations has only a small effect on the tight constraints HERA can achieve on parameters associated with the astrophysics of reionization.Comment: Accepted to ApJ, 18 pages, 17 Figures. Replacement matches accepted manuscrip

WSCLEAN: an implementation of a fast, generic wide-field imager for radio astronomy

Astronomical wide-field imaging of interferometric radio data is computationally expensive, especially for the large data volumes created by modern non-coplanar many-element arrays. We present a new wide-field interferometric imager that uses the w-stacking algorithm and can make use of the w-snapshot algorithm. The performance dependences of CASA's w-projection and our new imager are analysed and analytical functions are derived that describe the required computing cost for both imagers. On data from the Murchison Widefield Array, we find our new method to be an order of magnitude faster than w-projection, as well as being capable of full-sky imaging at full resolution and with correct polarization correction. We predict the computing costs for several other arrays and estimate that our imager is a factor of 2-12 faster, depending on the array configuration. We estimate the computing cost for imaging the lowfrequency Square Kilometre Array observations to be 60 PetaFLOPS with current techniques. We find that combining w-stacking with the w-snapshot algorithm does not significantly improve computing requirements over pure w-stacking. The source code of our new imager is publicly released

The Hydrogen Epoch of Reionization Array Dish I: Beam Pattern Measurements and Science Implications

The Hydrogen Epoch of Reionization Array (HERA) is a radio interferometer aiming to detect the power spectrum of 21 cm fluctuations from neutral hydrogen from the Epoch of Reionization (EOR). Drawing on lessons from the Murchison Widefield Array (MWA) and the Precision Array for Probing the Epoch of Reionization (PAPER), HERA is a hexagonal array of large (14 m diameter) dishes with suspended dipole feeds. Not only does the dish determine overall sensitivity, it affects the observed frequency structure of foregrounds in the interferometer. This is the first of a series of four papers characterizing the frequency and angular response of the dish with simulations and measurements. We focus in this paper on the angular response (i.e., power pattern), which sets the relative weighting between sky regions of high and low delay, and thus, apparent source frequency structure. We measure the angular response at 137 MHz using the ORBCOMM beam mapping system of Neben et al. We measure a collecting area of 93 m^2 in the optimal dish/feed configuration, implying HERA-320 should detect the EOR power spectrum at z~9 with a signal-to-noise ratio of 12.7 using a foreground avoidance approach with a single season of observations, and 74.3 using a foreground subtraction approach. Lastly we study the impact of these beam measurements on the distribution of foregrounds in Fourier space.Comment: 13 pages, 9 figures. Replaced to match accepted ApJ versio

WSClean : an implementation of a fast, generic wide-field imager for radio astronomy

This article has been accepted for publication in Monthly Notices of the Royal Astronomical Society. © 2014 The Authors. Published by Oxford University Press on behalf of the Royal Astronomical Society.Astronomical widefield imaging of interferometric radio data is computationally expensive, especially for the large data volumes created by modern non-coplanar many-element arrays. We present a new widefield interferometric imager that uses the w-stacking algorithm and can make use of the w-snapshot algorithm. The performance dependencies of CASA's w-projection and our new imager are analysed and analytical functions are derived that describe the required computing cost for both imagers. On data from the Murchison Widefield Array, we find our new method to be an order of magnitude faster than w-projection, as well as being capable of full-sky imaging at full resolution and with correct polarisation correction. We predict the computing costs for several other arrays and estimate that our imager is a factor of 2-12 faster, depending on the array configuration. We estimate the computing cost for imaging the low-frequency Square-Kilometre Array observations to be 60 PetaFLOPS with current techniques. We find that combining w-stacking with the w-snapshot algorithm does not significantly improve computing requirements over pure w-stacking. The source code of our new imager is publicly released.Peer reviewedFinal Published versio