The distribution of supermassive black holes in the nuclei of nearby galaxies

The growth of supermassive black holes by merging and accretion in hierarchical models of galaxy formation is studied by means of Monte Carlo simulations. A tight linear relation between masses of black holes and masses of bulges arises if the mass accreted by supermassive black holes scales linearly with the mass-forming stars and if the redshift evolution of mass accretion tracks closely that of star formation. Differences in redshift evolution between black hole accretion and star formation introduce a considerable scatter in this relation. A non-linear relation between black hole accretion and star formation results in a non-linear relation between masses of remnant black holes and masses of bulges. The relation of black hole mass to bulge luminosity observed in nearby galaxies and its scatter are reproduced reasonably well by models in which black hole accretion and star formation are linearly related but do not track each other in redshift. This suggests that a common mechanism determines the efficiency for black hole accretion and the efficiency for star formation, especially for bright bulges.     

On the evolution of the black hole: spheroid mass ratio

We present the results of a study which uses the 3C RR sample of radio-loud active galactic nuclei to investigate the evolution of the black hole:spheroid mass ratio in the most massive early-type galaxies from  0 < z < 2  . Radio-loud unification is exploited to obtain virial (linewidth) black hole mass estimates from the 3C RR quasars, and stellar mass estimates from the 3C RR radio galaxies, thereby providing black hole and stellar mass estimates for a single population of early-type galaxies. At low redshift  ( z ≲ 1)  , the 3C RR sample is consistent with a black hole:spheroid mass ratio of   M _bh/ M _sph≃ 0.002  , in good agreement with that observed locally for quiescent galaxies of similar stellar mass  ( M _sph≃ 5 × 10¹¹ M_⊙)  . However, over the redshift interval  0 < z < 2  the 3C RR black hole:spheroid mass ratio is found to evolve as   M _bh/ M _sph∝ (1 + z )^2.07±0.76  , reaching   M _bh/ M _sph≃ 0.008  by redshift   z ≃ 2  . This evolution is found to be inconsistent with the local black hole:spheroid mass ratio remaining constant at a moderately significant level (98 per cent). If confirmed, the detection of evolution in the 3C RR black hole:spheroid mass ratio further strengthens the evidence that, at least for massive early-type galaxies, the growth of the central supermassive black hole may be completed before that of the host spheroid.     

The evolution of massive black hole seeds

We investigate the evolution of high-redshift seed black hole masses at late times and their observational signatures. The massive black hole seeds studied here form at extremely high redshifts from the direct collapse of pre-galactic gas discs. Populating dark matter haloes with seeds formed in this way, we follow the mass assembly of these black holes to the present time using a Monte Carlo merger tree. Using this machinery, we predict the black hole mass function at high redshifts and at the present time, the integrated mass density of black holes and the luminosity function of accreting black holes as a function of redshift. These predictions are made for a set of three seed models with varying black hole formation efficiency. Given the accuracy of present observational constraints, all three models can be adequately fitted. Discrimination between the models appears predominantly at the low-mass end of the present-day black hole mass function which is not observationally well constrained. However, all our models predict that low surface brightness, bulgeless galaxies with large discs are least likely to be sites for the formation of massive seed black holes at high redshifts. The efficiency of seed formation at high redshifts has a direct influence on the black hole occupation fraction in galaxies at   z = 0  . This effect is more pronounced for low-mass galaxies. This is the key discriminant between the models studied here and the Population III remnant seed model. We find that there exist a population of low-mass galaxies that do not host nuclear black holes. Our prediction of the shape of the M _BH–σ relation at the low-mass end is in agreement with the recent observational determination from the census of low-mass galaxies in the Virgo cluster.     

Cosmological physics with black holes (and possibly white dwarfs)

The notion that microparsec-scale black holes can be used to probe gigaparsec-scale physics may seem counterintuitive, at first. Yet, the gravitational observatory LISA will detect cosmologically-distant coalescing pairs of massive black holes, accurately measure their luminosity distance and help identify an electromagnetic counterpart or a host galaxy. A wide variety of new black hole studies and a gravitational version of Hubble’s diagram become possible, if host galaxies are successfully identified. Furthermore, if dark energy is a manifestation of large-scale modified gravity, deviations from general relativistic expectations could become apparent in a gravitational signal propagated over cosmological scales, especially when compared to the electromagnetic signal from a same source. Finally, since inspirals of white dwarfs into massive black holes at cosmological distances may permit pre-merger localizations, we suggest that careful monitoring of these events and any associated electromagnetic counterpart could lead to high-precision cosmological measurements with LISA.     

