Cosmological evolution of the Fanaroff–Riley type II source population

By combining a model for the evolution of the radio luminosity of an individual source with the radio luminosity function, we perform a multidimensional Monte Carlo simulation to investigate the cosmological evolution of the Fanaroff–Riley type II (FR II) radio galaxy population by generating large artificial samples. The properties of FR II sources are required to evolve with redshift in the artificial samples to fit the observations. Either the maximum jet age or the maximum density of the jet environment or both evolve with redshift. We also study the distribution of FR II source properties as a function of redshift. From currently available data we cannot constrain the shape of the distribution of environment density or age, but jet power is found to follow a power-law distribution with an exponent of approximately −2. This power-law slope does not change with redshift out to   z = 0.6  . We also find the distribution of the pressure in the lobes of FR II sources to evolve with redshift up to   z ∼ 1.2  .     

A submillimetre survey of the star formation history of radio galaxies

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-μm 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 at a rate ∝(1+ z )³ 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.     

Giant radio sources

We present multifrequency Very Large Array (VLA) observations of two giant quasars, 0437−244 and 1025−229, from the Molonglo Complete Sample. These sources have well-defined FR II radio structure, possible one-sided jets, no significant depolarization between 1365 and 4935 MHz and low rotation measure (|RM|<20 rad m⁻²). The giant sources are defined to be those with overall projected size 1 Mpc. We have compiled a sample of about 50 known giant radio sources from the literature, and have compared some of their properties with a complete sample of 3CR radio sources of smaller sizes to investigate the evolution of giant sources, and test their consistency with the unified scheme for radio galaxies and quasars. We find an inverse correlation between the degree of core prominence and total radio luminosity, and show that the giant radio sources have similar core strengths to smaller sources of similar total luminosity. Hence their large sizes are unlikely to be caused by stronger nuclear activity. The degree of collinearity of the giant sources is also similar to that of the sample of smaller sources. The luminosity–size diagram shows that the giant sources are less luminous than our sample of smaller sized 3CR sources, consistent with evolutionary scenarios in which the giants have evolved from the smaller sources, losing energy as they expand to these large dimensions. For the smaller sources, radiative losses resulting from synchrotron radiation are more significant while for the giant sources the equipartition magnetic fields are smaller and inverse Compton loss owing to microwave background radiation is the dominant process. The radio properties of the giant radio galaxies and quasars are consistent with the unified scheme.     

Optical spectroscopy of two overlapping, flux-density-limited samples of radio sources in the North Ecliptic Cap, selected at 38 and 151 MHz

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 contain many of the lowest-luminosity z >2 radio galaxies so far discovered. 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 ₁₅₁≥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 ₃₈≥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, including a number of cases of probable gravitational lensing. 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 by comparing the results from traditional flux-density selection with our new method. 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.     

A sample of 6C radio sources designed to find objects at redshift z>4– III. Imaging and the radio galaxy K–z relation

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 (Ω_M=0.3 and Ω_Λ=0.7) , the most luminous radio sources appear to be associated with galaxies with a luminosity distribution with a high mean (≈5  L *), and a low dispersion ( σ ∼0.5 mag) which formed their stars at epochs corresponding to z ≳2.5 . This result is in line with recent submillimetre studies of high-redshift radio galaxies and the inferred ages of extremely red objects from faint radio samples.     

On the redshift cut-off for steep-spectrum radio sources

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. Though 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, incorporating 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     and     . We find that 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 probability ratios of     and     for spatially flat cosmologies with     and     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 and Jarvis & Rawlings, in line with the expectations of unified schemes.     

The Phoenix Deep Survey: 1.4-GHz source counts

We present the first results from the Phoenix Deep Survey, a multiwavelength survey of a 2° diameter region. Observations in the radio continuum at 1.4 GHz carried out with the Australia Telescope Compact Array are described. The catalogue of over 1000 radio sources compiled from these observations is analysed, and the source counts are presented. We model the observational source counts using a two-population model and published luminosity functions for these populations. Upon invoking luminosity and density evolution, we find that a luminosity evolution model best fits the radio observations, consistent with earlier work. The redshift distribution of the two galaxy populations investigated is also modelled and discussed.     

The relationship between star formation rate and radio synchrotron luminosity at 0 < z < 2

We probe the relationship between star formation rate (SFR) and radio synchrotron luminosity in galaxies at  0 < z < 2  within the northern Spitzer Wide-area Infrared Extragalactic survey (SWIRE) fields, in order to investigate some of the assumptions that go into calculating the star formation history of the Universe from deep radio observations. We present new 610-MHz Giant Metrewave Radio Telescope (GMRT) observations of the European Large-Area ISO Survey-North 2 (ELAIS-N2) field, and using this data, along with previous GMRT surveys carried out in the ELAIS-N1 (North 1) and Lockman Hole regions, we construct a sample of galaxies which have redshift and SFR information available from the SWIRE survey. We test whether the local relationship between SFR and radio luminosity is applicable to   z = 2  galaxies, and look for evolution in this relationship with both redshift and SFR in order to examine whether the physical processes which lead to synchrotron radiation have remained the same since the peak of star formation in the Universe. We find that the local calibration between radio luminosity and star formation can be successfully applied to radio-selected high-redshift, high-SFR galaxies, although we identify a small number of sources where this may not be the case; these sources show evidence for inaccurate estimations of their SFR, but there may also be some contribution from physical effects such as the recent onset of starburst activity, or suppression of the radio luminosity within these galaxies.     

