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 an analysis of Murchison Widefield Array radio telescope data from $\omega$ Cen, possibly a stripped dwarf spheroidal galaxy core captured by our Galaxy. Recent interpretations of Fermi-LAT $\gamma$-ray data by Brown {\it et al.} (2019) and Reynoso-Cordova {\it et al.} (2019) suggest that $\omega$ Cen may contain significant Dark Matter. We utilise their best-fit Dark Matter annihilation models, and an estimate of the magnetic field strength in $\omega$ Cen, to calculate the expected radio synchrotron signal from annihilation, and show that one can usefully rule out significant parts of the magnetic field - diffusion coefficient plane using the current observational limits. Improvement by a factor of 10-100 on these limits could constrain the models even more tightly.

The entire southern sky (Declination, $\delta <$ 30 deg) has been observed using the Murchison Widefield Array (MWA), which provides radio imaging of $\sim$2-arcmin resolution at low frequencies (72-231 MHz). This is the GaLactic and Extragalactic All-sky MWA (GLEAM) Survey, and we have previously used a combination of visual inspection, cross-checks against the literature, and internal matching to identify the 'brightest' radio-sources ($S_{\mathrm{151MHz}} >$ 4 Jy) in the extragalactic catalogue (Galactic latitude, $|b| >$ 10 deg). We refer to these 1,863 sources as the GLEAM 4-Jy (G4Jy) Sample, and use radio images (of $\leq$ 45-arcsec resolution), and multi-wavelength information, to assess their morphology and identify the galaxy that is hosting the radio emission (where appropriate). Details of how to access all of the overlays used for this work are available at https://github.com/svw26/G4Jy. Alongside this we conduct further checks against the literature, which we document in this paper for individual sources. Whilst the vast majority of the G4Jy Sample are active galactic nuclei with powerful radio-jets, we highlight that it also contains a nebula, two nearby, star-forming galaxies, a cluster relic, and a cluster halo. There are also three extended sources for which we are unable to infer the mechanism that gives rise to the low-frequency emission. In the G4Jy catalogue we provide mid-infrared identifications for 86% of the sources, and flag the remainder as: having an uncertain identification (129 sources), having a faint/uncharacterised mid-infrared host (126 sources), or it being inappropriate to specify a host (2 sources). For the subset of 129 sources, there is ambiguity concerning candidate host-galaxies, and this includes four sources (B0424$-$728, B0703$-$451, 3C 198, and 3C 403.1) where we question the existing identification.

The Murchison Widefield Array (MWA) has observed the entire southern sky (Declination, $\delta <$ 30 deg) at low radio-frequencies, over the range 72-231 MHz. These observations constitute the GaLactic and Extragalactic All-sky MWA (GLEAM) Survey, and we use the extragalactic catalogue (Galactic latitude, $|b| >$ 10 deg) to define the GLEAM 4-Jy (G4Jy) Sample. This is a complete sample of the 'brightest' radio-sources ($S_{\mathrm{151MHz}} >$ 4 Jy), the majority of which are active galactic nuclei with powerful radio-jets. Crucially, low-frequency observations allow the selection of such sources in an orientation-independent way (i.e. minimising the bias caused by Doppler boosting, inherent in high-frequency surveys). We then use higher-resolution radio images, and information at other wavelengths, to morphologically classify the brightest components in GLEAM. We also conduct cross-checks against the literature, and perform internal matching, in order to improve sample completeness (which is estimated to be $>$ 95.5%). This results in a catalogue of 1,863 sources, making the G4Jy Sample over 10 times larger than that of the revised Third Cambridge Catalogue of Radio Sources (3CRR; $S_{\mathrm{178MHz}} >$ 10.9 Jy). Of these G4Jy sources, 78 are resolved by the MWA (Phase-I) synthesised beam ($\sim$2 arcmin at 200 MHz), and we label 67% of the sample as 'single', 26% as 'double', 4% as 'triple', and 3% as having 'complex' morphology at $\sim$1 GHz (45-arcsec resolution). Alongside this, our value-added catalogue provides mid-infrared source associations (subject to 6-arcsec resolution at 3.4 micron) for the radio emission, as identified through visual inspection and thorough checks against the literature. As such, the G4Jy Sample can be used as a reliable training set for cross-identification via machine-learning algorithms. [Abstract abridged for arXiv submission.]

To make a power spectrum (PS) detection of the 21 cm signal from the Epoch of Reionisation (EoR), one must avoid/subtract bright foreground sources. Sources such as Fornax A present a modelling challenge due to spatial structures spanning from arc seconds up to a degree. We compare modelling with multi-scale (MS) CLEAN components to 'shapelets', an alternative set of basis functions. We introduce a new image-based shapelet modelling package, SHAMFI. We also introduce a new CUDA simulation code (WODEN) to generate point source, Gaussian, and shapelet components into visibilities. We test performance by modelling a simulation of Fornax A, peeling the model from simulated visibilities, and producing a residual PS. We find the shapelet method consistently subtracts large-angular-scale emission well, even when the angular-resolution of the data is changed. We find that when increasing the angular-resolution of the data, the MS CLEAN model worsens at large angular-scales. When testing on real MWA data, the expected improvement is not seen in real data because of the other dominating systematics still present. Through further simulation we find the expected differences to be lower than obtainable through current processing pipelines. We conclude shapelets are worthwhile for subtracting extended galaxies, and may prove essential for an EoR detection in the future, once other systematics have been addressed.

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.

We present a calibration component for the Murchison Widefield Array All-Sky Virtual Observatory (MWA ASVO) utilising a newly developed PostgreSQL database of calibration solutions. Since its inauguration in 2013, the MWA has recorded over thirty-four petabytes of data archived at the Pawsey Supercomputing Centre. According to the MWA Data Access policy, data become publicly available eighteen months after collection. Therefore, most of the archival data are now available to the public. Access to public data was provided in 2017 via the MWA ASVO interface, which allowed researchers worldwide to download MWA uncalibrated data in standard radio astronomy data formats (CASA measurement sets or UV FITS files). The addition of the MWA ASVO calibration feature opens a new, powerful avenue for researchers without a detailed knowledge of the MWA telescope and data processing to download calibrated visibility data and create images using standard radio-astronomy software packages. In order to populate the database with calibration solutions from the last six years we developed fully automated pipelines. A near-real-time pipeline has been used to process new calibration observations as soon as they are collected and upload calibration solutions to the database, which enables monitoring of the interferometric performance of the telescope. Based on this database we present an analysis of the stability of the MWA calibration solutions over long time intervals.

The search for emission from weakly interacting massive particle (WIMP) dark matter annihilation and decay has become a multi-pronged area of research not only targeting a diverse selection of astrophysical objects, but also taking advantage of the entire electromagnetic spectrum. The decay of WIMP particles into standard model particles has been suggested as a possible channel for synchrotron emission to be detected at low radio frequencies. Here, we present the stacking analysis of a sample of 33 dwarf spheroidal (dSph) galaxies with low-frequency (72 - 231 MHz) radio images from the GaLactic and Extragalactic All-sky Murchison Widefield Array (GLEAM) survey. We produce radial surface brightness profiles of images centred upon each dSph galaxy with background radio sources masked. We remove ten fields from the stacking due to contamination from either poorly subtracted, bright radio sources or strong background gradients across the field. The remaining 23 dSph galaxies are stacked in an attempt to obtain a statistical detection of any WIMP-induced synchrotron emission in these systems. We find that the stacked radial brightness profile does not exhibit a statistically significant detection above the 95% confidence level of $\sim$1.5 mJy beam$^{-1}$. This novel technique shows the potential of using low-frequency radio images to constrain fundamental properties of particle dark matter.

Long-term 17.6~GHz radio monitoring of the broad absorption line quasar, Mrk\,231, detected a strong flare in late 2017. This triggered four epochs of Very Long Baseline Array (VLBA) observations from 8.4~GHz to 43~GHz over a 10-week period as well as an X-ray observation with NuSTAR. This was the third campaign of VLBA monitoring that we have obtained. The 43~GHz VLBA was degraded in all epochs with only 7 of 10 antennas available in three epochs and 8 in the first epoch. However, useful results were obtained due to a fortuitous capturing of a complete short 100~mJy flare at 17.6~GHz: growth and decay. This provided useful constraints on the physical model of the ejected plasma that were not available in previous campaigns. We consider four classes of models, discrete ejections (both protonic and positronic) and jetted (protonic and positronic). The most viable model is a "dissipative bright knot" in a faint background leptonic jet with an energy flux $\sim10^{43}$ ergs/sec. Inverse Compton scattering calculations (based on these models) in the ambient quasar photon field explains the lack of a detectable increase in X-ray luminosity measured by NuSTAR. We show that the core (the bright knot) moves towards a nearby secondary at $\approx 0.97$c. The background jet is much fainter. Evidently, the high frequency VLBA core does not represent the point of origin of blazar jets, in general, and optical depth "core shift" estimates of jet points of origin can be misleading.

We present 1 to 10GHz radio continuum flux density, spectral index, polarisation and Rotation Measure (RM) images of the youngest known Galactic Supernova Remnant (SNR) G1.9+0.3, using observations from the Australia Telescope Compact Array (ATCA). We have conducted an expansion study spanning 8 epochs between 1984 and 2017, yielding results consistent with previous expansion studies of G1.9+0.3. We find a mean radio continuum expansion rate of ($0.78 \pm 0.09$) per cent year$^{-1}$ (or $\sim8900$ km s$^{-1}$ at an assumed distance of 8.5 kpc), although the expansion rate varies across the SNR perimeter. In the case of the most recent epoch between 2016 and 2017, we observe faster-than-expected expansion of the northern region. We find a global spectral index for G1.9+0.3 of $-0.81\pm0.02$ (76 MHz$-$10 GHz). Towards the northern region, however, the radio spectrum is observed to steepen significantly ($\sim -$1). Towards the two so called (east & west) "ears" of G1.9+0.3, we find very different RM values of 400-600 rad m$^{2}$ and 100-200 rad m$^{2}$ respectively. The fractional polarisation of the radio continuum emission reaches (19 $\pm$ 2)~per~cent, consistent with other, slightly older, SNRs such as Cas~A.

