Pulsars have been studied extensively over the last few decades and have proven instrumental in exploring a wide variety of physics. Discovering more pulsars emitting at low radio frequencies is crucial to further our understanding of spectral properties and emission mechanisms. The Murchison Widefield Array Voltage Capture System (MWA VCS) has been routinely used to study pulsars at low frequencies and discover new pulsars. The MWA VCS offers the unique opportunity of recording complex voltages from all individual antennas (tiles), which can be off-line beamformed or correlated/imaged at millisecond time resolution. Devising imaged-based methods for finding pulsar candidates, which can be verified in beamformed data, can accelerate the complete process and lead to more pulsar detections. Image-based searches for pulsar candidates can reduce the number of tied-array beams required, increasing compute resource efficiency. Despite a factor of $\sim$ 4 loss in sensitivity, searching for pulsar candidates in images from the MWA VCS, we can explore a larger parameter space, potentially leading to discoveries of pulsars missed by high-frequency surveys such as steep spectrum pulsars, exotic binary systems, or pulsars obscured in high-time resolution time series data by propagation effects. Image-based searches are also essential to probing parts of parameter space inaccessible to traditional beamformed searches with the MWA (e.g. at high dispersion measures). In this paper we describe the innovative approach and capability of dual-processing MWA VCS data, that is forming 1-s visibilities and sky images, finding pulsar candidates in these images, and verifying by forming tied-array beam. We developed and tested image-based methods of finding pulsar candidates, which are based on pulsar properties such as steep spectral index, polarisation and variability. The efficiency of these methodologies has been verified on known pulsars, and the main limitations explained in terms of sensitivity and low-frequency spectral turnover of some pulsars. No candidates were confirmed to be a new pulsar, but this new capability will now be applied to a larger subset of observations to accelerate pulsar discoveries with the MWA and potentially speed up future searches with the SKA-Low.

The Rapid ASKAP Continuum Survey (RACS) is the first large sky survey using the Australian Square Kilometre Array Pathfinder (ASKAP), covering the sky south of $+41^\circ$ declination. With ASKAP’s large, instantaneous field of view, ${\sim}31\,\mathrm{deg}^2$ , RACS observed the entire sky at a central frequency of 887.5 MHz using 903 individual pointings with 15 minute observations. This has resulted in the deepest radio survey of the full Southern sky to date at these frequencies. In this paper, we present the first Stokes I catalogue derived from the RACS survey. This catalogue was assembled from 799 tiles that could be convolved to a common resolution of $25^{\prime\prime}$ , covering a large contiguous region in the declination range $\delta=-80^{\circ}$ to $+30^\circ$ . The catalogue provides an important tool for both the preparation of future ASKAP surveys and for scientific research. It consists of $\sim$ 2.1 million sources and excludes the $|b|<5^{\circ}$ region around the Galactic plane. This provides a first extragalactic catalogue with ASKAP covering the majority of the sky ( $\delta<+30^{\circ}$ ). We describe the methods to obtain this catalogue from the initial RACS observations and discuss the verification of the data, to highlight its quality. Using simulations, we find this catalogue detects 95% of point sources at an integrated flux density of $\sim$ 5 mJy. Assuming a typical sky source distribution model, this suggests an overall 95% point source completeness at an integrated flux density $\sim$ 3 mJy. The catalogue will be available through the CSIRO ASKAP Science Data Archive (CASDA).