High-redshift galaxies, their active nuclei and central black holes

We demonstrate that the luminosity function of the recently detected population of actively star-forming galaxies at redshift z  = 3 and the B -band luminosity function of quasi-stellar objects (QSOs) at the same redshift can both be matched with the mass function of dark matter haloes predicted by standard variants of hierarchical cosmogonies for lifetimes of optically bright QSOs anywhere in the range 10⁶ to 10⁸ yr. There is a strong correlation between the lifetime and the required degree of non-linearity in the relation between black hole and halo mass. We suggest that the mass of supermassive black holes may be limited by the back-reaction of the emitted energy on the accretion flow in a self-gravitating disc. This would imply a relation of black hole to halo mass of the form M _bh ∝  v ⁵_halo ∝  M ^5/3_halo and a typical duration of the optically bright QSO phase of a few times 10⁷ yr. The high integrated mass density of black holes inferred from recent black hole mass estimates in nearby galaxies may indicate that the overall efficiency of supermassive black holes for producing blue light is smaller than previously assumed. We discuss three possible accretion modes with low optical emission efficiency: (i) accretion at far above the Eddington rate, (ii) accretion obscured by dust, and (iii) accretion below the critical rate leading to an advection-dominated accretion flow lasting for a Hubble time. We further argue that accretion with low optical efficiency might be closely related to the origin of the hard X-ray background and that the ionizing background might be progressively dominated by stars rather than QSOs at higher redshift.     

Is there an upper limit to black hole masses?

We make a case for the existence for ultra-massive black holes (UMBHs) in the Universe, but argue that there exists a likely upper limit to black hole (BH) masses of the order of   M ∼ 10¹⁰ M_⊙  . We show that there are three strong lines of argument that predicate the existence of UMBHs: (i) expected as a natural extension of the observed BH mass bulge luminosity relation, when extrapolated to the bulge luminosities of bright central galaxies in clusters; (ii) new predictions for the mass function of seed BHs at high redshifts predict that growth via accretion or merger-induced accretion inevitably leads to the existence of rare UMBHs at late times; (iii) the local mass function of BHs computed from the observed X-ray luminosity functions of active galactic nuclei predict the existence of a high-mass tail in the BH mass function at   z = 0  . Consistency between the optical and X-ray census of the local BH mass function requires an upper limit to BH masses. This consistent picture also predicts that the slope of the   M _bh–σ  relation will evolve with redshift at the high-mass end. Models of self-regulation that explain the co-evolution of the stellar component and nuclear BHs naturally provide such an upper limit. The combination of multiwavelength constraints predicts the existence of UMBHs and simultaneously provides an upper limit to their masses. The typical hosts for these local UMBHs are likely the bright, central cluster galaxies in the nearby Universe.     

Low-frequency gravitational waves from cosmological compact binaries

We consider gravitational waves emitted by various populations of compact binaries at cosmological distances. We use population synthesis models to characterize the properties of double neutron stars, double black holes and double white dwarf binaries, and white dwarf–neutron star, white dwarf–black hole and black hole–neutron star systems.
We use the observationally determined cosmic star formation history to reconstruct the redshift distribution of these sources and their merging rate evolution.
The gravitational signals emitted by each source during its early spiralling in phase add randomly to produce a stochastic background in the low-frequency band with spectral strain amplitude between ~10⁻¹⁸ and ~5×10⁻¹⁷ Hz^−1/2 at frequencies in the interval ~5×10⁻⁶–5×10⁻⁵ Hz.
The overall signal, which at frequencies above 10⁻⁴ Hz is largely dominated by double white dwarf systems, might be detectable with LISA in the frequency range 1–10 mHz and acts like a confusion-limited noise component, which might limit the LISA sensitivity at frequencies above 1 mHz.     

Heating of the intergalactic medium by primordial miniquasars

A simple analytical model is used to calculate the X-ray heating of the intergalactic medium (IGM) for a range of black hole masses. This process is efficient enough to decouple the spin temperature of the IGM from the cosmic microwave background (CMB) temperature and produce a differential brightness temperature of the order of ∼ 5–20 mK out to distances as large as a few comoving Mpc, depending on the redshift, black hole mass and lifetime. We explore the influence of two types of black holes, those with and without ionizing ultraviolet radiation. The results of the simple analytical model are compared to those of a full spherically symmetric radiative transfer code. Two simple scenarios are proposed for the formation and evolution of black hole mass density in the Universe. The first considers an intermediate mass black hole that form as an end-product of pop III stars, whereas the second considers supermassive black holes that form directly through the collapse of massive haloes with low spin parameter. These scenarios are shown not to violate any of the observational constraints, yet produce enough X-ray photons to decouple the spin temperature from that of the CMB. This is an important issue for future high-redshift 21-cm observations.     

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.     

