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.

We compute the spherically-averaged power spectrum from four seasons of data obtained for the Epoch of Reionisation (EoR) project observed with the Murchison Widefield Array (MWA). We measure the EoR power spectrum over $k= 0.07-3.0~h$Mpc$^{-1}$ at redshifts $z=6.5-8.7$. The largest aggregation of 110 hours on EoR0 high-band (3,340 observations), yields a lowest measurement of (43~mK)$^2$ = 1.8$\times$10$^3$ mK$^2$ at $k$=0.14~$h$Mpc$^{-1}$ and $z=6.5$ (2$\sigma$ thermal noise plus sample variance). Using the Real-Time System to calibrate and the CHIPS pipeline to estimate power spectra, we select the best observations from the central five pointings within the 2013--2016 observing seasons, observing three independent fields and in two frequency bands. This yields 13,591 2-minute snapshots (453 hours), based on a quality assurance metric that measures ionospheric activity. We perform another cut to remove poorly-calibrated data, based on power in the foreground-dominated and EoR-dominated regions of the two-dimensional power spectrum, reducing the set to 12,569 observations (419 hours). These data are processed in groups of 20 observations, to retain the capacity to identify poor data, and used to analyse the evolution and structure of the data over field, frequency, and data quality. We subsequently choose the cleanest 8,935 observations (298 hours of data) to form integrated power spectra over the different fields, pointings and redshift ranges.

The compact configuration of Phase II of the Murchison Widefield Array (MWA) consists of both a redundant subarray and pseudo-random baselines, offering unique opportunities to perform sky-model and redundant interferometric calibration. The highly redundant hexagonal cores give improved power spectrum sensitivity. In this paper, we present the analysis of nearly 40 hours of data targeting one of the MWA's EoR fields observed in 2016. We use both improved analysis techniques presented in Barry et al. (2019) as well as several additional techniques developed for this work, including data quality control methods and interferometric calibration approaches. We show the EoR power spectrum limits at redshift 6.5, 6.8 and 7.1 based on our deep analysis on this 40-hour data set. These limits span a range in $k$ space of $0.18$ $h$ $\mathrm{Mpc^{-1}}$ $<k<1.6$ $h$ $\mathrm{Mpc^{-1}}$, with a lowest measurement of $\Delta^2\leqslant2.39\times 10^3$ $\mathrm{mK}^2$ at $k=0.59$ $h$ $\mathrm{Mpc^{-1}}$ and $z=6.5$.

The Giant Radio Array for Neutrino Detection (GRAND) is a planned large-scale observatory of ultra-high-energy (UHE) cosmic particles, with energies exceeding 10^8 GeV. Its goal is to solve the long-standing mystery of the origin of UHE cosmic rays. To do this, GRAND will detect an unprecedented number of UHE cosmic rays and search for the undiscovered UHE neutrinos and gamma rays associated to them with unmatched sensitivity. GRAND will use large arrays of antennas to detect the radio emission coming from extensive air showers initiated by UHE particles in the atmosphere. Its design is modular: 20 separate, independent sub-arrays, each of 10 000 radio antennas deployed over 10 000 km^2. A staged construction plan will validate key detection techniques while achieving important science goals early. Here we present the science goals, detection strategy, preliminary design, performance goals, and construction plans for GRAND.

The overwhelming foreground contamination is one of the primary impediments to probing the Epoch of Reionization (EoR) through measuring the redshifted 21 cm signal. Among various foreground components, radio halos are less studied and their impacts on the EoR observations are still poorly understood. In this work, we employ the Press-Schechter formalism, merger-induced turbulent re-acceleration model, and the latest SKA1-Low layout configuration to simulate the SKA "observed" images of radio halos. We calculate the one-dimensional power spectra from simulated images and find that radio halos can be about $10^4$, $10^3$ and $10^{2.5}$ times more luminous than the EoR signal on scales of $0.1\,\text{Mpc}^{-1} < k < 2\,\text{Mpc}^{-1}$ in the 120-128, 154-162, and 192-200 MHz bands, respectively. By examining the two-dimensional power spectra inside properly defined EoR windows, we find that the power leaked by radio halos can still be significant, as the power ratios of radio halos to the EoR signal on scales of $0.5\,\text{Mpc}^{-1} \lesssim k \lesssim 1\,\text{Mpc}^{-1}$ can be up to about 230-800%, 18-95%, and 7-40% in the three bands, when the 68% uncertainties caused by the variation of the number density of bright radio halos are considered. Furthermore, we find that radio halos located inside the far side-lobes of the station beam can also impose strong contamination within the EoR window. In conclusion, we argue that radio halos are severe foreground sources and need serious treatments in future EoR experiments.

