Formation of supermassive black holes

Evidence shows that massive black holes reside in most local galaxies. Studies have also established a number of relations between the MBH mass and properties of the host galaxy such as bulge mass and velocity dispersion. These results suggest that central MBHs, while much less massive than the host (~0.1%), are linked to the evolution of galactic structure. In hierarchical cosmologies, a single big galaxy today can be traced back to the stage when it was split up in hundreds of smaller components. Did MBH seeds form with the same efficiency in small proto-galaxies, or did their formation had to await the buildup of substantial galaxies with deeper potential wells? I briefly review here some of the physical processes that are conducive to the evolution of the massive black hole population. I will discuss black hole formation processes for ‘seed’ black holes that are likely to place at early cosmic epochs, and possible observational tests of these scenarios.     

Measuring the accretion rate and kinetic luminosity functions of supermassive black holes

We derive accretion rate functions (ARFs) and kinetic luminosity functions (KLFs) for jet-launching supermassive black holes. The accretion rate as well as the kinetic power of an active galaxy is estimated from the radio emission of the jet. For compact low-power jets, we use the core radio emission while the jet power of high-power radio-loud quasars is estimated using the extended low-frequency emission to avoid beaming effects. We find that at low luminosities the ARF derived from the radio emission is in agreement with the measured bolometric luminosity function (BLF) of active galactic nucleus (AGN), i.e. all low-luminosity AGN launch strong jets. We present a simple model, inspired by the analogy between X-ray binaries (XRBs) and AGN, that can reproduce both the measured ARF of jet-emitting sources as well as the BLF. The model suggests that the break in power-law slope of the BLF is due to the inefficient accretion of strongly sub-Eddington sources. As our accretion measure is based on the jet power it also allows us to calculate the KLF and therefore the total kinetic power injected by jets into the ambient medium. We compare this with the kinetic power output from supernova remnants (SNRs) and XRBs, and determine its cosmological evolution.     

Mass-losing accretion discs around supermassive black holes

We study the effects of outflow/wind on the gravitational stability of accretion discs around supermassive black holes using a set of analytical steady-state solutions. Mass-loss rate by the outflow from the disc is assumed to be a power-law of the radial distance and the amount of the energy and the angular momentum which are carried away by the wind are parameterized phenomenologically. We show that the mass of the first clumps at the self-gravitating radius linearly decreases with the total mass-loss rate of the outflow. Except for the case of small viscosity and high accretion rate, generally, the self-gravitating radius increases as the amount of mass-loss by the outflow increases. Our solutions show that as more angular momentum is lost by the outflow, then reduction to the mass of the first clumps is more significant.     

Probing the nuclear activity with supermassive black holes

We report three new or updated techniques for probing the parameters of active galaxies based on the masses of their central black holes M_BH). First, we derived a near-IR analog of the bulge luminosity versus M_BH relationship. The low scatter makes it a promising new tool to study the black hole demographics. Next, we present relations between M_BH and the10 μm and 2-10 keV nuclear luminosity. They may help to study the M_BH evolution over wide redshift ranges. Finally, we measured M_BH in quasars from z ∼ 3.4 to z ∼ 0.3 to search directly for M_BH growth. Surprisingly, we found no evidence for growth implying that the majority of quasar host galaxies have undergone their last major merger at z ≥ 3. This revised version was published online in August 2006 with corrections to the Cover Date.     

X-ray reverberation around accreting black holes

Luminous accreting stellar mass and supermassive black holes produce power–law continuum X-ray emission from a compact central corona. Reverberation time lags occur due to light travel time delays between changes in the direct coronal emission and corresponding variations in its reflection from the accretion flow. Reverberation is detectable using light curves made in different X-ray energy bands, since the direct and reflected components have different spectral shapes. Larger, lower frequency, lags are also seen and are identified with propagation of fluctuations through the accretion flow and associated corona. We review the evidence for X-ray reverberation in active galactic nuclei and black hole X-ray binaries, showing how it can be best measured and how it may be modelled. The timescales and energy dependence of the high-frequency reverberation lags show that much of the signal is originating from very close to the black hole in some objects, within a few gravitational radii of the event horizon. We consider how these signals can be studied in the future to carry out X-ray reverberation mapping of the regions closest to black holes.     

Self-gravitating warped discs around supermassive black holes

We consider warped equilibrium configurations for stellar and gaseous discs in the Keplerian force field of a supermassive black hole, assuming that the self-gravity of the disc provides the only acting torques. Modelling the disc as a collection of concentric circular rings and computing the torques in the non-linear regime, we show that stable, strongly warped precessing equilibria are possible. These solutions exist for a wide range of disc-to-black-hole mass ratios   M _d/ M _bh  , can span large warp angles of up to  ±∼120°  , have inner and outer boundaries, and extend over a radial range of a factor of typically two to four. These equilibrium configurations obey a scaling relation such that in good approximation     where     is the (retrograde) precession frequency and Ω is a characteristic orbital frequency in the disc. Stability was determined using linear perturbation theory and, in a few cases, confirmed by numerical integration of the equations of motion. Most of the precessing equilibria are found to be stable, but some are unstable. The main result of this study is that highly warped discs near black holes can persist for long times without any persistent forcing other than by their self-gravity. The possible relevance of this to galactic nuclei is briefly discussed.     

Migration of supermassive black holes in galactic nuclei

The migration of central black holes in galactic nuclei through their encounters with galactic globular clusters is studied. The black hole moves in the field of the galactic bulge with a fixed potential. The dependences of the black-hole drift amplitude on orbital parameters of the globular cluster, its mass, and bulge parameters have been found. The drift amplitude of the central black hole can reach several parsecs in our Galaxy and several tens of parsecs in early-type (Sa) and late-type (Sc) spiral galaxies.     

