We have found an unexpected class of astronomical objects which have not previously been reported, 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, and do not seem to correspond to any known type of object. We speculate that they may represent a spherical shock wave from an extra-galactic transient event, or the outflow, or a remnant, from a radio galaxy viewed end-on.

We present the full panchromatic afterglow light curve data of GW170817, including new radio data as well as archival optical and X-ray data, between 0.5 and 940 days post-merger. By compiling all archival data, and reprocessing a subset of it, we have ensured that the panchromatic dataset is uniform and therefore immune to the differences in data processing or flux determination methods used by different groups. Simple power-law fits to the uniform afterglow light curve indicate a $t^{0.86\pm0.04}$ rise, a $t^{-1.90\pm0.12}$ decline, and a peak occurring at $155\pm4$ days. The afterglow is optically thin throughout its evolution, consistent with a single spectral index ($-0.569\pm0.002$) across all epochs. This gives a precise and updated estimate of the electron power-law index, $p=2.138\pm0.004$. By studying the diffuse X-ray emission from the host galaxy, we place a conservative upper limit on the hot ionized ISM density, $<$0.01 cm$^{-3}$, consistent with previous afterglow studies. Using the late-time afterglow data we rule out any long-lived neutron star remnant having magnetic field strength between 10$^{10.4}$ G and 10$^{16}$ G. Our fits to the afterglow data using an analytical model that includes VLBI proper motion from Mooley et al (2018), and a structured jet model that ignores the proper motion, indicates that the proper motion measurement needs to be considered while seeking an accurate estimate of the viewing angle.

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 Murchison Widefield Array (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$ rad m$^{-2}$ 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.

Accurate instrumental bandpass corrections are essential for the reliable interpretation of spectral lines from targeted and survey-mode observations with radio interferometers. Bandpass correction is typically performed by comparing measurements of a strong calibrator source to an assumed model, typically an isolated point source. The wide field-of-view and high sensitivity of modern interferometers means that additional sources are often detected in observations of calibrators. This can introduce errors into bandpass corrections and subsequently the target data if not properly accounted for. Focusing on the standard calibrator PKS B1934-638, we perform simulations to asses this effect by constructing a wide-field sky model. The cases of ASKAP (0.7-1.9 GHz), MeerKAT (UHF: 0.58-1.05 GHz; L-band: 0.87-1.67 GHz) and Band 2 (0.95-1.76 GHz) of SKA-MID are examined. The use of a central point source model during bandpass calibration is found to impart amplitude errors into spectra measured by the precursor instruments at the ~0.2-0.5% level dropping to ~0.01% in the case of SKA-MID. This manifests itself as ripples in the source spectrum, the behaviour of which is coupled to the distribution of the array baselines, the solution interval, the primary beam size, the hour-angle of the calibration scan, as well as the weights used when imaging the target. Calibration pipelines should routinely employ complete field models for standard calibrators to remove this potentially destructive contaminant from the data, a recommendation we validate by comparing our simulation results to a MeerKAT scan of PKS B1934-638, calibrated with and without our expanded sky model.

We present results from a search for a radio transient associated with the LIGO/Virgo source S190814bv, a likely neutron star-black hole (NSBH) merger, with the Australian Square Kilometre Array Pathfinder. We imaged a $30\,{\rm deg}^2$ field at $\Delta T$=2, 9 and 33 days post-merger at a frequency of 944\,MHz, comparing them to reference images from the Rapid ASKAP Continuum Survey observed 110 days prior to the event. Each epoch of our observations covers $89\%$ of the LIGO/Virgo localisation region. We conducted an untargeted search for radio transients in this field, resulting in 21 candidates. For one of these, \object[AT2019osy]{AT2019osy}, we performed multi-wavelength follow-up and ultimately ruled out the association with S190814bv. All other candidates are likely unrelated variables, but we cannot conclusively rule them out. We discuss our results in the context of model predictions for radio emission from neutron star-black hole mergers and place constrains on the circum-merger density and inclination angle of the merger. This survey is simultaneously the first large-scale radio follow-up of an NSBH merger, and the most sensitive widefield radio transients search to-date.

We identified a highly-polarized, steep-spectrum radio source in a deep image with the Australian Square Kilometre Array Pathfinder (ASKAP) telescope at 888 MHz. After considering and rejecting a stellar origin for this source, we discovered a new millisecond pulsar (MSP) using observations from the Parkes radio telescope. This pulsar has period 2.77 ms and dispersion measure 228.27 pc/cm**3. Although this pulsar does not yet appear to be particularly remarkable, the short spin period, wide profile and high dispersion measure do make it relatively hard to discover through traditional blind periodicity searches. Over the course of several weeks we see changes in the barycentric period of this pulsar that are consistent with orbital motion in a binary system, but the properties of any binary need to be confirmed by further observations. While even a deep ASKAP survey may not identify large numbers of new MSPs compared to the existing population, it would be competitive with existing all-sky surveys and could discover interesting new MSPs at high Galactic latitude without the need for computationally-expensive all-sky periodicity searches.

We present spectropolarimetric imaging observations of the solar corona at low frequencies (80 - 240 MHz) using the Murchison Widefield Array (MWA). These images are the first of their kind, and we introduce an algorithm to mitigate an instrumental artefact by which the total intensity signal contaminates the polarimetric images due to calibration errors. We then survey the range of circular polarization (Stokes V) features detected in over 100 observing runs near solar maximum during quiescent periods. First, we detect around 700 compact polarized sources across our dataset with polarization fractions ranging from less than 0.5% to nearly 100%. These sources exhibit a positive correlation between polarization fraction and total intensity, and we interpret them as a continuum of plasma emission noise storm (Type I burst) continua sources associated with active regions. Second, we report a characteristic "bullseye" structure observed for many low-latitude coronal holes in which a central polarized component is surrounded by a ring of the opposite sense. The central component does not match the sign expected from thermal bremsstrahlung emission, and we speculate that propagation effects or an alternative emission mechanism may be responsible. Third, we show that the large-scale polarimetric structure at our lowest frequencies is reasonably well-correlated with the line-of-sight (LOS) magnetic field component inferred from a global potential field source surface (PFSS) model. The boundaries between opposite circular polarization signs are generally aligned with polarity inversion lines in the model at a height roughly corresponding to that of the radio limb. This is not true at our highest frequencies, however, where the LOS magnetic field direction and polarization sign are often not straightforwardly correlated.

Fast Radio Bursts (FRBs) are brief radio emissions from distant astronomical sources. Some are known to repeat, but most are single bursts. Non-repeating FRB observations have had insufficient positional accuracy to localize them to an individual host galaxy. We report the interferometric localization of the single pulse FRB 180924 to a position 4 kpc from the center of a luminous galaxy at redshift 0.3214. The burst has not been observed to repeat. The properties of the burst and its host are markedly different from the only other accurately localized FRB source. The integrated electron column density along the line of sight closely matches models of the intergalactic medium, indicating that some FRBs are clean probes of the baryonic component of the cosmic web.

Active M-dwarfs are known to produce bursty radio emission, and multi-wavelength studies have shown that Solar-like magnetic activity occurs in these stars. However, coherent bursts from active M-dwarfs have often been difficult to interpret in the Solar activity paradigm. We present Australian Square Array Pathfinder (ASKAP) observations of UV Ceti at a central frequency of 888 MHz. We detect several periodic, coherent pulses occurring over a timescale consistent with the rotational period of UV Ceti. The properties of the pulsed emission show that they originate from the electron cyclotron maser instability, in a cavity $\sim7$ orders of magnitude less dense than the mean coronal density at the estimated source altitude. These results confirm that auroral activity can occur in active M-dwarfs, suggesting that these stars mark the beginning of the transition from Solar-like to auroral magnetospheric behaviour. These results demonstrate the capabilities of ASKAP for detecting polarized, coherent bursts from active stars and other systems.

In this paper we discuss the radio continuum and X-ray properties of the so-far poorly studied Galactic supernova remnant (SNR) G5.9+3.1. We present the radio spectral energy distribution (SED) of the Galactic SNR G5.9+3.1 obtained with the Murchison Widefield Array (MWA). Combining these new observations with the surveys at other radio continuum frequencies, we discuss the integrated radio continuum spectrum of this particular remnant. We have also analyzed an archival XMM-Newton observation, which represents the first detection of X-ray emission from this remnant. The SNR SED is very well explained by a simple power-law relation. The synchrotron radio spectral index of G5.9+3.1, is estimated to be 0.42$\pm$0.03 and the integrated flux density at 1GHz to be around 2.7Jy. Furthermore, we propose that the identified point radio source, located centrally inside the SNR shell, is most probably a compact remnant of the supernova explosion. The shell-like X-ray morphology of G5.9+3.1 as revealed by XMM-Newton broadly matches the spatial distribution of the radio emission, where the radio-bright eastern and western rims are also readily detected in the X-ray while the radio-weak northern and southern rims are weak or absent in the X-ray. Extracted MOS1+MOS2+PN spectra from the whole SNR as well as the north, east, and west rims of the SNR are fit successfully with an optically thin thermal plasma model in collisional ionization equilibrium with a column density N_H~0.80x$10^{22}$ cm$^{-2}$ and fitted temperatures spanning the range kT~0.14-0.23keV for all of the regions. The derived electron number densities n_e for the whole SNR and the rims are also roughly comparable (ranging from ~$0.20f^{-1/2}$ cm$^{-3}$ to ~$0.40f^{-1/2}$ cm$^{-3}$, where f is the volume filling factor). We also estimate the swept-up mass of the X-ray emitting plasma associated with G5.9+3.1 to be ~$46f^{-1/2}M_{\odot}$.

We present observations of linear polarisation in the southern radio lobe of Centaurus A, conducted during commissioning of the Australian Square Kilometre Array Pathfinder (ASKAP) telescope. We used 16 antennas to observe a 30 square degree region in a single 12 hour pointing over a 240 MHz band centred on 913 MHz. Our observations achieve an angular resolution of $26\times33$ arcseconds (480 parsecs), a maximum recoverable angular scale of 30 arcminutes, and a full-band sensitivity of 85 $\muup$Jy beam$^{-1}$. The resulting maps of polarisation and Faraday rotation are amongst the most detailed ever made for radio lobes, with of order 10$^5$ resolution elements covering the source. We describe several as-yet unreported observational features of the lobe, including its detailed peak Faraday depth structure, and intricate networks of depolarised filaments. These results demonstrate the exciting capabilities of ASKAP for widefield radio polarimetry.

We present new 0.6-10 GHz observations of the binary neutron star merger GW170817 covering the period up to 300 days post-merger, taken with the Karl G. Jansky Very Large Array, the Australia Telescope Compact Array, the Giant Metrewave Radio Telescope and the MeerKAT telescope. We use these data to precisely characterize the decay phase of the late-time radio light curve. We find that the temporal decay is consistent with a power-law slope of t^-2.2, and that the transition between the power-law rise and decay is relatively sharp. Such a slope cannot be produced by a quasi-isotropic (cocoon-dominated) outflow, but is instead the classic signature of a relativistic jet. This provides strong observational evidence that GW170817 produced a successful jet, and directly demonstrates the link between binary neutron star mergers and short-hard GRBs. Using simple analytical arguments, we derive constraints on the geometry and the jet opening angle of GW170817. These results are consistent with those from our companion Very Long Baseline Interferometry (VLBI) paper, reporting superluminal motion in GW170817.

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 (MWA; 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 MWA (GLEAM) survey. Our catalogue covers the Declination range $-17^{\circ}$ to $-37^{\circ}$ and 24 hours in Right Ascension, at a resolution of around 3 arcminutes. We detect 81 sources (including both a known pulsar and new pulsar candidate) with linearly-polarised flux densities in excess of 18 mJy across a survey area of approximately 6400 square degrees, corresponding to a surface density of 1 source per 79 square degrees. 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.

We present early results from a project to measure the sky-averaged (global), redshifted $21\,$cm signal from the Epoch of Reionisation (EoR), using the Murchison Widefield Array (MWA) telescope. Because interferometers are not sensitive to a spatially-invariant global average, they cannot be used to detect this signal using standard techniques. However, lunar occultation of the radio sky imprints a spatial structure on the global signal, allowing us to measure the average brightness temperature of the patch of sky immediately surrounding the Moon. In this paper we present one night of Moon observations with the MWA between 72 - 230 MHz and verify our techniques to extract the background sky temperature from measurements of the Moon's flux density. We improve upon previous work using the lunar occultation technique by using a more sophisticated model for reflected `earthshine' and by employing image differencing to remove imaging artefacts. We leave the Moon's (constant) radio brightness temperature as a free parameter in our fit to the data and as a result, measure $T_{\rm{moon}} = 180 \pm 12 $ K and a Galactic synchrotron spectral index of $-2.64\pm0.14$, at the position of the Moon. Finally, we evaluate the prospects of the lunar occultation technique for a global EoR detection and map out a way forward for future work with the MWA.

We derive the synchrotron distribution in the Milky Way disk from HII region absorption observations over -40{\deg} < l < 40{\deg} at six frequencies of 76.2, 83.8, 91.5, 99.2, 106.9, and 114.6 MHz with the GaLactic and Extragalactic All-sky Murchison widefield array survey (GLEAM). We develop a new method of emissivity calculation by taking advantage of the Haslam et al., (1981) map and known spectral indices, which enable us to simultaneously derive the emissivity and the optical depth of HII regions at each frequency. We show our derived synchrotron emissivities based on 152 absorption features of HII regions using both the method previously adopted in the literature and our improved method. We derive the synchrotron emissivity from HII regions to the Galactic edge along the line of sight and, for the first time, derive the emissivity from HII regions to the Sun. These results provide direct information on the distribution of the Galactic magnetic field and cosmic-ray electrons for future modelling.

We present results from the first all-sky radio survey in circular polarisation. The survey uses the Murchison Widefield Array (MWA) to cover 30900 sq. deg., over declinations south of +30$^{\circ}$ and north of -86$^{\circ}$ centred at 200 MHz (over a 169-231 MHz band). We achieve a spatial resolution of approx. 3' and a typical sensitivity of 3.0 mJy PSF$^{-1}$ over most of the survey region. We demonstrate a new leakage mitigation technique that reduces the leakage from total intensity into circular polarisation by an order of magnitude. In a blind survey of the imaged region, we detect 14 pulsars in circular polarisation above a 6$\sigma$ threshold. We also detect six transient sources associated with artificial satellites. A targeted survey of 2376 pulsars within the surveyed region yielded 33 detections above $4\sigma$. Looking specifically at pulsars previously detected at 200 MHz in total intensity, this represents a 35% detection rate. We also conducted a targeted survey of 2400 known flare stars, this resulted in two tentative detections above $4\sigma$. A similar targeted search for 1506 known exoplanets in the field yielded no detections above $4\sigma$. The success of the survey suggests that similar surveys at longer wavelength bands and of deeper fields are warranted.

We present 2-9 GHz radio observations of GW170817 covering the period 125-200 days post-merger, taken with the Australia Telescope Compact Array and the Karl G. Jansky Very Large Array. Our observations demonstrate that the radio afterglow peaked at $149\pm2$ days post-merger and is now declining in flux density. We see no evidence for evolution in the radio-only spectral index, which remains consistent with optically-thin synchrotron emission connecting the radio, optical, and X-ray regimes. The peak implies a total energy in the synchrotron-emitting component of a ${\rm few}\times 10^{50}\,$erg. The temporal decay rate is most consistent with mildly- or non-relativistic material and we do not see evidence for a very energetic off-axis jet, but we cannot distinguish between a lower-energy jet and more isotropic emission.

Like the magnetised planets in our Solar System, magnetised exoplanets should emit strongly at radio wavelengths. Radio emission directly traces the planetary magnetic fields and radio detections can place constraints on the physical parameters of these features. Large comparative studies of predicted radio emission characteristics for the known population of exoplanets help to identify what physical parameters could be key for producing bright, observable radio emission. Since the last comparative study, many thousands of exoplanets have been discovered. We report new estimates for the radio flux densities and maximum emission frequencies for the current population of known exoplanets orbiting pre-main sequence and main-sequence stars with spectral types F-M. The set of exoplanets predicted to produce observable radio emission are Hot Jupiters orbiting young stars. The youth of these system predicts strong stellar magnetic fields and/or dense winds, which are key for producing bright, observable radio emission. We use a new all-sky circular polarisation Murchison Widefield Array survey to place sensitive limits on 200 MHz emission from exoplanets, with $3\sigma$ values ranging from 4.0 - 45.0 mJy. Using a targeted Giant Metre Wave Radio Telescope observing campaign, we also report a $3\sigma$ upper limit of 4.5 mJy on the radio emission from V830 Tau b, the first Hot Jupiter to be discovered orbiting a pre-main sequence star. Our limit is the first to be reported for the low-frequency radio emission from this source.

We present Very Long Baseline Interferometry observations of a faint and low-luminosity ($L_{\rm 1.4 GHz} < 10^{27}~\mbox{W Hz}^{-1}$) Gigahertz-Peaked Spectrum (GPS) and Compact Steep Spectrum (CSS) sample. We select eight sources from deep radio observations that have radio spectra characteristic of a GPS or CSS source and an angular size of $\theta \lesssim 2$ arcsec, and detect six of them with the Australian Long Baseline Array. We determine their linear sizes, and model their radio spectra using Synchrotron Self Absorption (SSA) and Free Free Absorption (FFA) models. We derive statistical model ages, based on a fitted scaling relation, and spectral ages, based on the radio spectrum, which are generally consistent with the hypothesis that GPS and CSS sources are young and evolving. We resolve the morphology of one CSS source with a radio luminosity of $10^{25}~\mbox{W Hz}^{-1}$, and find what appear to be two hotspots spanning 1.7 kpc. We find that our sources follow the turnover-linear size relation, and that both homogenous SSA and an inhomogeneous FFA model can account for the spectra with observable turnovers. All but one of the FFA models do not require a spectral break to account for the radio spectrum, while all but one of the alternative SSA and power law models do require a spectral break to account for the radio spectrum. We conclude that our low-luminosity sample is similar to brighter samples in terms of their spectral shape, turnover frequencies, linear sizes, and ages, but cannot test for a difference in morphology.

We examine data from the Murchison Widefield Array (MWA) in the frequency range 72 -- 102 MHz for a field-of-view that serendipitously contained the interstellar object 'Oumuamua on 2017 November 28. Observations took place with time resolution of 0.5 s and frequency resolution of 10 kHz. %This observation was undertaken for another purpose but due to the MWA's extremely large field-of-view, 'Oumuamua was serendipitously observed simultaneously. Based on the interesting but highly unlikely suggestion that 'Oumuamua is an interstellar spacecraft, due to some unusual orbital and morphological characteristics, we examine our data for signals that might indicate the presence of intelligent life associated with 'Oumuamua. We searched our radio data for: 1) impulsive narrow-band signals; 2) persistent narrow-band signals; and 3) impulsive broadband signals. We found no such signals with non-terrestrial origins and make estimates of the upper limits on Equivalent Isotropic Radiated Power (EIRP) for these three cases of approximately 7 kW, 840 W, and 100 kW, respectively. These transmitter powers are well within the capabilities of human technologies, and are therefore plausible for alien civilizations. While the chances of positive detection in any given Search for Extraterrestrial Intelligence (SETI) experiment are vanishingly small, the characteristics of new generation telescopes such as the MWA (and in the future, the Square Kilometre Array) make certain classes of SETI experiment easy, or even a trivial by-product of astrophysical observations. This means that the future costs of SETI experiments are very low, allowing large target lists to partially balance the low probability of a positive detection.

We present a low-frequency, broadband polarization study of the FRII radio galaxy PKS J0636-2036 (z = 0.0551), using the Murchison Widefield Array (MWA) from 70 to 230 MHz. The northern and southern hotspots (separated by ~14.5' on the sky) are resolved by the MWA (3'.3 resolution) and both are detected in linear polarization across the full frequency range. A combination of Faraday rotation measure (RM) synthesis and broadband polarization model-fitting are used to constrain the Faraday depolarization properties of the source. For the integrated southern hotspot emission, two RM component models are strongly favoured over a single RM component, and the best-fitting model requires Faraday dispersions of approximately 0.7 and 1.2 rad/m$^2$ (with a mean RM of ~50 rad/m$^2$). High resolution imaging at 5" with the ATCA shows significant sub-structure in the southern hotspot and highlights some of the limitations in the polarization modelling of the MWA data. Based on the observed depolarization, combined with extrapolations of gas density scaling-relations for group environments, we estimate magnetic field strengths in the intergalactic medium between ~0.04 and 0.5 {\mu}G. We also comment on future prospects of detecting more polarized sources at low frequencies.

GW170817 is the first gravitational wave detection of a binary neutron star merger. It was accompanied by radiation across the electromagnetic spectrum and localized to the galaxy NGC 4993 at a distance of 40 Mpc. It has been proposed that the observed gamma-ray, X-ray and radio emission is due to an ultra-relativistic jet launched during the merger, directed away from our line of sight. The presence of such a jet is predicted from models positing neutron star mergers as the central engines driving short-hard gamma-ray bursts (SGRBs). Here we show that the radio light curve of GW170817 has no direct signature of an off-axis jet afterglow. While we cannot rule out the existence of a jet pointing elsewhere, the observed gamma-rays could not have originated from such a jet. Instead, the radio data requires a mildly relativistic wide-angle outflow moving towards us. This outflow could be the high velocity tail of the neutron-rich material dynamically ejected during the merger or a cocoon of material that breaks out when a jet transfers its energy to the dynamical ejecta. The cocoon scenario can explain the radio light curve of GW170817 as well as the gamma-rays and X-rays (possibly also ultraviolet and optical emission), and hence is the model most consistent with the observational data. Cocoons may be a ubiquitous phenomenon produced in neutron star mergers, giving rise to a heretofore unidentified population of radio, ultraviolet, X-ray and gamma-ray transients in the local universe.

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 program, it is also a testbed for wide range of future SKA activities ranging from hardware, software to data analysis. The key science programs 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 MWA (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 x 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.

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 afterglow cooled from approximately 6400K to 2100K over a 7-day 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 (about 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.

Gravitational waves have been detected from a binary neutron star merger event, GW170817. The detection of electromagnetic radiation from the same source has shown that the merger occurred in the outskirts of the galaxy NGC 4993, at a distance of 40 megaparsecs from Earth. We report the detection of a counterpart radio source that appears 16 days after the event, allowing us to diagnose the energetics and environment of the merger. The observed radio emission can be explained by either a collimated ultra-relativistic jet viewed off-axis, or a cocoon of mildly relativistic ejecta. Within 100 days of the merger, the radio light curves will distinguish between these models and very long baseline interferometry will have the capability to directly measure the angular velocity and geometry of the debris.

We have acquired radio continuum data between 70\,MHz and 48\,GHz for a sample of 19 southern starburst galaxies at moderate redshifts ($0.067 < z < 0.227$) with the aim of separating synchrotron and free-free emission components. Using a Bayesian framework we find the radio continuum is rarely characterised well by a single power law, instead often exhibiting low frequency turnovers below 500\,MHz, steepening at mid-to-high frequencies, and a flattening at high frequencies where free-free emission begins to dominate over the synchrotron emission. These higher order curvature components may be attributed to free-free absorption across multiple regions of star formation with varying optical depths. The decomposed synchrotron and free-free emission components in our sample of galaxies form strong correlations with the total-infrared bolometric luminosities. Finally, we find that without accounting for free-free absorption with turnovers between 90 to 500\,MHz the radio-continuum at low frequency ($\nu < 200$\,MHz) could be overestimated by upwards of a factor of twelve if a simple power law extrapolation is used from higher frequencies. The mean synchrotron spectral index of our sample is constrained to be $\alpha=-1.06$, which is steeper then the canonical value of $-0.8$ for normal galaxies. We suggest this may be caused by an intrinsically steeper cosmic ray distribution.

Synchrotron emission pervades the Galactic plane at low radio frequencies, originating from cosmic ray electrons interacting with the Galactic magnetic field. Using a low-frequency radio telescope, the Murchison Widefield Array (MWA), we measure the free-free absorption of this Galactic synchrotron emission by intervening HII regions along the line of sight. These absorption measurements allow us to calculate the Galactic cosmic-ray electron emissivity behind and in front of 47 detected HII regions in the region $250^\circ < l < 355^\circ$, $|b| < 2^\circ$. We find that all average emissivities between the HII regions and the Galactic edge along the line of sight ($\epsilon_b$) are in the range of 0.24$\,\,\sim\,\,$0.70$\,\,$K$\,\,$pc$^{-1}$ with a mean of 0.40$\,\,$K$\,\,$pc$^{-1}$ and a variance of 0.10$\,\,$K$\,\,$pc$^{-1}$ at 88$\,\,$MHz. Our best model, the Two-circle model, divides the Galactic disk into three regions using two circles centring on the Galactic centre. It shows a high emissivity region near the Galactic centre, a low emissivity region near the Galactic edge, and a medium emissivity region between these two regions, contrary to the trend found by previous studies.

We present techniques developed to calibrate and correct Murchison Widefield Array (MWA) 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 MWA provides enable instrumental calibration, removal of instrumental artefacts, and correction for ionospheric Faraday rotation through imaging techniques. With the demonstrated polarimetric capabilities of the MWA, 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 5049 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_\bot$ modes that do not affect the overall result for the large-scale cosmological signal of interest. The brightest ten 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.

Many radio galaxies show the presence of dense and dusty gas near the active nucleus. This can be traced by both 21cm HI absorption and soft X-ray absorption, offering new insight into the physical nature of the circumnuclear medium of these distant galaxies. To better understand this relationship, we investigate soft X-ray absorption as an indicator for the detection of associated HI absorption, as part of preparation for the First Large Absorption Survey in HI (FLASH) to be undertaken with the Australian Square Kilometre Array Pathfinder (ASKAP). We present the results of our pilot study using the Boolardy Engineering Test Array, a precursor to ASKAP, to search for new absorption detections in radio sources brighter than 1 Jy that also feature soft X-ray absorption. Based on this pilot survey, we detected HI absorption towards the radio source PKS 1657-298 at a redshift of z = 0.42. This source also features the highest X-ray absorption ratio of our pilot sample by a factor of 3, which is consistent with our general findings that X-ray absorption predicates the presence of dense neutral gas. By comparing the X-ray properties of AGN with and without detection of HI absorption at radio wavelengths, we find that X-ray hardness ratio and HI absorption optical depth are correlated at a statistical significance of 4.71{\sigma}. We conclude by considering the impact of these findings on future radio and X-ray absorption studies.

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

We present new radio continuum observations of NGC253 from the Murchison Widefield Array at frequencies between 76 and 227 MHz. We model the broadband radio spectral energy distribution for the total flux density of NGC253 between 76 MHz and 11 GHz. The spectrum is best described as a sum of central starburst and extended emission. The central component, corresponding to the inner 500pc of the starburst region of the galaxy, is best modelled as an internally free-free absorbed synchrotron plasma, with a turnover frequency around 230 MHz. The extended emission component of the NGC253 spectrum is best described as a synchrotron emission flattening at low radio frequencies. We find that 34% of the extended emission (outside the central starburst region) at 1 GHz becomes partially absorbed at low radio frequencies. Most of this flattening occurs in the western region of the SE halo, and may be indicative of synchrotron self-absorption of shock re-accelerated electrons or an intrinsic low-energy cut off of the electron distribution. Furthermore, we detect the large-scale synchrotron radio halo of NGC253 in our radio images. At 154 - 231 MHz the halo displays the well known X-shaped/horn-like structure, and extends out to ~8kpc in z-direction (from major axis).

We present a search for the synchrotron emission from the synchrotron cosmic web by cross correlating 180MHz radio images from the Murchison Widefield Array with tracers of large scale structure (LSS). We use two versions of the radio image covering $21.76\times 21.76$ degrees with point sources brighter than 0.05 Jy subtracted, with and without filtering of Galactic emission. As tracers of the LSS we use the Two-Micron-All-Sky-Survey (2MASS) and the Widefield InfraRed Explorer (WISE) redshift catalogues to produce galaxy number density maps. The cross correlation functions all show peak amplitudes at zero degrees, decreasing with varying slopes towards zero correlation over a range of one degree. The cross correlation signals include components from point source, Galactic, and extragalactic diffuse emission. We use models of the diffuse emission from smoothing the density maps with Gaussians of sizes 1-4 Mpc to find limits on the cosmic web components. From these models we find surface brightness 99.7 per cent upper limits in the range of 0.09-2.20 mJy beam$^{-1}$ (average beam size of 2.6 arcmin), corresponding to 0.01-0.30 mJy arcmin$^{-2}$. Assuming equipartition between energy densities of cosmic rays and the magnetic field, the flux density limits translate to magnetic field strength limits of 0.03-1.98 $\mu$G, depending heavily on the spectral index. We conclude that for a 3$\sigma$ detection of 0.1 $\mu$G magnetic field strengths via cross correlations, image depths of sub-mJy to sub-$\mu$Jy are necessary. We include discussion on the treatment and effect of extragalactic point sources and Galactic emission, and next steps for building on this work.

We have detected four flares from UV Ceti at 154 MHz using the Murchison Widefield Array. The flares have flux densities between 10--65 mJy --- a factor of 100 fainter than most flares in the literature at these frequencies --- and are only detected in polarization. The circular polarized fractions are limited to $>27$% at 3$\sigma$ confidence and two of the flares exhibit polarity reversal. We suggest that these flares occur periodically on a time scale consistent with the rotational period of UV Ceti. During the brightest observed flare, we also detect significant linear polarization with polarization fraction $>18$%. Averaging the data in 6-minute, 10 MHz frequency bins we find that the flux density of these flares does not vary over the 30 MHz bandwidth of the Murchison Widefield Array, however we cannot rule out finer time-frequency structure. Using the measured flux densities for the flares, we estimate brightness temperatures between $(10^{13}-10^{14})\,$K, indicative of a coherent emission mechanism. The brightness temperature and polarization characteristics point to the electron cyclotron maser mechanism. We also calculate the flare rates given our four observed flares and compare them to flare rates for the set of M dwarf stars with known 100--200 MHz flares. Our measurement is the first for flares with intensities $<100$ mJy at 100-200 MHz.

We have studied radio haloes and relics in nine merging galaxy clusters using the Murchison Widefield Array (MWA). The images used for this study were obtained from the GaLactic and Extragalactic All-sky MWA (GLEAM) Survey which was carried out at 5 frequencies, viz. 88, 118, 154, 188 and 215 MHz. We detect diffuse radio emission in 8 of these clusters. We have estimated the spectra of haloes and relics in these clusters over the frequency range 80-1400 MHz; the first such attempt to estimate their spectra at low frequencies. The spectra follow a power law with a mean value of $\alpha = -1.13\pm0.21$ for haloes and $\alpha = -1.2\pm0.19$ for relics where, $S \propto \nu^{\alpha}$. We reclassify two of the cluster sources as radio galaxies. The low frequency spectra are thus an independent means of confirming the nature of cluster sources. Five of the nine clusters host radio haloes. For the remaining four clusters, we place upper limits on the radio powers of possible haloes in them. These upper limits are a factor of 2-20 below those expected from the $L_{\rm X}-P_{\rm 1.4}$ relation. These limits are the lowest ever obtained and the implications of these limits to the hadronic model of halo emission are discussed.

We present a sample of 1,483 sources that display spectral peaks between 72 MHz and 1.4 GHz, selected from the GaLactic and Extragalactic All-sky Murchison Widefield Array (GLEAM) survey. The GLEAM survey is the widest fractional bandwidth all-sky survey to date, ideal for identifying peaked-spectrum sources at low radio frequencies. Our peaked-spectrum sources are the low frequency analogues of gigahertz-peaked spectrum (GPS) and compact-steep spectrum (CSS) sources, which have been hypothesized to be the precursors to massive radio galaxies. Our sample more than doubles the number of known peaked-spectrum candidates, and 95% of our sample have a newly characterized spectral peak. We highlight that some GPS sources peaking above 5 GHz have had multiple epochs of nuclear activity, and demonstrate the possibility of identifying high redshift ($z > 2$) galaxies via steep optically thin spectral indices and low observed peak frequencies. The distribution of the optically thick spectral indices of our sample is consistent with past GPS/CSS samples but with a large dispersion, suggesting that the spectral peak is a product of an inhomogeneous environment that is individualistic. We find no dependence of observed peak frequency with redshift, consistent with the peaked-spectrum sample comprising both local CSS sources and high-redshift GPS sources. The 5 GHz luminosity distribution lacks the brightest GPS and CSS sources of previous samples, implying that a convolution of source evolution and redshift influences the type of peaked-spectrum sources identified below 1 GHz. Finally, we discuss sources with optically thick spectral indices that exceed the synchrotron self-absorption limit.

Using the Boolardy Engineering Test Array of the Australian Square Kilometre Array Pathfinder (ASKAP BETA), we have carried out the first $z = 0 - 1$ survey for HI and OH absorption towards the gravitationally lensed quasars PKSB1830$-$211 and MGJ0414$+$0534. Although we detected all previously reported intervening systems towards PKSB1830$-$211, in the case of MGJ0414+0534 three systems were not found, indicating that the original identifications may have been confused with radio frequency interference. Given the sensitivity of our data, we find that our detection yield is consistent with the expected frequency of intervening HI systems estimated from previous surveys for 21-cm emission in nearby galaxies and $z \sim 3$ damped Lyman $\alpha$ absorbers. We find spectral variability in the $z = 0.886$ face-on spiral galaxy towards PKSB1830$-$211, from observations undertaken with the Westerbork Synthesis Radio Telescope in 1997/1998 and ASKAP BETA in 2014/2015. The HI equivalent width varies by a few per cent over approximately yearly time-scales. This long-term spectral variability is correlated between the north-east and south-west images of the core, and with the total flux density of the source, implying that it is observationally coupled to intrinsic changes in the quasar. The absence of any detectable variability in the ratio of HI associated with the two core images is in stark contrast to the behaviour previously seen in the molecular lines. We therefore infer that coherent opaque HI structures in this galaxy are larger than the parsec-scale molecular clouds found at mm-wavelengths.

We present a search for transient and highly variable sources at low radio frequencies (150-200 MHz) that explores long timescales of 1-3 years. We conducted this search by comparing the TIFR GMRT Sky Survey Alternative Data Release 1 (TGSS ADR1) and the GaLactic and Extragalactic All-sky Murchison Widefield Array (GLEAM) survey catalogues. To account for the different completeness thresholds in the individual surveys, we searched for compact GLEAM sources above a flux density limit of 100 mJy that were not present in the TGSS ADR1; and also for compact TGSS ADR1 sources above a flux density limit of 200 mJy that had no counterpart in GLEAM. From a total sample of 234 333 GLEAM sources and 275 612 TGSS ADR1 sources in the overlap region between the two surveys, there were 99658 GLEAM sources and 38 978 TGSS ADR sources that passed our flux density cutoff and compactness criteria. Analysis of these sources resulted in three candidate transient sources. Further analysis ruled out two candidates as imaging artefacts. We analyse the third candidate and show it is likely to be real, with a flux density of 182 +/- 26 mJy at 147.5 MHz. This gives a transient surface density of rho = (6.2 +/- 6) x 10-5 deg-2 . We present initial follow-up observations and discuss possible causes for this candidate. The small number of spurious sources from this search demonstrates the high reliability of these two new low-frequency radio catalogues.

Using the Murchison Widefield Array (MWA), the low-frequency Square Kilometre Array (SKA1 LOW) precursor located in Western Australia, we have completed the GaLactic and Extragalactic All-sky MWA (GLEAM) survey, and present the resulting extragalactic catalogue, utilising the first year of observations. The catalogue covers 24,831 square degrees, over declinations south of $+30^\circ$ and Galactic latitudes outside $10^\circ$ 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--231MHz, selected from a time- and frequency- integrated image centred at 200MHz, with a resolution of $\approx 2$'. Over the catalogued region, we estimate that the catalogue is 90% complete at 170mJy, and 50% complete at 55mJy, and large areas are complete at even lower flux density levels. Its reliability is 99.97% above the detection threshold of $5\sigma$, which itself is typically 50mJy. 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%. This paper presents details of the flagging, imaging, mosaicking, and source extraction/characterisation, as well as estimates of the completeness and reliability. All source measurements and images are available online (http://www.mwatelescope.org/science/gleam-survey). This is the first in a series of publications describing the GLEAM survey results.

The Detection of redshifted 21 cm emission from the epoch of reionization (EoR) is a challenging task owing to strong foregrounds that dominate the signal. In this paper, we propose a general method, based on the delay spectrum approach, to extract HI power spectra that is applicable to tracking observations using an imaging radio interferometer (Delay Spectrum with Imaging Arrays (DSIA)). Our method is based on modelling the HI signal taking into account the impact of wide field effects such as the $w$-term which are then used as appropriate weights in cross-correlating the measured visibilities. Our method is applicable to any radio interferometer that tracks a phase center and could be utilized for arrays such as MWA, LOFAR, GMRT, PAPER and HERA. In the literature the delay spectrum approach has been implemented for near-redundant baselines using drift scan observations. In this paper we explore the scheme for non-redundant tracking arrays, and this is the first application of delay spectrum methodology to such data to extract the HI signal. We analyze 3 hours of MWA tracking data on the EoR1 field. We present both 2-dimensional ($k_\parallel,k_\perp$) and 1-dimensional (k) power spectra from the analysis. Our results are in agreement with the findings of other pipelines developed to analyse the MWA EoR data.

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 $\Delta^2 \leq 2.7 \times 10^4$ mK$^2$ 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 deep polarimetric observations at 154 MHz with the Murchison Widefield Array (MWA), covering 625 deg^2 centered on RA=0 h, Dec=-27 deg. 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 deg^2) 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 deg - 8 deg in extent) in linear polarization clearly detectable in ~2 minute snapshots, which would remain undetectable by interferometers with minimum baseline lengths >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 rad m^-2 to 10 rad m^-2 with a large fraction peaking at ~+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.

We describe the performance of the Boolardy Engineering Test Array (BETA), the prototype for the Australian Square Kilometre Array Pathfinder telescope ASKAP. BETA is the first aperture synthesis radio telescope to use phased array feed technology, giving it the ability to electronically form up to nine dual-polarization 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 BETA's performance: sensitivity, beam characteristics, polarimetric properties and image quality. We summarise the astronomical science that it has produced and draw lessons from operating BETA that will be relevant to the commissioning and operation of the final ASKAP telescope.

Pulsars are the only cosmic radio sources known to be sufficiently compact to show diffractive interstellar scintillations. Images of the variance of radio signals in both time and frequency can be used to detect pulsars in large-scale continuum surveys using the next generation of synthesis radio telescopes. This technique allows a search over the full field of view while avoiding the need for expensive pixel-by-pixel high time resolution searches. We investigate the sensitivity of detecting pulsars in variance images. We show that variance images are most sensitive to pulsars whose scintillation time-scales and bandwidths are close to the subintegration time and channel bandwidth. Therefore, in order to maximise the detection of pulsars for a given radio continuum survey, it is essential to retain a high time and frequency resolution, allowing us to make variance images sensitive to pulsars with different scintillation properties. We demonstrate the technique with Murchision Widefield Array data and show that variance images can indeed lead to the detection of pulsars by distinguishing them from other radio sources.

Detection of the Epoch of Reionization HI 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 toward 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 approximately 1400 sq. deg. 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.

Infrared-faint radio sources (IFRS) are a class of radio-loud (RL) active galactic nuclei (AGN) at high redshifts (z > 1.7) that are characterised by their relative infrared faintness, resulting in enormous radio-to-infrared flux density ratios of up to several thousand. We aim to test the hypothesis that IFRS are young AGN, particularly GHz peaked-spectrum (GPS) and compact steep-spectrum (CSS) sources that have a low frequency turnover. We use the rich radio data set available for the Australia Telescope Large Area Survey fields, covering the frequency range between 150 MHz and 34 GHz with up to 19 wavebands from different telescopes, and build radio spectral energy distributions (SEDs) for 34 IFRS. We then study the radio properties of this class of object with respect to turnover, spectral index, and behaviour towards higher frequencies. We also present the highest-frequency radio observations of an IFRS, observed with the Plateau de Bure Interferometer at 105 GHz, and model the multi-wavelength and radio-far-infrared SED of this source. We find IFRS usually follow single power laws down to observed frequencies of around 150 MHz. Mostly, the radio SEDs are steep, but we also find ultra-steep SEDs. In particular, IFRS show statistically significantly steeper radio SEDs than the broader RL AGN population. Our analysis reveals that the fractions of GPS and CSS sources in the population of IFRS are consistent with the fractions in the broader RL AGN population. We find that at least 18% of IFRS contain young AGN, although the fraction might be significantly higher as suggested by the steep SEDs and the compact morphology of IFRS. The detailed multi-wavelength SED modelling of one IFRS shows that it is different from ordinary AGN, although it is consistent with a composite starburst-AGN model with a star formation rate of 170 solar masses per year.

We estimate spatial gradients in the ionosphere using the Global Positioning System (GPS) and GLONASS (Russian global navigation system) observations, utilising data from multiple GPS stations in the vicinity of Murchison Radio-astronomy Observatory (MRO). 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 (MWA). 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 (GNSS) stations than is currently available at the MRO.

We have compiled a catalogue of HII regions detected with the Murchison Widefield Array (MWA) between 72 and 231MHz. The multiple frequency bands provided by the MWA allow us identify the characteristic spectrum generated by the thermal Bremsstrahlung process in HII regions. We detect 302 HII regions between 260 < l < 340 and report on the positions, sizes, peak, integrated flux density, and spectral indices of these HII regions. By identifying the point at which HII 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 HII regions. This catalogue of HII regions represents the most extensive and uniform low frequency survey of HII regions in the Galaxy to date.

We present Murchison Widefield Array observations of the supernova remnant (SNR) 1987A between 72 and 230 MHz, representing the lowest frequency observations of the source to date. This large lever arm in frequency space constrains the properties of the circumstellar medium created by the progenitor of SNR 1987A when it was in its red supergiant phase. As of late-2013, the radio spectrum of SNR 1987A between 72 MHz and 8.64 GHz does not show any deviation from a non-thermal power-law with a spectral index of $-0.74 \pm 0.02$. This spectral index is consistent with that derived at higher frequencies, beneath 100 GHz, and with a shock in its adiabatic phase. A spectral turnover due to free-free absorption by the circumstellar medium has to occur below 72 MHz, which places upper limits on the optical depth of $\leq$ 0.1 at a reference frequency of 72 MHz, emission measure of $\lesssim$ 13,000 cm$^{-6}$ pc, and an electron density of $\lesssim$ 110 cm$^{-3}$. This upper limit on the electron density is consistent with the detection of prompt radio emission and models of the X-ray emission from the supernova. The electron density upper limit implies that some hydrodynamic simulations derived a red supergiant mass loss rate that is too high, or a wind velocity that is too low. The mass loss rate of $\sim 5 \times 10^{-6}$ $M_{\odot}$ yr$^{-1}$ and wind velocity of 10 km s$^{-1}$ obtained from optical observations are consistent with our upper limits, predicting a current turnover frequency due to free-free absorption between 5 and 60 MHz.

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.

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). Three hours 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 hours 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 $10^4$ mK on comoving scales $k\lesssim 0.5 h$Mpc$^{-1}$. This represents the first upper limits on the $21$ cm power spectrum fluctuations at redshifts $12\lesssim z \lesssim 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.

We report the discovery of new, high-velocity narrow-line components of the OH megamaser in IRAS 20100-4156. Results from the Australian Square Kilometre Array Pathfinder (ASKAP)'s Boolardy Engineering Test Array (BETA) and the Australia Telescope Compact Array (ATCA) provide two independent measurements of the OH megamaser spectrum. We found evidence for OH megamaser clumps at $-$409 and $-$562 km/s (blue-shifted) from the systemic velocity of the galaxy, in addition to the lines previously known. The presence of such high velocities in the molecular emission from IRAS 20100$-$4156 could be explained by a ~50 pc molecular ring enclosing an approximately 3.8 billion solar mass black hole. We also discuss two alternatives, i.e. that the narrow-line masers are dynamically coupled to the wind driven by the active galactic nucleus or they are associated with two separate galactic nuclei. The comparison between the BETA and ATCA spectra provides another scientific verification of ASKAP's BETA. Our data, combined with previous measurements of the source enabled us to study the variability of the source over a twenty-six year period. The flux density of the brightest OH maser components has reduced by more than a factor of two between 1988 and 2015, whereas a secondary narrow-line component has more than doubled in the same time. Plans for high-resolution VLBI follow-up of this source are discussed, as are prospects for discovering new OH megamasers during the ASKAP early science program.

The W50/SS433 system is an unusual Galactic outflow-driven object of debatable origin. We have used the Australia Telescope Compact Array (ATCA) to observe a new 198 pointing mosaic, covering $3^\circ \times 2^\circ$, and present the highest-sensitivity full-Stokes data of W50 to date using wide-field, wide-band imaging over a 2 GHz bandwidth centred at 2.1 GHz. We also present a complementary H$\alpha$ mosaic created using the Isaac Newton Telescope Photometric H$\alpha$ Survey of the Northern Galactic Plane (IPHAS). The magnetic structure of W50 is found to be consistent with the prevailing hypothesis that the nebula is a reanimated shell-like supernova remnant (SNR), that has been re-energised by the jets from SS433. We observe strong depolarization effects that correlate with diffuse H$\alpha$ emission, likely due to spatially-varying Faraday rotation measure (RM) fluctuations of $\ge48$ to 61 rad m$^{-2}$ on scales $\le4.5$ to 6 pc. We also report the discovery of numerous, faint, H$\alpha$ filaments that are unambiguously associated with the central region of W50. These thin filaments are suggestive of a SNR's shock emission, and almost all have a radio counterpart. Furthermore, an RM-gradient is detected across the central region of W50, which we interpret as a loop magnetic field with a symmetry axis offset by $\approx90^{\circ}$ to the east-west jet-alignment axis, and implying that the evolutionary processes of both the jets and the SNR must be coupled. A separate RM-gradient is associated with the termination shock in the Eastern ear, which we interpret as a ring-like field located where the shock of the jet interacts with the circumstellar medium. Future optical observations will be able to use the new H$\alpha$ filaments to probe the kinematics of the shell of W50, potentially allowing for a definitive experiment on W50's formation history.

The GaLactic and Extragalactic All-sky MWA survey (GLEAM) is a new relatively low resolution, contiguous 72-231 MHz survey of the entire sky south of declination +25 deg. In this paper, we outline one approach to determine the relative contribution of system noise, classical confusion and sidelobe confusion in GLEAM images. An understanding of the noise and confusion properties of GLEAM is essential if we are to fully exploit GLEAM data and improve the design of future low-frequency surveys. Our early results indicate that sidelobe confusion dominates over the entire frequency range, implying that enhancements in data processing have the potential to further reduce the noise.

This paper outlines how the new GaLactic and Extragalactic All-sky MWA Survey (GLEAM, Wayth et al. 2015), observed by the Murchison Widefield Array covering the frequency range 72 - 231 MHz, allows identification of a new large, complete, sample of more than 2000 bright extragalactic radio sources selected at 151 MHz. With a flux density limit of 4 Jy this sample is significantly larger than the canonical fully-complete sample, 3CRR (Laing, Riley & Longair 1983). In analysing this small bright subset of the GLEAM survey we are also providing a first user check of the GLEAM catalogue ahead of its public release (Hurley-Walker et al. in prep). Whilst significant work remains to fully characterise our new bright source sample, in time it will provide important constraints to evolutionary behaviour, across a wide redshift and intrinsic radio power range, as well as being highly complementary to results from targeted, small area surveys.

Low-frequency radio arrays are opening a new window for the study of the sky, both to study new phenomena and to better characterize known source classes. Being flat-spectrum sources, blazars are so far poorly studied at low radio frequencies. We characterize the spectral properties of the blazar population at low radio frequency compare the radio and high-energy properties of the gamma-ray blazar population, and search for radio counterparts of unidentified gamma-ray sources. We cross-correlated the 6,100 deg^2 Murchison Widefield Array Commissioning Survey catalogue with the Roma blazar catalogue, the third catalogue of active galactic nuclei detected by Fermi-LAT, and the unidentified members of the entire third catalogue of gamma-ray sources detected by \fermilat. When available, we also added high-frequency radio data from the Australia Telescope 20 GHz catalogue. We find low-frequency counterparts for 186 out of 517 (36%) blazars, 79 out of 174 (45%) gamma-ray blazars, and 8 out of 73 (11%) gamma-ray blazar candidates. The mean low-frequency (120--180 MHz) blazar spectral index is $\langle \alpha_\mathrm{low} \rangle=0.57\pm0.02$: blazar spectra are flatter than the rest of the population of low-frequency sources, but are steeper than at $\sim$GHz frequencies. Low-frequency radio flux density and gamma-ray energy flux display a mildly significant and broadly scattered correlation. Ten unidentified gamma-ray sources have a (probably fortuitous) positional match with low radio frequency sources. Low-frequency radio astronomy provides important information about sources with a flat radio spectrum and high energy. However, the relatively low sensitivity of the present surveys still misses a significant fraction of these objects. Upcoming deeper surveys, such as the GaLactic and Extragalactic All-Sky MWA (GLEAM) survey, will provide further insight into this population.

Experiments that pursue detection of signals from the Epoch of Reionization (EoR) are relying on spectral smoothness of source spectra at low frequencies. This article empirically explores the effect of foreground spectra on EoR experiments by measuring high-resolution full-polarization spectra for the 586 brightest unresolved sources in one of the MWA EoR fields using 45 h of observation. A novel peeling scheme is used to subtract 2500 sources from the visibilities with ionospheric and beam corrections, resulting in the deepest, confusion-limited MWA image so far. The resulting spectra are found to be affected by instrumental effects, which limit the constraints that can be set on source-intrinsic spectral structure. The sensitivity and power-spectrum of the spectra are analysed, and it is found that the spectra of residuals are dominated by PSF sidelobes from nearby undeconvolved sources. We release a catalogue describing the spectral parameters for each measured source.

In this paper we present observations, simulations, and analysis demonstrating the direct connection between the location of foreground emission on the sky and its location in cosmological power spectra from interferometric redshifted 21 cm experiments. We begin with a heuristic formalism for understanding the mapping of sky coordinates into the cylindrically averaged power spectra measurements used by 21 cm experiments, with a focus on the effects of the instrument beam response and the associated sidelobes. We then demonstrate this mapping by analyzing power spectra with both simulated and observed data from the Murchison Widefield Array. We find that removing a foreground model which includes sources in both the main field-of-view and the first sidelobes reduces the contamination in high k_parallel modes by several percent relative to a model which only includes sources in the main field-of-view, with the completeness of the foreground model setting the principal limitation on the amount of power removed. While small, a percent-level amount of foreground power is in itself more than enough to prevent recovery of any EoR signal from these modes. This result demonstrates that foreground subtraction for redshifted 21 cm experiments is truly a wide-field problem, and algorithms and simulations must extend beyond the main instrument field-of-view to potentially recover the full 21 cm power spectrum.

The Boolardy Engineering Test Array is a 6 x 12 m dish interferometer and the prototype of the Australian Square Kilometre Array Pathfinder (ASKAP), equipped with the first generation of ASKAP's phased array feed (PAF) receivers. These facilitate rapid wide-area imaging via the deployment of simultaneous multiple beams within a 30 square degree field of view. By cycling the array through 12 interleaved pointing positions and using 9 digitally formed beams we effectively mimic a traditional 1 hour x 108 pointing survey, covering 150 square degrees over 711 - 1015 MHz in 12 hours of observing time. Three such observations were executed over the course of a week. We verify the full bandwidth continuum imaging performance and stability of the system via self-consistency checks and comparisons to existing radio data. The combined three epoch image has arcminute resolution and a 1-sigma thermal noise level of 375 micro-Jy per beam, although the effective noise is a factor 3 higher due to residual sidelobe confusion. From this we derive a catalogue of 3,722 discrete radio components, using the 35 percent fractional bandwidth to measure in-band spectral indices for 1,037 of them. A search for transient events reveals one significantly variable source within the survey area. The survey covers approximately two-thirds of the Spitzer South Pole Telescope Deep Field. This pilot project demonstrates the viability and potential of using PAFs to rapidly and accurately survey the sky at radio wavelengths.

Detection of the cosmological neutral hydrogen signal from the Epoch of Reionization, and estimation of its basic physical parameters, is the principal scientific aim of many current low-frequency radio telescopes. Here we describe the Cosmological HI Power Spectrum Estimator (CHIPS), an algorithm developed and implemented with data from the Murchison Widefield Array (MWA), 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 maximising 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 hours of data, we set a 2$\sigma$ upper limit on the EoR dimensionless power at $k=0.05$~h.Mpc$^{-1}$ of $\Delta_k^2<7.6\times{10^4}$~mK$^2$ in the redshift range $z=[6.2-6.6]$, consistent with previous estimates.

We use observations from the Boolardy Engineering Test Array (BETA) of the Australian Square Kilometre Array Pathfinder (ASKAP) telescope to search for transient radio sources in the field around the intermittent pulsar PSR J1107-5907. The pulsar is thought to switch between an "off" state in which no emission is detectable, a weak state and a strong state. We ran three independent transient detection pipelines on two-minute snapshot images from a 13 hour BETA observation in order to 1) study the emission from the pulsar, 2) search for other transient emission from elsewhere in the image and 3) to compare the results from the different transient detection pipelines. The pulsar was easily detected as a transient source and, over the course of the observations, it switched into the strong state three times giving a typical timescale between the strong emission states of 3.7 hours. After the first switch it remained in the strong state for almost 40 minutes. The other strong states lasted less than 4 minutes. The second state change was confirmed using observations with the Parkes radio telescope. No other transient events were found and we place constraints on the surface density of such events on these timescales. The high sensitivity Parkes observations enabled us to detect individual bright pulses during the weak state and to study the strong state over a wide observing band. We conclude by showing that future transient surveys with ASKAP will have the potential to probe the intermittent pulsar population.

