Understanding the evolution of accretion activity is fundamental to our understanding of how galaxies form and evolve over the history of the Universe. We analyse a complete sample of 27 radio galaxies which includes both high-excitation (HEGs) and low excitation galaxies (LEGs), spanning a narrow redshift range of 0.9 < z < 1.1 and covering a factor of ~1000 in radio luminosity. Using data from the Spitzer Space Telescope combined with ground-based optical and near-infrared imaging, we show that the host galaxies have masses in the range of 10.7 < log (M /M_sun) < 12.0 with HEGs and LEGs exhibiting no difference in their mass distributions. We also find that HEGs accrete at significantly higher rates than LEGs, with the HEG/LEG division lying at an Eddington ratio of ~0.04, which is in excellent agreement with theoretical predictions of where the accretion rate becomes radiatively inefficient, thus supporting the idea of HEGs and LEGs being powered by different modes of accretion. Our study also shows that at least up to L_151MHz ~3x10^27 W /Hz /sr, HEGs and LEGs are indistinguishable in terms of their radio properties. From this result we infer that, at least for the lower radio luminosity range, another factor besides accretion rate must play an important role in the process of triggering jet activity.

We model a 21 cm intensity mapping survey in the redshift range 0.01<z<1.5 designed to simulate the skies as seen by future radio telescopes such as the Square Kilometre Array (SKA), including instrumental noise and Galactic foregrounds. In our pipeline, we remove the introduced Galactic foregrounds with a fast independent component analysis (fastica) technique. We present the power spectrum of the large-scale matter distribution, C(l), before and after the application of this foreground removal method and calculate the resulting systematic errors. We attempt to reduce systematics in the foreground subtraction by optimally masking the maps to remove high foregrounds in the Galactic plane. Our simulations show a certain level of bias remains in the power spectrum at all scales l<400. At large-scales l<30 this bias is particularly significant. We measure the impact of these systematic effects in two different ways: firstly we fit cosmological parameters to the broadband shape of the power spectrum and secondly we extract the position of the Baryon Acoustic Oscillations (BAO). In the first analysis, we find that the systematics introduce an significant shift in the best fit cosmological parameters at the 2 to 3 sigma level which depends on the masking and noise levels. However, cosmic distances can be recovered in an unbiased way after foreground removal at all simulated redshifts by fitting the BAOs in the power spectrum. We conclude that further advances in foreground removal are needed in order to recover unbiased information from the broadband shape of the power spectrum, however, intensity mapping experiments will be a powerful tool for mapping cosmic distances across a wide redshift range.

We predict the emission line luminosity functions (LFs) of the first 10 rotational transitions of CO in galaxies at redshift z=0 to z=10. This prediction relies on a recently presented simulation of the molecular cold gas content in ~3e7 evolving galaxies based on the Millennium Simulation. We combine this simulation with a model for the conversion between molecular mass and CO-line intensities, which incorporates the following mechanisms: (i) molecular gas is heated by the CMB, starbursts (SBs), and active galactic nuclei (AGNs); (ii) molecular clouds in dense or inclined galaxies can overlap; (iii) compact gas can attain a smooth distribution in the densest part of disks; (iv) CO-luminosities scale with metallicity changes between galaxies; (v) CO-luminosities are always detected against the CMB. We analyze the relative importance of these effects and predict the cosmic evolution of the CO-LFs. The most notable conclusion is that the detection of regular galaxies (i.e. no AGN, no massive SB) at high z>7 in CO-emission will be dramatically hindered by the weak contrast against the CMB, in contradiction to earlier claims that CMB-heating will ease the detection of high-redshift CO. The full simulation of extragalactic CO-lines and the predicted CO-LFs at any redshift can be accessed online, prior registration required} and they should be useful for the modeling of CO-line surveys with future telescopes, such as ALMA, the LMT, or the SKA.

We forecast the abilities of the Atacama Large Millimeter/submillimeter Array (ALMA) and the Square Kilometer Array (SKA) to detect CO and HI emission lines in galaxies at redshift z=3. A particular focus is set on Milky Way (MW) progenitors at z=3 for their detection within 24 h constitutes a key science goal of ALMA. The analysis relies on a semi-analytic model, which permits the construction of a MW progenitor sample by backtracking the cosmic history of all simulated present-day galaxies similar to the real MW. Results: (i) ALMA can best observe a MW at z=3 by looking at CO(3-2) emission. The probability of detecting a random model MW at 3-sigma in 24 h using 75 km/s channels is roughly 50%, and these odds can be increased by co-adding the CO(3-2) and CO(4-3) lines. These lines fall into ALMA band 3, which therefore represents the optimal choice towards MW detections at z=3. (ii) Higher CO transitions contained in the ALMA bands geq6 will be invisible, unless the considered MW progenitor coincidentally hosts a major starburst or an active black hole. (iii) The high-frequency array of SKA, fitted with 28.8 GHz receivers, would be a powerful instrument for observing CO(1-0) at z=3, able to detect nearly all simulated MWs in 24 h. (iv) HI detections in MWs at z=3 using the low-frequency array of SKA will be impossible in any reasonable observing time. (v) SKA will nonetheless be a supreme ha survey instrument through its enormous instantaneous field-of-view (FoV). A one year pointed HI survey with an assumed FoV of 410 sqdeg would reveal at least 10^5 galaxies at z=2.95-3.05. (vi) If the positions and redshifts of those galaxies are known from an optical/infrared spectroscopic survey, stacking allows the detection of HI at z=3 in less than 24 h.

An analysis of 44 GHz VLA observations of the z = 1.574 radio-loud quasar 3C318 has revealed emission from the redshifted J = 1 - 0 transition of the CO molecule and spatially resolved the 6.3 kpc radio jet associated with the quasar at 115 GHz rest-frame. The continuum-subtracted line emitter is spatially offset from the quasar nucleus by 0.33" (2.82 kpc in projection). This spatial offset has a significance of >8-sigma and, together with a previously published -400 km/s velocity offset measured in the J = 2 - 1 CO line relative to the systemic redshift of the quasar, rules out a circumnuclear starburst or molecular gas ring and suggests that the quasar host galaxy is either undergoing a major merger with a gas-rich galaxy or is otherwise a highly disrupted system. If the merger scenario is correct then the event may be in its early stages, acting as the trigger for both the young radio jets in the quasar and a starburst in the merging galaxy. The total molecular gas mass in the spatially offset line emitter as measured from the ground-state CO line M_H2 = 3.7 (+/-0.4) x 10^10 (alpha_CO/0.8) M_solar. Assuming that the line-emitter can be modelled as a rotating disk, an inclination-dependent upper limit is derived for its dynamical mass M_dyn sin^2(i) < 3.2 x 10^9 M_solar, suggesting that for M_H2 to remain less than M_dyn the inclination angle must be i < 16 degrees. The far infrared and CO luminosities of 246 extragalactic systems are collated from the literature for comparison. The high molecular gas content of 3C318 is consistent with that of the general population of high redshift quasars and sub-millimetre galaxies.

We present JVLA observations of the cold (CO (1-0)) molecular gas in IRAS F10214+4724, a lensed ULIRG at z=2.3 with an obscured active nucleus. The galaxy is spatially and spectrally well-resolved in the CO (1-0) emission line. A CO (1-0) counter-image is detected at the 3-sigma level. Five of the 42 km/s channels (with >5-sigma detections) are mapped back into the source plane and their total magnification posterior PDFs sampled. This reveals a roughly linear arrangement, tentatively a rotating disk. We derive a molecular gas mass of M_gas = 1.2 +- 0.2 x 10^10 M_sun, assuming a ULIRG L_{CO}-to-M_{gas} conversion ratio of \alpha = 0.8 M_sun / (K km/s pc^2) that agrees well with the derived range of \alpha = 0.3 - 1.3 for separate dynamical mass estimates at assumed inclinations of i = 90 - 30 degrees. Based on the AGN and CO (1-0) peak emission positions and the lens model, we predict a distortion of the CO Spectral Line Energy Distribution (SLED) where higher order J lines that may be partially excited by AGN heating will be preferentially lensed owing to their smaller solid angles and closer proximity to the AGN and therefore the cusp of the caustic. Comparison with other lensing inversion results shows that the narrow line region and AGN radio core in IRAS F10214+4724 are preferentially lensed by a factor >~ 3 and 11 respectively, relative to the molecular gas emission. This distorts the global continuum emission Spectral Energy Distribution (SED) and suggests caution in unsophisticated uses of IRAS F10214+4724 as an archetype high-redshift ULIRG. We explore two Large Velocity Gradient (LVG) models, incorporating spatial CO (1-0) and (3-2) information and present tentative evidence for an extended, low excitation cold gas component that implies that the total molecular gas mass in IRAS F10214+4724 is a factor >~2 greater than that calculated using spatially unresolved CO observations.

This is the first paper in a series that present a multi-wavelength analysis of the archetype Ultra-Luminous InfraRed Galaxy (ULIRG) IRAS FSC10214+4724, a gravitationally lensed, starburst/AGN at z=2.3. Here we present a new lens model and spatially-resolved radio data, as well as a deep HST F160W map. The lens modelling employs a Bayesian Markov Chain Monte Carlo algorithm with extended-source, forward ray-tracing. Using these high resolution HST, MERLIN and VLA maps, the algorithm allows us to constrain the level of distortion to the continuum spectral energy distribution resulting from emission components with differing magnification factors, due to their size and proximity to the caustic. Our lens model finds the narrow line region (NLR), and by proxy the active nucleus, is preferentially magnified. This supports previous claims that preferential magnification could mask the expected polycyclic aromatic hydrocarbon spectral features in the Spitzer mid-infrared spectrum which roughly trace the star-forming regions. Furthermore, we show the arc-to-counter-image flux ratio is not a good estimate of the magnification in this system, despite its common use in the IRAS FSC10214+4724 literature. Our lens modelling suggests magnifications of \mu ~ 15-20+-2 for the HST F814W, MERLIN 1.7 GHz and VLA 8 GHz maps, significantly lower than the canonical values of \mu = 50-100 often used for this system. Systematic errors such as the dark matter density slope and co-location of stellar and dark matter centroids dominate the uncertainties in the lens model at the 40 percent level.

We report 1.7 GHz Very Long Baseline Interferometry (VLBI) observations of IRAS F10214+4724, a lensed z=2.3 obscured quasar with prodigious star formation. We detect what we argue to be the obscured active nucleus with an effective angular resolution of < 50 pc at z = 2.3 . The S_{1.7} = 210 micro-Jy (9-\sigma) detection of this unresolved source is located within the HST rest-frame ultraviolet/optical arc, however, >~100 mas northward of the arc centre of curvature. This leads to a source plane inversion that places the European VLBI Network detection to within milli-arcseconds of the modelled cusp caustic, resulting in a very large magnification (\mu ~70), over an order of magnitude larger than the CO (1-0) derived magnification of a spatially resolved JVLA map, using the same lens model. We estimate the quasar bolometric luminosity from a number of independent techniques and with our X-ray modelling find evidence that the AGN may be close to Compton-thick, with an intrinsic bolometric luminosity log(L_{bol,QSO} / L_sun) = 11.34 +- 0.27 dex. We make the first black hole mass estimate of IRAS F10214+4724 and find log(M_{BH}/M_sun) = 8.36 +- 0.56 which suggests a low black hole accretion rate (\lambda = \dot{M} / \dot{M}_{Edd} ~ 3\pm^7_2 percent). We find evidence for a M_{BH}/M_{spheroid} ratio that is 1-2 orders of magnitude larger than that of submillimetre galaxies (SMGs) at z~2. At face value, this suggests IRAS F10214+4724 has undergone a different evolutionary path compared to SMGs at the same epoch. A primary result of this work is the demonstration that emission regions of differing size and position can undergo significantly different magnification boosts (> 1 dex) and therefore distort our view of high-redshift, gravitationally lensed galaxies.

We describe a 325-MHz survey, undertaken with the Giant Metrewave Radio Telescope (GMRT), which covers a large part of the three equatorial fields at 9, 12 and 14.5 h of right ascension from the Herschel-Astrophysical Terahertz Large Area Survey (H-ATLAS) in the area also covered by the Galaxy And Mass Assembly survey (GAMA). The full dataset, after some observed pointings were removed during the data reduction process, comprises 212 GMRT pointings covering ~90 deg^2 of sky. We have imaged and catalogued the data using a pipeline that automates the process of flagging, calibration, self-calibration and source detection for each of the survey pointings. The resulting images have resolutions of between 14 and 24 arcsec and minimum rms noise (away from bright sources) of ~1 mJy/beam, and the catalogue contains 5263 sources brighter than 5 sigma. We investigate the spectral indices of GMRT sources which are also detected at 1.4 GHz and find them to agree broadly with previously published results; there is no evidence for any flattening of the radio spectral index below S_1.4=10 mJy. This work adds to the large amount of available optical and infrared data in the H-ATLAS equatorial fields and will facilitate further study of the low-frequency radio properties of star formation and AGN activity in galaxies out to z~1.

LOFAR, the LOw-Frequency ARray, is a new-generation radio interferometer constructed in the north of the Netherlands and across europe. Utilizing a novel phased-array design, LOFAR covers the largely unexplored low-frequency range from 10-240 MHz and provides a number of unique observing capabilities. Spreading out from a core located near the village of Exloo in the northeast of the Netherlands, a total of 40 LOFAR stations are nearing completion. A further five stations have been deployed throughout Germany, and one station has been built in each of France, Sweden, and the UK. Digital beam-forming techniques make the LOFAR system agile and allow for rapid repointing of the telescope as well as the potential for multiple simultaneous observations. With its dense core array and long interferometric baselines, LOFAR achieves unparalleled sensitivity and angular resolution in the low-frequency radio regime. The LOFAR facilities are jointly operated by the International LOFAR Telescope (ILT) foundation, as an observatory open to the global astronomical community. LOFAR is one of the first radio observatories to feature automated processing pipelines to deliver fully calibrated science products to its user community. LOFAR's new capabilities, techniques and modus operandi make it an important pathfinder for the Square Kilometre Array (SKA). We give an overview of the LOFAR instrument, its major hardware and software components, and the core science objectives that have driven its design. In addition, we present a selection of new results from the commissioning phase of this new radio observatory.

We perform a stacking analysis of H-ATLAS data in order to obtain isothermal dust temperatures and rest-frame luminosities at 250um (L250), for 1599 radio sources over the H-ATLAS P1 GAMA area. The radio sample is generated using a combination of NVSS data and K-band UKIDSS-LAS data, over 0.01<z<0.8. The FIR properties of the sample are investigated as a function of L_1.4 GHz, z, projected radio-source size and radio spectral index. In order to search for stellar-mass dependent relations, we split the parent sample into those sources which are below and above 1.5 L_{K}^{*}. Correcting for stellar mass and redshift, we find no relation between the L250 and L_1.4GHz of radio AGN. This implies that a galaxy's nominal radio luminosity has little or no bearing on the SFR and/or dust mass content of the host system. The L250 of both the radio detected and non radio-detected galaxies (defined as those sources not detected at 1.4GHz but detected in the SDSS with r'<22) rises with increasing z. For mass and colour matched samples, sub-1.5 L_{K}^{*} and super-1.5 L_{K}^{*} radio-detected galaxies have 0.89\pm0.18 and 0.49\pm0.12 times the L250 of their non radio-detected counterparts. We explain these results in terms of the hotter, denser and richer halo environments massive radio-galaxies maintain and are embedded in. These results indicate that all massive radio galaxies (>1.5 L_{K}^{*}) may have systematically lower FIR luminosities (~25%) than their colour-matched non radio-detected counterparts. Compact radio sources (<30kpc) are associated with higher L250s and dust temperatures than their more extended (>30kpc) counterparts. The higher dust temperature suggests that this may be attributed to enhanced SFRs, but whether this is directly or indirectly due to radio activity (e.g. jet induced or merger-driven SF) is as yet unknown.

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)

Aims: This paper discusses the spectral occupancy for performing radio astronomy with the Low-Frequency Array (LOFAR), with a focus on imaging observations. Methods: We have analysed the radio-frequency interference (RFI) situation in two 24-h surveys with Dutch LOFAR stations, covering 30-78 MHz with low-band antennas and 115-163 MHz with high-band antennas. This is a subset of the full frequency range of LOFAR. The surveys have been observed with a 0.76 kHz / 1 s resolution. Results: We measured the RFI occupancy in the low and high frequency sets to be 1.8% and 3.2% respectively. These values are found to be representative values for the LOFAR radio environment. Between day and night, there is no significant difference in the radio environment. We find that lowering the current observational time and frequency resolutions of LOFAR results in a slight loss of flagging accuracy. At LOFAR's nominal resolution of 0.76 kHz and 1 s, the false-positives rate is about 0.5%. This rate increases approximately linearly when decreasing the data frequency resolution. Conclusions: Currently, by using an automated RFI detection strategy, the LOFAR radio environment poses no perceivable problems for sensitive observing. It remains to be seen if this is still true for very deep observations that integrate over tens of nights, but the situation looks promising. Reasons for the low impact of RFI are the high spectral and time resolution of LOFAR; accurate detection methods; strong filters and high receiver linearity; and the proximity of the antennas to the ground. We discuss some strategies that can be used once low-level RFI starts to become apparent. It is important that the frequency range of LOFAR remains free of broadband interference, such as DAB stations and windmills.

[Abridged] Sub-mm observations of the William Herschel Deep Field using LABOCA revealed possible counterparts for 2 X-ray absorbed QSOs. The aim here is to exploit EVLA imaging at 8.4 GHz to establish the QSOs as radio/sub-mm sources. The challenge in reducing the EVLA data was the presence of a strong 4C source in the field. A new calibration algorithm was applied to the data to subtract it. The resulting thermal noise limited radio map covers the 16'x16' Extended WHDF. It contains 41 sources above a 4-sigma limit, 17 of which have primary beam corrected flux. The radio observations show that the absorbed AGN with LABOCA detections are coincident with radio sources, confirming the tendency for X-ray absorbed AGN to be sub-mm bright. These sources show strong ultraviolet excess (UVX) suggesting the nuclear sightline is gas- but not dust-absorbed. Of the 3 remaining LABOCA sources within the ~5' half-power beam width, 1 is identified with a faint nuclear X-ray/radio source in a nearby galaxy, 1 with a faint radio source and 1 is unidentified in any other band. More generally, differential radio source counts are in good agreement with previous observations, showing at S<50 micro-Jy a significant excess over a pure AGN model. In the full area, of 10 sources fainter than this limit, 6 have optical counterparts of which 3 are UVX (i.e. likely QSOs) including the 2 absorbed quasar LABOCA sources. The other faint radio counterparts are not UVX but are only slightly less blue and likely to be star-forming/merging galaxies, predominantly at lower luminosities and redshifts. The 4 faint, optically unidentified radio sources may be either dust obscured QSOs or galaxies. These high-z obscured AGN and lower-z star-forming populations are thus the main candidates to explain the observed excess in faint source counts and hence the excess radio background found previously by the ARCADE2 experiment.

We present the Spitzer Extragalactic Representative Volume Survey (SERVS), an 18 square degrees medium-deep survey at 3.6 and 4.5 microns with the post-cryogenic Spitzer Space Telescope to ~2 microJy (AB=23.1) depth of five highly observed astronomical fields (ELAIS-N1, ELAIS-S1, Lockman Hole, Chandra Deep Field South and XMM-LSS). SERVS is designed to enable the study of galaxy evolution as a function of environment from z~5 to the present day, and is the first extragalactic survey both large enough and deep enough to put rare objects such as luminous quasars and galaxy clusters at z>1 into their cosmological context. SERVS is designed to overlap with several key surveys at optical, near- through far-infrared, submillimeter and radio wavelengths to provide an unprecedented view of the formation and evolution of massive galaxies. In this paper, we discuss the SERVS survey design, the data processing flow from image reduction and mosaicing to catalogs, as well as coverage of ancillary data from other surveys in the SERVS fields. We also highlight a variety of early science results from the survey.

We present new detections of the CO(5-4), CO(7-6), [CI](1-0) and [CI](2-1) molecular and atomic line transitions towards the unlensed, obscured quasar AMS12 (z=2.7672), observed with the IRAM PdBI. This is the first unlensed, high redshift source to have both [CI] transitions detected. Continuum measurements between 70 $\mu$m and 3 mm are used to constrain the FIR SED, and we find a best fit FIR luminosity of log[Lfir/Lsol] = 13.5+/-0.1, dust temperature T_d = 88+/-8 K and emissivity index {\beta} = 0.6+/-0.1. The highly-excited molecular gas probed by CO(3-2), (5-4) and (7-6), is modelled with large velocity gradient (LVG) models. The gas kinetic temperature T_g, density n(H2), and the characteristic size r0, are determined using the dust temperature from the FIR SED as a prior for the gas temperature. The best fitting parameters are T_g = 90+/-8 K, n(H2) = 10^(3.9+/-0.1) cm^(-3) and r0 = 0.8+/-0.04 kpc. The ratio of the [CI] lines gives a [CI] excitation temperature of 43+/-10 K, indicating the [CI] and the high-excitation CO are not in thermal equilibrium. The [CI] excitation temperature is below that of T_d and T_g of the high-excitation CO, perhaps because [CI] lies at a larger radius where there may also be a large reservoir of CO at a cooler temperature, perhaps detectable through the CO(1-0). Using the [CI](1-0) line we can estimate the strength of the CO(1-0) line and hence the gas mass. This suggests that a significant fraction (~30%) of the molecular gas is missed from the high-excitation line analysis. The Eddington limited black hole mass is found from the bolometric luminosity to be Mbh >~ 1.5x10^9 Msol. Along with the stellar mass of 3x10^11 Msol, these give a black hole - bulge mass ratio of Mbh/Mbulge >~ 0.005. This is in agreement with studies on the evolution of the Mbh/Mbulge relationship at high redshifts, which find a departure from the local value ~0.002.

Under the assumption that jets in active galactic nuclei are powered by accretion and the spin of the central supermassive black hole, we are able to reproduce the radio luminosity functions of high- and low-excitation galaxies. High-excitation galaxies are explained as high-accretion rate but very low spin objects, while low-excitation galaxies have low accretion rates and bimodal spin distributions, with approximately half of the population having maximal spins. At higher redshifts (z~1), the prevalence of high accretion rate objects means the typical spin was lower, while in the present day Universe is dominated by low accretion rate objects, with bimodal spin distributions.