We have detected 27 new supernova remnants (SNRs) using a new data release of the GLEAM survey from the Murchison Widefield Array (MWA) telescope, including the lowest surface-brightness SNR ever detected, G0.1-9.7. Our method uses spectral fitting to the radio continuum to derive spectral indices for 26/27 candidates, and our low-frequency observations probe a steeper-spectrum population than previously discovered. None of the candidates have coincident Wide-field Infrared Survey Explorer mid-IR emission, further showing that the emission is non-thermal. Using pulsar associations we derive physical properties for six candidate SNRs, finding G0.1-9.7 may be younger than 10kyr. 60% of the candidates subtend areas larger than 0.2deg$^{2}$ on the sky, compared to $<25$% of previously-detected SNRs. We also make the first detection of two SNRs in the Galactic longitude range $220^\circ-240^\circ$.

We examined the latest data release from the GaLactic and Extragalactic All-sky Murchison Widefield Array (GLEAM) survey covering $345^\circ < l < 60^\circ$, $180^\circ < l < 240^\circ$, using these data and that of the Widefield Infrared Survey Explorer to follow up proposed candidate Supernova Remnants from other sources. Of the 101 candidates proposed in the region, we are able to definitively confirm ten as SNRs, tentatively confirm two as SNRs, and reclassify five as Hii regions. A further two are detectable in our images but difficult to classify; the remaining 82 are undetectable in these data. We also investigated the 18 unclassified Multi-Array Galactic Plane Imaging Survey (MAGPIS) candidate SNRs, newly confirming three as SNRs, reclassifying two as Hii regions, and exploring the unusual spectra and morphology of two others.

This work makes available a further 2,860deg$^2$ of the GLEAM survey, covering half of the accessible Galactic Plane, across twenty frequency bands sampling $72-231$MHz, with resolution $4'-2'$. Unlike previous GLEAM data releases, we used multi-scale clean to better deconvolve large-scale Galactic structure. For the Galactic longitude ranges $345^\circ < l < 67^\circ$, $180^\circ < l < 240^\circ$, we provide a compact source catalogue of 22,037 components selected from a 60-MHz bandwidth image centred at 200-MHz, with RMS noise $\approx10-20$mJy beam$^{-1}$ and position accuracy better than $2"$. The catalogue has a completeness of 50% at $\approx120$mJy, and a reliability of 99.86%. It covers Galactic latitudes $1^\circ\leq|b|\leq10^\circ$ toward the Galactic Centre and $|b|\leq10^\circ$ for other regions, and is available from Vizier; images covering $|b|\leq10^\circ$ for all longitudes are made available on the GLEAM VO server and SkyView.

The Murchison Widefield Array (MWA) is an open access telescope dedicated to studying the low frequency (80$-$300 MHz) southern sky. Since beginning operations in mid 2013, the MWA has opened a new observational window in the southern hemisphere enabling many science areas. The driving science objectives of the original design were to observe 21\,cm radiation from the Epoch of Reionisation (EoR), explore the radio time domain, perform Galactic and extragalactic surveys, and monitor solar, heliospheric, and ionospheric phenomena. All together 60$+$ programs recorded 20,000 hours producing 146 papers to date. In 2016 the telescope underwent a major upgrade resulting in alternating compact and extended configurations. Other upgrades, including digital back-ends and a rapid-response triggering system, have been developed since the original array was commissioned. In this paper we review the major results from the prior operation of the MWA, and then discuss the new science paths enabled by the improved capabilities. We group these science opportunities by the four original science themes, but also include ideas for directions outside these categories.

We report the discovery of a radio transient in the spiral galaxy M81. The transient was detected in early 2015 as part of a two-year survey of M81 made up of 12 epochs using the Karl G. Jansky Very Large Array. While undetected on 2014 September 12, the source was first detected on 2015 January 2, from which point it remained visible at an approximately constant luminosity of $L_{R} = 1.5 \pm 0.1 \times 10^{24}$erg s$^{-1}$ Hz$^{-1}$ at the observing frequency of 6 GHz for at least 2 months. Assuming this is a synchrotron event with a rise-time between 2.6 and 112 days, the peak luminosity (at equipartition) corresponds to a minimum energy of $10^{44} \lesssim E_{min} \lesssim 10^{46}$erg and jet power of $P_{min} \sim 10^{39}$erg s$^{-1}$, which are higher than most known X-ray binaries. Given its longevity, lack of short-term radio variability, and the absence of any multi-wavelength counterpart (X-ray luminosity $L_{X} \lesssim 10^{36}$erg s$^{-1}$), it does not behave like known Galactic or extragalactic X-ray binaries. The M81 transient radio properties more closely resemble the unidentified radio transient 43.78+59.3 discovered in M82, which has been suggested to be a radio nebula associated with an accreting source similar to SS 433. One possibility is that both the new M81 transient and the M82 transient may be the birth of a short-lived radio bubble associated with a discrete accretion event similar to those observed from the ULX Holmberg II X-1. However, it is not possible to rule out other identifications including long-term supernova shockwave interactions with the surrounding medium from a faint supernova or a background active galaxy.

We examine the new Galactic supernova remnant (SNR) candidate, G23.11+0.18, as seen by the Murchison Widefield Array (MWA) radio telescope. We describe the morphology of the candidate and find a spectral index of -0.63+/-0.05 in the 70-170MHz domain. A coincident TeV gamma-ray detection in High-Energy Stereoscopic System (HESS) data supports the SNR nature of G23.11+0.18 and suggests that G23.11+0.18 is accelerating particles beyond TeV energies, thus making this object a promising new cosmic ray hadron source candidate. The remnant cannot be seen in current optical, infrared and X-ray data-sets. We do find, however, a dip in CO-traced molecular gas at a line-of-sight velocity of ~85 km/s, suggesting the existence of a G23.11+0.18 progenitor wind-blown bubble. Furthermore, the discovery of molecular gas clumps at a neighbouring velocity towards HESS J1832-085 adheres to the notion that a hadronic gamma-ray production mechanism is plausible towards the north of the remnant. Based on these morphological arguments, we propose an interstellar medium association for G23.11+0.18 at a kinematic distance of 4.6+/-0.8 kpc.

Surveys for transient and variable phenomena can be confounded by the presence of extrinsic variability such as refractive interstellar scintillation (RISS). We have developed an all-sky model for RISS which can predict variability on a variety of timescales, survey locations, and observing frequencies. The model makes use of Halpha intensity maps to probe the emission measure along the line of sight, convert this to a scattering measure, and finally a scintillation strength. The model uses previously developed and long understood physics along with (indirect) measurements of the electron content and distribution within the Milky Way. We develop a set of expectations that are useful in the planning of future surveys for transient and radio variability, and demonstrate that the 1-GHz sky is a poor predictor of the variable nature of the $100$-MHz sky. Interestingly, the correlation between the incidence of variability and Galactic latitude which has been seen at 1GHz, is reversed at 100MHz. We compare the predictions of our model to a low-frequency radio survey that was conducted with the Murchison Widefield Array, and find good qualitative agreement. We discuss the implications, current limitations, and future development of the model. The model has been implemented in a Python code and is available on GitHub/Zenodo.

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 report the discovery of a new Small Magellanic Cloud Pulsar Wind Nebula (PWN) at the edge of the Supernova Remnant (SNR)-DEM S5. The pulsar powered object has a cometary morphology similar to the Galactic PWN analogs PSR B1951+32 and 'the mouse'. It is travelling supersonically through the interstellar medium. We estimate the Pulsar kick velocity to be in the range of 700-2000 km/s for an age between 28-10 kyr. The radio spectral index for this SNR PWN pulsar system is flat (-0.29 $\pm$ 0.01) consistent with other similar objects. We infer that the putative pulsar has a radio spectral index of -1.8, which is typical for Galactic pulsars. We searched for dispersion measures (DMs) up to 1000 cm/pc^3 but found no convincing candidates with a S/N greater than 8. We produce a polarisation map for this PWN at 5500 MHz and find a mean fractional polarisation of P $\sim 23$ percent. The X-ray power-law spectrum (Gamma $\sim 2$) is indicative of non-thermal synchrotron emission as is expected from PWN-pulsar system. Finally, we detect DEM S5 in Infrared (IR) bands. Our IR photometric measurements strongly indicate the presence of shocked gas which is expected for SNRs. However, it is unusual to detect such IR emission in a SNR with a supersonic bow-shock PWN. We also find a low-velocity HI cloud of $\sim 107$ km/s which is possibly interacting with DEM S5. SNR DEM S5 is the first confirmed detection of a pulsar-powered bow shock nebula found outside the Galaxy.

We present Robbie: a general work-flow for the detection and characterization of radio variability and transient events in the image domain. Robbie is designed to work in a batch processing paradigm with a modular design so that components can be swapped out or upgraded to adapt to different input data, whilst retaining a consistent and coherent methodological approach. Robbie is based on commonly used and open software, and is encapsulated in a Makefile to aid portability and reproducibility. In this paper wedescribe the methodology behind Robbie, and demonstrate its use on real and simulated data. Robbie is available on GitHub.

The GaLactic and Extragalactic All-sky MWA survey (GLEAM) is a radio continuum survey at 72-231 MHz of the whole sky south of declination +30 deg, carried out with the Murchison Widefield Array (MWA). In this paper, we derive source counts from the GLEAM data at 200, 154, 118 and 88 MHz, to a flux density limit of 50, 80, 120 and 290 mJy respectively, correcting for ionospheric smearing, incompleteness and source blending. These counts are more accurate than other counts in the literature at similar frequencies as a result of the large area of sky covered and this survey's sensitivity to extended emission missed by other surveys. At S_154MHz > 0.5 Jy, there is no evidence of flattening in the average spectral index (alpha approx. -0.8 where S proportional to nu^alpha) towards the lower frequencies. We demonstrate that the SKA Design Study (SKADS) model by Wilman et al. (2008) significantly underpredicts the observed 154 MHz GLEAM counts, particularly at the bright end. Using deeper LOFAR counts and the SKADS model, we find that sidelobe confusion dominates the thermal noise and classical confusion at nu >~ 100 MHz due to both the limited CLEANing depth and undeconvolved sources outside the field-of-view. We show that we can approach the theoretical noise limit using a more efficient and automated CLEAN algorithm.

