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.     

M BH–σ relation in Sloan Digital Sky Survey flat-spectrum radio quasars

The relationship between the black hole mass and velocity dispersion indicated with [O  iii ] linewidth is investigated for a sample of 87 flat-spectrum radio quasars selected from the Sloan Digital Sky Survey Data Release 3 quasar catalogue. We found that the   M _BH−σ_[O III]  relation is different from the Tremaine et al. relation for nearby inactive galaxies, with a larger black hole mass at given velocity dispersion. There is no strong evidence of cosmology evolution in the   M _BH−σ_[O III]  relation up to   z ∼ 0.8  . A significant correlation between the [O  iii ] luminosity and broad-line region (BLR) luminosity is found. When transferring the [O  iii ] luminosity to narrow-line region (NLR) luminosity, the BLR luminosity is, on average, larger than the NLR one by about one order of magnitude. We found a strong correlation between the synchrotron peak luminosity and NLR luminosity, which implies a tight relation between the jet physics and accretion process.     

Radio Luminosity, Black Hole Mass and Eddington Ratio for Quasars from the Sloan Digital Sky Survey

We investigate the MBH-O'.relation for radio-loud quasars with redshift z＜0.83 in Data Release 3 of the Sloan Digital Sky Survey (SDSS).The sample consists of 3772 quasars with better models of the H/~ and [O III] lines and available radio luminosity,including 306 radio-loud quasars,3466 radio-quiet quasars with measured radio luminosity or upper-limit of radio luminosity (181 radio-quiet quasars with measured radio luminosity).The virial supermassive black hole mass (MBH) is calculated from the broad H/line,and the host stellar velocity dispersion (σ*) is traced by the core [O III] gaseous velocity dispersion.The radio luminosity and radio loudness are derived from the FIRST catalog.Our results are as follows:(1) For radio-quiet quasars,we confirm that there is no obvious deviation from the MBH-O".relation defined for inactive galaxies when the uncertainties in MBH and the luminosity bias are concerned.(2) We find that the radio-loud quasars deviate more from the MBH-σ.relation than do the radio-quiet quasars.This deviation is only partly due to a possible cosmological evolution of the MBH-σ* relation and the luminosity bias.(3) The radioluminosity is proportional to M128 0.23-0.16 BH(LBol/LEdd)1.29 0.31-0.24 for radio-quiet quasars and to M3.10 0.6.-0.70(LBol.LEdd)4.18 1.40-1.10 for radio-loud quasars.The weaker dependence of the radio luminosity on the mass and the Eddington ratio for radio-loud quasars shows that other physical effects would account for their radio luminosities,such as the spin of the black hole.     

Measuring the black hole masses of high-redshift quasars

A new technique is presented for determining the black hole masses of high-redshift quasars from optical spectroscopy. The new method utilizes the full-width at half-maximum (FWHM) of the low-ionization Mg  ii emission line and the correlation between the broad-line region (BLR) radius and the continuum luminosity at 3000 Å. Using archival ultraviolet (UV) spectra it is found that the correlation between BLR radius and 3000-Å luminosity is tighter than the established correlation with 5100-Å luminosity. Furthermore, it is found that the correlation between BLR radius and 3000-Å continuum luminosity is consistent with a relation of the form   R _BLR∝λ L ^1/2_λ  , as expected for a constant ionization parameter. Using a sample of objects with broad-line radii determined from reverberation mapping it is shown that the FWHM of Mg  ii and Hβ are consistent with following an exact one-to-one relation, as expected if both Hβ and Mg  ii are emitted at the same radius from the central ionizing source. The resulting virial black hole mass estimator based on rest-frame UV observables is shown to reproduce black hole mass measurements based on reverberation mapping to within a factor of 2.5 (1σ). Finally, the new UV black hole mass estimator is shown to produce identical results to the established optical (Hβ) estimator when applied to 128 intermediate-redshift  (0.3 < z < 0.9)  quasars drawn from the Large Bright Quasar Survey and the radio-selected Molonglo quasar sample. We therefore conclude that the new UV virial black hole mass estimator can be reliably used to estimate the black hole masses of quasars from   z ∼ 0.25  through to the peak epoch of quasar activity at   z ∼ 2.5  via optical spectroscopy alone.     

Modelling the cosmological co-evolution of supermassive black holes and galaxies – II. The clustering of quasars and their dark environment