FR II radio galaxies with z < 0.3 — I. Properties of jets, cores and hotspots

In previous papers we have discussed high-resolution observations of a large sample of powerful radio galaxies with z  < 0.3. Jets are detected in up to 80 per cent of the sample, and radio cores in nearly all the objects; in addition, we are able to resolve the hotspots in most sources. In this paper we present measurements of the radio properties of these components.   The prominences of the jets detected do not appear to be a function of radio luminosity, providing the clearest evidence yet that the reported low detection rate of jets in radio galaxies has been an artefact of low-sensitivity observations. We find a positive correlation between the total source length and core prominence in the narrow-line radio galaxies. We have found evidence for a relationship between hotspot size and total source size, but few other significant relationships between hotspot properties and those of the jets or lobes. We compare our measurements with those of Bridle et al., based on observations of a sample of quasars, and argue that the results are consistent with a modification of the unified model in which the broad-line radio galaxies are the low-luminosity counterparts of quasars, although the situation is complicated by contamination with low-excitation radio galaxies which appear to have radio properties different from the high-excitation objects. We discuss the classes of empirical model that can be fitted to the data set.     

HST and UKIRT imaging observations of z∼ 1 6C radio galaxies – II. Galaxy morphologies and the alignment effect

Powerful radio galaxies often display enhanced optical/ultraviolet emission regions, elongated and aligned with the radio jet axis. The aim of this series of papers is to investigate separately the effects of radio power and redshift on the alignment effect, together with other radio galaxy properties. In this second paper, we present a deeper analysis of the morphological properties of these systems, including both the host galaxies and their surrounding aligned emission.
The host galaxies of our 6C subsample are well described as de Vaucouleurs ellipticals, with typical scale sizes of  ∼10 kpc  . This is comparable to the host galaxies of low- z radio sources of similar powers, and also the more powerful 3CR sources at the same redshift. The contribution of nuclear point source emission is also comparable, regardless of radio power.
The 6C alignment effect is remarkably similar to that seen around more powerful 3CR sources at the same redshift in terms of extent and degree of alignment with the radio source axis, although it is generally less luminous. The bright, knotty features observed in the case of the z ∼ 1 3CR sources are far less frequent in our 6C subsample; neither do we observe such strong evidence for evolution in the strength of the alignment effect with radio source size/age. However, we do find a very strong link between the most extreme alignment effects and emission-line region properties indicative of shocks, regardless of source size/age or power. In general, the 6C alignment effect is still considerably stronger than that seen around lower redshift galaxies of similar radio powers. Cosmic epoch is clearly just as important a factor as radio power: although aligned emission is observed on smaller scales at lower redshifts, the processes which produce the most extreme features simply no longer occur, suggesting considerable evolution in the properties of the extended haloes surrounding the radio source.     

The 6C** sample of steep-spectrum radio sources – II. Redshift distribution and the space density of high-redshift radio galaxies

We use the 6C** sample to investigate the comoving 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 comoving space density of steep-spectrum radio sources beyond z ≳ 2.5, and rule out a steep decline.     

The SCUBA Half-Degree Extragalactic Survey (SHADES) – VIII. The nature of faint submillimetre galaxies in SHADES, SWIRE and SXDF surveys

We present the optical-to-submillimetre spectral energy distributions (SEDs) for 33 radio and mid-infrared (mid-IR) identified submillimetre galaxies discovered via the SHADES 850-μm SCUBA imaging in the Subaru- XMM Deep Field (SXDF). Optical data for the sources come from the 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-μm fluxes. We then fit 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 IR template fits are also used to resolve ambiguous identifications and cases of redshift aliasing. The redshift distribution obtained broadly matches previous results for submillimetre sources and on the SHADES SXDF field. Our template fitting finds that active galactic nuclei, while present in about 10 per cent of our 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 fitted 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 submillimetre 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.     

Radio properties of the Shapley Concentration — II. J1324–3138: a remnant of a radio galaxy in the Abell cluster A3556?

In this paper we present a detailed study of the radio galaxy J1324–3138, located at a projected distance of 2 arcmin from the centre of the Abell cluster of galaxies A3556, belonging to the core of the Shapley Concentration, at an average redshift of z  = 0.05. We have observed J1324–3138 over a wide range of frequencies: at 327 MHz (VLA), 843 MHz (MOST), and at 1376, 2382, 4790 and 8640 MHz (ATCA).   Our analysis suggests that J1324–3138 is a remnant of a tailed radio galaxy, in which the nuclear engine has switched off and the radio source is now at a late stage of its evolution, confined by the intracluster gas. The radio galaxy is not in pressure equilibrium with the external medium, as is often found for extended radio sources in clusters of galaxies. We favour the hypothesis that the lack of observed polarized radio emission in the source is a result of Faraday rotation by a foreground screen, i.e. the source is seen through a dense cluster gas, characterized by a random magnetic field.   An implication of the head–tail nature of the source is that J1324–3138 is moving away from the core of A3556 and that, possibly, a major merging event between the core of A3556 and the subgroup hosting J1324–3138 has already taken place.     