We report on a search for low-frequency radio variability in 944 bright (> 4Jy at 154 MHz) unresolved, extragalactic radio sources monitored monthly for several years with the Murchison Widefield Array. In the majority of sources we find very low levels of variability with typical modulation indices < 5%. We detect 15 candidate low frequency variables that show significant long term variability (>2.8 years) with time-averaged modulation indices M = 3.1 - 7.1%. With 7/15 of these variable sources having peaked spectral energy distributions, and only 5.7% of the overall sample having peaked spectra, we find an increase in the prevalence of variability in this spectral class. We conclude that the variability seen in this survey is most probably a consequence of refractive interstellar scintillation and that these objects must have the majority of their flux density contained within angular diameters less than 50 milli-arcsec (which we support with multi-wavelength data). At 154 MHz we demonstrate that interstellar scintillation time-scales become long (~decades) and have low modulation indices, whilst synchrotron driven variability can only produce dynamic changes on time-scales of hundreds of years, with flux density changes less than one milli-jansky (without relativistic boosting). From this work we infer that the low frequency extra-galactic southern sky, as seen by SKA-Low, will be non-variable on time-scales shorter than one year.

Powerful radio-galaxies feature heavily in our understanding of galaxy evolution. However, when it comes to studying their properties as a function of redshift and/or environment, the most-detailed studies tend to be limited by small-number statistics. During Focus Meeting 3, on "Radio Galaxies: Resolving the AGN phenomenon", SVW presented a new sample of nearly 2,000 of the brightest radio-sources in the southern sky (Dec. $<$ 30 deg). These were observed at low radio-frequencies as part of the GaLactic and Extragalactic All-sky MWA (GLEAM) Survey, which is a continuum survey conducted using the Murchison Widefield Array (MWA). This instrument is the precursor telescope for the low-frequency component of the Square Kilometre Array, and allows us to select radio galaxies in an orientation-independent way (i.e. minimising the bias caused by Doppler boosting, inherent in high-frequency surveys). Being brighter than 4 Jy at 151 MHz, we refer to these objects as the GLEAM 4-Jy (G4Jy) Sample. The G4Jy catalogue is close to being finalised, with SVW describing how multi-wavelength data have been used to determine the morphology of the radio emission, and identify the host galaxy. In addition, the MWA's excellent spectral-coverage and sensitivity to extended/diffuse emission were highlighted. Both of these aspects are important for understanding the physical mechanisms that take place within active galaxies, and how they interact with their environment.

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.

The Earth's ionosphere refracts radio waves incident on an interferometer, resulting in shifts to the measured positions of radio sources. We present a method to smoothly remove these shifts and restore sources to their reference positions, in both the catalogue and image domains. The method is applicable to instruments and ionospheric weather such that all antennas see the same ionosphere. The method is generalisable to repairing any sparsely-sampled vector field distortion to some input data. The code is available under the Academic Free License (https://opensource.org/licenses/AFL-3.0) from https://github.com/nhurleywalker/fits_warp

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.

The low-frequency sky may reveal some of the secrets yet to be discovered and until recently, molecules had never been detected within interstellar clouds at frequencies below 700MHz (Lovas, 2003) . Following the pilot survey towards the Galactic Centre at 103--133MHz with the Murchison Widefield Array (Tremblay, 2017), we surveyed 400 degree squared centered on the Orion KL nebula from 99--170MHz. Orion is a nearby region of active star formation and known to be a chemically rich environment. In this paper, we present tentative detections of nitric oxide and its isotopologues, singularly deuterated formic acid, molecular oxygen, and several unidentified transitions. The three identified molecules are particularly interesting as laboratory experiments have suggested these molecules are precursors to the formation of amines.

We detail tentative detections of low-frequency carbon radio recombination lines from within the Orion molecular cloud complex observed at 99-129 MHz. These tentative detections include one alpha transition and one beta transition over three locations and are located within the diffuse regions of dust observed in the infrared at 100umm, the Halpha emission detected in the optical, and the synchrotron radiation observed in the radio. With these observations, we are able to study the radiation mechanism transition from collisionally pumped to radiatively pumped within the HII regions of the Orion molecular cloud complex.

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.

In the era of the SKA precursors, telescopes are producing deeper, larger images of the sky on increasingly small time-scales. The greater size and volume of images place an increased demand on the software that we use to create catalogues, and so our source finding algorithms need to evolve accordingly. In this paper we discuss some of the logistical and technical challenges that result from the increased size and volume of images that are to be analysed, and demonstrate how the Aegean source finding package has evolved to address these challenges. In particular we address the issues of source finding on spatially correlated data, and on images in which the background, noise, and point spread function, vary across the sky. We also introduce the concept of forced or priorized fitting.

We present new radio and optical images of the nearest radio galaxy Centaurus A and its host galaxy NGC 5128. We focus our investigation on the northern transition region, where energy is transported from the ~5 kpc (~5 arcmin) scales of the Northern Inner Lobe (NIL) to the ~30 kpc (~30 arcmin) scales of the Northern Middle Lobe (NML). Our Murchison Widefield Array observations at 154 MHz and our Parkes radio telescope observations at 2.3 GHz show diffuse radio emission connecting the NIL to the NML, in agreement with previous Australia Telescope Compact Array observations at 1.4 GHz. Comparison of these radio data with our widefield optical emission line images show the relationship between the NML radio emission and the ionised filaments that extend north from the NIL, and reveal a new ionised filament to the east, possibly associated with a galactic wind. Our deep optical images show clear evidence for a bipolar outflow from the central galaxy extending to intermediate scales, despite the non-detection of a southern radio counterpart to the NML. Thus, our observational overview of Centaurus A reveals a number of features proposed to be associated with AGN feedback mechanisms, often cited as likely to have significant effects in galaxy evolution models. As one of the closest galaxies to us, Centaurus A therefore provides a unique laboratory to examine feedback mechanisms in detail.

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.

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 detail and present results from a pilot study to assess the feasibility of detecting molecular lines at low radio frequencies. We observed a 400 square degree region centred on the Galactic Centre with the Murchison Widefield Array (MWA) between 103 and 133\,MHz targeting 28 known molecular species that have significant transitions. The results of this survey yield tentative detections of nitric oxide (NO) and the mercapto radical (SH). Both of these molecules appear to be associated with evolved stars.

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.

This Rescaled Subset of the Alternative Data Release 1 to the Tata Institute of Fundamental Physics Giant Metrewave Radio Telescope Sky Survey (TGSS-RSADR1) modifies the initial data release of TGSS-ADR1 (Intema et al. 2017) to bring that catalogue to the same flux scale as the extragalactic catalogue from the GaLactic and Extragalactic All-sky Murchison Widefield Array survey (GLEAM: Wayth et al. 2015; Hurley-Walker et al. 2017). In this paper we motivate the derivation of correct and complementary flux density scales, introduce a methodology for correction based on radial basis functions, apply it to TGSS-ADR1, and create a modified catalogue, TGSS-RSADR1. This catalogue comprises 383,589 TGSS-ADR1 sources with updated flux density and flux density uncertainty values, and covers $\mathrm{Declination}\leq+30^\circ$, $|b|\geq10^\circ$, a sky area of 18,800 deg$^2$.

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.

The S-band Polarisation All-Sky Survey (S-PASS) has observed the entire southern sky using the 64-metre Parkes radio telescope at 2.3GHz with an effective bandwidth of 184MHz. The surveyed sky area covers all declinations $\delta\leq 0^\circ$. To analyse compact sources the survey data have been re-processed to produce a set of 107 Stokes $I$ maps with 10.75arcmin resolution and the large scale emission contribution filtered out. In this paper we use these Stokes $I$ images to create a total intensity southern-sky extragalactic source catalogue at 2.3GHz. The source catalogue contains 23,389 sources and covers a sky area of 16,600deg$^2$, excluding the Galactic plane for latitudes $|b|<10^\circ$. Approximately 8% of catalogued sources are resolved. S-PASS source positions are typically accurate to within 35arcsec. At a flux density of 225mJy the S-PASS source catalogue is more than 95% complete, and $\sim$94% of S-PASS sources brighter than 500mJy beam$^{-1}$ have a counterpart at lower frequencies.

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

Long term radio monitoring of the broad absorption line quasar, Mrk\,231, at 17.6~GHz detected a strong flare in 2015. This triggered four epochs of Very Long Baseline Array (VLBA) observations from 8.4 GHz to 43 GHz as well as three epochs of X-ray observations with NuSTAR and two with XMM over a 15 week period. Two ejected components were detected by the VLBA observations. A conservative lower bound on the apparent speed of the first ejection is attained by assuming that it was ejected when the flare began, $v_{\rm{app}}>3.15$c. Serendipitous far UV Hubble Space Telescope observations combined with our long term radio monitoring seem to indicate that episodes of relativistic ejections suppress flux that is emitted at wavelengths shortward of the peak of the far UV spectral energy distribution, similar to what has been observed in radio loud quasars. Episodes of strong jet production also seem to suppress the high ionization broad absorption line wind seen in weak jet states. We found a statistically significant increase ($\sim 25\%$) of the 3-12 keV flux during the radio flare relative to a quiescent radio state. This is explained by an ultra-fast ($\sim 0.06$c) X-ray absorbing photo-ionized wind that is significantly detected only in the low radio state (similar to Galactic black holes). Mrk~231 is becoming more radio loud. We found that the putative parsec scale radio lobe doubled in brightness in 9 years. Furthermore, large flares are more frequent with 3 major flares occurring at $\sim 2$ year intervals.

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.

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.

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 154 MHz Murchison Widefield Array imaging observations and variability information for a sample of pulsars. Over the declination range $-80^{\circ} < {\delta} < 10^{\circ}$ we detect 17 known pulsars with mean flux density greater than 0.3 Jy. We explore the variability properties of this sample on timescales of minutes to years. For three of these pulsars, PSR J0953+0755, PSR J0437-4715 and PSR J0630-2834 we observe interstellar scintillation and variability on timescales of greater than 2 minutes. One further pulsar, PSR J0034-0721, showed significant variability, the physical origins of which are difficult to determine. The dynamic spectra for PSR J0953+0755 and PSR J0437-4715 show discrete time and frequency structure consistent with diffractive interstellar scintillation and we present the scintillation bandwidth and timescales from these observations. The remaining pulsars within our sample were statistically non-variable. We also explore the spectral properties of this sample and find spectral curvature in pulsars PSR J0835-4510, PSR J1752-2806 and PSR J0437-4715.

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.

