Understanding the formation of stellar clusters requires following the interplay between gas and newly formed stars accurately. We therefore couple the magnetohydrodynamics code FLASH to the N-body code ph4 and the stellar evolution code SeBa using the Astrophysical Multipurpose Software Environment (AMUSE) to model stellar dynamics, evolution, and collisional N-body dynamics and the formation of binary and higher-order multiple systems, while implementing stellar feedback in the form of radiation, stellar winds and supernovae in FLASH. We use this novel numerical method to simulate the formation and early evolution of open clusters of $\sim 1000$ stars formed from clouds with a mass range of $10^3$-$10^5$ M$_\odot$. Analyzing the effects of stellar feedback on the gas and stars of the natal cluster, we find that our clusters are resilient to disruption, even in the presence of intense feedback. This can even slightly increase the amount of dense, Jeans unstable gas by sweeping up shells; thus, a stellar wind strong enough to trap its own H II region shows modest triggering of star formation. Our clusters are born moderately mass segregated, an effect enhanced by feedback, and retained after the ejection of their natal gas, in agreement with observations.

This white paper submitted for the 2020 Canadian Long-Range Planning process (LRP2020) presents the prospects for Canada and the Square Kilometre Array (SKA) from 2020-2030, focussing on the first phase of the project (SKA1) scheduled to begin construction early in the next decade. SKA1 will make transformational advances in our understanding of the Universe across a wide range of fields, and Canadians are poised to play leadership roles in several. Canadian key SKA technologies will ensure a good return on capital investment in addition to strong scientific returns, positioning Canadian astronomy for future opportunities well beyond 2030. We therefore advocate for Canada's continued scientific and technological engagement in the SKA from 2020-2030 through participation in the construction and operations phases of SKA1.

We describe a star cluster formation model that includes individual star formation from self-gravitating, magnetized gas, coupled to collisional stellar dynamics. The model uses the Astrophysical Multi-purpose Software Environment (AMUSE) to integrate an adaptive-mesh magnetohydrodynamics code (FLASH) with a fourth order Hermite N-body code (ph4), a stellar evolution code (SeBa), and a method for resolving binary evolution (multiples). This combination yields unique star formation simulations that allow us to study binaries formed dynamically from interactions with both other stars and dense, magnetized gas subject to stellar feedback during the birth and early evolution of stellar clusters. We find that for massive stars, our simulations are consistent with the observed dynamical binary fractions and mass ratios. However, our binary fraction drops well below observed values for lower mass stars, presumably due to unincluded binary formation during initial star formation. Further, we observe a build up of binaries near the hard-soft boundary that may be an important mechanism driving early cluster contraction.

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.

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.

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.

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.

The new broad-band capabilities of large radio interferometers such as the GMRT and JVLA allow for long-integration mosaic imaging observations to create ultra-deep full-polarization images of the sky over wide frequency ranges. Achieving rms sensitivities of order 1 $\mu$Jy, these observations explore the radio source population at flux densities well below the regime dominated by classical radio galaxies and Active Galactic Nuclei. We present initial results from radio sources revealed with deep mosaicking observations with the GMRT and JVLA at respectively 0.6 and 5 GHz, and evidence that the $\mu$Jy sensitivity level marks the transition to detection of polarized emission from a population of sources dominated by emission from magnetic fields in the disks of starburst and normal galaxies.

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.

The goal of this work is to determine the nature of the relation between morphology and accretion mode in radio galaxies, including environmental parameters. The CoNFIG extended catalogue (improved by new Ks-band identifications and estimated redshifts from UKIDSS, and spectral index measurements from new GMRT observations) is used to select a sub-sample of 206 radio galaxies with z<0.3 over a wide range of radio luminosity, which are morphology-classified using the Fanaroff-Riley (FR) classification of extended radio sources. For each galaxy, spectroscopic data are retrieved to determine the high/low excitation status of the source, related to its accretion mode. Environmental factors, such as the host galaxy luminosity and a richness factor are also computed, generally using SDSS data.

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.

EMU is a wide-field radio continuum survey planned for the new Australian Square Kilometre Array Pathfinder (ASKAP) telescope. The primary goal of EMU is to make a deep (rms ~ 10 microJy/beam) radio continuum survey of the entire Southern Sky at 1.3 GHz, extending as far North as +30 degrees declination, with a resolution of 10 arcsec. EMU is expected to detect and catalogue about 70 million galaxies, including typical star-forming galaxies up to z~1, powerful starbursts to even greater redshifts, and AGNs to the edge of the visible Universe. It will undoubtedly discover new classes of object. This paper defines the science goals and parameters of the survey, and describes the development of techniques necessary to maximise the science return from EMU.

