We have found an unexpected class of astronomical objects which have not previously been reported, in the Evolutionary Map of the Universe Pilot survey, using the Australian Square Kilometre Array Pathfinder telescope. The objects appear in radio images as circular edge-brightened discs about one arcmin diameter, and do not seem to correspond to any known type of object. We speculate that they may represent a spherical shock wave from an extra-galactic transient event, or the outflow, or a remnant, from a radio galaxy viewed end-on.

Modern radio telescopes are routinely reaching depths where normal starforming galaxies are the dominant observed population. Realising the potential of radio as a tracer of star formation and black hole activity over cosmic time involves achieving such depths over representative volumes, with radio forming part of a larger multiwavelength campaign. In pursuit of this we used the Karl G. Jansky Very Large Array (VLA) to image $\sim$5 deg$^{2}$ of the VIDEO/XMM-LSS extragalactic deep field at 1--2 GHz. We achieve a median depth of 16 $\mu$Jy beam$^{-1}$ with an angular resolution of 4.5\arcsec. Comparisons with existing radio observations of XMM-LSS showcase the improved survey speed of the upgraded VLA: we cover 2.5 times the area and increase the depth by $\sim$20\% in 40\% of the time. Direction-dependent calibration and wide-field imaging were required to suppress the error patterns from off-axis sources of even modest brightness. We derive a catalogue containing 5,762 sources from the final mosaic. Sub-band imaging provides in-band spectral indices for 3,458 (60\%) sources, with the average spectrum becoming flatter than the canonical synchrotron slope below 1 mJy. Positional and flux-density accuracy of the observations, and the differential source counts are in excellent agreement with those of existing measurements. A public release of the images and catalogue accompanies this article.

Tighter constraints on measurements of primordial non-Gaussianity will allow the differentiation of inflationary scenarios. The cosmic microwave background bispectrum-the standard method of measuring the local non-Gaussianity-is limited by cosmic variance. Therefore, it is sensible to investigate measurements of non-Gaussianity using the large-scale structure. This can be done by investigating the effects of non-Gaussianity on the power spectrum on large scales. In this study we forecast the constraints on the local primordial non-Gaussianity parameter $f_{\rm NL}$ that can be obtained with future radio surveys. We utilize the multi-tracer method which reduces the effect of cosmic variance and takes advantage of the multiple radio galaxy populations which are differently biased tracers of the same underlying dark matter distribution. Improvements on previous work include the use of observational bias and halo mass estimates, updated simulations and realistic photometric redshift expectations, thus producing more realistic forecasts. Combinations of SKA simulations and radio observations were used as well as different redshift ranges and redshift bin sizes. It was found that in the most realistic case the 1 - $\sigma$ error on $f_{\rm NL}$ falls within the range 4.07 and 6.58, rivalling the tightest constraints currently available.

The LOFAR Two-metre Sky Survey (LoTSS) will map the complete Northern sky and provide an excellent opportunity to study the distribution and evolution of the large-scale structure of the Universe. We study the completeness of the LoTSS first data release (DR1) and find a point-source completeness of 99 % above flux densities of 0.8 mJy and define a suite of quality cuts. We determine the count-in-cell statistics and differential source counts statistic and measure the angular two-point correlation function of the LoTSS radio sources. The counts-in-cell statistic reveals that the distribution of radio sources cannot be described by a spatial Poisson process. Instead, a good fit is provided by a compound Poisson distribution. The differential source counts are in good agreement with previous findings in deep fields at low radio frequencies. Simulated catalogues from the SKA design study sky match with our findings, in contrast to the more recent T-RECS models. For angular separations between 0.1 deg and 1 deg we find self-consistency amongst different estimates of the angular two-point correlation. The angular two-point correlation is $<10^{-2}$ at angular scales $> 1$ deg and up to the largest scales probed. We see indications of unidentified systematic issues at angular scales $>1$ deg, likely due to data calibration issues. Based on the distribution of photometric redshifts of LoTSS sources and the Planck 2018 best-fit cosmological model, the theoretically predicted angular two-point correlation between 0.1 deg and 1 deg agrees with the measured clustering for a subsample of radio sources with redshift information. We find agreement with the expectation of large-scale statistical isotropy of the radio sky at the per cent level and below angular separations of 1 deg the angular two-point correlation agrees with the expectation of the cosmological standard model. (abbreviated)

