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

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

We present a new redshift survey, the 2dF Quasar Dark Energy Survey pilot (2QDESp), which consists of ${\approx}10000$ quasars from ${\approx}150$ deg$^2$ of the southern sky, based on VST-ATLAS imaging and 2dF/AAOmega spectroscopy. Combining our optical photometry with the WISE (W1,W2) bands we can select essentially contamination free quasar samples with $0.8{<}z{<}2.5$ and $g{<}20.5$. At fainter magnitudes, optical UVX selection is still required to reach our $g{\approx}22.5$ limit. Using both these techniques we observed quasar redshifts at sky densities up to $90$ deg$^{-2}$. By comparing 2QDESp with other surveys (SDSS, 2QZ and 2SLAQ) we find that quasar clustering is approximately luminosity independent, with results for all four surveys consistent with a correlation scale of $r_{0}{=}6.1{\pm}0.1 \: h^{-1}$Mpc, despite their decade range in luminosity. We find a significant redshift dependence of clustering, particularly when BOSS data with $r_{0}{=}7.3{\pm}0.1 \: h^{-1}$Mpc are included at $z{\approx}2.4$. All quasars remain consistent with having a single host halo mass of ${\approx}2{\pm}1{\times}10^{12} \: h^{-1}M_\odot$. This result implies that either quasars do not radiate at a fixed fraction of the Eddington luminosity or AGN black hole and dark matter halo masses are weakly correlated. No significant evidence is found to support fainter, X-ray selected quasars at low redshift having larger halo masses as predicted by the `hot halo' mode AGN model of Fanidakis et al. 2013. Finally, although the combined quasar sample reaches an effective volume as large as that of the original SDSS LRG sample, we do not detect the BAO feature in these data.

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

We apply a cross-correlation technique to infer the $S>3$mJy radio luminosity function (RLF) from the NRAO VLA sky survey (NVSS) to $z\sim3.5$. We measure $\Sigma$ the over density of radio sources around spectroscopically confirmed quasars. $\Sigma$ is related to the space density of radio sources at the distance of the quasars and the clustering strength between the two samples, hence knowledge of one constrains the other. Under simple assumptions we find $\Phi\propto (1+z)^{3.7\pm0.7}$ out to $z\sim2$. Above this redshift the evolution slows and we constrain the evolution exponent to $<1.01$ ($2\sigma$). This behaviour is almost identical to that found by previous authors for the bright end of the RLF potentially indicating that we are looking at the same population. This suggests that the NVSS is dominated by a single population; most likely radio sources associated with high-excitation cold-mode accretion. Inversely, by adopting a previously modelled RLF we can constrain the clustering of high-redshift radio sources and find a clustering strength consistent with $r_0=15.0\pm 2.5$ Mpc up to $z\sim3.5$. This is inconsistent with quasars at low redshift and some measurements of the clustering of bright FRII sources. This behaviour is more consistent with the clustering of lower luminosity radio galaxies in the local universe. Our results indicate that the high-excitation systems dominating our sample are hosted in the most massive galaxies at all redshifts sampled.

Over the past 20years reverberation mapping has proved one of the most successful techniques for studying the local (<1pc) environment of super-massive black holes that drive active galactic nuclei. Key successes of reverberation mapping have been direct black hole mass estimates, the radius-luminosity relation for the Hbeta line and the calibration of single-epoch mass estimators commonly employed up to z~7. However, observing constraints mean that few studies have been successful at z>0.1, or for the more-luminous quasars that make up the majority of current spectroscopic samples, or for the rest-frame ultra-violet emission lines available in optical spectra of z>0.5 objects. Previously we described a technique for stacking cross correlations to obtain reverberation mapping results at high z. Here we present the first results from a campaign designed for this purpose. We construct stacked cross-correlation functions for the CIV and MgII lines and find a clear peak in both. We find the peak in the MgII correlation is at longer lags than CIV consistent with previous results at low redshift. For the CIV sample we are able to bin by luminosity and find evidence for increasing lags for more-luminous objects. This CIV radius-luminosity relation is consistent with previous studies but with a fraction of the observational cost.

Using a sample of 24264 optically selected AGNs from the SDSS DR7 database, we characterise how the profile of the [OIII] emission line relates to bolometric luminosity (L_Bol), Eddington ratio, radio loudness, radio luminosity (L_Rad) and optical class (Sy 1/2) to determine what drives the kinematics of this kpc-scale line emitting gas. Spectral stacking is used to characterise how the average [OIII] profile changes as a function of these variables. After accounting for the correlation between L_Bol and L_Rad, we report that L_Rad has the strongest influence on the [OIII] profile, with moderate radio luminosity AGNs (log(L_Rad)=23-25 W/Hz) having the broadest [OIII] profiles. When binned according to Eddington ratio, only AGNs in our highest bin (i.e., R_Edd>0.3) show any signs of having broadened [OIII] profiles, although the small numbers of such extreme AGNs mean we cannot rule out that other processes (e.g., radio jets) are responsible for this broadening. We find no significant difference between the [OIII] profiles of broad and narrow-line Seyfert 1s. The [OIII] profiles of Sy 1 and 2 AGNs show the same trends in terms of line width, but Sy 1 AGNs display a much stronger blue wing, which we interpret as evidence of outflowing ionised gas. Using multi-component line fitting we calculate the proportions of AGNs with broad [OIII] profiles, finding that moderate radio luminosity AGNs are roughly 5 times more likely to have extremely broad [OIII] lines (FWHM>1000km/s) compared to lower L_Rad AGNs. Our results are consistent with the most disturbed gas kinematics being induced by compact radio cores (rather than powerful radio jets), although broadened [OIII] lines are also present, but much rarer, in low L_Rad systems (<10^23 W/Hz). Our catalogue of multi-component fits is freely available as an online resource.

Reverberation mapping offers one of the best techniques for studying the inner regions of QSOs. It is based on cross-correlating continuum and emission-line light curves. New time-resolved optical surveys will produce well sampled light curves for many thousands of QSOs. We explore the potential of stacking samples to produce composite cross-correlations for groups of objects that have well sampled continuum light curves, but only a few (~2) emission-line measurements. This technique exploits current and future wide-field optical monitoring surveys (e.g. Pan-STARRS, LSST) and the multiplexing capability of multi-object spectrographs (e.g. 2dF, Hectospec) to significantly reduce the observational expense of reverberation mapping, in particular at high redshift (0.5 to 2.5). We demonstrate the technique using simulated QSO light curves and explore the biases involved when stacking cross-correlations in some simplified situations. We show that stacked cross correlations have smaller amplitude peaks compared to well sampled correlation functions as the mean flux of the emission light curve is poorly constrained. However, the position of the peak remains intact. We find there can be `kinks' in stacked correlation functions due to different measurements contributing to different parts of the correlation function. Using the Pan-STARRS Medium-Deep Survey (MDS) as a template we show that cross-correlation lags should be measurable in a sample size of 500 QSOs that have weekly photometric monitoring and two spectroscopic observations. Finally we apply the technique to a small sample (42) of QSOs that have light curves from the MDS. We find no indication of a peak in the stacked cross-correlation. A larger spectroscopic sample is required to produce robust reverberation lags.