This paper presents a new grid-based method for investigating the evolution of the steep-spectrum radio luminosity function, with the aim of quantifying the high-redshift cut-off suggested by previous work. To achieve this, the Combined EIS-NVSS Survey of Radio Sources (CENSORS) has been developed; this is a 1.4 GHz radio survey, containing 135 sources complete to a flux density of 7.2 mJy, selected from the NRAO VLA Sky Survey (NVSS) over 6 deg^2 of the ESO Imaging Survey (EIS) Patch D. The sample is currently 73% spectroscopically complete, with the remaining redshifts estimated via the K-z or I-z magnitude-redshift relation. CENSORS is combined with additional radio data from the Parkes All-Sky, Parkes Selected Regions, Hercules and VLA COSMOS samples to provide comprehensive coverage of the radio power vs. redshift plane. The redshift distributions of these samples, together with radio source count determinations, and measurements of the local luminosity function, provide the input to the fitting process. The modelling reveals clear declines, at > 3sigma significance, in comoving density at z > 0.7 for lower luminosity sources (log P = 25-26); these turnovers are still present at log P > 27, but move to z > 3, suggesting a luminosity-dependent evolution of the redshift turnover, similar to the `cosmic downsizing' seen for other AGN populations. These results are shown to be robust to the estimated redshift errors and to increases in the spectral index for the highest redshift sources. Analytic fits to the best-fitting steep spectrum grid are provided so that the results presented here can be easily accessed by the reader, as well as allowing plausible extrapolations outside of the regions covered by the input datasets

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.

We describe a survey in the ELAIS N2 region with the VLA at 43.4 GHz, carried out with 1627 independent snapshot observations in D-configuration and covering about 0.5 square degrees. One certain source is detected, a previously-catalogued flat-spectrum QSO at z=2.2. A few (<5) other sources may be present at about the 3sigma level, as determined from positions of source-like deflections coinciding with blue stellar objects, or with sources from lower-frequency surveys. Independently we show how all the source-like detections identified in the data can be used with a maximum-likelihood technique to constrain the 43-GHz source counts at a level of ~7 mJy. Previous estimates of the counts at 43 GHz, based on lower-frequency counts and spectral measurements, are consistent with these constraints, although the present results are suggestive of somewhat higher surface densities at the 7 mJy level. They do not provide direct evidence of intrusion of a previously unknown source population, although the several candidate sources need examination before such a population can be ruled out.

This paper focuses on a comparison of the space densities of FRI and FRII sources at different epochs, with a particular focus on FRI sources. First, we present the concluding steps in constructing the Combined NVSS-FIRST Galaxy catalogue (CoNFIG), including new VLA observations, optical identifications and redshift estimates. The final catalogue consists of 859 sources over 4 samples (CoNFIG-1, 2, 3 and 4 with flux density limits of S_1.4GHz = 1.3, 0.8, 0.2 and 0.05 Jy respectively). It is 95.7% complete in radio morphology classification and 74.3% of the sources have redshift data. Combining CoNFIG with complementary samples, the distribution and evolution of FRI and FRII sources are investigated. We find that FRI sources undergo mild evolution and that, at the same radio luminosity, FRI and FRII sources show similar space density enhancements in various redshift ranges, possibly implying a common evolution.

We focus on a comparison of the space densities of FRI and FRII extended radio sources at different epochs, and find that FRI and FRII sources show similar space density enhancements in various redshift ranges, possibly implying a common evolution.

We review the statistical properties of the main populations of radio sources, as emerging from radio and millimeter sky surveys. Recent determinations of local luminosity functions are presented and compared with earlier estimates still in widespread use. A number of unresolved issues are discussed. These include: the (possibly luminosity-dependent) decline of source space densities at high redshifts; the possible dichotomies between evolutionary properties of low- versus high-luminosity and of flat- versus steep-spectrum AGN-powered radio sources; and the nature of sources accounting for the upturn of source counts at sub-mJy levels. It is shown that straightforward extrapolations of evolutionary models, accounting for both the far-IR counts and redshift distributions of star-forming galaxies, match the radio source counts at flux-density levels of tens of microJy remarkably well. We consider the statistical properties of rare but physically very interesting classes of sources, such as GHz Peak Spectrum and ADAF/ADIOS sources, and radio afterglows of gamma-ray bursts. We also discuss the exploitation of large-area radio surveys to investigate large scale structure through studies of clustering and the Integrated Sachs-Wolfe effect. Finally we briefly describe the potential of the new and forthcoming generations of radio telescopes. A compendium of source counts at different frequencies is given in an appendix.

We use 3 years of data from the Swift/BAT survey to select a complete sample of X-ray blazars above 15 keV. This sample comprises 26 Flat-Spectrum Radio Quasars (FSRQs) and 12 BL Lac objects detected over a redshift range of 0.03<z<4.0. We use this sample to determine, for the first time in the 15--55 keV band, the evolution of blazars. We find that, contrary to the Seyfert-like AGNs detected by BAT, the population of blazars shows strong positive evolution. This evolution is comparable to the evolution of luminous optical QSOs and luminous X-ray selected AGNs. We also find evidence for an epoch-dependence of the evolution as determined previously for radio-quiet AGNs. We interpret both these findings as a strong link between accretion and jet activity. In our sample, the FSRQs evolve strongly, while our best-fit shows that BL Lacs might not evolve at all. The blazar population accounts for 10-20 % (depending on the evolution of the BL Lacs) of the Cosmic X-ray background (CXB) in the 15--55 keV band. We find that FSRQs can explain the entire CXB emission for energies above 500 keV solving the mystery of the generation of the MeV background. The evolution of luminous FSRQs shows a peak in redshift ($z_c$=4.3$\pm0.5$) which is larger than the one observed in QSOs and X-ray selected AGNs. We argue that FSRQs can be used as tracers of massive elliptical galaxies in the early Universe.

[ABRIDGED VERSION] The future of cm and m-wave astronomy lies with the Square Kilometre Array (SKA), a telescope under development by a consortium of 17 countries. The SKA will be 50 times more sensitive than any existing radio facility. A majority of the key science for the SKA will be addressed through large-area imaging of the Universe at frequencies from 300 MHz to a few GHz. The Australian SKA Pathfinder (ASKAP) is aimed squarely in this frequency range, and achieves instantaneous wide-area imaging through the development and deployment of phase-array feed systems on parabolic reflectors. This large field-of-view makes ASKAP an unprecedented synoptic telescope poised to achieve substantial advances in SKA key science. The central core of ASKAP will be located at the Murchison Radio Observatory in inland Western Australia, one of the most radio-quiet locations on the Earth and one of the sites selected by the international community as a potential location for the SKA. Following an introductory description of ASKAP, this document contains 7 chapters describing specific science programmes for ASKAP. The combination of location, technological innovation and scientific program will ensure that ASKAP will be a world-leading radio astronomy facility, closely aligned with the scientific and technical direction of the SKA. A brief summary chapter emphasizes the point, and considers discovery space.

The CoNFIG (Combined NVSS-FIRST Galaxies) sample is a new sample of 274 bright radio sources at 1.4 GHz. It was defined by selecting all sources with S_1.4GHz > 1.3 Jy from the NRAO-VLA Sky Survey (NVSS) in the North field of the Faint Images of the Radio Sky at Twenty-cm (FIRST) survey. New radio observations obtained with the VLA for 31 of the sources are presented. The sample has complete FRI/FRII morphology identification; optical identifications and redshifts are available for 80% and 89% of the sample respectively, yielding a mean redshift of ~0.71. One of the goals of this survey is to get better definitions of luminosity distributions and source counts of FRI/FRII sources separately, in order to determine the evolution of the luminosity function for each type of source. We present a preliminary analysis, showing that these data are an important step towards examining various evolutionary schemes for these objects and to confirm or correct the dual population unified scheme for radio AGN. Improving our understanding of radio galaxy evolution will give better insight into the role of AGN feedback in galaxy formation.

A picture has emerged connecting QSOs with Sub-Millimetre Galaxies (SMGs) through an evolutionary sequence in which forming galaxies are initially FIR-luminous but X-ray weak, similar to known SMGs. As the black hole and spheroid grow with time, the central QSO becomes powerful enough to terminate star formation and eject much of the fuel supply. The unobscured QSO activity subsequently declines to leave a quiescent spheroidal galaxy. Here I describe parallel investigations of space density, one for a sample of radio-loud QSOs (RQSOs), and a second for SMGs. Each class shows both cosmic down-sizing and a redshift cutoff. The coincidence in apparent epoch of creation is marked; if it does not prove a causal connection, it is at least circumstantial evidence that the foregoing sequence is correct.

The Combined NVSS-FIRST Galaxy (CoNFIG) survey was defined by selecting all sources with S_1.4GHz > 1.3Jy from the NRAO VLA Sky Survey (NVSS) in the north field of the Faint Images of the Radio Sky at Twenty-cm (FIRST) survey. We carried out FRI/FRII morphology classification from NVSS and FIRST survey data; to complete this process, new 8GHz VLA observations for 31 sources were obtained at 0.24 arcsec resolution. Optical identifications and redshift information were compiled for about ~80% of the 270 radio sources in the sample, the mean redshift being ~0.6. A major goal of this sample is a better definition of the individual luminosity distributions and source counts for FRI and FRII sources, in order to determine accurately the evolution of the luminosity function for both types. Amongst the aims are the issues of whether the two populations are really distinct, whether physical evolution schemes permit one type to evolve into the other, whether the dual-population unified scheme for radio AGN remains viable, and the role radio AGN - star-formation feedback mechanisms.

We explore the epoch dependence of number density and star-formation rate for submillimetre galaxies (SMGs) found at 850 um. The study uses a sample of 38 SMG in the GOODS-N field, for which cross-waveband identifications have been obtained for 35/38 members together with redshift measurements or estimates. A maximum-likelihood analysis is employed, along with the `single-source-survey' technique. We find a diminution in both space density and star formation rate at z > 3, closely mimicking the redshift cut-offs found for QSOs selected in different wavebands. The diminution in redshift is particularly marked, at a significance level too small to measure. The data further suggest, at a significance level of about 0.001, that two separately-evolving populations may be present, with distinct luminosity functions. These results parallel the different evolutionary behaviours of LIRGs and ULIRGs, and represent another manifestation of `cosmic down-sizing', suggesting that differential evolution extends to the most extreme star-forming galaxies.