Electron density irregularities in the ionosphere are known to be magnetically anisotropic, preferentially elongated along the lines of force. While many studies of their morphology have been undertaken by topside sounding and whistler measurements, it is only recently that detailed regional-scale reconstructions have become possible, enabled by the advent of widefield radio telescopes. Here we present a new approach for visualising and studying field-aligned irregularities (FAIs), which involves transforming interferometric measurements of TEC gradients onto a magnetic shell tangent plane. This removes the perspective distortion associated with the oblique viewing angle of the irregularities from the ground, facilitating the decomposition of dynamics along and across magnetic field lines. We apply this transformation to the dataset of Loi et al. [2015a], obtained on 15 October 2013 by the Murchison Widefield Array (MWA) radio telescope and displaying prominent FAIs. We study these FAIs in the new reference frame, quantifying field-aligned and field-transverse behaviour, examining time and altitude dependencies, and extending the analysis to FAIs on sub-array scales. We show that the inclination of the plane can be derived solely from the data, and verify that the best-fit value is consistent with the known magnetic inclination. The ability of the model to concentrate the fluctuations along a single spatial direction may find practical application to future calibration strategies for widefield interferometry, by providing a compact representation of FAI-induced distortions.

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.

We present a survey for transient and variable sources, on timescales from 28 seconds to $\sim$1 year, using the Murchison Widefield Array (MWA) at 182 MHz. Down to a detection threshold of 0.285 Jy, no transient candidates were identified, making this the most constraining low-frequency survey to date and placing a limit on the surface density of transients of $<4.1 \times 10^{-7}$ deg$^{-2}$ for the shortest timescale considered. At these frequencies, emission from Fast Radio Bursts (FRBs) is expected to be detectable in the shortest timescale images without any corrections for interstellar or intergalactic dispersion. At an FRB limiting flux density of 7980 Jy, we find a rate of $<$82 FRBs per sky per day for dispersion measures $<$700 pc cm$^{-3}$. Assuming a cosmological population of standard candles, our rate limits are consistent with the FRB rates obtained by Thornton et al. (2013) if they have a flat spectral slope. Finally, we conduct an initial variability survey of sources in the field with flux densities $\gtrsim$0.5 Jy and identify no sources with significant variability in their lightcurves. However, we note that substantial further work is required to fully characterise both the short term and low level variability within this field.

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.

Geomagnetically-aligned density structures with a range of sizes exist in the near-Earth plasma environment, including 10-100 km-wide VLF/HF wave-ducting structures. Their small diameters and modest density enhancements make them difficult to observe, and there is limited evidence for any of the several formation mechanisms proposed to date. We present a case study of an event on 26 August 2014 where a travelling ionospheric disturbance (TID) shortly precedes the formation of a complex collection of field-aligned ducts, using data obtained by the Murchison Widefield Array (MWA) radio telescope. Their spatiotemporal proximity leads us to suggest a causal interpretation. Geomagnetic conditions were quiet at the time, and no obvious triggers were noted. Growth of the structures proceeds rapidly, within 0.5 hr of the passage of the TID, attaining their peak prominence 1-2 hr later and persisting for several more hours until observations ended at local dawn. Analyses of the next two days show field-aligned structures to be preferentially detectable under quiet rather than active geomagnetic conditions. We used a raster scanning strategy facilitated by the speed of electronic beamforming to expand the quasi-instantaneous field of view of the MWA by a factor of three. These observations represent the broadest angular coverage of the ionosphere by a radio telescope to date.

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 report characterization results for an engineering prototype of a next-generation low-frequency radio astronomy array. This prototype, which we refer to as the Aperture Array Verification System 0.5 (AAVS0.5), is a sparse pseudo-random array of 16 log-periodic antennas designed for 70-450 MHz. It is co-located with the Murchison Widefield Array (MWA) at the Murchison Radioastronomy Observatory (MRO) near the Australian Square Kilometre Array (SKA) core site. We characterize the AAVS0.5 using two methods: in-situ radio interferometry with astronomical sources and an engineering approach based on detailed full-wave simulation. In-situ measurement of the small prototype array is challenging due to the dominance of the Galactic noise and the relatively weaker calibration sources easily accessible in the southern sky. The MWA, with its 128 "tiles" and up to 3 km baselines, enabled in-situ measurement via radio interferometry. We present array sensitivity and beam pattern characterization results and compare to detailed full-wave simulation. We discuss areas where differences between the two methods exist and offer possibilities for improvement. Our work demonstrates the value of the dual astronomy-simulation approach in upcoming SKA design work.

The Arcminute Microkelvin Imager Galactic Plane Survey (AMIGPS) provides mJy-sensitivity, arcminute-resolution interferometric images of the northern Galactic plane at $\approx$ 16 GHz. The first data release covered $76^{\circ} \lessapprox \ell \lessapprox 170^{\circ}$ between latitudes of $|b| \lessapprox 5^{\circ}$; here we present a second data release, extending the coverage to $53^{\circ} \lessapprox \ell \lessapprox 193^{\circ}$ and including high-latitude extensions to cover the Taurus and California giant molecular cloud regions, and the recently discovered large supernova remnant G159.6+7.3. The total coverage is now 1777 deg$^2$ and the catalogue contains 6509 sources. We also describe the improvements to the data processing pipeline which improves the positional and flux density accuracies of the survey.

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 observations of high-amplitude rapid (2 s) variability toward two bright, compact extragalactic radio sources out of several hundred of the brightest radio sources in one of the 30x30 deg MWA Epoch of Reionization fields using the Murchison Widefield Array (MWA) at 155 MHz. After rejecting intrinsic, instrumental, and ionospheric origins we consider the most likely explanation for this variability to be interplanetary scintillation (IPS), likely the result of a large coronal mass ejection propagating from the Sun. This is confirmed by roughly contemporaneous observations with the Ooty Radio Telescope. We see evidence for structure on spatial scales ranging from <1000 km to >1e6 km. The serendipitous night-time nature of these detections illustrates the new regime that the MWA has opened for IPS studies with sensitive night-time, wide-field, low-frequency observations. This regime complements traditional dedicated strategies for observing IPS and can be utilized in real-time to facilitate dedicated follow-up observations. At the same time, it allows large-scale surveys for compact (arcsec) structures in low-frequency radio sources despite the 2 arcmin resolution of the array.

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.

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.

We update the all-sky Planck catalogue of 1227 clusters and cluster candidates (PSZ1) published in March 2013, derived from Sunyaev-Zeldovich (SZ) effect detections using the first 15.5 months of Planck satellite observations. Addendum. We deliver an updated version of the PSZ1 catalogue, reporting the further confirmation of 86 Planck-discovered clusters. In total, the PSZ1 now contains 947 confirmed clusters, of which 214 were confirmed as newly discovered clusters through follow-up observations undertaken by the Planck Collaboration. The updated PSZ1 contains redshifts for 913 systems, of which 736 (~80.6%) are spectroscopic, and associated mass estimates derived from the Y_z mass proxy. We also provide a new SZ quality flag, derived from a novel artificial neural network classification of the SZ signal, for the remaining 280 candidates. Based on this assessment, the purity of the updated PSZ1 catalogue is estimated to be 94%. In this release, we provide the full updated catalogue and an additional readme file with further information on the Planck SZ detections.

We present an in-situ antenna characterization method and results for a "low-frequency" radio astronomy engineering prototype array, characterized over the 75-300 MHz frequency range. The presence of multiple cosmic radio sources, particularly the dominant Galactic noise, makes in-situ characterization at these frequencies challenging; however, it will be shown that high quality measurement is possible via radio interferometry techniques. This method is well-known in the radio astronomy community but seems less so in antenna measurement and wireless communications communities, although the measurement challenges involving multiple undesired sources in the antenna field-of-view bear some similarities. We discuss this approach and our results with the expectation that this principle may find greater application in related fields.

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

Recent observations with the Murchison Widefield Array at 185~MHz have serendipitously unveiled a heretofore unknown giant and relatively nearby ($z = 0.0178$) radio galaxy associated with NGC\,1534. The diffuse emission presented here is the first indication that NGC\,1534 is one of a rare class of objects (along with NGC\,5128 and NGC\,612) in which a galaxy with a prominent dust lane hosts radio emission on scales of $\sim$700\,kpc. We present details of the radio emission along with a detailed comparison with other radio galaxies with disks. NGC1534 is the lowest surface brightness radio galaxy known with an estimated scaled 1.4-GHz surface brightness of just 0.2\,mJy\,arcmin$^{-2}$. The radio lobes have one of the steepest spectral indices yet observed: $\alpha=-2.1\pm0.1$, and the core to lobe luminosity ratio is $<0.1$\%. We estimate the space density of this low brightness (dying) phase of radio galaxy evolution as $7\times10^{-7}$\,Mpc$^{-3}$ and argue that normal AGN cannot spend more than 6\% of their lifetime in this phase if they all go through the same cycle.

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

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.

