Simulating wide-field quasar surveys from the optical to near-infrared

A number of deep, wide-field, near-infrared (NIR) surveys employing new infrared cameras on 4-m class telescopes are about to commence. These surveys have the potential to determine the fraction of luminous dust-obscured quasars that may have eluded surveys undertaken at optical wavelengths. In order to understand the new observations it is essential to make accurate predictions of surface densities and number–redshift relations for unobscured quasars in the NIR based on information from surveys at shorter wavelengths. The accuracy of the predictions depends critically on a number of key components. The commonly used single power-law representation for quasar spectral energy distributions (SEDs) is inadequate and the use of an SED incorporating the upturn in continuum flux at  λ∼ 12 000 Å  is essential. The presence of quasar host galaxies is particularly important over the rest-frame wavelength interval  8000 < λ < 16 000 Å  and we provide an empirical determination of the magnitude distribution of host galaxies using a low-redshift sample of quasars from the Sloan Digital Sky Survey Data Release 3 quasar catalogue. A range of models for the dependence of host galaxy luminosity on quasar luminosity is investigated, along with the implications for the NIR surveys. Even adopting a conservative model for the behaviour of host galaxy luminosity the number counts for shallow surveys in the K band increase by a factor of 2. The degree of morphological selection applied to define candidate quasar samples in the NIR is found to be an important factor in determining the fraction of the quasar population included in such samples.     

On quasar host galaxies as tests of non-cosmological redshifts

Despite a general consensus in the astronomical community that all quasars are located at the distances implied by their redshifts, a number of observations still challenge this interpretation, possibly indicating that some subpopulation of quasars may harbour significant redshift components not related to the expansion of the Universe. It has been suggested that these objects may have been ejected from local galaxies and are likely to evolve into new galaxies themselves. Here, a test of such exotic scenarios is proposed, based on the spectral energy distribution of the galaxies hosting quasars with suspected ejection origin. Provided that the time-scales over which the ejected objects manifest themselves as quasars is short, one would in the framework of the ejection scenarios expect to find either no quasar host galaxy, a pseudo-host consisting of gas ionized by the quasar, or a host galaxy consisting of young stars only. It is argued that the spectral energy distributions corresponding to the last two options should differ significantly from that of most quasar host galaxies detected at low redshift so far, thus providing a potential test of the claimed existence of ejected quasars. A minimal implementation of this test, involving optical and near-infrared broad-band photometry, is suggested.     

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.     

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.     

Infrared photometry of z ∼ 1 3C quasars

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

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.     

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.     

A population of high-redshift type 2 quasars – I. Selection criteria and optical spectra

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

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.     

Microlensing of broad absorption line quasars

The physical nature of the material responsible for the high-velocity broad absorption line (BAL) features seen in a small fraction of quasar spectra has been the subject of debate since their discovery. This has been especially compounded by the lack of observational probes of the absorbing region. In this paper we examine the role of 'microlenses' in external galaxies in the observed variability in the profiles of BALs in multiply imaged quasars. Utilizing realistic models for both the BAL region and the action of an ensemble of microlensing masses, we demonstrate that stars at cosmological distances can provide an important probe of the physical state and structure of material at the heart of these complex systems. Applying these results to the macrolensed BAL quasar system H1413+117, the observed spectral variations are readily reproduced, but without the fine-tuning requirements of earlier studies which employ more simplistic models.     

Quasar–galaxy associations revisited

Gravitational lensing predicts an enhancement of the density of bright, distant quasi-stellar objects (QSOs) around foreground galaxies. We measure this QSO–galaxy correlation w _qg for two complete samples of radio-loud quasars, the southern 1 Jy and half-Jansky samples. The existence of a positive correlation between     quasars and     galaxies is confirmed at a     per cent significance level (>99.9 per cent if previous measurements on the northern hemisphere are included). A comparison with the results obtained for incomplete quasar catalogues (e.g. the Veron-Cetty and Veron compilation) suggests the existence of an 'identification bias', which spuriously increases the estimated amplitude of w _qg for incomplete samples. This effect may explain many of the surprisingly strong quasar–galaxy associations found in the literature. Nevertheless, the value of w _qg that we measure in our complete catalogues is still considerably higher than the predictions from weak lensing. We consider two effects which could help to explain this discrepancy: galactic dust extinction and strong lensing.     

Color-Redshift Relations and Photometric Redshift Estimations of Quasars in Large Sky Surveys

With a recently constructed composite quasar spectrum and the X2 minimization technique, we describe a general method for estimating the photometric redshifts of a large sample of quasars by deriving theoretical color-redshift relations and comparing the theoretical colors with the observed ones. We estimated the photometric redshifts from the 5-band SDSS photometric data of 18678 quasars in the first major data release of SDSS and compared them with their spectroscopic redshifts. The difference is less than 0.1 for 47% of the quasars and less than 0.2 for 68%. Based on the calculation of the theoretical color-color diagrams of stars, galaxies and quasars both on the SDSS system and on the BATC system, we expect that we would be able to select candidates of high redshift quasars more efficaciously with the latter than with the former, provided the BATC survey can detect objects with magnitudes fainter than 21.     

The dependence of quasar variability on black hole mass

