We describe system verification tests and early science results from the pulsar processor (PTUSE) developed for the newly-commissioned 64-dish SARAO MeerKAT radio telescope in South Africa. MeerKAT is a high-gain (~2.8 K/Jy) low-system temperature (~18 K at 20cm) radio array that currently operates from 580-1670 MHz and can produce tied-array beams suitable for pulsar observations. This paper presents results from the MeerTime Large Survey Project and commissioning tests with PTUSE. Highlights include observations of the double pulsar J0737-3039A, pulse profiles from 34 millisecond pulsars from a single 2.5h observation of the Globular cluster Terzan 5, the rotation measure of Ter5O, a 420-sigma giant pulse from the Large Magellanic Cloud pulsar PSR J0540-6919, and nulling identified in the slow pulsar PSR J0633-2015. One of the key design specifications for MeerKAT was absolute timing errors of less than 5 ns using their novel precise time system. Our timing of two bright millisecond pulsars confirm that MeerKAT delivers exceptional timing. PSR J2241-5236 exhibits a jitter limit of <4 ns per hour whilst timing of PSR J1909-3744 over almost 11 months yields an rms residual of 66 ns with only 4 min integrations. Our results confirm that the MeerKAT is an exceptional pulsar telescope. The array can be split into four separate sub-arrays to time over 1000 pulsars per day and the future deployment of S-band (1750-3500 MHz) receivers will further enhance its capabilities.

In this study we define a new sample of distant powerful radio galaxies in order to study their host-galaxy properties and provide targets for future observations of HI absorption with new radio telescopes and to understand the fuelling and feedback from such sources. We have cross-matched the Sydney University Molonglo Sky Survey (SUMSS) radio catalogue at 843 MHz with the VISTA Hemisphere Survey (VHS) near-infrared catalogue using the Likelihood Ratio technique. Photometric redshifts from the Dark Energy Survey are then used to assign redshifts to the radio source counterparts. We found a total of 249 radio sources with photometric redshifts over a 148 deg^2 region. By fitting the optical and near-infrared photometry with spectral synthesis models we determine the stellar mass and star-formation rates of the radio sources, finding typical stellar masses of 10^{11} - 10^{12}M$_{\odot}$ for the powerful high-redshift radio galaxies. We also find a population of low-mass blue galaxies. However, by comparing the derived star-formation rates to the radio luminosity, we suggest that these sources are false positives in our likelihood ratio analysis. Follow up, higher-resolution (<5 arcsec) radio imaging would help alleviate these mid-identifications, as the limiting factor in our cross-identifications is the low resolution (~45 arcsec) of the SUMSS radio imaging.

Using a Karl G. Jansky Very Large Array (VLA) L-band (1-2 GHz) survey covering $\sim100$\,deg$^2$ of the Stripe 82 field, we have obtained a catalogue of 2716 radio AGN. For these AGN, we investigate the impact of galaxy density on 1.4 GHz radio luminosity ($L_{1.4}$). We determine their close environment densities using the surface density parameter, $\Sigma_N,$ for $N=2$ and $N=5,$ which we bin by redshift to obtain a pseudo-3D galaxy density measure. Matching the radio AGN to sources without radio detections in terms of redshift, $K-$band magnitude and ($g-K$)-colour index, we obtain samples of control galaxies and determine whether radio AGN environments differ from this general population. Our results indicate that the environmental density of radio AGN and their radio luminosity are not correlated up to $z$ $\sim0.8$, over the luminosity range $10^{23} < (L_{1.4} / $W~Hz$^{-1}) < 10^{26}$. We also find that, when using a control sample matched in terms of redshift, $K-$band magnitude and colour, environments of radio AGN are similar to those of the control sample but with an excess of overdense regions in which radio AGN are more prevalent. Our results suggest that the $<1$ Mpc-scale galaxy environment plays some role in determining whether a galaxy produces a radio AGN. The jet power, however, does not correlate with environment. From this, we infer that secular processes e.g. accretion flows of cold gas to the central black-hole are more critical in fuelling radio AGN activity than radio jet power.