We describe the all-sky Planck catalogue of clusters and cluster candidates derived from Sunyaev--Zeldovich (SZ) effect detections using the first 15.5 months of Planck satellite observations. The catalogue contains 1227 entries, making it over six times the size of the Planck Early SZ (ESZ) sample and the largest SZ-selected catalogue to date. It contains 861 confirmed clusters, of which 178 have been confirmed as clusters, mostly through follow-up observations, and a further 683 are previously-known clusters. The remaining 366 have the status of cluster candidates, and we divide them into three classes according to the quality of evidence that they are likely to be true clusters. The Planck SZ catalogue is the deepest all-sky cluster catalogue, with redshifts up to about one, and spans the broadest cluster mass range from (0.1 to 1.6) 10^{15}Msun. Confirmation of cluster candidates through comparison with existing surveys or cluster catalogues is extensively described, as is the statistical characterization of the catalogue in terms of completeness and statistical reliability. The outputs of the validation process are provided as additional information. This gives, in particular, an ensemble of 813 cluster redshifts, and for all these Planck clusters we also include a mass estimated from a newly-proposed SZ-mass proxy. A refined measure of the SZ Compton parameter for the clusters with X-ray counter-parts is provided, as is an X-ray flux for all the Planck clusters not previously detected in X-ray surveys.

We present a search for transient and variable radio sources at 154 MHz with the Murchison Widefield Array 32-tile prototype. Fifty-one images were obtained that cover a field of view of 1430 deg^2 centred on Hydra A. The observations were obtained over three days in 2010 March and three days in 2011 April and May. The mean cadence of the observations was 26 minutes and there was additional temporal information on day and year timescales. We explore the variability of a sample of 105 low frequency radio sources within the field. Four bright (S > 6 Jy) candidate variable radio sources were identified that displayed low levels of short timescale variability (26 minutes). We conclude that this variability is likely caused by simplifications in the calibration strategy or ionospheric effects. On the timescale of one year we find two sources that show significant variability. We attribute this variability to either refractive scintillation or intrinsic variability. No radio transients were identified and we place an upper limit on the surface density of sources rho < 7.5 x 10^-5 deg^-2 with flux densities > 5.5 Jy, and characteristic timescales of both 26 minutes and one year.

Radio monitoring of the broad absorption line quasar (BALQSO) Mrk 231 from 13.9 GHz to 17.6 GHz detected a strong flat spectrum flare. Even though BALQSOs are typically weak radio sources, the 17.6 GHz flux density doubled in ~150 days, from ~135 mJy to ~270 mJy. It is demonstrated that the elapsed rise time in the quasar rest frame and the relative magnitude of the flare is typical of some of the stronger flares in blazars that are associated with the ejection of discrete components on parsec scales. The decay of a similar flare was found in a previous monitoring campaign at 22 GHz. We conclude that these flares are not rare and indicate the likely ejection of a new radio component that can be resolved from the core with Very Long Baseline Interferometry. The implication is that Mrk 231 seems to be a quasar in which the physical mechanism that produces the BAL wind is in tension with the emergence of a fledgling blazar.

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 have obtained deep SZ observations towards 15 of the apparently hottest XMM Cluster Survey (XCS) clusters that can be observed with the Arcminute Microkelvin Imager (AMI). We use a Bayesian analysis to quantify the significance of our SZ detections. We detect the SZ effect at high significance towards three of the clusters and at lower significance for a further two clusters. Towards the remaining ten clusters, no clear SZ signal was measured. We derive cluster parameters using the XCS mass estimates as a prior in our Bayesian analysis. For all AMI-detected clusters, we calculate large-scale mass and temperature estimates while for all undetected clusters we determine upper limits on these parameters. We find that the large- scale mean temperatures derived from our AMI SZ measurements (and the upper limits from null detections) are substantially lower than the XCS-based core-temperature estimates. For clusters detected in the SZ, the mean temperature is, on average, a factor of 1.4 lower than temperatures from the XCS. For clusters undetected in SZ, the average 68% upper limit on the mean temperature is a factor of 1.9 below the XCS temperature.

We present CARMA observations of a massive galaxy cluster discovered in the AMI blind SZ survey. Without knowledge of the cluster redshift a Bayesian analysis of the AMI, CARMA and joint AMI & CARMA uv-data is used to quantify the detection significance and parameterise both the physical and observational properties of the cluster whilst accounting for the statistics of primary CMB anisotropies, receiver noise and radio sources. The joint analysis of the AMI & CARMA uv-data was performed with two parametric physical cluster models: the {\beta}-model; and the model described in Olamaie et al. 2012 with the pressure profile fixed according to Arnaud et al. 2010. The cluster mass derived from these different models is comparable but our Bayesian evidences indicate a preference for the {\beta}-profile which we, therefore, use throughout our analysis. From the CARMA data alone we obtain a Bayesian probability of detection ratio of 12.8:1 when assuming that a cluster exists within our search area; alternatively assuming that Jenkins et al. 2001 accurately predicts the number of clusters as a function of mass and redshift, the Bayesian probability of detection is 0.29:1. From the analysis of the AMI or AMI & CARMA data the probability of detection ratio exceeds 4.5x10^3:1. Performing a joint analysis of the AMI & CARMA data with a physical cluster model we derive the total mass internal to r200 as MT,200 = 4.1x10^14Msun. Using a phenomenological {\beta}-model to quantify the temperature decrement as a function of angular distance we find a central SZ temperature decrement of 170{\mu}K in the AMI & CARMA data. The SZ decrement in the CARMA data is weaker than expected and we speculate that this is a consequence of the cluster morphology. In a forthcoming study we will assess the impact of cluster morphology on the SZ decrements that are observed with interferometers such as AMI and CARMA.

We present 16-GHz Sunyaev-Zel'dovich observations using the Arcminute Microkelvin Imager (AMI) and subsequent Bayesian analysis of six galaxy clusters at redshift ($z \approx 1$) chosen from an X-ray and Infrared selected sample from Culverhouse et al. (2010). In the subsequent analysis we use two cluster models, an isothermal \beta-model and a Dark Matter GNFW (DM-GNFW) model in order to derive a formal detection probability and the cluster parameters. We detect two clusters (CLJ1415+3612 & XMJ0830+5241) and measure their total masses out to a radius of 200 $\times$ the critical density at the respective cluster's redshift. For CLJ1415+3612 and XMJ0830+5241, we find M_{\mathrm{T},200} for each model, which agree with each other for each cluster. We also present maps before and after source subtraction of the entire sample and provide 1D and 2D posterior marginalised probability distributions for each fitted cluster profile parameter of the detected clusters. Using simulations which take into account the measured source environment from the AMI Large Array (LA), source confusion noise, CMB primordials, instrument noise, we estimate from low-radius X-ray data from Culverhouse et al. (2010), the detectability of each cluster in the sample and compare it with the result from the Small Array (SA) data. Furthermore, we discuss the validity of the assumptions of isothermality and constant gas mass fraction. We comment on the bias that these small-radius estimates introduce to large-radius SZ predictions. In addition, we follow-up the two detections with deep, single-pointed LA observations. We find a 3 sigma tentative decrement toward CLJ1415+3612 at high-resolution and a 5 sigma high-resolution decrement towards XMJ0830+5241.

The AMI Galactic Plane Survey (AMIGPS) is a large area survey of the outer Galactic plane to provide arcminute resolution images at milli-Jansky sensitivity in the centimetre-wave band. Here we present the first data release of the survey, consisting of 868 deg^2 of the Galactic plane, covering the area 76 deg \lessapprox l \lessapprox 170 deg between latitudes of |b| \lessapprox 5 deg, at a central frequency of 15.75 GHz (1.9 cm). We describe in detail the drift scan observations which have been used to construct the maps, including the techniques used for observing, mapping and source extraction, and summarise the properties of the finalized datasets. These observations constitute the most sensitive Galactic plane survey of large extent at centimetre-wave frequencies greater than 1.4 GHz.

A subset of ultraluminous X-ray sources (those with luminosities < 10^40 erg/s) are thought to be powered by the accretion of gas onto black holes with masses of ~5-20 M_solar, probably via an accretion disc. The X-ray and radio emission are coupled in such Galactic sources, with the radio emission originating in a relativistic jet thought to be launched from the innermost regions near the black hole, with the most powerful emission occurring when the rate of infalling matter approaches a theoretical maximum (the Eddington limit). Only four such maximal sources are known in the Milky Way, and the absorption of soft X-rays in the interstellar medium precludes determining the causal sequence of events that leads to the ejection of the jet. Here we report radio and X-ray observations of a bright new X-ray source whose peak luminosity can exceed 10^39 erg/s in the nearby galaxy, M31. The radio luminosity is extremely high and shows variability on a timescale of tens of minutes, arguing that the source is highly compact and powered by accretion close to the Eddington limit onto a stellar mass black hole. Continued radio and X-ray monitoring of such sources should reveal the causal relationship between the accretion flow and the powerful jet emission.

The Planck 28.5 GHz maps were searched for potential Anomalous Microwave Emission (AME) regions on the scale of $\sim3^{\circ}$ or smaller, and several new regions of interest were selected. Ancillary data at both lower and higher frequencies were used to construct spectral energy distributions (SEDs), which seem to confirm an excess consistent with spinning dust models. Here we present higher resolution observations of two of these new regions with the Arcminute Microkelvin Imager Small Array (AMI SA) between 14 and 18 GHz to test for the presence of a compact ($\sim$10 arcmin or smaller) component. For AME-G107.1+5.2, dominated by the {\sc Hii} region S140, we find evidence for the characteristic rising spectrum associated with the either the spinning dust mechanism for AME or an ultra/hyper-compact \textsc{Hii} region across the AMI frequency band, however for AME-G173.6+2.8 we find no evidence for AME on scales of $\sim 2-10$ arcmin.

