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

Measuring the properties of extragalactic magnetic fields through the effect of Faraday rotation provides a means to understand the origin and evolution of cosmic magnetism. Here we use data from the LOFAR Two-Metre Sky Survey (LoTSS) to calculate the Faraday rotation measure (RM) of close pairs of extragalactic radio sources. By considering the RM difference ($\Delta$RM) between physical pairs (e.g. double-lobed radio galaxies) and non-physical pairs (i.e. close projected sources on the sky), we statistically isolate the contribution of extragalactic magnetic fields to $\Delta$RM along the line of sight between non-physical pairs. From our analysis, we find no significant difference between the $\Delta$RM distributions of the physical and non-physical pairs, limiting the excess Faraday rotation contribution to $< 1.9$ rad/m$^2$ ($\sim$$95\%$ confidence). We use this limit with a simple model of an inhomogeneous universe to place an upper limit of 4 nG on the cosmological co-moving magnetic field strength on Mpc scales. We also compare the RM data with a more realistic suite of cosmological MHD simulations, that explore different magnetogenesis scenarios. Both magnetization of the large scale structure by astrophysical processes such as galactic and AGN outflows, and simple primordial scenarios with seed magnetic field strengths $< 0.5$ nG cannot be rejected by the current data; while stronger primordial fields or models with dynamo amplification in filaments are disfavoured.

Hybrid morphology radio sources (HyMoRS) are a rare group of radio galaxies in which differing Fanaroff & Riley morphologies (FR I/II) are observed for each of the two lobes. While they potentially provide insights into the formation of lobe structure, particle acceleration, and the FR dichotomy, previous work on HyMoRS has mainly been limited to low-resolution studies, searches for new candidates, and milliarcsecond-scale VLBI observations of the core region. In this paper, we use new multi-array configuration Very Large Array (VLA) observations between 1 and 8 GHz to determine the morphology of HyMoRS on arcsecond scales and perform the first well-resolved spectral study of these unusual sources. We find that while the apparent FR I lobe is centre-brightened, this is the result of a compact acceleration region resembling a hotspot with a spectrum more consistent with an FR II ("strong-flavour") jet. We find that the spectra of the apparent FR I lobes are not similar to their classical counterparts and are likely the result of line-of-sight mixing of plasma across a range of spectral ages. We consider possible mechanisms that could lead to the formation of HyMoRS under such conditions, including environment asymmetry and restarted sources, concluding through the use of simple modelling that HyMoRS are the result of orientation effects on intrinsically FR II sources with lobes non-parallel to the inner jet.

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.

Faraday rotation measures (RMs) of extragalactic radio sources provide information on line-of-sight magnetic fields, including contributions from our Galaxy, source environments, and the intergalactic medium (IGM). Looking at differences in RMs, $\Delta$RM, between adjacent sources on the sky can help isolate these different components. In this work, we classify adjacent polarized sources in the NVSS as random or physical pairs. We recompute and correct the uncertainties in the NVSS RM catalog, since these were significantly overestimated. Our sample contains 317 physical and 5111 random pairs, all with Galactic latitudes $|b|\ge20^{\circ}$, polarization fractions $\ge2\%$, and angular separations between $1.^{'},$ and $20^{'}$. We find an rms $\Delta$RM of $14.9\pm0.4\,$rad m$^{-2}$ and $4.6\pm1.1\,$rad m$^{-2}$ for random and physical pairs, respectively. This means polarized extragalactic sources that are close on the sky, but at different redshifts, have larger differences in RM than two components of one source. This difference of $\sim10\,$rad m$^{-2}$ is significant at $5\sigma$, and persists in different data subsamples. While there have been other statistical studies of $\Delta$RM between adjacent polarized sources, this is the first unambiguous demonstration that some of this RM difference must be extragalactic, thereby providing a firm upper limit on the RM contribution of the IGM. If the $\Delta$RMs originate local to the sources, then the local magnetic field difference between random sources is a factor of two larger than between components of one source. Alternatively, attributing the difference in $\Delta$RMs to the intervening IGM yields an upper limit on the IGM magnetic field strength of $40\,$nG.

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

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.