We apply two methods to estimate the 21~cm bispectrum from data taken within the Epoch of Reionisation (EoR) project of the Murchison Widefield Array (MWA). Using data acquired with the Phase II compact array allows a direct bispectrum estimate to be undertaken on the multiple redundantly-spaced triangles of antenna tiles, as well as an estimate based on data gridded to the $uv$-plane. The direct and gridded bispectrum estimators are applied to 21 hours of high-band (167--197~MHz; $z$=6.2--7.5) data from the 2016 and 2017 observing seasons. Analytic predictions for the bispectrum bias and variance for point source foregrounds are derived. We compare the output of these approaches, the foreground contribution to the signal, and future prospects for measuring the bispectra with redundant and non-redundant arrays. We find that some triangle configurations yield bispectrum estimates that are consistent with the expected noise level after 10 hours, while equilateral configurations are strongly foreground-dominated. Careful choice of triangle configurations may be made to reduce foreground bias that hinders power spectrum estimators, and the 21~cm bispectrum may be accessible in less time than the 21~cm power spectrum for some wave modes, with detections in hundreds of hours.

We explore methods for robust estimation of the 21 cm signal from the Epoch of Reionisation (EoR). A Kernel Density Estimator (KDE) is introduced for measuring the spatial temperature fluctuation power spectrum from the EoR. The KDE estimates the underlying probability distribution function of fluctuations as a function of spatial scale, and contains different systematic biases and errors to the typical approach to estimating the fluctuation power spectrum. Extraction of histograms of visibilities allows moments analysis to be used to discriminate foregrounds from 21 cm signal and thermal noise. We use the information available in the histograms, along with the statistical dis-similarity of foregrounds from two independent observing fields, to robustly separate foregrounds from cosmological signal, while making no assumptions about the Gaussianity of the signal. Using two independent observing fields to robustly discriminate signal from foregrounds is crucial for the analysis presented in this paper. We apply the techniques to 13 hours of Murchison Widefield Array (MWA) EoR data over two observing fields. We compare the output to that obtained with a comparative power spectrum estimation method, and demonstrate the reduced foreground contamination using this approach. Using the second moment obtained directly from the KDE distribution functions yields a factor of 2-3 improvement in power for k < 0.3hMpc^{-1} compared with a matched delay space power estimator, while weighting data by additional statistics does not offer significant improvement beyond that available for thermal noise-only weights.

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 Square Kilometre Array (SKA) will be the first low-frequency instrument with the capability to directly image the structures of the Epoch of Reionization (EoR). Indeed, deep imaging of the EoR over 5 targeted fields of 20 square degrees each has been selected as the highest priority science objective for SKA1, probably the most exciting scientific outcome to be expected for the SKA1-Low. Aiming at preparing for this highly challenging observation, we perform an extensive pre-selection of the 'quietest' and 'cleanest' candidate fields in the southern sky to be suited for deep imaging of the EoR using existing catalogs and observations over a broad band. The candidate fields should meet a number of strict criteria to avoid contamination from foreground structures and sources such as the Galactic plane, LMC/SMC, clusters of galaxies and bright radio sources. The candidate fields should also exhibit both the lowest average surface brightness and smallest variance to ensure uniformity and high quality deep imaging over the fields. We use two sets of field of view in our selection: the exactly same size of 20 square degrees as defined in current SKA1 documentation and 40 square degrees - an area of twice big as the EoR imaging fields. The choice of a larger area allows us to maximally reduce the environmental effect from radio structures/sources in the vicinity of the candidate field. Our selection yields a sample of 29 (9) `ideal' fields of 20 (40) square degrees in the southern sky that could be targeted for deep imaging of the EoR. We have launched a campaign to observe these candidate fields with some of the SKA pathfinders to set further constraints on these candidates. These pre-selected candidate fields can provide a valuable guidance for eventual selection of the five targeted fields for imaging the EoR structures with the future SKA1-low.

We choose the bright compact group HCG 62, which was found to exhibit both excess X-ray emission and high Fe abundance to the southwest of its core, as an example to study the impact of mergers on chemical enrichment in the intragroup medium. We first reanalyze the high-quality Chandra and XMM-Newton archive data to search for the evidence for additional SN II yields, which is expected as a direct result of the possible merger-induced starburst. We reveal that, similar to the Fe abundance, the Mg abundance also shows a high value in both the innermost region and the southwest substructure, forming a high-abundance plateau, meanwhile all the SN Ia and SN II yields show rather flat distributions in $>0.1r_{200}$ in favor of an early enrichment. Then we carry out a series of idealized numerical simulations to model the collision of two initially isolated galaxy groups by using the TreePM-SPH GADGET-3 code. We find that the observed X-ray emission and metal distributions, as well as the relative positions of the two bright central galaxies with reference to the X-ray peak, can be well reproduced in a major merger with a mass ratio of 3 when the merger-induced starburst is assumed. The `best-match' snapshot is pinpointed after the third pericentric passage when the southwest substructure is formed due to gas sloshing. By following the evolution of the simulated merging system, we conclude that the effects of such a major merger on chemical enrichment are mostly restricted within the core region when the final relaxed state is reached.