The Variables and Slow Transients Survey (VAST) on the Australian Square Kilometre Array Pathfinder (ASKAP) is designed to detect highly variable and transient radio sources on timescales from 5 s to $\sim\!5$ yr. In this paper, we present the survey description, observation strategy and initial results from the VAST Phase I Pilot Survey. This pilot survey consists of $\sim\!162$ h of observations conducted at a central frequency of 888 MHz between 2019 August and 2020 August, with a typical rms sensitivity of $0.24\ \mathrm{mJy\ beam}^{-1}$ and angular resolution of $12-20$ arcseconds. There are 113 fields, each of which was observed for 12 min integration time, with between 5 and 13 repeats, with cadences between 1 day and 8 months. The total area of the pilot survey footprint is 5 131 square degrees, covering six distinct regions of the sky. An initial search of two of these regions, totalling 1 646 square degrees, revealed 28 highly variable and/or transient sources. Seven of these are known pulsars, including the millisecond pulsar J2039–5617. Another seven are stars, four of which have no previously reported radio detection (SCR J0533–4257, LEHPM 2-783, UCAC3 89–412162 and 2MASS J22414436–6119311). Of the remaining 14 sources, two are active galactic nuclei, six are associated with galaxies and the other six have no multi-wavelength counterparts and are yet to be identified.

The GaLactic and Extragalactic All-sky Murchison Widefield Array (GLEAM) is a radio continuum survey at 76–227 MHz of the entire southern sky (Declination $<\!{+}30^{\circ}$ ) with an angular resolution of ${\approx}2$ arcmin. In this paper, we combine GLEAM data with optical spectroscopy from the 6dF Galaxy Survey to construct a sample of 1 590 local (median $z \approx 0.064$ ) radio sources with $S_{200\,\mathrm{MHz}} > 55$ mJy across an area of ${\approx}16\,700\,\mathrm{deg}^{2}$ . From the optical spectra, we identify the dominant physical process responsible for the radio emission from each galaxy: 73% are fuelled by an active galactic nucleus (AGN) and 27% by star formation. We present the local radio luminosity function for AGN and star-forming (SF) galaxies at 200 MHz and characterise the typical radio spectra of these two populations between 76 MHz and ${\sim}1$ GHz. For the AGN, the median spectral index between 200 MHz and ${\sim}1$ GHz, $\alpha_{\mathrm{high}}$ , is $-0.600 \pm 0.010$ (where $S \propto \nu^{\alpha}$ ) and the median spectral index within the GLEAM band, $\alpha_{\mathrm{low}}$ , is $-0.704 \pm 0.011$ . For the SF galaxies, the median value of $\alpha_{\mathrm{high}}$ is $-0.650 \pm 0.010$ and the median value of $\alpha_{\mathrm{low}}$ is $-0.596 \pm 0.015$ . Among the AGN population, flat-spectrum sources are more common at lower radio luminosity, suggesting the existence of a significant population of weak radio AGN that remain core-dominated even at low frequencies. However, around 4% of local radio AGN have ultra-steep radio spectra at low frequencies ( $\alpha_{\mathrm{low}} < -1.2$ ). These ultra-steep-spectrum sources span a wide range in radio luminosity, and further work is needed to clarify their nature.

We present the first Faraday rotation measure (RM) grid study of an individual low-mass cluster—the Fornax cluster—which is presently undergoing a series of mergers. Exploiting commissioning data for the POlarisation Sky Survey of the Universe’s Magnetism (POSSUM) covering a ${\sim}34$ square degree sky area using the Australian Square Kilometre Array Pathfinder (ASKAP), we achieve an RM grid density of ${\sim}25$ RMs per square degree from a 280-MHz band centred at 887 MHz, which is similar to expectations for forthcoming GHz-frequency ${\sim}3\pi$-steradian sky surveys. These data allow us to probe the extended magnetoionic structure of the cluster and its surroundings in unprecedented detail. We find that the scatter in the Faraday RM of confirmed background sources is increased by $16.8\pm2.4$ rad m−2 within 1$^\circ$ (360 kpc) projected distance to the cluster centre, which is 2–4 times larger than the spatial extent of the presently detectable X-ray-emitting intracluster medium (ICM). The mass of the Faraday-active plasma is larger than that of the X-ray-emitting ICM and exists in a density regime that broadly matches expectations for moderately dense components of the Warm-Hot Intergalactic Medium. We argue that forthcoming RM grids from both targeted and survey observations may be a singular probe of cosmic plasma in this regime. The morphology of the global Faraday depth enhancement is not uniform and isotropic but rather exhibits the classic morphology of an astrophysical bow shock on the southwest side of the main Fornax cluster, and an extended, swept-back wake on the northeastern side. Our favoured explanation for these phenomena is an ongoing merger between the main cluster and a subcluster to the southwest. The shock’s Mach angle and stand-off distance lead to a self-consistent transonic merger speed with Mach 1.06. The region hosting the Faraday depth enhancement also appears to show a decrement in both total and polarised radio emission compared to the broader field. We evaluate cosmic variance and free-free absorption by a pervasive cold dense gas surrounding NGC 1399 as possible causes but find both explanations unsatisfactory, warranting further observations. Generally, our study illustrates the scientific returns that can be expected from all-sky grids of discrete sources generated by forthcoming all-sky radio surveys.

