Radiatively driven mass accretion on to galactic nuclei by circumnuclear starbursts

We examine the physical processes of radiatively driven mass accretion on to galactic nuclei, owing to intensive radiation from circumnuclear starbursts. The radiation from a starburst not only causes the inner gas disc to contract via radition flux force, but also extracts angular momentum owing to relativistic radiation drag, thereby inducing an avalanche of the surface layer of the disc. To analyse such a mechanism, the radiation–hydrodynamical equations are solved, including the effects of the radiation drag force as well as the radiation flux force. As a result, it is found that the mass accretion rate owing to the radiative avalanche is given by M ˙ ( r )= η ( L _*/ c ²)( r / R )² (Δ R / R )(1 −  e ^−τ) at radius r , where the efficiency η ranges from 0.2 up to 1, L _* and R are respectively the bolometric luminosity and the radius of the starburst ring, Δ R is the extent of the emission regions, and τ is the face-on optical depth of the disc. In an optically thick regime, the rate depends upon neither the optical depth nor the surface mass density distribution of the disc. The present radiatively driven mass accretion may provide a physical mechanism which enables mass accretion from 100-pc scales down to ∼ parsec scales, and it may eventually be linked to advection-dominated viscous accretion on to a massive black hole. The radiation–hydrodynamical and self-gravitational instabilities of the disc are briefly discussed. In particular, the radiative acceleration possibly builds up a dusty wall, which 'shades' the nucleus in edge-on views. This provides another version of the model for the formation of an obscuring torus.     

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 detectability of distant Compton-thick obscured quasars

Chandra and XMM–Newton have resolved the     X-ray background (XRB) into point sources. Many of the fainter sources are obscured active galactic nuclei (AGN) with column densities in the range of     , some of which have quasar-like luminosities. According to obscuration models, the XRB above 8 keV is dominated by emission from Compton-thick AGN, with column densities exceeding     . Here, we consider whether Compton-thick quasars are detectable by Chandra and XMM–Newton by their direct (i.e. not scattered) X-ray emission. Detectability is optimized if the objects individually have a high luminosity and high redshift, so that the direct emission has a significant flux in the observed band. Using a simple galaxy formation model incorporating accreting black holes, in which quasars build most of their mass in a Compton-thick manner before expelling the obscuring matter, we predict that moderately deep 100-ks Chandra and XMM–Newton exposures may contain a handful of detectable Compton-thick quasars. Deep Ms or more Chandra images should contain     distant, optically faint, Compton-thick sources. In passing we show that radiation pressure can be as effective in expelling the obscuring gas as quasars winds, and yields a black hole mass proportional to the velocity dispersion of the host bulge to the fourth power.     

Feast and Famine: regulation of black hole growth in low-redshift galaxies

We analyse the observed distribution of Eddington ratios  ( L / L _Edd)  as a function of supermassive black hole mass for a large sample of nearby galaxies drawn from the Sloan Digital Sky Survey. We demonstrate that there are two distinct regimes of black hole growth in nearby galaxies. The first is associated with galaxies with significant star formation [   M _*/star formation rate (SFR) ∼  a Hubble time] in their central kiloparsec regions, and is characterized by a broad lognormal distribution of accretion rates peaked at a few per cent of the Eddington limit. In this regime, the Eddington ratio distribution is independent of the mass of the black hole and shows little dependence on the central stellar population of the galaxy. The second regime is associated with galaxies with old central stellar populations (   M _*/SFR ≫  a Hubble time), and is characterized by a power-law distribution function of Eddington ratios. In this regime, the time-averaged mass accretion rate on to black holes is proportional to the mass of stars in the galaxy bulge, with a constant of proportionality that depends on the mean stellar age of the stars. This result is once again independent of black hole mass. We show that both the slope of the power law and the decrease in the accretion rate on to black holes in old galaxies are consistent with population synthesis model predictions of the decline in stellar mass loss rates as a function of mean stellar age. Our results lead to a very simple picture of black hole growth in the local Universe. If the supply of cold gas in a galaxy bulge is plentiful, the black hole regulates its own growth at a rate that does not further depend on the properties of the interstellar medium. Once the gas runs out, black hole growth is regulated by the rate at which evolved stars lose their mass.     

Black hole demographics from the relation