Our goals are (i) to search for BAO and large-scale structure in current QSO survey data and (ii) to use these and simulation/forecast results to assess the science case for a new, >10x larger, QSO survey. We first combine the SDSS, 2QZ and 2SLAQ surveys to form a survey of ~60000 QSOs. We find a hint of a peak in the QSO 2-point correlation function, xi(s), at the same scale (~105h^-1 Mpc) as detected by Eisenstein et al (2005) in their sample of DR5 LRGs but only at low statistical significance. We then compare these data with QSO mock catalogues from the Hubble Volume simulation used by Hoyle et al (2002) and find that both routes give statistical error estimates that are consistent at ~100h^-1 Mpc scales. Mock catalogues are then used to estimate the nominal survey size needed for a 3-4 sigma detection of the BAO peak. We find that a redshift survey of ~250000 z<2.2 QSOs is required over ~3000 deg^2. This is further confirmed by static log-normal simulations where the BAO are clearly detectable in the QSO power spectrum and correlation function. The nominal survey would on its own produce the first detection of, for example, discontinuous dark energy evolution in the so far uncharted 1<z<2.2 redshift range. A survey with ~50% higher QSO sky densities and 50% bigger area will give an ~6sigma BAO detection, leading to an error ~60% of the size of the BOSS error on the dark energy evolution parameter, w_a. Another important aim for a QSO survey is to place new limits on primordial non-Gaussianity at large scales, testing tentative evidence we have found for the evolution of the linear form of the combined QSO xi(s) at z~1.6. Such a QSO survey will also determine the gravitational growth rate at z~1.6 via z-space distortions, allow lensing tomography via QSO magnification bias while also measuring the exact luminosity dependence of small-scale QSO clustering.