In this paper, we describe the system design and capabilities of the Australian Square Kilometre Array Pathfinder (ASKAP) radio telescope at the conclusion of its construction project and commencement of science operations. ASKAP is one of the first radio telescopes to deploy phased array feed (PAF) technology on a large scale, giving it an instantaneous field of view that covers $31\,\textrm{deg}^{2}$ at $800\,\textrm{MHz}$. As a two-dimensional array of 36$\times$12 m antennas, with baselines ranging from 22 m to 6 km, ASKAP also has excellent snapshot imaging capability and 10 arcsec resolution. This, combined with 288 MHz of instantaneous bandwidth and a unique third axis of rotation on each antenna, gives ASKAP the capability to create high dynamic range images of large sky areas very quickly. It is an excellent telescope for surveys between 700 and $1800\,\textrm{MHz}$ and is expected to facilitate great advances in our understanding of galaxy formation, cosmology, and radio transients while opening new parameter space for discovery of the unknown.

We have found a class of circular radio objects in the Evolutionary Map of the Universe Pilot Survey, using the Australian Square Kilometre Array Pathfinder telescope. The objects appear in radio images as circular edge-brightened discs, about one arcmin diameter, that are unlike other objects previously reported in the literature. We explore several possible mechanisms that might cause these objects, but none seems to be a compelling explanation.

The Rapid ASKAP Continuum Survey (RACS) is the first large-area survey to be conducted with the full 36-antenna Australian Square Kilometre Array Pathfinder (ASKAP) telescope. RACS will provide a shallow model of the ASKAP sky that will aid the calibration of future deep ASKAP surveys. RACS will cover the whole sky visible from the ASKAP site in Western Australia and will cover the full ASKAP band of 700–1800 MHz. The RACS images are generally deeper than the existing NRAO VLA Sky Survey and Sydney University Molonglo Sky Survey radio surveys and have better spatial resolution. All RACS survey products will be public, including radio images (with $\sim$ 15 arcsec resolution) and catalogues of about three million source components with spectral index and polarisation information. In this paper, we present a description of the RACS survey and the first data release of 903 images covering the sky south of declination $+41^\circ$ made over a 288-MHz band centred at 887.5 MHz.

The low-frequency linearly polarised radio source population is largely unexplored. However, a renaissance in low-frequency polarimetry has been enabled by pathfinder and precursor instruments for the Square Kilometre Array. In this second paper from the POlarised GaLactic and Extragalactic All-Sky MWA Survey-the POlarised GLEAM Survey, or POGS-we present the results from our all-sky MWA Phase I Faraday Rotation Measure survey. Our survey covers nearly the entire Southern sky in the Declination range $-82^\circ$ to $+30^\circ$ at a resolution between around three and seven arcminutes (depending on Declination) using data in the frequency range 169−231 MHz. We have performed two targeted searches: the first covering 25 489 square degrees of sky, searching for extragalactic polarised sources; the second covering the entire sky South of Declination $+30^\circ$, searching for known pulsars. We detect a total of 517 sources with 200 MHz linearly polarised flux densities between 9.9 mJy and 1.7 Jy, of which 33 are known radio pulsars. All sources in our catalogues have Faraday rotation measures in the range $-328.07$ to $+279.62$ rad m−2. The Faraday rotation measures are broadly consistent with results from higher-frequency surveys, but with typically more than an order of magnitude improvement in the precision, highlighting the power of low-frequency polarisation surveys to accurately study Galactic and extragalactic magnetic fields. We discuss the properties of our extragalactic and known-pulsar source population, how the sky distribution relates to Galactic features, and identify a handful of new pulsar candidates among our nominally extragalactic source population.

