[en] As observations of the Epoch of Reionization (EoR) in redshifted 21 cm emission begin, we assess the accuracy of the early catalog results from the Precision Array for Probing the Epoch of Reionization (PAPER) and the Murchison Wide-field Array (MWA). The MWA EoR approach derives much of its sensitivity from subtracting foregrounds to <1% precision, while the PAPER approach relies on the stability and symmetry of the primary beam. Both require an accurate flux calibration to set the amplitude of the measured power spectrum. The two instruments are very similar in resolution, sensitivity, sky coverage, and spectral range and have produced catalogs from nearly contemporaneous data. We use a Bayesian Markov Chain Monte Carlo fitting method to estimate that the two instruments are on the same flux scale to within 20% and find that the images are mostly in good agreement. We then investigate the source of the errors by comparing two overlapping MWA facets where we find that the differences are primarily related to an inaccurate model of the primary beam but also correlated errors in bright sources due to CLEAN. We conclude with suggestions for mitigating and better characterizing these effects.

[en] Subtraction of astrophysical foreground contamination from 'dirty' sky maps produced by simulated measurements of the Murchison Widefield Array (MWA) has been performed by fitting a third-order polynomial along the spectral dimension of each pixel in the data cubes. The simulations are the first to include the unavoidable instrumental effects of the frequency-dependent primary antenna beams and synthesized array beams. They recover the one-dimensional spherically binned input redshifted 21 cm power spectrum within ∼1% over the scales probed most sensitively by the MWA (0.01 ∼< k ∼< 1 Mpc^-1) and demonstrate that realistic instrumental effects will not mask the epoch of reionization signal. We find that the weighting function used to produce the dirty sky maps from the gridded visibility measurements is important to the success of the technique. Uniform weighting of the visibility measurements produces the best results, whereas natural weighting significantly worsens the foreground subtraction by coupling structure in the density of the visibility measurements to spectral structure in the dirty sky map data cube. The extremely dense uv-coverage of the MWA was found to be advantageous for this technique and produced very good results on scales corresponding to |u| ∼< 500λ in the uv-plane without any selective editing of the uv-coverage.

[en] The response of the antenna is a source of uncertainty in measurements with the Experiment to Detect the Global Epoch of Reionization Signature (EDGES). We aim to validate the electromagnetic beam model of the low-band (50–100 MHz) dipole antenna with comparisons between models and against data. We find that simulations of a simplified model of the antenna over an infinite perfectly conducting ground plane are, with one exception, robust to changes in the numerical electromagnetic solver code or algorithm. For simulations of the antenna with the actual finite ground plane and realistic soil properties, we find that two out of three numerical solvers agree well. Applying our analysis pipeline to a simulated drift-scan observation from an early EDGES low-band instrument that had a 10 m × 10 m ground plane, we find residual levels after fitting and removing a five-term foreground model from the simulated data binned in local sidereal time (LST) average about 250 mK with ±40 mK variation between numerical solvers. A similar analysis of the primary 30 m × 30 m sawtooth ground plane reduced the LST-averaged residuals to about 90 mK with ±10 mK between the two viable solvers. More broadly we show that larger ground planes generally perform better than smaller ground planes. Simulated data have a power that is within 4% of real observations, a limitation of net accuracy of the sky and beam models. We observe that residual spectral structures after foreground model fits match qualitatively between simulated data and observations, suggesting that the frequency dependence of the beam is reasonably represented by the models. We find that a soil conductivity of 0.02 S m⁻¹ and relative permittivity of 3.5 yield good agreement between simulated spectra and observations. This is consistent with the soil properties reported by Sutinjo et al. for the Murchison Radio-astronomy Observatory, where EDGES is located.

[en] We report constraints on the global 21 cm signal due to neutral hydrogen at redshifts $14.8⩾<!--\; ???\; -->z⩾<!--\; ???\; -->6.5$. We derive our constraints from low-foreground observations of the average sky brightness spectrum conducted with the EDGES High-band instrument between 2015 September 7 and October 26. Observations were calibrated by accounting for the effects of antenna beam chromaticity, antenna and ground losses, signal reflections, and receiver parameters. We evaluate the consistency between the spectrum and phenomenological models for the global 21 cm signal. For tanh-based representations of the ionization history during the epoch of reionization, we rule out, at $⩾<!--\; ???\; -->2\mathrm{\sigma <!--\; ??\; -->}$ significance, models with duration of up to $\mathrm{\Delta <!--\; ??\; -->}z=1$ at $z\approx <!--\; ???\; -->8.5$ and higher than $\mathrm{\Delta <!--\; ??\; -->}z=0.4$ across most of the observed redshift range under the usual assumption that the 21 cm spin temperature is much larger than the temperature of the cosmic microwave background during reionization. We also investigate a “cold” intergalactic medium (IGM) scenario that assumes perfect Lyα coupling of the 21 cm spin temperature to the temperature of the IGM, but that the latter is not heated by early stars or stellar remants. Under this assumption, we reject tanh-based reionization models of duration $\mathrm{\Delta <!--\; ??\; -->}z\lesssim <!--\; ???\; -->2$ over most of the observed redshift range. Finally, we explore and reject a broad range of Gaussian models for the 21 cm absorption feature expected in the First Light era. As an example, we reject 100 mK Gaussians with duration (full width at half maximum) $\mathrm{\Delta <!--\; ??\; -->}z⩽<!--\; ???\; -->4$ over the range $14.2⩾<!--\; ???\; -->z⩾<!--\; ???\; -->6.5$ at $⩾<!--\; ???\; -->2\mathrm{\sigma <!--\; ??\; -->}$ significance.