We analyse a sample of 32 galaxies for which a dynamical estimate of the mass of the hot stellar component, M _bulge, is available. For each of these galaxies, we calculate the mass of the central black hole, M _•, using the tight empirical correlation between M _• and bulge stellar velocity dispersion. The frequency function     is reasonably well described as a Gaussian with     and standard deviation ∼0.45; the implied mean ratio of black hole mass to bulge mass is a factor of ∼5 smaller than generally quoted in the literature. We present marginal evidence for a lower, average black hole mass fraction in more massive galaxies. The total mass density in black holes in the local Universe is estimated to be ∼     consistent with that inferred from high-redshift     active galactic nuclei.     

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.     

Supermassive black holes, star formation and downsizing of elliptical galaxies

The overabundance of Mg relative to Fe, observed in the nuclei of bright ellipticals, and its increase with galactic mass, poses a serious problem for all current models of galaxy formation. Here, we improve on the one-zone chemical evolution models for elliptical galaxies by taking into account positive feedback produced in the early stages of supermassive central black hole growth. We can account for both the observed correlation and the scatter if the observed anti-hierarchical behaviour of the AGN population couples to galaxy assembly and results in an enhancement of the star formation efficiency which is proportional to galactic mass. At low and intermediate galactic masses, however, a slower mode for star formation suffices to account for the observational properties.     

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.     

How galaxies lose their angular momentum

The processes are investigated by which gas loses its angular momentum during the protogalactic collapse phase, leading to disc galaxies that are too compact with respect to the observations. High-resolution N -body/SPH simulations in a cosmological context are presented including cold gas and dark matter (DM). A halo with quiet merging activity since redshift   z ∼ 3.8  and with a high-spin parameter is analysed that should be an ideal candidate for the formation of an extended galactic disc. We show that the gas and the DM have similar specific angular momenta until a merger event occurs at   z ∼ 2  with a mass ratio of 5:1. All the gas involved in the merger loses a substantial fraction of its specific angular momentum due to tidal torques and dynamical friction processes falls quickly into the centre. In contrast, gas infall through small subclumps or accretion does not lead to catastrophic angular momentum loss. In fact, a new extended disc begins to form from gas that was not involved in the 5:1 merger event and that falls in subsequently. We argue that the angular momentum problem of disc galaxy formation is a merger problem: in cold dark matter cosmology substantial mergers with mass ratios of 1:1 to 6:1 are expected to occur in almost all galaxies. We suggest that energetic feedback processes could in principle solve this problem, however only if the heating occurs at the time or shortly before the last substantial merger event. Good candidates for such a coordinated feedback would be a merger-triggered starburst or central black hole heating. If a large fraction of the low angular momentum gas would be ejected, late-type galaxies could form with a dominant extended disc component, resulting from late infall, a small bulge-to-disc ratio and a low baryon fraction, in agreement with observations.     

The host galaxy of a narrow-line Seyfert 1 galaxy, RE J1034+396, with X-ray quasi-periodic oscillations

Using simple stellar population synthesis, we model the bulge stellar contribution in the optical spectrum of a narrow-line Seyfert 1 galaxy, RE J1034+396. We find that its bulge stellar velocity dispersion is  67.7 ± 8 km s⁻¹  . The supermassive black hole (SMBH) mass is about  (1–4) × 10⁶ M_⊙  if it follows the well-known   M _BH–σ_*  relation found in quiescent galaxies. We also derive the SMBH mass from the Hβ second moment, which is consistent with that from its bulge stellar velocity dispersion. The SMBH mass of (1–4)  × 10⁶ M_⊙  implies that the X-ray quasi-periodic oscillation (QPO) of RE J1034+396 can be scaled to a high-frequency QPO at 27–108 Hz found in Galactic black hole binaries with a  10-M_⊙  black hole. With the mass distribution in different age stellar populations, we find that the mean specific star formation rate (SSFR) over the past 0.1 Gyr is  0.0163 ± 0.0011  Gyr⁻¹, the stellar mass in the logarithm is  10.155 ± 0.06  in units of solar mass and the current star formation rate is  0.23 ± 0.016 M_⊙ yr⁻¹  . For RE J1034+396, there is no relation between the Eddington ratio and the SSFR as suggested by Chen et al., despite a larger scatter in their relation. We also suggest that about 7.0 per cent of the total Hα luminosity and 50 per cent of the total [O  ii ] luminosity come from the star formation process.     

