Interferometric experiments designed to detect the highly redshifted 21-cm signal from neutral hydrogen are producing increasingly stringent constraints on the 21-cm power spectrum, but some k-modes remain systematics-dominated. Mutual coupling is a major systematic that must be overcome in order to detect the 21-cm signal, and simulations that reproduce effects seen in the data can guide strategies for mitigating mutual coupling. In this paper, we analyse 12 nights of data from the Hydrogen Epoch of Reionization Array and compare the data against simulations that include a computationally efficient and physically motivated semi-analytic treatment of mutual coupling. We find that simulated coupling features qualitatively agree with coupling features in the data; however, coupling features in the data are brighter than the simulated features, indicating the presence of additional coupling mechanisms not captured by our model. We explore the use of fringe-rate filters as mutual coupling mitigation tools and use our simulations to investigate the effects of mutual coupling on a simulated cosmological 21-cm power spectrum in a "worst case" scenario where the foregrounds are particularly bright. We find that mutual coupling contaminates a large portion of the "EoR Window", and the contamination is several orders-of-magnitude larger than our simulated cosmic signal across a wide range of cosmological Fourier modes. While our fiducial fringe-rate filtering strategy reduces mutual coupling by roughly a factor of 100 in power, a non-negligible amount of coupling cannot be excised with fringe-rate filters, so more sophisticated mitigation strategies are required.
The 21 cm transition from neutral Hydrogen promises to be the best observational probe of the Epoch of Reionization (EoR). This has led to the construction of low-frequency radio interferometric arrays, such as the Hydrogen Epoch of Reionization Array (HERA), aimed at systematically mapping this emission for the first time. Precision calibration, however, is a requirement in 21 cm radio observations. Due to the spatial compactness of HERA, the array is prone to the effects of mutual coupling, which inevitably lead to non-smooth calibration errors that contaminate the data. When unsmooth gains are used in calibration, intrinsically spectrally smooth foreground emission begins to contaminate the data in a way that can prohibit a clean detection of the cosmological EoR signal. In this paper, we show that the effects of mutual coupling on calibration quality can be reduced by applying custom time-domain filters to the data prior to calibration. We find that more robust calibration solutions are derived when filtering in this way, which reduces the observed foreground power leakage. Specifically, we find a reduction of foreground power leakage by 2 orders of magnitude at k|| approximate to 0 . 5 h Mpc(-1).
Tomas Cassanelli1,2,3,†, Calvin Leung4,5,6,∗,†, Pranav Sanghavi7,8,9,†, Juan Mena-Parra, Savannah Cary, Ryan Mckinven, Mohit Bhardwaj, Kiyoshi W. Masui, Daniele Michilli, Kevin Bandura, Shami Chatterjee, Jeffrey B. Peterson, Jane Kaczmarek, Chitrang Patel, Mubdi Rahman, Kaitlyn Shin, Keith Vanderlinde, Sabrina Berger, Charanjot Brar, P. J. Boyle, Daniela Breitman, Pragya Chawla, Alice P. Curtin, Matt Dobbs, Fengqiu Adam Dong, Emmanuel Fonseca, B. M. Gaensler, Adaeze Ibik, Victoria M. Kaspi, Khairy Kholoud, T. L. Landecker, Adam E. Lanman, Mattias Lazda, Hsiu-Hsien Lin, Jing Luo, Bradley W. Meyers, Nikola Milutinovic, Cherry Ng, Gavin Noble, Aaron B. Pearlman, Ue-Li Pen, Emily Petroff, Ziggy Pleunis, Brendan Quine, Masoud Rafiei-Ravandi, Andre Renard, Ketan R. Sand, Eve Schoen, Paul Scholz, Kendrick M. Smith, Ingrid Stairs, Shriharsh P. Tendulkar
Combining the visibilities measured by an interferometer to form a cosmological power spectrum is a complicated process in which the window functions play a crucial role. Indeed, in a delay-based analysis, the mapping between instrumental space, made of per-baseline delay spectra, 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 in this paper the impact of instrument characteristics and analysis choices on the power spectrum estimator by deriving its exact window functions outside of the delay approximation. Focusing on the Hydrogen Epoch of Reionization Array (HERA) as a case study, we find that observations made with long baselines tend to correspond to enhanced low- 𝑘 tails of the window functions, which facilitate foreground leakage outside the wedge, whilst the choice of bandwidth and frequency taper can help narrowing them down. With the help of simple test cases and more realistic visibility simulations, we show that, apart from tracing mode mixing, the window functions can accurately reconstruct the power spectrum estimator of simulated visibilities. We note that the window functions depend strongly on the chromaticity of the beam, 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- 𝑘 power to avoid foreground leakage. The window functions presented in this work 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 can be used in future analyses to confront theoretical models and data directly in cylindrical space.
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.
The compact configuration of Phase II of the Murchison Widefield Array (MWA) consists of both a redundant subarray and pseudo-random baselines, offering unique opportunities to perform sky-model and redundant interferometric calibration. The highly redundant hexagonal cores give improved power spectrum sensitivity. In this paper, we present the analysis of nearly 40 hr of data targeting one of the MWA's epoch of reionization (EoR) fields observed in 2016. We use both improved analysis techniques presented in Barry et al. and several additional techniques developed for this work, including data quality control methods and interferometric calibration approaches. We show the EoR power spectrum limits at redshift 6.5, 6.8, and 7.1 based on our deep analysis on this 40 hr data set. These limits span a range in k-space of 0.18 h Mpc(-1) < k < 1.6 h Mpc(-1), with a lowest measurement of Delta(2) <= 2.39 x 10(3) mK(2) at k = 0.59 h Mpc(-1) and z = 6.5.
Measurements of 21 cm Epoch of Reionization (EoR) structure are subject to systematics originating from both the analysis and the observation conditions. Using 2013 data from the Murchison Widefield Array (MWA), we show the importance of mitigating both sources of contamination. A direct comparison between results from Beardsley et al. and our updated analysis demonstrates new precision techniques, lowering analysis systematics by a factor of 2.8 in power. We then further lower systematics by excising observations contaminated by ultra-faint RFI, reducing by an additional factor of 3.8 in power for the zenith pointing. With this enhanced analysis precision and newly developed RFI mitigation, we calculate a noise-dominated upper limit on the EoR structure of Δ 2 ≤ 3.9 × 10 3 mK 2 at k = 0.20 h Mpc −1 and z = 7 using 21 hr of data, improving previous MWA limits by almost an order of magnitude.
Interferometric arrays seeking to measure the 21 cm signal from the epoch of reionization (EOR) must contend with overwhelmingly bright emission from foreground sources. Accurate recovery of the 21 cm signal will require precise calibration of the array, and several new avenues for calibration have been pursued in recent years, including methods using redundancy in the antenna configuration. The newly upgraded Phase II of Murchison Widefield Array (MWA) is the first interferometer that has large numbers of redundant baselines while retaining good instantaneous UV coverage. This array therefore provides a unique opportunity to compare redundant calibration with sky-model-based algorithms. In this paper, we present the first results from comparing both calibration approaches with MWA Phase II observations. For redundant calibration, we use the package OMNICAL and produce sky-based calibration solutions with the analysis package Fast Holographic Deconvolution (FHD). There are three principal results: (1) We report the success of OMNICAL on observations of ORBComm satellites, showing substantial agreement between redundant visibility measurements after calibration. (2) We directly compare OMNICAL calibration solutions with those from FHD and demonstrate that these two different calibration schemes give extremely similar results. (3) We explore improved calibration by combining OMNICAL and FHD. We evaluate these combined methods using power spectrum techniques developed for EOR analysis and find evidence for marginal improvements mitigating artifacts in the power spectrum. These results are likely limited by the signal-to-noise ratio in the 6 hr of data used, but they suggest future directions for combining these two calibration schemes.
The Murchison Widefield Array is designed to measure the fluctuations in the 21cm emission from neutral hydrogen during the Epoch of Reionisation. The new hex configuration is explicitly designed to test the predicted increase in sensitivity of redundant baselines. However the challenge of the new array is to understand calibration with the new configuration. We have developed two new pipelines to reduce the hex data, and will compare the results with previous datasets from the Phase 1 array. We have now processed 80 hours of data refining the data analysis through our two established Phase 1 pipelines. This proposal requests as much observing time as possible in semester 2017-A to (1) obtain a comparable hex dataset to test the sensitivity and systematic limits with redundant arrays,(2) establish the optimal observing strategy for an EoR detection, and (3) continue to explore observational strategies in the …