The low-frequency polarisation properties of radio sources are poorly studied, particularly in statistical samples. However, the new generation of low-frequency telescopes, such as the Murchison Widefield Array (the precursor for the low-frequency component of the Square Kilometre Array) offers an opportunity to probe the physics of radio sources at very low radio frequencies. In this paper, we present a catalogue of linearly polarised sources detected at 216 MHz, using data from the Galactic and Extragalactic All-sky Murchison Widefield Array survey. Our catalogue covers the Declination range –17° to –37° and 24 h in Right Ascension, at a resolution of around 3 arcminutes. We detect 81 sources (including both a known pulsar and a new pulsar candidate) with linearly polarised flux densities in excess of 18 mJy across a survey area of approximately 6 400 deg2, corresponding to a surface density of 1 source per 79 deg2. The level of Faraday rotation measured for our sources is broadly consistent with those recovered at higher frequencies, with typically more than an order of magnitude improvement in the uncertainty compared to higher-frequency measurements. However, our catalogue is likely incomplete at low Faraday rotation measures, due to our practice of excluding sources in the region where instrumental leakage appears. The majority of sources exhibit significant depolarisation compared to higher frequencies; however, a small sub-sample repolarise at 216 MHz. We also discuss the polarisation properties of four nearby, large-angular-scale radio galaxies, with a particular focus on the giant radio galaxy ESO 422–G028, in order to explain the striking differences in polarised morphology between 216 MHz and 1.4 GHz.

The discovery of the first electromagnetic counterpart to a gravitational wave signal has generated follow-up observations by over 50 facilities world-wide, ushering in the new era of multi-messenger astronomy. In this paper, we present follow-up observations of the gravitational wave event GW170817 and its electromagnetic counterpart SSS17a/DLT17ck (IAU label AT2017gfo) by 14 Australian telescopes and partner observatories as part of Australian-based and Australian-led research programs. We report early- to late-time multi-wavelength observations, including optical imaging and spectroscopy, mid-infrared imaging, radio imaging, and searches for fast radio bursts. Our optical spectra reveal that the transient source emission cooled from approximately 6 400 K to 2 100 K over a 7-d period and produced no significant optical emission lines. The spectral profiles, cooling rate, and photometric light curves are consistent with the expected outburst and subsequent processes of a binary neutron star merger. Star formation in the host galaxy probably ceased at least a Gyr ago, although there is evidence for a galaxy merger. Binary pulsars with short (100 Myr) decay times are therefore unlikely progenitors, but pulsars like PSR B1534+12 with its 2.7 Gyr coalescence time could produce such a merger. The displacement (~2.2 kpc) of the binary star system from the centre of the main galaxy is not unusual for stars in the host galaxy or stars originating in the merging galaxy, and therefore any constraints on the kick velocity imparted to the progenitor are poor.