Astronomers traditionally find steep, open-ended functions difficult: the log N - log S curve for example, or the probability-area-radius relation in making cross-waveband identifications. The debate over the existence of a redshift cutoff for RQSOs (radio-loud QSOs; Jarvis and Rawlings 2000, Wall et al. 2005) is a case in point: an instance in which a second-order effect involving radio-variable QSOs, initially selected by survey from a population with a steep source count, leads to discrepant results from different analyses. A redshift cutoff is present in optically-selected QSOs of the SDSS and in X-ray selected QSOs of XMM, ROSAT and CHANDRA. A similar cutoff for radio-loud QSOs is important to establish, because if present, no known type of obscuration can be responsible. Such a space-density diminution then defines an epoch of creation for galaxy hosts of massive black holes powering AGN, a datum for galaxy formation. I discuss the resolution of the debate: the issue is highly relevant to modern surveys and their follow-up observations.

We analyze the Parkes quarter-Jansky flat-spectrum sample of QSOs in terms of space density, including the redshift distribution, the radio luminosity function, and the evidence for a redshift cutoff. With regard to the luminosity function, we note the strong evolution in space density from the present day to epochs corresponding to redshifts ~ 1. We draw attention to a selection effect due to spread in spectral shape that may have misled other investigators to consider the apparent similarities in shape of luminosity functions in different redshift shells as evidence for luminosity evolution. To examine the evolution at redshifts beyond 3, we develop a model-independent method based on the V_max test using each object to predict expectation densities beyond z=3. With this we show that a diminution in space density at z > 3 is present at a significance level >4 sigma. We identify a severe bias in such determinations from using flux-density measurements at epochs significantly later than that of the finding survey. The form of the diminution is estimated, and is shown to be very similar to that found for QSOs selected in X-ray and optical wavebands. The diminution is also compared with the current estimates of star-formation evolution, with less conclusive results. In summary we suggest that the reionization epoch is little influenced by powerful flat-spectrum QSOs, and that dust obscuration does not play a major role in our view of the QSO population selected at radio, optical or X-ray wavelengths.