We report on deep Chandra observations of the nearby broad-line radio galaxy Pictor A, which we combine with new Australia Telescope Compact Array (ATCA) observations. The new X-ray data have a factor 4 more exposure than observations previously presented and span a 15-year time baseline, allowing a detailed study of the spatial, temporal and spectral properties of the AGN, jet, hotspot and lobes. We present evidence for further time variation of the jet, though the flare that we reported in previous work remains the most significantly detected time-varying feature. We also confirm previous tentative evidence for a faint counterjet. Based on the radio through X-ray spectrum of the jet and its detailed spatial structure, and on the properties of the counterjet, we argue that inverse-Compton models can be conclusively rejected, and propose that the X-ray emission from the jet is synchrotron emission from particles accelerated in the boundary layer of a relativistic jet. For the first time, we find evidence that the bright western hotspot is also time-varying in X-rays, and we connect this to the small-scale structure in the hotspot seen in high-resolution radio observations. The new data allow us to confirm that the spectrum of the lobes is in good agreement with the predictions of an inverse-Compton model and we show that the data favour models in which the filaments seen in the radio images are predominantly the result of spatial variation of magnetic fields in the presence of a relatively uniform electron distribution.

We report the discovery of a new 21-cm HI absorption system using commissioning data from the Boolardy Engineering Test Array of the Australian Square Kilometre Array Pathfinder (ASKAP). Using the 711.5 - 1015.5 MHz band of ASKAP we were able to conduct a blind search for the 21-cm line in a continuous redshift range between $z = 0.4$ and 1.0, which has, until now, remained largely unexplored. The absorption line is detected at $z = 0.44$ towards the GHz-peaked spectrum radio source PKS B1740$-$517 and demonstrates ASKAP's excellent capability for performing a future wide-field survey for HI absorption at these redshifts. Optical spectroscopy and imaging using the Gemini-South telescope indicates that the HI gas is intrinsic to the host galaxy of the radio source. The narrow OIII emission lines show clear double-peaked structure, indicating either large-scale outflow or rotation of the ionized gas. Archival data from the \emph{XMM-Newton} satellite exhibit an absorbed X-ray spectrum that is consistent with a high column density obscuring medium around the active galactic nucleus. The HI absorption profile is complex, with four distinct components ranging in width from 5 to 300 km s$^{-1}$ and fractional depths from 0.2 to 20 per cent. In addition to systemic HI gas, in a circumnuclear disc or ring structure aligned with the radio jet, we find evidence for a possible broad outflow of neutral gas moving at a radial velocity of $v \sim 300$ km s$^{-1}$. We infer that the expanding young radio source ($t_{\rm age} \approx 2500$ yr) is cocooned within a dense medium and may be driving circumnuclear neutral gas in an outflow of $\sim$ 1 $\mathrm{M}_{\odot}$ yr$^{-1}$.

