The evolution of star formation in quasar host galaxies

We have used far-infrared data from IRAS , Infrared Space Observatory ( ISO ), Spitzer Wide-Area Infrared Extragalactic (SWIRE), Submillimetre Common User Bolometer Array (SCUBA) and Max-Planck Millimetre Bolometer (MAMBO) to constrain statistically the mean far-infrared luminosities of quasars. Our quasar compilation at redshifts  0 < z < 6.5  and I -band luminosities  −20 < I _AB < −32  is the first to distinguish evolution from quasar luminosity dependence in such a study. We carefully cross-calibrate IRAS against Spitzer and ISO , finding evidence that IRAS 100-μm fluxes at <1 Jy are overestimated by ∼30 per cent. We find evidence for a correlation between star formation in quasar hosts and the quasar optical luminosities, varying as star formation rate (SFR)  ∝ L ^0.44±0.07_opt  at any fixed redshift below   z = 2  . We also find evidence for evolution of the mean SFR in quasar host galaxies, scaling as  (1 + z )^1.6±0.3  at   z < 2  for any fixed quasar I -band absolute magnitude fainter than −28. We find no evidence for any correlation between SFR and black hole mass at  0.5 < z < 4  . Our data are consistent with feedback from black hole accretion regulating stellar mass assembly at all redshifts.     

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.     

The prevalence of Fanaroff–Riley type I radio quasars

We present deep, multi-Very Large Array configuration radio images for a set of 18 quasars, having redshifts between 0.36 and 2.5, from the 7C quasar survey. Approximately one quarter of these quasars have Fanaroff–Riley type I (FR I) type twin-jet structures and the remainder are a broad range of wide angle tail, fat double, classical double, core-jet and hybrid sources. These images demonstrate that FR I quasars are prevalent in the Universe, rather than non-existent as had been suggested in the literature prior to the serendipitous discovery of the first FR I quasar a few years ago, the optically powerful 'radio-quiet' quasar E 1821+643.
Some of the FR I quasars have radio luminosities exceeding the traditional FR I/FR II break luminosity; however, we find no evidence for FR II quasars with luminosities significantly below the break. We consider whether the existence of such high-luminosity FR I structures is due to the increasingly inhomogeneous environments in the higher redshift Universe.     

The chemical evolution of the Universe – I. High column density absorbers

We construct a simple, robust model of the chemical evolution of galaxies from high to low redshift, and apply it to published observations of damped Lyman α quasar absorption line systems (DLAs). The elementary model assumes quiescent star formation and isolated galaxies (no interactions, mergers or gas flows). We consider the influence of dust and chemical gradients in the galaxies, and hence explore the selection effects in quasar surveys. We fit individual DLA systems to predict some observable properties of the absorbing galaxies, and also indicate the expected redshift behaviour of chemical element ratios involving nucleosynthetic time delays.
Despite its simplicity, our 'monolithic collapse' model gives a good account of the distribution and evolution of the metallicity and column density of DLAs, and of the evolution of the global star formation rate and gas density below redshifts z ∼3. However, from the comparison of DLA observations with our model, it is clear that star formation rates at higher redshifts ( z >3) are enhanced. Galaxy interactions and mergers, and gas flows very probably play a major role.     

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.     

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.     

Contribution of galaxies to the background hydrogen-ionizing flux

We estimate the evolution of the contribution of galaxies to the cosmic background flux at 912 Å by means of a semi-analytic model of galaxy formation and evolution. Such modelling has been quite successful in reproducing the optical properties of galaxies. We assume that high-redshift damped Lyman α systems are the progenitors of present-day galaxies, and we design a series of models that are consistent with the evolution of cosmic comoving emissivities in the available near-infrared, optical, ultraviolet and far-infrared bands along with the evolution of the neutral hydrogen content and average metallicity of damped Lyman α systems. We use these models to compute the galactic contribution to the Lyman-limit emissivity and background flux for 0 ≃  z  ≤ 4. We take into account the absorption of Lyman-limit photons by H  I and dust in the interstellar medium of the galaxies. We find that the background Lyman-limit flux due to galaxies might dominate (or be comparable to) the contribution from quasars at almost all redshifts if the absorption by H  I in the interstellar medium is neglected. Such H  I absorption would result in a severe diminishing of this flux — by almost three orders of magnitude at high redshifts and by one to two orders at z  ≃ 0. Though the resulting galaxy flux is completely negligible at high redshifts, it is comparable to the quasar flux at z  ≃ 0.     

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.     

Cosmological evolution and hierarchical galaxy formation

We calculate the rate at which dark matter haloes merge to form higher mass systems. Two complementary derivations using Press–Schechter theory are given, both of which result in the same equation for the formation rate. First, a derivation using the properties of the Brownian random walks within the framework of Press–Schechter theory is presented. We then use Bayes' theorem to obtain the same result from the standard Press–Schechter mass function. The rate obtained is shown to be in good agreement with results from Monte Carlo and N -body simulations. We illustrate the usefulness of this formula by calculating the expected cosmological evolution in the rate of star formation that is due to short-lived, merger-induced starbursts. The calculated evolution is well-matched to the observed evolution in ultraviolet luminosity density, in contrast to the lower rates of evolution that are derived from semi-analytic models that do not include a dominant contribution from starbursts. Hence we suggest that the bulk of the observed ultraviolet starlight at z >1 arises from merger-induced starbursts. Finally, we show that a simple merging-halo model can also account for the bulk of the observed evolution in the comoving quasar space density.     

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.     