Resonant tidal disruption in galactic nuclei

It has recently been shown by Rauch 38 Tremaine that the rate of angular momentum relaxation in nearly Keplerian star clusters is greatly increased by a process termed 'resonant relaxation'; it was also argued, via a series of scaling arguments, that tidal disruption of stars in galactic nuclei containing massive black holes could be noticeably enhanced by this process. We describe here the results of numerical simulations of resonant tidal disruption which quantitatively test the predictions made by Rauch 38 Tremaine. The simulation method is based on an N -body routine incorporating cloning of stars near the loss cone and a semirelativistic symplectic integration scheme. Normalized disruption rates for resonant and non-resonant nuclei are derived at orbital energies both above and below the critical energy, and the corresponding angular momentum distribution functions are found. The black hole mass above which resonant tidal disruption is quenched by relativistic precession is determined. We also briefly describe the discovery of chaos in the Wisdom–Holman symplectic integrator applied to highly eccentric orbits and propose a modified integration scheme that remains robust under these conditions. We find that resonant disruption rates exceed their non-resonant counterparts by an amount consistent with the predictions; in particular, we estimate the net tidal disruption rate for a fully resonant cluster to be about twice that of its non-resonant counterpart. No significant enhancement in rates is observed outside the critical radius. Relativistic quenching of the effect is found to occur for hole masses M  >  M _Q  = (8 ± 3) × 10⁷  M . The numerical results combined with the observed properties of galactic nuclei indicate that for most galaxies the resonant enhancement to tidal disruption rates will be very small.     

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 effect of a central supermassive black hole on gas fuelling

When a supermassive black hole exists in the centre of a galaxy, an additional inner Lindblad resonance (ILR) exists inside the usual ILRs. We study gas dynamics in a weakly barred potential with a central supermassive black hole by using 2D numerical simulations, and we investigate the effect of the additional ILR on the fuelling of gas into nuclear starburst regions or active galactic nuclei (AGNs). Our numerical results show that strong trailing spiral shocks are formed at the resonance region, and that the gas in the shock region is rapidly fuelled into a central region and makes a nuclear gas ring. As a result, a large amount of gas is concentrated in the nuclear region beyond the ILR in a dynamical time-scale.     

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.     

Reprocessed emission line profiles from dense clouds in geometrically thick accretion engines

The central engines of active galactic nuclei (AGN) contain cold, dense material as well as hot X-ray-emitting gas. The standard paradigm for the engine geometry is a cold thin disc sandwiched between hot X-ray coronae. Strong support for this geometry in Seyferts comes from the study of fluorescent iron line profiles, although the evidence is not ubiquitously airtight. The thin disc model of line profiles in AGN and in X-ray binaries should still be benchmarked against other plausible possibilities. One proposed alternative is an engine consisting of dense clouds embedded in an optically thin, geometrically thick X-ray-emitting engine. This model is also motivated by studies of geometrically thick engines such as advection-dominated accretion flows (ADAFs). Here we compute the reprocessed iron line profiles from dense clouds embedded in geometrically thick, optically thin X-ray-emitting discs near a Schwarzschild black hole. We consider a range of cloud distributions and disc solutions, including ADAFs, pure radial infall and bipolar outflows. We find that such models can reproduce line profiles similar to those from geometrically thin, optically thick discs and might help alleviate some of the problems encountered from the latter. Thus, independent of thin discs, thick disc engines can also exhibit iron line profiles if embedded dense clouds can survive long enough to reprocess radiation.     

The changing interstellar medium of massive elliptical galaxies and cosmic evolution of radio galaxies and quasars

The recently discovered apparent dramatic expansion in the effective radii of massive elliptical galaxies from   z ≃ 2  to ≃0.1 has been interpreted in terms of either galaxy mergers or the rapid loss of cold gas due to active galactic nuclei (AGN) feedback. In examining the latter case, we have quantified the extent of the expansion, which is uncertain observationally, in terms of the star formation parameters and time of the expulsion of the cold gas. In either case, the large global decrease in stellar density should translate into a major drop in the interstellar medium density and pressure with cosmic epoch. These cosmological changes are expected to have a major influence on the gas accretion mode, which will shift from 'cold' thin disc accretion at high redshifts towards 'hot' Bondi fed Advection Dominated Accretion Flow (ADAF) accretion at low redshifts. The decline of angular momentum inflow would then lead to a spin down of the black hole, for which we have calculated more precise time-scales; a value of about 0.2 Gyr is typical for a  10⁹ M_⊙  central black hole. These results have implications for the different cosmological evolutionary patterns found for the luminosity functions of powerful and weak radio galaxies.     

