Astronomical wide-field imaging of interferometric radio data is computationally expensive, especially for the large data volumes created by modern non-coplanar many-element arrays. We present a new ...wide-field interferometric imager that uses the w-stacking algorithm and can make use of the w-snapshot algorithm. The performance dependences of casa's w-projection and our new imager are analysed and analytical functions are derived that describe the required computing cost for both imagers. On data from the Murchison Widefield Array, we find our new method to be an order of magnitude faster than w-projection, as well as being capable of full-sky imaging at full resolution and with correct polarization correction. We predict the computing costs for several other arrays and estimate that our imager is a factor of 2–12 faster, depending on the array configuration. We estimate the computing cost for imaging the low-frequency Square Kilometre Array observations to be 60 PetaFLOPS with current techniques. We find that combining w-stacking with the w-snapshot algorithm does not significantly improve computing requirements over pure w-stacking. The source code of our new imager is publicly released.
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 × 103 mK2 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.
Abstract
We compute the spherically averaged power spectrum from four seasons of data obtained for the Epoch of Reionization (EoR) project observed with the Murchison Widefield Array (MWA). We ...measure the EoR power spectrum over k = 0.07–3.0 h Mpc−1 at redshifts $z$ = 6.5–8.7. The largest aggregation of 110 h on EoR0 high band (3340 observations), yields a lowest measurement of (43 mK)2 = 1.8 × 103 mK2 at k = 0.14 h Mpc−1 and $z$ = 6.5 (2σ thermal noise plus sample variance). Using the Real-Time System to calibrate and the CHIPS pipeline to estimate power spectra, we select the best observations from the central five pointings within the 2013–2016 observing seasons, observing three independent fields and in two frequency bands. This yields 13 591 2-min snapshots (453 h), based on a quality assurance metric that measures ionospheric activity. We perform another cut to remove poorly calibrated data, based on power in the foreground-dominated and EoR-dominated regions of the two-dimensional power spectrum, reducing the set to 12 569 observations (419 h). These data are processed in groups of 20 observations, to retain the capacity to identify poor data, and used to analyse the evolution and structure of the data over field, frequency, and data quality. We subsequently choose the cleanest 8935 observations (298 h of data) to form integrated power spectra over the different fields, pointings, and redshift ranges.
21 cm epoch of reionization (EoR) observations promise to transform our understanding of galaxy formation, but these observations are impossible without unprecedented levels of instrument ...calibration. We present end-to-end simulations of a full EoR power spectrum (PS) analysis including all of the major components of a real data processing pipeline: models of astrophysical foregrounds and EoR signal, frequency-dependent instrument effects, sky-based antenna calibration, and the full PS analysis. This study reveals that traditional sky-based per-frequency antenna calibration can only be implemented in EoR measurement analyses if the calibration model is unrealistically accurate. For reasonable levels of catalogue completeness, the calibration introduces contamination in otherwise foreground-free PS modes, precluding a PS measurement. We explore the origin of this contamination and potential mitigation techniques. We show that there is a strong joint constraint on the precision of the calibration catalogue and the inherent spectral smoothness of antennas, and that this has significant implications for the instrumental design of the SKA (Square Kilometre Array) and other future EoR observatories.
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 Δ2 ≤ 2.39 × 103 mK2 at k = 0.59 h Mpc−1 and z = 6.5.
ABSTRACT
Observations in the lowest Murchison Widefield Array (MWA) band between 75 and 100 MHz have the potential to constrain the distribution of neutral hydrogen in the intergalactic medium at ...redshift ∼13–17. Using 15 h 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 h of the cleanest data. Using these data, we obtain 2σ upper limits on the 21 cm power spectrum in the range of $0.1~ h~{\mathrm{ Mpc}}^{-1}\lessapprox k \lessapprox 1 ~ ~h~{\mathrm{ 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 $k=0.14 ~h~{\mathrm{ Mpc}}^{-1}$ and at z = 15.2 with a possibility of a few per cent of signal loss due to direction-independent calibration.
Detection of the epoch of reionization H i signal requires a precise understanding of the intervening galaxies and AGN, both for instrumental calibration and foreground removal. We present a ...catalogue of 7394 extragalactic sources at 182 MHz detected in the RA = 0 field of the Murchison Widefield Array Epoch of Reionization observation programme. Motivated by unprecedented requirements for precision and reliability we develop new methods for source finding and selection. We apply machine learning methods to self-consistently classify the relative reliability of 9490 source candidates. A subset of 7466 are selected based on reliability class and signal-to-noise ratio criteria. These are statistically cross-matched to four other radio surveys using both position and flux density information. We find 7369 sources to have confident matches, including 90 partially resolved sources that split into a total of 192 sub-components. An additional 25 unmatched sources are included as new radio detections. The catalogue sources have a median spectral index of −0.85. Spectral flattening is seen towards lower frequencies with a median of −0.71 predicted at 182 MHz. The astrometric error is 7 arcsec compared to a 2.3 arcmin beam FWHM. The resulting catalogue covers ∼1400 deg2 and is complete to approximately 80 mJy within half beam power. This provides the most reliable discrete source sky model available to date in the MWA EoR0 field for precision foreground subtraction.
PDF