3C 254: MERLIN observations of a highly asymmetric quasar

Multifrequency, high-resolution radio observations of the quasar 3C 254 using the Multi-Element Radio-linked Interferometer Network (MERLIN) are presented. The quasar has a highly asymmetric radio structure, with the eastern component of the double-lobed structure being much closer to the nucleus and significantly less polarized than the western one. However, the two lobes are more symmetric in their total flux densities. The observations show the detailed structure of the hotspots which are very different on opposite sides of the radio core, reveal no radio jet and suggest that the oppositely directed jets may be intrinsically asymmetric.     

Downsizing of supermassive black holes from the SDSS quasar survey – II. Extension to z∼ 4

Starting from the quasar sample of the Sloan Digital Sky Survey (SDSS) for which the C  iv line is observed, we use an analysis scheme to derive the z -dependence of the maximum mass of active black holes, which overcomes the problems related to the Malmquist bias. The same procedure is applied to the low-redshift sample of SDSS quasars for which Hβ measurements are available. Combining with the results from the previously studied Mg  ii sample, we find that the maximum mass of the quasar population increases as  (1 + z )^1.64±0.04  in the redshift range  0.1 z 4  , which includes the epoch of maximum quasar activity.     

ISO observations of the dusty quasar BR 1202−0725

We present mid- and far-infrared photometry of the high-redshift     dusty quasar BR 1202−0725. The quasar was detected in the near-infrared, at a flux level     consistent with an average radio-quiet quasar at its redshift. Only upper limits for the emission were obtained in the far-infrared. These upper limits, when combined with data from ground-based telescopes, are the first direct evidence for a turnover in the far-infrared emission, and hence confirm that a blackbody dominates the spectral energy distribution at far-infrared wavelengths. This blackbody is most probably cool dust, constrained to have a temperature below 80 K, for a β of 1.5.     

Constraining quasar host halo masses with the strength of nearby Lyα forest absorption

Using cosmological hydrodynamic simulations, we measure the mean transmitted flux in the Lyα forest for quasar sightlines that pass near a foreground quasar. We find that the trend of absorption with pixel quasar separation distance can be fitted using a simple power-law form including the usual correlation function parameters r ₀ and γ, so that     . From the simulations, we find the relation between r ₀ and quasar host mass, and formulate this as a way to estimate quasar host dark matter halo masses, quantifying uncertainties due to cosmological and IGM parameters, and redshift errors. With this method, we examine data for ∼9000 quasars from the Sloan Digital Sky Survey (SDSS) Data Release 5, assuming that the effect of ionizing radiation from quasars (the so-called transverse proximity effect) is unimportant (no evidence for it is seen in the data). We find that the best-fitting host halo mass for SDSS quasars with mean redshift z = 3 and absolute G -band magnitude −27.5 is  log  M /M_⊙= 12.68^+0.81_−0.67  . We also use the Lyman-Break Galaxy (LBG) and Lyα forest data of Adelberger et al. in a similar fashion to constrain the halo mass of LBGs to be  log₁₀  M /M_⊙= 11.41^+0.54_−0.59  , a factor of ∼20 lower than the bright quasars. In addition, we study the redshift distortions of the Lyα forest around quasars, using the simulations. We use the quadrupole to monopole ratio of the quasar Lyα forest correlation function as a measure of the squashing effect. We find its dependence on halo mass difficult to measure, but find that it may be useful for constraining cosmic geometry.     

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.     

Formation and evolution of elliptical galaxies and QSO activity

We present the results of a numerical code that combines multi-zone chemical evolution with 1D hydrodynamics to follow in detail the evolution and radial behaviour of gas and stars during the formation of elliptical galaxies. We use the model to explore the links between the evolution and formation of elliptical galaxies and QSO activity. The knowledge of the radial gas flows in the galaxy allows us to trace metallicity gradients, and, in particular, the formation of a high-metallicity core in ellipticals. The high-metallicity core is formed soon enough to explain the metal abundances inferred in high-redshift quasars. The star formation rate and the subsequent feedback regulate the episodes of wind, outflow and cooling flow, thus affecting the recycling of the gas and the chemical enrichment of the intergalactic medium. The evolution of the galaxy shows several stages, some of which are characterized by a complex flow pattern, with inflow in some regions and outflow in other regions. All models, however, exhibit during their late evolution a galactic wind at the outer boundary and, during their early evolution, an inflow towards the galactic nucleus. The characteristics of the inner inflow could explain the bolometric luminosity of a quasar lodged at the galactic centre as well as the evolution of the optical luminosity of 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.     

Black hole growth in hierarchical galaxy formation

We incorporate a model for black hole growth during galaxy mergers into the semi-analytical galaxy formation model based on ΛCDM proposed by Baugh et al. Our black hole model has one free parameter, which we set by matching the observed zero-point of the local correlation between black hole mass and bulge luminosity. We present predictions for the evolution with redshift of the relationships between black hole mass and bulge properties. Our simulations reproduce the evolution of the optical luminosity function of quasars. We study the demographics of the black hole population and address the issue of how black holes acquire their mass. We find that the direct accretion of cold gas during starbursts is an important growth mechanism for lower mass black holes and at high redshift. On the other hand, the re-assembly of pre-existing black hole mass into larger units via merging dominates the growth of more massive black holes at low redshift. This prediction could be tested by future gravitational wave experiments. As redshift decreases, progressively less massive black holes have the highest fractional growth rates, in line with recent claims of 'downsizing' in quasar activity.