The Oxford Dartmouth Thirty Degree Survey (ODTS) is a deep, wide, multi-band imaging survey designed to cover a total of 30 square degrees in BVRi'Z, with a subset of U and K band data, in four separate fields of 5-10 deg^2 centred at 00:18:24 +34:52, 09:09:45 +40:50, 13:40:00 +02:30 and 16:39:30 +45:24. Observations have been made using the Wide Field Camera on the 2.5-m Isaac Newton Telescope in La Palma to average limiting depths (5 sigma Vega, aperture magnitudes) of U=24.8, B=25.6, V=25.0, R=24.6, and i'=23.5, with observations taken in ideal conditions reaching the target depths of U=25.3, B=26.2, V=25.7, R=25.4, and i'=24.6. The INT Z band data was found to be severely effected by fringing and, consequently, is now being obtained at the MDM observatory in Arizona. A complementary K-band survey has also been carried out at MDM, reaching an average depth of K_{5\sigma}~18.5. At present, approximately 23 deg^2 of the ODTS have been observed, with 3.5 deg^2 of the K band survey completed. This paper details the survey goals, field selection, observation strategy and data reduction procedure, focusing on the photometric calibration and catalogue construction. Preliminary photometric redshifts have been obtained for a subsample of the objects with R <= 23. These results are presented alongside a brief description of the photometric redshift determination technique used. The median redshift of the survey is estimated to be z~0.7 from a combination of the ODTS photometric redshifts and comparison with the redshift distributions of other surveys. Finally, galaxy number counts for the ODTS are presented which are found to be in excellent agreement with previous studies.

A new sample of radio sources, with the designated name CENSORS (A Combined EIS-NVSS Survey Of Radio Sources), has been defined by combining the NRAO VLA Sky Survey (NVSS) at 1.4 GHz with the ESO Imaging Survey (EIS) Patch D, a 3 by 2 degree region of sky centred at 09 51 36.0, -21 00 00 (J2000). New radio observations of 199 NVSS radio sources with NVSS flux densities S(1.4GHz) > 7.8mJy are presented, and compared with the EIS I-band imaging observations which reach a depth of I~23; optical identifications are obtained for over two-thirds of the ~150 confirmed radio sources within the EIS field. The radio sources have a median linear size of 6 arcseconds, consistent with the trend for lower flux density radio sources to be less extended. Other radio source properties, such as the lobe flux density ratios, are consistent with those of brighter radio source samples. From the optical information, 30-40% of the sources are expected to lie at redshifts z >~ 1.5. One of the key goals of this survey is to accurately determine the high redshift evolution of the radio luminosity function. These radio sources are at the ideal flux density level to achieve this goal; at redshifts z~2 they have luminosities which are around the break of the luminosity function and so provide a much more accurate census of the radio source population at those redshifts than the existing studies of extreme, high radio power sources. Other survey goals include investigating the dual--population unification schemes for radio sources, studying the radio luminosity dependence of the evolution of radio source environments, and understanding the radio power dependence of the K-z relation for radio galaxies.