A sample of 1161 radio galaxies with $0.01 < z < 0.7$ and $10^{21} < L_{1.4~\rm GHz} / \textrm{W Hz}^{-1} < 10^{27}$ is selected from the Stripe 82 1-2 GHz Karl G. Jansky Very Large Array Snapshot Survey, which covers 100 square degrees and has a 1$\sigma$ noise level of 88 $\mu$Jy beam$^{-1}$. Optical spectra are used to classify these sources as high excitation and low excitation radio galaxies (HERGs and LERGs), resulting in 60 HERGs, 149 LERGs and 600 `probable LERGs'. The host galaxies of the LERGs have older stellar populations than those of the HERGs, in agreement with previous results in the literature. We find that the HERGs tend to have higher Eddington-scaled accretion rates than the LERGs but that there is some overlap between the two distributions. We show that the properties of the host galaxies vary continuously with accretion rate, with the most slowly accreting sources having the oldest stellar populations, consistent with the idea that these sources lack a supply of cold gas. We find that 84 per cent of our sample release more than 10 per cent of their accretion power in their jets, showing that mechanical AGN feedback is significantly underestimated in many hydrodynamical simulations. There is a scatter of $\sim 2$ dex in the fraction of the accreted AGN power deposited back into the ISM in mechanical form, showing that the assumption in many simulations that there is a direct scaling between accretion rate and radio-mode feedback does not necessarily hold. We also find that mechanical feedback is significant for many of the HERGs in our sample as well as the LERGs.

We have combined spectrosopic and photometric data from the Sloan Digital Sky Survey (SDSS) with $1.4$ GHz radio observations, conducted as part of the Stripe 82 $1-2$ GHz Snapshot Survey using the Karl G. Jansky Very Large Array (VLA), which covers $\sim100$ sq degrees, to a flux limit of 88 $\mu$Jy rms. Cross-matching the $11\,768$ radio source components with optical data via visual inspection results in a final sample of $4\,795$ cross-matched objects, of which $1\,996$ have spectroscopic redshifts and $2\,799$ objects have photometric redshifts. Three previously undiscovered Giant Radio Galaxies (GRGs) were found during the cross-matching process, which would have been missed using automated techniques. For the objects with spectroscopy we separate radio-loud Active Galactic Nuclei (AGN) and star-forming galaxies (SFGs) using three diagnostics and then further divide our radio-loud AGN into the HERG and LERG populations. A control matched sample of HERGs and LERGs, matched on stellar mass, redshift and radio luminosity, reveals that the host galaxies of LERGs are redder and more concentrated than HERGs. By combining with near-infrared data, we demonstrate that LERGs also follow a tight $K-z$ relationship. These results imply the LERG population are hosted by population of massive, passively evolving early-type galaxies. We go on to show that HERGs, LERGs, QSOs and star-forming galaxies in our sample all reside in different regions of a WISE colour-colour diagram. This cross-matched sample bridges the gap between previous `wide but shallow' and `deep but narrow' samples and will be useful for a number of future investigations.

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.

We present a catalogue of nearly 3,000 submillimetre sources detected at 850um over ~5 square degrees surveyed as part of the James Clerk Maxwell Telescope (JCMT) SCUBA-2 Cosmology Legacy Survey (S2CLS). This is the largest survey of its kind at 850um, probing a meaningful cosmic volume at the peak of star formation activity and increasing the sample size of submillimetre galaxies selected at 850um by an order of magnitude. We describe the wide 850um survey component of S2CLS, which covers the key extragalactic survey fields: UKIDSS-UDS, COSMOS, Akari-NEP, Extended Groth Strip, Lockman Hole North, SSA22 and GOODS-North. The average 1-sigma depth of S2CLS is 1.2 mJy/beam, approaching the SCUBA-2 850um confusion limit, which we determine to be ~0.8 mJy/beam. We measure the single dish 850um number counts to unprecedented accuracy, reducing the Poisson errors on the differential counts to approximately 4% at S_850~3mJy. With several independent fields, we investigate field-to-field variance, finding that the number counts on 0.5-1 degree scales are generally within 50% of the S2CLS mean for S_850>3mJy, with scatter consistent with the Poisson and estimated cosmic variance uncertainties, although there is a marginal (2-sigma) density enhancement in the GOODS-North field. The observed number counts are in reasonable agreement with recent phenomenological and semi-analytic models. Finally, the large solid angle of S2CLS allows us to measure the bright-end counts: at S_850>10mJy there are approximately ten sources per square degree, and we detect the distinctive up-turn in the number counts indicative of the detection of local sources of 850um emission and strongly lensed high-redshift galaxies. Here we describe the data collection and reduction procedures and present calibrated maps and a catalogue of sources; these are made publicly available.

