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.     

Radio imaging of the Subaru/XMM–Newton Deep Field – II. The 37 brightest radio sources

We study the 37 brightest radio sources in the Subaru/ XMM–Newton Deep Field. 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 'Fanaroff–Riley type I/type II (FR I/FR II)' luminosity divide. Using mid-infrared (mid-IR) ( Spitzer MIPS 24 μm) 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 FR I/FR II 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 FR I radio structure. For extended steep-spectrum radio sources, the fraction of objects with mid-IR excess drops dramatically below the FR I/FR II 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 per cent of the light, and the dust above and below the torus scatters ≳1 per cent of the light.     

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 Deep X-Ray Radio Blazar Survey (DXRBS) – II. New identifications

We have searched the archived, pointed ROSAT Position Sensitive Proportional Counter data for blazars by correlating the WGACAT X-ray data base with several publicly available radio catalogues, restricting our candidate list to serendipitous X-ray sources with a flat radio spectrum ( α _r≤0.70, where S _ν ∝ ν ^{− α} ). This makes up the Deep X-ray Radio Blazar Survey (DXRBS). Here we present new identifications and spectra for 106 sources, including 86 radio-loud quasars, 11 BL Lacertae objects, and nine narrow-line radio galaxies. Together with our previously published objects and already-known sources, our sample now contains 298 identified objects: 234 radio-loud quasars [181 flat-spectrum quasars: FSRQ ( α _r≤0.50) and 53 steep-spectrum quasars: SSRQ], 36 BL Lacs and 28 narrow-line radio galaxies. Redshift information is available for 96 per cent of these. Thus our selection technique is ∼90 per cent efficient at finding radio-loud quasars and BL Lacs. Reaching 5-GHz radio fluxes ∼50 mJy and 0.1–2.0 keV X-ray fluxes a few ×10⁻¹⁴ erg cm⁻² s⁻¹, DXRBS is the faintest and largest flat-spectrum radio sample with nearly complete (∼85 per cent) identification. We review the properties of the DXRBS blazar sample, including redshift distribution and coverage of the X-ray-radio–power plane for quasars and BL Lacs. Additionally, we touch upon the expanded multiwavelength view of blazars provided by DXRBS. By sampling for the first time the faint end of the radio and X-ray luminosity functions, this sample will allow us to investigate the blazar phenomenon and the validity of unified schemes down to relatively low powers.     

Ultraviolet Fe ii emission in z∼ 2 quasars

We present spectra of six luminous quasars at   z ∼ 2  , covering rest wavelengths 1600−3200 Å. The fluxes of the UV Fe  ii emission lines and Mg  ii λ2798 doublet, the line widths of Mg  ii and the 3000 Å luminosity were obtained from the spectra. These quantities were compared with those of low-redshift quasars at   z = 0.06–0.55  studied by Tsuzuki et al. In a plot of the Fe  ii (UV)/Mg  ii flux ratio as a function of the central black hole mass, Fe  ii (UV)/Mg  ii in our   z ∼ 2  quasars is systematically greater than in the low-redshift quasars. We confirmed that luminosity is not responsible for this excess. It is unclear whether this excess is caused by rich Fe abundance at   z ∼ 2  over low-redshift or by non-abundance effects such as high gas density, strong radiation field and high microturbulent velocity.     

Downsizing of supermassive black holes from the SDSS quasar survey

Starting from the ∼50 000 quasars of the Sloan Digital Sky Survey for which Mg  ii line width and 3000 Å monochromatic flux are available, we aim to study the dependence of the mass of active black holes on redshift. We focus on the observed distribution in the full width at half-maximum–nuclear luminosity plane, which can be reproduced at all redshifts assuming a limiting M _BH, a maximum Eddington ratio and a minimum luminosity (due to the survey flux limit). We study the z -dependence of the best-fitting parameters of assumed distributions at increasing redshift and find that the maximum mass of the quasar population evolves as  log ( M _BH(max)/M_⊙) ∼ 0.3 z + 9  , while the maximum Eddington ratio (∼0.45) is practically independent of cosmic time. These results are unaffected by the Malmquist bias.     

The active nuclei of z < 1.0 3CRR radio sources