The Arcminute Microkelvin Imager (AMI) is a telescope specifically designed for high sensitivity measurements of low-surface-brightness features at cm-wavelength and has unique, important capabilities. It consists of two interferometer arrays operating over 13.5-18 GHz that image structures on scales of 0.5-10 arcmin with very low systematics. The Small Array (AMI-SA; ten 3.7-m antennas) couples very well to Sunyaev-Zel'dovich features from galaxy clusters and to many Galactic features. The Large Array (AMI-LA; eight 13-m antennas) has a collecting area ten times that of the AMI-SA and longer baselines, crucially allowing the removal of the effects of confusing radio point sources from regions of low surface-brightness, extended emission. Moreover AMI provides fast, deep object surveying and allows monitoring of large numbers of objects. In this White Paper we review the new science - both Galactic and extragalactic - already achieved with AMI and outline the prospects for much more.

We present 16-GHz AMI SZ observations of 19 clusters with L_X >7x10^37 W (h50=1) selected from the LoCuS survey (0.142<z<0.295) and of A1758b, in the FoV of A1758a. We detect 17 clusters with 5-23sigma peak surface brightnesses. Cluster parameters are obtained using a Bayesian cluster analysis. We fit isothermal beta-models to our data and assume the clusters are virialized (with all the kinetic energy in gas internal energy). Our gas temperature, T_AMI, is derived from AMI SZ data, not from X-ray spectroscopy. Cluster parameters internal to r500 are derived assuming HSE. We find: (i) Different gNFW parameterizations yield significantly different parameter degeneracies. (ii) For h70 = 1, we find the virial radius r200 to be typically 1.6+/-0.1 Mpc and the total mass M_T(r200) typically to be 2.0-2.5xM_T(r500).(iii) Where we have found M_T X-ray (X) and weak-lensing (WL) values in the literature, there is good agreement between WL and AMI estimates (with M_{T,AMI}/M_{T,WL} =1.2^{+0.2}_{-0.3} and =1.0+/-0.1 for r500 and r200, respectively). In comparison, most Suzaku/Chandra estimates are higher than for AMI (with M_{T,X}/M_{T,AMI}=1.7+/-0.2 within r500), particularly for the stronger mergers.(iv) Comparison of T_AMI to T_X sheds light on high X-ray masses: even at large r, T_X can substantially exceed T_AMI in mergers. The use of these higher T_X values will give higher X-ray masses. We stress that large-r T_SZ and T_X data are scarce and must be increased. (v) Despite the paucity of data, there is an indication of a relation between merger activity and SZ ellipticity. (vi) At small radius (but away from any cooling flow) the SZ signal (and T_AMI) is less sensitive to ICM disturbance than the X-ray signal (and T_X) and, even at high r, mergers affect n^2-weighted X-ray data more than n-weighted SZ, implying significant shocking or clumping or both occur even in the outer parts of mergers.

We present 16-GHz observations using the Arcminute Microkelvin Imager (AMI) of 11 clusters with 7 x 10^{37}W < L_X < 11 x 10^{37}W (h_{50}=1.0) selected from the Local Cluster Substructure Survey (LoCuSS) and compare them to X-ray data. We use a fast, Bayesian cluster analysis to explore the high-dimensional parameter space of the cluster-plus-sources model and obtain robust cluster parameter estimates in the presence of radio point sources, receiver noise and primordial CMB anisotropy. Our analysis fits a spherical, isothermal beta-model to our data and assumes the cluster follows the theoretical mass-temperature relation. Large-scale cluster parameters internal to r_{500} are derived under the assumption of hydrostatic equilibrium. Posterior distributions for the large-scale parameters of 8 of our clusters are given; SZ effects towards Abell 1704 and Zw0857.9+2107 were not detected and our spherical beta-profile was found to be an inadequate fit to the decrement on our map for Abell 2409.

A comparison is presented of Sunyaev-Zeldovich measurements for 11 galaxy clusters as obtained by Planck and by the ground-based interferometer, the Arcminute Microkelvin Imager. Assuming a universal spherically-symmetric Generalised Navarro, Frenk & White (GNFW) model for the cluster gas pressure profile, we jointly constrain the integrated Compton-Y parameter (Y_500) and the scale radius (theta_500) of each cluster. Our resulting constraints in the Y_500-theta_500 2D parameter space derived from the two instruments overlap significantly for eight of the clusters, although, overall, there is a tendency for AMI to find the Sunyaev-Zeldovich signal to be smaller in angular size and fainter than Planck. Significant discrepancies exist for the three remaining clusters in the sample, namely A1413, A1914, and the newly-discovered Planck cluster PLCKESZ G139.59+24.18. The robustness of the analysis of both the Planck and AMI data is demonstrated through the use of detailed simulations, which also discount confusion from residual point (radio) sources and from diffuse astrophysical foregrounds as possible explanations for the discrepancies found. For a subset of our cluster sample, we have investigated the dependence of our results on the assumed pressure profile by repeating the analysis adopting the best-fitting GNFW profile shape which best matches X-ray observations. Adopting the best-fitting profile shape from the X-ray data does not, in general, resolve the discrepancies found in this subset of five clusters. Though based on a small sample, our results suggest that the adopted GNFW model may not be sufficiently flexible to describe clusters universally.

We present an interesting Sunyaev-Zel'dovich (SZ) detection in the first of the Arcminute Microkelvin Imager (AMI) 'blind', degree-square fields to have been observed down to our target sensitivity of 100{\mu}Jy/beam. In follow-up deep pointed observations the SZ effect is detected with a maximum peak decrement greater than 8 \times the thermal noise. No corresponding emission is visible in the ROSAT all-sky X-ray survey and no cluster is evident in the Palomar all-sky optical survey. Compared with existing SZ images of distant clusters, the extent is large (\approx 10') and complex; our analysis favours a model containing two clusters rather than a single cluster. Our Bayesian analysis is currently limited to modelling each cluster with an ellipsoidal or spherical beta-model, which do not do justice to this decrement. Fitting an ellipsoid to the deeper candidate we find the following. (a) Assuming that the Evrard et al. (2002) approximation to Press & Schechter (1974) correctly gives the number density of clusters as a function of mass and redshift, then, in the search area, the formal Bayesian probability ratio of the AMI detection of this cluster is 7.9 \times 10^4:1; alternatively assuming Jenkins et al. (2001) as the true prior, the formal Bayesian probability ratio of detection is 2.1 \times 10^5:1. (b) The cluster mass is MT,200 = 5.5+1.2\times 10^14h-1M\odot. (c) Abandoning a physical model with num- -1.3 70 ber density prior and instead simply modelling the SZ decrement using a phenomenological {\beta}-model of temperature decrement as a function of angular distance, we find a central SZ temperature decrement of -295+36 {\mu}K - this allows for CMB primary anisotropies, receiver -15 noise and radio sources. We are unsure if the cluster system we observe is a merging system or two separate clusters.

We present 16 GHz (1.9 cm) deep radio continuum observations made with the Arcminute Microkelvin Imager (AMI) of a sample of low-mass young stars driving jets. We combine these new data with archival information from an extensive literature search to examine spectral energy distributions (SEDs) for each source and calculate both the radio and sub-mm spectral indices in two different scenarios: (1) fixing the dust temperature (Td) according to evolutionary class; (2) allowing Td to vary. We use the results of this analysis to place constraints on the physical mechanisms responsible for the radio emission. From AMI data alone, as well as from model fitting to the full SED in both scenarios, we find that 80 per cent of the objects in this sample have spectral indices consistent with free-free emission. We find an average spectral index in both Td scenarios consistent with free-free emission. We examine correlations of the radio luminosity with bolometric luminosity, envelope mass, and outflow force and find that these data are consistent with the strong correlation with envelope mass seen in lower luminosity samples. We examine the errors associated with determining the radio luminosity and find that the dominant source of error is the uncertainty on the opacity index, beta. We examine the SEDs for variability in these young objects, and find evidence for possible radio flare events in the histories of L1551 IRS 5 and Serpens SMM 1.

We present an analysis of observations made with the Arcminute Microkelvin Imager (AMI) and the Canada-France-Hawaii Telescope (CFHT) of six galaxy clusters in a redshift range of 0.16--0.41. The cluster gas is modelled using the Sunyaev--Zel'dovich (SZ) data provided by AMI, while the total mass is modelled using the lensing data from the CFHT. In this paper, we: i) find very good agreement between SZ measurements (assuming large-scale virialisation and a gas-fraction prior) and lensing measurements of the total cluster masses out to r_200; ii) perform the first multiple-component weak-lensing analysis of A115; iii) confirm the unusual separation between the gas and mass components in A1914; iv) jointly analyse the SZ and lensing data for the relaxed cluster A611, confirming our use of a simulation-derived mass-temperature relation for parameterizing measurements of the SZ effect.

Most Sunyaev--Zel'dovich (SZ) and X-ray analyses of galaxy clusters try to constrain the cluster total mass and/or gas mass using parameterised models and assumptions of spherical symmetry and hydrostatic equilibrium. By numerically exploring the probability distributions of the cluster parameters given the simulated interferometric SZ data in the context of Bayesian methods, and assuming a beta-model for the electron number density we investigate the capability of this model and analysis to return the simulated cluster input quantities via three rameterisations. In parameterisation I we assume that the T is an input parameter. We find that parameterisation I can hardly constrain the cluster parameters. We then investigate parameterisations II and III in which fg(r200) replaces temperature as a main variable. In parameterisation II we relate M_T(r200) and T assuming hydrostatic equilibrium. We find that parameterisation II can constrain the cluster physical parameters but the temperature estimate is biased low. In parameterisation III, the virial theorem replaces the hydrostatic equilibrium assumption. We find that parameterisation III results in unbiased estimates of the cluster properties. We generate a second simulated cluster using a generalised NFW (GNFW) pressure profile and analyse it with an entropy based model to take into account the temperature gradient in our analysis and improve the cluster gas density distribution. This model also constrains the cluster physical parameters and the results show a radial decline in the gas temperature as expected. The mean cluster total mass estimates are also within 1 sigma from the simulated cluster true values. However, we find that for at least interferometric SZ analysis in practice at the present time, there is no differences in the AMI visibilities between the two models. This may of course change as the instruments improve.