Measurement of the evolution of both active galactic nuclei (AGN) and star-formation in galaxies underpins our understanding of galaxy evolution over cosmic time. Radio continuum observations can provide key information on these two processes, in particular via the mechanical feedback produced by radio jets in AGN, and via an unbiased dust-independent measurement of star-formation rates. In this paper we determine radio luminosity functions at 325 MHz for a sample of AGN and star-forming galaxies by matching a 138 deg sq. radio survey conducted with the Giant Metrewave Radio Telescope (GMRT), with optical imaging and redshifts from the Galaxy And Mass Assembly (GAMA) survey. We find that the radio luminosity function at 325 MHz for star-forming galaxies closely follows that measured at 1.4 GHz. By fitting the AGN radio luminosity function out to $z = 0.5$ as a double power law, and parametrizing the evolution as ${\Phi} \propto (1 + z)^{k}$ , we find evolution parameters of $k = 0.92 \pm 0.95$ assuming pure density evolution and $k = 2.13 \pm 1.96$ assuming pure luminosity evolution. We find that the Low Excitation Radio Galaxies are the dominant population in space density at lower luminosities. Comparing our 325 MHz observations with radio continuum imaging at 1.4 GHz, we determine separate radio luminosity functions for steep and flat-spectrum AGN, and show that the beamed population of flat-spectrum sources in our sample can be shifted in number density and luminosity to coincide with the unbeamed population of steep-spectrum sources, as is expected in the orientation based unification of AGN.

It has become apparent that active galactic nuclei (AGN) may have a significant impact on the growth and evolution of their host galaxies and vice versa but a detailed understanding of the interplay between these processes remains elusive. Deep radio surveys provide a powerful, obscuration-independent tool for measuring both star formation and AGN activity in high-redshift galaxies. Multiwavelength studies of deep radio fields show a composite population of star-forming galaxies and AGN, with the former dominating at the lowest flux densities (S$_{1.4\mathrm{GHz}}<$100~$\mu$Jy). The sensitivity and resolution of the SKA will allow us to identify, and separately trace, the total star formation in the bulges of individual high-redshift galaxies, the related nuclear activity and any star formation occurring on larger scales within a disc. We will therefore gain a detailed picture of the apparently simultaneous development of stellar populations and black holes in the redshift range where both star-formation and AGN activity peak (1$\leq$z$\leq$4). In this chapter we discuss the role of the SKA in studying the connection between AGN activity and galaxy evolution, and the most critical technical requirements for such of studies

Ultra Steep Spectrum (USS) radio sources are one of the efficient tracers of powerful High-z Radio Galaxies (HzRGs). In contrast to searches for powerful HzRGs from radio surveys of moderate depths, fainter USS samples derived from deeper radio surveys can be useful in finding HzRGs at even higher redshifts and in unveiling a population of obscured weaker radio-loud AGN at moderate redshifts. Using our 325 MHz GMRT observations (5-sigma ~ 800 microJy) and 1.4 GHz VLA observations (5-sigma ~ 80 - 100 microJy) available in two subfields (viz., VLA-VIMOS VLT Deep Survey (VLA-VVDS) and Subaru X-ray Deep Field (SXDF)) of the XMM-LSS field, we derive a large sample of 160 faint USS radio sources and characterize their nature. The optical, IR counterparts of our USS sample sources are searched using existing deep surveys, at respective wavelengths. We attempt to unveil the nature of our faint USS sources using diagnostic techniques based on mid-IR colors, flux ratios of radio to mid-IR, and radio luminosities. Redshift estimates are available for 86/116 (~ 74%) USS sources in the VLA-VVDS field and for 39/44 (~ 87%) USS sources in the SXDF fields with median values (z_median) ~ 1.18 and ~ 1.57, which are higher than that for non-USS radio sources (z_median non-USS ~ 0.99 and ~ 0.96), in the two subfields, respectively. The flux ratio of radio to mid-IR (S_1.4 GHz/S_3.6 micron) versus redshift diagnostic plot suggests that more than half of our USS sample sources distributed over z ~ 0.5 to 3.8 are likely to be hosted in obscured environments. A significant fraction (~ 26% in the VLA-VVDS and ~ 13% in the SXDF) of our USS sources without redshift estimates mostly remain unidentified in the existing optical, IR surveys, and exhibit high radio to mid-IR flux ratio limits similar to HzRGs, and thus, can be considered as potential HzRG candidates.