The Torpedo cluster, SPT-CL J0245-5302 (S0295) is a massive, merging cluster at a redshift of z = 0.300, which exhibits a strikingly similar morphology to the Bullet cluster 1E 0657-55.8 (z = 0.296), including a classic bow shock in the cluster's intra-cluster medium revealed by Chandra X-ray observations. We present Australia Telescope Compact Array data centred at 2.1 GHz and Murchison Widefield Array data at frequencies between 72 MHz and 231 MHz which we use to study the properties of the cluster. We characterise a number of discrete and diffuse radio sources in the cluster, including the detection of two previously unknown radio relics on the cluster periphery. The average spectral index of the diffuse emission between 70 MHz and 3.1 GHz is $\alpha=-1.63_{-0.10}^{+0.10}$ and a radio-derived Mach number for the shock in the west of the cluster is calculated as M = 2.04. The Torpedo cluster is thus as a double relic system at moderate redshift.

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

The 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 Giant Radio Array for Neutrino Detection (GRAND) aims at detecting ultra-high energy extraterrestrial neutrinos via the extensive air showers induced by the decay of tau leptons created in the interaction of neutrinos under the Earth's surface. Consisting of an array of $\sim10^5$ radio antennas deployed over $\sim 2\times10^5\,\rm {km}^2$, GRAND plans to reach, for the first time, an all-flavor sensitivity of $\sim1.5\times10^{-10} \,\rm GeV\, cm^{-2} \,s^{-1}\, sr^{-1}$ above $5\times10^{17}$ eV and a sub-degree angular resolution, beyond the reach of other planned detectors. We describe here preliminary designs and simulation results, plans for the ongoing, staged approach to the construction of GRAND, and the rich research program made possible by GRAND's design sensitivity and angular resolution.