We present the third data release from the Australia Telescope Large Area Survey (ATLAS). These data combine the observations at 1.4 GHz before and after upgrades to the Australia Telescope Compact Array reaching a sensitivity of 14 microJy/beam in 3.6 deg^2 over the Chandra Deep Field South (CDFS) and of 17 microJy/beam in 2.7 deg^2 over the European Large Area ISO Survey South 1 (ELAIS-S1). We used a variety of array configurations to maximise the uv coverage resulting in a resolution of 16 by 7 arcsec in CDFS and of 12 by 8 arcsec in ELAIS-S1. After correcting for peak bias and bandwidth smearing, we find a total of 3034 radio source components above 5 sigma in CDFS, of which 514 (17 per cent) are considered to be extended. The number of components detected above 5 sigma in ELAIS-S1 is 2084, of which 392 (19 per cent) are classified as extended. The catalogues include reliable spectral indices (delta alpha < 0.2) between 1.40 and 1.71 GHz for ~350 of the brightest components.

Refraction and diffraction of incoming radio waves by the ionosphere induce time variability in the angular positions, peak amplitudes and shapes of radio sources, potentially complicating the automated cross-matching and identification of transient and variable radio sources. In this work, we empirically assess the effects of the ionosphere on data taken by the Murchison Widefield Array (MWA) radio telescope. We directly examine 51 hours of data observed over 10 nights under quiet geomagnetic conditions (global storm index Kp < 2), analysing the behaviour of short-timescale angular position and peak flux density variations of around ten thousand unresolved sources. We find that while much of the variation in angular position can be attributed to ionospheric refraction, the characteristic displacements (10-20 arcsec) at 154 MHz are small enough that search radii of 1-2 arcmin should be sufficient for cross-matching under typical conditions. By examining bulk trends in amplitude variability, we place upper limits on the modulation index associated with ionospheric scintillation of 1-3% for the various nights. For sources fainter than ~1 Jy, this variation is below the image noise at typical MWA sensitivities. Our results demonstrate that the ionosphere is not a significant impediment to the goals of time-domain science with the MWA at 154 MHz.