In the current era of radio astronomy, continuum surveys observe a multitude of objects with complex morphologies and sizes, and are not limited to observing point sources. Typical radio source extraction software generates catalogues by using Gaussian components to form a model of the emission. This may not be well suited to complicated jet structures and extended emission, particularly in the era of interferometers with a high density of short baselines, that are sensitive to extended emission. In this paper, we investigate how the optically-motivated source detection package ProFound (Robotham et al. 2018) may be used to model radio emission of both complicated and point-like radio sources. We use a combination of observations and simulations to investigate how ProFound compares to other source extractor packages used for radio surveys. We find that ProFound can accurately recover both the flux densities of simulated Gaussian sources as well as extended radio galaxies. ProFound can create models that trace the complicated nature of these extended galaxies, which we show is not necessarily the case with other source extraction software. Our work suggests that our knowledge of the emission from extended radio objects may be both over or under-estimated using traditional software. We suggest that ProFound offers a useful alternative to the fitting of Gaussian components for generating catalogues from current and future radio surveys. Furthermore, ProFound's multi-wavelength capabilities will be useful in investigating radio sources in combination with multiwavelength data.

We present observations of the XMM Large-Scale Structure (XMM-LSS) field observed with the LOw Frequency ARray (LOFAR) at 120-168 MHz. Centred at a J2000 declination of $-4.5^{\circ}$, this is a challenging field to observe with LOFAR because of its low elevation with respect to the array. The low elevation of this field reduces the effective collecting area of the telescope, thereby reducing sensitivity. This low elevation also causes the primary beam to be elongated in the north-south direction, which can introduce side lobes in the synthesised beam in this direction. However the XMM-LSS field is a key field to study because of the wealth of ancillary information, encompassing most of the electromagnetic spectrum. The field was observed for a total of 12 hours from three four-hour LOFAR tracks using the Dutch array. The final image presented encompasses $\sim 27$ deg$^2$, which is the region of the observations with a $>$50\% primary beam response. Once combined, the observations reach a central rms of 280 $\mu$Jy beam$^{-1}$ at 144 MHz and have an angular resolution of $7.5 \times \ 8.5$". We present our catalogue of detected sources and investigate how our observations compare to previous radio observations. This includes investigating the flux scale calibration of these observations compared to previous measurements, the implied spectral indices of the sources, the observed source counts and corrections to obtain the true source counts, and finally the clustering of the observed radio sources.

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.

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.

By studying the large-scale structure of the bright high-redshift Lyman-break galaxy (LBG) population it is possible to gain an insight into the role of environment in galaxy formation physics in the early Universe. We measure the clustering of a sample of bright (-22.7<M_UV<-21.125) LBGs at z~6 and use a halo occupation distribution (HOD) model to measure their typical halo masses. We find that the clustering amplitude and corresponding HOD fits suggests that these sources are highly biased (b~8) objects in the densest regions of the high-redshift Universe. Coupled with the observed rapid evolution of the number density of these objects, our results suggest that the shape of high luminosity end of the luminosity function is related to feedback processes or dust obscuration in the early Universe - as opposed to a scenario where these sources are predominantly rare instances of the much more numerous M_UV ~ -19 population of galaxies caught in a particularly vigorous period of star formation. There is a slight tension between the number densities and clustering measurements, which we interpret this as a signal that a refinement of the model halo bias relation at high redshifts or the incorporation of quasi-linear effects may be needed for future attempts at modelling the clustering and number counts. Finally, the difference in number density between the fields (UltraVISTA has a surface density ~1.8 times greater than UDS) is shown to be consistent with the cosmic variance implied by the clustering measurements.

