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.     

The cluster environments of the z ∼ 1 3CR radio galaxies

An analysis of the environments around a sample of 28 3CR radio galaxies with redshifts 0.6< z <1.8 is presented, based primarily upon K -band images down to K ∼20 taken using the UK Infrared Telescope (UKIRT). A net overdensity of K -band galaxies is found in the fields of the radio galaxies, with the mean excess counts being comparable to that expected for clusters of Abell Class 0 richness. A sharp peak is found in the angular cross-correlation amplitude centred on the radio galaxies that, for reasonable assumptions about the luminosity function of the galaxies, corresponds to a spatial cross-correlation amplitude between those determined for low-redshift Abell Class 0 and 1 clusters.
These data are complemented by J -band images also from UKIRT, and by optical images from the Hubble Space Telescope . The fields of the lower redshift ( z ≲0.9) radio galaxies in the sample generally show well-defined near-infrared colour–magnitude relations with little scatter, indicating a significant number of galaxies at the redshift of the radio galaxy; the relations involving colours at shorter wavelengths than the 4000 Å break show considerably greater scatter, suggesting that many of the cluster galaxies have low levels of recent or on-going star formation. At higher redshifts the colour–magnitude sequences are less prominent owing to the increased field galaxy contribution at faint magnitudes, but there is a statistical excess of galaxies with the very red infrared colours ( J − K ≳1.75) expected of old cluster galaxies at these redshifts.
Although these results are appropriate for the mean of all of the radio galaxy fields, there exist large field-to-field variations in the richness of the environments. Many, but certainly not all, powerful z ∼1 radio galaxies lie in (proto)cluster environments.     

Near-infrared luminosity function and colours of dwarf galaxies in the Coma cluster

We present K -band observations of the low-luminosity galaxies in the Coma cluster, which are responsible for the steep upturn in the optical luminosity function at M _R∼−16, discovered recently. The main results of this study are as follows.
(i) The optical–near-infrared colours of these galaxies imply that they are dwarf spheroidal galaxies. The median B − K colour for galaxies with −19.3< M_K <−16.3 is 3.6 mag.
(ii) The K -band luminosity function in the Coma cluster is not well constrained, because of the uncertainties due to the field-to-field variance of the background. However, within the estimated large errors, this is consistent with the R -band luminosity function, shifted by ∼3 mag.
(iii) Many of the cluster dwarfs lie in a region of the B − K versus B − R colour–colour diagram where background galaxies are rare ( B − K <5; 1.2< B − R <1.6). Local dwarf spheroidal galaxies lie in this region too. This suggests that a better measurement of the K -band cluster luminosity can be made if the field-to-field variance of the background can be measured as a function of colour, even if it is large.
(iv) If we assume that none of the galaxies in the region of the B − K versus B − R plane given in (iii) in our cluster fields are background, and that all the cluster galaxies with 15.5< K <18.5 lie in this region of the plane, then we measure α=−1.41^+0.34_−0.37 for −19.3< M_K −16.3, where α is the logarithmic slope of the luminosity function. The uncertainties in this number come from counting statistics.     

The ages of quasar host galaxies

We present the results of fitting deep off-nuclear optical spectra of radio-quiet quasars, radio-loud quasars and radio galaxies at z ≃0.2 with evolutionary synthesis models of galaxy evolution. Our aim was to determine the age of the dynamically dominant stellar populations in the host galaxies of these three classes of powerful active galactic nuclei (AGN). Some of our spectra display residual nuclear contamination at the shortest wavelengths, but the detailed quality of the fits longward of the 4000-Å break provides unequivocal proof, if further proof were needed, that quasars lie in massive galaxies with (at least at z ≃0.2) evolved stellar populations. By fitting a two-component model we have separated the very blue (starburst and/or AGN contamination) from the redder underlying spectral energy distribution, and find that the hosts of all three classes of AGN are dominated by old stars of age 8–14 Gyr. If the blue component is attributed to young stars, we find that, at most, 1 per cent of the visible baryonic mass of these galaxies is involved in star formation activity at the epoch of observation, at least over the region sampled by our spectroscopic observations. These results strongly support the conclusion reached by McLure et al. that the host galaxies of luminous quasars are massive ellipticals which have formed by the epoch of peak quasar activity at z ≃2.5.     

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.     