We present broadband observations and spectral modeling of PKS B0008-421, and identify it as an extreme gigahertz-peaked spectrum (GPS) source. PKS B0008-421 is characterized by the steepest known spectral slope below the turnover, close to the theoretical limit of synchrotron self-absorption, and the smallest known spectral width of any GPS source. Spectral coverage of the source spans from 0.118 to 22 GHz, which includes data from the Murchison Widefield Array and the wide bandpass receivers on the Australia Telescope Compact Array. We have implemented a Bayesian inference model fitting routine to fit the data with various absorption models. We find that without the inclusion of a high-frequency exponential break the absorption models can not accurately fit the data, with significant deviations above and below the peak in the radio spectrum. The addition of a high-frequency break provides acceptable spectral fits for the inhomogeneous free-free absorption and double-component synchrotron self-absorption models, with the inhomogeneous free-free absorption model statistically favored. The requirement of a high-frequency spectral break implies that the source has ceased injecting fresh particles. Additional support for the inhomogeneous free-free absorption model as being responsible for the turnover in the spectrum is given by the consistency between the physical parameters derived from the model fit and the implications of the exponential spectral break, such as the necessity of the source being surrounded by a dense ambient medium to maintain the peak frequency near the gigahertz region. The discovery of PKS B0008-421 suggests that the next generation of low radio frequency surveys could reveal a large population of GPS sources that have ceased activity, and that a portion of the ultra-steep spectrum source population could be composed of these GPS sources in a relic phase.

