Search

The locations of Ly-$\alpha$ emitting galaxies (LAEs) at the end of the Epoch of Reionisation (EoR) are expected to correlate with regions of ionised hydrogen, traced by the redshifted 21~cm hyperfine line. Mapping the neutral hydrogen around regions with detected and localised LAEs offers an avenue to constrain the brightness temperature of the Universe within the EoR by providing an expectation for the spatial distribution of the gas, thereby providing prior information unavailable to power spectrum measurements. We use a test set of 12 hours of observations from the Murchison Widefield Array (MWA) in extended array configuration, to constrain the neutral hydrogen signature of 58 LAEs, detected with the Subaru Hypersuprime Cam in the \textit{Silverrush} survey, centred on $z$=6.58. We assume that detectable emitters reside in the centre of ionised HII bubbles during the end of reionization, and predict the redshifted neutral hydrogen signal corresponding to the remaining neutral regions using a set of different ionised bubble radii. A prewhitening matched filter detector is introduced to assess detectability. We demonstrate the ability to detect, or place limits upon, the amplitude of brightness temperature fluctuations, and the characteristic HII bubble size. With our limited data, we constrain the brightness temperature of neutral hydrogen to $\Delta{\rm T}_B, Comment: 9 pages, 6 figures, accepted for publication in MNRAS

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 acknowledge the support of the Australian Renewable Energy Agency (ARENA), the Australian Centre for Advanced Photovoltaics (ACAP), and the ANFF ACT Node in carrying out this research. H.T.N. acknowledges the fellowship support of the Australian Centre for Advanced Photovoltaics. T.P.W. is the recipient of an Australian Research Council Future Fellowship (Project No. FT180100302) funded by the Australian Government.

We present a broadband map of polarized diffuse emission at 167-198 MHz developed from data from the Murchison Widefield Array (MWA). The map is designed to improve visibility simulation and precision calibration for 21 cm Epoch of Reionization (EoR) experiments. It covers a large swath - 11,000 sq. deg. - of the Southern Hemisphere sky in all four Stokes parameters and captures emission on angular scales of 1 to 9 degrees. The band-averaged diffuse structure is predominantly unpolarized but has significant linearly polarized structure near RA = 0 h. We evaluate the accuracy of the map by combining it with the GLEAM catalog and simulating an observation from the MWA, demonstrating that the accuracy of the short baselines (6.1-50 wavelengths) now approaches the accuracy of the longer baselines typically used for EoR calibration. We discuss how to use the map for visibility simulation for a variety of interferometric arrays. The map has potential to improve calibration accuracy for experiments such as the Hydrogen Epoch of Reionization Array (HERA) and the forthcoming Square Kilometre Array (SKA) as well as the MWA., Comment: Accepted for publication in MNRAS

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

The key challenge in the observation of the redshifted 21-cm signal from cosmic reionization is its separation from the much brighter foreground emission. Such separation relies on the different spectral properties of the two components, although, in real life, the foreground intrinsic spectrum is often corrupted by the instrumental response, inducing systematic effects that can further jeopardize the measurement of the 21-cm signal. In this paper, we use Gaussian Process Regression to model both foreground emission and instrumental systematics in $\sim 2$ hours of data from the Hydrogen Epoch of Reionization Array. We find that a simple co-variance model with three components matches the data well, giving a residual power spectrum with white noise properties. These consist of an "intrinsic" and instrumentally corrupted component with a coherence-scale of 20 MHz and 2.4 MHz respectively (dominating the line of sight power spectrum over scales $k_{\parallel} \le 0.2$ h cMpc$^{-1}$) and a baseline dependent periodic signal with a period of $\sim 1$ MHz (dominating over $k_{\parallel} \sim 0.4 - 0.8$h cMpc$^{-1}$) which should be distinguishable from the 21-cm EoR signal whose typical coherence-scales is $\sim 0.8$ MHz., 15 pages, 15 figures, 1 table, Accepted to MNRAS

Observations in the lowest MWA band between $75-100$ MHz have the potential to constrain the distribution of neutral hydrogen in the intergalactic medium at redshift $\sim 13-17$. Using 15 hours of MWA data, we analyse systematics in this band such as radio-frequency interference (RFI), ionospheric and wide field effects. By updating the position of point sources, we mitigate the direction independent calibration error due to ionospheric offsets. Our calibration strategy is optimized for the lowest frequency bands by reducing the number of direction dependent calibrators and taking into account radio sources within a wider field of view. We remove data polluted by systematics based on the RFI occupancy and ionospheric conditions, finally selecting 5.5 hours of the cleanest data. Using these data, we obtain two sigma upper limits on the 21 cm power spectrum in the range of $0.1\lessapprox k \lessapprox 1 ~\rm ~h~Mpc^{-1}$ and at $z$=14.2, 15.2 and 16.5, with the lowest limit being $6.3\times 10^6 ~\rm mK^2$ at $\rm k=0.14 \rm ~h~Mpc^{-1}$ and at $z=15.2$ with a possibility of a few \% of signal loss due to direction independent calibration., Comment: 17 pages, 18 figures, accepted to MNRAS

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

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

In 21 cm cosmology, precision calibration is key to the separation of the neutral hydrogen signal from very bright but spectrally-smooth astrophysical foregrounds. The Hydrogen Epoch of Reionization Array (HERA), an interferometer specialized for 21 cm cosmology and now under construction in South Africa, was designed to be largely calibrated using the self-consistency of repeated measurements of the same interferometric modes. This technique, known as "redundant-baseline calibration" resolves most of the internal degrees of freedom in the calibration problem. It assumes, however, on antenna elements with identical primary beams placed precisely on a redundant grid. In this work, we review the detailed implementation of the algorithms enabling redundant-baseline calibration and report results with HERA data. We quantify the effects of real-world non-redundancy and how they compare to the idealized scenario in which redundant measurements differ only in their noise realizations. Finally, we study how non-redundancy can produce spurious temporal structure in our calibration solutions--both in data and in simulations--and present strategies for mitigating that structure., Comment: 24 Pages, 19 Figures. Updated to match the accepted MNRAS version