The LBDS Hercules sample of mJy radio sources at 1.4 GHz – I. Multicolour photometry

The results are presented of an extensive programme of optical and infrared imaging of radio sources in a complete subsample of the Leiden–Berkeley Deep Survey. The LBDS Hercules sample consists of 72 sources observed at 1.4 GHz, with flux densities S _1.41.0 mJy, in a 1.2 deg² region of Hercules. This sample is almost completely identified in the g , r , i and K bands, with some additional data available at J and H . The magnitude distributions peak at r ≃22 mag, K ≃16 mag and extend down to r ≃26 mag, K ≃21 mag. The K -band magnitude distributions for the radio galaxies and quasars are compared with those of other radio surveys. At S _1.4 GHz≲1 Jy, the K -band distribution does not change significantly with radio flux density. The sources span a broad range of colours, with several being extremely red ( r − K ≳6). Though small, this is the most optically complete sample of mJy radio sources available at 1.4 GHz, and is ideally suited for studying the evolution of the radio luminosity function out to high redshifts.     

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.     

Optical and near-infrared imaging of faint Gigahertz Peaked Spectrum sources

A sample of 47 faint Gigahertz Peaked Spectrum (GPS) radio sources selected from the Westerbork Northern Sky Survey (WENSS) has been imaged in the optical and near-infrared, resulting in an identification fraction of 87 per cent. The R  −  I R  −  K colours of the faint optical counterparts are as expected for passively evolving elliptical galaxies, assuming that they follow the R -band Hubble diagram as determined for radio-bright GPS galaxies. We find evidence that the radio spectral properties of the GPS quasars are different from those of GPS galaxies. The observed distribution of radio spectral peak frequencies for GPS sources optically identified with bright stellar objects (presumably quasars) is shifted compared with GPS sources identified with faint or extended optical objects (presumably galaxies), in the sense that a GPS quasar is likely to have a higher peak frequency than a GPS galaxy. This means that the true peak frequency distribution is different for the GPS galaxies and quasars, because the sample selection effects are independent of optical identification. The correlation between peak frequency and redshift that has been suggested for bright sources has not been found in this sample; no correlation exists between R magnitude (and therefore redshift) and peak frequency for the GPS galaxies. We therefore believe that the claimed correlation is actually caused by the dependence of the peak frequency on optical host, because the GPS galaxies are generally at lower redshifts than the quasars. The difference in the peak frequency distributions of the GPS galaxies and quasars is further evidence against the hypothesis that they form a single class of object.     

Upper limits on K-band polarization in three high-redshift radio galaxies: 53W091, 3C 441 and MRC 0156−252

We present the results of K -band imaging polarimetry of three radio galaxies, including the very red and apparently old z =1.55 galaxy 53W091. We find weak evidence for polarization in components of 3C 441 and in the south-eastern companion of 53W091, but no evidence of significant polarization in 53W091 itself. We also find strong evidence that MRC 0156−252 is unpolarized. We present upper limits for the K -band polarization of all three sources. For 53W091, the lack of significant K -band polarization provides further confidence that its red R − K colour can be attributed to a mature stellar population, consistent with the detailed analyses of its ultraviolet spectral-energy distribution which indicate a minimum age of 2–3.5 Gyr.     

The K-band Hubble diagram of submillimetre galaxies and hyperluminous galaxies