Quantifying how the baryonic matter traces the underlying dark matter distribution is key to both understanding galaxy formation and our ability to constrain the cosmological model. Using the cross-correlation function of radio and near-infrared galaxies, we present a large-scale clustering analysis of radio galaxies to z ~ 2.2. We measure the angular auto-correlation function of Ks < 23.5 galaxies in the VIDEO-XMM3 field with photometric redshifts out to z = 4 using VIDEO and CFHTLS photometry in the near-infrared and optical. We then use the cross-correlation function of these sources with 766 radio sources at S_1.4 > 90 {\mu}Jy to infer linear bias of radio galaxies in four redshift bins. We find that the bias evolves from b = 0.57 +/- 0.06 at z ~ 0.3 to 8.55 +/- 3.11 at z ~ 2.2. Furthermore, we separate the radio sources into subsamples to determine how the bias is dependent on the radio luminosity, and find a bias which is significantly higher than predicted by the simulations of Wilman et al., and consistent with the lower luminosity but more abundant FR-I population having a similar bias to the highly luminous but rare FR-IIs. Our results are suggestive of a higher mass, particularly for FR-I sources than assumed in simulations, especially towards higher redshift.

Ultra Steep Spectrum (USS) radio sources are one of the efficient tracers of High Redshift Radio Galaxies (HzRGs). To search for HzRGs candidates, we investigate properties of a large sample of faint USS sources derived from our deep 325 MHz GMRT observations combined with 1.4 GHz VLA data on the two subfields (i.e., VLA-VIMOS VLT Deep Survey (VVDS) and Subaru X-ray Deep Field (SXDF)) in the XMM-LSS field. The available redshift estimates show that majority of our USS sample sources are at higher redshifts with the median redshifts ~ 1.18 and ~ 1.57 in the VLA-VVDS and SXDF fields. In the VLA-VVDS field, ~ 20% of USS sources lack the redshift estimates as well as the detection in the deep optical, IR surveys, and thus these sources may be considered as potential high-z candidates. The radio luminosity distributions suggest that a substantial fraction (~ 40%) of our USS sample sources are radio-loud sources, distributed over redshifts ~ 0.5 to 4.

It has been speculated that low luminosity radio-loud AGN have the potential to serve as an important source of AGN feedback, and may be responsible for suppressing star-formation activity in massive elliptical galaxies at late times. As such the cosmic evolution of these sources is vitally important to understanding the significance of such AGN feedback processes and their influence on the global star-formation history of the universe. In this paper we present a new investigation of the evolution of faint radio sources out to $z{\sim}2.5$. We combine a 1 square degree VLA radio survey, complete to a depth of 100 $\mu$Jy, with accurate 10 band photometric redshifts from the VIDEO and CFHTLS surveys. The results indicate that the radio population experiences mild positive evolution out to $z{\sim}1.2$ increasing their space density by a factor of $\sim$3, consistent with results of several previous studies. Beyond $z$=1.2 there is evidence of a slowing down of this evolution. Star-forming galaxies drive the more rapid evolution at low redshifts, $z{<}$1.2, while more slowly evolving AGN populations dominate at higher redshifts resulting in a decline in the evolution of the radio luminosity function at $z{>}$1.2. The evolution is best fit by pure luminosity evolution with star-forming galaxies evolving as $(1+z)^{2.47\pm0.12}$ and AGN as $(1+z)^{1.18\pm0.21}$.

In the lead-up to the Square Kilometre Array (SKA) project, several next-generation radio telescopes and upgrades are already being built around the world. These include APERTIF (The Netherlands), ASKAP (Australia), eMERLIN (UK), VLA (USA), e-EVN (based in Europe), LOFAR (The Netherlands), Meerkat (South Africa), and the Murchison Widefield Array (MWA). Each of these new instruments has different strengths, and coordination of surveys between them can help maximise the science from each of them. A radio continuum survey is being planned on each of them with the primary science objective of understanding the formation and evolution of galaxies over cosmic time, and the cosmological parameters and large-scale structures which drive it. In pursuit of this objective, the different teams are developing a variety of new techniques, and refining existing ones. Here we describe these projects, their science goals, and the technical challenges which are being addressed to maximise the science return.

In this paper we describe the first data release of the the Visible and Infrared Survey Telescope for Astronomy (VISTA) Deep Extragalactic Observations (VIDEO) survey. VIDEO is a ~12degree^2 survey in the near-infrared Z,Y,J,H and K_s bands, specifically designed to enable the evolution of galaxies and large structures to be traced as a function of both epoch and environment from the present day out to z=4, and active galactic nuclei (AGN) and the most massive galaxies up to and into the epoch of reionization. With its depth and area, VIDEO will be able to fully explore the period in the Universe where AGN and starburst activity were at their peak and the first galaxy clusters were beginning to virialize. VIDEO therefore offers a unique data set with which to investigate the interplay between AGN, starbursts and environment, and the role of feedback at a time when it was potentially most crucial. We provide data over the VIDEO-XMM3 tile, which also covers the Canada-France-Hawaii-Telescope Legacy Survey Deep-1 field (CFHTLS-D1). The released VIDEO data reach a 5-sigma AB-magnitude depth of Z=25.7, Y=24.5, J=24.4, H=24.1 and K_s=23.8 in 2 arcsec diameter apertures (the full depth of Y=24.6 will be reached within the full integration time in future releases). The data are compared to previous surveys over this field and we find good astrometric agreement with the Two-Micron All Sky Survey, and source counts in agreement with the recently released UltraVISTA survey data. The addition of the VIDEO data to the CFHTLS-D1 optical data increases the accuracy of photometric redshifts and significantly reduces the fraction of catastrophic outliers over the redshift range 0<z<1 from 5.8 to 3.1 per cent in the absence of an i-band luminosity prior. (Truncated Abstract)