In this paper, we describe our efforts towards the development of a real-time radio imaging correlator for the Long-Wavelength Array station in Sevilleta, New Mexico. We briefly discuss the direct-imaging algorithm and present the architecture of the GPU implementation. We describe the code-level modifications carried out for one of the modules in the algorithm that improves GPU-memory management and highlight the performance improvements achieved through it. We emphasize our ongoing efforts in tuning the overall run-time duration of the correlator which in turn is expected to increase the operating bandwidth in order to address the demands of wide-band capability for radio transient science.

The experimental efforts to detect the redshifted 21 cm signal from the Epoch of Reionization (EoR) are limited predominantly by the chromatic instrumental systematic effect. The delay spectrum methodology for 21 cm power spectrum measurements brought new attention to the critical impact of an antenna's chromaticity on the viability of making this measurement. This methodology established a straightforward relationship between time-domain response of an instrument and the power spectrum modes accessible to a 21 cm EoR experiment. We examine the performance of a prototype of the Hydrogen Epoch of Reionization Array (HERA) array element that is currently observing in Karoo desert, South Africa. We present a mathematical framework to derive the beam integrated frequency response of a HERA prototype element in reception from the return loss measurements between 100-200 MHz and determined the extent of additional foreground contamination in the delay space. The measurement reveals excess spectral structures in comparison to the simulation studies of the HERA element. Combined with the HERA data analysis pipeline that incorporates inverse covariance weighting in optimal quadratic estimation of power spectrum, we find that in spite of its departure from the simulated response, HERA prototype element satisfies the necessary criteria posed by the foreground attenuation limits and potentially can measure the power spectrum at spatial modes as low as $k_{\parallel} > 0.1h$~Mpc$^{-1}$. The work highlights a straightforward method for directly measuring an instrument response and assessing its impact on 21 cm EoR power spectrum measurements for future experiments that will use reflector-type antenna., 13 pages, 10 figures, Submitted to APJ

Characterizing the epoch of reionization (EoR) at $z\gtrsim 6$ via the redshifted 21 cm line of neutral Hydrogen (HI) is critical to modern astrophysics and cosmology, and thus a key science goal of many current and planned low-frequency radio telescopes. The primary challenge to detecting this signal is the overwhelmingly bright foreground emission at these frequencies, placing stringent requirements on the knowledge of the instruments and inaccuracies in analyses. Results from these experiments have largely been limited not by thermal sensitivity but by systematics, particularly caused by the inability to calibrate the instrument to high accuracy. The interferometric bispectrum phase is immune to antenna-based calibration and errors therein, and presents an independent alternative to detect the EoR HI fluctuations while largely avoiding calibration systematics. Here, we provide a demonstration of this technique on a subset of data from the Hydrogen Epoch of Reionization Array (HERA) to place approximate constraints on the brightness temperature of the intergalactic medium (IGM). From this limited data, at $z=7.7$ we infer "$1\sigma$" upper limits on the IGM brightness temperature to be $\le 316$ "pseudo" mK at $\kappa_\parallel=0.33$ "pseudo" $h$ Mpc$^{-1}$ (data-limited) and $\le 1000$ "pseudo" mK at $\kappa_\parallel=0.875$ "pseudo" $h$ Mpc$^{-1}$ (noise-limited). The "pseudo" units denote only an approximate and not an exact correspondence to the actual distance scales and brightness temperatures. By propagating models in parallel to the data analysis, we confirm that the dynamic range required to separate the cosmic HI signal from the foregrounds is similar to that in standard approaches, and the power spectrum of the bispectrum phase is still data-limited (at $\gtrsim 10^6$ dynamic range) indicating scope for further improvement in sensitivity as the array build-out continues., Comment: 22 pages, 12 figures (including sub-figures). Published in PhRvD. Abstract may be slightly abridged compared to the actual manuscript due to length limitations on arXiv

Precise instrumental calibration is of crucial importance to 21-cm cosmology experiments. The Murchison Widefield Array's (MWA) Phase II compact configuration offers us opportunities for both redundant calibration and sky-based calibration algorithms; using the two in tandem is a potential approach to mitigate calibration errors caused by inaccurate sky models. The MWA Epoch of Reionization (EoR) experiment targets three patches of the sky (dubbed EoR0, EoR1, and EoR2) with deep observations. Previous work in \cite{Li_2018} and \cite{Wenyang_2019} studied the effect of tandem calibration on the EoR0 field and found that it yielded no significant improvement in the power spectrum over sky-based calibration alone. In this work, we apply similar techniques to the EoR1 field and find a distinct result: the improvements in the power spectrum from tandem calibration are significant. To understand this result, we analyze both the calibration solutions themselves and the effects on the power spectrum over three nights of EoR1 observations. We conclude that the presence of the bright radio galaxy Fornax A in EoR1 degrades the performance of sky-based calibration, which in turn enables redundant calibration to have a larger impact. These results suggest that redundant calibration can indeed mitigate some level of model-incompleteness error., Comment: 20 pages, 11 figures, accepted for publication in PASA

We analyze data from the Hydrogen Epoch of Reionization Array. This is the third in a series of papers on the closure phase delay-spectrum technique designed to detect the HI 21cm emission from cosmic reionization. We present the details of the data and models employed in the power spectral analysis, and discuss limitations to the process. We compare images and visibility spectra made with HERA data, to parallel quantities generated from sky models based on the GLEAM survey, incorporating the HERA telescope model. We find reasonable agreement between images made from HERA data, with those generated from the models, down to the confusion level. For the visibility spectra, there is broad agreement between model and data across the full band of $\sim 80$MHz. However, models with only GLEAM sources do not reproduce a roughly sinusoidal spectral structure at the tens of percent level seen in the observed visibility spectra on scales $\sim 10$ MHz on 29 m baselines. We find that this structure is likely due to diffuse Galactic emission, predominantly the Galactic plane, filling the far sidelobes of the antenna primary beam. We show that our current knowledge of the frequency dependence of the diffuse sky radio emission, and the primary beam at large zenith angles, is inadequate to provide an accurate reproduction of the diffuse structure in the models. We discuss implications due to this missing structure in the models, including calibration, and in the search for the HI 21cm signal, as well as possible mitigation techniques., Comment: 17 pages, 10 figures, accepted for the ApJ Supplement