We combine Chandra and XMM–Newton X-ray data from our previous papers with new X-ray observations and with Spitzer mid-infrared (mid-IR) data in order to study the nature of the nuclei of radio galaxies and radio-loud quasars with   z < 1.0  from the 3CRR sample. The significant increase in sample size over our previous work, the reduction of bias in the sample as a result of new observations and the availability of more mid-IR data allow us to show conclusively that almost all objects classed as low-excitation radio galaxies in optical spectroscopic studies lack a radiatively efficient active nucleus. We show that the distribution of absorbing columns in the narrow-line radio galaxies differs from the population of X-ray-selected radio-quiet type 2 quasars and from that in local Seyfert 2s. We comment on the current evidence for the nature of the soft X-ray component in radio-galaxy nuclear spectra, concluding that a jet origin for this component is very hard to evade. Finally, we discuss the recently discovered 'fundamental plane' of black hole activity, showing that care must be taken when placing radio-loud active galactic nucleus (AGN) on such diagnostic diagrams.     

A search for distant radio-loud quasars in the CLASS survey: three new radio-selected quasars at z > 4

We report on the search for distant radio-loud quasars in the Cosmic Lens All Sky Survey (CLASS) of flat spectrum radio sources with S _5GHz>30 mJy . Unresolved optical counterparts were selected from APM scans of POSS-I plates, with e <19.0 and red o − e >2.0 colours, in an effective area of ∼6400 deg². Four sources were found to be quasars with z >4 , of which one was previously known. This sample bridges the gap between the strong radio surveys with S _5GHz>200 mJy and the samples of radio-weak quasars that can be generated via radio observations of optically selected quasars. In addition, four new quasars at z >3 have been found. The selection criteria result in a success-rate of ∼1:7 for radio-loud quasars at z >4 , which is a significant improvement over previous studies. This search yields a surface density of 1 per 1600 deg², which is about a factor of ∼15 lower than that found in a similar search for radio-quiet quasars at z >4 . The study presented here is strongly biased against quasars beyond z >4.5 , since the e -passband of the POSS-I only samples the spectra shortward of 1200 Å at these redshifts.     

The central engines of radio-quiet quasars

Two rival hypotheses have been proposed for the origin of the compact radio flux observed in radio-quiet quasars (RQQs). It has been suggested that the radio emission in these objects, typically some two or three orders of magnitude less powerful than in radio-loud quasars (RLQs), either represents emission from a circumnuclear starburst or is produced by radio jets with bulk kinetic powers ∼ 10³ times lower than those of RLQs with similar luminosity ratios in other wavebands. We describe the results of high-resolution (∼pc-scale) radio-imaging observations of a sample of 12 RQQs using the Very Long Baseline Array (VLBA). We find strong evidence for jet-producing central engines in eight members of our sample.     

Optical spectroscopy of faint gigahertz peaked‐spectrum sources

We present spectroscopic observations of a sample of faint gigahertz peaked‐spectrum (GPS) radio sources drawn from the Westerbork Northern Sky Survey (WENSS). Redshifts have been determined for 19 (40 per cent) of the objects. The optical spectra of the GPS sources identified with low‐redshift galaxies show deep stellar absorption features. This confirms previous suggestions that their optical light is not significantly contaminated by active galactic nucleus-related emission, but is dominated by a population of old (>9 Gyr) and metal-rich (>0.2 [Fe/H]) stars, justifying the use of these (probably) young radio sources as probes of galaxy evolution. The optical spectra of GPS sources identified with quasars are indistinguishable from those of flat-spectrum quasars, and clearly different from the spectra of compact steep‐spectrum (CSS) quasars. The redshift distribution of the GPS quasars in our radio-faint sample is comparable to that of the bright samples presented in the literature, peaking at z ∼2–3. It is unlikely that a significant population of low-redshift GPS quasars is missed as a result of selection effects in our sample. We therefore claim that there is a genuine difference between the redshift distributions of GPS galaxies and quasars, which, because it is present in both the radio-faint and bright samples, cannot be caused by a redshift–luminosity degeneracy. It is therefore unlikely that the GPS quasars and galaxies are unified by orientation, unless the quasar opening angle is a strong function of redshift. We suggest that the GPS quasars and galaxies are unrelated populations and just happen to have identical observed radio spectral properties, and hypothesize that GPS quasars are a subclass of flat-spectrum quasars.     

First results from the Isaac Newton Telescope Wide Angle Survey: the quasar survey

