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

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.

Infrared imaging of the colliding-wind binary Apep has revealed a spectacular dust plume with complicated internal dynamics that challenges standard colliding-wind binary physics. Such challenges can be potentially resolved if a rapidly-rotating Wolf-Rayet star is located at the heart of the system, implicating Apep as a Galactic progenitor system to long-duration gamma-ray bursts. One of the difficulties in interpreting the dynamics of Apep is that the spectral composition of the stars in the system was unclear. Here we present visual to near-infrared spectra that demonstrate that the central component of Apep is composed of two classical Wolf-Rayet stars of carbon- (WC8) and nitrogen-sequence (WN4-6b) subtypes. We argue that such an assignment represents the strongest case of a classical WR+WR binary system in the Milky Way. The terminal line-of-sight wind velocities of the WC8 and WN4-6b stars are measured to be $2100 \pm 200$ and $3500 \pm 100$ km s$^{-1}$, respectively. If the mass-loss rate of the two stars are typical for their spectral class, the momentum ratio of the colliding winds is expected to be $\approx$ 0.4. Since the expansion velocity of the dust plume is significantly smaller than either of the measured terminal velocities, we explore the suggestion that one of the Wolf-Rayet winds is anisotropic. We can recover a shock-compressed wind velocity consistent with the observed dust expansion velocity if the WC8 star produces a significantly slow equatorial wind with a velocity of $\approx$530 km s$^{-1}$. Such slow wind speeds can be driven by near-critical rotation of a Wolf-Rayet star.

We present LOFAR observations at 150 MHz of the borderline FRI/FRII giant radio galaxy NGC 6251. This paper presents the most sensitive and highest-resolution images of NGC 6251 at these frequencies to date, revealing for the first time a low-surface-brightness extension to the northern lobe, and a possible backflow associated with the southern lobe. The integrated spectra of components of NGC 6251 are consistent with previous measurements at higher frequencies, similar to results from other LOFAR studies of nearby radio galaxies. We find the outer structures of NGC 6251 to be either at equipartition or slightly electron dominated, similar to those of FRII sources rather than FRIs; but this conclusion remains tentative because of uncertainties associated with the geometry and the extrapolation of X-ray measurements to determine the external pressure distribution on the scale of the outer lobes. We place lower limits on the ages of the extension of the northern lobe and the backflow of the southern lobe of $t \gtrsim 250$ Myr and $t \gtrsim 210$ Myr respectively. We present the first detection of polarisation at 150 MHz in NGC 6251. Taking advantage of the high Faraday resolution of LOFAR, we place an upper limit on the magnetic field in the group of $B < 0.2 (\Lambda_B / 10 {\rm kpc})^{-0.5} \mu$G for a coherence scale of $\Lambda_B < 60 {\rm kpc}$ and $B < 13 \mu$G for $\Lambda_B = 240$ kpc.

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.