The Murchison Widefield Array (MWA), located in Western Australia, is one of the low-frequency precursors of the international Square Kilometre Array (SKA) project. In addition to pursuing its own ambitious science programme, it is also a testbed for wide range of future SKA activities ranging from hardware, software to data analysis. The key science programmes for the MWA and SKA require very high dynamic ranges, which challenges calibration and imaging systems. Correct calibration of the instrument and accurate measurements of source flux densities and polarisations require precise characterisation of the telescope’s primary beam. Recent results from the MWA GaLactic Extragalactic All-sky Murchison Widefield Array (GLEAM) survey show that the previously implemented Average Embedded Element (AEE) model still leaves residual polarisations errors of up to 10–20% in Stokes Q. We present a new simulation-based Full Embedded Element (FEE) model which is the most rigorous realisation yet of the MWA’s primary beam model. It enables efficient calculation of the MWA beam response in arbitrary directions without necessity of spatial interpolation. In the new model, every dipole in the MWA tile (4 × 4 bow-tie dipoles) is simulated separately, taking into account all mutual coupling, ground screen, and soil effects, and therefore accounts for the different properties of the individual dipoles within a tile. We have applied the FEE beam model to GLEAM observations at 200–231 MHz and used false Stokes parameter leakage as a metric to compare the models. We have determined that the FEE model reduced the magnitude and declination-dependent behaviour of false polarisation in Stokes Q and V while retaining low levels of false polarisation in Stokes U.

We present techniques developed to calibrate and correct Murchison Widefield Array low-frequency (72–300 MHz) radio observations for polarimetry. The extremely wide field-of-view, excellent instantaneous (u, v)-coverage and sensitivity to degree-scale structure that the Murchison Widefield Array provides enable instrumental calibration, removal of instrumental artefacts, and correction for ionospheric Faraday rotation through imaging techniques. With the demonstrated polarimetric capabilities of the Murchison Widefield Array, we discuss future directions for polarimetric science at low frequencies to answer outstanding questions relating to polarised source counts, source depolarisation, pulsar science, low-mass stars, exoplanets, the nature of the interstellar and intergalactic media, and the solar environment.

The current generation of experiments aiming to detect the neutral hydrogen signal from the Epoch of Reionisation (EoR) is likely to be limited by systematic effects associated with removing foreground sources from target fields. In this paper, we develop a model for the compact foreground sources in one of the target fields of the MWA’s EoR key science experiment: the ‘EoR1’ field. The model is based on both the MWA’s GLEAM survey and GMRT 150 MHz data from the TGSS survey, the latter providing higher angular resolution and better astrometric accuracy for compact sources than is available from the MWA alone. The model contains 5 049 sources, some of which have complicated morphology in MWA data, Fornax A being the most complex. The higher resolution data show that 13% of sources that appear point-like to the MWA have complicated morphology such as double and quad structure, with a typical separation of 33 arcsec. We derive an analytic expression for the error introduced into the EoR two-dimensional power spectrum due to peeling close double sources as single point sources and show that for the measured source properties, the error in the power spectrum is confined to high k⊥ modes that do not affect the overall result for the large-scale cosmological signal of interest. The brightest 10 mis-modelled sources in the field contribute 90% of the power bias in the data, suggesting that it is most critical to improve the models of the brightest sources. With this hybrid model, we reprocess data from the EoR1 field and show a maximum of 8% improved calibration accuracy and a factor of two reduction in residual power in k-space from peeling these sources. Implications for future EoR experiments including the SKA are discussed in relation to the improvements obtained.

We present low-frequency spectral energy distributions of 60 known radio pulsars observed with the Murchison Widefield Array telescope. We searched the GaLactic and Extragalactic All-sky Murchison Widefield Array survey images for 200-MHz continuum radio emission at the position of all pulsars in the Australia Telescope National Facility (ATNF) pulsar catalogue. For the 60 confirmed detections, we have measured flux densities in 20 × 8 MHz bands between 72 and 231 MHz. We compare our results to existing measurements and show that the Murchison Widefield Array flux densities are in good agreement.

We describe the performance of the Boolardy Engineering Test Array, the prototype for the Australian Square Kilometre Array Pathfinder telescope. Boolardy Engineering Test Array is the first aperture synthesis radio telescope to use phased array feed technology, giving it the ability to electronically form up to nine dual-polarisation beams. We report the methods developed for forming and measuring the beams, and the adaptations that have been made to the traditional calibration and imaging procedures in order to allow BETA to function as a multi-beam aperture synthesis telescope. We describe the commissioning of the instrument and present details of Boolardy Engineering Test Array’s performance: sensitivity, beam characteristics, polarimetric properties, and image quality. We summarise the astronomical science that it has produced and draw lessons from operating Boolardy Engineering Test Array that will be relevant to the commissioning and operation of the final Australian Square Kilometre Array Path telescope.