We present spectroscopic and eleven-band photometric redshifts for galaxies in the 100-uJy Subaru/XMM-Newton Deep Field radio source sample. We find good agreement between our redshift distribution and that predicted by the SKA Simulated Skies project. We find no correlation between K-band magnitude and radio flux, but show that sources with 1.4-GHz flux densities below ~1mJy are fainter in the near-infrared than brighter radio sources at the same redshift, and we discuss the implications of this result for spectroscopically-incomplete samples where the K-z relation has been used to estimate redshifts. We use the infrared--radio correlation to separate our sample into radio-loud and radio-quiet objects and show that only radio-loud hosts have spectral energy distributions consistent with predominantly old stellar populations, although the fraction of objects displaying such properties is a decreasing function of radio luminosity. We calculate the 1.4-GHz radio luminosity function (RLF) in redshift bins to z=4 and find that the space density of radio sources increases with lookback time to z~2, with a more rapid increase for more powerful sources. We demonstrate that radio-loud and radio-quiet sources of the same radio luminosity evolve very differently. Radio-quiet sources display strong evolution to z~2 while radio-loud AGNs below the break in the radio luminosity function evolve more modestly and show hints of a decline in their space density at z>1, with this decline occurring later for lower-luminosity objects. If the radio luminosities of these sources are a function of their black hole spins then slowly-rotating black holes must have a plentiful fuel supply for longer, perhaps because they have yet to encounter the major merger that will spin them up and use the remaining gas in a major burst of star formation.

We use results from simulations of the production of magnetohydrodynamic jets around black holes to derive the cosmic spin history of the most massive black holes. We assume that the efficiency of jet production is a monotonic function of spin a, as given by the simulations, and that the accretion flow geometry is similarly thick for quasars accreting close to the Eddington ratio and for low-excitation radio galaxies accreting at very small Eddington rates. We use the ratio of the comoving densities of the jet power and the radiated accretion power associated with supermassive black holes with Mbh>~10^8 Msol to estimate the cosmic history of the characteristic spin a. The evolution of this ratio, which increases with decreasing z, is consistent with a picture where the z~0 active galactic nuclei have typically higher spins than those at z~2 (with typical values a~0.35-0.95 and a~0.0-0.25 respectively). We discuss the implications in terms of the relative importance of accretion and mergers in the growth of supermassive black holes with Mbh>~10^8 Msol.

We review how the Square Kilometre Array (SKA) will address fundamental questions in cosmology, focussing on its use for neutral Hydrogen (HI) surveys. A key enabler of its unique capabilities will be large (but smart) receptors in the form of aperture arrays. We outline the likely contributions of Phase-1 of the SKA (SKA1), Phase-2 SKA (SKA2) and pathfinding activities (SKA0). We emphasise the important role of cross-correlation between SKA HI results and those at other wavebands such as: surveys for objects in the EoR with VISTA and the SKA itself; and huge optical and near-infrared redshift surveys, such as those with HETDEX and Euclid. We note that the SKA will contribute in other ways to cosmology, e.g. through gravitational lensing and $H_{0}$ studies.

We review the current status of the Square Kilometre Array (SKA) by outlining the science drivers for its Phase-1 (SKA1) and setting out the timeline for the key decisions and milestones on the way to the planned start of its construction in 2016. We explain how Phase-2 SKA (SKA2) will transform the research scope of the SKA infrastructure, placing it amongst the great astronomical observatories and survey instruments of the future, and opening up new areas of discovery, many beyond the confines of conventional astronomy.

The benefits for the e-MERLIN and EVN arrays of using antennae at the satellite communication station at Goonhilly in Cornwall are discussed. The location of this site - new to astronomy - will provide an almost equal distribution of long baselines in the east-west- and north-south directions, and opens up the possibility to get significantly improved observations of equatorial fields with e-MERLIN. These additional baselines will improve the sensitivity on a set of critical spatial scales and will increase the angular resolution of e-MERLIN by a factor of two. e-MERLIN observations, including many allocated under the e-MERLIN Legacy programme, will benefit from the enhanced angular resolution and imaging capability especially for sources close to or below the celestial equator (where ESO facilities such as ALMA will operate) of including the Goonhilly telescopes. Furthermore, the baselines formed between Goonhilly and the existing stations will close the gap between the baselines of e-MERLIN and those of the European VLBI Network (EVN) and therefore enhance the legacy value of e-MERLIN datasets.

[Abridged] Observations of molecular gas at all redshifts are critical for measuring the cosmic evolution in molecular gas density and understanding the star-formation history of the Universe. The 12CO molecule (J=1-0 transition = 115.27 GHz) is the best proxy for extragalactic H2, which is the gas reservoir from which star formation occurs, and has been detected out to z~6. Typically, redshifted high-J lines are observed at mm-wavelengths, the most commonly targeted systems exhibiting high SFRs (e.g. submm galaxies), and far-IR-bright QSOs. While the most luminous objects are the most readily observed, detections of more typical galaxies with modest SFRs are essential for completing the picture. ALMA will be revolutionary in terms of increasing the detection rate and pushing the sensitivity limit down to include such galaxies, however the limited FoV when observing at such high frequencies makes it difficult to use ALMA for studies of the large-scale structure traced out by molecular gas in galaxies. This article introduces a strategy for a systematic search for molecular gas during the EoR (z~7 and above), capitalizing on the fact that the J=1-0 transition of 12CO enters the upper bands of cm-wave instruments at high-z. The FoV advantage gained by observing at such frequencies, coupled with modern broadband correlators allows significant cosmological volumes to be probed on reasonable timescales. In this article we present an overview of our future observing programme which has been awarded 6,500 hours as one of the Large Survey Projects for MeerKAT, the forthcoming South African SKA pathfinder instrument. Its large FoV and correlator bandwidth, and high-sensitivity provide unprecedented survey speed for such work. An existing astrophysical simulation is coupled with instrumental considerations to demonstrate the feasibility of such observations and predict detection rates.

A consortium of universities has recently been formed with the goal of using the decommissioned telecommunications infrastructure at the Goonhilly Earth Station in Cornwall, UK, for astronomical purposes. One particular goal is the introduction of one or more of the ~30-metre parabolic antennas into the existing e-MERLIN radio interferometer. This article introduces this scheme and presents some simulations which quantify the improvements that would be brought to the e-MERLIN system. These include an approximate doubling of the spatial resolution of the array, an increase in its N-S extent with strong implications for imaging the most well-studied equatorial fields, accessible to ESO facilities including ALMA. It also increases the overlap between the e-MERLIN array and the European VLBI Network. We also discuss briefly some niche science areas in which an e-MERLIN array which included a receptor at Goonhilly would be potentially world-leading, in addition to enhancing the existing potential of e-MERLIN in its role as a Square Kilometer Array pathfinder instrument.

We use recent progress in simulating the production of magnetohydrodynamic jets around black holes to derive the cosmic spin history of the most massive black holes, with masses >~10^8 Msol. Assuming the jet efficiency depends on spin a, we can approximately reproduce the observed `radio loudness' of quasars and the local radio luminosity function. Using the X-ray luminosity function and the local mass function of supermassive black holes, SMBHs we can reproduce the individual radio luminosity functions of radio sources showing high- and low-excitation narrow emission lines. The data favour spin distributions that are bimodal, with one component around spin zero and the other close to maximal spin. In the low-excitation galaxies, the two components have similar amplitudes. For the high-excitation galaxies, the amplitude of the high-spin peak is typically much smaller than that of the low-spin peak. A bimodality should be seen in the radio loudness of quasars. We predict that the low-excitation galaxies are dominated by SMBHs with masses >~10^8 Msol, down to radio luminosity densities ~10^21 W Hz-1 sr-1 at 1.4~GHz. Our model is also able to predict the radio luminosity function at z=1, and predicts it to be dominated by high-excitation galaxies above luminosity densities >~10^26 W Hz-1 sr-1, in full agreement with the observations. From our parametrisation and using the best fitting jet efficiencies there is marginal evidence for evolution in spin: the mean spin increases slightly from <a>~0.25 at z=1 to <a>~0.35 at z=0, and the fraction of SMBHs with a>=0.5 increases from 0.16+-0.03 at z=1 to 0.24+-0.09 at z=0. Our results are in excellent agreement with the mean radiative efficiency of quasars, as well as recent cosmological simulations. We discuss the implications in terms of accretion and SMBH mergers, and galactic black holes (Abridged).

We use new mid-infrared (mid-IR) photometry from the Spitzer Space Telescope to study the relations between low-frequency radio luminosity density L_151MHz, mid-IR (12um rest-frame) luminosity L_12um, and optical-emission-line ([OII]) luminosity L_[OII], for a complete sample of z~1 radio galaxies from the 3CRR, 6CE, 6C*, 7CRS and TOOT00 surveys. The narrow redshift span of our sample (0.9<z<1.1) means that it is unbiased to evolutionary effects. We find evidence that these three quantities are positively correlated. The scaling between L_12um and L_[OII] is similar to that seen in other AGN samples, consistent with both L_12um and L_[OII] tracing accretion rate. We show that the positive correlation between L_12um and L_151MHz implies that there is a genuine lack of objects with low values of L_12um at high values of L_151MHz. Given that L_12um traces accretion rate, while L_151MHz traces jet power, this can be understood in terms of a minimum accretion rate being necessary to produce a given jet power. This implies that there is a maximum efficiency with which accreted energy can be chanelled into jet power and that this efficiency is of order unity.

Gravitational lensing is a powerful astrophysical and cosmological probe and is particularly valuable at submillimeter wavelengths for the study of the statistical and individual properties of dusty starforming galaxies. However the identification of gravitational lenses is often time-intensive, involving the sifting of large volumes of imaging or spectroscopic data to find few candidates. We used early data from the Herschel Astrophysical Terahertz Large Area Survey to demonstrate that wide-area submillimeter surveys can simply and easily detect strong gravitational lensing events, with close to 100% efficiency.

We present a parameterized model of the intra-cluster medium that is suitable for jointly analysing pointed observations of the Sunyaev-Zel'dovich (SZ) effect and X-ray emission in galaxy clusters. The model is based on assumptions of hydrostatic equilibrium, the Navarro, Frenk and White (NFW) model for the dark matter, and a softened power law profile for the gas entropy. We test this entropy-based model against high and low signal-to-noise mock observations of a relaxed and recently-merged cluster from N-body/hydrodynamic simulations, using Bayesian hyper-parameters to optimise the relative statistical weighting of the mock SZ and X-ray data. We find that it accurately reproduces both the global values of the cluster temperature, total mass and gas mass fraction (fgas), as well as the radial dependencies of these quantities outside of the core (r > kpc). For reference we also provide a comparison with results from the single isothermal beta model. We confirm previous results that the single isothermal beta model can result in significant biases in derived cluster properties.

We use the Herschel-ATLAS science demonstration data to investigate the star-formation properties of radio-selected galaxies in the GAMA-9h field as a function of radio luminosity and redshift. Radio selection at the lowest radio luminosities, as expected, selects mostly starburst galaxies. At higher radio luminosities, where the population is dominated by AGN, we find that some individual objects are associated with high far-infrared luminosities. However, the far-infrared properties of the radio-loud population are statistically indistinguishable from those of a comparison population of radio-quiet galaxies matched in redshift and K-band absolute magnitude. There is thus no evidence that the host galaxies of these largely low-luminosity (Fanaroff-Riley class I), and presumably low-excitation, AGN, as a population, have particularly unusual star-formation histories. Models in which the AGN activity in higher-luminosity, high-excitation radio galaxies is triggered by major mergers would predict a luminosity-dependent effect that is not seen in our data (which only span a limited range in radio luminosity) but which may well be detectable with the full Herschel-ATLAS dataset.

We present an analysis of observations of the doubly-lensed double hotspot in the giant radio galaxy J0816+5003 from MERLIN, MDM, WIYN, WHT, UKIRT and the VLA. The images of the two hotspot components span a factor of two in radius on one side of the lensing galaxy at impact parameters of less than 500pc. Hence we measure the slope of the lensing potential over a large range in radius, made possible by significant improvement in the accuracy of registration of the radio and optical frame and higher resolution imaging data than previously available. We also infer the lens and source redshifts to be 0.332 and > 1 respectively. Purely on the basis of lens modelling, and independently of stellar velocity dispersion measurements, we find the potential to be very close to isothermal.

The SKA at mid and low frequencies will be constructed in two distinct phases, the first being a subset of the second. This document defines the main scientific goals and baseline technical concept for the SKA Phase 1 (SKA_1). The major science goals for SKA_1 will be to study the history and role of neutral Hydrogen in the Universe from the dark ages to the present-day, and to employ pulsars as probes of fundamental physics. The baseline technical concept of SKA_1 will include a sparse aperture array operating at frequencies up to 450 MHz, and an array of dishes, initially operating at frequencies up to 3 GHz but capable of 10 GHz in terms of antenna surface accuracy. An associated Advanced Instrumentation Program (AIP) allows further development of new technologies currently under investigation. Construction will take place in 2016-2019 at a total capital cost of 350M\texteuro, including an element for contingency. The cost estimates of the SKA_1 telescope are now the subject of a more detailed and thorough costing exercise led by the SKA Project Development Office (SPDO). The 350 M\texteuro total for SKA_1 is a cost-constrained cap; an additional contingency is to reduce the overall scope of the project. The design of the SKA_1 is expected to evolve as the major cost estimates are refined, in particular the infrastructure costs at the two sites. The SKA_1 facility will represent a major step forward in terms of sensitivity, survey speed, image fidelity, temporal resolution and field-of-view. It will open up new areas of discovery space and demonstrate the science and technology underpinning the SKA Phase 2 (SKA_2).

We present deep spectroscopic data for a 24-object subsample of our full 41-object z~0.5 radio galaxy sample in order to investigate the evolution of the Fundamental Plane of radio galaxies. We find that the low-luminosity, FRI-type, radio galaxies in our sample are consistent with the local Fundamental Plane of radio galaxies defined by Bettoni et al. when corrected for simple passive evolution of their stellar populations. However, we find that the higher luminosity, FRII-type radio galaxies are inconsistent with the local Fundamental Plane if only passive evolution is considered, and find evidence for a rotation in the Fundamental Plane at z~0.5 when compared with the local relation. We show that neither passive evolution, nor a mass-dependent evolution in the mass-to-light ratio, nor an evolution in the size of the host galaxies can, by themselves, plausibly explain the observed tilt. However, we suggest that some combination of all three effects, with size evolution as the dominant factor, may be sufficient to explain the difference between the planes. We also find evidence for a correlation between host galaxy velocity dispersion and radio luminosity at the 97% significance level within our sub-sample, although further observations are required in order to determine whether this is different for the FRI and FRII radio sources. Assuming that the M_BH - sigma relation still holds at z~0.5, this implies that radio luminosity scales with black hole mass, in agreement with previous studies.

To investigate the poorly constrained sub-mm counts and spectral properties of blazars we searched for these in the Herschel-ATLAS (H-ATLAS) science demostration phase (SDP) survey catalog. We cross-matched 500$\mu$m sources brighter than 50 mJy with the FIRST radio catalogue. We found two blazars, both previously known. Our study is among the first blind blazar searches at sub-mm wavelengths, i.e., in the spectral regime where little is still known about the blazar SEDs, but where the synchrotron peak of the most luminous blazars is expected to occur. Our early results are consistent with educated extrapolations of lower frequency counts and question indications of substantial spectral curvature downwards and of spectral upturns at mm wavelengths. One of the two blazars is identified with a Fermi/LAT $\gamma$-ray source and a WMAP source. The physical parameters of the two blazars are briefly discussed.These observations demonstrate that the H-ATLAS survey will provide key information about the physics of blazars and their contribution to sub-mm counts.

We present the results of our investigation into the stellar populations of 24 radio galaxies at z~0.5 drawn from four complete, low-frequency selected radio surveys. We use the strength of the 4000A break as an indicator of recent star formation, and compare this with radio luminosity, optical spectral classification and morphological classification. We find evidence of different star formation histories for high- and low-luminosity radio sources; our group of low radio luminosity sources (typically FRI-type sources) has systematically older stellar populations than the higher radio luminosity group. Our sample is also fairly well divided by optical spectral classification. We find that galaxies classified as having low excitation spectra (LEGs) possess older stellar populations than high excitation line objects (HEGs), with the HEGs showing evidence for recent star formation. We also investigate the link between radio morphology, as used by Owen & Laing (1989), and the stellar populations. We find that there is a preference for the "fat-double" sources to have older stellar populations than the "classical double" sources, although this is also linked to these sources lying predominantly in the LEG and HEG categories respectively. These results are consistent with the hypothesis that HEGs are powered by accretion of cold gas, which could be supplied, for example, by recent mergers, secular instabilities, or filamentary cold flows. These processes could also trigger star formation in the host galaxy. The host galaxies of the LEGs do not show evidence for recent star formation and an influx of cold gas, and are consistent with being powered by the accretion of the hot phase of the inter-stellar medium.

We present a sky simulation of the atomic HI emission line and the first ten CO rotational emission lines of molecular gas in galaxies beyond the Milky Way. The simulated sky field has a comoving diameter of 500/h Mpc, hence the actual field-of-view depends on the (user-defined) maximal redshift zmax; e.g. for zmax=10, the field of view yields ~4x4 sqdeg. For all galaxies, we estimate the line fluxes, line profiles, and angular sizes of the HI and CO emission lines. The galaxy sample is complete for galaxies with cold hydrogen masses above 10^8 Msun. This sky simulation builds on a semi-analytic model of the cosmic evolution of galaxies in a Lambda-cold dark matter (LCDM) cosmology. The evolving CDM-distribution was adopted from the Millennium Simulation, an N-body CDM-simulation in a cubic box with a side length of 500/h Mpc. This side length limits the coherence scale of our sky simulation: it is long enough to allow the extraction of the baryon acoustic oscillations (BAOs) in the galaxy power spectrum, yet the position and amplitude of the first acoustic peak will be imperfectly defined. This sky simulation is a tangible aid to the design and operation of future telescopes, such the SKA, the LMT, and ALMA. The results presented in this paper have been restricted to a graphical representation of the simulated sky and fundamental dN/dz-analyzes for peak flux density limited and total flux limited surveys of HI and CO. A key prediction is that HI will be harder to detect at redshifts z>2 than predicted by a no-evolution model. The future verification or falsification of this prediction will allow us to qualify the semi-analytic models.

To exploit synergies between the Herschel Space Observatory and next generation radio facilities, we have extended the semi-empirical extragalactic radio continuum simulation of Wilman et al. (2008) to the mid- and far-infrared. Here we describe the assignment of infrared spectral energy distributions (SEDs) to the star-forming galaxies and active galactic nuclei, using Spitzer 24, 70 and 160 micron and SCUBA 850 micron survey results as the main constraints. Star-forming galaxies dominate the source counts, and a model in which their far-infrared-radio correlation and infrared SED assignment procedure are invariant with redshift underpredicts the observed 24 and 70 micron source counts. The 70 micron deficit can be eliminated if the star-forming galaxies undergo stronger luminosity evolution than originally assumed for the radio simulation, a requirement which may be partially ascribed to known non-linearity in the far-infrared--radio correlation at low luminosity if it evolves with redshift. At 24 micron, the shortfall is reduced if the star-forming galaxies develop SEDs with cooler dust and correspondingly stronger Polycyclic Aromatic Hydrocarbon (PAH) emission features with increasing redshift at a given far-infrared luminosity, but this trend may reverse at z>1 in order not to overproduce the sub-mm source counts. The resulting model compares favourably with recent BLAST results and we have extended the simulation database to aid the interpretation of Herschel surveys. Such comparisons may also facilitate further model refinement and revised predictions for the SKA and its precursors.

The current status of the HI simulation efforts is presented, in which a self consistent simulation path is described and basic equations to calculate array sensitivities are given. There is a summary of the SKA Design Study (SKADS) sky simulation and a method for implementing it into the array simulator is presented. A short overview of HI sensitivity requirements is discussed and expected results for a simulated HI survey are presented.

We present a range of new observations of the "broad absorption line radio galaxy" 4C +72.26 (z = 3.5), including sensitive rest-frame ultraviolet integral field spectroscopy using the Gemini/GMOS-N instrument and Subaru/CISCO K-band imaging and spectroscopy. We show that 4C +72.26 is a system of two vigorously star-forming galaxies superimposed along the line of sight separated by ~1300 +/- 200 km/s in velocity, with each demonstrating spectroscopically resolved absorption lines. The most active star-forming galaxy also hosts the accreting supermassive black hole which powers the extended radio source. We conclude that the star formation is unlikely to have been induced by a shock caused by the passage of the radio jet, and instead propose that a collision is a more probable trigger for the star formation. Despite the massive starburst, the UV-mid-infrared spectral energy distribution suggests that the pre-existing stellar population comprises ~10^12 Msolar of stellar mass, with the current burst only contributing a further ~2%, suggesting that 4C +72.26 has already assembled most of its final stellar mass.

We present the results of a search for galaxy clusters in Subaru-XMM Deep Field. We reach a depth for a total cluster flux in the 0.5-2 keV band of 2x10^{-15} ergs cm^{-2} s^{-1} over one of the widest XMM-Newton contiguous raster surveys, covering an area of 1.3 square degrees. Cluster candidates are identified through a wavelet detection of extended X-ray emission. The red sequence technique allows us to identify 57 cluster candidates. We report on the progress with the cluster spectroscopic follow-up and derive their properties based on the X-ray luminosity and cluster scaling relations. In addition, 3 sources are identified as X-ray counterparts of radio lobes, and in 3 further sources, X-ray counterpart of radio lobes provides a significant fraction of the total flux of the source. In the area covered by NIR data, our identification success rate achieves 86%. We detect a number of radio galaxies within our groups and for a luminosity-limited sample of radio galaxies we compute halo occupation statistics using a marked cluster mass function. We compare the cluster detection statistics in the SXDF with the predictions of concordance cosmology and current knowledge of the X-ray cluster properties, concluding that a reduction of concordance sigma_8 value by 5% is required in order to match the prediction of the model and the data. This conclusion still needs verification through the completion of cluster follow-up.

The S3-Tools are a set of Python-based routines and interfaces whose purpose is to provide user-friendly access to the SKA Simulated Skies (S3) set of simulations, an effort led by the University of Oxford in the framework of the European Union's SKADS program (http://www.skads-eu.org). The databases built from the S3 simulations are hosted by the Oxford e-Research Center (OeRC), and can be accessed through a web portal at http://s-cubed.physics.ox.ac.uk. This paper focuses on the practical steps involved to make radio images from the S3-SEX and S3-SAX simulations using the S3-Map tool and should be taken as a broad overview. For a more complete description, the interested reader should look up the user's guide. The output images can then be used as input to instrument simulators, e.g. to assess technical designs and observational strategies for the SKA and SKA pathfinders.

The Herschel ATLAS is the largest open-time key project that will be carried out on the Herschel Space Observatory. It will survey 510 square degrees of the extragalactic sky, four times larger than all the other Herschel surveys combined, in five far-infrared and submillimetre bands. We describe the survey, the complementary multi-wavelength datasets that will be combined with the Herschel data, and the six major science programmes we are undertaking. Using new models based on a previous submillimetre survey of galaxies, we present predictions of the properties of the ATLAS sources in other wavebands.

We present observations at 1.2 mm with MAMBO-II of a sample of z>~2 radio-intermediate obscured quasars, as well as CO observations of two sources with the Plateau de Bure Interferometer. Five out of 21 sources (24%) are detected at a significance of >=3sigma. Stacking all sources leads to a statistical detection of <S_1.2mm>= 0.96+-0.11 mJy and stacking only the non-detections also yields a statistical detection, with <S_1.2mm>= 0.51+-0.13 mJy. This corresponds to a typical far-infrared luminosity L_FIR~4x10^12 Lsol. If the far-infrared luminosity is powered entirely by star-formation, and not by AGN-heated dust, then the characteristic inferred star-formation rate is ~700 Msol yr-1. This far-infrared luminosity implies a dust mass of M_dust~3x10^8 Msol. We estimate that such large dust masses on kpc scales can plausibly cause the obscuration of the quasars. We present dust SEDs for our sample and derive a mean SED for our sample. This mean SED is not well fitted by clumpy torus models, unless additional extinction and far-infrared re-emission due to cool dust are included. There is a hint that the host galaxies of obscured quasars must have higher far-infrared luminosities and cool-dust masses and are therefore often found at an earlier evolutionary phase than those of unobscured quasars. For one source at z=2.767, we detect the CO(3-2) transition, with S_CO Delta nu=630+-50 mJy km s-1, corresponding to L_CO(3-2)= 3.2x10^7 Lsol, or L'_CO(3-2)=2.4x10^10 K km s-1 pc2. For another source at z=4.17, the lack of detection of the CO(4-3) line yields a limit of L'_CO(4-3)<1x10^10 K km s-1 pc2. Molecular gas masses, gas depletion timescales and gas-to-dust ratios are estimated (Abridged).

Infrared spectroscopy of the H-alpha emission lines of a sub-sample of 19 high-redshift (0.8 < z < 2.3) Molonglo quasars, selected at 408 MHz, is presented. These emission lines are fitted with composite models of broad and narrow emission, which include combinations of classical broad-line regions of fast-moving gas clouds lying outside the quasar nucleus, and/or a theoretical model of emission from an optically-thick, flattened, rotating accretion disk. All bar one of the nineteen sources are found to have emission consistent with the presence of an optically-emitting accretion disk, with the exception appearing to display complex emission including at least three broad components. Ten of the quasars have strong Bayesian evidence for broad-line emission arising from an accretion disk together with a standard broad-line region, selected in preference to a model with two simple broad lines. Thus the best explanation for the complexity required to fit the broad H-alpha lines in this sample is optical emission from an accretion disk in addition to a region of fast-moving clouds. We derive estimates of the angle between the rotation axis of the accretion disk and the line of sight. A weak correlation is found between the accretion disk angle and the logarithm of the low-frequency radio luminosity. This is direct, albeit tenuous, evidence for the receding torus model. Velocity shifts of the broad H-alpha components are analysed and the results found to be consistent with a two-component model comprising one single-peaked broad line emitted at the same redshift as the narrow lines, and emission from an accretion disk which appears to be preferentially redshifted with respect to the narrow lines for high-redshift sources and blueshifted relative to the narrow lines for low-redshift sources.

We present optical spectroscopy, near-infrared (mostly K-band) and radio (151-MHz and 1.4-GHz) imaging of the first complete region (TOOT00) of the TexOx-1000 (TOOT) redshift survey of radio sources. The 0.0015-sr (~ 5 deg^2) TOOT00 region is selected from pointed observations of the Cambridge Low-Frequency Survey Telescope at 151 MHz at a flux density limit of ~= 100 mJy, ~ 5-times fainter than the 7C Redshift Survey (7CRS), and contains 47 radio sources. We have obtained 40 spectroscopic redshifts (~ 85% completeness). Adding redshifts estimated for the 7 other cases yields a median redshift z_med ~ 1.25. We find a significant population of objects with FRI-like radio structures at radio luminosities above both the low-redshift FRI/II break and the break in the radio luminosity function. The redshift distribution and sub-populations of TOOT00 are broadly consistent with extrapolations from the 7CRS/6CE/3CRR datasets underlying the SKADS Simulated Skies Semi-Empirical Extragalactic Database, S^3-SEX.

We present the first results from the largest deep extragalactic millimetre-wavelength survey undertaken to date. These results are derived from maps covering over 0.7 deg^2, made at 1.1mm, using the AzTEC continuum camera mounted on the James Clerk Maxwell Telescope. The maps were made in the two fields originally targeted at 0.85mm with SCUBA in the SHADES project, namely the Lockman Hole East (mapped to a depth of 0.9-1.3 mJy rms) and the Subaru XMM Deep Field (1.0-1.7 mJy rms). The wealth of existing and forthcoming deep multi-frequency data in these two fields will allow the bright mm source population revealed by these images to be explored in detail in subsequent papers. Here we present the maps themselves, a catalogue of 114 high-significance sub-millimetre galaxy detections, and a thorough statistical analysis leading to the most robust determination to date of the 1.1mm source number counts. Through careful comparison, we find that both the COSMOS and GOODS North fields, also imaged with AzTEC, contain an excess of mm sources over the new 1.1mm source-count baseline established here. In particular, our new AzTEC/SHADES results indicate that very luminous high-redshift dust enshrouded starbursts (S_{1.1} > 3 mJy) are 25-50% less common than would have been inferred from these smaller surveys, thus highlighting the potential roles of cosmic variance and clustering in such measurements. We compare number count predictions from recent models of the evolving mm/sub-mm source population to these SMG surveys, which provide important constraints for the ongoing refinement of semi-analytic and hydrodynamical models of galaxy formation, and find that all recent models over-predict the number of bright sub-millimetre galaxies found in this survey.

We use two independent methods to forecast the dark energy measurements achievable by combining future galaxy redshift surveys based on the radio HI emission line with Cosmic Microwave Background (CMB) data from the {\sl Planck} satellite. In the first method we focus on the `standard ruler' provided by the baryon acoustic oscillation (BAO) length scale. In the second method we utilize additional information encoded in the galaxy power spectrum including galaxy bias from velocity-space distortions and the growth of cosmic structure. We find that a radio synthesis array with about 10 per cent of the collecting area of the Square Kilometre Array (SKA), equipped with a wide ($10-100 ~ {\rm deg}^2$) field-of-view, would have the capacity to perform a $20{,}000 ~ {\rm deg}^2$ redshift survey to a maximum redshift $z_{\rm max} \sim 0.8$ and thereby produce dark energy measurements that are competitive with surveys likely to be undertaken by optical telescopes around 2015. There would then be powerful arguments for adding collecting area to such a `Phase-1' SKA because of the square-law scaling of survey speed with telescope sensitivity for HI surveys, compared to the linear scaling for optical redshift surveys. The full SKA telescope should, by performing a $20{,}000 ~ {\rm deg}^2$ HI redshift survey to $z_{\rm max} \sim 2$ around 2020, yield an accurate measurement of cosmological parameters independent of CMB datasets. Combining CMB ({\sl Planck}) and galaxy power spectrum (SKA) measurements will drive errors in the dark energy equation-of-state parameter $w$ well below the 1 per cent level. The major systematic uncertainty in these forecasts is the lack of direct information about the mass function of high-redshift HI-emitting galaxies.