Implications for unified schemes from submillimetre and far-infrared follow-up of radio-selected samples

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-μm) 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 ₈₅₀ of low-frequency-selected radio sources is correlated with the bolometric luminosity L _bol of their quasar nuclei via an approximate scaling relation   L ₈₅₀∝ L ^0.7±0.2_bol  . Secondly, 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. Thirdly, 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 ₈₅₀∝ L ^1/2_bol  which is within the errors of that determined here for radio-selected quasars, and consistent with that inferred for radio-quiet quasars.     

The quasar fraction in low-frequency-selected complete samples and implications for unified schemes

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 [O  ii ] emission-line luminosity of log₁₀( L _{[O  ii ]}/W)≳35 [or radio luminosity log₁₀( L ₁₅₁/W Hz⁻¹ sr⁻¹)≳ 26.5], 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 θ _trans≈53°. 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 θ _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.     

Radio sources in the 2dF Galaxy Redshift Survey – II. Local radio luminosity functions for AGN and star-forming galaxies at 1.4 GHz

We have cross-matched the 1.4-GHz NRAO VLA Sky Survey (NVSS) with the first 210 fields observed in the 2dF Galaxy Redshift Survey (2dFGRS), covering an effective area of 325 deg² (about 20 per cent of the final 2dFGRS area). This yields a set of optical spectra of 912 candidate NVSS counterparts, of which we identify 757 as genuine radio identifications – the largest and most homogeneous set of radio source spectra ever obtained. The 2dFGRS radio sources span the redshift range     to 0.438, and are a mixture of active galaxies (60 per cent) and star-forming galaxies (40 per cent). About 25 per cent of the 2dFGRS radio sources are spatially resolved by NVSS, and the sample includes three giant radio galaxies with projected linear size greater than 1 Mpc. The high quality of the 2dF spectra means we can usually distinguish unambiguously between AGN and star-forming galaxies. We make a new determination of the local radio luminosity function at 1.4 GHz for both active and star-forming galaxies, and derive a local star formation density of         .     

Clustering of galaxies around radio quasars at 0.5≤z≤0.8

We have observed the galaxy environments around a sample of 21 radio-loud, steep-spectrum quasars at 0.5≤ z ≤0.82, spanning several orders of magnitude in radio luminosity. The observations also include background control fields used to obtain the excess number of galaxies in each quasar field. The galaxy excess was quantified using the spatial galaxy–quasar correlation amplitude, B _gq, and an Abell-type measurement, N _0.5. A few quasars are found in relatively rich clusters, but on average, they seem to prefer galaxy groups or clusters of approximately Abell class 0. We have combined our sample with literature samples extending down to z ≈0.2 and covering the same range in radio luminosity. By using the Spearman statistic to disentangle redshift and luminosity dependences, we detect a weak, but significant, positive correlation between the richness of the quasar environment and the radio luminosity of the quasar. However, we do not find any epoch dependence in B _gq, as has previously been reported for radio quasars and galaxies. We discuss the radio luminosity–cluster richness link and possible explanations for the weak correlation that is seen.     

X-ray variability in a deep, flux-limited sample of QSOs

We present an analysis of X-ray variability in a flux-limited sample of quasi-stellar objects (QSOs). Selected from our deep ROSAT survey, these QSOs span a wide range in redshift (0.1< z <3.2) and are typically very faint, so we have developed a method to constrain the amplitude of variability in ensembles of low signal-to-noise ratio light curves. We find evidence for trends in this variability amplitude with both redshift and luminosity. The mean variability amplitude declines sharply with luminosity, as seen in local active galactic nuclei (AGN), but with some suggestion of an upturn for the most powerful sources. We find tentative evidence that this is caused by redshift evolution, since the high-redshift QSOs ( z >0.5) do not show the anticorrelation with luminosity seen in local AGN. We speculate on the implications of these results for physical models of AGN and their evolution. Finally, we find evidence for X-ray variability in an object classified as a narrow-emission-line galaxy, suggesting the presence of an AGN.     

1–5 μm imaging of 3CRR galaxies: the K–z relation and the geometry of the torus

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

The far-infrared–radio relationship at high and low redshift

We use the results of the SCUBA Local Universe Galaxy Survey, a submillimetre (submm) survey of galaxies in the nearby Universe, to investigate the relationship between the far-infrared (FIR)–submm and radio emission of galaxies at both low and high redshift. At low redshift we show that the correlation between radio and FIR emission is much stronger than the correlation between radio and submm emission, which is evidence that massive stars are the source of both the FIR and radio emission. At high redshift we show that the submm sources detected by SCUBA are brighter sources of radio emission than are predicted from the properties of galaxies in the local Universe. We discuss possible reasons for the cosmic evolution of the relationship between radio and FIR emission.