Acceleration of black hole universe

Recently, Zhang slightly modified the standard big bang theory and developed a new cosmological model called black hole universe, which is consistent with Mach’s principle, governed by Einstein’s general theory of relativity, and able to explain all observations of the universe. Previous studies accounted for the origin, structure, evolution, expansion, and cosmic microwave background radiation of the black hole universe, which grew from a star-like black hole with several solar masses through a supermassive black hole with billions of solar masses to the present state with hundred billion-trillions of solar masses by accreting ambient matter and merging with other black holes. This paper investigates acceleration of the black hole universe and provides an alternative explanation for the redshift and luminosity distance measurements of type Ia supernovae. The results indicate that the black hole universe accelerates its expansion when it accretes the ambient matter in an increasing rate. In other words, i.e., when the second-order derivative of the mass of the black hole universe with respect to the time is positive $\ddot{M}(t) > 0$ . For a constant deceleration parameter $q = -M(t) \ddot{M}(t)/\dot{M}(t) \sim-0.6$ , we can perfectly explain the type Ia supernova measurements with the reduced chi-square to be very close to unity, χ _red～1.0012. The expansion and acceleration of black hole universe are driven by external energy.     

The imprint of massive black hole formation models on the LISA data stream

The formation, merging and accretion history of massive black holes (MBHs) along the hierarchical build-up of cosmic structures leaves a unique imprint on the background of gravitational waves (GWs) at mHz frequencies. We study here, by means of dedicated simulations of black hole build-up, the possibility of constraining different models of black hole cosmic evolution using future GW space-borne missions, such as LISA . We consider two main scenarios for black hole formation, namely, one where seeds are light (  ≃10² M_⊙  , remnant of Population III stars) and one where seeds are heavy (  ≳10⁴ M_⊙  , direct collapse). In all the models we have investigated, MBH binary coalescences do not produce a stochastic GW background, but rather, a set of individual resolved events. Detection of several hundreds merging events in a 3-yr LISA mission will be the sign of a heavy seed scenario with efficient formation of black hole seeds in a large fraction of high-redshift haloes. At the other extreme, a low event rate, about a few tens in 3 yr, is peculiar of scenarios where either the seeds are light, and many coalescences do not fall into the LISA band, or seeds are massive, but rare. In this case a decisive diagnostic is provided by the shape of the mass distribution of detected events. Light binaries  ( m < 10⁴ M_⊙)  are predicted in a fairly large number in Population III remnant models, but are totally absent in direct collapse models. Finally, a further, helpful diagnostic of black hole formation models lies in the distribution of the mass ratios in binary coalescences. While heavy seed models predict that most of the detected events involve equal-mass binaries, in the case of light seeds, mass ratios are equally distributed in the range 0.1–1.     

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.     

High-frequency X-ray variability as a mass estimator of stellar and supermassive black holes

There is increasing evidence that supermassive black holes in active galactic nuclei (AGN) are scaled-up versions of Galactic black holes. We show that the amplitude of high-frequency X-ray variability in the hard spectral state is inversely proportional to the black hole mass over eight orders of magnitude. We have analysed all available hard-state data from RXTE of seven Galactic black holes. Their power density spectra change dramatically from observation to observation, except for the high-frequency (≳10 Hz) tail, which seems to have a universal shape, roughly represented by a power law of index −2. The amplitude of the tail,   C _M   (extrapolated to 1 Hz), remains approximately constant for a given source, regardless of the luminosity, unlike the break or quasi-periodic oscillation frequencies, which are usually strongly correlated with luminosity. Comparison with a moderate-luminosity sample of AGN shows that the amplitude of the tail is a simple function of black hole mass,   C _M = C / M   , where   C ≈ 1.25 M_⊙ Hz⁻¹  . This makes   C _M   a robust estimator of the black hole mass which is easy to apply to low- to moderate-luminosity supermassive black holes. The high-frequency tail with its universal shape is an invariant feature of a black hole and, possibly, an imprint of the last stable orbit.     

Joint cosmological formation of QSOs and bulge-dominated galaxies