We report the first results of an observational programme designed to determine the luminosity density of high-redshift quasars     quasars) using deep multicolour CCD data. We report the discovery and spectra of three     high-redshift     quasars, including one with     . At     , this is the fourth highest redshift quasar currently published. Using these preliminary results we derive an estimate of the         quasar space density in the redshift range     of     . When completed, the survey will provide a firm constraint on the contribution to the ionizing UV background in the redshift range     from quasars by determining the faint-end slope of the quasar luminosity function. The survey uses imaging data taken with the 2.5-m Isaac Newton Telescope as part of the Public Isaac Newton Group Wide Field Survey (WFS). This initial sample of objects is taken from two fields of effective area ∼12.5 deg² from the final ∼100 deg².     

The radio galaxy 3C 356 and clues to the trigger mechanisms for powerful radio sources

We present deep near-infrared images, taken with the 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 active galactic nuclei, 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.     

An 80-kpc Lyα halo around a high-redshift type-2 quasi-stellar object

We announce the discovery of an extended emission-line region associated with a high-redshift type-2 quasi-stellar object (QSO). The halo, which was discovered in our new wide-field narrow-band survey, resides at   z = 2.85  in the Spitzer First Look Survey region and is extended over ∼80 kpc. Deep very long baseline interferometry (VLBI) observations imply that approximately 50 per cent of the radio emission is extended on scales >200 pc. The inferred active galactic nuclei (AGN) luminosity is sufficient to ionize the extended halo, and the optical emission is consistent with being triggered coevally with the radio source. The Lyα halo is as luminous as those found around high-redshift radio galaxies; however, the active nucleus is several orders of magnitude less luminous at radio wavelengths than those Fanarof–Riley type II (FRIIs) more commonly associated with extended emission-line regions. AMS05 appears to be a high-redshift analogue to the radio-quiet quasar E1821+643 which is core dominated, but which also exhibits extended Fanarof–Riley type I (FRI)-like structure and contains an optically powerful AGN. We also find evidence for more quiescent kinematics in the Lyα emission line in the outer regions of the halo, reminiscent of the haloes around the more powerful FRIIs. The optical to mid-infrared spectral energy distribution is well described by a combination of an obscured QSO  ( L _bol∼ 3.4 ± 0.2 × 10¹³ L_⊙)  and a 1.4 Gyr old simple stellar population with mass  ∼3.9 ± 0.3 × 10¹¹ M_⊙  .     

Hubble Space Telescope imaging survey of sub-mJy star-forming galaxies – I. Morphologies at z ∼ 0.2

We present the first results of our Hubble Space Telescope HST WFPC2 F814W snapshot imaging survey, targeting virtually all sub-mJy decimetric radio-selected star-forming galaxies. The radio selection at ∼1 GHz is free from extinction effects and the radio luminosities are largely unaffected by AGN contamination, making these galaxies ideal tracers of the cosmic star formation history. A subsample of four targets is presented here, selected at 1.4 GHz from the spectroscopically homogenous and complete samples of Benn et al. and Hopkins et al. The redshifts are confined to a narrow range around z ∼0.2, to avoid differential evolution, with a radio luminosity close to L ∗ where the galaxies dominate the comoving volume-averaged star formation rate. We find clearly disturbed morphologies resembling those of ultraluminous infrared galaxies, indicating that galaxy interactions may be the dominant mechanism for triggering star formation at these epochs. The morphologies are also clearly different from those of coeval quasars and radio galaxies, as found in star-forming galaxies selected at other wavelengths. This may prove challenging for models that propose direct causal links between AGN evolution and the cosmic star formation history at these epochs. The asymmetries are typically much larger than seen in the Canada–France Redshift Survey at similar redshifts, optical luminosities and H α -derived star formation rates, indicating the possible existence of an obscuration-related morphological bias in such samples.     

A comparative HST imaging study of the host galaxies of radio-quiet quasars, radio-loud quasars and radio galaxies – I

We present the first results from a major HST WFPC2 imaging study aimed at providing the first statistically meaningful comparison of the morphologies, luminosities, scalelengths and colours of the host galaxies of radio-quiet quasars, radio-loud quasars and radio galaxies. We describe the design of this study and present the images that have been obtained for the first half of our 33-source sample. We find that the hosts of all three classes of luminous AGN are massive elliptical galaxies, with scalelengths ≃10 kpc, and R − K colours consistent with mature stellar populations. Most importantly, this is first unambiguous evidence that, just like radio-loud quasars, essentially all radio-quiet quasars brighter than M _R =−24 reside in massive ellipticals. This result removes the possibility that radio 'loudness' is directly linked to host galaxy morphology, but is however in excellent accord with the black hole/spheroid mass correlation recently highlighted by Magorrian et al. We apply the relations given by Magorrian et al. to infer the expected Eddington luminosity of the putative black hole at the centre of each of the spheroidal host galaxies we have uncovered. Comparison with the actual nuclear R -band luminosities suggests that the black holes in most of these galaxies are radiating at a few per cent of the Eddington luminosity; the brightest host galaxies in our low- z sample are capable of hosting quasars with M _R ≃− 28, comparable to the most luminous quasars at z ≃3. Finally, we discuss our host-derived black hole masses in the context of the radio luminosity:black hole mass correlation recently uncovered for nearby galaxies by Franceschini et al., and consider the resulting implications for the physical origin of radio loudness.     

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.     