In order to investigate the dependence of quasar variability on fundamental physical parameters like black hole mass, we have matched quasars from the Quasar Equatorial Survey Team, Phase 1 (QUEST1) variability survey with broad-lined objects from the Sloan Digital Sky Survey. The matched sample contains ≈100 quasars, and the Sloan spectra are used to estimate black hole masses and bolometric luminosities. Variability amplitudes are measured from the QUEST1 light curves. We find that black hole mass correlates with several measures of the variability amplitude at the 99 per cent significance level or better. The correlation does not appear to be caused by obvious selection effects inherent to flux-limited quasar samples, host galaxy contamination or other well-known correlations between quasar variability and luminosity/redshift. We evaluate variability as a function of rest-frame time lag using structure functions and find further support for the variability–black hole mass correlation. The correlation is strongest for time lags of the order of a few months up to the QUEST1 maximum temporal resolution of ≈2 yr, and may provide important clues for understanding the long-standing problem of the origin of quasar optical variability. We discuss whether our result is a manifestation of a relation between characteristic variability time-scale and black hole mass, where the variability time-scale is typical for accretion disc thermal time-scales, but find little support for this. Our favoured explanation is that more massive black holes have larger variability amplitudes, and we highlight the need for larger samples with more complete temporal sampling to test the robustness of this result.     

The far-infrared–submillimetre spectral energy distribution of high-redshift quasars

We combine photometric observations of high-redshift     quasars, obtained at submillimetre to millimetre wavelengths, to obtain a mean far-infrared (rest-frame) spectral energy distribution (SED) of the thermal emission from dust, parametrized by a single temperature ( T ) and power-law emissivity index ( β ). The best-fitting values are     and     . Our method exploits the redshift spread of this set of quasars, which allows us to sample the SED at a larger number of rest wavelengths than is possible for a single object: the wavelength range extends down to ∼60 μm, and therefore samples the turnover in the greybody curve for these temperatures. This parametrization is of use to any studies that extrapolate from a flux at a single wavelength, for example to infer dust masses and far-infrared luminosities.
We interpret the cool, submillimetre component as arising from dust heated by star formation in the host galaxy of the quasar, although we do not exclude the presence of dust heated directly by the active galactic nucleus (AGN). Applying the mean SED to the data, we derive consistent star formation rates ∼1000 M_⊙ yr⁻¹ and dust masses ∼10⁹ M_⊙, and investigate a simple scheme of AGN and host galaxy co-evolution to account for these quantities. The time-scale for formation of the host galaxy is     , and the luminous quasar phase occurs towards the end of this period, just before the reservoir of cold gas is depleted. Given the youth of the Universe at     (1.6 Gyr), the coexistence of a massive black hole and a luminous starburst at high redshifts is a powerful constraint on models of quasar host galaxy formation.     

Far-infrared-loud quasars – I. Disturbed and quiescent quasars in the PG survey

We use host galaxy imaging studies of the PG quasar survey to compare the far-infrared (FIR) properties of quasars with disturbed and undisturbed host galaxies. By using survival analysis, we show that the quasars with disturbed host galaxies, with morphologies classified from a homogeneous data set, have a 60-μm luminosity distribution that is different from that of those with undisturbed hosts with >97 per cent confidence. For morphological classifications using an inhomogeneous data set, including HST data for some objects, this confidence rises to >99 per cent confidence. The mean 60-μm luminosity of the disturbed-host quasars is several times greater than that of the undisturbed-host quasars. However, possible biases in the PG survey might affect these conclusions. Our results are interpreted as supporting the idea that quasars are related to at least some ultraluminous infrared galaxies. We discuss the implications of this result for studies of quasar and galaxy evolution.     

Red synchrotron jets in Parkes quasars

We present model fits to spectral energy distributions in the optical and near-infrared of >100 flat-spectrum radio quasars from the Parkes Half-Jansky Flat-spectrum Sample. We find that ∼40 per cent of the sources have power-law spectral energy distributions (SEDs), while a similar number show evidence for two primary components: a blue power law and optical synchrotron emission. The blue power law is similar to the dominant component observed in the spectra of optically selected quasars. There is strong evidence that the synchrotron component has a turnover in the ultraviolet–optical rest frame of the spectrum. In the remaining sources, it is likely that the synchrotron peaks at longer wavelengths. This mixture of two components is supported by optical polarization measurements in a subgroup of the sources. The sources with power-law SEDs show evidence for an excess number of red power-law slopes compared with optically selected quasars. There are additional spectral components in some of the sources, such as dust and the underlying galaxy, which have not been considered here.     

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.     

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.     

The galaxy environment of a quasar at z= 1.226: a possible cluster merger

We have conducted ultra-deep optical and deep near-infrared observations of a field around the z =1.226 radio-quiet quasar 104420.8+055739 from the Clowes–Campusano LQG of 18 quasars at z ∼1.3, in search of associated galaxy clustering. Galaxies at these redshifts are distinguished by their extremely red colours, with I − K >3.75, and we find a factor ∼11 overdensity of such galaxies in a 2.25×2.25 arcmin² field centred on the quasar. In particular, we find 15–18 galaxies that have colours consistent with being a population of passively evolving massive ellipticals at the quasar redshift. They form 'fingers' in the V − K K , I − K K colour–magnitude plots at V − K ≃6.9, I − K ≃4.3 comparable to the red sequences observed in other z ≃1.2 clusters. We find suggestive evidence for substructure among the red sequence galaxies in the K image, in the form of two compact groups, 40 arcsec to the north, and 60 arcsec to the south-east of the quasar. An examination of the wider optical images indicates that this substructure is significant, and that the clustering extends to form a large-scale structure 2–3  h ⁻¹ Mpc across. We find evidence for a high (≳50 per cent) fraction of blue galaxies in this system, in the form of 15–20 'red outlier' galaxies with I − K >3.75 and V − I <2.00, which we suggest are dusty, star-forming galaxies at the quasar redshift. Within 30 arcsec of the quasar we find a concentration of blue ( V − I <1) galaxies in a band that bisects the two groups of red sequence galaxies. This band of blue galaxies is presumed to correspond to a region of enhanced star formation. We explain this distribution of galaxies as the early stages of a cluster merger which has triggered both the star formation and the quasar.     

Near-infrared spectroscopy of powerful compact steep-spectrum radio sources

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