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.
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.
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.
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.)
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.
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.
ABSTRACT 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.
ABSTRACT
CXO J133815.6+043255 is an ultraluminous X-ray source (ULX) with ultraviolet, optical, and radio counterparts located 10 kpc away from the nucleus of the galaxy NGC 5252. Optical ...spectroscopic studies indicate that the ULX is kinematically associated with NGC 5252; yet, the compactness of its radio emission could not rule out the possibility that the ULX is a background blazar. We present follow-up VLBA radio observations that are able to resolve the compact radio emission of the ULX into two components, making the blazar scenario very unlikely. The east component is extended at 4.4 GHz and its detection also at 7.6 GHz reveals a steep spectral index. The west component is only detected at 4.4 GHz, is not firmly resolved, and has a flatter spectral index. Considering that the west component hosts the radio core, we constrain the black hole mass of the ULX to 103.5 < MBH ≲2 × 106 M⊙ and its Eddington ratio to ∼10−3. The ULX is thus most likely powered by an intermediate-mass black hole or low-mass AGN. Our results constitute the first discovery of a multicomponent radio jet in a ULX and possible intermediate-mass black hole.
ABSTRACT We present deep polarimetric observations at 154 MHz with the Murchison Widefield Array (MWA), covering 625 deg2 centered on = 0hand δ = −27°. The sensitivity available in our deep ...observations allows an in-band, frequency-dependent analysis of polarized structure for the first time at long wavelengths. Our analysis suggests that the polarized structures are dominated by intrinsic emission but may also have a foreground Faraday screen component. At these wavelengths, the compactness of the MWA baseline distribution provides excellent snapshot sensitivity to large-scale structure. The observations are sensitive to diffuse polarized emission at ∼54′ resolution with a sensitivity of 5.9 mJy beam−1 and compact polarized sources at ∼2 4 resolution with a sensitivity of 2.3 mJy beam−1 for a subset (400 deg2) of this field. The sensitivity allows the effect of ionospheric Faraday rotation to be spatially and temporally measured directly from the diffuse polarized background. Our observations reveal large-scale structures (∼1°-8° in extent) in linear polarization clearly detectable in ∼2 minute snapshots, which would remain undetectable by interferometers with minimum baseline lengths of >110 m at 154 MHz. The brightness temperature of these structures is on average 4 K in polarized intensity, peaking at 11 K. Rotation measure synthesis reveals that the structures have Faraday depths ranging from −2 to 10 rad m−2 with a large fraction peaking at approximately +1 rad m−2. We estimate a distance of 51 20 pc to the polarized emission based on measurements of the in-field pulsar J2330-2005. We detect four extragalactic linearly polarized point sources within the field in our compact source survey. Based on the known polarized source population at 1.4 GHz and non-detections at 154 MHz, we estimate an upper limit on the depolarization ratio of 0.08 from 1.4 GHz to 154 MHz.