M. E. Bell

wsclean: an implementation of a fast, generic wide-field imager for radio astronomy

Astronomical widefield 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 widefield interferometric imager that uses the w-stacking algorithm and can make use of the w-snapshot algorithm. The performance dependencies of CASAs 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 polarisation 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.

GaLactic and Extragalactic All-sky Murchison Widefield Array (GLEAM) survey – I. A low-frequency extragalactic catalogue

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. 1 This is the first in a series of publications describing the GLEAM survey results.

GLEAM: The GaLactic and Extragalactic All-Sky MWA Survey

© Astronomical Society of Australia 2015; published by Cambridge University Press. This is an Open Access article distributed in accordance with the terms of the Creative Commons Attribution (CC BY 4.0) license, which permits others to distribute, remix, adapt and build upon this work, for commercial use, provided the original work is properly cited. See: http://creativecommons.org/licenses/by/4.0/

The Murchison Widefield Array Commissioning Survey: A Low-Frequency Catalogue of 14 110 Compact Radio Sources over 6 100 Square Degrees

We present the results of an approximately 6,100 square degree 104--196MHz radio sky survey performed with the Murchison Widefield Array during instrument commissioning between 2012 September and 2012 December: the Murchison Widefield Array Commissioning Survey (MWACS). The data were taken as meridian drift scans with two different 32-antenna sub-arrays that were available during the commissioning period. The survey covers approximately 20.5 h < Right Ascension (RA) < 8.5 h, -58 deg < Declination (Dec) < -14 deg over three frequency bands centred on 119, 150 and 180 MHz, with image resolutions of 6--3 arcmin. The catalogue has 3-arcmin angular resolution and a typical noise level of 40 mJy/beam, with reduced sensitivity near the field boundaries and bright sources. We describe the data reduction strategy, based upon mosaiced snapshots, flux density calibration and source-finding method. We present a catalogue of flux density and spectral index measurements for 14,110 sources, extracted from the mosaic, 1,247 of which are sub-components of complexes of sources.

Calibrating high-precision Faraday rotation measurements for LOFAR and the next generation of low-frequency radio telescopes

Faraday rotation measurements using the current and next generation of low-frequency radio telescopes will provide a powerful probe of astronomical magnetic fields. However, achieving the full potential of these measurements requires accurate removal of the time-variable ionospheric Faraday rotation contribution. We present ionFR, a code that calculates the amount of ionospheric Faraday rotation for a specific epoch, geographic location, and line-of-sight. ionFR uses a number of publicly available, GPS-derived total electron content maps and the most recent release of the International Geomagnetic Reference Field. We describe applications of this code for the calibration of radio polarimetric observations, and demonstrate the high accuracy of its modeled ionospheric Faraday rotations using LOFAR pulsar observations. These show that we can accurately determine some of the highest-precision pulsar rotation measures ever achieved. Precision rotation measures can be used to monitor rotation measure variations - either intrinsic or due to the changing line-of-sight through the interstellar medium. This calibration is particularly important for nearby sources, where the ionosphere can contribute a significant fraction of the observed rotation measure. We also discuss planned improvements to ionFR, as well as the importance of ionospheric Faraday rotation calibration for the emerging generation of low-frequency radio telescopes, such as the SKA and its pathfinders.

Radio Continuum Surveys with Square Kilometre Array Pathfinders

In the lead-up to the Square Kilometre Array (SKA) project, several next-generation radio telescopes and upgrades are already being built around the world. These include APERTIF (The Netherlands), ASKAP (Australia), e-MERLIN (UK), VLA (USA), e-EVN (based in Europe), LOFAR (The Netherlands), MeerKAT (South Africa), and the Murchison Widefield Array. Each of these new instruments has different strengths, and coordination of surveys between them can help maximise the science from each of them. A radio continuum survey is being planned on each of them with the primary science objective of understanding the formation and evolution of galaxies over cosmic time, and the cosmological parameters and large-scale structures which drive it. In pursuit of this objective, the different teams are developing a variety of new techniques, and refining existing ones. To achieve these exciting scientific goals, many technical challenges must be addressed by the survey instruments. Given the limited resources of the global radio-astronomical community, it is essential that we pool our skills and knowledge. We do not have sufficient resources to enjoy the luxury of re-inventing wheels. We face significant challenges in calibration, imaging, source extraction and measurement, classification and cross-identification, redshift determination, stacking, and data-intensive research. As these instruments extend the observational parameters, we will face further unexpected challenges in calibration, imaging, and interpretation. If we are to realise the full scientific potential of these expensive instruments, it is essential that we devote enough resources and careful study to understanding the instrumental effects and how they will affect the data. We have established an SKA Radio Continuum Survey working group, whose prime role is to maximise science from these instruments by ensuring we share resources and expertise across the projects. Here we describe these projects, their science goals, and the technical challenges which are being addressed to maximise the science return.

M 87 at metre wavelengths: the LOFAR picture