A simple way to classify supermassive black holes

We propose a classification of supermassive black holes (SMBHs) based on their efficiency in the conversion of infalling mass in emitted radiation. We use a theoretical model that assumes a conservation of angular momentum between the gas falling inside the hole and the photons emitted outwards, and suggests the existence of the scaling relation M_•‐R_eσ³, where M_• is the mass of the central SMBH, whereas R_e and σ are the effective radius and velocity dispersion of the host galaxies (bulges), respectively. We apply our model on a data set of 57 galaxies of different morphological types and with M_• measurements, obtained through the analysis of Spitzer /IRAC 3.6‐µ m images. In order to find the best fit of the corresponding scaling law, we use the FITEXY routine to perform a least‐squares regression of M_• on R_eσ³ for the considered sample of galaxies. Our analysis shows that the relation is tight and our theoretical model allows to easily estimate the efficiency of mass conversion into radiation of the central SMBHs. Finally we propose a new appealing way to classify the SMBHs in terms of this parameter. (© 2014 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Ultra-high energy neutrino fluxes from supermassive AGN black holes

We compute the ultra-high energy (UHE) neutrino fluxes from plausible accreting supermassive black holes closely linking to the 377 active galactic nuclei (AGNs). They have well-determined black hole masses collected from the literature. The neutrinos are produced via simple or modified URCA processes, even after the neutrino trapping, in superdense proto-matter medium. The resulting fluxes are ranging from: (1) (quark reactions)— $J^{q}_{\nu\varepsilon}/(\varepsilon_{d}\ \mathrm{erg}\,\mathrm{cm}^{-2}\,\mathrm{s}^{-1}\,\mathrm{sr}^{-1})\simeq8.29\times 10^{-16}$ to 3.18×10^?4, with the average $\overline{J}^{q}_{\nu\varepsilon}\simeq5.53\times 10^{-10}\varepsilon_{d}\ \mathrm{erg}\,\mathrm{cm}^{-2}\,\mathrm{s}^{-1}\,\mathrm{sr}^{-1}$ , where ε _d ～10^?12 is the opening parameter; (2) (pionic reactions)— $J^{\pi}_{\nu\varepsilon} \simeq0.112J^{q}_{\nu\varepsilon}$ , with the average $J^{\pi}_{\nu\varepsilon} \simeq3.66\times 10^{-11}\varepsilon_{d}\ \mathrm{erg}\,\mathrm{cm}^{-2}\,\mathrm{s}^{-1}\,\mathrm{sr}^{-1}$ ; and (3) (modified URCA processes)— $J^{URCA}_{\nu\varepsilon}\simeq7.39\times10^{-11} J^{q}_{\nu\varepsilon}$ , with the average $\overline{J}^{URCA}_{\nu\varepsilon} \simeq2.41\times10^{-20} \varepsilon_{d}\ \mathrm{erg}\,\mathrm{cm}^{-2}\,\mathrm{s}^{-1}\,\mathrm{sr}^{-1}$ . We conclude that the AGNs are favored as promising pure neutrino sources, because the computed neutrino fluxes are highly beamed along the plane of accretion disk, peaked at high energies and collimated in smaller opening angle θ～ε _d .     

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.     

Radial self-gravity of accretion disc around supermassive black holes

This paper studies the properties of self-gravity of accretion discs around supermassive black holes. With integration on the thin disc configuration, this paper has calculated the radial and vertical components of self-gravity of accretion discs. The discussion mainly concentrates on the radial component, and the results are briefly as follows: for accretion discs around supermassive black holes (M10 ⁸–10¹⁰ M). At the distance where (R/R _g)10⁵–10⁴, the radial component of self-gravity dominates over the central one where the dynamical structure of the accretion discs completely differs from that of Keplerian disc. A turbulence driven by radial self-gravity instability as a kind of energy source is proposed. This paper has two criteria for the comparison of magnitude between the self-gravity of accretion discs and the gravity of the central black hole, from which an analytic estimation for the outer radius of the accretion discs has been derived. The results of this paper may be used to explain the accretion discs of quasars and AGNs.This research was supported by the National Natural Science Foundation of China.     

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.     

Gravitational instability of discs with dissipative coronae around supermassive black holes

We study the dynamical structure of a self-gravitating disc with coronae around a supermassive black hole. Assuming that the magnetorotational instability responsible for generating the turbulent stresses inside the disc is also the source for a magnetically dominated corona, a fraction of the power released when the disc matter accretes is transported to and dissipated in the corona. This has a major effect on the structure of the disc and its gravitational (in)stability according to our analytical and self-consistent solutions. We determine the radius where the disc crosses the inner radius of gravitational instability and forms the first stars. Not only the location of this radius which may extend to very large distances from the central black hole, but also the mass of the first stars highly depends on the input parameters, notably the viscosity coefficient, the mass of the central object and the accretion rate. For accretion discs around quasi-stellar objects (QSOs) and the Galactic Centre, we determine the self-gravitating radius and the mass of the first clumps. Comparing the cases with a corona and without a corona for typical discs around QSOs or the Galactic Centre, when the viscosity coefficient is around 0.3, we show that the self-gravitating radius decreases by a factor of approximately 2, but the mass of the fragments increases with more or less the same factor. The existence of a corona implies a more gravitationally unstable disc according to our results. The effect of a corona on the instability of the disc is more effective when the viscosity coefficient increases.