Galaxies in the early Universe were more compact and contained more molecular gas than today. In this paper, we revisit the relation between these empirical findings, and we quantitatively predict the cosmic evolution of the surface densities of atomic (HI) and molecular (H2) hydrogen in regular galaxies. Our method uses a pressure-based model for the H2/HI-ratio of the Interstellar Medium, applied to ~3*10^7 virtual galaxies in the Millennium Simulation. We predict that, on average, the HI-surface density of these galaxies saturates at Sigma_HI<10 Msun/pc^2 at all redshifts (z), while H2-surface densities evolve dramatically as Sigma_H2(1+z)^2.4. This scaling is dominated by a (1+z)^2 surface brightness scaling originating from the (1+z)^-1 size scaling of galaxies at high z. Current measurements of Sigma_H2 at high z, derived from CO-observations, tend to have even higher values, which can be quantitatively explained by a selection bias towards merging systems. However, despite the consistency between our high-z predictions and the sparse empirical data, we emphasize that the empirical data potentially suffer from serious selection biases and that the semi-analytic models remain in many regards uncertain. As a case study, we investigate the cosmic evolution of simulated galaxies, which resemble the Milky Way at z=0. We explicitly predict their HI- and H2-distribution at z=1.5, corresponding to the CO-detected galaxy BzK-21000, and at z=3, corresponding to the primary science goal of the Atacama Large Millimeter/submillimeter Array (ALMA).

We use a pressure-based model for splitting cold hydrogen into its atomic (HI) and molecular (H2)components to tackle the co-evolution of HI, H2, and star formation rates (SFR) in ~3e7 simulated galaxies in the Millennium simulation. The main prediction is that galaxies contained similar amounts of HI at redshift z=1-5 than today, but substantially more H2, in quantitative agreement with the strong molecular line emission already detected in a few high redshift galaxies and approximately consistent with inferences from studies of the damped Lyman-alpha absorbers seen in the spectra of quasars. The cosmic H2/HI-ratio is predicted to evolve monotonically as Omega(H2)/Omega(HI) (1+z)^1.6. This decline of the H2/HI-ratio as a function of cosmic time is driven by the growth of galactic disks and the progressive reduction of the mean cold gas pressure. Finally, a comparison between the evolutions of HI, H2, and SFRs reveals two distinct cosmic epochs of star formation: an early epoch (z>3), driven by the evolution of Omega(HI+H2), and a late epoch (z<3), driven by the evolution of Omega(H2)/Omega(HI).

We present a simulation of the cosmic evolution of the atomic and molecular phases of the cold hydrogen gas in about 3e7 galaxies, obtained by post-processing the virtual galaxy catalog produced by (De Lucia et al. 2007) on the Millennium Simulation of cosmic structure (Springel et al. 2005). Our method uses a set of physical prescriptions to assign neutral atomic hydrogen (HI) and molecular hydrogen (H2) to galaxies, based on their total cold gas masses and a few additional galaxy properties. These prescriptions are specially designed for large cosmological simulations, where, given current computational limitations, individual galaxies can only be represented by simplistic model-objects with a few global properties. Our recipes allow us to (i) split total cold gas masses between HI, H2, and Helium, (ii) assign realistic sizes to both the HI- and H2-disks, and (iii) evaluate the corresponding velocity profiles and shapes of the characteristic radio emission lines. The results presented in this paper include the local HI- and H2-mass functions, the CO-luminosity function, the cold gas mass--diameter relation, and the Tully-Fisher relation (TFR), which all match recent observational data from the local Universe. We also present high-redshift predictions of cold gas diameters and the TFR, both of which appear to evolve markedly with redshift.

We revisit the mass ratio Rmol between molecular hydrogen (H2) and atomic hydrogen (HI) in different galaxies from a phenomenological and theoretical viewpoint. First, the local H2-mass function (MF) is estimated from the local CO-luminosity function (LF) of the FCRAO Extragalactic CO-Survey, adopting a variable CO-to-H2 conversion fitted to nearby observations. This implies an average H2-density Omega_H2=(6.9+-2.7) 10^5/h and Omega_H2/Omega_HI=0.26+-0.11 in the local Universe. Second, we investigate the correlations between Rmol and global galaxy properties in a sample of 245 local galaxies. Based on these correlations we introduce four phenomenological models for Rmol, which we apply to estimate H2-masses for each HI-galaxy in the HIPASS catalogue. The resulting H2-MFs (one for each model for Rmol) are compared to the reference H2-MF derived from the CO-LF, thus allowing us to determine the Bayesian evidence of each model and to identify a clear best model, in which, for spiral galaxies, Rmol negatively correlates with both galaxy Hubble type and total gas mass. Third, we derive a theoretical model for Rmol for regular galaxies based on an expression for their axially symmetric pressure profile dictating the degree of molecularization. This model is quantitatively similar to the best phenomenological one at redshift z=0, and hence represents a consistent generalization while providing a physical explanation for the dependence of Rmol on global galaxy properties. Applying the best phenomenological model for Rmol to the HIPASS sample, we derive the first integral cold gas-MF (HI+H2+helium) of the local Universe.

In this letter we report the discovery of a z=4.88 radio galaxy discovered with a new technique which does not rely on pre-selection of a sample based on radio properties such as steep-spectral index or small angular size. This radio galaxy was discovered in the Elais-N2 field and has a spectral index of alpha = 0.75, i.e. not ultra-steep spectrum. It also has a luminosity consistent with being drawn from the break of the radio luminosity function and can therefore be considered as a typical radio galaxy. Using the Spitzer-SWIRE data over this field we find that the host galaxy is consistent with being similarly massive to the lower redshift powerful radio galaxies (~1-3L*). We note however, that at z=4.88 the H-alpha line is redshifted into the IRAC 3.6micron filter and some of the flux in this band may be due to this rather than stellar continuum emission. The discovery of such a distant radio source from our initial spectroscopic observations demonstrate the promise of our survey for finding the most distant radio sources.

The radio properties of 11 obscured `radio-intermediate' quasars at redshifts z>~2 have been investigated using the European Very-Long-Baseline-Interferometry Network (EVN) at 1.66 GHz. A sensitivity of ~25 micro Jy per 14x17 mas2 beam was achieved, and in 7 out of 11 sources unresolved radio emission was securely detected. The detected radio emission of each source accounts for ~30-100 % of the total source flux density. The physical extent of this emission is ~<150 pc, and the derived properties indicate that this emission originates from an active galactic nucleus (AGN). The missing flux density is difficult to account for by star-formation alone, so radio components associated with jets of physical size >~150 pc, and ~< 40 kpc are likely to be present in most of the sources. Amongst the observed sample steep, flat, gigahertz-peaked and compact-steep spectrum sources are all present. Hence, as well as extended and compact jets, examples of beamed jets are also inferred, suggesting that in these sources, the obscuration must be due to dust in the host galaxy, rather than the torus invoked by the unified schemes. Comparing the total to core (~< 150 pc) radio luminosities of this sample with different types of AGN suggests that this sample of z >~2 radio-intermediate obscured quasars shows radio properties that are more similar to those of the high-radio-luminosity end of the low-redshift radio-quiet quasar population than those of FR I radio galaxies. This conclusion may reflect intrinsic differences, but could be strongly influenced by the increasing effect of inverse-Compton cooling of extended radio jets at high redshift.

We outline the case for a comprehensive wide and deep survey ultimately targeted at obtaining 21-cm HI line emission spectroscopic observations of more than a billion galaxies to redshift z=1.5 and greater over half the sky. This survey provides a database of galaxy redshifts, HI gas masses, and galaxy rotation curves that would enable a wide range of science, including fundamental cosmology and studies of Dark Energy. This science requires the next generation of radio arrays, which are being designed under the umbrella of the Square Kilometer Array (SKA) project. We present a science roadmap, extending to 2020 and beyond, that would enable this ambitious survey. We also place this survey in the context of other multi-wavelength surveys.

We discuss the optical properties, X-ray detections, and Active Galactic Nucleus (AGN) populations of four clusters at z ~ 1 in the Subaru-XMM Deep Field (SXDF). The velocity distribution and plausible extended X-ray detections are examined, as well as the number of X-ray point sources and radio sources associated with the clusters. We find that the two clusters that appear virialised and have an extended X-ray detection contain few, if any, AGN, whereas the two pre-virialised clusters have a large AGN population. This constitutes evidence that the AGN fraction in clusters is linked to the clusters' evolutionary stage. The number of X-ray AGN in the pre-virialised clusters is consistent with an overdensity of factor ~ 200; the radio AGN appear to be clustered with a factor of three to six higher. The median K-band luminosities of L_K = 1.7 +/- 0.7 L* for the X-ray sources and L_K = 2.3 +/- 0.1 L* for the radio sources support the theory that these AGN are triggered by galaxy interaction and merging events in sub-groups with low internal velocity distributions, which make up the cluster environment in a pre-virialisation evolutionary stage.

We report here on the discovery of a red quasar, J004929.4+351025.7 at a redshift of z = 2.48, situated within a large Lya emission-line halo. The radio spectral energy distribution implies that the radio jets were triggered < 10^4 years prior to the time at which the object is observed, suggesting that the jet triggering of the active galactic nucleus is recent. The loosely biconical structure of the emission-line halo suggests that it is ionised by photons emitted by the central quasar nucleus and that the central nucleus is obscured by a dusty torus with Av ~ 3.0. The large spatial extent of the Lya halo relative to the radio emission means this could only have occurred if the radio jets emerged from an already established highly-accreting black hole. This suggests that the radio-jet triggering is delayed with respect to the onset of accretion activity on to the central supermassive black hole.

We have developed a semi-empirical simulation of the extragalactic radio continuum sky suitable for aiding the design of next generation radio interferometers such as the Square Kilometre Array (SKA). The emphasis is on modelling the large-scale cosmological distribution of radio sources rather than the internal details of individual galaxies. Here we provide a description of the simulation to accompany the online release of a catalogue of 320 million simulated radio sources. The simulation covers 20x20 deg^2 - a plausible upper limit to the instantaneous field of view attainable with future (e.g. SKA) aperture array technologies - out to redshift z=20, and down to flux density limits of 10 nJy at 151, 610 MHz, 1.4, 4.86 and 18 GHz. Five distinct source types are included: radio-quiet AGN, radio-loud AGN of the FRI and FRII structural classes, and star-forming galaxies, the latter split into populations of quiescent and starbursting galaxies. In our semi-empirical approach, the simulated sources are drawn from observed (or extrapolated) luminosity functions and grafted onto an underlying dark matter density field with biases which reflect their measured large-scale clustering. A numerical Press-Schechter-style filtering of the density field is used to identify and populate clusters of galaxies. Radio source structures are built from point source and elliptical sub-components, and for FRI and FRII sources an orientation-based unification and beaming model is used to partition flux between the core and extended lobes and hotspots. The simulation output can be post-processed to achieve more complete agreement with observational data in the years ahead, with the aim of using these 'idealised skies' in telescope simulators to optimise the design of the SKA itself (abridged).

We study the 37 brightest radio sources in the Subaru/XMM-Newton Deep Field (SXDF). We have spectroscopic redshifts for 24 of 37 objects and photometric redshifts for the remainder, yielding a median redshift z_med for the whole sample of z_med ~= 1.1 and a median radio luminosity close to the `FRI/FRII' luminosity divide. Using mid-IR (Spitzer MIPS 24 um) data we expect to trace nuclear accretion activity, even if it is obscured at optical wavelengths, unless the obscuring column is extreme. Our results suggest that above the FRI/FRII radio luminosity break most of the radio sources are associated with objects that have excess mid-IR emission, only some of which are broad-line objects, although there is one clear low-accretion-rate object with an FRI radio structure. For extended steep-spectrum radio sources, the fraction of objects with mid-IR excess drops dramatically below the FRI/FRII luminosity break, although there exists at least one high-accretion-rate `radio-quiet' QSO. We have therefore shown that the strong link between radio luminosity (or radio structure) and accretion properties, well known at z ~ 0.1, persists to z ~ 1. Investigation of mid-IR and blue excesses shows that they are correlated as predicted by a model in which, when significant accretion exists, a torus of dust absorbs ~30% of the light, and the dust above and below the torus scatters >~1% of the light.

We present the optical-to-submm spectral energy distributions for 33 radio & mid-IR identified submillimetre galaxies discovered via the SHADES 850 micron SCUBA imaging in the Subaru-XMM Deep Field (SXDF). Optical data for the sources comes from the Subaru-XMM Deep Field (SXDF) and mid- and far-IR fluxes from SWIRE. We obtain photometric redshift estimates for our sources using optical and IRAC 3.6 and 4.5 micron fluxes. We then fit spectral energy distribution (SED) templates to the longer wavelength data to determine the nature of the far-IR emission that dominates the bolometric luminosity of these sources. The infrared template fits are also used to resolve ambiguous identifications and cases of redshift aliasing. The redshift distribution obtained broadly matches previous results for submm sources and on the SHADES SXDF field. Our template fitting finds that AGN, present in about 10% of sources, do not contribute significantly to their bolometric luminosity. Dust heating by starbursts, with either Arp220 or M82 type SEDs, appears to be responsible for the luminosity in most sources (23/33 are fitted by Arp220 templates, 2/33 by the warmer M82 templates). 8/33 sources, in contrast, are fit by a cooler cirrus dust template, suggesting that cold dust has a role in some of these highly luminous objects. Three of our sources appear to have multiple identifications or components at the same redshift, but we find no statistical evidence that close associations are common among our SHADES sources. Examination of rest-frame K -band luminosity suggests that 'downsizing' is underway in the submm galaxy population, with lower redshift systems lying in lower mass host galaxies. Of our 33 identifications six are found to be of lower reliability but their exclusion would not significantly alter our conclusions.

We have analysed the 24um properties of a radio-selected sample in the Subaru-XMM/Newton Deep Field in order to explore the behaviour of the FIR/radio relation at high redshifts. Statistically, the correlation is described by q24, the ratio between the observed flux densities at 24um and 1.4GHz, respectively. Using 24um data results in considerably more scatter in the correlation than previous work using data at 60-70um. Nevertheless, we do observe a steady correlation as a function of redshift, up to z~3.5, suggesting its validity back to primeval times. We find q24 = 0.30 +/- 0.56 for the observed and q24 = 0.71 +/- 0.47 for the k-corrected radio sample, based on sources with 300uJy < S(1.4GHz) < 3.2mJy and 24um detections. A suitable k-correction given by a M82-like mid-IR template suggests no extreme silicate absorption in the bulk of our radio sample. Using thresholds in q24 to identify radio-excess sources, we have been able to characterise the transition from radio-loud AGN to star-forming galaxies and radio-quiet AGN at faint (<1mJy) radio flux densities. Our results are in broad agreement with previous studies which show a dominant radio-loud AGN population at >1mJy. The rest-frame U-B colours of the expected radio-excess population have redder distribution than those that follow the correlation. This is therefore a promising way to select obscured Type-2 AGN, with a radio loud nature, missed by deep X-ray observations. Spectroscopic follow-up of these sources is required to fully test this method.