We present a brief overview of the ongoing Complete Optical and Radio Absorption Line System (CORALS) survey. We review the results of our survey for z>2.2 damped Lyman alpha systems (DLAs) and present the preliminary results of an extended survey to lower redshifts (0.7<z<1.5). Although final observations are still pending, our preliminary statistics do not find an excess of absorption systems selected via MgII towards optically faint QSOs. We also briefly review the results from some spin-off projects from CORALS including the low alpha/Fe abundances of multiple DLAs and the possible excess of absorbers within 3000 km/s of the QSO redshift.

We present optical spectra and redshift measurements for 178 flat-spectrum objects from the Parkes quarter-Jansky flat-spectrum sample. These spectra were obtained in order to compile a complete sample of quasars for use in a study of quasar evolution. We present a composite optical spectrum made from the subset of 109 quasars that have flux densities in the range 0.25Jy < S(2.7GHz)< 0.5Jy, and make a comparison with a composite for radio-quiet QSOs from the Large Bright Quasar Survey. Our large sample of radio-loud quasars allows us to strengthen previous reports that the Ly-alpha and CIV emission lines have larger equivalent width in radio-loud quasars than radio-quiet QSOs to greater than the 3-sigma level. However we see no significant difference in the equivalent widths of CIII] or MgII. We also show that the flux decrements across the Lyman-alpha line (D_A) measured from these spectra show the same trend with redshift as for optically selected QSOs.

We quantify the angular clustering of radio galaxies in the NVSS and FIRST radio surveys using the two-point correlation function and the moments of counts-in-cells - both important points of comparison with theory. These investigations consistently demonstrate that the slope of the correlation function for radio galaxies agrees with that for optically-selected galaxies, gamma = 1.8. We describe how to disentangle the imprint of galaxy clustering from the two observational problems: resolution of radio galaxies into multiple components and gradients in source surface density induced by difficulties in processing "snapshot" radio observations (significant in both surveys below 15 mJy). This study disagrees in some respects with previous analyses of the angular clustering of radio galaxies.

We are in motion against the cosmic backdrop. This motion is evidenced by the systematic temperature shift - or dipole anisotropy - observed in the Cosmic Microwave Background radiation (CMB). Because of the Doppler effect, the temperature of the CMB is 0.1 per cent higher in our direction of motion through the Universe. If our standard cosmological understanding is correct, this dipole should also be present as an enhancement in the surface density of distant galaxies. The main obstacle in finding this signal is the very uneven distribution of nearby galaxies in the Local Supercluster, which drowns out the small cosmological imprint. Here we report the first detection of the expected dipole anisotropy in the galaxy distribution, in a survey of galaxies detected in radio waves. Radio galaxies are mostly located at cosmological distances, so the contamination from nearby clusters should be small. With local radio sources removed, we find a dipole anisotropy in the radio galaxy distribution in the same direction as the CMB, close to the expected amplitude. This result is confirmation of the standard cosmological interpretation of the CMB.

We present a list of z_abs ~ z_em (Delta v < 3000 km/s) DLAs discovered during the CORALS survey for absorbers in a radio-selected QSO sample. On the assumption that z_abs ~ z_em DLAs are neither ejecta from the QSO, nor part of the host galaxy itself, we use the number density statistics of these DLAs to investigate galaxy clustering near the QSO redshift. We find that n(z) of z_abs \~z_em DLAs in our radio-selected QSO sample is ~4 times larger than the number density of intervening DLAs, implying an excess of galaxies near the QSO. This result is further supported with the inclusion of the radio-loud QSOs in the FIRST survey, although the total number of DLAs is still small (4) and the result is currently only significant at approximately the 2 sigma level. Since all of the z_abs ~ z_em DLAs we identify in CORALS are found towards optically bright (B < 20) QSOs, there is no strong evidence (based on these limited statistics) that this population suffers from a severe dust bias. We compare our results with those from an optically-selected, radio-quiet QSO sample in order to determine whether there is evidence for an excess of galaxies around radio-loud versus radio-quiet QSOs. We find that the n(z) of z_abs ~ z_em DLAs towards radio-quiet QSOs is in agreement with the number density of intervening absorbers. This result, although currently limited by the small number statistics of our survey, supports the conclusion that radio-loud quasars are found preferentially in rich galaxy environments at high redshift. Finally, we report that one of the new z_abs ~ z_em DLAs discovered by CORALS has some residual flux in the base of the Lya trough which may be due to Lya emission, either from star formation in the DLA galaxy or from gas photoionised by the QSO.

We quantify the angular distribution of radio sources in the NVSS by measuring the two-point angular correlation function. By careful consideration of the resolution of radio galaxies into multiple components, we are able to determine both the galaxy angular clustering and the size distribution of giant radio galaxies. The slope of the correlation function for radio galaxies agrees with that for other classes of galaxy, gamma = 1.8, with a 3D correlation length r_0 = 6 h^-1 Mpc (under certain assumptions). Calibration problems in the survey prevent clustering analysis below 10 mJy. About 7 per cent of radio galaxies are resolved by NVSS into multiple components, with a power-law size distribution. Our work calls into question previous analyses and interpretations of the angular correlation function from radio surveys.