We cross match the NVSS and FIRST surveys with three large photometric catalogues of luminous red galaxies (LRGs) to define radio-loud samples. These have median redshifts 0.35, 0.55 and 0.68 and, by matching rest-frame optical and radio properties, we construct uniform samples across the three surveys. This paper is concerned with the clustering properties of these samples derived from the angular correlation function. The primary aim is to characterise any evolution in the clustering amplitude of radio galaxies bellow z~0.68. We find no evidence for evolution in the large-scale (~1-50h^{-1}Mpc) clustering amplitude. Our radio galaxy autocorrelations are consistent with previous findings indicating little-to-no evolution in the redshift range 0.68 to 0 (~6Gyr of time). We also cross correlate radio galaxies with the parent LRG samples to increase the precision of our results and again find no evidence for evolution. Our results are inconsistent with a long-lived model for the clustering evolution that assumes radio sources randomly sample the LRG population. A model where the halo mass is constant with redshift is consistent with the data. This is similar to QSOs that have clustering amplitudes consistent with a single halo mass at all redshifts. Given that the brightest radio sources show stronger evolution in space density compared to fainter radio sources we restrict our samples to include only objects with L>10^{26}W/Hz and repeat the analysis. Again we find no evidence for evolution in the comoving correlation amplitude. These radio sources appear to inhabit the same mass halos as fainter radio galaxies (~9x10^{13}h^{-1}M_\odot). These halos are ~twice as massive as those of the general LRG population and ~30 times as massive as optical AGN/QSOs.

QSOs from SDSS, 2QZ and 2SLAQ covering an order of magnitude in luminosity at fixed redshift exhibit similar amplitudes of clustering. In addition, QSO clustering evolution at z>0.5 is well fitted by a model that assumes a fixed host halo mass, implying that QSOs may occur in a relatively narrow range of halo and BH mass. We argue that the slow evolution of early-type galaxies out to z~1-2 may also provide support for a slow evolution of QSO host BH masses. The result would mean that if high-z QSOs radiate at Eddington rates then low-z SyI must radiate at ~100x less than Eddington. We conclude that models where QSOs radiate at L_Edd require M_BH and M_halo to be decoupled to circumvent the clustering results. While single BH mass and flickering models fit the z>0.5 clustering results, they appear to be rejected by the z~0, M_BH-L relation from reverberation mapping. We find that the inclusion of z<0.5 QSO clustering data improves the fit of a long-lived QSO model and suggest that the predictions of a PLE model for QSO BH masses agree reasonably with UV-bump and reverberation estimates (abridged).

We use the SDSS, along with the NVSS and the WENSS to define a sample of 746 radio-loud quasars and measure their 330MHz to 1.4GHz spectral indexes. Following previous authors we take the spectral index as an indicator of the orientation towards the quasars with more pole-on sources having flatter spectral indexes. We use this proxy for orientation to investigate the effect observing angle may have on optical spectra. Quasars with flatter spectral indexes tend to be brighter. However, we find no indication of reddening in steep-spectrum objects to indicate obscuration by a torus as a possible explanation. Nor do we find increased redddening in the flat-spectrum sources which could imply a contribution from jet-related emission. We reproduce a previously-described anti-correlation between the width of the MgII line and radio spectral index indicating a disk-like geometry for the MgII BLR. In contrast to previous authors we find no such correlation for the CIV line suggesting a more isotropic high-ionisation BLR. Both the [OII] and [OIII] narrow lines have more flux in steep spectrum sources while the [OII]/[OIII] flux ratio is lower in these sources. To describe both of these effects we propose a simple geometric model in which the NLR exists primarily on the surface of optically thick clouds facing the active nucleus and the NLR is stratified such that higher-ionisation lines are found preferentially closer to the nucleus. Quantitatively we find that orientation may effect the observed strength of narrow lines, as well as ratios between lines, by a factor of ~2. These findings have implications for the use of narrow emission lines to estimate bolometric luminosities, as well as comparisons between narrow line luminosity functions for type 1s and type 2s, and the potential of emission-line diagnostic diagrams as an accurate tool with which to distinguish classes of AGN.

We perform an extensive analysis of the Civ line in three large spectroscopic surveys of quasars. Differing approaches for fitting the Civ line can be found in the literature, and we compare the most common methods to highlight the relative systematics associated with each. We develop a line fitting procedure and apply it to the Civ line in spectra from the SDSS, 2QZ and 2SLAQ surveys. Our results are compared with a previous study of the Mgii line in the same sample. Civ tends to be broader than the Mgii line in spectra that have both, and the average ratio between the lines is consistent with a simplistic model for a photoionised, virialised and stratified broad-line region. There is a statistically significant correlation between the widths of the Civ and Mgii lines. However, the correlation is weak, and the scatter around a best fit is only marginally less than the full dynamic range of line widths. Motivated by previous work, we examine the dispersion in the distribution of Civ line widths. We find that the dispersion in Civ line widths is essentially independent of both redshift and luminosity. This is in stark contrast to the Mgii line, which shows a strong luminosity dependence. Finally we consider our results in terms of their implications for virial black hole mass estimation. The inconsistency between Mgii and Civ line widths in single spectra, combined with the differing behaviour of the Mgii and Civ line width distributions, indicates that there must be an inconsistency between Mgii and Civ virial mass estimators. Furthermore, the level of intrinsic dispersion in Mgii and Civ line widths contributes less dynamic range to virial mass estimates than the error associated with the estimates. The indication is that the line width term in these UV virial mass estimators may be essentially irrelevant with respect to the typical uncertainty on a mass estimate.