We discuss the equatorial imaging benefits that arise from the addition of the 25-metre dish at Chilbolton to the e-MERLIN array. Its inclusion considerably enhances the capabilities of e-MERLIN on and below the equator. This will become particularly important in the era of ALMA and other upcoming southern hemisphere facilities. We present simulated observations of point sources in the equatorial region of the sky which is the target area for many existing sky surveys. We find that the additional baselines created by the inclusion of the Chilbolton dish favourably adjust the beam shape of e-MERLIN to a more compact and circular shape, with significantly reduced sidelobe structure. Putting aside the benefits of increased collecting area, the modified beam shape has implications for more rapidly reaching a given completeness limit for equatorial surveys.

We present mid-infrared observations of 18 sources from a sample of 21 z~2 radio-intermediate obscured quasars. The mid-infrared spectra of the sources are continuum dominated, and 12 sources show deep silicate absorption with tau_9.7~1-2. Combining mid-infrared and optical spectra, we achieve 86% spectroscopic completeness which allows us to confirm that most (63 +14 -22 %) z~2 radio-intermediate quasars are obscured. The new spectra also prove that many high-redshift type-2 quasars do not show any rest-frame ultraviolet emission lines. From the 18 individual mid-infrared spectra, we classify most of the sources into three subsamples: those with hints of the 7.7 and 6.2 micron polyaromatic hydrocarbons (3/18 sources show PAHs, subsample A), those with an excess of emission around 8 micron but no hint of the 6.2 micron PAH (7/18 cases, subsample B) and pure-continuum sources with no visible excess (4/18 sources, subsample C). The remaining 4/18 sources have spectra that are featureless or too noisy for any features to be visible. In subsample A, averaging the spectra leads to a statistical detection of both 6.2 and 7.7 micron PAHs over the continuum, with the strength of the 7.7 micron PAH comparable to that of submillimetre-selected galaxies (SMGs) at similar redshifts. These sources are in a phase of coeval growth of a supermassive black hole and a host galaxy.

We study the 37 brightest radio sources in the Subaru/XMM-Newton Deep Field (SXDF). Using mid-IR (Spitzer MIPS 24 micron) data we expect to trace nuclear accretion activity, even if it is obscured at optical wavelengths, unless the obscuring column is extreme. Our results suggest that above the `FRI/FRII' radio luminosity break most of the radio sources are associated with objects that have excess mid-IR emission, only some of which are broad-line objects, although there is one clear low-accretion-rate FRI. The fraction of objects with mid-IR excess drops dramatically below the FRI/FRII break, although there exists at least one high-accretion-rate QSO. Investigation of mid-IR and blue excesses shows that they are correlated as predicted by a model in which a torus of dust absorbs ~30% of the light, and the dust above and below the torus scatters >~1% of the light.

We have obtained deep optical spectroscopic data of the highest-redshift cluster candidate (z ~ 1.4, CVB13) selected by Van Breukelen et al. (2006) in a photometric optical/infrared catalogue of the Subaru XMM-Newton Deep Field. The data, which comprise 104 targeted galaxies, were taken with the DEep Imaging Multi-Object Spectrograph (DEIMOS) on the Keck 2 telescope and yielded 31 secure redshifts in the range 1.25 < z < 1.54 within a 7' x 4' field centred on CVB13. Instead of one massive cluster at z = 1.4, we find evidence for three projected structures at z = 1.40, z = 1.45, and z = 1.48. The most statistically robust of these structures, at z = 1.454, has six spectroscopically confirmed galaxies. Its total mass is estimated at > 10^14 M_sun and it may therefore be termed a cluster. There is an X-ray source at the cluster position which is marginally spatially resolved but whose X-ray spectrum is too hard to be thermal cluster emission. Its origin could be the summed X-ray emission from active galaxies in, and projected onto, the cluster. Serendipitously we have discovered a cluster at z = 1.28 with a mass of > 10^14 M_sun at the same position on the sky, comprising six spectroscopically confirmed cluster galaxies and at least one additional radio source. The selection of CVB13 for the cluster catalogue was evidently aided by the superposition of other, presumably lower-mass, structures, whereas the single cluster at z = 1.28 contained too few galaxies to be isolated by the same algorithm. Given the complicated nature of such structures, caution must be employed when measuring the mass function of putative high-redshift clusters with photometric techniques alone.

We present deep, multi-VLA-configuration radio images for a set of 18 quasars, having redshifts between 0.36 and 2.5, from the 7C quasar survey. Approximately one quarter of these quasars have FRI-type twin-jet structures and the remainder are a broad range of wide angle tail, fat double, classical double, core-jet and hybrid sources. These images demonstrate that FRI quasars are prevalent in the universe, rather than non-existent as had been suggested in the literature prior to the serendipitous discovery of the first FRI quasar a few years ago, the optically powerful "radio quiet" quasar E1821+643. Some of the FRI quasars have radio luminosities exceeding the traditional FRI / FRII break luminosity, however we find no evidence for FRII quasars with luminosities significantly below the break. We consider whether the existence of such high luminosity FRI structures is due to the increasingly inhomogeneous environments in the higher redshift universe.

We present multi-object spectroscopy of galaxies in the immediate (Mpc-scale) environments of four low-power (L_1.4 GHz < 10^25 W/Hz) radio galaxies at z~0.5, selected from the Subaru/XMM-Newton Deep Field. We use the spectra to calculate velocity dispersions and central redshifts of the groups the radio galaxies inhabit, and combined with XMM-Newton (0.3-10 keV) X-ray observations investigate the L_X--sigma_v and T_X--sigma_v scaling relationships. All the radio galaxies reside in moderately rich groups -- intermediate environments between poor groups and rich clusters, with remarkably similar X-ray properties. We concentrate our discussion on our best statistical example that we interpret as a low-power (FRI) source triggered within a sub-group, which in turn is interacting with a nearby group of galaxies, containing the bulk of the X-ray emission for the system -- a basic scenario which can be compared to more powerful radio sources at both high (z>4) and low (z<0.1) redshifts. This suggests that galaxy-galaxy interactions triggered by group mergers may play an important role in the life-cycle of radio galaxies at all epochs and luminosities.

We analyse the first publicly released deep field of the UKIDSS Deep eXtragalactic Survey (DXS) to identify candidate galaxy over-densities at z~1 across ~1 sq. degree in the ELAIS-N1 field. Using I-K, J-K and K-3.6um colours we identify and spectroscopically follow-up five candidate structures with Gemini/GMOS and confirm they are all true over-densities with between five and nineteen members each. Surprisingly, all five structures lie in a narrow redshift range at z=0.89+/-0.01, although they are spread across 30Mpc on the sky. We also find a more distant over-density at z=1.09 in one of the spectroscopic survey regions. These five over-dense regions lying in a narrow redshift range indicate the presence of a supercluster in this field and by comparing with mock cluster catalogs from N-body simulations we discuss the likely properties of this structure. Overall, we show that the properties of this supercluster are similar to the well-studied Shapley and Hercules superclusters at lower redshift.

We present the redshift distribution of the SHADES galaxy population based on the rest-frame radio-mm-FIR colours of 120 robustly detected 850um sources in the Lockman Hole East (LH) and Subaru XMM-Newton Deep Field (SXDF). The redshift distribution derived from the full SED information is shown to be narrower than that determined from the radio-submm spectral index, as more photometric bands contribute to a higher redshift accuracy. The redshift distribution of sources derived from at least two photometric bands peaks at z ~ 2.4 and has a near-Gaussian distribution, with 50 per cent (interquartile range) of sources at z=1.8-3.1. We find a statistically-significant difference between the measured redshift distributions in the two fields; the SXDF peaking at a slightly lower redshift (median z ~ 2.2) than the LH (median z ~ 2.7), which we attribute to the noise-properties of the radio observations. We demonstrate however that there could also be field-to-field variations that are consistent with the measured differences in the redshift distributions, and hence, that the incomplete area observed by SHADES with SCUBA, despite being the largest sub-mm survey to date, may still be too small to fully characterize the bright sub-mm galaxy population. Finally we present a brief comparison with the predicted, or assumed, redshift distributions of sub-mm galaxy formation and evolution models, and we derive the contribution of these SHADES sources and the general sub-mm galaxy population to the star formation-rate density at different epochs.

Determining an accurate position for a submm galaxy (SMG) is the crucial step that enables us to move from the basic properties of an SMG sample - source counts and 2-D clustering - to an assessment of their detailed, multi-wavelength properties, their contribution to the history of cosmic star formation and their links with present-day galaxy populations. In this paper, we identify robust radio and/or IR counterparts, and hence accurate positions, for over two thirds of the SCUBA HAlf-Degree Extragalactic Survey (SHADES) Source Catalogue, presenting optical, 24-um and radio images of each SMG. Observed trends in identification rate have given no strong rationale for pruning the sample. Uncertainties in submm position are found to be consistent with theoretical expectations, with no evidence for significant additional sources of error. Employing the submm/radio redshift indicator, via a parameterisation appropriate for radio-identified SMGs with spectroscopic redshifts, yields a median redshift of 2.8 for the radio-identified subset of SHADES, somewhat higher than the median spectroscopic redshift. We present a diagnostic colour-colour plot, exploiting Spitzer photometry, in which we identify regions commensurate with SMGs at very high redshift. Finally, we find that significantly more SMGs have multiple robust counterparts than would be expected by chance, indicative of physical associations. These multiple systems are most common amongst the brightest SMGs and are typically separated by 2-6", or 15-50/(sin i) kpc at z ~ 2, consistent with early bursts seen in merger simulations.

Using mid-infrared and radio selection criteria, we pre-select a sample of candidate high-redshift type-2 quasars in the Subaru XMM-Newton Deep Field (SXDF). To filter out starburst contaminants, we use a bayesian method to fit the spectral energy distributions (SEDs) between 24-microns and B-band, obtain photometric redshifts, and identify the best candidates for high-z type-2 quasars. This leaves us with 12 z_phot >= 1.7 type-2 quasar candidates in an area ~0.8 deg^2, of which only two have secure X-ray detections. The two detected sources have estimated column densities N_H~2 & 3x10^27 m^-2, i.e. heavily obscured but Compton-thin quasars. Given the large bolometric luminosities and redshifts of the undetected objects, the lack of X-ray detections suggests extreme absorbing columns N_H >= 10^28 m^-2 are typical. We have found evidence for a population of ``Compton-thick'' high-redshift type-2 quasars, at least comparable to, and probably larger than the type-1 quasar population, although spectroscopic confirmation of their AGN nature is important.

Current measurements of the large-scale structure of galaxies are able to place an $\sim 0.5 ~ \rm eV$ upper limit on the absolute mass scale of neutrinos. An order-of-magnitude improvement in raw sensitivity, together with an insensitivity to systematic effects, is needed to reach the lowest value allowed by particle physics experiments. We consider the prospects of determining both the neutrino mass scale and the number of of massive neutrinos using future redshift surveys, specifically those undertaken with the Square Kilometre Array (SKA), with and without additional constraints from the upcoming Planck CMB experiment. If the sum of the neutrino masses $\sum m_i \gtsimeq 0.25 ~ \rm eV$ then the imprint of neutrinos on large-scale structure (LSS) should be enough, on its own, to establish the neutrino mass scale and, considered alongside CMB constraints, it will also determine the number of massive neutrinos $N_{\nu}$, and hence the mass hierarchy. If $\sum m_i \sim 0.05 ~ \rm eV$, at the bottom end of the allowed range, then a combination of LSS, CMB and particle physics constraints should be able to determine $\sum m_i$, $N_{\nu}$ and the hierarchy. If $\sum m_i$ is in the specific range $0.1-0.25 ~ \rm eV$, then a combination of LSS, CMB and particle physics experiments should determine $\sum m_i$, but not $N_{\nu}$ or the hierarchy. Once an SKA-like LSS survey is available there are good prospects of obtaining a full understanding of the conventional neutrino sector, and a chance of finding evidence for sterile neutrinos.

We use the 6C** sample to investigate the co-moving space density of powerful, steep-spectrum radio sources. This sample, consisting of 68 objects, has virtually complete K-band photometry and spectroscopic redshifts for 32 per cent of the sources. In order to find its complete redshift distribution, we develop a method of redshift estimation based on the K-z diagram of the 3CRR, 6CE, 6C* and 7CRS radio galaxies. Based on this method, we derive redshift probability density functions for all the optically identified sources in the 6C** sample. Using a combination of spectroscopic and estimated redshifts, we select the most radio luminous sources in the sample. Their redshift distribution is then compared with the predictions of the radio luminosity function of Jarvis et al. We find that, within the uncertainties associated with the estimation method, the data are consistent with a constant co-moving space density of steep-spectrum radio sources beyond z > 2.5, and rule out a steep decline.

Using mid-infrared and radio criteria, we select a sample of candidate z~2 obscured quasars. Optical spectroscopy confirms about half of these as type-2 quasars, and modelling the population suggests 50-80% of the quasars are obscured. We find some flat radio spectrum type-2 quasars, and tentative evidence for obscuration unrelated to the torus. Using a similar sample, we also find evidence for a significant fraction of Compton-thick quasars.

3C 318 is a z=1.574 radio-loud quasar. The small physical size of its radio jets indicate that these jets were triggered relatively recently. In addition to the ultraviolet continuum emission being reddened by dust, detections with IRAS and SCUBA show it to have an exceptionally high far-infrared luminosity. We present CO(2-1) observations of 3C 318 made with the IRAM Plateau de Bure Interferometer. We detect CO(2-1) emission with a FWHM=200 km/s at a signal-to-noise ratio of 5.4. There is evidence for positional (~ 20 kpc) and velocity (~ -400 km/s) offsets between the molecular gas and the quasar which may be due to the quasar experiencing a major merger. The mass of molecular gas inferred from our observations is M_{H_2}=(3.0 +/- 0.6) x 10^{10} M_sun. This molecular gas mass is comparable to that in sub-mm-selected galaxies at similar redshifts. The large molecular gas mass is consistent with the primary source of heating for the cool dust in this quasar to be massive star formation with a star formation rate of 1700 M_sun/yr and a gas depletion timescale of 20 Myr. Our observations support the idea that star formation episodes and jet triggering can be synchronised.

We present basic observational data on the 6C** sample. This is a new sample of radio sources drawn from the 151 MHz 6C survey, which was filtered with radio criteria chosen to optimize the chances of finding radio galaxies at z > 4. The filtering criteria are a steep-spectral index and a small angular size. The final sample consists of 68 sources from a region of sky covering 0.421 sr. We present VLA radio maps, and the results of K-band imaging and optical spectroscopy. Near-infrared counterparts are identified for 66 of the 68 sources, down to a 3-sigma limiting magnitude of K ~ 22 mag in a 3-arcsec aperture. Eight of these identifications are spatially compact, implying an unresolved nuclear source. The K-magnitude distribution peaks at a median K=18.7 mag, and is found to be statistically indistinguishable from that of the similarly selected 6C* sample, implying that the redshift distribution could extend to z > 4. Redshifts determined from spectroscopy are available for 22 (32 per cent) of the sources, over the range of 0.2 < z < 3.3 . We measure 15 of these, whereas the other 7 were previously known. Six sources are at z > 2.5. Four sources show broad emission lines in their spectra and are classified as quasars. Three of these show also an unresolved K-band identification. Eleven sources fail to show any distinctive emission and/or absorption features in their spectra. We suggest that these could be (i) in the so-called `redshift desert' region of 1.2 < z < 1.8, or (ii) at a greater redshift, but feature weak emission line spectra.

We study the 36 brightest radio sources in the Subaru/XMM-Newton Deep Field SXDF). Using MIPS 24 micron data from Spitzer we expect to trace accretion activity, even if it is hidden at optical wavelengths, unless the obscuring column is extreme. Our results suggest that in the decade or so below the break in the radio luminosity function that at least half, and potentially nearly all, radio sources are associated with accreting quasar-like objects. This is not true at lower radio luminosities where the quasar-like fraction approaches zero once compact sources are excluded.

We describe deep radio imaging at 1.4-GHz of the 1.3 square degree Subaru/XMM-Newton Deep Field (SXDF), made with the Very Large Array in B and C configurations. We present a radio map of the entire field, and a catalogue of 505 sources covering 0.8 square degrees to a peak flux density limit of 100 microJy. Robust optical identifications are provided for 90% of the sources, and suggested IDs are presented for all but 14 (of which 7 are optically blank, and 7 are close to bright contaminating objects). We show that the optical properties of the radio sources do not change with flux density, suggesting that AGNs continue to contribute significantly at faint flux densities. We test this assertion by cross-correlating our radio catalogue with the X-ray source catalogue and conclude that radio-quiet AGNs become a significant population at flux densities below 300 microJy, and may dominate the population responsible for the flattening of the radio source counts if a significant fraction of them are Compton-thick.

We present multi-frequency radio observations of a sample of z~2 obscured (type-2) quasars in the Spitzer extragalactic First Look Survey area. We combine the public data at 1.4 GHz, used in the selection of these sources, with new observations at 610 MHz (GMRT) and at 4.9 GHz (VLA). We find the sample includes sources with steep, flat and gigahertz-peaked spectra. There are no strong correlations between the presence or absence of emission lines in the optical spectra and the radio spectral properties of the sample. However, there are no secure flat-spectrum type-2 quasars with narrow emission lines which would be problematic for unified schemes. Most of the population have straight radio spectra with spectral index alpha~1 as is expected for developed, potentially FRI-like, jets in which continous injection of relativistic electrons is accompanied by inverse-Compton losses against the cosmic microwave background.

We present the maps, source catalogue and number counts of the largest, most complete and unbiased extragalactic submillimetre survey ever undertaken: the 850-micron SCUBA HAlf Degree Extragalactic Survey (SHADES). Using the Submillimetre Common-User Bolometer Array (SCUBA) on the James Clerk Maxwell Telescope (JCMT), SHADES mapped two separate regions of sky: the Subaru/XMM-Newton Deep Field (SXDF) and the Lockman Hole East (LH). These SCUBA maps cover 720 square arcmin with an RMS noise level of about 2 mJy and have uncovered >100 submillimetre galaxies. In order to ensure the utmost robustness of the resulting source catalogue, data reduction was independently carried out by four sub-groups drawn within the SHADES team, providing an unprecedented degree of reliability with respect to other SCUBA catalogues available from the literature. Individual source lists from the four groups were combined to produce a robust 120-object SHADES catalogue; an invaluable resource for follow-up campaigns aiming to study the properties of a complete and consistent sample of submillimetre galaxies. We present differential and integral source number counts of submillimetre galaxies and find that the differential counts are better fit with a broken power-law or a Schechter function than with a single power-law; the SHADES data alone significantly show that a break is required at several mJy, although the precise position of the break is not well constrained. We also find that an 850-micron survey complete down to 2 mJy would resolve 20-30 per cent of the Far-IR background into point sources. [abridged]

We present the first cluster catalogue extracted from the UKIRT Infrared Deep Sky Survey Early Data Release. The catalogue is created using UKIDSS Ultra Deep Survey infrared J and K data combined with 3.6 micro-m and 4.5 micro-m Spitzer bands and optical BVRi'z' imaging from the Subaru Telescope over 0.5 square degrees in the Subaru XMM-Newton Deep Field. We have created a new cluster-detection algorithm, based on the Friends-Of-Friends and Voronoi Tessellation methods, which utilises probability distribution functions derived from a photometric redshift analysis. We employ mock catalogues to understand the selection effects and contamination associated with the algorithm. The cluster catalogue contains 13 clusters at redshifts 0.61 <= z <= 1.39 with luminosities 10 L* <~ L_tot <~ 50 L*, corresponding to masses 5 x 10^13 M_sun <~ M_cluster <~ 3 x 10^14 M_sun for (M/M_sun) / (L/L_sun) = 75h. The measured sky surface density of ~ 10 deg^-2 for high-redshift (z=0.5-1.5), massive (>10^14 M_sun) clusters is precisely in line with theoretical predictions presented by Kneissl et al. (2001).

We construct a sample of extremely red objects (EROs) within the UKIDSS Ultra Deep Survey by combining the Early Data Release with optical data from the Subaru/XMM-Newton Deep Field. We find a total of 3715 objects over 2013 sq. arcmin with R-K>5.3 and K<=20.3, which is a higher surface density than found by previous studies. This is partly due to our ability to use a small aperture in which to measure colours, but is also the result of a genuine overdensity of objects compared to other fields. We separate our sample into passively-evolving and dusty star-forming galaxies using their RJK colours and investigate their radio properties using a deep radio map. The dusty population has a higher fraction of individually-detected radio sources and a higher mean radio flux density among the undetected objects, but the passive population has a higher fraction of bright radio sources, suggesting that AGNs are more prevalent among the passive ERO population.

We discuss the relative merits of mid-infrared and X-ray selection of type-2 quasars. We describe the mid-infrared, near-infrared and radio selection criteria used to find a population of redshift z~2 type-2 quasars which we previously argued suggests that most supermassive black hole growth in the Universe is obscured (Martinez-Sansigre et al., 2005). We present the optical spectra obtained from the William Herschel Telescope, and we compare the narrow emission line luminosity, radio luminosity and maximum size of jets to those of objects from radio-selected samples. This analysis suggests that these are genuine radio-quiet type-2 quasars, albeit the radio-bright end of this population. We also discuss the possibility of two different types of quasar obscuration, which could explain how the ~2-3:1 ratio of type-2 to type-1 quasars preferred by modelling our population can be reconciled with the ~1:1 ratio predicted by unified schemes.

We present an alternative method to X-ray surveys for hunting down the high-redshift type-2 quasar population, using Spitzer and VLA data on the Spitzer First Look Survey. By demanding objects to be bright at 24 microns but faint at 3.6 microns, and combining this with a radio criterion, we find 21 type-2 radio-quiet quasar candidates at the epoch at which the quasar activity peaked. Optical spectroscopy with the WHT confirmed 10 of these objects to be type-2s with 1.4 < z < 4.2 while the rest are blank. There is no evidence for contamination in our sample, and we postulate that our 11 blank-spectrum candidates are obscured by kpc-scale dust as opposed to dust from a torus around the accretion disk. By carefully modelling our selection criteria, we conclude that, at high redshift, 50-80 % of the supermassive black hole growth is obscured by dust.