The tidal disruption rate in dense galactic cusps containing a supermassive binary black hole

We consider the problem of tidal disruption of stars in the centre of a galaxy containing a supermassive binary black hole with unequal masses. We assume that over the separation distance between the black holes, the gravitational potential is dominated by the more massive primary black hole. Also, we assume that the number density of stars is concentric with the primary black hole and has a power-law cusp. We show that the bulk of stars with a small angular-momentum component normal to the black hole binary orbit can reach a small value of total angular momentum through secular evolution in the gravitational field of the binary, and hence they can be tidally disrupted by the larger black hole. This effect is analogous to the so-called Kozai effect well known in celestial mechanics. We develop an analytical theory for the secular evolution of the stellar orbits and calculate the rate of tidal disruption. We compare our analytical theory with a simple numerical model and find very good agreement.
Our results show that for a primary black hole mass of  ∼10⁶–10⁷ M_⊙  , the black hole mass-ratio   q > 10⁻²  , cusp size ∼1 pc, the tidal disruption rate can be as large as  ∼10⁻²–1 M_⊙ yr⁻¹  . This is at least 10²–10⁴ times larger than estimated for the case of a single supermassive black hole. The duration of the phase of enhanced tidal disruption is determined by the dynamical-friction time-scale, and it is rather short: ∼10⁵ yr. The dependence of the tidal disruption rate on the mass ratio, and on the size of the cusp, is also discussed.     

The origin of the density distribution of disc galaxies: a new problem for the standard model of disc formation

We present new models for the formation of disc galaxies that improve upon previous models by following the detailed accretion and cooling of the baryonic mass, and by using realistic distributions of specific angular momentum. Under the assumption of detailed angular momentum conservation, the discs that form have density distributions that are more centrally concentrated than an exponential. We examine the influence of star formation, bulge formation, and feedback on the outcome of the surface brightness distributions of the stars. Low angular momentum haloes yield disc galaxies with a significant bulge component and with a stellar disc that is close to exponential, in good agreement with observations. High angular momentum haloes, on the other hand, produce stellar discs that are much more concentrated than an exponential, in clear conflict with observations. At large radii, the models reveal distinct truncation radii in both the stars and the cold gas. The stellar truncation radii result from our implementation of star formation threshold densities, and are in excellent agreement with observations. The truncation radii in the density distribution of the cold gas reflect the maximum specific angular momentum of the gas that has cooled. We find that these truncation radii occur at H  i surface densities of roughly 1 M_⊙ pc⁻², in conflict with observations. We examine various modifications to our models, including feedback, viscosity, and dark matter haloes with constant-density cores, but show that the models consistently fail to produce bulge less discs with exponential surface brightness profiles. This signals a new problem for the standard model of disc formation: if the baryonic component of the protogalaxies out of which disc galaxies form has the same angular momentum distribution as the dark matter, discs are too compact.     

Cooling flows as a calorimeter of active galactic nucleus mechanical power

The assumption that radiative cooling of gas in the centres of galaxy clusters is approximately balanced by energy input from a central supermassive black hole implies that the observed X-ray luminosity of the cooling flow region sets a lower limit on active galactic nucleus (AGN) mechanical power. The conversion efficiency of the mechanical power of the AGN into gas heating is uncertain, but we argue that it can be high even in the absence of strong shocks. These arguments inevitably lead to the conclusion that the time-averaged mechanical power of AGNs in cooling flows is much higher than the bolometric luminosity of these objects observed currently.
The energy balance between cooling losses and AGN mechanical power requires some feedback mechanism. We consider a toy model in which the accretion rate on to a black hole is set by the classic Bondi formula. Application of this model to the best studied case of M87 suggests that accretion proceeds at approximately the Bondi rate down to a few gravitational radii with most of the power (at the level of a few per cent of the rest mass) being carried away by an outflow.