We detect and characterise extended, diffuse radio emission from galaxy clusters at 168 MHz within the Epoch of Reionization 0-hour field; a $45^\circ \times 45^\circ$ region of the southern sky centred on R. A. $= 0^\circ$, decl. $=-27^\circ$. We detect 31 diffuse radio sources; 3 of which are newly detected haloes in Abell 0141, Abell 2811, and Abell S1121; 2 newly detected relics in Abell 0033 and Abell 2751; 5 new halo candidates and a further 5 new relic candidates. Further, we detect a new phoenix candidate in Abell 2556 as well as 2 candidate dead radio galaxies at the centres of Abell 0122 and Abell S1136 likely associated with the brightest cluster galaxies. Beyond this we find 2 clusters with unclassifiable, diffuse steep-spectrum emission as well as a candidate double relic system associated with RXC J2351.0-1934. We present measured source properties such as their integrated flux densities, spectral indices, and sizes where possible. We find several of the diffuse sources to be ultra-steep including the halo in Abell 0141 which has $\alpha \leq -2.1 \pm 0.1$ making it one of the steepest halos known. Finally, we compare our sample of haloes with previously detected haloes and revisit established scaling relations of the radio halo power ($P_{1.4}$) with the cluster X-ray luminosity ($L_{\mathrm{X}}$) and mass ($M_{500}$). We find consistent fitting parameters for assumed power law relationships, and find that the $P_{1.4}-L_{\mathrm{X}}$ has less raw scatter than the corresponding $P_{1.4}-M_{500}$ despite inhomogeneous $L_{\mathrm{X}}$ measurements. These scaling relation properties are consistent with a sample of only non-cool core clusters hosting radio haloes.

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.

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 a catalog of 624 radio sources detected around the North Celestial Pole (NCP) with the 21 Centimeter Array (21CMA), a radio interferometer dedicated to the statistical measurement of the epoch of reionization (EoR). The data are taken from a 12 h observation made on 2013 April 13, with a frequency coverage from 75 to 175 MHz and an angular resolution of ~ 4 arcmin. The catalog includes flux densities at eight sub-bands across the 21CMA bandwidth and provides the in-band spectral indices for the detected sources. To reduce the complexity of interferometric imaging from the so-called "w" term and ionospheric effects, the present analysis are restricted to the east-west baselines within 1500 m only. The 624 radio sources are found within 5 degrees around the NCP down to ~ 0.1 Jy. Our source counts are compared, and also exhibit a good agreement, with deep low-frequency observations made recently with the GMRT and MWA. In particular, for fainter radio sources below ~ 1 Jy, we find a flattening trend of source counts towards lower frequencies. While the thermal noise (~0.4 mJy) is well controlled to below the confusion limit, the dynamical range (~10^4) and sensitivity of current 21CMA imaging is largely limited by calibration and deconvolution errors, especially the grating lobes of very bright sources, such as 3C061.1, in the NCP field which result from the regular spacings of the 21CMA. We note that particular attention should be paid to the extended sources, and their modeling and removals may constitute a large technical challenge for current EoR experiments. Our analysis may serve as a useful guide to design of next generation low-frequency interferometers like the Square Kilometre Array.

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.

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.

Detection and mitigation of radio frequency interference (RFI) is the first and also the key step for data processing in radio observations, especially for ongoing low frequency radio experiments towards the detection of the cosmic dawn and epoch of reionization (EoR). In this paper we demonstrate the technique and efficiency of RFI identification and mitigation for the 21 Centimeter Array (21CMA), a radio interferometer dedicated to the statistical measurement of EoR. For terrestrial, man-made RFI, we concentrate mainly on a statistical approach by identifying and then excising non-Gaussian signatures, in the sense that the extremely weak cosmic signal is actually buried under thermal and therefore Gaussian noise. We also introduce the so-called visibility correlation coefficient instead of conventional visibility, which allows a further suppression of rapidly time-varying RFI. Finally, we briefly discuss removals of the sky RFI, the leakage of sidelobes from off-field strong radio sources with time-invariant power and a featureless spectrum. It turns out that state of the art technique should allow us to detect and mitigate RFI to a satisfactory level in present low frequency interferometer observations such as those acquired with the 21CMA, and the accuracy and efficiency can be greatly improved with the employment of low-cost, high-speed computing facilities for data acquisition and processing.

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

The redshifted 21cm signal of neutral hydrogen from the epoch of reionization (EoR) is extremely weak and its first detection is therefore expected to be statistical with first-generation low-frequency radio interferometers. In this letter we propose a method to extract the angular power spectrum of EoR from the visibility correlation coefficients p_{ij}(u,v), instead of the visibilities V_{ij}(u,v) measured directly by radio interferometers in conventional algorithm. The visibility correlation coefficients are defined as p_{ij}(u,v)=V_{ij}(u,v)/\sqrt{|V_{ii}||V_{jj}|} by introducing the auto-correlation terms V_{ii} and V_{jj} such that the angular power spectrum C_{\ell} can be obtained through C_{\ell}=4pi^2T_0^2<|p_{ij}(u,v)|^2>, independently of the primary beams of antennas. This also removes partially the influence of receiver gains in the measurement of C_{\ell} because the amplitudes of the gains cancel each other out in the statistical average operation of <|p_{ij}(u,v)|^2>.We use the average system temperature T_0 as a calibrator of C_{\ell}, which is dominated by the Milky Way and extragalactic sources in our interested frequency range below 200 MHz. Finally we demonstrate the feasibility of the novel method using the simulated sky maps as targets and the 21 CentiMeter Array (21CMA) as interferometer.