We present first results from the study of the TOOT00 region consisting of 47 radio sources brighter than 100 mJy at 151 MHz. We have 81% spectroscopic redshift completeness. From the K-z diagram we deduce that the host galaxies are similar to ~3 L* passively evolved elliptical galaxies and thus estimate the redshifts of the 9 sources without a secure spectroscopic redshift yielding a median redshift of 1.287. Above the RLF break we have a quasar fraction f ~ 0.3 although the quasars appear reddened; below the RLF break f -> 0 if we exclude flat-spectrum radio sources. We present a histogram of the number of TOOT00 radio sources versus their redshift which looks broadly like the Willott et al. (2001) prediction for TOOT, although the observed ratio of high to low redshift objects is somewhat lower than the prediction.

Supermassive black holes underwent periods of exponential growth during which we seem them as quasars in the distant Universe. The summed emission from these quasars generates the cosmic X-ray background, the spectrum of which has been used to argue that most black-hole growth is obscured. There are clear examples of obscured black-hole growth in the form of `type-2' quasars, but their numbers are fewer than expected from modelling of the X-ray background. Here we report on the direct detection of a population of distant type-2 quasars which is at least comparable in size to the well-known unobscured type-1 population. We selected objects that have mid-infrared and radio emissions characteristic of quasars, but which are faint at near-infrared and optical wavelengths. This population is responsible for most of the black hole growth in the young Universe and, throughout cosmic history, black-hole growth occurs in the dusty, gas-rich centres of active galaxies.

We extend our previous analysis which used generalized luminosity functions (GLFs) to predict the number of quasars and galaxies in low-radio-frequency-selected samples as a function of redshift, radio luminosity, narrow-emission-line luminosity and type of unified scheme. Our extended analysis incorporates the observed submillimetre (850 micron) flux densities of radio sources, employs a new method which allows us to deal with non detections, and focuses on the high-luminosity population. First, we conclude that the submillimetre luminosity L_{850} of low-frequency-selected radio sources is correlated with the bolometric luminosity L_{Bol} of their quasar nuclei via an approximate scaling relation L_{850} \propto L_{Bol}^{0.7 \pm 0.2}. Second, we conclude that there is quantitative evidence for a receding-torus-like physical process for the high-luminosity population within a two-population unified scheme for radio sources; this evidence comes from the fact that radio quasars are brighter in both narrow emission lines and submillimetre luminosity than radio galaxies matched in radio luminosity and redshift. Third, we note that the combination of a receding-torus-like scheme and the assumption that the observed submillimetre emission is dominated by quasar-heated dust yields a scaling relation L_{850} \propto L_{Bol}^{0.5} which is within the errors of that determined here for radio-selected quasars, and consistent with that inferred for radio-quiet quasars by Willott, Rawlings & Grimes (2003).

We present a clustering analysis of the Texas-Oxford NVSS Structure (TONS) radio galaxy redshift survey. This complete flux-limited survey consists of 268 radio galaxies with spectroscopic redshifts in three separate regions of the sky covering a total of 165 deg^2. By going to faint radio flux densities (s_1.4>3 mJy) but imposing relatively bright optical limits (E R 19.5), the TONS sample is optimised for looking at the clustering properties of low luminosity radio galaxies in a region of moderate (0 < z < 0.5) redshifts. We use the two point correlation function to determine the clustering strength of the combined TONS08 and TONS12 sub-samples and find a clustering strength of r_0(z)=8.7+/-1.6 Mpc (h=0.7). If we assume growth of structure by linear theory and that the median redshift is 0.3, this corresponds to r_0(0)=11.0+/-2.0 Mpc which is consistent with the clustering strength of the underlying host galaxies (~ 2.5 Lstar ellipticals) of the TONS radio galaxy population.

We present early results from the spectroscopic follow-up of a sample of candidate obscured AGN selected in the mid-infrared from the Spitzer First Look Survey. Our selection allows a direct comparison of the numbers of obscured and unobscured AGN at a given luminosity for the first time, and shows that the ratio of obscured to unobscured AGN at infrared luminosities corresponding to low luminosity quasars is ~1:1 at z~0.5. Most of our optically-faint candidate obscured AGN have the high-ionization, narrow-line spectra expected from type-2 AGN. A composite spectrum shows evidence for Balmer absorption lines, indicating recent star-formation activity in the host galaxies. There is tentative evidence for a decrease in the obscured AGN fraction with increasing AGN luminosity.

Current surveys may be on the verge of measuring the baryonic oscillations in the galaxy power spectrum which are clearly seen imprinted on the Cosmic Microwave Background. It has recently been proposed that these oscillations allow a `standard ruler' method of probing the equation of state of dark energy. In this paper we present a new calculation of the number of galaxies future adio telescopes will detect in surveys of the sky in neutral Hydrogen (HI). We estimate the likely statistical errors if the standard ruler method were to be applied to such surveys. We emphasise uncertainties in our calculations, and pinpoint the most important features of future HI surveys if they are to provide new constraints on dark energy via baryonic oscillations. Designs of future radio telescopes are required to have a large bandwidth (characterised by \beta$, the ratio of the instantaneous bandwidth to the bandwidth required by survey) and to have the widest instantaneous (1.4 GHz) field of view ($FOV$) possible. Given the expected sensitivity of a future Square Kilometre Array (SKA), given that half of its collecting area will be concentrated in a core of diameter $\sim 5 ~ \rm km$, and given a reasonable survey duration ($T_0$ $\sim$ 1 yr), we show that there will be negligible shot noise on a power spectrum derived from HI galaxies out to redshift $z \simeq 1.5$. To access the largest cosmic volume possible by surveying all the sky available, we argue that $\beta$, $T_0$ and $FOV$ must obey the relation $\beta FOV T_0 \gtsimeq 10 \rm deg^2$ yr. An $\sim$1-yr SKA survey would then contain $\gtsimeq 10^{9} (f_{\rm sky}/0.5)$ HI galaxies and provide constraints on the dark-energy parameter $w$ of order $\Delta w \simeq 0.01 (f_{\rm sky}/0.5)^{-0.5}$.

The relationships between supermassive black holes and the properties of their associated dark-matter halos imply that outflows from accreting black holes provide a feedback mechanism regulating galaxy formation. Accreting black holes with weak or undetectable radio jets (radio-quiet quasars) outnumber those with powerful jets (radio-loud quasars) by a factor ~10-100, so powerful-jet outflows are often neglected. However, whenever powerful jets are triggered, there is a dramatic (factor >100) step-function increase in the efficiency of feedback. We use a feedback model, together with the measured space density of flat-spectrum radio-loud quasars, to show that a powerful-jet episode probably occurred in every protocluster in the Universe. Before jet triggering, there was time for gravitational collapse to create many (~10-100) surrounding protogalaxies massive enough to host radio-quiet quasars. After triggering, the powerful jet pushes back and heats ionized gas so that it cannot fall onto these protogalaxies and cool. Once neutral/molecular gas reservoirs become exhausted, there is a synchronized shut down in both star-formation and black-hole activity throughout the protocluster. These considerations imply that radio-loud quasars have a profound influence on the evolution of all the galaxies seen in clusters today.

The present-day Universe is seemingly dominated by dark energy and dark matter, but mapping the normal (baryonic) content remains vital for both astrophysics - understanding how galaxies form - and astro-particle physics - inferring properties of the dark components. The Square Kilometre Array (SKA) will provide the only means of studying the cosmic evolution of neutral Hydrogen (HI) which, alongside information on star formation from the radio continuum, is needed to understand how stars formed from gas within dark-matter over-densities and the roles of gas accretion and galaxy merging. `All hemisphere' HI redshift surveys to redshift 1.5 are feasible with wide-field-of-view realizations of the SKA and, by measuring the galaxy power spectrum in exquisite detail, will allow the first precise studies of the equation-of-state of dark energy. The SKA will be capable of other uniquely powerful cosmological studies including the measurement of the dark-matter power spectrum using weak gravitational lensing, and the precise measurement of H0 using extragalactic water masers. The SKA is likely to become the premier dark-energy-measuring machine, bringing breakthroughs in cosmology beyond those likely to be made possible by combining CMB (e.g. Planck), optical (e.g. LSST, SNAP) and other early-21st-century datasets.

The epoch of reionization (EoR) sets a fundamental benchmark in cosmic structure formation, corresponding to the formation of the first luminous objects that act to ionize the neutral intergalactic medium (IGM). Recent observations at near-IR and radio wavelengths imply that we are finally probing into this key epoch of galaxy formation at $z \ge 6$. The Square Kilometer Array (SKA) will provide critical insight into the EoR in a number of ways. First, the ability of the SKA to image the neutral IGM in 21cm emission is a truly unique probe of the process of reionization, and is recognized as the next necessary and fundamental step in our study of the evolution of large scale structure and cosmic reionization. Second, study of HI 21cm absorption toward the first radio loud objects probes small to intermediate scale structure in the neutral 'cosmic web', as well as HI in the first collapsed structures (proto-disks and mini-halos). And third, the incomparable sensitivity of the SKA allows for the study of the molecular gas, dust, and star formation activity in the first galaxies, as well as the radio continuum emission from the first accreting massive black holes. Such objects will be obscured at optical wavelengths due to absorption by the neutral IGM.

We argue that the Square Kilometre Array has the potential to make both redshift (HI) surveys and radio continuum surveys that will revolutionize cosmological studies, provided that it has sufficient instantaneous field-of-view that these surveys can cover a hemisphere in a timescale ~1 yr. Adopting this assumption, we focus on two key experiments which will yield fundamental new measurements in cosmology, characterizing the properties of the mysterious dark energy which dominates the dynamics of today's Universe. Experiment I will map out ~10^9 HI galaxies to redshift z~1.5, providing the premier measurement of the clustering power spectrum of galaxies: accurately delineating the acoustic oscillations and the `turnover'. Experiment II will quantify the cosmic shear distortion of ~10^10 radio continuum sources, determining a precise power spectrum of the dark matter, and its growth as a function of cosmic epoch. We contrast the performance of the SKA in precision cosmology with that of other facilities which will, probably or possibly, be available on a similar timescale. We conclude that data from the SKA will yield transformational science as the direct result of four key features: (i) the immense cosmic volumes probed, exceeding future optical redshift surveys by more than an order of magnitude; (ii) well-controlled systematic effects such as the narrow `k-space window function' for Experiment I and the accurately-known `point-spread function' for Experiment II; (iii) the ability to measure with high precision large-scale modes in the clustering power spectra, for which nuisance effects such as non-linear structure growth, peculiar velocities and `galaxy bias' are minimised; and (iv) different degeneracies between key parameters to those which are inherent in the CMB.

