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.     

Effects of Kerr space–time on spectral features from X-ray illuminated accretion discs

We performed detailed calculations of the relativistic effects acting on both the reflection continuum and the iron line from accretion discs around rotating black holes. Fully relativistic transfer of both illuminating and reprocessed photons has been considered in Kerr space–time. We calculated overall spectra, line profiles and integral quantities, and present their dependences on the black hole angular momentum. We show that the observed EW of the lines is substantially enlarged when the black hole rotates rapidly and/or the source of illumination is near above the hole. Therefore, such calculations provide a way to distinguish between different models of the central source.     

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.     

Quasi-viscous accretion flow – I. Equilibrium conditions and asymptotic behaviour

In a novel approach to studying viscous accretion flows, viscosity has been introduced as a perturbative effect, involving a first-order correction in the α-viscosity parameter. This method reduces the problem of solving a second-order non-linear differential equation (Navier–Stokes equation) to that of an effective first-order equation. Viscosity breaks down the invariance of the equilibrium conditions for stationary inflow and outflow solutions, and distinguishes accretion from wind. Under a dynamical systems classification, the only feasible critical points of this 'quasi-viscous' flow are saddle points and spirals. On large spatial scales of the disc, where a linearized and radially propagating time-dependent perturbation is known to cause a secular instability, the velocity evolution equation of the quasi-viscous flow has been transformed to bear a formal closeness with Schrödinger's equation with a repulsive potential. Compatible with the transport of angular momentum to the outer regions of the disc, a viscosity-limited length-scale has been defined for the full spatial extent over which the accretion process would be viable.     

A comparative HST imaging study of the host galaxies of radio-quiet quasars, radio-loud quasars and radio galaxies – I

We present the first results from a major HST WFPC2 imaging study aimed at providing the first statistically meaningful comparison of the morphologies, luminosities, scalelengths and colours of the host galaxies of radio-quiet quasars, radio-loud quasars and radio galaxies. We describe the design of this study and present the images that have been obtained for the first half of our 33-source sample. We find that the hosts of all three classes of luminous AGN are massive elliptical galaxies, with scalelengths ≃10 kpc, and R − K colours consistent with mature stellar populations. Most importantly, this is first unambiguous evidence that, just like radio-loud quasars, essentially all radio-quiet quasars brighter than M _R =−24 reside in massive ellipticals. This result removes the possibility that radio 'loudness' is directly linked to host galaxy morphology, but is however in excellent accord with the black hole/spheroid mass correlation recently highlighted by Magorrian et al. We apply the relations given by Magorrian et al. to infer the expected Eddington luminosity of the putative black hole at the centre of each of the spheroidal host galaxies we have uncovered. Comparison with the actual nuclear R -band luminosities suggests that the black holes in most of these galaxies are radiating at a few per cent of the Eddington luminosity; the brightest host galaxies in our low- z sample are capable of hosting quasars with M _R ≃− 28, comparable to the most luminous quasars at z ≃3. Finally, we discuss our host-derived black hole masses in the context of the radio luminosity:black hole mass correlation recently uncovered for nearby galaxies by Franceschini et al., and consider the resulting implications for the physical origin of radio loudness.     

Shock oscillation model for quasi-periodic oscillations in stellar mass and supermassive black holes

We numerically examine centrifugally supported shock waves in 2D rotating accretion flows around a stellar mass  (10 M_⊙)  and a supermassive  (10⁶ M_⊙)  black holes over a wide range of input accretion rates of     . The resultant 2D shocks are unstable with time and the luminosities show quasi-periodic oscillations (QPOs) with modulations of a factor of 2–3 and with periods of a tenth of a second to several hours, depending on the black hole masses. The shock oscillation model may explain the intermediate frequency QPOs with 1–10 Hz observed in the stellar mass black hole candidates and also suggest the existence of QPOs with the period of hours in active galactic nuclei. When the accretion rate     is low, the luminosity increases in proportion to the accretion rate. However, when     greatly exceeds the Eddington critical rate     , the luminosity is insensitive to the accretion rate and is kept constantly around  ∼3 L _E  . On the other hand, the mass-outflow rate     increases in proportion to     and it amounts to about a few per cent of the input mass-flow rate.     

Hydrodynamical non-radiative accretion flows in two dimensions

Two-dimensional (axially symmetric) numerical hydrodynamical calculations of accretion flows that cannot cool through emission of radiation are presented. The calculations begin from an equilibrium configuration consisting of a thick torus with constant specific angular momentum. Accretion is induced by the addition of a small anomalous azimuthal shear stress which is characterized by a function ν . We study the flows generated as the amplitude and form of ν are varied. A spherical polar grid which spans more than two orders of magnitude in radius is used to resolve the flow over a wide range of spatial scales. We find that convection in the inner regions produces significant outward mass motions that carry away both the energy liberated by and a large fraction of the mass participating in the accretion flow. Although the instantaneous structure of the flow is complex and dominated by convective eddies, long-time averages of the dynamical variables show remarkable correspondence to certain steady-state solutions. The two-dimensional structure of the time-averaged flow is marginally stable to the Høiland criterion, indicating that convection is efficient. Near the equatorial plane, the radial profiles of the time-averaged variables are power laws with an index that depends on the radial scaling of the shear stress. A stress in which ν ∝ r ^1/2 recovers the widely studied self-similar solution corresponding to an ' α -disc'. We find that, regardless of the adiabatic index of the gas, or the form or magnitude of the shear stress, the mass inflow rate is a strongly increasing function of radius, and is everywhere nearly exactly balanced by mass outflow. The net mass accretion rate through the disc is only a fraction of the rate at which mass is supplied to the inflow at large radii, and is given by the local, viscous accretion rate associated with the flow properties near the central object.