Presence of spectral curvature due to synchrotron self-absorption (SSA) of extragalactic radio sources may break down the spectral smoothness feature-the premise that bright radio foreground can be successfully removed in the 21cm experiments of searching for the epoch of reionization (EOR).We present a quantitative estimate of the effect on the measurement of the angular power spectrum of the low-frequency sky,incorporating a phenomenological model,characterized by the fraction of radio sources with turnover frequencies in 100-1000MHz range and a broken power law for the spectral transition around turnover frequencies nu_m,into the simulated radio sources over a small sky area of 10x10 deg^2.We compare statistically the changes in their residual maps with/without inclusion of the SSA after the bright sources of S_150MHz>=100mJy are excised and the best-fitted polynomials in frequency domain on each pixel are further subtracted.It has been shown that the effect of SSA on the detection of EOR depends sensitively on the spectral profiles of radio sources around the turnover frequencies: A hard transition model described by the broken power law with the turnover of spectral index at nu_m would leave pronounced imprints on the residual background and cause serious confusion with the EOR signal.However,the spectral signatures on the angular power spectrum of extragalactic foreground generated by a soft transition model,in which the rise and fall power laws of spectral distribution around nu_m are connected through a smooth transition spanning >=200 MHz in characteristic width,can be fitted and consequently subtracted by employment of polynomials to an acceptable degree(delta T<1mK).As this latter scenario seems to be favored by both theoretical expectation and radio spectral observations,we conclude that the influence of SSA on the 21cm experiments is probably very minor.

In this work we simulate the $50-200$ MHz radio sky that is constrained in the field of view ($5^{\circ}$ radius) of the 21 Centimeter Array (21CMA), by carrying out Monte-Carlo simulations to model redshifted cosmological reionization signals and strong contaminating foregrounds, including emissions from our Galaxy, galaxy clusters, and extragalactic point sources. As an improvement of previous works, we consider in detail not only random variations of morphological and spectroscopic parameters within the ranges allowed by multi-band observations, but also evolution of radio halos in galaxy clusters, assuming that relativistic electrons are re-accelerated in the ICM in merger events and lose energy via both synchrotron emission and inverse Compton scattering with CMB photons. By introducing a new approach designed on the basis of independent component analysis (ICA) and wavelet detection algorithm, we prove that, with a cumulative observation of one month with the 21CMA array, about $80\%$ of galaxy clusters with central brightness temperatures of $> 10~{\rm K}$ at 65 MHz can be safely identified and separated from the overwhelmingly bright foreground. We find that the morphological and spectroscopic distortions are extremely small as compared to the input simulated clusters, and the reduced $\chi^2$ of brightness temperature profiles and spectra are controlled to be $\lesssim 0.5$ and $\lesssim 1.3$, respectively. These results robustly indicate that in the near future a sample of dozens of bright galaxy clusters will be disentangled from the foreground in 21CMA observations, the study of which will greatly improve our knowledge about cluster merger rates, electron acceleration mechanisms in cluster radio halos, and magnetic field in the ICM.

Measurement of the 21cm hyperfine transition of neutral hydrogen provides a unique probe of the epoch of reionization and the Dark Ages. Three major mechanisms are believed to dominate the radiation process: emission from neutral hydrogen surrounding the ionized bubbles of first galaxies and/or quasars, emission from neutral hydrogen inside minihalos, and absorption of diffuse neutral hydrogen against the cosmic microwave background. In the present work, by simply combining the existing analytic models for the three mechanisms, we investigate the contribution of cross-correlation between these three components to the total 21cm angular power spectrum, in the sense that neutral hydrogen associated with different radiation processes traces the large-scale structures of underlying density perturbations. While the overall 21cm power spectrum remains almost unchanged with the inclusion of the cross-correlations, the cross-correlation may play a key role in the determination of the 21cm power spectrum during the transition of 21cm radiation from emission-dominated phase to absorption-dominated phase at redshift z~20. A significant suppression in the 21cm angular power spectrum during this transition is anticipated as the result of negative contribution of the cross-correlation between the absorption of diffuse neutral hydrogen and the emission components. Therefore, an accurate prediction of the cosmic 21cm power spectrum should take the cross-correlation into account especially at the transition phase.