(Abridged) We present the first results from the Complete Optical and Radio Absorption Line System (CORALS) survey. We have compiled a homogeneous sample of radio-selected QSOs from the Parkes Catalogue and searched for damped Lyman alpha systems (DLAs) towards every target, irrespective of its optical magnitude. This approach circumvents selection effects -- particularly from intervening dust -- which have long been suspected to affect DLA surveys in optically-selected, magnitude-limited QSO samples. The CORALS data set consists of 66 z_em > 2.2 QSOs in which 22 DLAs with absorption redshifts 1.8 < z_abs < z_em have been identified over a total redshift interval Delta z = 55.46. In this first paper of the CORALS series we describe the sample, present intermediate resolution spectroscopy and determine the population statistics of DLAs. We deduce a value of the neutral gas mass density traced by DLAs (expressed as a fraction of the closure density) log Omega_DLA h = -2.59^{+0.17}_{-0.24}, and a number density of DLAs per unit redshift n(z) = 0.31^{+0.09}_{-0.08}, both at a mean redshift <z> = 2.37. Taking into account the errors, we conclude that dust-induced bias in previous surveys may have led to an underestimate of these quantities by at most a factor of two and we have not uncovered a previously unrecognised population of high column density DLAs in front of faint QSOs.

We present the first results from the Complete Optical and Radio Absorption Line System (CORALS) survey for Damped Lyman Alpha (DLA) systems based on a radio-selected sample of quasars. We find that the number density and neutral gas content of CORALS DLAs is in good agreement with previous surveys. These results indicate that dust obscuration has probably not been responsible for a significant selection bias in the past, although the statistics still permit an underestimate of n(z) and Omega_DLA by up to a factor of about 2.

Observations of the extraordinarily bright optical afterglow (OA) of GRB 991208 started 2.1 d after the event. The flux decay constant of the OA in the R-band is -2.30 +/- 0.07 up to 5 d, which is very likely due to the jet effect, and after that it is followed by a much steeper decay with constant -3.2 +/- 0.2, the fastest one ever seen in a GRB OA. A negative detection in several all-sky films taken simultaneously to the event implies either a previous additional break prior to 2 d after the occurrence of the GRB (as expected from the jet effect). The existence of a second break might indicate a steepening in the electron spectrum or the superposition of two events. Once the afterglow emission vanished, contribution of a bright underlying SN is found, but the light curve is not sufficiently well sampled to rule out a dust echo explanation. Our determination of z = 0.706 indicates that GRB 991208 is at 3.7 Gpc, implying an isotropic energy release of 1.15 x 10E53 erg which may be relaxed by beaming by a factor > 100. Precise astrometry indicates that the GRB coincides within 0.2" with the host galaxy, thus given support to a massive star origin. The absolute magnitude is M_B = -18.2, well below the knee of the galaxy luminosity function and we derive a star-forming rate of 11.5 +/- 7.1 Mo/yr. The quasi-simultaneous broad-band photometric spectral energy distribution of the afterglow is determined 3.5 day after the burst (Dec 12.0) implying a cooling frequency below the optical band, i.e. supporting a jet model with p = -2.30 as the index of the power-law electron distribution.