We present the K -band Hubble diagrams ( K – z relations) of submillimetre-selected galaxies and hyperluminous galaxies (HLIRGs). We report the discovery of a remarkably tight K – z relation of HLIRGs, indistinguishable from that of the most luminous radio galaxies. Like radio galaxies, the HLIRG K – z relation at   z ≲ 3  is consistent with a passively evolving ∼3 L _* instantaneous starburst starting from a redshift of   z ∼ 10  . In contrast, many submillimetre-selected galaxies are ≳2 mag fainter, and the population has a much larger dispersion. We argue that dust obscuration and/or a larger mass range may be responsible for this scatter. The galaxies so far proved to be hyperluminous may have been biased towards higher AGN bolometric contributions than submillimetre-selected galaxies due to the 60-μm selection of some, so the location on the K – z relation may be related to the presence of the most massive active galactic nucleus. Alternatively, a particular host galaxy mass range may be responsible for both extreme star formation and the most massive active nuclei.     

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.     

Red companions to a z=2.15 radio-loud quasar

We have conducted observations of the environment around the z =2.15 radio-loud quasar 1550–269 in search of distant galaxies associated either with it or the z =2.09 C  iv absorber along its line of sight. Such objects will be distinguished by their red colours, R − K >4.5. We find five such objects in a 1.5 arcmin² field around the quasar, with typical K ' magnitudes of ∼20.4 and no detected R -band emission. We also find a sixth object with K =19.6±0.3, and undetected at R , just two arcsec from the quasar. The nature of all these objects is currently unclear, and will remain so until we have determined their redshifts. We suggest that it is likely that they are associated with either the quasar or the C  iv absorber, in which case their properties might be similar to those of the z =2.38 red Ly α emitting galaxies discovered by Francis et al. The small separation between the quasar and the brightest of our objects suggests that it may be the galaxy responsible for the C  iv metal line absorption system. The closeness to the quasar and the red colour might have precluded similar objects from being uncovered in previous searches for emission from C  iv and damped absorbers.     

Discovery of a single faint AGN in a large sample of z > 5 Lyman break galaxies

As part of a large spectroscopic survey of   z > 5  Lyman break galaxies (LBGs), we have identified a single source which is clearly hosting an active galactic nucleus (AGN). Out of a sample of more than 50 spectroscopically confirmed R -band dropout galaxies at   z ∼ 5  and above, only J104048.6−115550.2 at   z = 5.44  shows evidence for a high ionization potential emission line indicating the presence of a hard ionizing continuum from an AGN. Like most objects in our sample the rest-frame-UV spectrum shows the UV continuum breaking across a Lyα line. Uniquely within this sample of LBGs, emission from N  v is also detected, a clear signature of AGN photoionization. The object is spatially resolved in Hubble Space Telescope ( HST ) imaging. This, and the comparatively high Lyα/N  v flux ratio indicates that the majority of the Lyα (and the UV continuum longward of it) originates from stellar photoionization, a product of the ongoing starburst in the LBG. Even without the AGN emission, this object would have been photometrically selected and spectroscopically confirmed as a Lyman break in our survey. The measured optical flux  ( I _AB= 26.1)  is therefore an upper limit to that from the AGN and is of order 100 times fainter than the majority of known quasars at these redshifts. The detection of a single object in our survey volume is consistent with the best current models of high redshift AGN luminosity function, providing a substantial fraction of such AGN is found within luminous starbursting galaxies. We discuss the cosmological implications of this discovery.     

Optical off-nuclear spectra of quasar hosts and radio galaxies