Dark matter haloes in which galaxies reside are likely to have a significant impact on their evolution. We investigate the link between dark matter haloes and their constituent galaxies by measuring the angular two-point correlation function of radio sources, using recently released 3 GHz imaging over $\sim 2 \ \mathrm{deg}^2$ of the COSMOS field. We split the radio source population into Star Forming Galaxies (SFGs) and Active Galactic Nuclei (AGN), and further separate the AGN into radiatively efficient and inefficient accreters. Restricting our analysis to $z<1$, we find SFGs have a bias, $b = 1.5 ^{+0.1}_{-0.2}$, at a median redshift of $z=0.62$. On the other hand, AGN are significantly more strongly clustered with $b = 2.1\pm 0.2$ at a median redshift of 0.7. This supports the idea that AGN are hosted by more massive haloes than SFGs. We also find low-accretion rate AGN are more clustered ($b = 2.9 \pm 0.3$) than high-accretion rate AGN ($b = 1.8^{+0.4}_{-0.5}$) at the same redshift ($z \sim 0.7$), suggesting that low-accretion rate AGN reside in higher mass haloes. This supports previous evidence that the relatively hot gas that inhabits the most massive haloes is unable to be easily accreted by the central AGN, causing them to be inefficient. We also find evidence that low-accretion rate AGN appear to reside in halo masses of $M_{h} \sim 3-4 \times 10^{13}h^{-1}$M$_{\odot}$ at all redshifts. On the other hand, the efficient accreters reside in haloes of $M_{h} \sim 1-2 \times 10^{13}h^{-1}$M$_{\odot}$ at low redshift but can reside in relatively lower mass haloes at higher redshifts. This could be due to the increased prevalence of cold gas in lower mass haloes at $z \ge 1$ compared to $z<1$.

The MIGHTEE large survey project will survey four of the most well-studied extragalactic deep fields, totalling 20 square degrees to $\mu$Jy sensitivity at Giga-Hertz frequencies, as well as an ultra-deep image of a single ~1 square degree MeerKAT pointing. The observations will provide radio continuum, spectral line and polarisation information. As such, MIGHTEE, along with the excellent multi-wavelength data already available in these deep fields, will allow a range of science to be achieved. Specifically, MIGHTEE is designed to significantly enhance our understanding of, (i) the evolution of AGN and star-formation activity over cosmic time, as a function of stellar mass and environment, free of dust obscuration; (ii) the evolution of neutral hydrogen in the Universe and how this neutral gas eventually turns into stars after moving through the molecular phase, and how efficiently this can fuel AGN activity; (iii) the properties of cosmic magnetic fields and how they evolve in clusters, filaments and galaxies. MIGHTEE will reach similar depth to the planned SKA all-sky survey, and thus will provide a pilot to the cosmology experiments that will be carried out by the SKA over a much larger survey volume.

In these proceedings we highlight the primary scientific goals and design of the WEAVE-LOFAR survey, which will use the new WEAVE spectrograph on the 4.2m William Herschel Telescope to provide the primary source of spectroscopic information for the LOFAR Surveys Key Science Project. Beginning in 2018, WEAVE-LOFAR will generate more than 10$^6$ R=5000 365-960 nm spectra of low-frequency selected radio sources, across three tiers designed to efficiently sample the redshift-luminosity plane, and produce a data set of enormous legacy value. The radio frequency selection, combined with the high multiplex and throughput of the WEAVE spectrograph, make obtaining redshifts in this way very efficient, and we expect that the redshift success rate will approach 100 per cent at $z < 1$. This unprecedented spectroscopic sample - which will be complemented by an integral field component - will be transformational in key areas, including studying the star formation history of the Universe, the role of accretion and AGN-driven feedback, properties of the epoch of reionisation, cosmology, cluster haloes and relics, as well as the nature of radio galaxies and protoclusters. Each topic will be addressed in unprecedented detail, and with the most reliable source classifications and redshift information in existence.