We use semi-analytic modelling on top of the Millennium simulation to study the joint formation of galaxies and their embedded supermassive black holes. Our goal is to test scenarios in which black hole accretion and quasar activity are triggered by galaxy mergers, and to constrain different models for the light curves associated with individual quasar events. In the present work, we focus on studying the spatial distribution of simulated quasars. At all luminosities, we find that the simulated quasar two-point correlation function is fit well by a single power law in the range  0.5 ≲ r ≲ 20  h ⁻¹ Mpc  , but its normalization is a strong function of redshift. When we select only quasars with luminosities within the range typically accessible by today's quasar surveys, their clustering strength depends only weakly on luminosity, in agreement with observations. This holds independently of the assumed light-curve model, since bright quasars are black holes accreting close to the Eddington limit, and are hosted by dark matter haloes with a narrow mass range of a few  10¹²  h ⁻¹ M_⊙  . Therefore, the clustering of bright quasars cannot be used to disentangle light-curve models, but such a discrimination would become possible if the observational samples can be pushed to significantly fainter limits. Overall, our clustering results for the simulated quasar population agree rather well with observations, lending support to the conjecture that galaxy mergers could be the main physical process responsible for triggering black hole accretion and quasar activity.     

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.     

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.     

The quasar epoch and the stellar ages of early-type galaxies

We investigate the hypothesis that quasars formed together with the stellar populations of early-type galaxies. This hypothesis – in conjunction with the stellar ages of early-type galaxies from population synthesis models, the relation of black hole mass to bulge velocity dispersion, and the velocity dispersion distribution of spheroids from the Sloan Digital Sky Survey – completely determines the cosmic accretion history of supermassive black holes and the redshift evolution of the characteristic luminosity. On the other hand, the precise shape of the luminosity function of quasars depends on the light curve of quasars and – in the optical, but not so much in X-rays – on the covering factor of the dust surrounding the active nucleus. We find a plausible set of assumptions for which the coeval formation of supermassive black holes and elliptical galaxies is in good agreement with the observed B -band and X-ray luminosity functions of quasars.     

Radio-quiet quasar environments at 0.5 ≤ z ≤ 0.8

We quantify the galaxy environments around a sample of 0.5≤ z ≤0.8 radio-quiet quasars using the amplitude of the spatial galaxy–quasar correlation function, B _gq. The quasars exist in a wide variety of environments; some sources are located in clusters as rich as Abell class 1–2 clusters, whereas others exist in environments comparable to the field. We find that, on average, the quasars prefer poorer clusters of ≈Abell class 0, which suggests that quasars are biased tracers of mass compared with galaxies. The mean B _gq for the sample is found to be indistinguishable from the mean amplitude for a sample of radio-loud quasars matched in redshift and optical luminosity. These observations are consistent with recent studies of the hosts of radio-quiet quasars at low to intermediate redshifts, and suggest that the mechanism for the production of powerful radio jets in radio-loud quasars is controlled by processes deep within the active galactic nucleus itself, and is unrelated to the nature of the hosts or their environments.     

A unified model for the evolution of galaxies and quasars

We incorporate a simple scheme for the growth of supermassive black holes into semi-analytic models that follow the formation and evolution of galaxies in a cold dark matter-dominated Universe. We assume that supermassive black holes are formed and fuelled during major mergers. If two galaxies of comparable mass merge, their central black holes coalesce and a few per cent of the gas in the merger remnant is accreted by the new black hole over a time-scale of a few times 10⁷ yr. With these simple assumptions, our model not only fits many aspects of the observed evolution of galaxies, but also reproduces quantitatively the observed relation between bulge luminosity and black hole mass in nearby galaxies, the strong evolution of the quasar population with redshift, and the relation between the luminosities of nearby quasars and those of their host galaxies. The strong decline in the number density of quasars from z ∼2 to z =0 is a result of the combination of three effects: (i) a decrease in the merging rate; (ii) a decrease in the amount of cold gas available to fuel black holes, and (iii) an increase in the time-scale for gas accretion. The predicted decline in the total content of cold gas in galaxies is consistent with that inferred from observations of damped Ly α systems. Our results strongly suggest that the evolution of supermassive black holes, quasars and starburst galaxies is inextricably linked to the hierarchical build-up of galaxies.     

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 quasar-frame velocity distribution of narrow C iv absorbers

We report on a survey for narrow (full widths at half-minimum <600 km s⁻¹) C  iv absorption lines in a sample of bright quasars at redshifts  1.8 ≤ z < 2.25  in the Sloan Digital Sky Survey. Our main goal is to understand the relationship of narrow C  iv absorbers to quasar outflows and, more generally, to quasar environments. We determine velocity zero-points using the broad Mg  ii emission line, and then measure the absorbers' quasar-frame velocity distribution. We examine the distribution of lines arising in quasar outflows by subtracting model fits to the contributions from cosmologically intervening absorbers and absorption due to the quasar host galaxy or cluster environment. We find that a substantial number (  ≥43 ± 6  per cent) of absorbers with   W ^λ1548₀ > 0.3  Å in the velocity range  +750 ≲ v ≲+ 12 000  km s⁻¹ are intrinsic to the active galactic nucleus outflow. This 'outflow fraction' peaks near   v =+2000  km s⁻¹ with a value of   f _outflow≃ 0.81 ± 0.13  . At velocities below   v ≈+ 2000  km s⁻¹, the incidence of outflowing systems drops, possibly due to geometric effects or to the over-ionization of gas that is nearer the accretion disc. Furthermore, we find that outflow absorbers are on average broader and stronger than cosmologically intervening systems. Finally, we find that ∼14 per cent of the quasars in our sample exhibit narrow, outflowing C  iv absorption with   W ^λ1548₀ > 0.3  Å, slightly larger than that for broad absorption line systems.     