In this paper we discuss the demographics of the radio blazar population: (i) what are their parent (`unbeamed') sources and (ii) what magnitude and/or type of evolution have they undergone ? The discussion is based on models of radio source evolution and beaming based on a `dual population' unification paradigm. These models, developed from radio blazar properties in bright samples, predict blazar demographic trends at the lower flux-density levels; samples from deep mJy-level surveys (e.g. NVSS and FIRST) may now provide direct tests of these predictions.

Radio spectra of radio galaxies are often ascribed a simple power law form ($S_{\nu} \propto \nu^{\alpha}$) with the spectral index, $\alpha$, being of order -0.7. However, all radio galaxies deviate from this simple power law behaviour. In this paper we derive simple expressions for the average rest-frame spectra of FRI and FRII radio galaxies. These will be used to describe the spectral curvature of the parent (FRI and FRII) populations in models of radio source evolution.

It has been claimed by Taylor et al. that the low-redshift end of the K-z relation for radio galaxies is too bright by about half a magnitude due to contributions from the obscured quasar nuclei. Such a result has major implications for the use of the K-band Hubble diagram in understanding the cosmological evolution of radio galaxies. In this paper we present 1-5 micron imaging data of a nearly-complete sample of low-redshift radio galaxies; this approach allows us to accurately determine the strengths of any unresolved nuclear components in the galaxies. We detect nuclear sources in five targets, whose broad-band colours are consistent with reddened quasar spectra. In all five cases the ratio of the inferred intrinsic near-infrared luminosity to the narrow-line luminosity is typical of quasars. We find a correlation between the inferred nuclear extinction and core-to-lobe ratio, which places constraints on the geometry of the torus. We find evidence for a shift of the K-z relation to fainter magnitudes, but by a much smaller amount (~0.1 mag) than Taylor et al. determined. Under the assumption that the nuclear sources in radio galaxies have the same intrinsic near-infrared spectra as quasars, our multi-wavelength images allow us to limit any possible shift to less than 0.3 magnitudes.

We present spectra for a sample of radio sources from the FIRST survey, and use them to define the form of the redshift distribution of radio sources at mJy levels.We targeted 365 sources and obtained 46 redshifts (13 per cent of the sample). We find that our sample is complete in redshift measurement to R $\sim 18.6$, corresponding to $z\sim 0.2$. Early-type galaxies represent the largest subset (45 per cent) of the sample and have redshifts $0.15\la z \la 0.5$ ; late-type galaxies make up 15 per cent of the sample and have redshifts $0.05\la z \la 0.2$; starbursting galaxies are a small fraction ($\sim 6$ per cent), and are very nearby ($z\la 0.05$). Some 9 per cent of the population have Seyfert1/quasar-type spectra, all at $z\ga 0.8$, and there are 4 per cent are Seyfert2 type galaxies at intermediate redshifts ($z\sim 0.2$). Using our measurements and data from the Phoenix survey, we obtain an estimate for $N(z)$ at $S_{1.4 \rm {GHz}}\ge 1$ mJy and compare this with model predictions. At variance with previous conclusions, we find that the population of starbursting objects makes up $\la 5 $ per cent of the radio population at S $\sim 1$ mJy.

We report on two programs to address differential evolution between the radio-loud and radio-quiet quasar populations at high (z>4) redshift. Both programs entail studying the radio properties of optically-selected quasars. First, we have observed 32 optically-selected, high-redshift (z>4) quasars with the VLA at 6 cm (5 GHz). These sources comprise a statistically complete and well-understood sample. We detect four quasars above our 3-sigma limit of ~0.15 mJy, which is sufficiently sensitive to detect all radio-loud quasars at the probed redshift range. Second, we have correlated 134 z>4 quasars, comprising all such sources that we are aware of as of mid-1999, with FIRST and NVSS. These two recent 1.4 GHz VLA sky surveys reach 3-sigma limits of approximately 0.6 mJy and 1.4 mJy respectively. We identify a total of 15 z>4 quasars, of which six were not previously known to be radio-loud. The depth of these surveys does not reach the radio-loud/radio-quiet demarcation luminosity density (L(1.4 GHz) = 10^32.5 h(50)^(-2) ergs/s/Hz) at the redshift range considered; this correlation therefore only provides a lower limit to the radio-loud fraction of quasars at high-redshift. The two programs together identify eight new radio-loud quasars at z>4, a significant increase over the seven currently in the published literature. We find no evidence for radio-loud fraction depending on optical luminosity for -25 > M_B > -28 at z~2, or for -26>M_B>-28 at z>4. Our results also show no evolution in the radio-loud fraction between z~2 and z>4 (-26>M_B>-28).

A new method of deriving the cosmological evolution of radio AGN is described which is based on the premise of unified models for quasars and radio galaxies.

We discuss how different theoretical predictions for the variance $\sigma^2$ of the distribution of radio sources can be matched to measurements from the FIRST survey at different flux limits. The predictions are given by the integration of models for the angular correlation function $w(\theta)$ for three different functional forms of the redshift distribution $N(z)$, different spatial correlation functions and by different evolutions of the bias $b(z)$ with redshift. We also consider the two cases of open and flat Universes. Although the predicted $w(\theta)$ show substantial differences due to differences in the $N(z)$'s, these differences are not significant compared to the uncertainties in the current observations. It turns out that the best fit is provided by models with constant biasing at all times, although the difference between models with epoch-independent bias and models with bias that evolves linearly with redshift is not very large. All models with strong evolution of bias with epoch are ruled out. As a further step we directly calculated $w_{obs}(\theta)$ at 3mJy from the catalogue and matched it with our models for the angular correlation function in the hypothesis that the clustering signal comes from two different populations, namely AGN-powered sources and starbursting galaxies. The results are consistent with a scenario for hierarchical clustering where the fainter starbursting galaxies trace the mass at all epochs, while brighter AGN's are strongly biased, with $b(z)$ evolving linearly with redshift, as suggested by some theories of galaxy formation and evolution.

We show that a dual-population unification scheme provides a successful paradigm with which to describe the evolution and beaming of all bright extragalactic radio sources. The paradigm consists of two intrinsic radio-source populations, based on the two distinct radio-galaxy morphologies of Fanaroff-Riley classes I and II. These represent the `unbeamed' or `side-on' parent populations of steep radio spectra; the `beamed' source types including flat-spectrum quasars and BL Lac objects, arise through the random alignment of their radio-axis to our line-of-sight where Doppler-beaming of the relativistic radio jets produces highly anisotropic radio emission.

Extragalactic radio sources are a unique cosmological probe in that they trace large-scale structure on scales inaccessible to other wavelengths. However as radio survey data is inherently 2D, the redshift distribution, N(z), is necessary to derive spatial information. To obtain this distribution either we measure thousands of radio source redshifts to directly determine N(z) or we derive N(z) from statistical analyses of radio source count and identification data. In this paper we show how the dual-population unification scheme can be incorporated into a rigorous statistical analysis of radio source count data, with the result that our simple parametric evolution and beaming model revises previous estimations of N(z), specifically at low flux densites. This revision is particularly pertinent given that the new generation of radio surveys extend to milli-jansky flux density levels: sampling source densities high enough to reveal spatial structure. In turn, these new radio surveys will provide potent tests which will refine our model.

We present the population data as predicted by our space density analysis at the SUMSS frequency of 843 MHz. This data demonstrates the potential of the SUMSS survey to trace in detail the change-over in radio source populations from `radio monsters' at high flux densities to a local starbursting population at the survey limit. The exact form of the observed change-over will be used to refine our models of radio source evolution.

We investigate the large-scale clustering of radio sources in the FIRST 1.4-GHz survey by analysing the distribution function (counts in cells). We select a reliable sample from the the FIRST catalogue, paying particular attention to the problem of how to define single radio sources from the multiple components listed. We also consider the incompleteness of the catalogue. We estimate the angular two-point correlation function $w(\theta)$, the variance $\Psi_2$, and skewness $\Psi_3$ of the distribution for the various sub-samples chosen on different criteria. Both $w(\theta)$ and $\Psi_2$ show power-law behaviour with an amplitude corresponding a spatial correlation length of $r_0 \sim 10 h^{-1}$Mpc. We detect significant skewness in the distribution, the first such detection in radio surveys. This skewness is found to be related to the variance through $\Psi_3=S_3(\Psi_2)^{\alpha}$, with $\alpha=1.9\pm 0.1$, consistent with the non-linear gravitational growth of perturbations from primordial Gaussian initial conditions. We show that the amplitude of variance and skewness are consistent with realistic models of galaxy clustering.

We have obtained complete redshift information for a sample of 442 radio-loud quasars with flux densities S(11cm) >= 0.25Jy. These come from a completely-identified sample of 878 flat-spectrum radio sources, so there is no optical magnitude limit. With these quasars, therefore, we can map out the entire quasar epoch, unhindered by optical selection effects and any intervening dust. The rapid decline in true space density above redshifts of about 2.5 is confirmed. It agrees closely with that obtained from optically-selected samples, indicating that these are also not significantly affected by dust. The evolution of the quasar space density is similar to that deduced for the global star formation rate, implying that the quasar epoch may also be that of star formation. This sample of radio-selected quasars can now be used in a variety of other studies, including unbiased searches for damped Ly-alpha absorbers, studies of the phenomenon of ``red quasars'', and searches for high-redshift molecular absorption systems.

(Abridged Abstract) We calculate the expected amplitude of the dipole and higher spherical harmonics in the angular distribution of radio galaxies. The median redshift of radio sources in existing catalogues is z=1, which allows us to study large scale structure on scales between those accessible to present optical and infrared surveys, and that of the Cosmic Microwave Background (CMB). The dipole is due to 2 effects which turn out to be of comparable magnitude: (i) our motion with respect to the CMB, and (ii) large scale structure, parameterised here by a family of Cold Dark Matter power-spectra. We make specific predictions for the Green Bank (87GB) and Parkes-MIT-NRAO (PMN) catalogues. For these relatively sparse catalogues both the motion and large scale structure dipole effects are expected to be smaller than the Poisson shot-noise. However, we detect dipole and higher harmonics in the combined 87GB-PMN catalogue which are far larger than expected. We attribute this to a 2 % flux mismatch between the two catalogues. We also investigate the existence and extent of the Supergalactic Plane in the above catalogues. In a strip of +- 10 deg of the standard Supergalactic equator, we find a 3-sigma detection in PMN, but only 1-sigma in 87GB.

We investigate the large-scale clustering of radio sources in the Green Bank and Parkes-MIT-NRAO 4.85 GHz surveys by measuring the angular two-point correlation function w(\theta). Excluding contaminated areas, the two surveys together cover 70 per cent of the whole sky. We find both surveys to be reasonably complete above 50 mJy. On the basis of previous studies, the radio sources are galaxies and radio-loud quasars lying at redshifts up to z \sim 4, with a median redshift z \sim 1. This provides the opportunity to probe large-scale structures in a volume far larger than that within the reach of present optical and infrared surveys. We detect a clustering signal w(\theta) \approx 0.01 for \theta = 1\degr. By assuming an evolving power-law spatial correlation function in comoving coordinates \xi(r_c,z) = ( r_c / r_0 )^{-\gamma} (1+z)^{\gamma-(3+\epsilon)}, where \gamma = 1.8, and the redshift distribution N(z) of the radio galaxies, we constrain the r_0--\epsilon parameter space. For `stable clustering' (\epsilon = 0), we find the correlation length r_0 \approx 18 Mpc/h, larger than the value for nearby normal galaxies and comparable to the cluster-cluster correlation length.