In order to study the "Cosmic Dawn" and the Epoch of Reionization with 21 cm tomography, we need to statistically separate the cosmological signal from foregrounds known to be orders of magnitude brighter. Over the last few years, we have learned much about the role our telescopes play in creating a putatively foreground-free region called the "EoR window." In this work, we examine how an interferometer's effects can be taken into account in a way that allows for the rigorous estimation of 21 cm power spectra from interferometric maps while mitigating foreground contamination and thus increasing sensitivity. This requires a precise understanding of the statistical relationship between the maps we make and the underlying true sky. While some of these calculations would be computationally infeasible if performed exactly, we explore several well-controlled approximations that make mapmaking and the calculation of map statistics much faster, especially for compact and highly redundant interferometers designed specifically for 21 cm cosmology. We demonstrate the utility of these methods and the parametrized trade-offs between accuracy and speed using one such telescope, the upcoming Hydrogen Epoch of Reionization Array, as a case study.

In this white paper, we lay out a US roadmap for high-redshift 21 cm cosmology (30 < z < 6) in the 2020s. Beginning with the currently-funded HERA and MWA Phase II projects and advancing through the decade with a coordinated program of small-scale instrumentation, software, and analysis projects targeting technology development, this roadmap incorporates our current best understanding of the systematics confronting 21 cm cosmology into a plan for overcoming them, enabling next-generation, mid-scale 21 cm arrays to be proposed late in the decade. Submitted for consideration by the Astro2020 Decadal Survey Program Panel for Radio, Millimeter, and Submillimeter Observations from the Ground as a Medium-Sized Project.

Radio frequency interference (RFI) is an ever-present limiting factor among radio telescopes even in the most remote observing locations. When looking to retain the maximum amount of sensitivity and reduce contamination for Epoch of Reionization studies, the identification and removal of RFI is especially important. In addition to improved RFI identification, we must also take into account computational efficiency of the RFI-Identification algorithm as radio interferometer arrays such as the Hydrogen Epoch of Reionization Array (HERA) grow larger in number of receivers. To address this, we present a deep fully convolutional neural network (DFCN) that is comprehensive in its use of interferometric data, where both amplitude and phase information are used jointly for identifying RFI. We train the network using simulated HERA visibilities containing mock RFI, yielding a known ‘ground truth’ data set for evaluating the accuracy of various RFI algorithms. Evaluation of the DFCN model is performed on observations from the 67 dish build-out, HERA-67, and achieves a data throughput of 1.6 × 105 HERA time-ordered 1024 channelled visibilities per hour per GPU. We determine that relative to an amplitude only network including visibility phase adds important adjacent time–frequency context which increases discrimination between RFI and non-RFI. The inclusion of phase when predicting achieves a recall of 0.81, precision of 0.58, and F2 score of 0.75 as applied to our HERA-67 observations.

Spectral structures due to the instrument response is the current limiting factor for the experiments attempting to detect the redshifted 21 cm signal from the Epoch of Reionization (EoR). Recent advances in the delay spectrum methodology for measuring the redshifted 21 cm EoR power spectrum brought new attention to the impact of an antenna’s frequency response on the viability of making this challenging measurement. The delay spectrum methodology provides a somewhat straightforward relationship between the time-domain response of an instrument that can be directly measured and the power spectrum modes accessible to a 21 cm EoR experiment. In this paper, we derive the explicit relationship between antenna reflection coefficient (S11) measurements made by a Vector Network Analyzer (VNA) and the extent of additional foreground contaminations in delay space. In the light of this mathematical framework, we examine the chromaticity of a prototype antenna element that will constitute the Hydrogen Epoch of Reionization Array (HERA) between 100 and 200 MHz. These reflectometry measurements exhibit additional structures relative to electromagnetic simulations, but we find that even without any further design improvement, such an antenna element will support measuring spatial k modes with line-of-sight components of k∥ > 0.2h Mpc− 1. We also find that when combined with the powerful inverse covariance weighting method used in optimal quadratic estimation of redshifted 21 cm power spectra the HERA prototype elements can successfully measure the power spectrum at spatial modes as low as k∥ > 0.1h Mpc− 1. This work represents a major step toward understanding the HERA antenna element and highlights a straightforward method for characterizing instrument response for future experiments designed to detect the 21 cm EoR power spectrum.

Development of the hardware, data analysis, and simulation techniques for large compact radio arrays dedicated to mapping the 21 cm line of neutral hydrogen gas has proven to be more difficult than imagined twenty years ago when such telescopes were first proposed. Despite tremendous technical and methodological advances, there are several outstanding questions on how to optimally calibrate and analyze such data. On the positive side, it has become clear that the outstanding issues are purely technical in nature and can be solved with sufficient development activity. Such activity will enable science across redshifts, from early galaxy evolution in the pre-reionization era to dark energy evolution at low redshift.