The Square Kilometer Array (SKA) represents the next major, and natural, step in radio astronomical facilities, providing two orders of magnitude increase in collecting area over existing telescopes. In a series of meetings, starting in Groningen, the Netherlands (August 2002) and culminating in a `science retreat' in Leiden (November 2003), the SKA International Science Advisory Committee (ISAC), conceived of, and carried-out, a complete revision of the SKA science case (to appear in New Astronomy Reviews). This preface includes: (i) general introductory material, (ii) summaries of the key science programs, and (iii) a detailed listing of standards and assumptions used in the revised science case.

In this paper we investigate how the Square Kilometre Array (SKA) can aid in determining the evolutionary history of active galactic nuclei (AGN) from redshifts z= 0 ->6. Given the vast collecting area of the SKA, it will be sensitive to both `radio-loud' AGN and the much more abundant `radio-quiet' AGN, namely the radio-quiet quasars and their `Type-II' counterparts, out to the highest redshifts. Not only will the SKA detect these sources but it will also often be able to measure their redshifts via the Hydrogen 21 cm line in emission and/or absorption. We construct a complete radio luminosity function (RLF) for AGN, combining the most recent determinations for powerful radio sources with an estimate of the RLF for radio-quiet objects using the hard X-ray luminosity function of Ueda et al., including both Type-I and Type-II AGN. We use this complete RLF to determine the optimal design of the SKA for investigating the accretion history of the Universe for which it is likely to be a uniquely powerful instrument.

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.

We have analyzed the optical and radio properties of a sample of 3CR FR I radio galaxies which have Hubble Space Telescope imaging capable of detecting optical cores. The jet powers of the FR I radio galaxies are estimated from their low-frequency radio luminosities, and the optical core luminosity is taken as an upper limit on the emission from any unobscured accretion disc. We argue that if the accretion discs in these sources are assumed to be advection dominated accretion flows (ADAFs), or adiabatic inflow-outflow solution (ADIOS) flows, then the Blandford-Znajek mechanism (Blandford & Znajek 1977) provides insufficient power to explain the high radio luminosities of at least a third, and perhaps all, of the sample. We suggest instead that a significant fraction (the 'high-jet-power' third), and perhaps most, of the 3CR FR I radio galaxies have normal accretion discs, but that their optical cores can be hidden, with any HST-detected optical synchrotron emission coming from jets on scales larger than the obscuring material. A normal accretion disc hypothesis, at least for the high-jet-power third of the 3CR FR Is, explains why narrow-line luminosity correlates with radio luminosity. It also explains why one object in the sample (3C 386) has an observed broad-line nucleus. We conclude that there is no evidence to suggest that there is a difference in accretion mode between FR I and FR II radio galaxies.

The Universe became fully reionized, and observable optically, at a time corresponding to redshift z ~ 6.5, so it is only by studying the HI and molecular absorption lines against higher-redshift, radio-loud sources that one can hope to make detailed studies of the earliest stages of galaxy formation. At present no targets for such studies are known. In these proceedings we describe a survey which is underway to find radio-loud quasars at z > 6.5, and present broad-band SEDs of our most promising candidates.

We present submillimetre (SCUBA) photometry of a sample of eight high-redshift (2.5 < z < 3.5) radio quasars from two redshift surveys: the TEXOX-1000, or TOOT, survey; and the 7C Redshift Survey (7CRS). Unlike the powerful high-redshift radio sources observed previously in the submillimetre, these radio sources are typical of those dominating the radio luminosity density of the population. We detect just two of the TOOT/7CRS targets at 850 microns, and one of these detections is probably due to synchrotron emission rather than dust. The population represented by the other six objects is detected in a statistical sense with their average 850-micron flux density implying that they are similar to low-redshift, far-infrared-luminous quasars undergoing at most moderate (< 200 solar masses per year) starbursts. By considering all the SCUBA data available for radio sources, we conclude that positive correlations between rest-frame far-infrared luminosity, 151-MHz luminosity and redshift, although likely to be present, are hard to interpret because of subtle selection and classification biases, small number statistics and uncertainties concerning synchrotron contamination and K-correction. We argue that there is not yet any compelling evidence for significant differences in the submillimetre properties of radio-loud and radio-quiet quasars at high redshift.

We present a new radio sample, 6C** designed to find radio galaxies at z > 4 and discuss some of its near-infrared imaging follow-up results.

We present new results on the clustering and three-dimensional distribution of radio galaxies from the Texas-Oxford NVSS Structure (TONS) survey. The TONS survey was constructed to look at the distribution of radio galaxies in a region of moderate (0 < z < 0.5) redshifts by matching NVSS sources with objects in APM catalogues to obtain a sample of optically bright (R < 19.5), radio faint (1.4-GHz flux density S_1.4 > 3 mJy) radio galaxies over large areas on the sky. We find that redshift spikes, which represent large concentrations of radio galaxies which trace (~ 100 Mpc^3) super-structures are a common phenomena in these surveys. Under the assumption of quasi-linear structure formation theory and a canonical radio galaxy bias, the structures represent ~ 4-5 sigma peaks in the primordial density field and their expected number is low. The most plausible explanation for these low probabilities is an increase in the radio galaxy bias with redshift. To investigate potential mechanisms which have triggered the radio activity in these galaxies - and hence may account for an increase in the bias of this population, we performed imaging studies of the cluster environment of the radio galaxies in super-structure regions. Preliminary results show that these radio galaxies may reside preferentially at the edges of rich clusters. If radio galaxies are preferentially triggered as they fall towards rich clusters then they would effectively adopt the cluster bias.

The results of analysis of HST I-band imaging of a sample of 41 radio galaxies spanning three orders of magnitude in radio luminosity at redshift z~0.5 are presented. The full sample is drawn from four complete, low-frequency selected radio samples with progressively fainter flux-density limits (3CRR, 6CE, 7CRS and the new TexOx-1000 sample). Modelling of the HST imaging data shows that the host galaxies have surface-brightness distributions consistent with those expected for classic ellipticals (Sersic parameter, beta~0.25), with beta in the range 0.17<beta<0.30, and a mean of <beta>=0.23\pm0.01. The host-galaxy luminosities are comparable with those of galaxies drawn from the bright end of the local cluster galaxy luminosity function, spanning the range 0.7L^{\star}<L<10L^{\star}, with a mean of 3.2\pm0.3 L^{\star}. In addition, the radio galaxies are shown to follow a Kormendy relation indistinguishable from that of powerful low-redshift radio galaxies. Combining our new results with those in the literature it is found that the scalelengths and Kormendy relations of 3C-class radio galaxies do not vary significantly over the redshift range 0.0<z<0.8, providing no evidence for dynamical evolution of this class of host galaxy within this redshift interval. Converting the host-galaxy luminosities into black-hole mass estimates predicts that the radio galaxies harbour central black holes with masses in the range 10^{8.1}\Msun <M_{bh}< 10^{9.5}\Msun, with a geometric mean of <M_{bh}>=10^{8.87\pm 0.04}\Msun. Finally, a significant (~3 sigma) correlation is found between black-hole mass and 151-MHz radio luminosity for those objects in the sample with either high-excitation nuclear spectra or classical double radio structures. (abridged)

We use a principal components analysis of radio-selected (3CRR, 6CE and 7CRS) AGN datasets to define two parameters related to low-frequency (151 MHz) radio luminosity L_151 and [OIII] luminosity L_[OIII]: a parameter alpha encoding the L_151 - L_[OIII] correlation and a parameter beta encoding scatter about this correlation. We describe methods for constructing generalized luminosity functions (GLFs) based on alpha, beta, redshift and schemes for unifying quasars and radio galaxies. Luminosity-dependent unified schemes (e.g. a receding-torus scheme) have been invoked to explain the low quasar to radio galaxy fraction at low alpha and the differences in emission-line luminosities of radio quasars and radio galaxies. With the constraints of the 3CRR, 6CE and 7CRS datasets and radio source counts, our GLF approach was used to determine whether a receding-torus-like scheme is required if there are two populations of radio sources: one at low alpha, consisting of 'starved AGN'; the other at high alpha consisting of 'Eddington-tuned AGN'. Because of the overlap between these two populations and the effects of the beta parameter, schemes with or without a receding torus can produce a low quasar fraction at low alpha and differences in [OIII] luminosity between radio galaxies and quasars. The receding torus may be a physical process important in one or more populations of radio sources, but this is not yet proved either by the quasar fraction or the emission-line properties of radio-selected samples.

The results of a sub-millimeter survey of SDSS broad CIV absorption line quasars is discussed. It is found that the sub-millimeter flux distribution of BAL quasars is similar to that of non-BAL quasars. This is consistent with the idea that all quasars contain broad absorption line regions, but only a fraction of them are visible along our line-of-sight. The observations are inconsistent with BAL quasars being observed at a special evolutionary epoch co-inciding with a high star-formation rate and dust mass.

The Universe became fully reionized, and observable optically, at a time corresponding to redshift z~6.5, so it is only by studying the HI and molecular absorption lines against higher-redshift, radio-loud sources that one can hope to make detailed studies of the earliest stages of galaxy formation. At present no targets for such studies are known. In these proceedings we describe a survey which is underway to find radio-loud quasars at z > 6.5.

We have obtained near-infrared spectroscopy of a small sample of powerful Compact Steep-Spectrum (CSS) radio sources mainly, but not exclusively, from the 3CR sample. We find no differences between the distributions in the equivalent width and luminosity of the [OIII] 5007A line for our sample and other larger, presumably older, high-redshift 3C objects, suggesting that the underlying quasar luminosity remains roughly constant as quasars age. We also find a possible broad line in 3C241, adding to recent evidence for broad lines in some radio galaxies.

We have carried out the first systematic survey of the sub-millimeter properties of broad absorption line (BAL) quasars. 30 BAL quasars drawn from a homogeneously selected sample from the Sloan Digital Sky Survey at redshifts 2<z<2.6 were observed with the SCUBA array at the JCMT to a typical rms sensitivity of 2.5 mJy. Eight quasars were detected at > 2 sigma significance, four of which are at > 3 sigma significance. The far-infrared luminosities of these quasars are > 10^{13} L_solar. There is no correlation of sub-millimeter flux with either the strength of the broad absorption feature or with absolute magnitude in our sample. We compare the sub-millimeter flux distribution of the BAL quasar sample with that of a sample of quasars which do not show BAL features in their optical spectra and find that the two are indistinguishable. BAL quasars do not have higher sub-millimeter luminosities than non-BAL quasars. These findings are consistent with the hypothesis that all quasars would contain a BAL if viewed along a certain line-of-sight. The data are inconsistent with a model in which the BAL phenomenon indicates a special evolutionary stage which co-incides with a large dust mass in the host galaxy and a high sub-millimeter luminosity. Our work provides constraints on alternative evolutionary explanations of BAL quasars.

The results of Hubble Space Telescope and UKIRT imaging observations are presented for a sample of 11 6C radio galaxies with redshifts 0.85 < z < 1.5. The observations of the 6C sources reveal a variety of different features, similar to those observed around the higher luminosity of the aligned emission appears less extreme in the case of the 6C radio galaxies. For both samples, the aligned emission clearly cannot be explained by a single emission mechanism; line emission and related nebular continuum emission, however, often provide a significant contribution to the aligned emission.

We present unequivocal evidence for a huge (~80 x 100 x 100 Mpc^3) super-structure at redshift z=0.27 in the 3-D distribution of radio galaxies from the TONS08 sample, confirming tentative evidence for such a structure from the 7C redshift survey (7CRS). A second, newly discovered super-structure is also less securely found at redshift 0.35 (of dimensions ~100 x 100 x 100 Mpc^3). We present full observational details on the TONS08 sample which was constructed to probe structures in the redshift range 0 < z < 0.5 by matching NVSS sources with objects in APM catalogues to obtain a sample of optically bright (E ~ R < 19.83), radio faint (1.4-GHz flux density S_1.4 > 3 mJy) radio galaxies in the same 25 deg^2 area as part-II of the 7CRS. Out of the total sample size of 84 radio galaxies, at least 25 are associated with the two ~100 Mpc-scale super-structures. We use quasi-linear structure formation theory to estimate the number of such structures expected in the TONS08 volume if the canonical value for radio galaxy bias is assumed. Under this assumption, the structures represent ~ 4-5 sigma peaks in the primordial density field and their expected number is low (~10^-2 - 10^-4). Fortunately, there are several plausible explanations (many of which are testable) for these low probabilities in the form of potential mechanisms for boosting the bias on large scales. These include: the association of radio galaxies with highly biased rich clusters in super-structures, enhanced triggering by group/group mergers, and enhanced triggering and/or redshift space distortion in collapsing systems as the growth of super-structures moves into the non-linear regime (abridged).

The next generation of surveys of extragalactic radio sources will be dominated by different types of objects than the jetted-AGN that dominate surveys like 3C, 6C and 7C. Before radio astronomy becomes concentrated on the new types of object, it is vital that we understand the cosmological importance of the jetted-AGN that have been studied for many years. I argue that, as observational manifestations of Eddington-tuned outflow events, these objects may have more significance for galaxy/cluster formation and evolution than is typically appreciated. Outstanding problems in galaxy formation may be solved by cosmological simulations in which this type of outflow, as well as other types associated with weaker jets, are properly taken into account. I will highlight areas of ignorance which are currently hindering attempts to do this.

Radio galaxies are excellent at tracing large-scale structure due to their high bias. We present new results from the TONS08 radio galaxy redshift survey. We find unequivocal evidence for a huge (at least 80 x 80 x 100 Mpc^3) super-structure at redshift z=0.27, confirming tentative evidence for such a structure from the 7C redshift survey (7CRS). A second, newly discovered super-structure is also tentatively found at redshift 0.35 (of dimensions at least 100 x 100 x 100 Mpc^3). Out of the total sample size of 84 radio galaxies, at least 25 are associated with the two super-structures. We use quasi-linear structure formation theory to estimate the number of such structures expected in the TONS08 volume if the canonical value for radio galaxy bias is assumed. Under this assumption, the structures represent ~ 4-5 sigma peaks in the primordial density field and their expected number is low (10^{-2}-10^{-4}). Fortunately, there are several plausible explanations (many of which are testable) for these low probabilities in the form of potential mechanisms for boosting the bias on large scales. These include: the association of radio galaxies with highly biased rich clusters in super-structures, enhanced triggering by group/group mergers, and enhanced triggering and/or redshift space distortion in collapsing systems as the growth of super-structures moves into the non-linear regime. Similar structures could have been missed in previous surveys because of the effects of Poisson-sampling fluctuations.

The TexOx-1000 (TOOT) radio source redshift survey is designed to find and study typical radio-loud active galaxies to high redshift. They are typical in the same sense that L* galaxies are typical of galaxies in the optical. Previous surveys have only included the most luminous, rare objects at and beyond the peak of activity at z~2, but in going a factor of 100 fainter than the 3C survey, and in assembling a large sample, TOOT probes for the first time the objects that dominate the radio luminosity density of the universe at high redshift. Here we describe the current status of the TOOT survey and draw preliminary conclusions about the redshift distribution of the radio sources. So far, ~520 of the 1000 radio sources have redshifts, with ~440 of those in well-defined, complete, sub-regions of the survey. For these we find a median redshift of z=1, but the measured redshift distribution has a deficit of objects with z~2, when compared to predictions based on extrapolating luminosity functions constrained by higher-flux-density samples. These are the more luminous objects that usually show emission lines, and which should not be missed in the survey unless they are heavily reddened. The deficit may be real, but it would not be too surprising to find a population of faint, reddened radio sources at z~2-3 among the TOOT sources yet to have accurate redshifts.

The TexOx Cluster (TOC) Survey uses overdensities of radiosources in the NVSS to trace clusters of galaxies. The links between radiosources and rich environments make this a powerful way to find clusters which may potentially be overlooked by other selection techniques. By including constraints from optical surveys, TOC is an extremely efficient way to find clusters at high redshift. One such field, TOC J0233.3+3021, contains at least one galaxy cluster (at z ~ 1.4) and has been detected using the Sunyaev-Zel'dovich (SZ) effect. Even in targetted deep optical observations, however, distinguishing the cluster galaxies from the background is difficult, especially given the tendency of TOC to select fields containing multiple structures at different redshifts.

The sources discovered in deep hard X-ray surveys with 2-8 keV fluxes of 10^-14 erg cm^-2 s^-1 make up the bulk of the X-ray background at these energies. We present here detailed multi-wavelength observations of three such sources from the ELAIS Deep X-ray Survey. The observations include sensitive near-infrared spectroscopy with the Subaru Telescope and X-ray spectral information from the Chandra X-ray Observatory. The sources observed all have optical-to-near-IR colours redder than an unobscured quasar and comprise a reddened quasar, a radio galaxy and an optically-obscured AGN. The reddened quasar is at a redshift z=2.61 and shows a very large X-ray absorbing column of N_H approx 3.10^23 cm^-2. This contrasts with the relatively small amount of dust reddening, implying a gas-to-dust ratio along the line-of-sight a hundred times greater than that of the Milky Way. The radio galaxy at z=1.57 shows only narrow emission lines, but has a surprisingly soft X-ray spectrum. The softness of this spectrum either indicates an unusually low gas-to-dust ratio for the absorbing medium or X-ray emission related to the young radio source. The host galaxy is extremely red (R-K=6.4) and its optical/near-IR spectrum is best fit by a strongly reddened (A_V~2) starburst. The third X-ray source discussed is also extremely red (R-K=6.1) and lies in a close grouping of three other R-K>6 galaxies. No emission or absorption lines were detected from this object, but its redshift (and that of one of the nearby galaxies) are constrained by SED-fitting to be just greater than z=1. The extremely red colours of these two galaxies can be accounted for by old stellar populations. These observations illustrate the diverse properties of hard X-ray selected AGN.

We present K-band imaging of all 49 radio galaxies in the 7C-I and 7C-II regions of the 7C Redshift Survey (7CRS). The low-frequency (151 MHz) selected 7CRS sample contains all sources with flux-densities S_151 > 0.5 Jy in three regions of the sky. We combine the K-band magnitudes of the 7CRS radio galaxies with those from the 3CRR, 6CE and 6C* samples to investigate the nature of the relationship between K-magnitude and redshift and whether there is any dependence upon radio luminosity. We find that radio galaxies appear to belong to a homogeneous population which formed the bulk of their stars at high redshifts (z_f>5) and evolved passively from then until they reach a mean present-day luminosity of 3L*. We find a significant difference between the K-magnitudes of the 7CRS and 3CRR radio galaxies with the 7CRS galaxies being \~0.55 mag fainter at all redshifts. The cause of this weak correlation between stellar and radio luminosities probably lies in mutual correlations of these properties with the central black hole mass. We compare the evolution-corrected host luminosities at a constant radio luminosity and find that the typical host luminosity (mass) increases by approximately 1L* from z~2 to z~0.5 which, although a much smaller factor than predicted by semi-analytic models of galaxy formation, is in line with results on optically-selected quasars. Our study has therefore revealed that the small dispersion in stellar luminosity of radio galaxies around 3L* includes subtle but significant differences between the host galaxies of extreme- and moderate-power radio sources at fixed redshift, and between those of high- and low-redshift radio sources at fixed radio luminosity.

We present an up-to-date critique of the physical basis for the spectral ageing method. We find that the number of cases where this method may be meaningfully applied to deduce the ages of classical double radio sources is small indeed. This critique is much more than merely a re-expression of anxieties about the calibration of spectral ageing (which have been articulated by others in the past).

(abridged) The results of deep long-slit optical spectroscopy for a sample of 8 6C radio galaxies at z~1 are presented. Emission line ratios are derived for emission lines with rest-frame wavelengths of 1500 - 4500Ang, and the kinematic propertiess of the emission line gas are derived from an analysis of the 2-d structure of thee [OII]3727 Ang emission line at ~5Ang spectral resolution. In general, the 6C spectra display many characteristics similar to those of more powerful 3CR sources at the same redshifts. The emission line region gas kinematics are more extreme for the smaller radio sources in the sample, which often display distorted velocity profiles. The ionization state of the emission line region also varies with radio size: the spectra of large radio sources (>120kpc) are consistent with photoionization by an obscured AGN, whilst smaller (<120kpc) sources typically exist in a lower ionization state and have spectra which are better explained by additional ionization due to shocks associated with the expanding radio source. The kinematic and ionization properties of the 6C radio galaxies are clearly linked. As for the 3CR sources, smaller radio sources also typically possess more extensive emission line regions, with enhanced emission line luminosities. It is clear that the best interpretation of the spectra of radio sources requires a combination of ionization mechanisms.

The results of a spectroscopic analysis of 3CR and 6C radio galaxies at redshift z~1 are contrasted with the properties of lower redshift radio galaxies, chosen to be matched in radio luminosity to the 6C sources studied at z~1, thus enabling the P-z degeneracy to be broken. Partial rank correlations and principal component analysis have been used to determine which of z and P are the critical parameters underlying the observed variation of the ionization state andd kinematics of the emission line gas. [OII]/H-beta is shown to be a useful ionization mechanism diagnostic. Statistical analysis of the data shows that the ionization state of the emission line gas is strongly correlated with radio power, once the effects of other parameters are removed. No dependence of ionization state on z is observed, implying that the ionization state of the emission line gas is solely a function of the AGN properties rather than the hostt galaxy and/or environment. Statistical analysis of the kinematic properties of the emission line gas shows that these are strongly correlated independently withh both P and z. The correlation with redshift is the stronger of the two, suggesting that host galaxy composition or environment may play a role in producing the less extreme gas kinematics observed in the emission line regions of low redshift galaxies. For both the ionization and kinematic properties of thee galaxies, the independent correlations observed with radio size are strongest. Radio source age is a determining factor for the extended emission line regions.

We present optical spectroscopy of all 49 radio galaxies in the 7C-I and 7C-II regions of the 7C Redshift Survey (7CRS). The low-frequency (151 MHz) selected 7CRS sample contains all sources with flux-densities S_151 > 0.5 Jy in three regions of the sky; 7C-I and 7C-II were chosen to overlap with the 5C6 and 5C7 surveys respectively, and cover a total sky area of 0.013 sr. The sample has been completely identified and spectroscopy of the quasars and broad-lined radio galaxies has been presented in Willott et al. (1998). Only seven of the radio galaxies do not have redshift determinations from the spectroscopy, giving a redshift completeness for the sample of >90%. The median redshift of the 7CRS is 1.1. We present a composite 0.2<z<0.8 7CRS radio galaxy spectrum and investigate the strengths of the 4000 Angstrom breaks in these radio galaxies. We find an anti-correlation between the 4000 Angstrom break strength and emission line luminosity, indicating that departures from old elliptical galaxy continuum shapes are most likely due to non-stellar emission associated with the active nucleus.

We present deep near-infrared images, taken with Subaru Telescope, of the region around the z=1.08 radio source 3C 356 which show it to be associated with a poor cluster of galaxies. We discuss evidence that this cluster comprises two subclusters traced by the two galaxies previously proposed as identifications for 3C 356, which both seem to harbour AGN, and which have the disturbed morphologies expected if they underwent an interpenetrating collision at the time the radio jets were triggered. We explain the high luminosity and temperature of the diffuse X-ray emission from this system as the result of shock-heating of intracluster gas by the merger of two galaxy groups. Taken together with the results on other well-studied powerful radio sources, we suggest that the key ingredient for triggering a powerful radio source, at least at epochs corresponding to z~1, is a galaxy--galaxy interaction which can be orchestrated by the merger of their parent subclusters. This provides an explanation for the rapid decline in the number density of powerful radio sources since z~1. We argue that attempts to use distant radio-selected clusters to trace the formation and evolution of the general cluster population must address ways in which X-ray properties can be influenced by the radio source, both directly, by mechanisms such as inverse-Compton scattering, and indirectly, by the fact that the radio source may be preferentially triggered at a specific time during the formation of the cluster.

We present submillimetre photometry of 11 3CR and 6CE radio quasars of similar 151-MHz radio luminosity and redshifts to the radio galaxies studied at 1.3<z<2 by Archibald et al. (2001). We detect all bar one quasar at 850 microns, and five quasars are confirmed as dusty hyperluminous objects by detections at shorter wavelengths. Our observations reveal a clear difference between the submillimetre luminosity distributions of the radio quasars and a matched sample of radio galaxies: the quasars are ~5-times brighter than the radio galaxies and >~2-times brighter accounting for possible synchrotron contamination. This difference is in quantitative agreement with a receding-torus unified scheme for radio sources in which the torus opening angle depends on quasar optical luminosity, provided that there is a close relationship between the optical and submillimetre luminosities. The implication is that quasar-heated dust dominates the submillimetre emission for powerful radio quasars at z~1.5. We find a significant anti-correlation between submillimetre/far-infrared luminosity L_FIR and radio source age in the sense that hyperluminous quasars tend to be associated with young (<10^7 yr) sources. This suggests that the processes controlling L_FIR are synchronized with the jet-triggering event and that L_FIR declines on a ~10^7 yr timescale. There is evidence for a small amount of obscuration in the hyperluminous quasars from reddening and associated or broad absorption lines. We speculate that shortly after jet triggering, dust is more widely distributed than at later times and is heated by the quasar nucleus and/or a short-lived synchronized starburst.

Using a new radio sample, 6C* designed to find radio galaxies at z > 4 along with the complete 3CRR and 6CE sample we extend the radio galaxy K-z relation to z~4.5. The 6C* K-z data significantly improve delineation of the K-z relation for radio galaxies at high redshift (z >2). In a spatially flat universe with a cosmological constant (Omega_M=0.3 and Omega_Lambda=0.7), the most luminous radio sources appear to be associated with galaxies with a luminosity distribution with a high mean (~5L*), and a low dispersion (~0.5 mag) which formed their stars at epochs corresponding to z >~ 2.5.

We present a deep image of the optically-powerful quasar E1821+643 at 18cm made with the Very Large Array (VLA). This image reveals radio emission, over 280 kpc in extent, elongated way beyond the quasar's host galaxy. Its radio structure has decreasing surface brightness with increasing distance from the bright core, characteristic of FRI sources (Fanaroff & Riley 1974). Its radio luminosity at 5GHz falls in the classification for `radio-quiet' quasars (it is only 10^23.9 W/Hz/sr; see e.g. Kellermann et al 1994). Its radio luminosity at 151MHz (which is 10^25.3 W/Hz/sr) is at the transition luminosity observed to separate FRIs and FRIIs. Hitherto, no optically-powerful quasar had been found to have a conventional FRI radio structure. For searches at low-frequency this is unsurprising given current sensitivity and plausible radio spectral indices for radio-quiet quasars. We demonstrate the inevitability of the extent of any FRqI radio structures being seriously under-estimated by existing targetted follow-up observations of other optically-selected quasars, which are typically short exposures of z > 0.3 objects, and discuss the implications for the purported radio bimodality in quasars. The nature of the inner arcsec-scale jet in E1821+643, together with its large-scale radio structure, suggest that the jet-axis in this quasar is precessing (cf. Galactic jet sources such as SS433). A possible explanation for this is that its central engine is a binary whose black holes have yet to coalesce. The ubiquity of precession in `radio-quiet' quasars, perhaps as a means of reducing the observable radio luminosity expected in highly-accreting systems, remains to be established.

We present some initial results from the TexOx (Texas-Oxford) Cluster (TOC) survey - a new method of selecting distant galaxy clusters. The cosmic evolution of the radio source population suggests that some massive clusters at high redshift will contain several radio-loud AGN. We searched for extreme over-densities at ~mJy levels in 7' x 7' boxes within the NVSS radio catalogue, covering a large (~1100 square degree) sky area. We have acquired optical images for ~130 cluster candidates, and followed up a subset of these with the VLA, and with Calar Alto near-IR imaging. Ryle Telescope observations have yielded at least one Sunyaev-Zel'dovich (SZ) detection of a massive z~1 system. Spectroscopic follow-up with 8-m class telescopes is in progress.

We present the results of radio, optical and near-infrared observations of the field of TOC J0233.3+3021, a cluster of milliJansky radiosources from the TexOx Cluster survey. In an observation of this field with the Ryle Telescope (RT) at 15 GHz, we measure a decrement in the cosmic microwave background (CMB) of $-675 \pm 95 \mu$Jy on the RT's $\approx$ 0.65 k$\lambda$ baseline. Using optical and infrared imaging with the McDonald 2.7-m Smith Reflector, Calar Alto 3.5-m telescope and UKIRT, we identify the host galaxies of five of the radiosources and measure magnitudes of $R \approx 24$, $J \approx 20$, $K \approx 18$. The CMB decrement is consistent with the Sunyaev-Zel'dovich (SZ) effect of a massive cluster of galaxies, which if modelled as a spherical King profile of core radius $\theta_C = 20^{\prime\prime}$ has a central temperature decrement of $900 \mu$K. The magnitudes and colours of the galaxies are consistent with those of old ellipticals at $z \sim 1$. We therefore conclude that TOC J0233.3+3021 is a massive, high redshift cluster. These observations add to the growing evidence for a significant population of massive clusters at high redshift, and demonstrate the effectiveness of combining searches for AGN `signposts' to clusters with the redshift-independence of the SZ effect.

We use three samples (3CRR, 6CE and 6C*) selected at low radio frequency to constrain the cosmic evolution in the radio luminosity function (RLF) for the `most luminous' steep-spectrum radio sources. Although intrinsically rare, such sources give the largest possible baseline in redshift for the complete flux-density-limited samples currently available. Using parametric models to describe the RLF which incorporate distributions in radio spectral shape and linear size as well as the usual luminosity and redshift, we find that the data are consistent with a constant comoving space density between z~2.5 and z~4.5. We find this model is favoured over a model with similar evolutionary behaviour to that of optically-selected quasars (i.e. a roughly Gaussian distribution in redshift) with a probability ratio of ~25:1 and ~100:1 for spatially-flat cosmologies with Omega_Lambda = 0 and Omega_Lambda = 0.7 respectively. Within the uncertainties, this evolutionary behaviour may be reconciled with the shallow decline preferred for the comoving space density of flat-spectrum sources by Dunlop & Peacock (1990) and Jarvis & Rawlings (2000), in line with the expectations of Unified Schemes.

(Abridged) This is the second in a series of three papers which present and interpret basic observational data on the 6C* 151-MHz radio sample: a low-frequency selected sample which exploits filtering criteria based on radio properties (steep spectral index and small angular size) to find radio sources at redshift z > 4 within a 0.133sr patch of sky. We present results of a programme of optical spectroscopy which has yielded redshifts in the range 0.5 < z < 4.4 for the 29 sources in the sample, all but six of which are secure. We find that the fil tering criteria used for 6C* are very effective in excluding the low-redshift, low-luminosity radio sources: the median redshift of 6C* is z~1.9 compared to z~1.1 for a complete sample matched in 151-MHz flux density. By combining the emission-line dataset for the 6C* radio sources with those for the 3CRR, 6CE and 7CRS samples we establish that z > 1.75 radio galaxies follow a rough proportionality between Lyalpha- and 151 MHz-luminosity which, like similar correlations seen in samples of lower-redshift radio sources, are indicative of a primary link between the power in the source of the photoionising photons (most likely a hidden quasar nucleus) and the power carried by the radio jets. We argue that radio sources modify their environments and that the range of emission-line properties seen is determined more by the range of source age than by the range in ambient environment. This is in accord with the idea that all high-redshift, high-luminosity radio sources are triggered in similar environments, presumably recently collapsed massive structures.

In this paper, the third and final of a series, we present complete K-band imaging and some complementary I-band imaging of the filtered 6C* sample. We find no systematic differences between the K-z relation of 6C* radio galaxies and those from complete samples, so the near-infrared properties of luminous radio galaxies are not obviously biased by the additional 6C* radio selection criteria (steep spectral index and small angular size). The 6C* K-z data significantly improve delineation of the K-z relation for radio galaxies at high-redshift (z >2). Accounting for non-stellar contamination, and for correlations between radio luminosity and stellar mass, we find little support for previous claims that the underlying scatter in the stellar luminosity of radio galaxies increases significantly at z >2. In a particular spatially-flat universe with a cosmological constant, the most luminous radio sources appear to be associated with galaxies with a luminosity distribution with a high mean (~5 Lstar), and a low dispersion (sigma ~ 0.5 mag) which formed their stars at epochs corresponding to z >~2.5. This result is in line with recent sub-mm studies of high-redshift radio galaxies and the inferred ages of extremely red objects from faint radio samples.

At least half the radio galaxies at z>1 in the 7C Redshift Survey have extremely red colours (R-K>5), consistent with stellar populations which formed at high redshift (z>5). We discuss the implications of this for the evolution of massive galaxies in general and for the fraction of near-IR-selected EROs which host AGN, a result which is now being tested by deep, hard X-ray surveys. The conclusion is that many massive galaxies undergo at least two active phases: one at z~5 when the black hole and stellar bulge formed and another at z~1-2 when activity is triggered by an event such as an interaction or merger.

Using a new radio sample, 6C* designed to find radio galaxies at z > 4 along with the complete 3CRR and 6CE sample we extend the radio galaxy K-z relation to z~4.5. The 6C* K-z data significantly improve delineation of the K-z relation for radio galaxies at high redshift (z > 2). Accounting for non-stellar contamination, and for correlations between radio luminosity and estimates of stellar mass, we find little support for previous claims that the underlying scatter in the stellar luminosity of radio galaxies increases significantly at z > 2. This indicates that we are not probing into the formation epoch until at least z > 3.

We present the results of the first major systematic submillimetre survey of radio galaxies spanning the redshift range 1 < z < 5. The primary aim of this work is to elucidate the star-formation history of this sub-class of elliptical galaxies by tracing the cosmological evolution of dust mass. Using SCUBA on the JCMT we have obtained 850-micron photometry of 47 radio galaxies to a consistent rms depth of 1 mJy, and have detected dust emission in 14 cases. The radio galaxy targets have been selected from a series of low-frequency radio surveys of increasing depth (3CRR, 6CE, etc), in order to allow us to separate the effects of increasing redshift and increasing radio power on submillimetre luminosity. Although the dynamic range of our study is inevitably small, we find clear evidence that the typical submillimetre luminosity (and hence dust mass) of a powerful radio galaxy is a strongly increasing function of redshift; the detection rate rises from 15 per cent at z < 2.5 to 75 per cent at z > 2.5, and the average submillimetre luminosity rises as (1+z)^3 out to z~4. Moreover our extensive sample allows us to argue that this behaviour is not driven by underlying correlations with other radio galaxy properties such as radio power, radio spectral index, or radio source size/age. Although radio selection may introduce other more subtle biases, the redshift distribution of our detected objects is in fact consistent with the most recent estimates of the redshift distribution of comparably bright submillimetre sources discovered in blank field surveys. The evolution of submillimetre luminosity found here for radio galaxies may thus be representative of massive ellipticals in general.