Relativistic beaming effects in the spectra of cores and hotspots in radio galaxies and quasars

We present the results of multiwavelength observations of cores and hotspots, at L , C , X and U bands with the Very Large Array, of a matched sample of radio galaxies and quasars selected from the Molonglo Reference Catalogue . We use these observations to determine the spectra of cores and hotspots, and test the unified scheme for radio galaxies and quasars. Radio cores have been detected at all wavelengths in all of the quasars in our sample, whereas only ∼50 per cent of the galaxies have cores detected in at least one of the wavelengths . The degree of core prominence in this sample is consistent with the unified scheme for radio galaxies and quasars. A comparison of the distributions of the two-point spectral index of the cores in our sample of lobe-dominated quasars, with the distributions in a matched sample of core-dominated quasars, shows that the distributions for these two classes are significantly different, and this is consistent with the expectations of the unified scheme. The difference in the spectral indices of the two hotspots on opposite sides is also significantly larger for quasars than for radio galaxies, as is expected in the unified scheme. We also investigate the relationship between the spectral index of the hotspots and the redshift or luminosity for our sample of sources.     

Chandra measurements of the X-ray core and cluster of 3C 220.1

We report results of an 18-ks exposure with the ACIS instrument on Chandra of the powerful z =0.62 radio galaxy 3C 220.1. The X-ray emission separates into cluster gas of emission-weighted kT ∼5 keV , 0.7–12 keV luminosity (to a radius of 45 arcsec) of 5.6×10⁴⁴ erg s⁻¹ and unresolved emission (coincident with the radio core). While the extended X-ray emission is clearly thermal in nature, a straightforward cooling-flow model, even in conjunction with a point-source component, is a poor fit to the radial profile of the X-ray emission. This is despite the fact that the measured properties of the gas suggest a massive cooling flow of ∼130 M_⊙ yr⁻¹, and the data show weak evidence for a temperature gradient. The central unresolved X-ray emission has a power-law spectral energy index α ∼0.7 and 0.7–12 keV luminosity of 10⁴⁵ erg s⁻¹, and any intrinsic absorption is relatively small. The two-point spectrum of the core emission between radio and X-ray energies has α _rx=0.75 . Since this is a flatter spectrum than seen in other sources where the X-ray emission is presumed to be radio-related, regions close to the active galactic nucleus (AGN) in this source may dominate the central X-ray output, as is believed to be the case for lobe-dominated quasars. Simple unification models would be challenged if this were found to be the case for a large fraction of high-power radio galaxies.     

A simple model for the evolution of supermassive black holes and the quasar population

An empirically motivated model is presented for accretion-dominated growth of supermassive black holes (SMBH) in galaxies, and the implications are studied for the evolution of the quasar population in the Universe. We investigate the core aspects of the quasar population, including space density evolution, evolution of the characteristic luminosity, plausible minimum masses of quasars, the mass function of SMBH and their formation epoch distribution. Our model suggests that the characteristic luminosity in the quasar luminosity function arises primarily as a consequence of a characteristic mass scale above which there is a systematic separation between the black hole and the halo merging rates. At lower mass scales, black hole merging closely tracks the merging of dark haloes. When combined with a declining efficiency of black hole formation with redshift, the model can reproduce the quasar luminosity function over a wide range of redshifts. The observed space density evolution of quasars is well described by formation rates of SMBH above  ∼10⁸  M_⊙  . The inferred mass density of SMBH agrees with that found independently from estimates of the SMBH mass function derived empirically from the quasar luminosity function.     

High-redshift obscured quasars: radio emission at sub-kiloparsec scales

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 μJy per 17 × 14 mas² beam  was achieved, and in seven out of 11 sources unresolved radio emission was securely detected. The detected radio emission of each source accounts for ∼30–100 per cent 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 Fanaroff–Riley type I (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.