We present a study of internal reflection and cross-coupling systematics in Phase I of the Hydrogen Epoch of Reionization Array (HERA). In a companion paper, we outlined the mathematical formalism for such systematics and presented algorithms for modeling and removing them from the data. In this work, we apply these techniques to data from HERA’s first observing season as a method demonstration. The data show evidence for systematics that, without removal, would hinder a detection of the 21 cm power spectrum for the targeted Epoch of Reionization (EoR) line-of-sight modes in the range 0.2 h -1 Mpc-1 < {k}\parallel < 0.5 h -1 Mpc-1. In particular, we find evidence for nonnegligible amounts of spectral structure in the raw autocorrelations that overlaps with the EoR window and is suggestive of complex instrumental effects. Through systematic modeling on a single night of data, we find we can recover these modes in the power spectrum down to the integrated noise floor, achieving a dynamic range in the EoR window of 106 in power (mK2 units) with respect to the bright galactic foreground signal. Future work with deeper integrations will help determine whether these systematics can continue to be mitigated down to EoR levels. For future observing seasons, HERA will have upgraded analog and digital hardware to better control these systematics in the field.

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 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" dataset 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$\times 10^{5}$ HERA time-ordered 1024 channeled 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 $F_{2}$ score of 0.75 as applied to our HERA-67 observations., Comment: 11 pages, 7 figures

The paper outlines research and development needed over the next decade to enable three-dimension mapping of the Universe using the 21 cm line of neutral hydrogen.

PUMA is a proposal for an ultra-wideband, low-resolution and transit interferometric radio telescope operating at $200-1100\,\mathrm{MHz}$. Its design is driven by six science goals which span three science themes: the physics of dark energy (measuring the expansion history and growth of the universe up to $z=6$), the physics of inflation (constraining primordial non-Gaussianity and primordial features) and the transient radio sky (detecting one million fast radio bursts and following up SKA-discovered pulsars). We propose two array configurations composed of hexagonally close-packed 6m dish arrangements with 50% fill factor. The initial 5,000 element 'petite array' is scientifically compelling, and can act as a demonstrator and a stepping stone to the full 32,000 element 'full array'. Viewed as a 21cm intensity mapping telescope, the program has the noise equivalent of a traditional spectroscopic galaxy survey comprised of 0.6 and 2.5 billion galaxies at a comoving wavenumber of $k=0.5\,h\mathrm{Mpc}^{-1}$ spanning the redshift range $z = 0.3 - 6$ for the petite and full configurations, respectively. At redshifts beyond $z=2$, the 21cm technique is a uniquely powerful way of mapping the universe, while the low-redshift range will allow for numerous cross-correlations with existing and upcoming surveys. This program is enabled by the development of ultra-wideband radio feeds, cost-effective dish construction methods, commodity radio-frequency electronics driven by the telecommunication industry and the emergence of sufficient computing power to facilitate real-time signal processing that exploits the full potential of massive radio arrays. The project has an estimated construction cost of 55 and 330 million FY19 USD for the petite and full array configurations. Including R&D, design, operations and science analysis, the cost rises to 125 and 600 million FY19 USD, respectively.

The central challenge in 21~cm cosmology is isolating the cosmological signal from bright foregrounds. Many separation techniques rely on the accurate knowledge of the sky and the instrumental response, including the antenna primary beam. For drift-scan telescopes such as the Hydrogen Epoch of Reionization Array \citep[HERA, ][]{DeBoer2017} that do not move, primary beam characterization is particularly challenging because standard beam-calibration routines do not apply \citep{Cornwell2005} and current techniques require accurate source catalogs at the telescope resolution. We present an extension of the method from \citet{Pober2012} where they use beam symmetries to create a network of overlapping source tracks that break the degeneracy between source flux density and beam response and allow their simultaneous estimation. We fit the beam response of our instrument using early HERA observations and find that our results agree well with electromagnetic simulations down to a -20~dB level in power relative to peak gain for sources with high signal-to-noise ratio. In addition, we construct a source catalog with 90 sources down to a flux density of 1.4~Jy at 151~MHz., 22 pages, 22 figures, accepted for publication in ApJ

The detection of the Epoch of Reionization (EoR) delay power spectrum using a "foreground avoidance method" highly depends on the instrument chromaticity. The systematic effects induced by the radio-telescope spread the foreground signal in the delay domain, which contaminates the EoR window theoretically observable. Applied to the Hydrogen Epoch of Reionization Array (HERA), this paper combines detailed electromagnetic and electrical simulations in order to model the chromatic effects of the instrument, and quantify its frequency and time responses. In particular, the effects of the analogue receiver, transmission cables, and mutual coupling are included. These simulations are able to accurately predict the intensity of the reflections occurring in the 150-m cable which links the antenna to the back-end. They also show that electromagnetic waves can propagate from one dish to another one through large sections of the array due to mutual coupling. The simulated system time response is attenuated by a factor $10^{4}$ after a characteristic delay which depends on the size of the array and on the antenna position. Ultimately, the system response is attenuated by a factor $10^{5}$ after 1400 ns because of the reflections in the cable, which corresponds to characterizable ${k_\parallel}$-modes above 0.7 $h\;\rm{Mpc}^{-1}$ at 150 MHz. Thus, this new study shows that the detection of the EoR signal with HERA Phase I will be more challenging than expected. On the other hand, it improves our understanding of the telescope, which is essential to mitigate the instrument chromaticity., Comment: 13 pages, 17 figures - Submitted to MNRAS - 2nd revision

The Murchison Widefield Array (MWA) is an open access telescope dedicated to studying the low frequency (80$-$300 MHz) southern sky. Since beginning operations in mid 2013, the MWA has opened a new observational window in the southern hemisphere enabling many science areas. The driving science objectives of the original design were to observe 21\,cm radiation from the Epoch of Reionisation (EoR), explore the radio time domain, perform Galactic and extragalactic surveys, and monitor solar, heliospheric, and ionospheric phenomena. All together 60$+$ programs recorded 20,000 hours producing 146 papers to date. In 2016 the telescope underwent a major upgrade resulting in alternating compact and extended configurations. Other upgrades, including digital back-ends and a rapid-response triggering system, have been developed since the original array was commissioned. In this paper we review the major results from the prior operation of the MWA, and then discuss the new science paths enabled by the improved capabilities. We group these science opportunities by the four original science themes, but also include ideas for directions outside these categories., Comment: 38 pages, 12 figures, accepted for publication in PASA