We present the final spectroscopic QSO catalogue from the 2dF-SDSS LRG and QSO (2SLAQ) Survey. This is a deep, 18<g<21.85 (extinction corrected), sample aimed at probing in detail the faint end of the broad line AGN luminosity distribution at z<2.6. The candidate QSOs were selected from SDSS photometry and observed spectroscopically with the 2dF spectrograph on the Anglo-Australian Telescope. This sample covers an area of 191.9 deg^2 and contains new spectra of 16326 objects, of which 8764 are QSOs, and 7623 are newly discovered (the remainder were previously identified by the 2QZ and SDSS surveys). The full QSO sample (including objects previously observed in the SDSS and 2QZ surveys) contains 12702 QSOs. The new 2SLAQ spectroscopic data set also contains 2343 Galactic stars, including 362 white dwarfs, and 2924 narrow emission line galaxies with a median redshift of z=0.22. We present detailed completeness estimates for the survey, based on modelling of QSO colours, including host galaxy contributions. This calculation shows that at g~21.85 QSO colours are significantly affected by the presence of a host galaxy up to redshift z~1 in the SDSS ugriz bands. In particular we see a significant reddening of the objects in g-i towards fainter g-band magnitudes. This reddening is consistent with the QSO host galaxies being dominated by a stellar population of age at least 2-3 Gyr. The full catalogue, including completeness estimates, is available on-line at http://www.2slaq.info/

We measure the width of the MgII $\lambda2799$ line in quasar spectra from the SDSS, 2QZ and 2SLAQ surveys and, by invoking an unnormalised virial mass estimator, relate the scatter in line width to the scatter in mass in the underlying black hole population. We find conclusive evidence for a trend such that there is less scatter in line width, and hence black hole mass, in more luminous objects. However, the most luminous objects in our sample show such a low degree of scatter in line width that, when combined with measures for the intrinsic scatter in the radius-luminosity relation for the broad-line region in active galaxies, an inconsistency arises in the virial technique for estimating black hole masses. This analysis implies that, at least for the most luminous quasars, either there is little-to-no intrinsic scatter in the radius-luminosity relation or the MgII broad emission line region is not totally dominated by virial velocities. Finally we exploit the measured scatter in line widths to constrain models for the velocity field of the broad-line region. We show that the lack of scatter in broad line-widths for luminous quasars is inconsistent with a pure planar/disk-like geometry for the broad-line region... (abridged)

We investigate the relation between the mass of super-massive black holes (Mbh) in QSOs and the mass of the dark matter halos hosting them (Mdh). We measure the widths of broad emission lines (Mgii lambda 2798, Civ lambda 1549) from QSO composite spectra as a function of redshift. These widths are then used to determine virial black hole mass estimates. We compare our virial black hole mass estimates to dark matter halo masses for QSO hosts derived by Croom et al. (2005) based on measurements of QSO clustering. This enables us to trace the Mbh-Mdh relation over the redshift range z=0.5 to 2.5. We calculate the mean zero-point of the Mbh-Mdh relation to be Mbh=10^(8.4+/-0.2)Msun for an Mdh=10^(12.5)Msun. These data are then compared with several models connecting Mbh and Mdh as well as recent hydrodynamical simulations of galaxy evolution. We note that the flux limited nature of QSO samples can cause a Malmquist-type bias in the measured zero-point of the Mbh-Mdh relation. The magnitude of this bias depends on the scatter in the Mbh-Mdh relation, and we reevaluate the zero-point assuming three published values for this scatter. (abridged)

We present a spectroscopic survey of almost 15,000 candidate intermediate-redshift Luminous Red Galaxies (LRGs) brighter than i=19.8, observed with 2dF on the Anglo-Australian Telescope. The targets were selected photometrically from the Sloan Digital Sky Survey (SDSS) and lie along two narrow equatorial strips covering 180 sq deg. Reliable redshifts were obtained for 92% of the targets and the selection is very efficient: over 90% have redshifts between 0.45 and 0.8. More than 80% of the ~11,000 red galaxies have pure absorption-line spectra consistent with a passively-evolving old stellar population. The redshift, photometric and spatial distributions of the LRGs are described. The 2SLAQ data will be released publicly from mid-2006, providing a powerful resource for observational cosmology and the study of galaxy evolution.