We made two new sensitive (rms noise sigma_n ~ 1 microJy/beam) high resolution (theta = 3.0" and theta = 0.66" FWHM) S--band (2 < nu < 4 GHz) images covering a single JVLA primary beam (FWHM ~ 14') centered on J2000 RA = 10 46, Dec = 59 01 in the Lockman Hole. These images yielded a catalog of 792 radio sources, 97.7 +/- 0.8% of which have infrared counterparts stronger than S ~ 2 microJy at lambda = 4.5 micron. About 91% of the radio sources found in our previously published, comparably sensitive low resolution (theta = 8" FWHM) image covering the same area were also detected at 0.66" resolution, so most radio sources with S_3GHz >~ 5 microJy have angular structure phi <~ 0.66". The ratios of peak brightness in the 0.66" and 3" images have a distribution indicating that most microJy radio sources are quite compact, with a median Gaussian angular diameter <phi> = 0.3" +/- 0.1" FWHM and an rms scatter sigma_phi <~ 0.3" of individual sizes. Most of our microJy radio sources obey the tight far-infrared/radio correlation, indicating that they are powered by star formation. The median effective angular radius enclosing half the light emitted by an exponential disk is <rho_e> ~ <phi>/2.43 ~ 0.12", so the median effective radius of star-forming galaxies at redshifts z~1 is <r_e> ~ 1.0 kpc.

We summarize the radio synchrotron background workshop that took place July 19-21, 2017 at the University of Richmond. This first scientific meeting dedicated to the topic was convened because current measurements of the diffuse radio monopole reveal a surface brightness that is several times higher than can be straightforwardly explained by known Galactic and extragalactic sources and processes, rendering it by far the least well understood photon background at present. It was the conclusion of a majority of the participants that the radio monopole level is at or near that reported by the ARCADE 2 experiment and inferred from several absolutely calibrated zero level lower frequency radio measurements, and unanimously agreed that the production of this level of surface brightness, if confirmed, represents a major outstanding question in astrophysics. The workshop reached a consensus on the next priorities for investigations of the radio synchrotron background.

We present the largest ever sample of radio polarization properties for $z>4$ sources, with 14 sources having significant polarization detections. Using wideband data from the Karl G. Jansky Very Large Array, we obtained the rest-frame total intensity and polarization properties of 37 radio sources, nine of which have spectroscopic redshifts in the range $1 \le z \le 1.4$, with the other 28 having spectroscopic redshifts in the range $3.5 \le z \le 6.21$. Fits are performed for the Stokes $I$ and fractional polarization spectra, and Faraday rotation measures are derived using Rotation measure synthesis and $QU$ fitting. Using archival data of 476 polarized sources, we compare high redshift ($z>3$) source properties to a $15\,$GHz rest-frame luminosity matched sample of low redshift ($z<3$) sources to investigate if the polarization properties of radio sources at high redshifts are intrinsically different than those at low redshift. We find a mean of the rotation measure absolute values, corrected for Galactic rotation, of $50 \pm 22\,$rad m$^{-2}$ for $z>3$ sources and $57 \pm 4\,$rad m$^{-2}$ for $z<3$. Although there is some indication of lower intrinsic rotation measures at high-$z$ possibly due to higher depolarization from the high density environments, using several statistical tests we detect no significant difference between low and high redshift sources. Larger samples are necessary to determine any true physical difference.

We set limits on the presence of the synchrotron cosmic web through the cross-correlation of the 2.3 GHz S-PASS survey with a model of the local cosmic web derived from constrained magnetohydrodynamic (MHD) simulations. The MHD simulation assumes cosmologically seeded magnetic fields amplified during large-scale structure formation, and a population of relativistic electrons/positrons from proton-proton collisions within the intergalactic medium. We set a model-dependent 3$\sigma$ upper limit on the synchrotron surface brightness of 0.16 mJy arcmin$^{-2}$ at 2.3 GHz in filaments. Extrapolating from magnetic field maps created from the simulation, we infer an upper limit (density-weighted) magnetic field of 0.03 (0.13) $\mu$G in filaments at the current epoch, and a limit on the primordial magnetic field (PMF) of B$_{PMF}$~1.0 nG.

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

We present the results of the first, deep ALMA imaging covering the full 4.5 sq arcmin of the Hubble Ultra Deep Field (HUDF) as previously imaged with WFC3/IR on HST. Using a mosaic of 45 pointings, we have obtained a homogeneous 1.3mm image of the HUDF, achieving an rms sensitivity of 35 microJy, at a resolution of 0.7 arcsec. From an initial list of ~50 >3.5sigma peaks, a rigorous analysis confirms 16 sources with flux densities S(1.3) > 120 microJy. All of these have secure galaxy counterparts with robust redshifts (<z> = 2.15), and 12 are also detected at 6GHz in new deep JVLA imaging. Due to the wealth of supporting data in this unique field, the physical properties of the ALMA sources are well constrained, including their stellar masses (M*) and UV+FIR star-formation rates (SFR). Our results show that stellar mass is the best predictor of SFR in the high-z Universe; indeed at z > 2 our ALMA sample contains 7 of the 9 galaxies in the HUDF with M* > 2 x 10^10 Msun and we detect only one galaxy at z > 3.5, reflecting the rapid drop-off of high-mass galaxies with increasing redshift. The detections, coupled with stacking, allow us to probe the redshift/mass distribution of the 1.3-mm background down to S(1.3) ~ 10 micro-Jy. We find strong evidence for a steep `main sequence' for star-forming galaxies at z ~ 2, with SFR \propto M* and a mean specific SFR = 2.2 /Gyr. Moreover, we find that ~85% of total star formation at z ~ 2 is enshrouded in dust, with ~65% of all star formation at this epoch occurring in high-mass galaxies (M* > 2 x 10^10 Msun), for which the average obscured:unobscured SF ratio is ~200. Finally, we combine our new ALMA results with the existing HST data to revisit the cosmic evolution of star-formation rate density; we find that this peaks at z ~ 2.5, and that the star-forming Universe transits from primarily unobscured to primarily obscured thereafter at z ~ 4.

This is the second of two papers describing the observations and source catalogues derived from sensitive 3-GHz images of the Lockman Hole North using the Karl G. Jansky Very Large Array. We describe the reduction and cataloguing process, which yielded an image with 8 arcsec resolution and instrumental noise of \sigma$_{\rm n}$=1.0 \mu Jy beam$^{-1}$ rms (before primary beam corrections) and a catalogue of 558 sources detected above 5\sigma$_{\rm n}$. We include details of how we estimate source spectral indices across the 2-GHz VLA bandwidth, finding a median index of -0.76$\pm$0.04. Stacking of source spectra reveals a flattening of spectral index with decreasing flux density. We present a source count derived from the catalogue. We show a traditional count estimate compared with a completely independent estimate made via a P(D) confusion analysis, and find very good agreement. Cross-matches of the catalogue with X-ray, optical, infrared, radio, and redshift catalogues are also presented. The X-ray, optical and infrared data, as well as AGN selection criteria allow us to classify 10 per cent as radio-loud AGN, 28 percent as radio-quiet AGN, and 58 per cent as star-forming galaxies, with only 4 per cent unclassified.

This is the first of two papers describing the observations and cataloguing of deep 3-GHz observations of the Lockman Hole North using the Karl G. Jansky Very Large Array. The aim of this paper is to investigate, through the use of simulated images, the uncertainties and accuracy of source-finding routines, as well as to quantify systematic effects due to resolution, such as source confusion and source size. While these effects are not new, this work is intended as a particular case study that can be scaled and translated to other surveys. We use the simulations to derive uncertainties in the fitted parameters, as well as bias corrections for the actual catalogue (presented in Paper 2). We compare two different source-finding routines, OBIT and AEGEAN, and two different effective resolutions, 8 and 2.75 arcsec. We find that the two routines perform comparably well, with OBIT being slightly better at de-blending sources, but slightly worse at fitting resolved sources. We show that 30 to 70 per cent of sources are missed or fit inaccurately once the source size becomes larger than the beam, possibly explaining source count errors in high-resolution surveys. We also investigate the effect of blending, finding that any sources with separations smaller than the beam size are fit as single sources. We show that the use of machine-learning techniques can correctly identify blended sources up to 90 per cent of the time, and prior-driven fitting can lead to a 70 per cent improvement in the number of de-blended sources.

Astronomy at or below the 'survey threshold' has expanded significantly since the publication of the original 'Science with the Square Kilometer Array' in 1999 and its update in 2004. The techniques in this regime may be broadly (but far from exclusively) defined as 'confusion' or 'P(D)' analyses (analyses of one-point statistics), and 'stacking', accounting for the flux-density distribution of noise-limited images co-added at the positions of objects detected/isolated in a different waveband. Here we discuss the relevant issues, present some examples of recent analyses, and consider some of the consequences for the design and use of surveys with the SKA and its pathfinders.

We present a study of diffuse extragalactic radio emission at $1.75\,$GHz from part of the ELAIS-S1 field using the Australia Telescope Compact Array. The resulting mosaic is $2.46\,$deg$^2$, with a roughly constant noise region of $0.61\,$deg$^2$ used for analysis. The image has a beam size of $150 \times60\,$arcsec and instrumental $\langle\sigma_{\rm n}\rangle= (52\pm5)\, \mu$Jy beam$^{-1}$. Using point-source models from the ATLAS survey, we subtract the discrete emission in this field for $S \ge 150\, \mu$Jy beam$^{-1}$. Comparison of the source-subtracted probability distribution, or \pd, with the predicted distribution from unsubtracted discrete emission and noise, yields an excess of $(76 \pm 23) \, \mu$Jy beam$^{-1}$. Taking this as an upper limit on any extended emission we constrain several models of extended source counts, assuming $\Omega_{\rm source} \le 2\,$arcmin. The best-fitting models yield temperatures of the radio background from extended emission of $T_{\rm b}=(10\pm7) \,$mK, giving an upper limit on the total temperature at $1.75\,$GHz of $(73\pm10)\,$mK. Further modelling shows that our data are inconsistent with the reported excess temperature of ARCADE2 to a source-count limit of $1\, \mu$Jy. Our new data close a loop-hole in the previous constraints, because of the possibility of extended emission being resolved out at higher resolution. Additionally, we look at a model of cluster halo emission and two WIMP dark matter annihilation source-count models, and discuss general constraints on any predicted counts from such sources. Finally, we report the derived integral count at $1.4\,$GHz using the deepest discrete count plus our new extended-emission limits, providing numbers that can be used for planning future ultra-deep surveys.

Radio source counts constrain galaxy populations and evolution, as well as the global star formation history. However, there is considerable disagreement among the published 1.4-GHz source counts below 100 microJy. Here we present a statistical method for estimating the microJy and even sub-microJy source count using new deep wide-band 3-GHz data in the Lockman Hole from the Karl G. Jansky Very Large Array (VLA). We analyzed the confusion amplitude distribution P(D), which provides a fresh approach in the form of a more robust model, with a comprehensive error analysis. We tested this method on a large-scale simulation, incorporating clustering and finite source sizes. We discuss in detail our statistical methods for fitting using Monte Carlo Markov chains, handling correlations, and systematic errors from the use of wide-band radio interferometric data. We demonstrated that the source count can be constrained down to 50 nJy, a factor of 20 below the rms confusion. We found the differential source count near 10 microJy to have a slope of -1.7, decreasing to about -1.4 at fainter flux densities. At 3GHz the rms confusion in an 8arcsec FWHM beam is ~ 1.2 microJy/beam, and a radio background temperature ~ 14mK. Our counts are broadly consistent with published evolutionary models. With these results we were also able to constrain the peak of the Euclidean normalized differential source count of any possible new radio populations that would contribute to the cosmic radio background down to 50 nJy.

We used the Karl G. Jansky Very Large Array (VLA) to image one primary beam area at 3 GHz with 8 arcsec FWHM resolution and 1.0 microJy/beam rms noise near the pointing center. The P(D) distribution from the central 10 arcmin of this confusion-limited image constrains the count of discrete sources in the 1 < S(microJy/beam) < 10 range. At this level the brightness-weighted differential count S^2 n(S) is converging rapidly, as predicted by evolutionary models in which the faintest radio sources are star-forming galaxies; and ~96$% of the background originating in galaxies has been resolved into discrete sources. About 63% of the radio background is produced by AGNs, and the remaining 37% comes from star-forming galaxies that obey the far-infrared (FIR) / radio correlation and account for most of the FIR background at lambda = 160 microns. Our new data confirm that radio sources powered by AGNs and star formation evolve at about the same rate, a result consistent with AGN feedback and the rough correlation of black hole and bulge stellar masses. The confusion at centimeter wavelengths is low enough that neither the planned SKA nor its pathfinder ASKAP EMU survey should be confusion limited, and the ultimate source detection limit imposed by "natural" confusion is < 0.01 microJy at 1.4 GHz. If discrete sources dominate the bright extragalactic background reported by ARCADE2 at 3.3 GHz, they cannot be located in or near galaxies and most are < 0.03 microJy at 1.4 GHz.

We examine the brightness of the Cosmic Radio Background (CRB) by comparing the contribution from individual source counts to absolute measurements. We use a compilation of radio counts to estimate the contribution of detected sources to the CRB in several different frequency bands.We apply a Monte Carlo Markov Chain technique to estimate the brightness values and uncertainties, paying attention to various sources of systematic error. We compare our results to absolute measurements from the ARCADE 2 experiment. At v = 150 MHz, 325 MHz, 408 MHz, 610 MHz, 1.4 GHz, 4.8 GHz, and 8.4 GHz our calculated contributions to the background sky temperature are 18, 2.8, 1.6, 0.71, 0.11, 0.0032, 0.0059 K, respectively. If the ARCADE 2 measurements are correct and come from sources, then there must be an additional population of radio galaxies, fainter than where current data are probing. More specifically, the Euclidean-normalized counts at 1.4 GHz have to have an additional bump below about 10 {\mu}Jy.