We present limits on the 21cm power spectrum from the Epoch of Reionization (EoR) using data from the 64 antenna configuration of the Donald C. Backer Precision Array for Probing the Epoch of Reionization (PAPER) analyzed through a power spectrum pipeline independent from previous PAPER analyses. Previously reported results from PAPER have been found to contain significant signal loss (Cheng et al. 2018, arxiv:1810.05175). Several lossy steps from previous PAPER pipelines have not been included in this analysis, namely: delay-based foreground filtering, optimal fringe-rate filtering, and empirical covariance-based estimators. Steps which remain in common with previous analyses include redundant calibration and local sidereal time (LST) binning. The power spectra reported here are effectively the result of applying a linear Fourier transform analysis to the calibrated, LST binned data. This analysis also uses more data than previous publications, including the complete available redshift range of $z \sim 7.5$ to $11$. In previous PAPER analyses, many power spectrum measurements were found to be detections of noncosmological power at levels of significance ranging from two to hundreds of times the theoretical noise. Here, excess power is examined using redundancy between baselines and power spectrum jackknives. The upper limits we find on the 21cm power spectrum from reionization are ($1500$ mK)$^{2}$, ($1900$ mK)$^{2}$, ($280$ mK)$^{2}$, ($200$ mK)$^{2}$, ($380$ mK)$^{2}$, ($300$ mK)$^{2}$ at redshifts $z=10.87,\ 9.93,\ 8.68,\ 8.37,\ 8.13,$ and $7.48$, respectively. For reasons described in Cheng et al. 2018 (arxiv:1810.05175), these limits supersede all previous PAPER results (Ali et al. 2018, arxiv:1502.06016)., Comment: 28 Pages, 17 Pages, Accepted to APJ

Foreground power dominates the measurements of interferometers that seek a statistical detection of highly-redshifted HI emission from the Epoch of Reionization (EoR). The chromaticity of the instrument creates a boundary in the Fourier transform of frequency (proportional to $k_\parallel$) between spectrally smooth emission, characteristic of the strong synchrotron foreground (the "wedge"), and the spectrally structured emission from HI in the EoR (the "EoR window"). Faraday rotation can inject spectral structure into otherwise smooth polarized foreground emission, which through instrument effects or miscalibration could possibly pollute the EoR window. Using data from the HERA 19-element commissioning array, we investigate the polarization response of this new instrument in the power spectrum domain. We perform a simple image-based calibration based on the unpolarized diffuse emission of the Global Sky Model, and show that it achieves qualitative redundancy between the nominally-redundant baselines of the array and reasonable amplitude accuracy. We construct power spectra of all fully polarized coherencies in all pseudo-Stokes parameters. We compare to simulations based on an unpolarized diffuse sky model and detailed electromagnetic simulations of the dish and feed, confirming that in Stokes I, the calibration does not add significant spectral structure beyond the expected level. Further, this calibration is stable over the 8 days of observations considered. Excess power is seen in the power spectra of the linear polarization Stokes parameters which is not easily attributable to leakage via the primary beam, and results from some combination of residual calibration errors and actual polarized emission. Stokes V is found to be highly discrepant from the expectation of zero power, strongly pointing to the need for more accurate polarized calibration., Comment: 21 pages, 12 figures, submitted to ApJ

Interferometric arrays seeking to measure the 21 cm signal from the Epoch of Reionization 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 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 signal-to-noise in the six hours of data used, but suggest future directions for combining these two calibration schemes., 20 pages, 11 figures. Accepted to ApJ

Roughly 40 billion M dwarfs in our galaxy host at least one small planet in the habitable zone (HZ). The stellar ultraviolet (UV) radiation from M dwarfs is strong and highly variable, and impacts planetary atmospheric loss, composition and habitability. These effects are amplified by the extreme proximity of their HZs (0.1–0.4 AU). Knowing the UV environments of M dwarf planets will be crucial to understanding their atmospheric composition and a key parameter in discriminating between biological and abiotic sources for observed biosignatures. The Star-Planet Activity Research CubeSat (SPARCS) will be a 6U CubeSat devoted to photometric monitoring of M stars in the far-UV and near-UV, measuring the time-dependent spectral slope, intensity and evolution of low-mass star high-energy radiation.