In this paper, we explore how the forthcoming generation of large-scale radio continuum surveys, with the inclusion of some degree of redshift information, can constrain cosmological parameters. By cross-matching these radio surveys with shallow optical to near-infrared surveys, we can essentially separate the source distribution into a low- and a high-redshift sample, thus providing a constraint on the evolution of cosmological parameters such as those related to dark energy. We examine two radio surveys, the Evolutionary Map of the Universe (EMU) and the Westerbork Observations of the Deep APERTIF Northern sky (WODAN). A crucial advantage is their combined potential to provide a deep, full-sky survey. The surveys used for the cross-identifications are SkyMapper and SDSS, for the southern and northern skies, respectively. We concentrate on the galaxy clustering angular power spectrum as our benchmark observable, and find that the possibility of including such low redshift information yields major improvements in the determination of cosmological parameters. With this approach, and provided a good knowledge of the galaxy bias evolution, we are able to put strict constraints on the dark energy parameters, i.e. w_0=-0.9+/-0.041 and w_a=-0.24+/-0.13, with type Ia supernovae and CMB priors (with a one-parameter bias in this case); this corresponds to a Figure of Merit (FoM) > 600, which is twice better than what is obtained by using only the cross-identified sources and greater than four time better than the case without any redshift information at all.

In this paper we investigate the performance of the likelihood ratio method as a tool for identifying optical and infrared counterparts to proposed radio continuum surveys with SKA precursor and pathfinder telescopes. We present a comparison of the infrared counterparts identified by the likelihood ratio in the VISTA Deep Extragalactic Observations (VIDEO) survey to radio observations with 6, 10 and 15 arcsec resolution. We cross-match a deep radio catalogue consisting of radio sources with peak flux density $>$ 60 $\mu$Jy with deep near-infrared data limited to $K_{\mathrm{s}}\lesssim$ 22.6. Comparing the infrared counterparts from this procedure to those obtained when cross-matching a set of simulated lower resolution radio catalogues indicates that degrading the resolution from 6 arcsec to 10 and 15 arcsec decreases the completeness of the cross-matched catalogue by approximately 3 and 7 percent respectively. When matching against shallower infrared data, comparable to that achieved by the VISTA Hemisphere Survey, the fraction of radio sources with reliably identified counterparts drops from $\sim$89%, at $K_{\mathrm{s}}\lesssim$22.6, to 47% with $K_{\mathrm{s}}\lesssim$20.0. Decreasing the resolution at this shallower infrared limit does not result in any further decrease in the completeness produced by the likelihood ratio matching procedure. However, we note that radio continuum surveys with the MeerKAT and eventually the SKA, will require long baselines in order to ensure that the resulting maps are not limited by instrumental confusion noise.

We investigate the cosmic evolution of low luminosity ($L_{\rm{1.4GHz}}<10^{25}\rm{W~Hz^{-1}sr^{-1}}$) radio sources in the XMM Large Scale Structure survey field (XMM-LSS). We match low frequency selected (610~MHz) radio sources in the XMM-LSS field with near infrared $K$-band observations over the same field from the UKIRT Infrared Deep Sky Survey. We use both the mean $V/V_{\rm{max}}$ statistic and the radio luminosity function of these matched sources to quantify the evolution of the co-moving space density of the low luminosity radio sources in our sample. Our results indicate that the low luminosity sources evolve differently to their high luminosity counterparts out to a redshift of z$\sim$0.8. The derived luminosity function is consistent with an increase in the co-moving space density of low luminosity sources by a factor of $\sim$1.5 at z=0.8. We show that the use of the $K-z$ diagram for the radio source population, although coarser than a full photometric redshift analysis, produces consistent results with previous studies using $\sim >10$ band photometry. This offers a promising method for conducting similar analyses over the whole sky with future near- and mid-infrared surveys.