We compare first order (refractive) ionospheric effects seen by the Murchison Widefield Array (MWA) with the ionosphere as inferred from Global Positioning System (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 Center for Orbit Determination in Europe (CODE), using data from globally distributed GPS receivers. However, for the more accurate local ionosphere estimates required 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 Differential Code Biases (DCBs). The receiver DCBs need to be estimated for any non-CODE GPS station used for ionosphere modelling, a requirement for establishing dense GPS networks in arbitrary locations in the vicinity of the MWA. In this work, single GPS station-based ionospheric modelling is performed at a time resolution of 10 minutes. Also the receiver DCBs are estimated for selected Geoscience Australia (GA) GPS receivers, located at Murchison Radio Observatory (MRO1), Yarragadee (YAR3), Mount Magnet (MTMA) and Wiluna (WILU). The ionospheric gradients estimated from GPS are compared with the ionospheric gradients inferred from radio source position shifts observed with the 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.

We confirm our recent prediction of the "pitchfork" foreground signature in power spectra of high-redshift 21 cm measurements where the interferometer is sensitive to large-scale structure on all baselines. This is due to the inherent response of a wide-field instrument and is characterized by enhanced power from foreground emission in Fourier modes adjacent to those considered to be the most sensitive to the cosmological H I signal. In our recent paper, many signatures from the simulation that predicted this feature were validated against Murchison Widefield Array (MWA) data, but this key pitchfork signature was close to the noise level. In this paper, we improve the data sensitivity through the coherent averaging of 12 independent snapshots with identical instrument settings and provide the first confirmation of the prediction with a signal-to-noise ratio > 10. This wide-field effect can be mitigated by careful antenna designs that suppress sensitivity near the horizon. Simple models for antenna apertures that have been proposed for future instruments such as the Hydrogen Epoch of Reionization Array and the Square Kilometre Array indicate they should suppress foreground leakage from the pitchfork by ~40 dB relative to the MWA and significantly increase the likelihood of cosmological signal detection in these critical Fourier modes in the three-dimensional power spectrum.

New images from the Hubble Space Telescope of the FRII radio galaxy Pictor A reveal a previously undiscovered tidal tail, as well as a number of jet knots coinciding with a known X-ray and radio jet. The tidal tail is approximately 5" wide (3 kpc projected), starting 18" (12 kpc) from the center of Pictor A, and extends more than 90" (60 kpc). The knots are part of a jet observed to be about 4' (160 kpc) long, extending to a bright hotspot. These images are the first optical detections of this jet, and by extracting knot flux densities through three filters we set constraints on emission models. While the radio and optical flux densities are usually explained by synchrotron emission, there are several emission mechanisms which might be used to explain the X-ray flux densities. Our data rule out Doppler boosted inverse Compton scattering as a source of the high energy emission. Instead, we find that the observed emission can be well described by synchrotron emission from electrons with a low energy index ($p\sim2$) that dominates the radio band, while a high energy index ($p\sim3$) is needed for the X-ray band and the transition occurs in the optical/infrared band. This model is consistent with a continuous electron injection scenario.

We present HI imaging of the galaxy group IC 1459 carried out with six antennas of the Australian SKA Pathfinder equipped with phased-array feeds. We detect and resolve HI in eleven galaxies down to a column density of $\sim10^{20}$ cm$^{-2}$ inside a ~6 deg$^2$ field and with a resolution of ~1 arcmin on the sky and ~8 km/s in velocity. We present HI images, velocity fields and integrated spectra of all detections, and highlight the discovery of three HI clouds -- two in the proximity of the galaxy IC 5270 and one close to NGC 7418. Each cloud has an HI mass of $10^9$ M$_\odot$ and accounts for ~15% of the HI associated with its host galaxy. Available images at ultraviolet, optical and infrared wavelengths do not reveal any clear stellar counterpart of any of the clouds, suggesting that they are not gas-rich dwarf neighbours of IC 5270 and NGC 7418. Using Parkes data we find evidence of additional extended, low-column-density HI emission around IC 5270, indicating that the clouds are the tip of the iceberg of a larger system of gas surrounding this galaxy. This result adds to the body of evidence on the presence of intra-group gas within the IC 1459 group. Altogether, the HI found outside galaxies in this group amounts to several times $10^9$ M$_\odot$, at least 10% of the HI contained inside galaxies. This suggests a substantial flow of gas in and out of galaxies during the several billion years of the group's evolution.

Low-frequency, wide field-of-view (FoV) radio telescopes such as the Murchison Widefield Array (MWA) enable the ionosphere to be sampled at high spatial completeness. We present the results of the first power spectrum analysis of ionospheric fluctuations in MWA data, where we examined the position offsets of radio sources appearing in two datasets. The refractive shifts in the positions of celestial sources are proportional to spatial gradients in the electron column density transverse to the line of sight. These can be used to probe plasma structures and waves in the ionosphere. The regional (10-100 km) scales probed by the MWA, determined by the size of its FoV and the spatial density of radio sources (typically thousands in a single FoV), complement the global (100-1000 km) scales of GPS studies and local (0.01-1 km) scales of radar scattering measurements. Our data exhibit a range of complex structures and waves. Some fluctuations have the characteristics of travelling ionospheric disturbances (TIDs), while others take the form of narrow, slowly-drifting bands aligned along the Earth's magnetic field.

The separation of the faint cosmological background signal from bright astrophysical foregrounds remains one of the most daunting challenges of mapping the high-redshift intergalactic medium with the redshifted 21 cm line of neutral hydrogen. Advances in mapping and modeling of diffuse and point source foregrounds have improved subtraction accuracy, but no subtraction scheme is perfect. Precisely quantifying the errors and error correlations due to missubtracted foregrounds allows for both the rigorous analysis of the 21 cm power spectrum and for the maximal isolation of the "EoR window" from foreground contamination. We present a method to infer the covariance of foreground residuals from the data itself in contrast to previous attempts at a priori modeling. We demonstrate our method by setting limits on the power spectrum using a 3 h integration from the 128-tile Murchison Widefield Array. Observing between 167 and 198 MHz, we find at 95% confidence a best limit of Delta^2(k) < 3.7 x 10^4 mK^2 at comoving scale k = 0.18 hMpc^-1 and at z = 6.8, consistent with existing limits.

We have carried out multiwavelength observations of the near-by ($z=0.046$) rich, merging galaxy cluster Abell 3376 with the Murchison Widefield Array (MWA). As a part of the GaLactic and Extragalactic All-sky MWA survey (GLEAM), this cluster was observed at 88, 118, 154, 188 and 215 MHz. The known radio relics, towards the eastern and western peripheries of the cluster, were detected at all the frequencies. The relics, with a linear extent of $\sim$ 1 Mpc each, are separated by $\sim$ 2 Mpc. Combining the current observations with those in the literature, we have obtained the spectra of these relics over the frequency range 80 -- 1400 MHz. The spectra follow power laws, with $\alpha$ = $-1.17\pm0.06$ and $-1.37\pm0.08$ for the west and east relics, respectively ($S \propto \nu^{\alpha}$). Assuming the break frequency to be near the lower end of the spectrum we estimate the age of the relics to be $\sim$ 0.4 Gyr. No diffuse radio emission from the central regions of the cluster (halo) was detected. The upper limit on the radio power of any possible halo that might be present in the cluster is a factor of 35 lower than that expected from the radio power and X-ray luminosity correlation for cluster halos. From this we conclude that the cluster halo is very extended ($>$ 500 kpc) and/or most of the radio emission from the halo has decayed. The current limit on the halo radio power is a factor of ten lower than the existing upper limits with possible implications for models of halo formation.

GLEAM, the GaLactic and Extragalactic All-sky MWA survey, is a survey of the entire radio sky south of declination +25 deg at frequencies between 72 and 231 MHz, made with the Murchison Widefield Array (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 x 2.2/cos(DEC+26.7) arcmin with sensitivity to structures up to ~10 deg 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.

Detection of 21~cm emission of HI from the epoch of reionization, at redshifts z>6, is limited primarily by foreground emission. We investigate the signatures of wide-field measurements and an all-sky foreground model using the delay spectrum technique that maps the measurements to foreground object locations through signal delays between antenna pairs. We demonstrate interferometric measurements are inherently sensitive to all scales, including the largest angular scales, owing to the nature of wide-field measurements. These wide-field effects are generic to all observations but antenna shapes impact their amplitudes substantially. A dish-shaped antenna yields the most desirable features from a foreground contamination viewpoint, relative to a dipole or a phased array. Comparing data from recent Murchison Widefield Array observations, we demonstrate that the foreground signatures that have the largest impact on the HI signal arise from power received far away from the primary field of view. We identify diffuse emission near the horizon as a significant contributing factor, even on wide antenna spacings that usually represent structures on small scales. For signals entering through the primary field of view, compact emission dominates the foreground contamination. These two mechanisms imprint a characteristic "pitchfork" signature on the "foreground wedge" in Fourier delay space. Based on these results, we propose that selective down-weighting of data based on antenna spacing and time can mitigate foreground contamination substantially by a factor ~100 with negligible loss of sensitivity.

Ionization of the Earth's atmosphere by sunlight forms a complex, multi-layered plasma environment within the Earth's magnetosphere, the innermost layers being the ionosphere and plasmasphere. The plasmasphere is believed to be embedded with cylindrical density structures (ducts) aligned along the Earth's magnetic field, but direct evidence for these remains scarce. Here we report the first direct wide-angle observation of an extensive array of field-aligned ducts bridging the upper ionosphere and inner plasmasphere, using a novel ground-based imaging technique. We establish their heights and motions by feature-tracking and parallax analysis. The structures are strikingly organized, appearing as regularly-spaced, alternating tubes of overdensities and underdensities strongly aligned with the Earth's magnetic field. These findings represent the first direct visual evidence for the existence of such structures.

The Murchison Widefield Array (MWA) is a Square Kilometre Array (SKA) Precursor. The telescope is located at the Murchison Radio--astronomy Observatory (MRO) in Western Australia (WA). The MWA consists of 4096 dipoles arranged into 128 dual polarisation aperture arrays forming a connected element interferometer that cross-correlates signals from all 256 inputs. A hybrid approach to the correlation task is employed, with some processing stages being performed by bespoke hardware, based on Field Programmable Gate Arrays (FPGAs), and others by Graphics Processing Units (GPUs) housed in general purpose rack mounted servers. The correlation capability required is approximately 8 TFLOPS (Tera FLoating point Operations Per Second). The MWA has commenced operations and the correlator is generating 8.3 TB/day of correlation products, that are subsequently transferred 700 km from the MRO to Perth (WA) in real-time for storage and offline processing. In this paper we outline the correlator design, signal path, and processing elements and present the data format for the internal and external interfaces.

The Murchison Widefield Array (MWA) is a new low-frequency interferometric radio telescope built in Western Australia at one of the locations of the future Square Kilometre Array (SKA). We describe the automated radio-frequency interference (RFI) detection strategy implemented for the MWA, which is based on the AOFlagger platform, and present 72-231-MHz RFI statistics from 10 observing nights. RFI detection removes 1.1% of the data. RFI from digital TV (DTV) is observed 3% of the time due to occasional ionospheric or atmospheric propagation. After RFI detection and excision, almost all data can be calibrated and imaged without further RFI mitigation efforts, including observations within the FM and DTV bands. The results are compared to a previously published Low-Frequency Array (LOFAR) RFI survey. The remote location of the MWA results in a substantially cleaner RFI environment compared to LOFAR's radio environment, but adequate detection of RFI is still required before data can be analysed. We include specific recommendations designed to make the SKA more robust to RFI, including: the availability of sufficient computing power for RFI detection; accounting for RFI in the receiver design; a smooth band-pass response; and the capability of RFI detection at high time and frequency resolution (second and kHz-scale respectively).

This paper offers an electromagnetic, more specifically array theory, perspective on understanding strong instrumental polarization effects for planar low-frequency "aperture arrays" with the Murchison Widefield Array (MWA) as an example. A long-standing issue that has been seen here is significant instrumental Stokes leakage after calibration, particularly in Stokes Q at high frequencies. A simple model that accounts for inter-element mutual coupling is presented which explains the prominence of Q leakage seen when the array is scanned away from zenith in the principal planes. On these planes, the model predicts current imbalance in the X (E-W) and Y (N-S) dipoles and hence the Q leakage. Although helpful in concept, we find that this model is inadequate to explain the full details of the observation data. This finding motivates further experimentation with more rigorous models that account for both mutual coupling and embedded element patterns. Two more rigorous models are discussed: the "full" and "average" embedded element patterns. The viability of the "full" model is demonstrated by simulating current MWA practice of using a Hertzian dipole model as a Jones matrix estimate. We find that these results replicate the observed Q leakage to approximately 2 to 5%. Finally, we offer more direct indication for the level of improvement expected from upgrading the Jones matrix estimate with more rigorous models. Using the "average" embedded pattern as an estimate for the "full" model, we find that Q leakage of a few percent is achievable.

The lenticular galaxy ESO 243-49 hosts the ultraluminous X-ray source HLX-1, the best candidate intermediate mass black hole (IMBH) currently known. The environments of IMBHs remain unknown, however the proposed candidates include the nuclei of dwarf galaxies or globular clusters. Evidence at optical wavelengths points at HLX-1 being the remnant of an accreted dwarf galaxy. Here we report the Australia Telescope Compact Array radio observations of HI emission in and around ESO 243-49 searching for signatures of a recent merger event. No HI line emission is detected in ESO 243-49 with a 5$\sigma$ upper limit on the HI gas mass of a few $10^8 M_{\odot}$. A likely reason for this non-detection is the cluster environment depleting ESO 243-49's HI gas reservoir. The upper limit is consistent with an interpretation of HLX-1 as a dwarf satellite of ESO 243-49, however more sensitive observations are required for a detection. We detect ~$5 \times 10^8 M_{\odot}$ of HI gas in the peculiar spiral galaxy AM 0108-462, located at a projected distance of ~170 kpc from ESO 243-49. This amount of HI gas is ~10 times less than in spiral galaxies with similar optical and near-infrared properties in the field, strengthening the conclusion that the cluster environment indeed depletes the HI gas reservoir of these two galaxies. Here we also report observations of AM 0108-462 in several optical and near-infrared bands using the Magellan 6.5 m telescopes, and archival X-ray and ultraviolet observations with XMM-Newton and Swift. These data combined with the HI line data suggest it is likely that AM 0108-462 is experiencing a merger event.

We present new low-frequency observations of the nearby radio galaxy Fornax A at 154 MHz with the Murchison Widefield Array, microwave flux-density measurements obtained from WMAP and Planck data, and gamma-ray flux densities obtained from Fermi data. We also compile a comprehensive list of previously published images and flux-density measurements at radio, microwave and X-ray energies. A detailed analysis of the spectrum of Fornax A between 154 MHz and 1510 MHz reveals that both radio lobes have a similar spatially-averaged spectral index, and that there exists a steep-spectrum bridge of diffuse emission between the lobes. Taking the spectral index of both lobes to be the same, we model the spectral energy distribution of Fornax A across an energy range spanning eighteen orders of magnitude, to investigate the origin of the X-ray and gamma-ray emission. A standard leptonic model for the production of both the X-rays and gamma-rays by inverse-Compton scattering does not fit the multi-wavelength observations. Our results best support a scenario where the X-rays are produced by inverse-Compton scattering and the gamma-rays are produced primarily by hadronic processes confined to the filamentary structures of the Fornax A lobes.

We present follow-up radio observations of ESO 243-49 HLX-1 from 2012 using the Australia Telescope Compact Array (ATCA) and the Karl G. Jansky Very Large Array (VLA). We report the detection of radio emission at the location of HLX-1 during its hard X-ray state using the ATCA. Assuming that the `Fundamental Plane' of accreting black holes is applicable, we provide an independent estimate of the black hole mass of $M_{\rm{BH}}\leq2.8^{+7.5}_{-2.1} \times 10^{6}$ M$_{\odot}$ at 90% confidence. However, we argue that the detected radio emission is likely to be Doppler-boosted and our mass estimate is an upper limit. We discuss other possible origins of the radio emission such as being due to a radio nebula, star formation, or later interaction of the flares with the large-scale environment. None of these were found adequate. The VLA observations were carried out during the X-ray outburst. However, no new radio flare was detected, possibly due to a sparse time sampling. The deepest, combined VLA data suggests a variable radio source and we briefly discuss the properties of the previously detected flares and compare them with microquasars and active galactic nuclei.

We present the results of a survey for low frequency radio emission from 17 known exoplanetary systems with the Murchison Widefield Array. This sample includes 13 systems that have not previously been targeted with radio observations. We detected no radio emission at 154 MHz, and put 3 sigma upper limits in the range 15.2-112.5 mJy on this emission. We also searched for circularly polarised emission and made no detections, obtaining 3 sigma upper limits in the range 3.4-49.9 mJy. These are comparable with the best low frequency radio limits in the existing literature and translate to luminosity limits of between 1.2 x 10^14 W and 1.4 x 10^17 W if the emission is assumed to be 100% circularly polarised. These are the first results from a larger program to systematically search for exoplanetary emission with the MWA.

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.

The Murchison Widefield Array (MWA) is a new low frequency interferometric radio telescope, operating in the remote Murchison Radio Observatory in Western Australia. In this paper we present the first MWA observations of the well known radio relics in Abell 3667 (A3667) between 120 and 226 MHz. We clearly detect the radio relics in A3667 and present flux estimates and spectral indices for these features. The average spectral index of the north-west (NW) and south-east (SE) relics is -0.9 +/- 0.1 between 120 and 1400 MHz. We are able to resolve spatial variation in the spectral index of the NW relic from -1.7 to -0.4, which is consistent with results found at higher frequencies.

This paper describes the system architecture of a newly constructed radio telescope - the Boolardy Engineering Test Array, which is a prototype of the Australian Square Kilometre Array Pathfinder telescope. Phased array feed technology is used to form multiple simultaneous beams per antenna, providing astronomers with unprecedented survey speed. The test array described here is a 6-antenna interferometer, fitted with prototype signal processing hardware capable of forming at least 9 dual-polarisation beams simultaneously, allowing several square degrees to be imaged in a single pointed observation. The main purpose of the test array is to develop beamforming and wide-field calibration methods for use with the full telescope, but it will also be capable of limited early science demonstrations.

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 origins of radio halos, relics and bridges is still unclear, however galaxy cluster mergers seems to be an important factor. We clearly detect the North-West (NW) and South-East (SE) 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. How- ever, 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.

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 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 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.

This is the first of two papers describing the second data release (DR2) of the Australia Telescope Large Area Survey (ATLAS) at 1.4 GHz, which comprises deep wide-field observations in total intensity, linear polarization, and circular polarization over the Chandra Deep Field-South and European Large Area Infrared Space Observatory Survey-South 1 regions. DR2 improves upon the first data release by maintaining consistent data reductions across the two regions, including polarization analysis, and including differential number counts in total intensity and linear polarization. Typical DR2 sensitivities across the mosaicked multi-pointing images are 30 uJy/beam at approximately 12"x6" resolution over a combined area of 6.4 square degrees. In this paper we present detailed descriptions of our data reduction and analysis procedures, including corrections for instrumental effects such as positional variations in image sensitivity, bandwidth smearing with a non-circular beam, and polarization leakage, and application of the BLOBCAT source extractor. We present the DR2 images and catalogues of components (discrete regions of radio emission) and sources (groups of physically associated radio components). We describe new analytic methods to account for resolution bias and Eddington bias when constructing differential number counts of radio components.

This thesis reports on the application of new wide-field Very Long Baseline Interferometry (VLBI) imaging techniques using real data for the first time. These techniques are used to target three specific science areas: (i) a sub-parsec-scale study of compact radio sources in nearby starburst galaxies, (ii) a study of jet interactions in active radio galaxies, and (iii) an unbiased study of the sub-arcsecond, 90 cm sky.

We report the detection of molecular CO(1-0) gas in F00183-7111, one of the most extreme Ultra-Luminous Infrared Galaxies known, with the Australia Telescope Compact Array. We measure a redshift of 0.3292 for F00183-7111 from the CO(1-0) line and estimate the mass of the molecular gas in 00183 to be 1 $\times$ 10$^{10}$ M$_{\odot}$. We find that F00183-7111 is predominately powered by the AGN and only $\sim$14 per cent of the total luminosity is contributed by star-formation (SFR $\sim$220 M$_{\odot}$ yr$^{-1}$). We also present an optical image of F00183-7111, which shows an extension to the East. We searched for star-formation in this extension using radio continuum observations but do not detect any. This suggests that the star-formation is likely to be predominately nuclear. These observations provide additional support for a model in which the radio emission from ULIRGs is powered by an intense burst of star-formation and by a radio-loud AGN embedded in its nucleus, both triggered by a merger of gas-rich galaxies.

Relativistic jets are streams of plasma moving at appreciable fractions of the speed of light. They have been observed from stellar mass black holes ($\sim$3$-$20 solar masses, M$_\odot$) as well as supermassive black holes ($\sim$10$^6$$-$10$^9$ M$_\odot$) found in the centres of most galaxies. Jets should also be produced by intermediate mass black holes ($\sim$10$^2$$-$10$^5$ M$_\odot$), although evidence for this third class of black hole has until recently been weak. We report the detection of transient radio emission at the location of the intermediate mass black hole candidate ESO 243-49 HLX-1, which is consistent with a discrete jet ejection event. These observations also allow us to refine the mass estimate of the black hole to be between $\sim$9 $\times$10$^{3}$ M$_\odot$ and $\sim$9 $\times$10$^{4}$ M$_\odot$.

The results of multi-epoch observations of the southern starburst galaxy, NGC 253, with the Australian Long Baseline Array (LBA) at 2.3 GHz are presented. As with previous radio interferometric observations of this galaxy, no new sources were discovered. By combining the results of this survey with Very Large Array observations at higher frequencies from the literature, spectra were derived and a free-free absorption model was fitted of 20 known sources in NGC 253. The results were found to be consistent with previous studies. The supernova remnant, 5.48-43.3, was imaged with the highest sensitivity and resolution to date, revealing a two-lobed morphology. Comparisons with previous observations of similar resolution give an upper limit of 10,000 km/s for the expansion speed of this remnant. We derive a supernova rate of <0.2 per yr for the inner 300 pc using a model that improves on previous methods by incorporating an improved radio supernova peak luminosity distribution and by making use of multi-wavelength radio data spanning 21 years. A star formation rate of SFR(M>5 M_solar)<4.9 M_solar per yr was also estimated using the standard relation between supernova and star-formation rates. Our improved estimates of supernova and star-formation rates are consistent with studies at other wavelengths. The results of our study point to the possible existence of a small population of undetected supernova remnants, suggesting a low rate of radio supernovae production in NGC 253.

The Murchison Widefield Array (MWA) is a new low frequency interferomeric radio telescope. The MWA is the low frequency precursor to the Square Kilometre Array (SKA) and is the first of three SKA precursors to be operational, supporting a varied science mission ranging from the attempted detection of the Epoch of Reionisation to the monitoring of solar flares and space weather. We explore the possibility that the MWA can be used for the purposes of Space Situational Awareness (SSA). In particular we propose that the MWA can be used as an element of a passive radar facility operating in the frequency range 87.5 - 108 MHz (the commercial FM broadcast band). In this scenario the MWA can be considered the receiving element in a bi-static radar configuration, with FM broadcast stations serving as non-cooperative transmitters. The FM broadcasts propagate into space, are reflected off debris in Earth orbit, and are received at the MWA. The imaging capabilities of the MWA can be used to simultaneously detect multiple pieces of space debris, image their positions on the sky as a function of time, and provide tracking data that can be used to determine orbital parameters. Such a capability would be a valuable addition to Australian and global SSA assets, in terms of southern and eastern hemispheric coverage. We provide a feasibility assessment of this proposal, based on simple calculations and electromagnetic simulations that shows the detection of sub-metre size debris should be possible (debris radius of >0.5 m to ~1000 km altitude). We also present a proof-of-concept set of observations that demonstrate the feasibility of the proposal, based on the detection and tracking of the International Space Station via reflected FM broadcast signals originating in south-west Western Australia. These observations broadly validate our calculations and simulations.

We present a Stokes I, Q and U survey at 189 MHz with the Murchison Widefield Array 32-element prototype covering 2400 square degrees. The survey has a 15.6 arcmin angular resolution and achieves a noise level of 15 mJy/beam. We demonstrate a novel interferometric data analysis that involves calibration of drift scan data, integration through the co-addition of warped snapshot images and deconvolution of the point spread function through forward modeling. We present a point source catalogue down to a flux limit of 4 Jy. We detect polarization from only one of the sources, PMN J0351-2744, at a level of 1.8 \pm 0.4%, whereas the remaining sources have a polarization fraction below 2%. Compared to a reported average value of 7% at 1.4 GHz, the polarization fraction of compact sources significantly decreases at low frequencies. We find a wealth of diffuse polarized emission across a large area of the survey with a maximum peak of ~13 K, primarily with positive rotation measure values smaller than +10 rad/m^2. The small values observed indicate that the emission is likely to have a local origin (closer than a few hundred parsecs). There is a large sky area at 2^h30^m where the diffuse polarized emission rms is fainter than 1 K. Within this area of low Galactic polarization we characterize the foreground properties in a cold sky patch at $(\alpha,\delta) = (4^h,-27^\circ.6)$ in terms of three dimensional power spectra

The extreme ULIRG F00183-7111 has recently been found to have a radio-loud AGN with jets in its centre, representing an extreme example of the class of radio-loud AGNs buried within dusty star-forming galaxies. This source appears to be a rare example of a ULIRG glimpsed in the (presumably) brief period as it changes from "quasar mode" to "radio mode" activity. Such transition stages probably account for many of the high-redshift radio-galaxies and extreme high-redshift ULIRGs, and so this object at the relatively low redshift of 0.328 offers a rare opportunity to study this class of objects in detail. We have also detected the CO signal from this galaxy with the ATCA, and here describe the implications of this detection for future ULIRG studies.

The Ultra Luminous X-ray (ULX) source HLX-1 in the galaxy ESO 243-49 has an observed maximum unabsorbed X-ray luminosity of 1.3e42 erg/s (0.2-10.0 keV). From the conservative assumption that this value exceeds the Eddington limit by at most a factor of 10, the minimum mass is then 500 solar masses. The X-ray luminosity varies by a factor of 40 with an apparent recurrence timescale of approximately one year. This X-ray variability is associated with spectral state transitions similar to those seen in black hole X-ray binaries. Here we discuss our recent modelling of all the X-ray data for HLX-1 and show that it supports the idea that this ULX is powered by sub- and near Eddington accretion onto an intermediate mass black hole. We also present evidence for transient radio emission which is consistent with a discrete jet ejection event as well as comment on the nature of the environment around HLX-1 in light of recent Hubble Space Telescope photometry.

We analyse new Chandra observations of the nearest (D=4 Mpc) Seyfert 2 active galaxy, Circinus, and match them to pre-existing radio, infrared and optical data to study the kpc-scale emission. The proximity of Circinus allows us to observe in striking detail the structure of the radio lobes, revealing for the first time edge-brightened emission both in X-rays and radio. After considering various other possible scenarios, we show that this extended emission in Circinus is most likely caused by a jet-driven outflow, which is driving shells of strongly shocked gas into the halo of the host galaxy. In this context, we estimate Mach numbers M=2.7-3.6 and M=2.8-5.3 for the W and E shells respectively. We derive temperatures of 0.74 (+0.06, -0.05) keV and 0.8-1.8 keV for the W and E shells, and an expansion velocity of ~900-950 km/s. We estimate that the total energy (thermal and kinetic) involved in creating both shells is ~2x10^55 erg, and their age is ~10^6 years. Comparing these results with those we previously obtained for Centaurus A, NGC 3801 and Mrk 6, we show that these parameters scale approximately with the radio power of the parent AGN. The spatial coincidence between the X-ray and edge-brightened radio emission in Circinus resembles the morphology of some SNR shocks. This parallel has been expected for AGN, but has never been observed before. We investigate what underlying mechanisms both types of systems may have in common, arguing that, in Circinus, the edge-brightening in the shells may be accounted for by a B field enhancement caused by shock compression, but do not preclude some local particle acceleration. These results can be extrapolated to other low-power systems, particularly those with late type hosts.

Star forming galaxies are thought to dominate the sub-mJy radio population, but recent work has shown that low luminosity AGN can still make a significant contribution to the faint radio source population. Spectral indices are an important tool for understanding the emission mechanism of the faint radio sources. We have observed the extended Chandra Deep Field South at 5.5 GHz using a mosaic of 42 pointings with the Australia Telescope Compact Array (ATCA). Our image reaches an almost uniform sensitivity of ~12 microJy rms over 0.25 deg^2 with a restoring beam of 4.9 x 2.0 arcsec, making it one of the deepest 6cm surveys to date. We present the 5.5 GHz catalogue and source counts from this field. We take advantage of the large amounts of ancillary data in this field to study the 1.4 to 5.5 GHz spectral indices of the sub-mJy population. For the full 5.5 GHz selected sample we find a flat median spectral index, alpha_med = -0.40, which is consistent with previous results. However, the spectral index appears to steepen at the faintest flux density levels (S_{5.5 GHz} < 0.1 mJy), where alpha_med = -0.68. We performed stacking analysis of the faint 1.4 GHz selected sample (40 < S_{1.4 GHz} < 200 microJy) and also find a steep average spectral index, alpha = -0.8, consistent with synchrotron emission. We find a weak trend of steepening spectral index with redshift. Several young AGN candidates are identified using spectral indices, suggesting Gigahertz Peaked Spectrum (GPS) sources are as common in the mJy population as they are at Jy levels.

F00183-7111 is one of the most extreme Ultra-Luminous Infrared Galaxies known. Here we present a VLBI image which shows that F00183-7111 is powered by a combination of a radio-loud Active Galactic Nucleus surrounded by vigorous starburst activity. Although already radio-loud, the quasar jets are only 1.7 kpc long, boring through the dense gas and starburst activity that confine them. We appear to be witnessing this remarkable source in the brief transition period between merging starburst and radio-loud "quasar-mode" accretion.

EMU is a wide-field radio continuum survey planned for the new Australian Square Kilometre Array Pathfinder (ASKAP) telescope. The primary goal of EMU is to make a deep (rms ~ 10 microJy/beam) radio continuum survey of the entire Southern Sky at 1.3 GHz, extending as far North as +30 degrees declination, with a resolution of 10 arcsec. EMU is expected to detect and catalogue about 70 million galaxies, including typical star-forming galaxies up to z~1, powerful starbursts to even greater redshifts, and AGNs to the edge of the visible Universe. It will undoubtedly discover new classes of object. This paper defines the science goals and parameters of the survey, and describes the development of techniques necessary to maximise the science return from EMU.

The Parkes pulsar data archive currently provides access to 144044 data files obtained from observations carried out at the Parkes observatory since the year 1991. Around 10^5 files are from surveys of the sky, the remainder are observations of 775 individual pulsars and their corresponding calibration signals. Survey observations are included from the Parkes 70cm and the Swinburne Intermediate Latitude surveys. Individual pulsar observations are included from young pulsar timing projects, the Parkes Pulsar Timing Array and from the PULSE@Parkes outreach program. The data files and access methods are compatible with Virtual Observatory protocols. This paper describes the data currently stored in the archive and presents ways in which these data can be searched and downloaded.

Here we describe the Compact Array Broadband Backend (CABB) and present first results obtained with the upgraded Australia Telescope Compact Array (ATCA). The 16-fold increase in observing bandwidth, from 2 x 128 MHz to 2 x 2048 MHz, high bit sampling, and addition of 16 zoom windows (each divided into a further 2048 channels) provide major improvements for all ATCA observations. The benefits of the new system are: (1) hugely increased radio continuum and polarization sensitivity as well as image fidelity, (2) substantially improved capability to search for and map emission and absorption lines over large velocity ranges, (3) simultaneous multi-line and continuum observations, (4) increased sensitivity, survey speed and dynamic range due to high-bit sampling, and (5) high velocity resolution, while maintaining full polarization output. The new CABB system encourages all observers to make use of both spectral line and continuum data to achieve their full potential. Given the dramatic increase of the ATCA capabilities in all bands (ranging from 1.1 to 105 GHz) CABB enables scientific projects that were not feasible before the upgrade, such as simultaneous observations of multiple spectral lines, on-the-fly mapping, fast follow-up of radio transients (e.g., the radio afterglow of new supernovae) and maser observations at high velocity resolution and full polarization. The first science results presented here include wide-band spectra, high dynamic-range images, and polarization measurements, highlighting the increased capability and discovery potential of the ATCA.

One of our closest neighbours, the Andromeda Galaxy (M31) has been the subject of numerous large area studies across the entire spectrum, but so far full-disk radio surveys have been conducted only at low resolution. The new wide-field capabilities of the DiFX software correlator present the possibility of imaging the entire primary beam of a VLBI array, thus enabling a high resolution wide-field study of the entire galaxy. Using the VLBA and EVN, pilot observations of M31 have been carried out with the aim of using these new wide-field techniques to characterise the population of compact components at VLBI resolution both within and behind one of our nearest neighbours. This contribution describes the observations carried out, the preliminary processing and first results.

Wide-field surveys are a commonly-used method for studying thousands of objects simultaneously, to investigate, e.g., the joint evolution of star-forming galaxies and active galactic nuclei. VLBI observations can yield valuable input to such studies because they are able to identify AGN. However, VLBI observations of large swaths of the sky are impractical using standard methods, because the fields of view of VLBI observations are of the order of 10" or less. We have embarked on a project to carry out Very Long Baseline Array (VLBA) observations of all 96 known radio sources in one of the best-studied areas in the sky, the Chandra Deep Field South (CDFS). The challenge was to develop methods which could significantly reduce the amount of observing (and post-processing) time. We have developed an extension to the DiFX software correlator which allows one to correlate hundreds of positions within the primary beams. This extension enabled us to target many sources, at full resolution and high sensitivity, using only a small amount of observing time. The combination of wide fields-of-view and high sensitivity across the field in this survey is unprecedented. We have observed with the VLBA a single pointing containing the Chandra Deep Field South, in which 96 radio sources were known from previous observations with the ATCA. From our input sample, 20 were detected with the VLBA. The majority of objects have flux densities in agreement with arcsec-scale observations, implying that their radio emission comes from very small regions. One VLBI-detected object had earlier been classified as a star-forming galaxy. Comparing the VLBI detections to sources found in sensitive, co-located X-ray observations we find that X-ray detections are not a good indicator for VLBI detections. Wide-field VLBI survey science is now coming of age.

Infrared-faint radio sources (IFRS) are objects that have flux densities of several mJy at 1.4GHz, but that are invisible at 3.6um when using sensitive Spitzer observations with uJy sensitivities. Their nature is unclear and difficult to investigate since they are only visible in the radio. High-resolution radio images and comprehensive spectral coverage can yield constraints on the emission mechanisms of IFRS and can give hints to similarities with known objects. We imaged a sample of 17 IFRS at 4.8GHz and 8.6GHz with the Australia Telescope Compact Array to determine the structures on arcsecond scales. We added radio data from other observing projects and from the literature to obtain broad-band radio spectra. We find that the sources in our sample are either resolved out at the higher frequencies or are compact at resolutions of a few arcsec, which implies that they are smaller than a typical galaxy. The spectra of IFRS are remarkably steep, with a median spectral index of -1.4 and a prominent lack of spectral indices larger than -0.7. We also find that, given the IR non-detections, the ratio of 1.4GHz flux density to 3.6um flux density is very high, and this puts them into the same regime as high-redshift radio galaxies. The evidence that IFRS are predominantly high-redshift sources driven by active galactic nuclei (AGN) is strong, even though not all IFRS may be caused by the same phenomenon. Compared to the rare and painstakingly collected high-redshift radio galaxies, IFRS appear to be much more abundant, but less luminous, AGN-driven galaxies at similar cosmological distances.

VLBI observations are a reliable method to identify AGN, since they require high brightness temperatures for a detection to be made. However, because of the tiny fields of view it is unpractical to carry out VLBI observations of many sources using conventional methods. We used an extension of the DiFX software correlator to image with high sensitivity 96 sources in the Chandra Deep Field South, using only 9h of observing time with the VLBA. We detected 20 sources, 8 of which had not been identified as AGN at any other wavelength, despite the comprehensive coverage of this field. The lack of X-ray counterparts to 1/3 of the VLBI-detected sources, despite the sensitivity of co-located X-ray data, demonstrates that X-ray observations cannot be solely relied upon when searching for AGN activity. Surprisingly, we find that sources classified as type 1 QSOs using X-ray data are always detected, in contrast to the 10% radio-loud objects which are found in optically-selected QSOs. We present the continuation of this project with the goal to image 1450 sources in the Lockman Hole/XMM region.

A Chandra X-ray imaging observation of the jet in Pictor A showed a feature that appears to be a flare that faded between 2000 and 2002. The feature was not detected in a follow-up observation in 2009. The jet itself is over 150 kpc long and a kpc wide, so finding year-long variability is surprising. Assuming a synchrotron origin of the observed high-energy photons and a minimum energy condition for the outflow, the synchrotron loss time of the X-ray emitting electrons is of order 1200 yr, which is much longer than the observed variability timescale. This leads to the possibility that the variable X-ray emission arises from a very small sub-volume of the jet, characterized by magnetic field that is substantially larger than the average over the jet.

The first phase of the ATLAS (Australia Telescope Large Area Survey) project surveyed a total 7 square degrees down to 30 micro Jy rms at 1.4 GHz and is the largest sensitive radio survey ever attempted. We report on the scientific achievements of ATLAS to date and plans to extend the project as a path finder for the proposed EMU (Evolutionary map of the Universe) project which has been designed to use ASKAP (Australian Square Kilometre Array Pathfinder).

High angular resolution VLBI observations of Centaurus A have been undertaken that allow access to a wide field-of-view, encompassing both the well-studied pc-scale jet and the inner part of the kpc-scale jet. The VLBI observations have detected compact regions of synchrotron emission in the kpc-scale jet that coincide with three stationary features identified from previous VLA monitoring observations. Each of these stationary features is associated with strong localised X-ray emission. The VLBI results strengthen arguments made by previous authors suggesting that the stationary features may be the result of stellar objects or gas clouds traversing the jet flow, intercepting the jet and causing strong shocks. The VLBI data show that the most strongly shocked regions in these features are resolved but have extents no larger than a few pc, reducing the required mass of the typical intercepting object by a factor of ~10 relative to previous estimates, making explanations based on high mass loss stars or low density gas clouds more plausible.

The Low Frequency Array (LOFAR) is a new generation of electronic radio telescope based on aperture array technology. The telescope is being developed by ASTRON, and currently being rolled out across the Netherlands and other countries in Europe. I present the current status of the project, and its relation to high resolution instruments such as the European VLBI Network (EVN) and the Very Long Baseline Array (VLBA). In particular, I present recent VLBI results at 327 MHz associated with: (i) a shallow survey based on VLBA archive data and (ii) a deep, wide-field Global VLBI survey centred on two in-beam calibrators, B0218+357 and J0226+3421. The results suggest that there will be no shortage of relatively bright primary calibrators that remain unresolved by LOFAR even on the longest European baselines. The sky density of fainter in-beam calibrators should also be more than adequate to permit the generation of high fidelity images over a large fraction of the sky, especially in the high-band observing band (120-240 MHz). Extending LOFAR via international stations to baseline lengths of several thousand kilometres is certainly practical and should significantly enhance the scientific output and capabilities of the array.

Wide-field VLBI observations of the nearby starburst galaxy NGC 4945, obtained with the Australian LBA, have produced 2.3 GHz images over two epochs with a maximum spatial resolution of 0.3 pc. 15 sources were detected, 13 of which correspond to sources identified in higher frequency (3 cm and 12 mm) ATCA images. Four of the sources are resolved into shell-like structures ranging between 1.1 to 2.1 pc in diameter. From these data the spectra of 13 compact radio sources in NGC 4945 were modelled; nine were found to be consistent with free-free absorbed power laws and four with a simple power law spectrum. The free-free opacity is highest toward the nucleus but varies significantly throughout the nuclear region, implying that the overall structure of the ionised medium is clumpy. Of the 13 sources, 10 have steep intrinsic spectra associated with synchrotron emission from supernova remnants, the remaining sources have flat intrinsic spectra which may be associated with thermal radio emission. A non-thermal source with a jet-like morphology is detected ~1" from the assumed location of the AGN. A type II supernova (SN) rate upper limit of 15.3/yr is determined for the inner 250 pc region of the galaxy at the 95% confidence level, based on the lack of detection of new sources in observations spanning 1.9 years and a simple model for the evolution of supernova remnants. A type II SN rate of >0.1(v/1e4)/yr is implied from estimates of supernova remnant source counts, sizes and expansion rates, where v is the radial expansion velocity of the supernova remnant in km/s. A SFR of 2.4(v/1e4)<SFR(M>=5Msun)<370 Msun/yr has been estimated directly from the SN rate limits and is of the same order of magnitude as rates determined from integrated FIR and radio luminosities.

Images made with the VLBA have resolved the region in a nearby radio galaxy, Pictor A, where the relativistic jet that originates at the nucleus terminates in an interaction with the intergalactic medium, a so-called radio galaxy hot spot. This image provides the highest spatial resolution view of such an object to date (16 pc), more than three times better than previous VLBI observations of similar objects. The north-west Pictor A hot spot is resolved into a complex set of compact components, seen to coincide with the bright part of the hot spot imaged at arcsecond-scale resolution with the VLA. In addition to a comparison with VLA data, we compare our VLBA results with data from the HST and Chandra telescopes, as well as new Spitzer data. The presence of pc-scale components in the hot spot, identifying regions containing strong shocks in the fluid flow, leads us to explore the suggestion that they represent sites of synchrotron X-ray production, contributing to the integrated X-ray flux of the hot spot, along with X-rays from synchrotron self-Compton scattering. This scenario provides a natural explanation for the radio morphology of the hot spot and its integrated X-ray emission, leading to very different predictions for the higher energy X-ray spectrum compared to previous studies. From the sizes of the individual pc-scale components and their angular spread, we estimate that the jet width at the hot spot is in the range 70 - 700 pc, which is comparable to similar estimates in PKS 2153-69, 3C 205, and 4C 41.17. The lower limit in this range arises from the suggestion that the jet may dither in its direction as it passes through hot spot backflow material close to the jet termination point, creating a "dentist drill" effect on the inside of a cavity 700 pc in diameter.

We report on the first wide-field, very long baseline interferometry (VLBI) survey at 90 cm. The survey area consists of two overlapping 28 deg^2 fields centred on the quasar J0226+3421 and the gravitational lens B0218+357. A total of 618 sources were targeted in these fields, based on identifications from Westerbork Northern Sky Survey (WENSS) data. Of these sources, 272 had flux densities that, if unresolved, would fall above the sensitivity limit of the VLBI observations. A total of 27 sources were detected as far as 2 arcdegrees from the phase centre. The results of the survey suggest that at least 10% of moderately faint (S~100 mJy) sources found at 90 cm contain compact components smaller than ~0.1 to 0.3 arcsec and stronger than 10% of their total flux densities. A ~90 mJy source was detected in the VLBI data that was not seen in the WENSS and NRAO VLA Sky Survey (NVSS) data and may be a transient or highly variable source that has been serendipitously detected. This survey is the first systematic (and non-biased), deep, high-resolution survey of the low-frequency radio sky. It is also the widest field of view VLBI survey with a single pointing to date, exceeding the total survey area of previous higher frequency surveys by two orders of magnitude. These initial results suggest that new low frequency telescopes, such as LOFAR, should detect many compact radio sources and that plans to extend these arrays to baselines of several thousand kilometres are warranted.

We report on a recent 90 cm wide-field VLBI survey of two 3.1 deg^2 fields using the VLBA, Westerbork and Jodrell Bank telescopes. In-beam calibration was used to calibrate each field, the process was simplified by imaging the calibrators in DIFMAP and transferring the calibration solutions to AIPS using the newly developed DIFMAP task - cordump. We detected and imaged 13 out of the 141 sources originally detected by the low resolution (54") WENSS survey of the same two fields. The sources were detected at 7-12 sigma levels above the image noise, had total flux densities ranging between 85-1640 mJy and were between 16'-58' from the phase centre of each field. This is the first systematic (and non-biased), deep, high resolution survey of the low frequency radio sky. These initial results suggest that new instruments such as LOFAR should detect many compact radio sources and that plans to extend these arrays to baselines of several thousand kilometres are warranted.

Wide-field, very long baseline interferometry (VLBI) observations of the nearby starburst galaxy NGC 253, obtained with the Australian Long Baseline Array (LBA), have produced a 2.3 GHz image with a maximum angular resolution of 15 mas (0.3 pc). Six sources were detected, all corresponding to sources identified in higher frequency (>5 GHz) VLA images. One of the sources, supernova remnant 5.48-43.3, is resolved into a shell-like structure approximately 90 mas (1.7 pc) in diameter. From these data and data from the literature, the spectra of 20 compact radio sources in NGC 253 were modelled and found to be consistent with free-free absorbed power laws. Broadly, the free-free opacity is highest toward the nucleus but varies significantly throughout the nuclear region (tau_0 ~ 1->20), implying that the overall structure of the ionised medium is clumpy. Of the 20 sources, nine have flat intrinsic spectra associated with thermal radio emission and the remaining 11 have steep intrinsic spectra, associated with synchrotron emission from supernova remnants. A supernova rate upper limit of 2.4 yr^-1 is determined for the inner 320 pc region of the galaxy at the 95% confidence level, based on the lack of detection of new sources in observations spanning almost 17 years and a simple model for the evolution of supernova remnants. A supernova rate of >0.14 (v/10^4) yr^-1 is implied from estimates of supernova remnant source counts, sizes and expansion rates, where v is the radial expansion velocity of the supernova remnant in km s^-1. A star formation rate of 3.4 (v/10^4) < SFR(M<=5Msun) < 59 Msun yr^-1 has been estimated directly from the supernova rate limits and is of the same order of magnitude as rates determined from integrated FIR and radio luminosities.