Author(s): Ali, ZS; Parsons, AR; Zheng, H; Pober, JC; Liu, A; Aguirre, JE; Bradley, RF; Bernardi, G; Carilli, CL; Cheng, C; Deboer, DR; Dexter, MR; Grobbelaar, J; Horrell, J; Jacobs, DC; Klima, P; Macmahon, DHE; Maree, M; Moore, DF; Razavi, N; Stefan, II; Walbrugh, WP; Walker, A | Abstract: In this erratum, we retract the upper limits on the 21 cm power spectrum presented in the published article. The published article reported an upper limit on δ 212(k ) of (22.4 mK)2 at z= 8.4 in the range 0.15 l k l 0.5h Mpc-1. This analysis underestimated the level of signal loss, or attenuation of the target cosmological 21 cm signal associated with the chosen power spectrum estimator, and also underestimated the statistical error on those estimates. A revised result, with a new analysis, is presented in M. K. Kolopanis et al. (2018, in preparation). Below, we briefly summarize the errors in the original analysis and how they are corrected. For an indepth analysis and discussion of the errors, we refer the reader to Cheng et al. (2018). Signal loss was expected in the original analysis because the covariance matrices, C, used to weight the un-normalized bandpower estimates, qa, in (Formula Presented) were empirically estimated from a time-averaged finite ensemble of the data, x, such that (Formula Presented). While the true covariance C leads to an inherently unbiased lossless estimator of the power spectrum, using an empirically estimated Ĉ can lead to signal loss. Specifically, weighting data by an empirically estimated covariance carries the risk of overfitting and downweighting EoR fluctuations that are coupled to the data. In Cheng et al. (2018), it is shown that these couplings are especially strong in the fringe-rate filtered PAPER-64 data set. The first and most impactful error relates to the method by which signal loss was estimated. To assess signal loss from the empirically estimated covariance matrix, different realizations of mock cosmological signals e of known amplitudes are added to the original data to form a new data vector, (Formula Presented). New covariance matrices, (Formula Presented) are used to estimate un-normalized bandpowers (Formula Presented), which can be written as (Formula Presented). The normalized power estimate can then be compared to the known injected power in e to estimate signal loss. The key error in the previous analysis was to assume that, since e was statistically independent of x, that the two middle crossterms in Equation (2) would average to zero in an ensemble. However, as shown in Cheng et al. (2018) and Switzer et al. (2015), these cross-terms can contain significant negative power because Ĉr contains information that correlates the two vectors. Ignoring these cross-terms leads to a significant underestimate of signal loss. As a result, we presented negligible signal loss in our original analysis, when in fact approximately 99.99% of the signal was removed (Cheng et al. 2018). Correcting for the actual signal loss is the biggest factor revising the upper limit on 21 D2 . The second mistake made in the original analysis was to underestimate the statistical errors in the reported power spectrum estimates. The original analysis used a bootstrap resampling technique on power spectral measurements over the baseline and time axes. However, fringe-rate filtering introduces significant correlations in the data along the time axis. As is discussed in Cheng et al. (2018), bootstrapping across correlated samples can result in a significant underestimate of the variation in the data if the number of resamplings is not equal to the number of independent samples in the data, as in the case of the original analysis. The error bars associated with this oversampling were underestimated by approximately a factor of 2 (in mK). The revised analysis in M. K. Kolopanis et al. (2018, in preparation) only applies bootstrap resampling across the baseline axis to avoid this problem. The mistake in estimating the statistical errors should have become apparent when comparing results to our theoretical thermal noise sensitivity. Unfortunately, a third miscalculation was made in estimating the thermal noise sensitivity. As detailed in Cheng et al. (2018), this miscalculation stemmed from numerous small mismatches between the idealized analysis pipeline used to estimate sensitivity and the actual analysis applied to the data. As a result, our estimated thermal noise sensitivity was approximately a factor of 3 low (in mK), leading to the mistaken impression that our error bars were consistent with the level of thermal noise. In summary, we retract the power spectrum results shown in Figures 18 and 20 in the published article. Results that relied on the original limits, including those presented in Figure 21, are retracted. Additionally, the companion paper to the original manuscript, Pober et al. (2015), used the original limits to place constraints on the spin temperature of the intergalactic medium (IGM) at z=8.4. Our revised limits do not place significant constraints on the IGM temperature, and the results of Figure 4 from Pober et al. (2015) should be disregarded. However, we note that their analysis would still be relevant should a future experiment place constraints on the 21 cm signal similar to those claimed in the published article. An updated analysis of this same data set is presented in M. K. Kolopanis et al. (2018, in preparation), where these revised results are put into context with measurements at other redshifts.

Observations of the 21cm Epoch of Reionization (EoR) signal are dominated by Galactic and extragalactic foregrounds. The need for foreground removal has led to the development of two main techniques, often referred to as "foreground avoidance" and "foreground subtraction." Avoidance is associated with filtering foregrounds in Fourier space, while subtraction uses an explicit foreground model that is removed. Using 1088 hours of data from the 64-element PAPER array, we demonstrate that subtraction of a foreground model prior to delay-space foreground filtering results in a modest but measurable improvement of the performance of the filter. This proof-of-concept result shows that improvement stems from the reduced dynamic range requirements needed for the foreground filter: subtraction of a foreground model reduces the total foreground power, so for a fixed dynamic range, the filter can push towards fainter limits. We also find that the choice of window function used in the foreground filter can have an appreciable affect on the performance near the edges of the observing band. We demonstrate these effects using a smaller 3 hour sampling of data from the MWA, and find that the hybrid filtering and subtraction removal approach provides similar improvements across the band as seen in the case with PAPER-64., Comment: 20 pages, 12 figures