We present deep 1.8 cm (16 GHz) radio continuum imaging of seven young stellar objects in the Taurus molecular cloud. These objects have previously been extensively studied in the sub-mm to NIR range and their SEDs modelled to provide reliable physical and geometrical parametres.We use this new data to constrain the properties of the long-wavelength tail of the greybody spectrum, which is expected to be dominated by emission from large dust grains in the protostellar disk. We find spectra consistent with the opacity indices expected for such a population, with an average opacity index of beta = 0.26+/-0.22 indicating grain growth within the disks. We use spectra fitted jointly to radio and sub-mm data to separate the contributions from thermal dust and radio emission at 1.8 cm and derive disk masses directly from the cm-wave dust contribution. We find that disk masses derived from these flux densities under assumptions consistent with the literature are systematically higher than those calculated from sub-mm data, and meet the criteria for giant planet formation in a number of cases.

The Planck Early Release Compact Source Catalogue includes 26 sources with no obvious matches in other radio catalogues (of primarily extragalactic sources). Here we present observations made with the Arcminute Microkelvin Imager Small Array (AMI SA) at 15.75 GHz of the eight of the unmatched sources at declination > +10 degrees. Of the eight, four are detected and are associated with known objects. The other four are not detected with the AMI SA, and are thought to be spurious.

We present deep radio continuum observations of the cores identified as deeply embedded young stellar objects in the Serpens molecular cloud by the Spitzer c2d programme at a wavelength of 1.8cm with the Arcminute Microkelvin Imager Large Array (AMI-LA). These observations have a resolution of ~30arcsec and an average sensitivity of 19microJy/beam. The targets are predominantly Class I sources, and we find the detection rate for Class I objects in this sample to be low (18%) compared to that of Class 0 objects (67%), consistent with previous works. For detected objects we examine correlations of radio luminosity with bolometric luminosity and envelope mass and find that these data support correlations found by previous samples, but do not show any indiction of the evolutionary divide hinted at by similar data from the Perseus molecular cloud when comparing radio luminosity with envelope mass. We conclude that envelope mass provides a better indicator for radio luminosity than bolometric luminosity, based on the distribution of deviations from the two correlations. Combining these new data with archival 3.6cm flux densities we also examine the spectral indices of these objects and find an average spectral index of 0.53+/-1.14, consistent with the canonical value for a partially optically thick spherical or collimated stellar wind. However, we caution that possible inter-epoch variability limits the usefulness of this value, and such variability is supported by our identification of a possible flare in the radio history of Serpens SMM 1.

We present the first all-sky sample of galaxy clusters detected blindly by the Planck satellite through the Sunyaev-Zeldovich (SZ) effect from its six highest frequencies. This early SZ (ESZ) sample is comprised of 189 candidates, which have a high signal-to-noise ratio ranging from 6 to 29. Its high reliability (purity above 95%) is further ensured by an extensive validation process based on Planck internal quality assessments and by external cross-identification and follow-up observations. Planck provides the first measured SZ signal for about 80% of the 169 previously-known ESZ clusters. Planck furthermore releases 30 new cluster candidates, amongst which 20 meet the ESZ signal-to-noise selection criterion. At the submission date, twelve of the 20 ESZ candidates were confirmed as new clusters, with eleven confirmed using XMM-Newton snapshot observations, most of them with disturbed morphologies and low luminosities. The ESZ clusters are mostly at moderate redshifts (86% with z below 0.3) and span more than a decade in mass, up to the rarest and most massive clusters with masses above 10^15 Msol.

We present follow-up observations of two galaxy clusters detected blindly via the Sunyaev-Zel'dovich (SZ) effect and released in the Planck Early Release Compact Source Catalogue. We use the Arcminute Microkelvin Imager, a dual-array 14-18 GHz radio interferometer. After radio source subtraction, we find a SZ decrement of integrated flux density -1.08+/-0.10 mJy toward PLCKESZ G121.11+57.01, and improve the position measurement of the cluster, finding the centre to be RA 12 59 36.4, Dec +60 04 46.8, to an accuracy of 20 arcseconds. The region of PLCKESZ G115.71+17.52 contains strong extended emission, so we are unable to confirm the presence of this cluster via the SZ effect.

We present deep radio continuum observations of the cores identified as deeply embedded young stellar objects in the Perseus molecular cloud by the Spitzer c2d programme at a wavelength of 1.8 cm with the Arcminute Microkelvin Imager Large Array (AMI-LA). We detect 72% of Class 0 objects from this sample and 31% of Class I objects. No starless cores are detected. We use the flux densities measured from these data to improve constraints on the correlations between radio luminosity and bolometric luminosity, infrared luminosity and, where measured, outflow force. We discuss the differing behaviour of these objects as a function of protostellar class and investigate the differences in radio emission as a function of core mass. Two of four possible very low luminosity objects (VeLLOs) are detected at 1.8 cm.

We have observed an area of approximatley 27 deg^2 to an rms noise level of less than 0.2 mJy at 15.7 GHz, using the Arcminute Microkelvin Imager Large Array. These observations constitute the most sensitive radio-source survey of any extent (greater than approximately 0.2 deg^2) above 1.4 GHz. This paper presents the techniques employed for observing, mapping and source extraction. We have used a systematic procedure for extracting information and producing source catalogues, from maps with varying noise and uv-coverage. We have performed simulations to test our mapping and source-extraction procedures, and developed methods for identifying extended, overlapping and spurious sources in noisy images. In an accompanying paper, AMI Consortium: Davies et al. 2010, the first results from the 10C survey, including the deep 15.7-GHz source count, are presented.

The first results from the Tenth Cambridge (10C) Survey of Radio Sources, carried out using the AMI Large Array (LA) at an observing frequency of 15.7 GHz, are presented. The survey fields cover an area of approximately 27 sq. degrees to a flux-density completeness of 1 mJy. Results for some deeper areas, covering approximately 12 sq. degrees, wholly contained within the total areas and complete to 0.5 mJy, are also presented. The completeness for both areas is estimated to be at least 93 per cent. The source catalogue contains 1897 entries and is available at www.mrao.cam.ac.uk/surveys/10C. It has been combined with that of the 9C Survey to calculate the 15.7-GHz source counts. A broken power law is found to provide a good parameterisation of the differential count between 0.5 mJy and 1 Jy. The measured count has been compared to that predicted by de Zotti et al. (2005). The model displays good agreement with the data at the highest flux densities but under-predicts the integrated count between 0.5 mJy and 1 Jy by about 30 per cent. Entries from the source catalogue have been matched to those contained in the catalogues of NVSS and FIRST (both of which have observing frequencies of 1.4 GHz). This matching provides evidence for a shift in the typical 1.4-to-15.7-GHz spectral index of the 15.7-GHz-selected source population with decreasing flux density towards sub-mJy levels - the spectra tend to become less steep. Automated methods for detecting extended sources have been applied to the data; approximately 5 per cent of the sources are found to be extended relative to the LA synthesised beam of approximately 30 arcsec. Investigations using higher-resolution data showed that most of the genuinely extended sources at 16 GHz are classical doubles, although some nearby galaxies and twin-jet sources were also identified.

We present 13.9-18.2 GHz observations of the Sunyaev-Zel'dovich (SZ) effect towards Abell 2146 using the Arcminute Microkelvin Imager (AMI). The cluster is detected with a peak SNR ratio of 13 sigma in the radio source subtracted map. Comparison of the SZ and X-ray images suggests that they both have extended regions which lie approximately perpendicular to one another, with their emission peaks significantly displaced. These features indicate non-uniformities in the distributions of the gas temperature and pressure, indicative of a cluster merger. We use a Bayesian cluster analysis to explore the high-dimensional parameter space of the cluster-plus-sources model to obtain cluster parameter estimates in the presence of radio point sources, receiver noise and primordial CMB anisotropy; the probability of SZ + CMB primordial structure + radio sources + receiver noise to CMB + radio sources + receiver noise is 3 x 10^{6}:1. We compare the results from three different cluster models. Our preferred model exploits the observation that the gas fractions do not appear to vary greatly between clusters. Given the relative masses of the two merging systems in Abell 2146, the mean gas temperature can be deduced from the virial theorem (assuming all of the kinetic energy is in the form of internal gas energy) without being affected significantly by the merger event, provided the primary cluster was virialized before the merger. In this model we fit a simple spherical isothermal beta-model, despite the inadequacy of this model for a merging system like Abell 2146, and assume the cluster follows the mass-temperature relation of a virialized, singular, isothermal sphere. We note that this model avoids inferring large-scale cluster parameters internal to r_200 under the widely used assumption of hydrostatic equilibrium. We find that at r_200 M_T= 4.1 \pm 0.5 x 10^{14} h^{-1}M_sun and T=4.5 \pm 0.5 keV.

We present observations from the Small Array of the Arcminute Microkelvin Imager (AMI) of eight high X-ray luminosity galaxy cluster systems selected from the Local Cluster Substructure Survey (LoCuSS) sample.We detect the Sunyaev-Zel'dovich (SZ) effect in seven of these clusters. With the assumptions that galaxy clusters are isothermal, have a density profile described by a spherical b -model and obey the theoretical M-T relation, we are able to derive cluster parameters at r200 from our SZ data. With the additional assumption of hydrostatic equilibrium we are able to derive parameters at r500. We present posterior probability distributions for cluster parameters such as mass, radius and temperature (TSZ, MT). Combining our sample with that of AMI Consortium: Rodr'\iguez-Gonz\'alvez et al. (2011) and using large-radius X-ray temperature estimates (TX) from Chandra and Suzaku observations, we find that there is reasonable correspondence between TX and TSZ,MT values at low TX, but that for clusters with TX above around 6keV the correspondence breaks down with TX exceeding TSZ, MT; we stress that this finding is based on just ten clusters.