Context. M87 is a giant elliptical galaxy located in the centre of the Virgo cluster, which harbours a supermassive black hole of mass 6.4×109 M, whose activity is responsible for the extended (80 kpc) radio lobes that surround the galaxy. The energy generated by matter falling onto the central black hole is ejected and transferred to the intra-cluster medium via a relativistic jet and morphologically complex systems of buoyant bubbles, which rise towards the edges of the extended halo. Aims. To place constraints on past activity cycles of the active nucleus, images of M 87 were produced at low radio frequencies never explored before at these high spatial resolution and dynamic range. To disentangle different synchrotron models and place constraints on source magnetic field, age and energetics, we also performed a detailed spectral analysis of M 87 extended radio-halo. Methods. We present the first observations made with the new Low-Frequency Array (LOFAR) of M 87 at frequencies down to 20 MHz. Three observations were conducted, at 15−30 MHz, 30−77 MHz and 116−162 MHz. We used these observations together with archival data to produce a low-frequency spectral index map and to perform a spectral analysis in the wide frequency range 30 MHz–10 GHz. Results. We do not find any sign of new extended emissions; on the contrary the source appears well confined by the high pressure of the intracluster medium. A continuous injection of relativistic electrons is the model that best fits our data, and provides a scenario in which the lobes are still supplied by fresh relativistic particles from the active galactic nuclei. We suggest that the discrepancy between the low-frequency radiospectral slope in the core and in the halo implies a strong adiabatic expansion of the plasma as soon as it leaves the core area. The extended halo has an equipartition magnetic field strength of 10 μG, which increases to 13 μG in the zones where the particle flows are more active. The continuous injection model for synchrotron ageing provides an age for the halo of 40 Myr, which in turn provides a jet kinetic power of 6−10 × 1044 erg s−1.

First LOFAR observations at very low frequencies of cluster-scale non-thermal emission: the case of Abell 2256

Abell 2256 is one of the best known examples of a galaxy cluster hosting large-scale diffuse radio emission that is unrelated to individual galaxies. It contains both a giant radio halo and a relic, as well as a number of head-tail sources and smaller diffuse steep-spectrum radio sources. The origin of radio halos and relics is still being debated, but over the last years it has become clear that the presence of these radio sources is closely related to galaxy cluster merger events. Here we present the results from the first LOFAR low band antenna (LBA) observations of Abell 2256 between 18 and 67 MHz. To our knowledge, the image presented in this paper at 63 MHz is the deepest ever obtained at frequencies below 100 MHz in general. Both the radio halo and the giant relic are detected in the image at 63 MHz, and the diffuse radio emission remains visible at frequencies as low as 20 MHz. The observations confirm the presence of a previously claimed ultra-steep spectrum source to the west of the cluster center with a spectral index of -2.3 +/- 0.4 between 63 and 153 MHz. The steep spectrum suggests that this source is an old part of a head-tail radio source in the cluster. For the radio relic we find an integrated spectral index of -0.81 +/- 0.03, after removing the flux contribution from the other sources. This is relatively flat which could indicate that the efficiency of particle acceleration at the shock substantially changed in the last similar to 0.1 Gyr due to an increase of the shock Mach number. In an alternative scenario, particles are re-accelerated by some mechanism in the downstream region of the shock, resulting in the relatively flat integrated radio spectrum. In the radio halo region we find indications of low-frequency spectral steepening which may suggest that relativistic particles are accelerated in a rather inhomogeneous turbulent region.

Investigating the disc-jet coupling in accreting compact objects using the black hole candidate Swift J1753.5-0127

In studies of accreting black holes in binary systems, empirical relations have been proposed to quantify the coupling between accretion processes and ejection mechanisms. These processes are probed, respectively, by means of X-ray and radio/optical-infrared observations. The relations predict, given certain accretion conditions, the expected energy output in the form of a jet. We investigated this coupling by studying the black hole candidate Swift J1753.5−0127, via multiwavelength-coordinated observations over a period of ∼4 yr. We present the results of our campaign showing that, all along the outburst, the source features a jet that is fainter than expected from the empirical correlation between the radio and the X-ray luminosities in a hard spectral state. Because the jet is so weak in this system the near-infrared emission is, unusually for this state and luminosity, dominated by thermal emission from the accretion disc. We briefly discuss the importance and the implications of a precise determination of both the slope and the normalization of the correlations, listing some possible parameters that broad-band jet models should take into account to explain the population of sources characterized by a dim jet. We also investigate whether our data can give any hint on the nature of the compact object in the system, since its mass has not been dynamically measured.

The Low-Frequency Environment of the Murchison Widefield Array: Radio-Frequency Interference Analysis and Mitigation

This is the Accepted Manuscript version of the following article: A. R. Offringa, et al., “The low-frequency environment of the Murchison Widefield Array: radio-frequency interference analysis and mitigation”, Publications of the Astronomical Society of Australia, Vol. 32, March 2015. The final published version is available at: https://doi.org/10.1017/pasa.2015.7