We present optical (∼3200 to ∼9000 Å) off-nuclear spectra of 26 powerful active galaxies in the redshift range 0.1≤ z ≤0.3, obtained with the Mayall and William Herschel 4-m class telescopes. The sample consists of radio-quiet quasars, radio-loud quasars (all with −23≥ M _V ≥−26) and radio galaxies of Fanaroff–Riley Type II (with extended radio luminosities and spectral indices comparable to those of the radio-loud quasars). The spectra were all taken approximately 5 arcsec off-nucleus, with offsets carefully selected so as to maximize the amount of galaxy light falling into the slit, whilst simultaneously minimizing the amount of scattered nuclear light. The majority of the resulting spectra appear to be dominated by the integrated stellar continuum of the underlying galaxies rather than by light from the non-stellar processes occurring in the active nuclei, and in many cases a 4000-Å break feature can be identified. The individual spectra are described in detail, and the importance of the various spectral components is discussed. Stellar population synthesis modelling of the spectra will follow in a subsequent paper.     

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.     

Optical and infrared diagnostics of SDSS galaxies in the SWIRE survey

We present the rest-frame optical and infrared colours of a complete sample of  1114 z < 0.3  galaxies from the Spitzer Wide-Area Infrared Extragalactic (SWIRE) Legacy Survey and the Sloan Digital Sky Survey (SDSS). We discuss the optical and infrared colours of our sample and analyse in detail the contribution of dusty star-forming galaxies and active galactic nuclei (AGN) to optically selected red sequence galaxies.
We propose that the optical  ( g − r )  colour and infrared  log( L ₂₄/ L _3.6)  colour of galaxies in our sample are determined primarily by a bulge-to-disc ratio. The  ( g − r )  colour is found to be sensitive to the bulge-to-disc ratio for disc-dominated galaxies, whereas the  log( L ₂₄/ L _3.6)  colour is more sensitive for bulge-dominated systems.
We identify ∼18 per cent (195 sources) of our sample as having red optical colours and infrared excess. Typically, the infrared luminosities of these galaxies are found to be at the high end of star-forming galaxies with blue optical colours. Using emission-line diagnostic diagrams, 78 are found to have an AGN contribution and 117 are identified as star-forming systems. The red  ( g − r )  colour of the star-forming galaxies could be explained by extinction. However, their high optical luminosities cannot. We conclude that they have a significant bulge component.
The number densities of optically red star-forming galaxies are found to correspond to ∼13 per cent of the total number density of our sample. In addition, these systems contribute ∼13 per cent of the total optical luminosity density, and 28 per cent of the total infrared luminosity density of our SWIRE/SDSS sample. These objects may reduce the need for 'dry mergers'.     

HST Planetary Camera images of quasar host galaxies

We present Hubble Space Telescope ( HST ) images of seven low-redshift quasars (six taken with the Planetary Camera, one with the Wide Field Camera). These complete the sample of 14 quasars observed by the Faint Object Camera Investigation Definition Team (FOC IDT). Following subtraction of the quasar nuclear light, host galaxies can be seen in all seven cases. A combination of the optical morphology and luminosity profiles of the residual host galaxies and the results of 2D cross-correlation model fitting implies that five of the objects have elliptical host galaxies and two have disc host galaxies. The luminosities vary from slightly fainter than L ^* to about 1.3 mag brighter than L ^*.   We discuss the properties of the complete sample of 14 quasars. Nine of the objects appear to have elliptical host galaxies (all six of the radio-loud quasars in the sample as well as three radio-quiet quasars). Two further radio-quiet quasars appear to lie in disc galaxies. The other three objects (radio-quiet, ultraluminous infrared quasars) all lie in violently interacting systems. The sample as a whole has an average luminosity about 0.8 mag brighter than L ^*, although the radio-loud objects have hosts on average 0.7 mag brighter than the radio-quiet objects.   We compare our results with those from HST imaging of quasars by other authors. Taken together, our observations are in broad agreement with those of Bahcall et al. Radio-loud quasars appear to lie in luminous elliptical galaxies whereas radio-quiet quasars are found to lie in either elliptical or spiral hosts. Host galaxy luminosities (of radio-quiet and radio-loud quasars) are much brighter than would be expected if they followed a Schechter luminosity function.     

Halo masses for optically selected and for radio-loud AGN from clustering and galaxy–galaxy lensing