We estimate spatial gradients in the ionosphere using the Global Positioning System and GLONASS (Russian global navigation system) observations, utilising data from multiple Global Positioning System stations in the vicinity of Murchison Radio-astronomy Observatory. In previous work, the ionosphere was characterised using a single-station to model the ionosphere as a single layer of fixed height and this was compared with ionospheric data derived from radio astronomy observations obtained from the Murchison Widefield Array. Having made improvements to our data quality (via cycle slip detection and repair) and incorporating data from the GLONASS system, we now present a multi-station approach. These two developments significantly improve our modelling of the ionosphere. We also explore the effects of a variable-height model. We conclude that modelling the small-scale features in the ionosphere that have been observed with the MWA will require a much denser network of Global Navigation Satellite System stations than is currently available at the Murchison Radio-astronomy Observatory.

We have compiled a catalogue of H ii regions detected with the Murchison Widefield Array between 72 and 231 MHz. The multiple frequency bands provided by the Murchison Widefield Array allow us identify the characteristic spectrum generated by the thermal Bremsstrahlung process in H ii regions. We detect 306 H ii regions between 260° < l < 340° and report on the positions, sizes, peak, integrated flux density, and spectral indices of these H ii regions. By identifying the point at which H ii regions transition from the optically thin to thick regime, we derive the physical properties including the electron density, ionised gas mass, and ionising photon flux, towards 61 H ii regions. This catalogue of H ii regions represents the most extensive and uniform low frequency survey of H ii regions in the Galaxy to date.

We compare first-order (refractive) ionospheric effects seen by the MWA with the ionosphere as inferred from GPS data. The first-order ionosphere manifests itself as a bulk position shift of the observed sources across an MWA field of view. These effects can be computed from global ionosphere maps provided by GPS analysis centres, namely the CODE. However, for precision radio astronomy applications, data from local GPS networks needs to be incorporated into ionospheric modelling. For GPS observations, the ionospheric parameters are biased by GPS receiver instrument delays, among other effects, also known as receiver DCBs. The receiver DCBs need to be estimated for any non-CODE GPS station used for ionosphere modelling. In this work, single GPS station-based ionospheric modelling is performed at a time resolution of 10 min. Also the receiver DCBs are estimated for selected Geoscience Australia GPS receivers, located at Murchison Radio Observatory, Yarragadee, Mount Magnet and Wiluna. The ionospheric gradients estimated from GPS are compared with that inferred from MWA. The ionospheric gradients at all the GPS stations show a correlation with the gradients observed with the MWA. The ionosphere estimates obtained using GPS measurements show promise in terms of providing calibration information for the MWA.

GLEAM, the GaLactic and Extragalactic All-sky MWA survey, is a survey of the entire radio sky south of declination + 25° at frequencies between 72 and 231 MHz, made with the MWA using a drift scan method that makes efficient use of the MWA’s very large field-of-view. We present the observation details, imaging strategies, and theoretical sensitivity for GLEAM. The survey ran for two years, the first year using 40-kHz frequency resolution and 0.5-s time resolution; the second year using 10-kHz frequency resolution and 2 s time resolution. The resulting image resolution and sensitivity depends on observing frequency, sky pointing, and image weighting scheme. At 154 MHz, the image resolution is approximately 2.5 × 2.2/cos (δ + 26.7°) arcmin with sensitivity to structures up to ~ 10° in angular size. We provide tables to calculate the expected thermal noise for GLEAM mosaics depending on pointing and frequency and discuss limitations to achieving theoretical noise in Stokes I images. We discuss challenges, and their solutions, that arise for GLEAM including ionospheric effects on source positions and linearly polarised emission, and the instrumental polarisation effects inherent to the MWA’s primary beam.