Low frequency ($\nu\lesssim 150\,$MHz) stellar radio emission is expected to originate in the outer corona at heights comparable to and larger than the stellar radius. Such emission from the Sun has been used to study coronal structure, mass ejections, space-weather conditions around the planets (Schwenn 2006). Searches for low-frequency emission from other stars have only detected a single active flare-star (Lynch et al. 2017) that is not representative of the wider stellar population. Here we report the detection of low-frequency radio emission from a quiescent star, GJ 1151--- a member of the most common stellar type (red dwarf or spectral class M) in the Galaxy. The characteristics of the emission are similar to those of planetary auroral emissions (Zarka 1998) (e.g. Jupiter's decametric emission), suggesting a coronal structure dominated by a global magnetosphere with low plasma density. Our results show that large-scale currents that power radio aurorae operate over a vast range of mass and atmospheric composition, ranging from terrestrial planets to main-sequence stars. The Poynting flux required to produce the observed radio emission cannot be generated by GJ 1151's slow rotation, but can originate in a sub-Alfv\'{e}nic interaction of its magnetospheric plasma with a short-period exoplanet. The emission properties are consistent with theoretical expectations (Zarka 2007; Lanza 2009; Saur et al. 2013; Turnpenney et al. 2018) for interaction with an Earth-size planet in a $\sim 1-5$ day-long orbit.

The recent detection of circularly polarized, long-duration (>8 hr) low-frequency (~150 MHz) radio emission from the M4.5 dwarf GJ 1151 has been interpreted as arising from a star-planet interaction via the electron cyclotron maser instability. The existence or parameters of the proposed planets have not been determined. Using 20 new HARPS-N observations, we put 99th-percentile upper limits on the mass of any close companion to GJ 1151 at Msini < 5.6 M earth. With no stellar, brown dwarf, or giant planet companion likely in a close orbit, our data are consistent with detected radio emission emerging from a magnetic interaction between a short-period terrestrial-mass planet and GJ 1151.

Serendipitously, we have rediscovered radio emission from the binary system Ross 867 (M4.5V) and Ross 868 (M3.5V) while inspecting archival Giant Metrewave Radio Telescope (GMRT) observations. The binary system consists of two M-dwarf stars that share common characteristics such as spectral type, astrometric parameters, age and emission at infrared, optical and X-rays frequencies. The GMRT data at 610 MHz taken on July 2011 shows that the radio emission from Ross 867 is polarized and highly variable on hour time scales with a peak flux of 10.4 $\pm$ 0.7 mJy/beam. Additionally, after reviewing archival data from several observatories (VLA, GMRT, JVLA and LOFAR), we confirm that although both stars are likely coeval, only Ross 867 has been detected, while Ross 868 remains undetected at radio wavelengths. As they have a a large orbital separation, this binary stellar system provides a coeval laboratory to examine and constrain the stellar properties linked to radio activity in M dwarfs. We speculate that the observed difference in radio activity between the dwarfs could be due to vastly different magnetic field topologies or that Ross 867 has an intrinsically different dynamo.

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

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.

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.

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.

Obscuration of quasars by accreted gas and dust, or dusty intervening galaxies, can cause active galactic nuclei (AGN) to be missed in optically-selected surveys. Radio observations can overcome this dust bias. In particular, radio surveys searching for HI absorption inform us on how the AGN can impact on the cold neutral gas medium within the host galaxy, or the population of intervening galaxies through the observed line of sight gas kinematics. We present the results of a HI absorption line survey at $0.4 < z < 1$ towards 34 obscured quasars with the Australian SKA Pathfinder (ASKAP) commissioning array. We detect three HI absorption lines, with one of these systems previously unknown. Through optical follow-up for two sources, we find that in all detections the HI gas is associated with the AGN, and hence that these AGN are obscured by material within their host galaxies. Most of our sample are compact, and in addition, are either gigahertz peaked spectrum (GPS), or steep spectrum (CSS) sources, both thought to represent young or recently re-triggered radio AGN. The radio spectral energy distribution classifications for our sample agree with galaxy evolution models in which the obscured AGN has only recently become active. Our associated HI detection rate for GPS and compact SS sources matches those of other surveys towards such sources. We also find shallow and asymmetric HI absorption features, which agrees with previous findings that the cold neutral medium in compact radio galaxies is typically kinematically disturbed by the AGN.

We have examined the giant radio galaxy 3C~236 using LOFAR at 143 MHz down to an angular resolution of 7", in combination with observations at higher frequencies. We have used the low frequency data to derive spectral index maps with the highest resolution yet at these low frequencies. We confirm a previous detection of an inner hotspot in the north-west lobe and for the first time observe that the south-east lobe hotspot is in fact a triple hotspot, which may point to an intermittent source activity. Also, the spectral index map of 3C 236 shows that the spectral steepening at the inner region of the northern lobe is prominent at low frequencies. The outer regions of both lobes show spectral flattening, in contrast with previous high frequency studies. We derive spectral age estimates for the lobes, as well as particle densities of the IGM at various locations. We propose that the morphological differences between the lobes are driven by variations in the ambient medium density as well as the source activity history.

Extragalactic peaked-spectrum radio sources are thought to be the progenitors of larger, radio-loud active galactic nuclei (AGN). Synchrotron self-absorption (SSA) has often been identified as the cause of their spectral peak. The identification of new megahertz-peaked spectrum sources from the GaLactic and Extragalactic All-sky Murchison Widefield Array (GLEAM) survey provides an opportunity to test how radio sources with spectral peaks below 1 GHz fit within this evolutionary picture. We observed six peaked-spectrum sources selected from the GLEAM survey, three that have spectral characteristics which violate SSA and three that have spectral peaks below 230 MHz, with the Very Long Baseline Array at 1.55 and 4.96 GHz. We present milliarcsecond resolution images of each source and constrain their morphology, linear size, luminosity, and magnetic field strength. Of the sources that are resolved by our study, the sources that violate SSA appear to be compact doubles, while the sources with peak frequencies below 230 MHz have core-jet features. We find that all of our sources are smaller than expected from SSA by factors of >20. We also find that component magnetic field strengths calculated from SSA are likely inaccurate, differing by factors of >5 from equipartition estimates. The calculated equipartition magnetic field strengths more closely resemble estimates from previously studied gigahertz-peaked spectrum sources. Exploring a model of the interaction between jets and the interstellar medium, we demonstrate that free-free absorption (FFA) can accurately describe the linear sizes and peak frequencies of our sources. Our findings support the theory that there is a fraction of peaked-spectrum sources whose spectral peaks are best modelled by FFA, implying our understanding of the early stages of radio AGN is incomplete.

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

The LOFAR Two-metre Sky Survey (LoTSS) is a sensitive, high-resolution 120-168 MHz survey of the Northern sky. The LoTSS First Data Release (DR1) presents 424 square degrees of radio continuum observations over the HETDEX Spring Field (10h45m00s $<$ right ascension $<$ 15h30m00s and 45$^\circ$00$'$00$'$ $<$ declination $<$ 57$^\circ$00$'$00$''$) with a median sensitivity of 71$\mu$Jy/beam and a resolution of 6$''$. In this paper we present photometric redshifts (photo-$z$) for 94.4% of optical sources over this region that are detected in the Panoramic Survey Telescope and Rapid Response System (Pan-STARRS) 3$\pi$ steradian survey. Combining the Pan-STARRS optical data with mid-infrared photometry from the Wide-field Infrared Survey Explorer, we estimate photo-$z$s using a novel hybrid photometric redshift methodology optimised to produce the best possible performance for the diverse sample of radio continuum selected sources. For the radio-continuum detected population, we find an overall scatter in the photo-$z$ of 3.9% and an outlier fraction ($\left | z_{\rm{phot}} - z_{\rm{spec}} \right | / (1+z_{\rm{spec}}) > 0.15$) of 7.9%. We also find that, at a given redshift, there is no strong trend in photo-$z$ quality as a function of radio luminosity. However there are strong trends as a function of redshift for a given radio luminosity, a result of selection effects in the spectroscopic sample and/or intrinsic evolution within the radio source population. Additionally, for the sample of sources in the LoTSS First Data Release with optical counterparts, we present rest-frame optical and mid-infrared magnitudes based on template fits to the consensus photometric (or spectroscopic when available) redshift.

Context: New generation low-frequency telescopes are exploring a new parameter space in terms of depth and resolution. The data taken with these interferometers, for example with the LOw Frequency ARray (LOFAR), are often calibrated in a low signal-to-noise ratio regime and the removal of critical systematic effects is challenging. The process requires an understanding of their origin and properties. Aim: In this paper we describe the major systematic effects inherent to next generation low-frequency telescopes, such as LOFAR. With this knowledge, we introduce a data processing pipeline that is able to isolate and correct these systematic effects. The pipeline will be used to calibrate calibrator observations as the first step of a full data reduction process. Methods: We processed two LOFAR observations of the calibrator 3C196: the first using the Low Band Antenna (LBA) system at 42-66 MHz and the second using the High Band Antenna (HBA) system at 115-189 MHz. Results: We were able to isolate and correct for the effects of clock drift, polarisation misalignment, ionospheric delay, Faraday rotation, ionospheric scintillation, beam shape, and bandpass. The designed calibration strategy produced the deepest image to date at 54 MHz. The image has been used to confirm that the spectral energy distribution of the average radio source population tends to flatten at low frequencies. Conclusions: We prove that LOFAR systematic effects can be described by a relatively small number of parameters. Furthermore, the identification of these parameters is fundamental to reducing the degrees of freedom when the calibration is carried out on fields that are not dominated by a strong calibrator.

We constructed a sample of 23,344 radio-loud active galactic nuclei (RLAGN) from the catalogue derived from the LOFAR Two-Metre Sky Survey (LoTSS) survey of the HETDEX Spring field. Although separating AGN from star-forming galaxies remains challenging, the combination of spectroscopic and photometric techniques we used gives us one of the largest available samples of candidate RLAGN. We used the sample, combined with recently developed analytical models, to investigate the lifetime distribution of RLAGN. We show that large or giant powerful RLAGN are probably the old tail of the general RLAGN population, but that the low-luminosity RLAGN candidates in our sample, many of which have sizes $<100$ kpc, either require a very different lifetime distribution or have different jet physics from the more powerful objects. We then used analytical models to develop a method of estimating jet kinetic powers for our candidate objects and constructed a jet kinetic luminosity function based on these estimates. These values can be compared to observational quantities, such as the integrated radiative luminosity of groups and clusters, and to the predictions from models of RLAGN feedback in galaxy formation and evolution. In particular, we show that RLAGN in the local Universe are able to supply all the energy required per comoving unit volume to counterbalance X-ray radiative losses from groups and clusters and thus prevent the hot gas from cooling. Our computation of the kinetic luminosity density of local RLAGN is in good agreement with other recent observational estimates and with models of galaxy formation.

Cosmological simulations predict that an intergalactic magnetic field (IGMF) pervades the large scale structure (LSS) of the Universe. Measuring the IGMF is important to determine its origin (i.e. primordial or otherwise). Using data from the LOFAR Two Metre Sky Survey (LoTSS), we present the Faraday rotation measure (RM) and depolarisation properties of the giant radio galaxy J1235+5317, at a redshift of $z = 0.34$ and 3.38 Mpc in size. We find a mean RM difference between the lobes of $2.5\pm0.1$ rad/m$^2$ , in addition to small scale RM variations of ~0.1 rad/m$^2$ . From a catalogue of LSS filaments based on optical spectroscopic observations in the local universe, we find an excess of filaments intersecting the line of sight to only one of the lobes. Associating the entire RM difference to these LSS filaments leads to a gas density-weighted IGMF strength of ~0.3 {\mu}G. However, direct comparison with cosmological simulations of the RM contribution from LSS filaments gives a low probability (~5%) for an RM contribution as large as 2.5 rad/m$^2$ , for the case of IGMF strengths of 10 to 50 nG. It is likely that variations in the RM from the Milky Way (on 11' scales) contribute significantly to the mean RM difference, and a denser RM grid is required to better constrain this contribution. In general, this work demonstrates the potential of the LOFAR telescope to probe the weak signature of the IGMF. Future studies, with thousands of sources with high accuracy RMs from LoTSS, will enable more stringent constraints on the nature of the IGMF.

Determining the star-forming properties of radio-quiet quasars is important for understanding the co-evolution of star formation and black hole accretion. Here, we present the detection of the gravitationally-lensed radio-quiet quasars SDSS J1055+4628, SDSS J1313+5151 and SBS 1520+530 at 144 MHz that fall in the HETDEX Spring Field targeted in the LOFAR Two-metre Sky Survey (LoTSS) first full data release. We compare their radio and far-infrared luminosities relative to the radio-infrared correlation and find that their radio luminosities can be explained by star formation. The implied star formation rates derived from their radio and infrared luminosities are between 20 and 300 $\rm{M_{\odot}~yr^{-1}}$. These detections represent the first study of gravitationally lensed sources with LOFAR, opening a new frequency window for investigating the star-forming properties of high-redshift quasars at radio wavelengths. We consider the implications for future data releases and estimate that many of the objects in our parent sample will be detected during LoTSS, significantly increasing the fraction of gravitationally lensed radio-quiet quasars with radio detections.

The LOFAR Two-metre Sky Survey (LoTSS) is an ongoing sensitive, high-resolution 120-168 MHz survey of the Northern sky with diverse and ambitious science goals. Many of the scientific objectives of LoTSS rely upon, or are enhanced by, the association or separation of the sometimes incorrectly catalogued radio components into distinct radio sources, and the identification and characterisation of the optical counterparts to these sources. Here we present the source associations and optical and/or IR identifications for sources in the first data release, which are made using a combination of statistical techniques and visual association and identification. We document in detail the colour- and magnitude-dependent likelihood ratio method used for statistical identification as well as the Zooniverse project, called LOFAR Galaxy Zoo, used for the visual classification. We describe the process used to select which of these two different methods is most appropriate for each LoTSS source. The final LoTSS-DR1-IDs value-added catalogue presented contains 318,520 radio sources, of which 231,716 (73%) have optical and/or IR identifications in Pan-STARRS and WISE. The value-added catalogue is available online at https://lofar-surveys.org/, as part of this data release.

The LOFAR Two-metre Sky Survey (LoTSS) is an ongoing sensitive, high-resolution 120-168MHz survey of the entire northern sky for which observations are now 20% complete. We present our first full-quality public data release. For this data release 424 square degrees, or 2% of the eventual coverage, in the region of the HETDEX Spring Field (right ascension 10h45m00s to 15h30m00s and declination 45$^\circ$00$'$00$''$ to 57$^\circ$00$'$00$''$) were mapped using a fully automated direction-dependent calibration and imaging pipeline that we developed. A total of 325,694 sources are detected with a signal of at least five times the noise, and the source density is a factor of $\sim 10$ higher than the most sensitive existing very wide-area radio-continuum surveys. The median sensitivity is S$_{\rm 144 MHz} = 71\,\mu$Jy beam$^{-1}$ and the point-source completeness is 90% at an integrated flux density of 0.45mJy. The resolution of the images is 6$''$ and the positional accuracy is within 0.2$''$. This data release consists of a catalogue containing location, flux, and shape estimates together with 58 mosaic images that cover the catalogued area. In this paper we provide an overview of the data release with a focus on the processing of the LOFAR data and the characteristics of the resulting images. In two accompanying papers we provide the radio source associations and deblending and, where possible, the optical identifications of the radio sources together with the photometric redshifts and properties of the host galaxies. These data release papers are published together with a further $\sim$20 articles that highlight the scientific potential of LoTSS.

Historically, the blazar population has been poorly understood at low frequencies because survey sensitivity and angular resolution limitations have made it difficult to identify megahertz counterparts. We used the LOFAR Two-Metre Sky Survey (LoTSS) first data release value-added catalogue (LDR1) to study blazars in the low-frequency regime with unprecedented sensitivity and resolution. We identified radio counterparts to all $98$ known sources from the Third \textit{Fermi}-LAT Point Source Catalogue (3FGL) or Roma-BZCAT Multi-frequency Catalogue of Blazars ($5^{\mathrm{th}}$ edition) that fall within the LDR1 footprint. Only the 3FGL unidentified $\gamma$-ray sources (UGS) could not be firmly associated with an LDR1 source; this was due to source confusion. We examined the redshift and radio luminosity distributions of our sample, finding flat-spectrum radio quasars (FSRQs) to be more distant and more luminous than BL Lacertae objects (BL Lacs) on average. Blazars are known to have flat spectra in the gigahertz regime but we found this to extend down to $144$ MHz, where the radio spectral index, $\alpha$, of our sample is $-0.17 \pm 0.14$. For BL Lacs, $\alpha = -0.13 \pm 0.16$ and for FSRQs, $\alpha = -0.15 \pm 0.17$. We also investigated the radio-to-$\gamma$-ray connection for the $30$ $\gamma$-ray-detected sources in our sample. We find Pearson's correlation coefficient is $0.45$ ($p = 0.069$). This tentative correlation and the flatness of the spectral index suggest that the beamed core emission contributes to the low-frequency flux density. We compare our sample distribution with that of the full LDR1 on colour-colour diagrams, and we use this information to identify possible radio counterparts to two of the four UGS within the LDR1 field. We will refine our results as LoTSS continues.

The radio-loud/radio-quiet (RL/RQ) dichotomy in quasars is still an open question. Although it is thought that accretion onto supermassive black holes in the centre the host galaxies of quasars is responsible for some radio continuum emission, there is still a debate as to whether star formation or active galactic nuclei (AGN) activity dominate the radio continuum luminosity. To date, radio emission in quasars has been investigated almost exclusively using high-frequency observations in which the Doppler boosting might have an important effect on the measured radio luminosity, whereas extended structures, best observed at low radio frequencies, are not affected by the Doppler enhancement. We used a sample of quasars selected by their optical spectra in conjunction with sensitive and high-resolution low-frequency radio data provided by the LOw Frequency ARray (LOFAR) as part of the LOFAR Two-Metre Sky Survey (LoTSS) to investigate their radio properties using the radio loudness parameter ($\mathcal{R} = \frac{L_{\mathrm{144-MHz}}}{L_{\mathrm{i\,band}}}$). The examination of the radio continuum emission and RL/RQ dichotomy in quasars exhibits that quasars show a wide continuum of radio properties (i.e. no clear bimodality in the distribution of $\mathcal{R}$). Radio continuum emission at low frequencies in low-luminosity quasars is consistent with being dominated by star formation. We see a significant albeit weak dependency of $\mathcal{R}$ on the source nuclear parameters. For the first time, we are able to resolve radio morphologies of a considerable number of quasars. All these crucial results highlight the impact of the deep and high-resolution low-frequency radio surveys that foreshadow the compelling science cases for the Square Kilometre Array (SKA).

The massive evolved Wolf-Rayet stars sometimes occur in colliding-wind binary systems in which dust plumes are formed as a result of the collision of stellar winds. These structures are known to encode the parameters of the binary orbit and winds. Here, we report observations of a previously undiscovered Wolf-Rayet system, 2XMM J160050.7-514245, with a spectroscopically determined wind speed of $\approx$3400 km s$^{-1}$. In the thermal infrared, the system is adorned with a prominent $\approx$12$''$ spiral dust plume, revealed by proper motion studies to be expanding at only $\approx$570 km s$^{-1}$. As the dust and gas appear coeval, these observations are inconsistent with existing models of the dynamics of such colliding wind systems. We propose that this contradiction can be resolved if the system is capable of launching extremely anisotropic winds. Near-critical stellar rotation is known to drive such winds, suggesting this Wolf-Rayet system as a potential Galactic progenitor system for long-duration gamma-ray bursts.

Nearly everything we know about extrasolar planets to date comes from optical astronomy. While exoplanetary aurorae are predicted to be bright at low radio frequencies (< 1 GHz), we consider the effect of an exoplanet transit on radio emission from the host star. As radio emission from solar-like stars is concentrated in active regions, a planet occulting a starspot can cause a disproportionately deep transit which should be detectable with major radio arrays currently under development, such as the Square Kilometre Array (SKA). We calculate the radiometric sensitivity of the SKA stages and components, finding that SKA2-Mid can expect to detect transits around the very nearest solar-like stars and many cool dwarfs. The shape of this radiometric light curve will be affected by scintillation and lensing from the planet's magnetosphere and thereby encode magnetospheric parameters. Furthermore, these transits will also probe the distribution of stellar activity across a star's surface, and will help scrub out contamination from stellar activity on exoplanet transmission spectra and radial velocity spectra. This radio window on exoplanets and their host stars is therefore a valuable complement to existing optical tools.

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 shape of low-frequency radio continuum spectra of normal galaxies is not well understood, the key question being the role of physical processes such as thermal absorption in shaping them. In this work we take advantage of the LOFAR Multifrequency Snapshot Sky Survey (MSSS) to investigate such spectra for a large sample of nearby star-forming galaxies. Using the measured 150MHz flux densities from the LOFAR MSSS survey and literature flux densities at various frequencies we have obtained integrated radio spectra for 106 galaxies. The spectra are explained through the use of a three-dimensional model of galaxy radio emission, and radiation transfer dependent on the galaxy viewing angle and absorption processes. Spectra of our galaxies are generally flatter at lower compared to higher frequencies but as there is no tendency for the highly inclined galaxies to have more flattened low-frequency spectra, we argue that the observed flattening is not due to thermal absorption, contradicting the suggestion of Israel & Mahoney (1990). According to our modelled radio maps for M51-like galaxies, the free-free absorption effects can be seen only below 30MHz and in the global spectra just below 20MHz, while in the spectra of starburst galaxies, like M82, the flattening due to absorption is instead visible up to higher frequencies of about 150MHz. Locally, within galactic disks, the absorption effects are distinctly visible in M51-like galaxies as spectral flattening around 100-200MHz in the face-on objects, and as turnovers in the edge-on ones, while in M82-like galaxies there are strong turnovers at frequencies above 700MHz, regardless of viewing angle. Our modelling suggests that the weak spectral flattening observed in the nearby galaxies studied here results principally from synchrotron spectral curvature due to cosmic ray energy losses and propagation effects.

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.

PKS 1718$-$649 is one of the closest and most comprehensively studied candidates of a young active galactic nucleus (AGN) that is still embedded in its optical host galaxy. The compact radio structure, with a maximal extent of a few parsecs, makes it a member of the group of compact symmetric objects (CSO). Its environment imposes a turnover of the radio synchrotron spectrum towards lower frequencies, also classifying PKS 1718$-$649 as gigahertz-peaked radio spectrum (GPS) source. Its close proximity has allowed the first detection of extended X-ray emission in a GPS/CSO source with Chandra that is for the most part unrelated to nuclear feedback. However, not much is known about the nature of this emission. By co-adding all archival Chandra data and complementing these datasets with the large effective area of XMM-Newton, we are able to study the detailed physics of the environment of PKS 1718$-$649. Not only can we confirm that the bulk of the $\lesssim$kiloparsec-scale environment emits in the soft X-rays, but we also identify the emitting gas to form a hot, collisionally ionized medium. While the feedback of the central AGN still seems to be constrained to the inner few parsecs, we argue that supernovae are capable of producing the observed large-scale X-ray emission at a rate inferred from its estimated star formation rate.

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 order to determine the causes of kinematic asymmetry in the H$\alpha$ gas in the SAMI Galaxy Survey sample, we investigate the comparative influences of environment and intrinsic properties of galaxies on perturbation. We use spatially resolved H$\alpha$ velocity fields from the SAMI Galaxy Survey to quantify kinematic asymmetry ($\overline{v_{asym}}$) in nearby galaxies and environmental and stellar mass data from the GAMA survey. {We find that local environment, measured as distance to nearest neighbour, is inversely correlated with kinematic asymmetry for galaxies with $\mathrm{\log(M_*/M_\odot)}>10.0$, but there is no significant correlation for galaxies with $\mathrm{\log(M_*/M_\odot)}<10.0$. Moreover, low mass galaxies ($\mathrm{\log(M_*/M_\odot)}<9.0$) have greater kinematic asymmetry at all separations, suggesting a different physical source of asymmetry is important in low mass galaxies.} We propose that secular effects derived from gas fraction and gas mass may be the primary causes of asymmetry in low mass galaxies. High gas fraction is linked to high $\frac{\sigma_{m}}{V}$ (where $\sigma_m$ is H$\alpha$ velocity dispersion and $V$ the rotation velocity), which is strongly correlated with $\overline{v_{asym}}$, and galaxies with $\log(M_*/M_\odot)<9.0$ have offset $\overline{\frac{\sigma_{m}}{V}}$ from the rest of the sample. Further, asymmetry as a fraction of dispersion decreases for galaxies with $\log(M_*/M_\odot)<9.0$. Gas mass and asymmetry are also inversely correlated in our sample. We propose that low gas masses in dwarf galaxies may lead to asymmetric distribution of gas clouds, leading to increased relative turbulence.

We investigate the Tully-Fisher Relation (TFR) for a morphologically and kine- matically diverse sample of galaxies from the SAMI Galaxy Survey using 2 dimensional spatially resolved Halpha velocity maps and find a well defined relation across the stellar mass range of 8.0 < log(M*) < 11.5. We use an adaptation of kinemetry to parametrise the kinematic Halpha asymmetry of all galaxies in the sample, and find a correlation between scatter (i.e. residuals off the TFR) and asymmetry. This effect is pronounced at low stellar mass, corresponding to the inverse relationship between stellar mass and kinematic asymmetry found in previous work. For galaxies with log(M*) < 9.5, 25 +/- 3% are scattered below the root mean square (RMS) of the TFR, whereas for galaxies with log(M*) > 9.5 the fraction is 10 +/- 1% We use 'simulated slits' to directly compare our results with those from long slit spectroscopy and find that aligning slits with the photometric, rather than the kinematic, position angle, increases global scatter below the TFR. Further, kinematic asymmetry is correlated with misalignment between the photometric and kinematic position angles. This work demonstrates the value of 2D spatially resolved kinematics for accurate TFR studies; integral field spectroscopy reduces the underestimation of rotation velocity that can occur from slit positioning off the kinematic axis.

We present high sensitivity ($\sigma_P \simeq 0.6\,$mJy) polarimetric observations in seven bands, from $2.1$ to $38\,$GHz, of a complete sample of $104$ compact extragalactic radio sources brighter than $200\,$mJy at $20\,$GHz. Polarization measurements in six bands, in the range $5.5-38\,$GHz, for $53$ of these objects were reported by \citet{Galluzzi2017}. We have added new measurements in the same six bands for another 51 sources and measurements at $2.1\,$GHz for the full sample of $104$ sources. Also, the previous measurements at $18$, $24$, $33$ and $38\,$GHz were re-calibrated using the updated model for the flux density absolute calibrator, PKS1934-638, not available for the earlier analysis. The observations, carried out with the Australia Telescope Compact Array (ATCA), achieved a $90\%$ detection rate (at $5\sigma$) in polarization. $89$ of our sources have a counterpart in the $72$ to $231\,$MHz GLEAM survey \citep{HurleyWalker2017}, providing an unparalleled spectral coverage of $2.7$ decades of frequency for these sources. While the total intensity data from $5.5$ to $38\,$GHz could be interpreted in terms of single component emission, a joint analysis of more extended total intensity spectra presented here, and of the polarization spectra, reveals that over $90\%$ of our sources show clear indications of at least two emission components. We interpret this as an evidence of recurrent activity. Our high sensitivity polarimetry has allowed a $5\,\sigma$ detection of the weak circular polarization for $\sim 38\%$ of the dataset, and a deeper estimate of $20\,$GHz polarization source counts than has been possible so far.

We report the first astrophysical application of the technique of wide-field Interplanetary Scintillation (IPS) with the Murchison Widefield Array (MWA). This powerful technique allows us to identify and measure sub-arcsecond compact components in low-frequency radio sources across large areas of sky without the need for long-baseline interferometry or ionospheric calibration. We present the results of a five-minute observation of a 30x30 sq. deg MWA field at 162 MHz with 0.5 second time resolution. Of the 2550 continuum sources detected in this field, 302 (12 per cent) show rapid fluctuations caused by IPS. We find that at least 32% of bright low-frequency radio sources contain a sub-arcsec compact component that contributes over 40% of the total flux density. Perhaps surprisingly, peaked-spectrum radio sources are the dominant population among the strongly-scintillating, low-frequency sources in our sample. While gamma-ray AGN are generally compact, flat-spectrum radio sources at higher frequencies, the 162 MHz properties of many of the Fermi blazars in our field are consistent with a compact component embedded within more extended low-frequency emission. The detection of a known pulsar in our field shows that the wide-field IPS technique is at the threshold of sensitivity needed to detect new pulsars using image plane analysis, and scaling the current MWA sensitivity to that expected for SKA-low implies that large IPS-based pulsar searches will be feasible with SKA. Calibration strategies for the SKA require a better knowledge of the space density of compact sources at low radio frequencies, which IPS observations can now provide.

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.

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

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

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

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

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

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

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.

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.

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.

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.

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

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

Using the new wideband capabilities of the Australia Telescope Compact Array (ATCA), we obtain spectra for PKS 1718-649, a well-known gigahertz-peaked spectrum radio source. The observations, between approximately 1 and 10 GHz over three epochs spanning approximately 21 months, reveal variability both above the spectral peak at ~3 GHz and below the peak. The combination of the low and high frequency variability cannot be easily explained using a single absorption mechanism, such as free-free absorption or synchrotron self-absorption. We find that the PKS 1718-649 spectrum and its variability are best explained by variations in the free-free optical depth on our line-of-sight to the radio source at low frequencies (below the spectral peak) and the adiabatic expansion of the radio source itself at high frequencies (above the spectral peak). The optical depth variations are found to be plausible when X-ray continuum absorption variability seen in samples of Active Galactic Nuclei is considered. We find that the cause of the peaked spectrum in PKS 1718-649 is most likely due to free-free absorption. In agreement with previous studies, we find that the spectrum at each epoch of observation is best fit by a free-free absorption model characterised by a power-law distribution of free-free absorbing clouds. This agreement is extended to frequencies below the 1 GHz lower limit of the ATCA by considering new observations with Parkes at 725 MHz and 199 MHz observations with the newly operational Murchison Widefield Array. These lower frequency observations argue against families of absorption models (both free-free and synchrotron self-absorption) that are based on simple homogenous structures.

M giants are among the longest-period pulsating stars which is why their studies were traditionally restricted to analyses of low-precision visual observations, and more recently, accurate ground-based data. Here we present an overview of M giant variability on a wide range of time-scales (hours to years), based on analysis of thirteen quarters of Kepler long-cadence observations (one point per every 29.4 minutes), with a total time-span of over 1000 days. About two-thirds of the sample stars have been selected from the ASAS-North survey of the Kepler field, with the rest supplemented from a randomly chosen M giant control sample. We first describe the correction of the light curves from different quarters, which was found to be essential. We use Fourier analysis to calculate multiple frequencies for all stars in the sample. Over 50 stars show a relatively strong signal with a period equal to the Kepler-year and a characteristic phase dependence across the whole field-of-view. We interpret this as a so far unidentified systematic effect in the Kepler data. We discuss the presence of regular patterns in the distribution of multiple periodicities and amplitudes. In the period-amplitude plane we find that it is possible to distinguish between solar-like oscillations and larger amplitude pulsations which are characteristic for Mira/SR stars. This may indicate the region of the transition between two types of oscillations as we move upward along the giant branch.

Detailed knowledge of the mesopheric sodium layer characteristics is crucial to estimate and optimize the performance of Laser Guide Star (LGS) assisted Adaptive Optics (AO) systems. In this paper, we present an analysis of two sets of data on the mesospheric sodium layer. The first set comes from a laser experiment that was carried out at Cerro Tololo to monitor the abundance and altitude of the mesospheric sodium in 2001, during six runs covering a period of one year. This data is used to derive the mesospheric sodium column density, the sodium layer thickness and the temporal behavior of the sodium layer mean altitude. The second set of data was gathered during the first year of the Gemini MCAO System (GeMS) commissioning and operations. GeMS uses five LGS to measure and compensate for atmospheric distortions. Analysis of the LGS wavefront sensor data provides information about the sodium photon return and the spot elongation seen by the WFS. All these parameters show large variations on a yearly, nightly and hourly basis, affecting the LGS brightness, shape and mean altitude. The sodium photon return varies by a factor of three to four over a year, and can change by a factor of two over a night. In addition, the comparison of the photon returns obtained in 2001 with those measured a decade later using GeMS shows a significant difference in laser format efficiencies. We find that the temporal power spectrum of the sodium mean altitude follows a linear trend, in good agreement with the results reported by Pfrommer & Hickson (2010).

We present detailed analysis of the transient X-ray source 2XMMi J003833.3+402133 detected by XMM-Newton in January 2008 during a survey of M 31. The X-ray spectrum is well fitted by either a steep power law plus a blackbody model or a double blackbody model. Prior observations with XMM-Newton, Chandra, Swift and ROSAT spanning 1991 to 2007, as well as an additional Swift observation in 2011, all failed to detect this source. No counterpart was detected in deep optical imaging with the Canada France Hawaii Telescope down to a 3sigma lower limit of g = 26.5 mag. This source has previously been identified as a black hole X-ray binary in M 31. While this remains a possibility, the transient behaviour, X-ray spectrum, and lack of an optical counterpart are equally consistent with a magnetar interpretation for 2XMMi J003833.3+402133. The derived luminosity and blackbody emitting radius at the distance of M 31 argue against an extragalactic location, implying that if it is indeed a magnetar it is located within the Milky Way but 22deg out of the plane. The high Galactic latitude could be explained if 2XMMi J003833.3+402133 were an old magnetar, or if its progenitor was a runaway star that traveled away from the plane prior to going supernova.