We compute two-point correlation functions and measure the shear signal due to galaxy–galaxy lensing for 80 000 optically identified and 5700 radio-loud active galactic nuclei (AGN) from Data Release 4 of the Sloan Digital Sky Survey. Halo occupation models are used to estimate halo masses and satellite fractions for these two types of AGN. The large sample size allows us to separate AGN according to the stellar mass of their host galaxies. We study how the halo masses of optical and radio AGN differ from those of the parent population at fixed   M _*  . Halo masses deduced from clustering and from lensing agree satisfactorily. Radio AGN are found in more massive haloes than optical AGN: in our samples, their mean halo masses are  1.6 × 10¹³  and  8 × 10¹¹  h ⁻¹ M_⊙  , respectively. Optical AGN follow the same relation between stellar mass and halo mass as galaxies selected without regard to nuclear properties, but radio-loud AGN deviate significantly from this relation. The dark matter haloes of radio-loud AGN are about twice as massive as those of control galaxies of the same stellar mass. This boost is independent of radio luminosity, and persists even when our analysis is restricted to field galaxies. The large-scale gaseous environment of the galaxy clearly plays a crucial role in producing observable radio emission. The dark matter halo masses that we derive for the AGN in our two samples are in good agreement with recent models in which feedback from radio AGN becomes dominant in haloes where gas cools quasi-statically.     

Environmental dependence of active galactic nuclei activity in the supercluster A901/2

We present XMM data for the supercluster A901/2, at   z ∼ 0.17  , which is combined with deep imaging and 17-band photometric redshifts (from the COMBO-17 survey), two degree field (2dF) spectra and Spitzer 24 μm data, to identify active galactic nuclei (AGN) in the supercluster. The 90 ksec XMM image contains 139 point sources, of which 11 are identified as supercluster AGN with   L _{X(0.5−7.5 keV)} > 1.7 × 10⁴¹ erg cm⁻² s⁻¹  . The host galaxies have   M _R < −20  and only two of eight sources with spectra could have been identified as AGN by the detected optical emission lines. Using a large sample of 795 supercluster galaxies, we define control samples of massive galaxies with no detected AGN. The local environments of the AGN and control samples differ at ≳98 per cent significance. The AGN host galaxies lie predominantly in areas of moderate projected galaxy density and with more local blue galaxies than the control sample, with the exception of one very bright type I AGN very near the centre of a cluster. These environments are similar to, but not limited to, cluster outskirts and blue groups. Despite the large number of potential host galaxies, no AGN are found in regions with the highest galaxy density (excluding some cluster cores where emission from the intra-cluster medium obscures moderate luminosity AGN). AGN are also absent from the areas with lowest galaxy density. We conclude that the prevalence of cluster AGN is linked to their environment.     

Morphology and surface brightness evolution of z∼1.1 radio galaxies

We use K '-band (2.1-μm) imaging to investigate the angular size and morphology of 10 6C radio galaxies, at redshifts 1≤ z ≤1.4. Two radio galaxies appear to be undergoing mergers, another contains, within a single envelope, 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 radio-luminous 3CR galaxies at similar redshifts. This would indicate that the higher mean K -band luminosity of the 3CR galaxies reflects a difference in the size of the host galaxies, and not solely a difference in the power of the active nuclei.
The size–luminosity relation of the z ∼1.1 6C galaxies indicates a 1.0–1.6 mag enhancement of their 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 ∼1.1 show a comparable enhancement in surface brightness. The mean radius of the 6C galaxies suggests that they evolve into ellipticals of L ∼ L * luminosity, and is consistent with their low-redshift counterparts being relatively small FRII galaxies ∼25 times lower in radio luminosity, or small FRI galaxies ∼1000 times lower in radio luminosity. Hence the 6C radio galaxies appear to undergo as much optical and radio evolution as the 3CR galaxies.