We present RIJHK imaging of seven radio galaxies from the 7C Redshift Survey (7CRS) which lack strong emission lines and we use these data to investigate their spectral energy distributions (SEDs) with models which constrain their redshifts. Six of these seven galaxies have extremely red colours (R-K>5.5) and we find that almost all of them lie in the redshift range 1<z<2. We also present near-infrared spectroscopy of these galaxies which demonstrate that their SEDs are not dominated by emission lines, although tentative lines, consistent with H-alpha at z=1.45 and z=1.61, are found in two objects. Although the red colours of the 7CRS galaxies can formally be explained by stellar populations which are either very old or young and heavily reddened, independent evidence favours the former hypothesis. At z~1.5 at least 1/4 of powerful radio jets are triggered in massive (>L*) galaxies which formed the bulk of their stars several Gyr earlier, that is at epochs corresponding to redshifts z>5. If a similar fraction of all z~1.5 radio galaxies are old, then extrapolation of the radio luminosity function shows that, depending on the radio source lifetimes, between 10-100% of the near-IR selected extremely red object (ERO) population undergo a radio outburst at epochs corresponding to 1<z<2. An ERO found serendipitously in the field of one of the 7CRS radio sources appears to be a radio-quiet analogue of the 7CRS EROs with an emission line likely to be [OII] at z=1.20. The implication is that some of the most massive elliptical galaxies formed the bulk of their stars at z>5 and these objects probably undergo at least two periods of AGN activity: one at high redshift during which the black hole forms and another one at an epoch corresponding to z~1.5.

This is the second of two papers presenting basic observational data on the 6CE sample of extragalactic radio sources. It presents the results of optical spectroscopy which has yielded virtually complete redshift information for the 6CE sample: 56 of the 59 sample members have spectroscopic redshifts which, with the exception of seven cases, are secure. The redshift distribution N(z) is fairly flat over the redshift range 0 < z < 2, with a median redshift of 1.1 and a high-redshift tail reaching to z = 3.4. The highest-redshift (z > 1.75) members of the 6CE sample have similar optical spectra and a tight (less than one dex) spread in narrow-Lyalpha emission line luminosity.

We measure the radio luminosity function (RLF) of steep-spectrum radio sources using three redshift surveys of flux-limited samples selected at low (151 & 178 MHz) radio frequency, low-frequency source counts and the local RLF. The redshift surveys used are the new 7C Redshift Survey (7CRS) and the brighter 3CRR and 6CE surveys totalling 356 sources with virtually complete redshift information. This yields unprecedented coverage of the radio luminosity versus z plane for steep-spectrum sources, and hence the most accurate measurements of the steep-spectrum RLF yet made. We find that a simple dual-population model for the RLF fits the data well, requiring differential density evolution (with z) for the two populations. The low-luminosity population can be associated with radio galaxies with weak emission lines, and includes sources with both FRI and FRII radio structures; its comoving space density $\rho$ rises by about one dex between z~0 and z~1 but cannot yet be meaningfully constrained at higher redshifts. The high-luminosity population can be associated with FRII radio galaxies and quasars with strong emission lines; its $\rho$ rises by nearly three dex between z~0 and z~2. These results mirror the situation seen in X-ray and optically-selected AGN. The integrated radio luminosity density of the combination of the two populations is controlled by the value of $\rho$ at the low-luminosity end of the RLF of the high-luminosity population, a quantity which has been directly measured at z~1 by the 7CRS. We argue that robust determination of this quantity at higher redshifts requires a new redshift survey based on a large (~1000 source) sample about five times fainter than the 7CRS.

I review some recent cosmological studies based on redshift surveys of radio sources selected at low frequencies. The accretion rate onto the central black hole is identified as the basis of a crude physical division of the low-frequency population into two sub-populations: the first comprises `Eddington-tuned' (high accretion rate) quasars and their torus-hidden counterparts; the second comprises `starved quasars' like M87. There exist remarkable similarities between the shapes and evolutionary behaviours of the luminosity functions of radio sources and radio-quiet quasars; all luminous AGN are $\sim300-$times rarer at epochs corresponding to $z \sim 0$ than at $z \sim 2.5$. I argue that any evidence that quasars were intrinsically rarer at $z \sim 5$ than at $z \sim 2.5$ is as yet both tentative and indirect. A simple calculation suggests that the radio source population has been over-looked as a potentially important contributor to the entropy budget of the Universe. A recent sub-mm survey of radio sources is used to demonstrate a connection between the events which trigger jets and intense bursts of star formation, and a close link between the histories of star formation and AGN activity is proposed. I discuss the aims and methods of future large redshift surveys of radio sources, emphasising the importance of dovetailing these with the development of robust physical models for radio sources and their epoch-dependent environments.

We use data from the Parkes Half-Jansky Flat-Spectrum (PHJFS) sample (Drinkwater et al. 1997) to constrain the cosmic evolution in the co-moving space density of radio sources in the top decade of the flat-spectrum radio luminosity function (RLF). A consistent picture for the high-redshift evolution is achieved using both simple parametric models, which are the first to allow for distributions in both radio luminosity and spectral index, and variants of the V / V_max test, some of which incorporate the effects of radio spectral curvature. For the most luminous flat-spectrum objects, the PHJFS sample is extremely similar to that used by Shaver et al. (1996, 1998) to argue for an abrupt `redshift cut-off': a decrease by a factor ~30 in space density between a peak redshift z~2.5 and z~5. Our analysis finds that the observable co-moving volume is too small to make definitive statements about any redshift cut-off for the most luminous flat-spectrum sources, although both constant-space density (no cut-off) models and models with cut-offs as abrupt as those envisaged by Shaver et al. are outside the 90% confidence region. The inference that the decline in space density is most likely to be gradual, by a factor ~4 between z~2.5 and z~5, is in accordance with previous work on the RLF by Dunlop & Peacock (1990), but different to the abrupt decline favoured by studies of optically-selected quasars. Dust obscuration provides one explanation for this difference. We show that a significant fraction of the most radio-luminous flat-spectrum objects are Giga-Hertz Peaked Spectrum (GPS) rather than Doppler-Boosted (DB) sources, complicating any interpretation of the redshift cut-off.

We present JHKL' photometry of a complete sample of steep-spectrum radio-loud quasars from the revised 3CR catalogue in the redshift range 0.65 < z < 1.20. After correcting for contributions from emission lines and the host galaxies, we investigate their spectral energy distributions (SEDs) around 1 micron. About 75% of the quasars are tightly grouped in the plane of optical spectral index, alpha_lo, versus near-infrared spectral index, alpha_hi, with the median value of alpha_lo close to the canonical value, and the median alpha_hi slightly flatter. We conclude that the fraction of moderately-obscured, red quasars decreases with increasing radio power, in accordance with the `receding torus' model which can also explain the relatively flat median near-infrared spectra of the 3CR quasars. Two of the red quasars have inverted infrared spectral indices, and we suggest that their unusual SEDs might result from a combination of dust-scattered and transmitted quasar light.

Low-frequency radio surveys are ideal for selecting orientation-independent samples of extragalactic sources because the sample members are selected by virtue of their isotropic steep-spectrum extended emission. We use the new 7C Redshift Survey along with the brighter 3CRR and 6C samples to investigate the fraction of objects with observed broad emission lines - the `quasar fraction' - as a function of redshift and of radio and narrow emission line luminosity. We find that the quasar fraction is more strongly dependent upon luminosity (both narrow line and radio) than it is on redshift. Above a narrow [OII] emission line luminosity of log L_[OII] > 35 W (or radio luminosity log L_151 > 26.5 W/Hz/sr), the quasar fraction is virtually independent of redshift and luminosity; this is consistent with a simple unified scheme with an obscuring torus with a half-opening angle theta_trans approx 53 degrees. For objects with less luminous narrow lines, the quasar fraction is lower. We show that this is not due to the difficulty of detecting lower-luminosity broad emission lines in a less luminous, but otherwise similar, quasar population. We discuss evidence which supports at least two probable physical causes for the drop in quasar fraction at low luminosity: (i) a gradual decrease in theta_trans and/or a gradual increase in the fraction of lightly-reddened (0 < A(V) < 5) lines-of-sight with decreasing quasar luminosity; and (ii) the emergence of a distinct second population of low luminosity radio sources which, like M87, lack a well-fed quasar nucleus and may well lack a thick obscuring torus.

We compare two temporal properties of classical double radio sources: i) radiative lifetimes of synchrotron-emitting particles and ii) dynamical source ages. We discuss how these can be quite discrepant from one another, rendering use of the traditional spectral ageing method inappropriate: we contend that spectral ages give meaningful estimates of dynamical ages only when these ages are << 10^7 years. In juxtaposing the fleeting radiative lifetimes with source ages which are significantly longer, a refinement of the paradigm for radio source evolution is required. The changing spectra along lobes are explained, not predominantly by synchrotron ageing but, by gentle gradients in a magnetic field mediated by a low-gamma matrix which illuminates an energy-distribution of particles, controlled largely by classical synchrotron loss in the high magnetic field of the hotspot. The energy in the particles is an order of magnitude higher than that inferred from the minimum-energy estimate, implying that the jet-power is of the same order as the accretion luminosity produced by the quasar central engine. This refined paradigm points to a resolution of the findings of Rudnick et al (1994) and Katz-Stone & Rudnick (1994) that both the Jaffe-Perola and Kardashev-Pacholczyk model spectra are invariably poor descriptions of the curved spectral shape of lobe emission, and indeed that for Cygnus A all regions of the lobes are characterised by a `universal spectrum'. [abridged]

We present near-infrared spectroscopy and imaging of the compact steep- spectrum radio source 3C 318 which shows it to be a quasar at redshift z=1.574 (the z=0.752 value previously reported is incorrect). 3C 318 is an IRAS, ISO and SCUBA source so its new redshift makes it the most intrinsically luminous far-infrared (FIR) source in the 3C catalogue (there is no evidence of strong gravitational lensing effects). Its bolometric luminosity greatly exceeds the 10^13 solar luminosity level above which an object is said to be hyperluminous. Its spectral energy distribution (SED) requires that the quasar heats the dust responsible for the FIR flux, as is believed to be the case in other hyperluminous galaxies, and contributes (at the >10% level) to the heating of the dust responsible for the sub-mm emission. We cannot determine whether a starburst makes an important contribution to the heating of the coolest dust, so evidence for a high star-formation rate (SFR) is circumstantial being based on the high dust mass required by its sub-mm detection. We show that the current sub-mm and FIR data available for high-z radio galaxies are consistent with SEDs similar to that of 3C 318. This indicates that at least some of this population may be detected in the sub-mm because of dust heated by the quasar nucleus, and that interpreting sub-mm detection as evidence for very high SFRs may not always be valid. We show that the 3C318 quasar is slightly reddened (A(V) ~ 0.5), the most likely cause of which is SMC-type dust in the host galaxy. If very distant radio galaxies are reddened in a similar way then we show that only slightly greater amounts of dust could obscure the quasars in these sources, thus explaining the low fraction of quasars amongst the very high redshift (z > 3) objects in low-frequency radio-selected samples.

We describe the pitfalls encountered in deducing from classical double radio source observables (luminosity, spectral index, redshift and linear size) the essential nature of how these objects evolve. We discuss the key role played by hotspots in governing the energy distribution of the lobes they feed, and subsequent spectral evolution. We present images obtained using the new 74 MHz receivers on the VLA and discuss constraints which these enforce on models of the backflow and ages in classical doubles.

We discuss a unifying explanation for many `trends with redshift' of radio galaxies which includes the relevance of their ages (time since their jet triggering event), and the marked dependence of their ages on redshift due to the selection effect of imposing a flux-limit. We briefly describe some important benefits which this `youth-redshift degeneracy' brings.

We use three samples (3CRR, 6CE and 6C*) to investigate the radio luminosity function (RLF) for the `most powerful' low-frequency selected radio sources. We find that the data are well fitted by a model with a constant co-moving space density at high redshift as well as by one with a declining co-moving space density above some particular redshift. This behaviour is very similar to that inferred for steep-spectrum radio quasars by Willott et al (1998) in line with the expectations of Unified Schemes. We conclude that there is as yet no evidence for a `redshift cut-off' in the co-moving space densities of powerful classical double radio sources, and rule out a cut-off at z < 2.5.

We describe the inter-dependence of four properties of classical double radio sources - spectral index, linear size, luminosity and redshift - from an extensive study based on spectroscopically-identified complete samples. We use these relationships to discuss aspects of strategies for searching for radio galaxies at extreme redshifts, in the context of possible capabilities of the new generation of proposed radio telescopes.

We have used narrow emission line data from the new 7C Redshift Survey to investigate correlations between the narrow-line luminosities and the radio properties of radio galaxies and steep-spectrum quasars. The 7C Redshift Survey is a low-frequency (151 MHz) selected sample with a flux-density limit about 25-times fainter than the 3CRR sample. By combining these samples, we can for the first time distinguish whether the correlations present are controlled by 151 MHz radio luminosity L_151 or redshift z. We find unequivocal evidence that the dominant effect is a strong positive correlation between narrow line luminosity L_NLR and L_151, of the form L_NLR proportional to L_151 ^ 0.79 +/- 0.04. Correlations of L_NLR with redshift or radio properties, such as linear size or 151 MHz (rest-frame) spectral index, are either much weaker or absent. We use simple assumptions to estimate the total bulk kinetic power Q of the jets in FRII radio sources, and confirm the underlying proportionality between jet power and narrow line luminosity first discussed by Rawlings & Saunders (1991). We make the assumption that the main energy input to the narrow line region is photoionisation by the quasar accretion disc, and relate Q to the disc luminosity, Q_phot. We find that 0.05 < Q / Q_phot < 1 so that the jet power is within about an order of magnitude of the accretion disc luminosity. The most powerful radio sources are accreting at rates close to the Eddington limit of supermassive black holes (~ 10^9 - 10^10 solar masses), whilst lower power sources are accreting at sub-Eddington rates.

We present the results of optical spectroscopy of two flux-density-limited samples of radio sources selected at frequencies of 38 and 151 MHz in the same region around the North Ecliptic Cap, the 8C-NEC and 7C-III samples respectively. Both samples are selected at flux density levels ~20 times fainter than samples based on the 3C catalogue. They are amongst the first low-frequency selected samples with no spectral or angular size selection for which almost complete redshift information has been obtained and they will therefore provide a valuable resource for understanding the cosmic evolution of radio sources and their hosts and environments. The 151-MHz 7C-III sample is selected to have S_151 >=0.5 Jy and is the more spectroscopically complete; out of 54 radio sources fairly reliable redshifts have been obtained for 44 objects. The 8C sample has a flux limit of S_38 >=1.3 Jy and contains 58 sources of which 46 have fairly reliable redshifts. We discuss possible biases in the observed redshift distribution, and some interesting individual objects. Using the 8C-NEC and 7C-III samples in conjunction, we form the first sample selected on low-frequency flux in the rest-frame of the source, rather than the usual selection on flux density in the observed frame. This allows us to remove the bias associated with an increasing rest-frame selection frequency with redshift. We investigate the difference this selection makes to correlations of radio source properties with redshift and luminosity. We show in particular that flux-density-based selection leads to an overestimate of the steepness of the correlation of radio source size with redshift. (abridged)

Further imaging observations of a sample of radio sources in the North Ecliptic Cap are presented and a number of new identifications are made. Using redshifts from spectroscopic data presented in a companion paper (Lacy et al.\ 1999b), the photometric properties of the galaxies in the sample are discussed. It is shown that: (1) out to at least z~0.6 radio galaxies are good standard candles irrespective of radio luminosity; (2) for 0.6~<z~<1 a large fraction of the sample have magnitudes and colours consistent with a non-evolving giant elliptical, and (3) at higher redshifts, where the R-band samples the rest-frame UV flux, most objects have less UV luminosity than expected if they form their stellar populations at a constant rate from a high redshift to $z\sim 1$ in unobscured star-forming regions (assuming an Einstein -- de Sitter cosmology). The consequences of these observations are briefly discussed.

