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.
We present a survey for transient and variable sources, on time-scales from 28 s to ∼1 yr, using the Murchison Widefield Array (MWA) at 182 MHz. Down to a detection threshold of 0.285 Jy, no ...transient candidates were identified, making this the most constraining low-frequency survey to date and placing a limit on the surface density of transients of <4.1 × 10−7 deg−2 for the shortest time-scale considered. At these frequencies, emission from Fast Radio Bursts (FRBs) is expected to be detectable in the shortest time-scale images without any corrections for interstellar or intergalactic dispersion. At an FRB limiting flux density of 7980 Jy, we find a rate of <82 FRBs per sky per day for dispersion measures <700 pc cm−3. Assuming a cosmological population of standard candles, our rate limits are consistent with the FRB rates obtained by Thornton et al. if they have a flat spectral slope. Finally, we conduct an initial variability survey of sources in the field with flux densities ≳0.5 Jy and identify no sources with significant variability in their light curves. However, we note that substantial further work is required to fully characterize both the short-term and low-level variability within this field.
ABSTRACT 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 mK2 at k = 0.27 h Mpc−1 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.
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.
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 origin of radio haloes, relics and bridges is still unclear, however galaxy cluster merger seems to be an important factor. We clearly detect the north-west (NW) and south-east 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. However, 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.
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−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.
Using the Murchison Widefield Array (MWA), the low-frequency Square Kilometre Array precursor located in Western Australia, we have completed the GaLactic and Extragalactic All-sky MWA (GLEAM) ...survey, and present the resulting extragalactic catalogue, utilizing the first year of observations. The catalogue covers 24 831 square degrees, over declinations south of +30... and Galactic latitudes outside 10... of the Galactic plane, excluding some areas such as the Magellanic Clouds. It contains 307 455 radio sources with 20 separate flux density measurements across 72-231 MHz, selected from a time- and frequency-integrated image centred at 200 MHz, with a resolution of ...2 arcmin. Over the catalogued region, we estimate that the catalogue is 90 per cent complete at 170 mJy, and 50 per cent complete at 55 mJy, and large areas are complete at even lower flux density levels. Its reliability is 99.97 per cent above the detection threshold of 5..., which itself is typically 50 mJy. These observations constitute the widest fractional bandwidth and largest sky area survey at radio frequencies to date, and calibrate the low-frequency flux density scale of the southern sky to better than 10 per cent. This paper presents details of the flagging, imaging, mosaicking and source extraction/characterization, as well as estimates of the completeness and reliability. All source measurements and images are available online. This is the first in a series of publications describing the GLEAM survey results. (ProQuest: ... denotes formulae/symbols omitted.)
ABSTRACT Detection of the cosmological neutral hydrogen signal from the Epoch of Reionization (EoR) and estimation of its basic physical parameters are principal scientific aims of many current ...low-frequency radio telescopes. Here we describe the Cosmological H i Power Spectrum Estimator (CHIPS), an algorithm developed and implemented with data from the Murchison Widefield Array, 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 maximizing 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 hr of data, we set a 2 upper limit on the EoR dimensionless power at Mpc−1 of mK2 in the redshift range z = 6.2-6.6, consistent with previous estimates.
We present a search for transient and variable radio sources at 154 MHz with the Murchison Widefield Array 32-tile prototype. 51 images were obtained that cover a field of view of 1430 deg2 centred ...on Hydra A. The observations were obtained over three days in 2010 March and three days in 2011 April and May. The mean cadence of the observations was 26 min and there was additional temporal information on day and year time-scales. We explore the variability of a sample of 105 low-frequency radio sources within the field. Four bright (S > 6 Jy) candidate variable radio sources were identified that displayed low levels of short time-scale variability (26 min). We conclude that this variability is likely caused by simplifications in the calibration strategy or ionospheric effects. On the time-scale of 1 yr we find two sources that show significant variability. We attribute this variability to either refractive scintillation or intrinsic variability. No radio transients were identified and we place an upper limit on the surface density of sources ρ < 7.5 × 10−5 deg−2 with flux densities >5.5 Jy, and characteristic time-scales of both 26 min and 1 yr.
We analyse a 154 MHz image made from a 12 h observation with the Murchison Widefield Array (MWA) to determine the noise contribution and behaviour of the source counts down to 30 mJy. The MWA image ...has a bandwidth of 30.72 MHz, a field-of-view within the half-power contour of the primary beam of 570 deg2, a resolution of 2.3 arcmin and contains 13 458 sources above 5σ. The rms noise in the centre of the image is 4–5 mJy beam−1. The MWA counts are in excellent agreement with counts from other instruments and are the most precise ever derived in the flux density range 30–200 mJy due to the sky area covered. Using the deepest available source count data, we find that the MWA image is affected by sidelobe confusion noise at the ≈3.5 mJy beam−1 level, due to incompletely peeled and out-of-image sources, and classical confusion becomes apparent at ≈1.7 mJy beam−1. This work highlights that (i) further improvements in ionospheric calibration and deconvolution imaging techniques would be required to probe to the classical confusion limit and (ii) the shape of low-frequency source counts, including any flattening towards lower flux densities, must be determined from deeper ≈150 MHz surveys as it cannot be directly inferred from higher frequency data.