The optical luminosity evolution function of quasars. I. A statistical evolution function that includes the large redshift range above redshift 2.2

We have redetermined the parameters in the Heisler-Ostriker quasar evolution formula so as to make it applicable to a larger redshift range. The form of the formula is appropriate because, physically, any variation in the mass of the central black hole and the accretion rate will cause a luminosity evolution while the birth of quasars will imply a luminosity-dependent density evolution. With our new parameters, we can 1) reproduce the quasar magnitude count down to magnitude 22.0, 2) reproduce the quasar redshift distribution in intervals of ΔZ = 0.1 and 3) substantially lessen the excess count for Z > 2 redicted by previous formulae.     

A Study of the Black Hole Mass and Accretion Rate Distributions of AGNs

The redshift, central black hole mass and accretion rate are important parameters when studying the AGN evolution. The central black hole masses for 172 quasars and Seyfert galaxies are calculated in this paper using the reverberation mapping method. The distributions of central black hole masses, redshifts and the Eddington accretion rates are analyzed, to verify the transition from the quasar to the Seyfert galaxy in the course of evolution.     

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.     

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.     

Fuelling quasars with hot gas

We consider a model for quasar formation in which massive black holes are formed and fuelled largely by the accretion of hot gas during the process of galaxy formation. In standard hierarchical collapse models, objects about the size of normal galaxies and larger form a dense hot atmosphere when they collapse. We show that if such an atmosphere forms a nearly 'maximal' cooling flow, then a central black hole can accrete at close to its Eddington limit. This leads to exponential growth of a seed black hole, resulting in a quasar in some cases. In this model, the first quasars form soon after the first collapses to produce hot gas. The hot gas is depleted as time progresses, mostly by cooling, so that the accretion rate eventually falls below the threshold for advection-dominated accretion, at which stage radiative efficiency plummets and any quasar turns off. A simple implementation of this model, incorporated into a semi-analytical model for galaxy formation, overproduces quasars when compared with observed luminosity functions, but is consistent with models of the X-ray background, which indicate that most accretion is obscured. It produces few quasars at high redshift owing to the lack of time needed to grow massive black holes. Quasar fuelling by hot gas provides a minimum level, sufficient to power most quasars at redshifts between one and two, to which other sources of fuel can be added. The results are sensitive to feedback effects, such as might result from radio jets and other outflows.     

Radio loudness of high-redshift (z≥3) quasars from the Sloan Digital Sky Survey

Using a homogenous sample of 1962 quasars with redshift 3.0≤z≲4.5 drawn from the Sloan Digital Sky Survey (SDSS), we study the relationships between radio loudness, virial black hole (BH) mass and Eddington ratio (accretion rate relative to the Eddington limit). For the radio-detected objects, we find a significant (>99.5 per cent) anticorrelation between radio loudness R parameter and BH mass, consistent with previous studies of low-redshift radio-loud quasars. The truly radio-loud quasars (R>30) are found to be confined to M _BH≲10¹⁰ M_⊙ within our sample. We also find that R is only weakly correlated with Eddington ratio L _bol/L _Edd. Combined with previous results on the low-redshift R−L _bol/L _Edd relation, this result indicates no strong L _bol/L _Edd dependence of R at L _bol/L _Edd≳10⁻². On the other hand, the large scatter in these relationships suggests that other physical properties such as BH spin and quasar clustering must also play an important role in quasar radio emission.     

Unification of Radio Galaxies and their Accretion Jet Properties

We investigate the relation between black hole mass, M _bh, and jet power, Q _jet, for a sample of BL Lacs and radio quasars. We find that BL Lacs are separated from radio quasars by the FR I/II dividing line in M _bh–Q _jet plane, which strongly supports the unification scheme of FR I/BL Lac and FR II/radio quasar. The Eddington ratio distribution of BL Lacs and radio quasars exhibits a bimodal nature with a rough division at L _bol/L _Edd~0.01, which imply that they may have different accretion modes. We calculate the jet power extracted from advection-dominated accretion flow (ADAF), and find that it requires dimensionless angular momentum of black hole j???0.9???0.99 to reproduce the dividing line between FR I/II or BL Lac/radio quasar if dimensionless accretion rate $\dot{m}=0.01$ is adopted, which is required by the above bimodal distribution of Eddington ratios. Our results suggest that black holes in radio galaxies are rapidly spinning.