We perform deep 1.8 cm radio continuum imaging towards thirteen protostellar regions selected from the Spitzer c2d small clouds and cores programme at high resolution (25") in order to detect and quantify the cm-wave emission from deeply embedded young protostars. Within these regions we detect fifteen compact radio sources which we identify as radio protostars including two probable new detections. The sample is in general of low bolometric luminosity and contains several of the newly detected VeLLO sources. We determine the 1.8 cm radio luminosity to bolometric luminosity correlation, L_rad -L_bol, for the sample and discuss the nature of the radio emission in terms of the available sources of ionized gas. We also investigate the L_rad-L_IR correlation and suggest that radio flux density may be used as a proxy for the internal luminosity of low luminosity protostars.

We present observations using the Small Array of the Arcminute Microkelvin Imager (AMI; 14-18 GHz) of four Abell and three MACS clusters spanning 0.171-0.686 in redshift. We detect Sunyaev-Zel'dovich (SZ) signals in five of these without any attempt at source subtraction, although strong source contamination is present. With radio-source measurements from high-resolution observations, and under the assumptions of spherical $\beta$-model, isothermality and hydrostatic equilibrium, a Bayesian analysis of the data in the visibility plane detects extended SZ decrements in all seven clusters over and above receiver noise, radio sources and primary CMB imprints. Bayesian evidence ratios range from 10^{11}:1 to 10^{43}:1 for six of the clusters and 3000:1 for one with substantially less data than the others. We present posterior probability distributions for, e.g., total mass and gas fraction averaged over radii internal to which the mean overdensity is 1000, 500 and 200, r_200 being the virial radius. Reaching r_200 involves some extrapolation for the nearer clusters but not for the more-distant ones. We find that our estimates of gas fraction are low (compared with most in the literature) and decrease with increasing radius. These results appear to be consistent with the notion that gas temperature in fact falls with distance (away from near the cluster centre) out to the virial radius.

We report new cm-wave measurements at five frequencies between 15 and 18GHz of the continuum emission from the reportedly anomalous "region 4" of the nearby galaxy NGC6946. We find that the emission in this frequency range is significantly in excess of that measured at 8.5GHz, but has a spectrum from 15-18GHz consistent with optically thin free-free emission from a compact HII region. In combination with previously published data we fit four emission models containing different continuum components using the Bayesian spectrum analysis package radiospec. These fits show that, in combination with data at other frequencies, a model with a spinning dust component is slightly preferred to those that possess better-established emission mechanisms.

We present 25 arcsecond resolution radio images of five Lynds Dark Nebulae (L675, L944, L1103, L1111 & L1246) at 16 GHz made with the Arcminute Microkelvin Imager (AMI) Large Array. These objects were previously observed with the AMI Small Array to have an excess of emission at microwave frequencies relative to lower frequency radio data. In L675 we find a flat spectrum compact radio counterpart to the 850 micron emission seen with SCUBA and suggest that it is cm-wave emission from a previously unknown deeply embedded young protostar. In the case of L1246 the cm-wave emission is spatially correlated with 8 micron emission seen with Spitzer. Since the MIR emission is present only in Spitzer band 4 we suggest that it arises from a population of PAH molecules, which also give rise to the cm-wave emission through spinning dust emission.

Observations at 14.2 to 17.9 GHz made with the AMI Small Array towards fourteen Lynds Dark Nebulae with a resolution of 2' are reported. These sources are selected from the SCUBA observations of Visser et al. (2001) as small angular diameter clouds well matched to the synthesized beam of the AMI Small Array. Comparison of the AMI observations with radio observations at lower frequencies with matched uv-plane coverage is made, in order to search for any anomalous excess emission which can be attributed to spinning dust. Possible emission from spinning dust is identified as a source within a 2' radius of the Scuba position of the Lynds dark nebula, exhibiting an excess with respect to lower frequency radio emission. We find five sources which show a possible spinning dust component in their spectra. These sources have rising spectral indices in the frequency range 14.2--17.9 GHz. Of these five one has already been reported, L1111, we report one new definite detection, L675, and three new probable detections (L944, L1103 and L1246). The relative certainty of these detections is assessed on the basis of three criteria: the extent of the emission, the coincidence of the emission with the Scuba position and the likelihood of alternative explanations for the excess. Extended microwave emission makes the likelihood of the anomalous emission arising as a consequence of a radio counterpart to a protostar or a proto-planetary disk unlikely. We use a 2' radius in order to be consistent with the IRAS identifications of dark nebulae (Parker 1988), and our third criterion is used in the case of L1103 where a high flux density at 850 microns relative to the FIR data suggests a more complicated emission spectrum.

We present follow-up observations of 97 point sources from the Wilkinson Microwave Anisotropy Probe (WMAP) 3-year data, contained within the New Extragalactic WMAP Point Source (NEWPS) catalogue between declinations of -4 and +60 degrees; the sources form a flux-density-limited sample complete to 1.1 Jy (approximately 5 sigma) at 33 GHz. Our observations were made at 16 GHz using the Arcminute Microkelvin Imager (AMI) and at 33 GHz with the Very Small Array (VSA). 94 of the sources have reliable, simultaneous -- typically a few minutes apart -- observations with both telescopes. The spectra between 13.9 and 33.75 GHz are very different from those of bright sources at low frequency: 44 per cent have rising spectra (alpha < 0.0), where flux density is proportional to frequency^-alpha, and 93 per cent have spectra with alpha < 0.5; the median spectral index is 0.04. For the brighter sources, the agreement between VSA and WMAP 33-GHz flux densities averaged over sources is very good. However, for the fainter sources, the VSA tends to measure lower values for the flux densities than WMAP. We suggest that the main cause of this effect is Eddington bias arising from variability.

Using the Arcminute Microkelvin Imager (AMI) at 16 GHz and the Very Small Array (VSA) at 33 GHz to make follow-up observations of sources in the New Extragalactic WMAP Point Source Catalogue, we have investigated the flux-density variability in a complete sample of 97 sources over timescales of a few months to approximately 1.5 years. We find that 53 per cent of the 93 sources, for which we have multiple observations, are variable, at the 99 per cent confidence level, above the flux density calibration uncertainties of approximately 4 per cent at 16 GHz; the fraction of sources having varied by more than 20 per cent is 15 per cent at 16 GHz and 20 per cent at 33 GHz. Not only is this common occurrence of variability at high frequency of interest for source physics, but strategies for coping with source contamination in CMB work must take this variability into account. There is no strong evidence of a correlation between variability and flux density for the sample as a whole. Using a maximum-likelihood method, we calculate the correlation in the variability at the two frequencies in a subset of sources classified as variable from both the AMI and VSA data and find the Pearson product-moment correlation coefficient to be very high (0.955 +/- 0.034). We also find the degree of variability at 16 GHz (0.202 +/- 0.028) to be very similar to that at 33 GHz (0.224 +/- 0.039). Finally, we have investigated the relationship between variability and spectral index, alpha_13.9^33.75 (where S proportional to nu^-alpha), and find a significant difference in the spectral indices of the variable sources (-0.06 +/- 0.05) and non-variable sources (0.13 +/- 0.04).

We present observations between 1.4 and 18 GHz confirming that G64.5+0.9 is new Galactic shell supernova remnant, using the Very Large Array and the Arcminute Microkelvin Imager. The remnant is a shell ~8 arcmin in diameter with a spectral index of alpha = 0.47 +/- 0.03 (with alpha defined such that flux density S varies with frequency nu as S proportional to nu to the power of -alpha). There is also emission near the centre of the shell, ~1 arcmin in extent, with a spectral index of alpha = 0.81 +/- 0.02. We do not find any evidence for spectral breaks for either source within our frequency range. The nature of the central object is unclear and requires further investigation, but we argue that is most unlikely to be extragalactic. It is difficult to avoid the conclusion that it is associated with the shell, although its spectrum is very unlike that of known pulsar wind nebulae.

We present observations between 14.2 and 17.9 GHz of 12 reported supernova remnants (SNRs) made with the Arcminute Microkelvin Imager Small Array (AMI SA). In conjunction with data from the literature at lower radio frequencies, we determine spectra of these objects. For well-studied SNRs (Cas A, Tycho's SNR, 3C58 and the Crab Nebula), the results are in good agreement with spectra based on previous results. For the less well-studied remnants the AMI SA observations provide higher-frequency radio observations than previously available, and better constrain their radio spectra. The AMI SA results confirm a spectral turnover at ~11 GHz for the filled-centre remnant G74.9+1.2. We also see a possible steepening of the spectrum of the filled-centre remnant G54.1+0.3 within the AMI SA frequency band compared with lower frequencies. We confirm that G84.9+0.5, which had previously been identified as a SNR, is rather an HII region and has a flat radio spectrum.

We present observations of the Lynds' dark nebula LDN 1111 made at microwave frequencies between 14.6 and 17.2 GHz with the Arcminute Microkelvin Imager (AMI). We find emission in this frequency band in excess of a thermal free--free spectrum extrapolated from data at 1.4 GHz with matched uv-coverage. This excess is > 15 sigma above the predicted emission. We fit the measured spectrum using the spinning dust model of Drain & Lazarian (1998a) and find the best fitting model parameters agree well with those derived from Scuba data for this object by Visser et al. (2001).

The Arcminute Microkelvin Imager is a pair of interferometer arrays operating with six frequency channels spanning 13.9-18.2 GHz, with very high sensitivity to angular scales 30''-10'. The telescope is aimed principally at Sunyaev-Zel'dovich imaging of clusters of galaxies. We discuss the design of the telescope and describe and explain its electronic and mechanical systems.

We present observations between 14.2 and 17.9 GHz of sixteen Galactic HII regions made with the Arcminute Microkelvin Imager (AMI). In conjunction with data from the literature at lower radio frequencies we investigate the possibility of a spinning dust component in the spectra of these objects. We conclude that there is no significant evidence for spinning dust towards these sources and measure an average spectral index of 0.15+/-0.07 between 1.4 and 17.9 GHz for the sample.