[en] Unaccounted for systematics from foregrounds and instruments can severely limit the sensitivity of current experiments from detecting redshifted 21 cm signals from the Epoch of Reionization (EoR). Upcoming experiments are faced with a challenge to deliver more collecting area per antenna element without degrading the data with systematics. This paper and its companions show that dishes are viable for achieving this balance using the Hydrogen Epoch of Reionization Array (HERA) as an example. Here, we specifically identify spectral systematics associated with the antenna power pattern as a significant detriment to all EoR experiments which causes the already bright foreground power to leak well beyond ideal limits and contaminate the otherwise clean EoR signal modes. A primary source of this chromaticity is reflections in the antenna-feed assembly and between structures in neighboring antennas. Using precise foreground simulations taking wide-field effects into account, we provide a generic framework to set cosmologically motivated design specifications on these reflections to prevent further EoR signal degradation. We show that HERA will not be impeded by such spectral systematics and demonstrate that even in a conservative scenario that does not perform removal of foregrounds, HERA will detect the EoR signal in line-of-sight k-modes, $k_{\parallel <!--\; ???\; -->}\gtrsim <!--\; ???\; -->0.2h$ Mpc⁻¹, with high significance. Under these conditions, all baselines in a 19-element HERA layout are capable of detecting EoR over a substantial observing window on the sky.

[en] We use measurements from the Experiment to Detect the Global EoR Signature (EDGES) to determine scale and zero-level corrections to the diffuse radio surveys by Guzmán et al. at 45 MHz and by Landecker & Wielebinski at 150 MHz. We find that the map of Guzmán et al. requires a scale correction of 1.076 ± 0.034 (2σ) and a zero-level correction of −160 ± 78 K (2σ) to best-fit the EDGES data. For the map of Landecker & Wielebinski, the scale correction is 1.112 ± 0.023 (2σ) and the zero-level correction is 0.7 ± 6.0 K (2σ). The correction uncertainties are dominated by systematic effects, of which the most significant are uncertainty in the calibration of the EDGES receivers, antenna pointing, and tropospheric and ionospheric effects. We propagate the correction uncertainties to estimate the uncertainties in the corrected maps themselves and find that the 2σ uncertainty in the map brightness temperature is in the range 3.2%–7.5% for the map of Guzmán et al. and 2.1%–9.0% for the map of Landecker & Wielebinski, with the largest percentage uncertainties occurring at high Galactic latitudes. The corrected maps could be used to improve existing diffuse low-frequency radio sky models, which are essential tools in analyses of cosmological 21 cm observations, as well as to investigate the existence of a radio monopole excess above the cosmic microwave background and known Galactic and extragalactic contributions.

[en] 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 a 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 _∥ ≲ 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.

[en] 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 and the Precision Array for Probing the EOR, HERA is a hexagonal array of large (14 m diameter) dishes with suspended dipole feeds. The dish not only determines overall sensitivity, but also 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. In this paper, we focus 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² in the optimal dish/feed configuration, implying that 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. Finally, we study the impact of these beam measurements on the distribution of foregrounds in Fourier space.

[en] The Murchison Widefield Array (MWA) has collected hundreds of hours of Epoch of Reionization (EoR) data and now faces the challenge of overcoming foreground and systematic contamination to reduce the data to a cosmological measurement. We introduce several novel analysis techniques, such as cable reflection calibration, hyper-resolution gridding kernels, diffuse foreground model subtraction, and quality control methods. Each change to the analysis pipeline is tested against a two-dimensional power spectrum figure of merit to demonstrate improvement. We incorporate the new techniques into a deep integration of 32 hours of MWA data. This data set is used to place a systematic-limited upper limit on the cosmological power spectrum of ${\mathrm{\Delta <!--\; ??\; -->}}^2⩽<!--\; ???\; -->2.7\times <!--\; ??\; -->{10}^4$ mK² at k = 0.27 h Mpc⁻¹ and z = 7.1, consistent with other published limits, and a modest improvement (factor of 1.4) over previous MWA results. From this deep analysis, we have identified a list of improvements to be made to our EoR data analysis strategies. These improvements will be implemented in the future and detailed in upcoming publications.

[en] We present the 21 cm power spectrum analysis approach of the Murchison Widefield Array Epoch of Reionization project. In this paper, we compare the outputs of multiple pipelines for the purpose of validating statistical limits cosmological hydrogen at redshifts between 6 and 12. Multiple independent data calibration and reduction pipelines are used to make power spectrum limits on a fiducial night of data. Comparing the outputs of imaging and power spectrum stages highlights differences in calibration, foreground subtraction, and power spectrum calculation. The power spectra found using these different methods span a space defined by the various tradeoffs between speed, accuracy, and systematic control. Lessons learned from comparing the pipelines range from the algorithmic to the prosaically mundane; all demonstrate the many pitfalls of neglecting reproducibility. We briefly discuss the way these different methods attempt to handle the question of evaluating a significant detection in the presence of foregrounds.