Older and more recent pieces of observational evidence suggest a strong connection between QSOs and galaxies; in particular, the recently discovered correlation between black hole and galactic bulge masses suggests that QSO activity is directly connected to the formation of galactic bulges. The cosmological problem of QSO formation is analysed in the framework of an analytical model for galaxy formation; for the first time a joint comparison with galaxy and QSO observables is performed. In this model it is assumed that the same physical variable that determines galaxy morphology is able to modulate the mass of the black hole responsible for QSO activity. Both halo spin and the occurrence of a major merger are considered as candidates for this role. The predictions of the model are compared with available data for the type-dependent galaxy mass functions, the star formation history of elliptical galaxies, the QSO luminosity function and its evolution (including the obscured objects contributing to the hard-X-ray background), the mass function of dormant black holes and the distribution of black hole-to-bulge mass ratios. A good agreement with observations is obtained if the halo spin modulates the efficiency of black hole formation, and if the galactic haloes at z =0 have shone in an inverted order with respect to the hierarchical one (i.e., stars and black holes in bigger galactic haloes have formed before those in smaller ones). This inversion of hierarchical order for galaxy formation, which reconciles galaxy formation with QSO evolution, is consistent with many pieces of observational evidence.     

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 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.     

The runaway instability of thick discs around black holes I. The constant angular momentum case

We present results from a numerical study of the runaway instability of thick discs around black holes. This instability is an important issue for most models of cosmic gamma-ray bursts, where the central engine responsible for the initial energy release is such a system consisting of a thick disc surrounding a black hole. We have carried out a comprehensive number of time-dependent simulations aimed at exploring the appearance of the instability. Our study has been performed using a fully relativistic hydrodynamics code. The general relativistic hydrodynamic equations are formulated as a hyperbolic flux-conservative system and solved using a suitable Godunov-type scheme. We build a series of constant angular momentum discs around a Schwarzschild black hole. Furthermore, the self-gravity of the disc is neglected and the evolution of the central black hole is assumed to be that of a sequence of exact Schwarzschild black holes of varying mass. The black hole mass increase is thus determined by the mass accretion rate across the event horizon. In agreement with previous studies based on stationary models, we find that by allowing the mass of the black hole to grow the disc becomes unstable. Our hydrodynamical simulations show that for all disc-to-hole mass ratios considered (between 1 and 0.05), the runaway instability appears very fast on a dynamical time-scale of a few orbital periods, typically a few 10 ms and never exceeding 1 s for our particular choice of the mass of the black hole (2.5 M) and a large range of mass fluxes  ( m 10^-3 M s^-1)  . The implications of our results in the context of gamma-ray bursts are briefly discussed.     

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 growth of supermassive black holes in pseudo-bulges, classical bulges and elliptical galaxies

Using results from structural analysis of a sample of nearly 1000 local galaxies from the Sloan Digital Sky Survey, we estimate how the mass in central black holes is distributed amongst elliptical galaxies, classical bulges and pseudo-bulges, and investigate the relation between their stellar masses and central stellar velocity dispersion σ. Assuming a single relation between elliptical galaxy/bulge mass, M _Bulge, and central black hole mass, M _BH, we find that  55⁺⁸₋₄  per cent of the mass in black holes in the local universe is in the centres of elliptical galaxies,  41⁺⁴₋₂  per cent in classical bulges and  4^+0.9_−0.4  per cent in pseudo-bulges. We find that ellipticals, classical bulges and pseudo-bulges follow different relations between their stellar masses and σ, and the most significant offset occurs for pseudo-bulges in barred galaxies. This structural dissimilarity leads to discrepant black hole masses if single   M _BH– M _Bulge  and   M _BH–σ  relations are used. Adopting relations from the literature, we find that the   M _BH–σ  relation yields an estimate of the total mass density in black holes that is roughly 55 per cent larger than if the   M _BH– M _Bulge  relation is used.     

Fe Kα line profile in AGNs: average shape and redshift dependence

We collect a sample of 381(271 type Ⅰ)active galactic nuclei(AGNs)with XMM-Newton observations for an analysis of the dependence of Fe Kα profiles on redshifts to test the potential cosmological evolution of spins of supermassive black holes(SMBHs).The sample spans a redshift range of z = 0.0008-4.76,which allows us to rebin the sample into 7 redshift groups.Phenomenological analysis of the Fe profile suggested that the line width(σ)of the narrow or broad Fe line does not show significant changes in redshift range z < 0.3.Using a physical model,we significantly detect a narrow Fe Kα line at 6.4 keV with an average equivalent width(EW)= 160 eV except for the two largest redshift bins.The EW of the Fe line does not show significant changes.We also find a relativistic broad line in three redshift bins(z < 0.116,0.056 < z < 0.12 and 0.12 < z < 0.3)with an average EW = 522eV.