Structure imprinted in foreground extragalactic point sources by ionospheric refraction has the potential to contaminate Epoch of Reionisation (EoR) power spectra of the 21~cm emission line of neutral hydrogen. The alteration of the spatial and spectral structure of foreground measurements due to total electron content (TEC) gradients in the ionosphere create a departure from the expected sky signal. We present a general framework for understanding the signatures of ionospheric behaviour in the two-dimensional (2D) neutral hydrogen power spectrum measured by a low-frequency radio interferometer. Two primary classes of ionospheric behaviour are considered, corresponding to dominant modes observed in Murchison Widefield Array (MWA) EoR data; namely, anisotropic structured wave behaviour, and isotropic turbulence. Analytic predictions for power spectrum bias due to this contamination are computed, and compared with simulations. We then apply the ionospheric metric described in Jordan et al. (2017) to study the impact of ionospheric structure on MWA data, by dividing MWA EoR datasets into classes with good and poor ionospheric conditions, using sets of matched 30-minute observations from 2014 September. The results are compared with the analytic and simulated predictions, demonstrating the observed bias in the power spectrum when the ionosphere is active (displays coherent structures or isotropic turbulence). The analysis demonstrates that unless ionospheric activity can be quantified and corrected, active data should not be included in EoR analysis in order to avoid systematic biases in cosmological power spectra. When data are corrected with a model formed from the calibration information, bias reduces below the expected 21~cm signal level. Data are considered `quiet' when the median measured source position offsets are less than 10-15~arcseconds., Comment: 19 pages, 11 figures, Accepted for publication in the Astrophysical Journal

We present first results from radio observations with the Murchison Widefield Array seeking to constrain the power spectrum of 21 cm brightness temperature fluctuations between the redshifts of 11.6 and 17.9 (113 and 75 MHz). 3 h of observations were conducted over two nights with significantly different levels of ionospheric activity. We use these data to assess the impact of systematic errors at low frequency, including the ionosphere and radio-frequency interference, on a power spectrum measurement. We find that after the 1–3 h of integration presented here, our measurements at the Murchison Radio Observatory are not limited by RFI, even within the FM band, and that the ionosphere does not appear to affect the level of power in the modes that we expect to be sensitive to cosmology. Power spectrum detections, inconsistent with noise, due to fine spectral structure imprinted on the foregrounds by reflections in the signal-chain, occupy the spatial Fourier modes where we would otherwise be most sensitive to the cosmological signal. We are able to reduce this contamination using calibration solutions derived from autocorrelations so that we achieve an sensitivity of 104 mK on comoving scales k ≲ 0.5 h Mpc[superscript −1]. This represents the first upper limits on the 21 cm power spectrum fluctuations at redshifts 12 ≲ z ≲ 18 but is still limited by calibration systematics. While calibration improvements may allow us to further remove this contamination, our results emphasize that future experiments should consider carefully the existence of and their ability to calibrate out any spectral structure within the EoR window., National Science Foundation (U.S.) (AST-0457585), National Science Foundation (U.S.) (AST-0821321), National Science Foundation (U.S.) (AST-1105835), National Science Foundation (U.S.) (AST-1410719), National Science Foundation (U.S.) (AST-1410484), National Science Foundation (U.S.) (AST-1411622), National Science Foundation (U.S.) (AST-1440343), National Science Foundation (U.S.) (1122374), National Science Foundation (U.S.) (PHY-0835713), National Science Foundation (U.S.) (CAREER-0847753), National Science Foundation (U.S.) (AST-0908884), United States. Air Force Office of Scientific Research (FA9550-0510247)

We present the first Murchison Widefield Array observations of the well-known cluster of galaxies Abell 3667 (A3667) between 105 and 241 MHz. A3667 is one of the best known examples of a galaxy cluster hosting a double radio relic and has been reported to contain a faint radio halo and bridge. The origins of radio halos, relics and bridges is still unclear, however galaxy cluster mergers seems to be an important factor. We clearly detect the North-West (NW) and South-East (SE) radio relics in A3667 and find an integrated flux density at 149 MHz of 28.1 +/- 1.7 and 2.4 +/- 0.1 Jy, respectively, with an average spectral index, between 120 and 1400 MHz, of -0.9 +/- 0.1 for both relics. We find evidence of a spatial variation in the spectral index across the NW relic steepening towards the centre of the cluster, which indicates an ageing electron population. These properties are consistent with higher frequency observations. We detect emission that could be associated with a radio halo and bridge. How- ever, due to the presence of poorly sampled large-scale Galactic emission and blended point sources we are unable to verify the exact nature of these features., 17 pages, 11 figures

We present a baseline sensitivity analysis of the Hydrogen Epoch of Reionization Array (HERA) and its build-out stages to one-point statistics (variance, skewness, and kurtosis) of redshifted 21 cm intensity fluctuation from the Epoch of Reionization (EoR) based on realistic mock observations. By developing a full-sky 21 cm lightcone model, taking into account the proper field of view and frequency bandwidth, utilising a realistic measurement scheme, and assuming perfect foreground removal, we show that HERA will be able to recover statistics of the sky model with high sensitivity by averaging over measurements from multiple fields. All build-out stages will be able to detect variance, while skewness and kurtosis should be detectable for HERA128 and larger. We identify sample variance as the limiting constraint of the variance measurement while skewness and kurtosis measurements will be primarily limited by thermal noise. The sensitivity can be improved by performing frequency binning and windowing. In addition, we find that strong sample variance fluctuation in the kurtosis measured from an individual field of observation indicates the present of outlying cold or hot regions in the underlying fluctuations, a feature that can potentially be used as an EoR bubble indicator. This manuscript is altered from the originally published paper to reflect corrections in the erratum., Comment: 14 pages, 10 figures, 2 tables, altered from the originally published paper to reflect corrections in the erratum

The External Calibrator for Hydrogen Observatories (ECHO) is an ongoing project to develop precision calibration sources for wide-field, low frequency, radio arrays such as the Murchison Widefield Array and the Hydrogen Epoch of Reionization Array which target 21cm radiation from primordial times. The accuracy of existing antenna beam models have been found to be a key factor limiting the realized sensitivity of the analysis pipelines. The precision required motivates us to go beyond applying a single model to each antenna in the array characterization and instead measure the response of the as-built system. We have developed a system for mapping the response pattern with a drone-mounted transmitter. Here we report a first series of tests, including our first test maps of known calibration standards and a second test mapping an MWA tile.

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., 13 pages, 9 figures. Replaced to match accepted ApJ version

The Epoch of Reionization (EoR) is an uncharted era in our Universe's history during which the birth of the first stars and galaxies led to the ionization of neutral hydrogen in the intergalactic medium. There are many experiments investigating the EoR by tracing the 21cm line of neutral hydrogen. Because this signal is very faint and difficult to isolate, it is crucial to develop analysis techniques that maximize sensitivity and suppress contaminants in data. It is also imperative to understand the trade-offs between different analysis methods and their effects on power spectrum estimates. Specifically, with a statistical power spectrum detection in HERA's foreseeable future, it has become increasingly important to understand how certain analysis choices can lead to the loss of the EoR signal. In this paper, we focus on signal loss associated with power spectrum estimation. We describe the origin of this loss using both toy models and data taken by the 64-element configuration of the Donald C. Backer Precision Array for Probing the Epoch of Reionization (PAPER). In particular, we highlight how detailed investigations of signal loss have led to a revised, higher 21cm power spectrum upper limit from PAPER-64. Additionally, we summarize errors associated with power spectrum error estimation that were previously unaccounted for. We focus on a subset of PAPER-64 data in this paper; revised power spectrum limits from the PAPER experiment are presented in a forthcoming paper by Kolopanis et al. (in prep.) and supersede results from previously published PAPER analyses., 25 pages, 18 figures, Accepted by ApJ

We are developing the Precision Array for Probing the Epoch of Reionization (PAPER) to detect 21cm emission from the early Universe, when the first stars and galaxies were forming. We describe the overall experiment strategy and architecture and summarize two PAPER deployments: a 4-antenna array in the low-RFI environment of Western Australia and an 8-antenna array at our prototyping site in Green Bank, WV. From these activities we report on system performance, including primary beam model verification, dependence of system gain on ambient temperature, measurements of receiver and overall system temperatures, and characterization of the RFI environment at each deployment site. We present an all-sky map synthesized between 139 MHz and 174 MHz using data from both arrays that reaches down to 80 mJy (4.9 K, for a beam size of 2.15e-5 steradians at 154 MHz), with a 10 mJy (620 mK) thermal noise level that indicates what would be achievable with better foreground subtraction. We calculate angular power spectra ($C_\ell$) in a cold patch and determine them to be dominated by point sources, but with contributions from galactic synchrotron emission at lower radio frequencies and angular wavemodes. Although the cosmic variance of foregrounds dominates errors in these power spectra, we measure a thermal noise level of 310 mK at $\ell=100$ for a 1.46-MHz band centered at 164.5 MHz. This sensitivity level is approximately three orders of magnitude in temperature above the level of the fluctuations in 21cm emission associated with reionization.

Experiments that pursue detection of signals from the Epoch of Reionization (EoR) are relying on spectral smoothness of source spectra at low frequencies. This article empirically explores the effect of foreground spectra on EoR experiments by measuring high-resolution full-polarization spectra for the 586 brightest unresolved sources in one of the Murchison Widefield Array (MWA) EoR fields using 45 h of observation. A novel peeling scheme is used to subtract 2500 sources from the visibilities with ionospheric and beam corrections, resulting in the deepest, confusion-limited MWA image so far. The resulting spectra are found to be affected by instrumental effects, which limit the constraints that can be set on source-intrinsic spectral structure. The sensitivity and power-spectrum of the spectra are analysed, and it is found that the spectra of residuals are dominated by point spread function sidelobes from nearby undeconvolved sources. We release a catalogue describing the spectral parameters for each measured source.

Detection of the cosmological neutral hydrogen signal from the Epoch of Reionization, and estimation of its basic physical parameters, is the principal scientific aim of many current low-frequency radio telescopes. Here we describe the Cosmological HI Power Spectrum Estimator (CHIPS), an algorithm developed and implemented with data from the Murchison Widefield Array (MWA), to compute the two-dimensional and spherically-averaged power spectrum of brightness temperature fluctuations. The principal motivations for CHIPS are the application of realistic instrumental and foreground models to form the optimal estimator, thereby maximising the likelihood of unbiased signal estimation, and allowing a full covariant understanding of the outputs. CHIPS employs an inverse-covariance weighting of the data through the maximum likelihood estimator, thereby allowing use of the full parameter space for signal estimation ("foreground suppression"). We describe the motivation for the algorithm, implementation, application to real and simulated data, and early outputs. Upon application to a set of 3 hours of data, we set a 2$\sigma$ upper limit on the EoR dimensionless power at $k=0.05$~h.Mpc$^{-1}$ of $\Delta_k^2, Comment: 34 pages, accepted for publication in ApJ

The Hydrogen Epoch of Reionization Array (HERA) is a staged experiment to measure 21 cm emission from the primordial intergalactic medium (IGM) throughout cosmic reionization ($z=6-12$), and to explore earlier epochs of our Cosmic Dawn ($z\sim30$). During these epochs, early stars and black holes heated and ionized the IGM, introducing fluctuations in 21 cm emission. HERA is designed to characterize the evolution of the 21 cm power spectrum to constrain the timing and morphology of reionization, the properties of the first galaxies, the evolution of large-scale structure, and the early sources of heating. The full HERA instrument will be a 350-element interferometer in South Africa consisting of 14-m parabolic dishes observing from 50 to 250 MHz. Currently, 19 dishes have been deployed on site and the next 18 are under construction. HERA has been designated as an SKA Precursor instrument. In this paper, we summarize HERA's scientific context and provide forecasts for its key science results. After reviewing the current state of the art in foreground mitigation, we use the delay-spectrum technique to motivate high-level performance requirements for the HERA instrument. Next, we present the HERA instrument design, along with the subsystem specifications that ensure that HERA meets its performance requirements. Finally, we summarize the schedule and status of the project. We conclude by suggesting that, given the realities of foreground contamination, current-generation 21 cm instruments are approaching their sensitivity limits. HERA is designed to bring both the sensitivity and the precision to deliver its primary science on the basis of proven foreground filtering techniques, while developing new subtraction techniques to unlock new capabilities. The result will be a major step toward realizing the widely recognized scientific potential of 21 cm cosmology., Comment: 26 pages, 24 figures, 2 tables

Current-generation low frequency interferometers constructed with the objective of detecting the high-redshift 21 cm background, aim to generate power spectra of the brightness-temperature contrast of neutral hydrogen in primordial intergalactic medium. Two-dimensional power spectra (power in Fourier modes parallel and perpendicular to the line of sight) formed from interferometric visibilities have been shown to delineate a boundary between spectrally-smooth foregrounds (known as the wedge) and spectrally-structured 21 cm background emission (the EoR-window). However, polarized foregrounds are known to possess spectral structure due to Faraday rotation, which can leak into the EoR window. In this work, we create and analyze 2D power spectra from the PAPER-32 imaging array in Stokes I, Q, U and V. These allow us to observe and diagnose systematic effects in our calibration at high signal-to-noise within the Fourier space most relevant to EoR experiments. We observe well-defined windows in the Stokes visibilities, with Stokes Q, U and V power spectra sharing a similar wedge shape to that seen in Stokes I. With modest polarization calibration, we see no evidence that polarization calibration errors move power outside the wedge in any Stokes visibility, to the noise levels attained. Deeper integrations will be required to confirm that this behavior persists to the depth required for EoR detection., 10 pages, 7 figures, accepted by ApJ

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., Comment: accepted to ApJ

The Murchison Widefield Array (MWA) is a Square Kilometre Array (SKA) Precursor. The telescope is located at the Murchison Radio--astronomy Observatory (MRO) in Western Australia (WA). The MWA consists of 4096 dipoles arranged into 128 dual polarisation aperture arrays forming a connected element interferometer that cross-correlates signals from all 256 inputs. A hybrid approach to the correlation task is employed, with some processing stages being performed by bespoke hardware, based on Field Programmable Gate Arrays (FPGAs), and others by Graphics Processing Units (GPUs) housed in general purpose rack mounted servers. The correlation capability required is approximately 8 TFLOPS (Tera FLoating point Operations Per Second). The MWA has commenced operations and the correlator is generating 8.3 TB/day of correlation products, that are subsequently transferred 700 km from the MRO to Perth (WA) in real-time for storage and offline processing. In this paper we outline the correlator design, signal path, and processing elements and present the data format for the internal and external interfaces., Comment: 17 pages, 9 figures. Accepted for publication in PASA. Some figures altered to meet astro-ph submission requirements

In this paper, we report new limits on 21cm emission from cosmic reionization based on a 135-day observing campaign with a 64-element deployment of the Donald C. Backer Precision Array for Probing the Epoch of Reionization (PAPER) in South Africa. This work extends the work presented in Parsons et al. (2014) with more collecting area, a longer observing period, improved redundancy-based calibration, optimal fringe-rate filtering, and improved power-spectral analysis using optimal quadratic estimators. The result is a new $2\sigma$ upper limit on $\Delta^{2}(k)$ of (22.4 mK)$^2$ in the range $0.15 < k < 0.5h\ {\rm Mpc}^{-1}$ at $z = 8.4$. This represents a three-fold improvement over the previous best upper limit. As we discuss in more depth in a forthcoming paper (Pober et al. 2015, in prep), this upper limit supports and extends previous evidence against extremely cold reionization scenarios. We conclude with a discussion of implications for future 21cm reionization experiments, including the newly funded Hydrogen Epoch of Reionization Array (HERA). $\textbf{The limits presented in this paper have been retracted: The erratum can be found in Appendix A.}$, Comment: Updated manuscript with an erratum in appendix A

We present constraints on both the kinetic temperature of the intergalactic medium (IGM) at z=8.4, and on models for heating the IGM at high-redshift with X-ray emission from the first collapsed objects. These constraints are derived using a semi-analytic method to explore the new measurements of the 21 cm power spectrum from the Donald C. Backer Precision Array for Probing the Epoch of Reionization (PAPER), which were presented in a companion paper, Ali et al. (2015). Twenty-one cm power spectra with amplitudes of hundreds of mK^2 can be generically produced if the kinetic temperature of the IGM is significantly below the temperature of the Cosmic Microwave Background (CMB); as such, the new results from PAPER place lower limits on the IGM temperature at z=8.4. Allowing for the unknown ionization state of the IGM, our measurements find the IGM temperature to be above ~5 K for neutral fractions between 10% and 85%, above ~7 K for neutral fractions between 15% and 80%, or above ~10 K for neutral fractions between 30% and 70%. We also calculate the heating of the IGM that would be provided by the observed high redshift galaxy population, and find that for most models, these galaxies are sufficient to bring the IGM temperature above our lower limits. However, there are significant ranges of parameter space that could produce a signal ruled out by the PAPER measurements; models with a steep drop-off in the star formation rate density at high redshifts or with relatively low values for the X-ray to star formation rate efficiency of high redshift galaxies are generally disfavored. The PAPER measurements are consistent with (but do not constrain) a hydrogen spin temperature above the CMB temperature, a situation which we find to be generally predicted if galaxies fainter than the current detection limits of optical/NIR surveys are included in calculations of X-ray heating., Comment: companion paper to Ali et al. (2015), ApJ 809, 61; matches version accepted to ApJ; 11 pages, 7 figures

The Murchison Widefield Array (MWA) is a new low-frequency interferometric radio telescope built in Western Australia at one of the locations of the future Square Kilometre Array (SKA). We describe the automated radio-frequency interference (RFI) detection strategy implemented for the MWA, which is based on the AOFlagger platform, and present 72-231-MHz RFI statistics from 10 observing nights. RFI detection removes 1.1% of the data. RFI from digital TV (DTV) is observed 3% of the time due to occasional ionospheric or atmospheric propagation. After RFI detection and excision, almost all data can be calibrated and imaged without further RFI mitigation efforts, including observations within the FM and DTV bands. The results are compared to a previously published Low-Frequency Array (LOFAR) RFI survey. The remote location of the MWA results in a substantially cleaner RFI environment compared to LOFAR's radio environment, but adequate detection of RFI is still required before data can be analysed. We include specific recommendations designed to make the SKA more robust to RFI, including: the availability of sufficient computing power for RFI detection; accounting for RFI in the receiver design; a smooth band-pass response; and the capability of RFI detection at high time and frequency resolution (second and kHz-scale respectively)., Comment: Accepted for publication in PASA

Detection of 21~cm emission of HI from the epoch of reionization, at redshifts z>6, is limited primarily by foreground emission. We investigate the signatures of wide-field measurements and an all-sky foreground model using the delay spectrum technique that maps the measurements to foreground object locations through signal delays between antenna pairs. We demonstrate interferometric measurements are inherently sensitive to all scales, including the largest angular scales, owing to the nature of wide-field measurements. These wide-field effects are generic to all observations but antenna shapes impact their amplitudes substantially. A dish-shaped antenna yields the most desirable features from a foreground contamination viewpoint, relative to a dipole or a phased array. Comparing data from recent Murchison Widefield Array observations, we demonstrate that the foreground signatures that have the largest impact on the HI signal arise from power received far away from the primary field of view. We identify diffuse emission near the horizon as a significant contributing factor, even on wide antenna spacings that usually represent structures on small scales. For signals entering through the primary field of view, compact emission dominates the foreground contamination. These two mechanisms imprint a characteristic "pitchfork" signature on the "foreground wedge" in Fourier delay space. Based on these results, we propose that selective down-weighting of data based on antenna spacing and time can mitigate foreground contamination substantially by a factor ~100 with negligible loss of sensitivity., Published in ApJ

Polarized foreground emission is a potential contaminant of attempts to measure the fluctuation power spectrum of highly redshifted 21 cm H i emission from the epoch of reionization. Using the Donald C. Backer Precision Array for Probing the Epoch of Reionization, we present limits on the observed power spectra of all four Stokes parameters in two frequency bands, centered at 126 MHz (z = 10.3) and 164 MHz (z = 7.66), for a three-month observing campaign of a deployment involving 32 antennas, for which results on the unpolarized power spectrum have been reported at z = 7.7 (by Parsons et al.) and at (by Jacobs et al.). The power spectra in this paper are processed in the same way as by those authors, and show no definitive detection of polarized power. This nondetection is consistent with what is known about polarized sources, combined with the suppression of polarized power by fluctuations in the ionospheric rotation measure, which can strongly affect Stokes Q and U. We are able to show that the net effect of polarized leakage is a negligible contribution at the levels of the limits reported by Parsons et al. and Jacobs et al.