Radio galaxies can be seen out to very high redshifts, where in principle they can serve as probes of the early evolution of the Universe. Here we show that for any model of radio-galaxy evolution in which the luminosity decreases with time after an initial rapid increase (that is, essentially all reasonable models), all observable high-redshift radio-galaxies must be seen when the lobes are less than 10^7 years old. This means that high-redshift radio galaxies can be used as a high-time-resolution probe of evolution in the early Universe. Moreover, this result helps to explain many observed trends of radio-galaxy properties with redshift [(i) the `alignment effect' of optical emission along radio-jet axes, (ii) the increased distortion in radio structure, (iii) the decrease in physical sizes, (iv) the increase in radio depolarisation, and (v) the increase in dust emission] without needing to invoke explanations based on cosmology or strong evolution of the surrounding intergalactic medium with cosmic time, thereby avoiding conflict with current theories of structure formation.

We present an optically-based study of the alignment between the radio axes and the optical major axes of eight z~0.7 radio galaxies in a 7C sample. The radio galaxies in this sample are ~20-times less radio luminous than 3C galaxies at the same redshift, and are significantly less radio-luminous than any other well-defined samples studied to date. Using Nordic Optical Telescope images taken in good seeing conditions at rest-frame wavelengths just longward of the 4000A break, we find a statistically significant alignment effect in the 7C sample. Furthermore, in two cases where the aligned components are well separated from the host we have been able to confirm spectroscopically that they are indeed at the same redshift as the radio galaxy. However, a quantitative analysis of the alignment in this sample and in a corresponding 3C sample from HST archival data indicates that the percentage of aligned flux may be lower and of smaller spatial scale in the 7C sample. Our study suggests that alignments on the 50-kpc scale are probably closely related to the radio luminosity, whereas those on the 15 kpc scale are not. We discuss these results in the context of popular models for the alignment effect.

We present the results of a programme of thermal-IR imaging of nineteen z~1 radio galaxies from the 3CR and 3CRR samples. We detect emission at L' (3.8um) from four objects; in each case the emission is unresolved at 1" resolution. Fifteen radio galaxies remain undetected to sensitive limits of L'~15.5. Using these data in tandem with archived HST data and near-IR spectroscopy we show that 3 of the detected `radio galaxies' (3C22, 3C41, and 3C65) harbour quasars reddened by Av<5. Correcting for this reddening 3C22 and 3C41 are very similar to coeval 3C quasars, whilst 3C65 seems unusually underluminous. The fourth radio galaxy detection (3C265) is a more highly obscured (Av~15) but otherwise typical quasar which previously has been evident only in scattered light. We determine the fraction of dust-reddened quasars at z~1 to be 28(+25)(-13)% at 90% confidence. On the assumption that the undetected radio galaxies harbour quasars similar to those in 3C22, 3C41 and 3C265 (as seems reasonable given their similar narrow emission line luminosities) we deduce extinctions of Av>15 towards their nuclei. The contributions of reddened quasar nuclei to the total K-band light ranges from ~0 per cent for the non-detections, through ~10 per cent for 3C265 to ~80 per cent for 3C22 and 3C41. Correcting for these effects does not remove the previously reported differences between the K magnitudes of 3C and 6C radio galaxies, so contamination by reddened quasar nuclei is not a serious problem for drawing cosmological conclusions from the K-z relation for radio galaxies. We discuss these results in the context of the `receding torus' model which predicts a small fraction of lightly-reddened quasars in samples of high radio luminosity sources. We also examine the likely future importance of thermal-IR imaging in the study of distant powerful radio sources.

We describe the selection of candidate radio-loud quasars obtained by cross-matching radio source positions from the low-frequency (151 MHz) 7C survey with optical positions from five pairs of EO POSS-I plates scanned with the Cambridge Automatic Plate-measuring Machine (APM). The sky region studied is centred at RA 10 h 28 m, Dec +41 and covers 0.057 sr. We present VLA observations of the quasar candidates, and tabulate various properties derived from the radio maps. We discuss the selection criteria of the resulting `7CQ' sample of radio-loud quasars. The 70 confirmed quasars, and some fraction of the 36 unconfirmed candidates, constitute a filtered sample with the following selection criteria: 151-MHz flux density S151 > 100 mJy; POSS-I E-plate magnitude E approx R < 20; and POSS-I colour (O - E) < 1.8; the effective area of the survey drops significantly below S151 approx 200 mJy. We argue that the colour criterion excludes few if any quasars, but note, on the basis of recent work by Willott et al. (1998b), that the E magnitude limit probably excludes more than 50 per cent of the radio-loud quasars.

We present deep optical images of the PC1643+4631 field obtained at the WHT. This field contains two quasars at redshifts z=3.79 & 3.83 and a cosmic microwave background (CMB) decrement detected with the Ryle Telescope. The images are in U,G,V,R and I filters, and are complete to 25th magnitude in R and G and to 25.5 in U. The isophotal galaxy counts are consistent with the results of Metcalde et al. (1996), Hogg et al. (1997), and others. We find an excess of robust high-redshift Ly-break galaxy candidates with 25.0<R<25.5 compared with the mean number found in the fields studied by Steidel et al. -we expect 7 but find 16 - but we do not find that the galaxies are concentrated in the direction of the CMB decrement. However, we are still not sure of the distance to the system causing the CMB decrement. We have also used our images to compare the commonly used object-finding algorithms of FOCAS and SExtractor: we find FOCAS the more efficient at detecting faint objects and the better at dealing with composite objects, whereas SExtractor's morphological classification is more reliable, especially for faint objects near the resolution limit. More generally, we have also compared the flux lost using isophotal apertures on a real image with that on a noise-only image: recovery of artificial galaxies from the noise-only image significantly overestimates the flux lost from the galaxies, and we find that the corrections made using this technique suffer a systematic error of some 0.4 magnitudes.

Deep radio observations at 5 and 8 GHz are presented of a complete sample of 11 radio galaxies with redshifts 0.85 < z < 1.5, selected from the 6C sample of Eales. The radio data, taken using the VLA in A, B and C array configurations, provide a best angular resolution of 0.25'' and reach an rms noise level of order 20 microJy. Radio core candidates are detected in 8 of the 11 sources. 9 of the 11 sources display deviations from `standard double radio source' morphologies, with multiple hotspots or a hotspot withdrawn from the leading edge of the radio emission. The sources are typically polarised at the 5 to 15% level. The mean rotation measures of the individual lobes are less than 50 rad/m^2 but strong asymmetries between the two lobes indicate that the Faraday rotation does not have a Milky Way origin; rather, the distant 6C radio sources lie in a relatively dense clumpy environment. The sources are compared with the more radio powerful distant 3CR radio galaxies and with low redshift sources. The ratio of core to extended radio flux is almost independent of the radio source size and only weakly inversely correlated with the total radio power. This indicates that the high radio luminosity of the most powerful radio sources must originate in the AGN, in contrast to the suggestion that they are so luminous only due to confinement by a dense surrounding environment. Environmental effects play a secondary role. (Abridged)

A systematic survey measuring the submillimetre continuum luminosity in radio galaxies between redshifts z ~ 0.1 - 5 is currently in progress. The first results from observations with the bolometer array SCUBA on the JCMT suggest a trend of increasing submillimetre luminosity with redshift out to z ~ 4. Assuming the continuum emission at 850um is dominated by thermal radiation from dust heated by young, massive stars, the straightforward interpretation of these data implies that the host galaxies of powerful radio sources, presumably ellipticals or their progenitors, exhibit increased star formation activity with increasing redshift. However a severe bias may distort the true picture since only a subset of 30 galaxies have been observed, which represent the most luminous radio sources (P_(151 MHz) > 10^27.5 W Hz^-1 sr^-1), whilst the complete sample covers ~4 decades in radio power (10^25 - 10^29 W Hz^-1 sr^-1). This on-going observational programme continues to improve coverage of the radio power-redshift plane in an attempt to remove this potential bias towards higher AGN activity and hence understand more fully the relationship between high-z AGN, their level of starformation, and the evolutionary status of their host galaxies.

We present a study of the trends in luminosity, linear size, spectral index, and redshift of classical double radio sources from three complete samples selected at successively fainter low radio-frequency flux-limits. We have been able to decouple the effects of the tight correlation between redshift and luminosity (inherent in any single flux-limited sample) which have hitherto hindered interpretation of the relationships between these four source properties. The major trends found are that (i) spectral indices increase with linear size, (ii) rest-frame spectral indices have a stronger dependence on luminosity than on redshift except at high (GHz) frequencies, and that (iii) the linear sizes are smaller at higher redshifts. We reproduce the observed dependences in a model for radio sources (born throughout cosmic time according to a radio-source birth function) whose lobes are fed with a synchrotron-emitting population (whose energy distribution is governed by compact hotspots), and which suffer inverse Compton, synchrotron and adiabatic expansion losses. In simulating the basic observed dependences, we find that there is no need to invoke any systematic change in the environments of these objects with redshift if the consequences of imposing a survey flux-limit on our simulated datasets are properly included in the model. We present evidence that for a radio survey there is an unavoidable `youth--redshift degeneracy', even though radio sources are short-lived relative to the age of the Universe; it is imperative to take this into account in studies which seemingly reveal correlations of source properties with redshift such as the `alignment effect'.

We present a complete sample of 24 radio-loud quasars (RLQs) from the new 7C Redshift Survey. Every quasar with a low-frequency (151 MHz) radio flux-density S_151 > 0.5 Jy in two regions of the sky covering 0.013 sr is included; 23 of these have sufficient extended flux to meet the selection criteria, 18 of these have steep radio spectra (hereafter denoted as SSQs). The key advantage of this sample over most samples of RLQs is the lack of an optical magnitude limit. By combining the 7C and 3CRR samples, we have investigated the properties of RLQs as a function of redshift z and radio luminosity L_151. We derive the radio luminosity function (RLF) of RLQs and find that the data are well fitted by a single power-law with slope alpha_1=1.9. We find that there must be a break in the RLQ RLF at log_10(L_151 / W Hz^-1 sr^-1) < 27, in order for the models to be consistent with the 7C and 6C source counts. The z-dependence of the RLF follows a one-tailed gaussian which peaks at z=1.7. We find no evidence for a decline in the co-moving space density of RLQs at higher redshifts. A positive correlation between the radio and optical luminosities of SSQs is observed, confirming a result of Serjeant et al. (1998). We are able to rule out this correlation being due to selection effects or biases in our combined sample. The radio-optical correlation and best-fit model RLF enable us to estimate the distribution of optical magnitudes of quasars in samples selected at low radio frequencies. We conclude that for samples with S_151 < 1 Jy one must use optical data significantly deeper than the POSS-I limit (R approx 20), in order to avoid severe incompleteness.

We use new near-infrared spectroscopy and our published optical spectroscopy of the gravitationally-lensed Seyfert-2 galaxy F10214+4724 to study both the links between the starburst and AGN in this object and the properties of the inner narrow-line clouds. The UV spectrum is consistent with a compact, moderately- reddened starburst providing about half the UV light. Spectroscopy of the Halpha /[NII] line blend has enabled us to distinguish emission from the narrow-line region of the Seyfert-2 and a moderately-reddened emission line region which we argue is associated with the starburst. Estimates of the star formation rate from the UV continuum flux and the Halpha flux are broadly consistent. We can explain the unusual emission line properties of F10214+4724 in terms of conventional models for nearby Seyfert-2 galaxies if lensing is preferentially magnifying the side of the inner narrow-line region between the AGN and the observer, and the other side is both less magnified and partially obscured by the torus. The hydrogen densities of clouds in this region are high enough to make the Balmer lines optically thick and to suppress forbidden emission lines with low critical densities. We have deduced the column density of both ionised and neutral gas in the narrow-line clouds, and the density of the ionised gas. Using these we have been able to estimate the mass of the inner narrow-line clouds to be ~ 1 solar mass, and show that the gas:dust ratio NH/E(B-V) in these clouds must be ~1.3x10^{27}m^{-2}mag^{-1}, significantly higher than in the Milky Way. The cloud properties are consistent with the those of the warm absorbers seen in the X-ray spectra of Seyfert-1 galaxies. Our results favour models in which narrow-line clouds start close to the nucleus and flow out.

Multi-wavelength imaging and spectroscopy of the z=0.708 radio galaxy 3C441 and a red aligned optical/infrared component are used to show that the most striking aspect of the radio-optical ``alignment effect'' in this object is due to the interaction of the radio jet with a companion galaxy in the same group or cluster. The stellar population of the red aligned continuum component is predominately old, but with a small post-starburst population superposed, and it is surrounded by a low surface-brightness halo, possibly a face-on spiral disc. The [OIII]500.7/[OII]372.7 emission line ratio changes dramatically from one side of the component to the other, with the low-ionisation material apparently having passed through the bow shock of the radio source and been compressed. A simple model for the interaction is used to explain the velocity shifts in the emission line gas, and to predict that the ISM of the interacting galaxy is likely to escape once the radio source bow shock has passed though. We also discuss another, much fainter, aligned component, and the sub-arcsecond scale alignment of the radio source host galaxy. Finally we comment on the implications of our explanation of 3C441 for theories of the alignment effect.

The z=2.286 IRAS galaxy F10214+4724 remains one of the most luminous galaxies in the Universe, despite its gravitational lens magnification. We present optical and near-infrared spectra of F10214+4724, with clear evidence for three distinct components: lines of width ~1000 km/s from a Seyfert-II nucleus; <~200 km/s lines which are likely to be associated with star formation; and a broad ~4000 km/s CIII] 1909ang emission line which is blue-shifted by ~1000 km/s with respect to the Seyfert-II lines. Our study of the Seyfert-II component leads to several new results, including: (i) From the double-peaked structure in the Ly alpha line, and the lack of Ly beta, we argue that the Ly alpha photons have emerged through a neutral column of N_H ~ 2.5 x 10^{25}/m^2, possibly located within the AGN narrow-line region as argued in several high redshift radiogalaxies. (ii) The resonant O VI 1032,1036ang doublet (previously identified as Ly beta) is in an optically thick (1:1) ratio. At face value this implies an an extreme density (n_e ~ 10^{17}/m^3) more typical of broad line region clouds. However, we attribute this instead to the damping wings of Ly beta from the resonant absorption. (iii) A tentative detection of HeII 1086 suggests little extinction in the rest-frame ultraviolet.

Although radio-quiet quasars (RQQs) constitute >90% of optically-identified quasar samples their radio properties are only poorly understood. In this paper we present the results of a multi-frequency VLA study of 27 low-redshift RQQs. We detect radio emission from 20 objects, half of which are unresolved (< 0.24''). In cases where significant structure can be resolved, double, triple and linear radio sources on scales of a few kpc are found. The radio emission (typically) has a steep spectrum (alpha ~ 0.7), and high brightness temperatures (T_B > 10^5 K) are measured in some of the radio components. The RQQs form a natural extension to the radio luminosity - absolute magnitude distribution of nearby Seyfert 1s. We conclude that a significant fraction of the radio emission in RQQs originates in a compact nuclear source directly associated with the quasar. There are no significant differences between the radio properties of RQQs with elliptical hosts and those in disc galaxies within the current sample.

We present preliminary results of a search for galaxies at z>4 through Lyman-limit imaging of the fields of known high-redshift radio-galaxies. Objects were selected by means of their broad-band colours, and spectroscopy of candidate objects in one of the fields has been performed through multi-slit spectroscopy at the 4.2m William Herschel Telescope. These spectra show some of the first z>4 galaxies to be identified using the Lyman break technique.

Using complete samples of steep-spectrum quasars, we present evidence for a correlation between radio and optical luminosity which is not caused by selection effects, nor caused by an orientation dependence (such as relativistic beaming), nor a byproduct of cosmic evolution. We argue that this rules out models of jet formation in which there are no parameters in common with the production of the optical continuum. This is arguably the most direct evidence to date for a close link between accretion onto a black hole and the fuelling of relativistic jets. The correlation also provides a natural explanation for the presence of aligned optical/radio structures in only the most radio luminous high-redshift galaxies.

We present the results of new rest-frame far-IR observations of the z = 4.25 radio galaxy, 8C1435+635, which not only confirm that it contains an enormous quantity of dust (as first inferred from its mm-wave detection by Ivison 1995) but also allow the first meaningful constraints to be placed on the mass of this dust and associated gas. The new measurements consist of: (i) clear detections of sub-mm continuum emission at 450 and 850um obtained with the new sub-mm bolometer array, SCUBA, on the JCMT, (ii) continuum upper limits at 350, 750 and 175um obtained with SCUBA and the PHT far-IR camera aboard ISO, and (iii) a sensitive upper limit on the CO(4-3) line flux obtained with the IRAM 30-m MRT. The resulting rest-frame 33- to 238-um continuum coverage allows us to deduce that 2E8 Msun of dust at a temperature of 40 +/- 5 K is responsible for the observed mm/sub-mm emission. Using our CO upper limit, which constrains M(H2)/M(dust) to <950, we go on to calculate robust limits on the total gas reserves (H2 + HI) which are thereby constrained to between 4E10 and 1.2E12 Msun. The sub-mm properties of 8C1435+635 are thus strikingly similar to those of the z = 3.80 radio galaxy, 4C41.17, the only other high-redshift galaxy detected to date at sub-mm wavelengths whose properties appear not to be exaggerated by gravitational lensing (Dunlop et al. 1994; Hughes, Dunlop & Rawlings 1997). The inferred gas masses of both objects are sufficiently large to suggest that the formative starbursts of massive elliptical galaxies are still in progress at around z = 4. Observations of complete samples of radio galaxies spanning a range of redshifts and radio luminosities will be required to determine if the spectacular far-IR properties of 8C1435+635 and 4C41.17 are primarily due to their extreme redshifts or their extreme radio luminosities.

We describe the selection of a sample of 34 radio sources from the 6C survey (Hales, Baldwin & Warner 1993) from a region of sky covering 0.133 sr. The selection criteria for this sample, hereafter called 6C*, were chosen to optimise the chances of finding radio galaxies at redshift z > 4. Optical follow-up observations have already led to the discovery of the most distant known radio galaxy at z = 4.41 (Rawlings et al. 1996). We present VLA radio maps and derive radio spectra for all the 6C* objects.

We use $K-band imaging to investigate the angular size and morphology of 10 6C radio galaxies, at redshifts $1\leq z\leq 1.4$. Two appear to be undergoing mergers, another contains two intensity peaks aligned with the radio jets, while the other seven appear consistent with being normal ellipticals in the K-band. Intrinsic half-light radii are estimated from the areas of each radio galaxy image above a series of thresholds. The 6C galaxy radii are found to be significantly smaller than those of the more radioluminous 3CR galaxies at similar redshifts. This would indicate that the higher mean K-band luminosity of the 3CR galaxies results from a difference in the size of the host galaxies, and not solely from a difference in the power of the active nuclei. The size-luminosity relation of the $z\sim 1.1$ 6C galaxies indicates a 1.0--1.8 mag enhancement of the rest-frame R-band surface brightness relative to either local ellipticals of the same size or FRII radio galaxies at $z<0.2$. The 3CR galaxies at $z\sim 1.1$ show a comparable enhancement in surface brightness. The mean radius of the 6C galaxies suggests that they evolve into ellipticals of $L\sim L^*$ luminosity, and is consistent with their low redshift counterparts being relatively small FRII galaxies, a factor $\sim 25$ lower in radio luminosity, or small FRI galaxies a factor of $\sim 1000$ lower in radio luminosity. Hence the 6C radio galaxies may undergo at least as much optical and radio evolution as the 3CR galaxies.

We present new results of a study of the sub-millimetre continuum emission from a sample of 9 radio galaxies and 4 radio-quiet quasars at redshifts $z = 0.75 - 4.26$. The observations were made at 800$\mu m$, using the single-element bolometer UKT14 on the James Clerk Maxwell Telescope (JCMT), reaching a typical r.m.s. sensitivity of $\sigma_{rms} \sim 4$ mJy and represent some of the deepest submillimetre extragalactic measurements made to date. Three detections were achieved, of which two are secure (4C41.17, Dunlop et al. (1994) and H1413$+$117) and one (53W002) is tentative, whilst comparable upper-limits were obtained for 7 of the 10 remaining sources. We use these data as the motivation for a detailed discussion of the conversion from submillimetre and millimetre continuum fluxes to dust/gas masses and star formation rates at high-redshift, and determine these quantities from our own and other data on high-redshift radio-galaxies and quasars. In particular we have investigated the impact of the four main sources of uncertainty in deriving physical quantities from such data, namely i) potential contamination by galactic cirrus, ii) uncertainty in the value of the dust rest-frequency mass-absorption coefficient, iii) difficulty in constraining the dust temperature, and iv) estimation of the appropriate gas:dust ratio in high-redshift objects. Our discussion emphasises how important it will be to quantify and, where possible, minimize such uncertainties (via, for example, appropriate observational strategies) in order to fully capitalize on the ten-fold improvement in sensitivity offered by the imminent arrival of the next generation of bolometer arrays, such as SCUBA on the JCMT.

We present details of a new sample of radio-loud quasars drawn from 0.013 sr of the 7C Redshift Survey. This sample is small (21 quasars) but complete in that every object with an unresolved nucleus and/or broad emission lines with S(151MHz) > 0.5 Jy has been discovered. The dependence of the quasar fraction with redshift and radio luminosity is investigated, providing new evidence supporting the unification of radio-loud quasars and powerful radio galaxies. This 7C sample is compared with optically-selected quasars, in order to determine whether there are systematic biases in the different selection techniques. There are no lightly reddened (Av approx. 1) quasars in our sample amongst the 14 with z < 2. The discovery of a reddened quasar at z = 2.034 and its implications are discussed. A tight correlation between radio luminosity and optical/near infrared continuum luminosity for a subset of the sample is also found.

We use the K-z relation for radio galaxies to illustrate why it has proved difficult to obtain definitive cosmological results from studies based entirely on catalogues of the brightest radio sources, e.g. 3C. To improve on this situation we have been undertaking redshift surveys of complete samples drawn from the fainter 6C and 7C radio catalogues. We describe these surveys, and illustrate the new studies they are allowing. We also discuss our `filtered' 6C redshift surveys: these have led to the discovery of a radio galaxy at z=4.4, and are sensitive to similar objects at higher redshift provided the space density of these objects is not declining too rapidly with z. There is currently no direct evidence for a sharp decline in the space density of radio galaxies for z > 4, a result only barely consistent with the observed decline of flat-spectrum radio quasars.

We present preliminary results from HST imaging of radio-loud quasar hosts, covering a ~x100 range in radio luminosity but in a narrow redshift range (0.5 to 0.65). The sample was selected from our new, spectroscopically complete 7C survey and the 3CRR catalogue. Despite the very large radio luminosity range, the host luminosities are only weakly correlated (if at all) with radio power, perhaps reflecting a predominance of purely central engine processes in the formation of radio jets, and hence perhaps also in the radio-loud/-quiet dichotomy at these redshifts. The results also contradict naive expectations from several quasar formation theories, but the host magnitudes support radio-loud Unified Schemes.

Although radio-quiet quasars (RQQs), which constitute the majority of optically-identified quasar samples, are by no means radio silent the properties of their radio emission are only poorly understood. We present the results of a multi-frequency VLA study of 27 low-redshift RQQs. In general, we find that the properties of the radio sources in RQQs are consistent with them being weak, small-scale (~1 kpc) jets similar to those observed in nearby Seyfert galaxies. We conclude that a significant fraction of the radio emission in RQQs is directly associated with the central engine and is not a result of stellar processes in the surrounding galaxy. There appears to be no difference between the radio properties of RQQs in elliptical and disc galaxies, implying that the relationship between the host galaxy and the `radio loudness' of the active nucleus is not straightforward.

We have obtained redshifts and infrared images for a sample of faint B2/6C radio sources whose fluxes are about six times fainter than those of sources in the bright 3C sample. We now have unambiguous redshifts for 90% of the sources, making this the first faint radio sample with such complete redshift information. We find that the infrared Hubble diagrams (K versus z) of the 3C sample and the B2/6C sample are similar at a low redshift (z < 0.6), but by a redshift about 1 the 6C/B2 galaxies are on average about 0.6 mags fainter in the K-band than 3C galaxies at the same redshift. This suggests that the bright K-magnitudes of 3C galaxies at a redshift of about 1 are not the result of stellar evolution, but of a correlation between infrared and radio luminosity. We also find that the infrared stuctures of B2/6C galaxies at z=1 are less clearly aligned with their radio structures than 3C galaxies at this redshift, implying that the strength of the alignment effect depends on radio luminosity. Finally, above a redshift of 2 we find that the dispersion in the K-z relation of the B2/6C sample is about 2.7 times greater than that at low redshift, a result which is expected if at these redshifts we are beginning to penetrate into the epoch in which radio galaxies formed.

In a companion paper (Paper I: Jones et al. 1996) we report the discovery of a cosmic microwave background decrement, indicative of a distant cluster with mass ~10^{15} solar masses, towards the quasar pair PC1643+4631A&B (z=3.79,3.83, separation 198''). To search for the cluster responsible, we have obtained R-, J- and K-band images of the field and have also carried out optical spectroscopy of selected objects in it. No such cluster is evident in these images. Assuming the cluster causing the decrement is similar to massive clusters already known, our magnitude limits imply that it must lie about or beyond z=1. This provides independent support for the X-ray-based distance argument of Paper I. The cluster must gravitationally lens objects behind it; for a cluster z around 1-2, the Einstein-ring radius for sources at z ~= 3.8 is ~ 100''. Simple modelling, producing simultaneously the S-Z effect and the lensing, shows that the source positions of quasars A and B lie within ~10'' of each other and may indeed be coincident. The two quasar spectra are found to be remarkably similar apart from their one-percent redshift difference. Assuming A and B are images of a single quasar, we present a possible explanation of this difference.

Recent investigations of how the median size of extragalactic radio sources change with redshift have produced inconsistent results. Eales compared the radio and optical properties of a bright 3C and faint 6C sample and concluded that $D\propto(1+z)^{-1.1\pm0.5}$ ($\Omega_0 = 0$), with $D$ being the median size of the radio sources at a given epoch and z the redshift. Oort, Katgert, and Windhorst, on the other hand, from a comparison of the properties of a number of radio samples, found much stronger evolution, with $D\propto(1+z)^{-3.3 \pm0.5}$. In this paper we attempt to resolve the difference. We have repeated the analysis of Eales using the virtually complete redshift information that now exists for the 6C sample. Confining our analysis to FR2 sources, which we argue is the best-understood class of radio sources and the least likely to be affected by selection effects, we find $D\propto(1+z)^{-1.2\pm0.5}$ ($\Omega_0 = 0$) and $D\propto(1+z)^{-1.7\pm0.4}$ ($\Omega_0 = 1$). Our complete redshift information allows us to gain insight into our result by plotting a radio luminosity-size (P-D) diagram for the 6C sample. The most obvious difference between the 3C and 6C P-D diagrams is the clump of sources in the 6C diagram at $D\sim 100 kpc, P_{151}\sim 5x10^{27} WHz^{-1}sr^{-1}$. These clump sources have similar sizes to the emission-line regions found around high-redshift radio galaxies, suggesting that the presence of dense line-emitting gas around high-redshift radio galaxies is responsible for the size evolution. We show that this explanation can quantitatively explain the observed size evolution, as long as there is either little X-ray emitting gas around these objects or, if there is, it is distributed in a similar way to the emission-line gas: highly anisotropic and inhomogeneous.

The extreme bolometric luminosity of IRAS F10214+4724, and in particular the huge mass in molecular gas inferred from the CO line fluxes have led to suggestions that this is a giant galaxy in the process of formation. An arc-like structure and the closeness of a second object suggest, however, that gravitational lensing may be responsible for its anomalously high luminosity and mass. In this paper, we use an optical spectrum taken in conditions of 0.7-arcsec seeing to provide further evidence that F10214+4724 is a gravitationally lensed system. We measure tentative redshifts of 0.896 and 0.899 for galaxies projected $\approx 1$ and $\approx 3$ arcsec from IRAS F10214+4724 respectively. Identifying the former as the lensing galaxy we obtain a mass:light ratio consistent with those derived for other lenses, and find that its luminosity is consistent with the velocity dispersion deduced from the radius of the Einstein ring. If lensed, our models suggest magnification by a factor $\sim 10$, and hence an intrinsic bolometric luminosity for F10214+4724 similar to local ULIRGs.