The outer parts of standard steady-state accretion discs around quasi-stellar objects (QSOs) are prone to self-gravity, and they might be expected to fragment into stars rather than feed the central black hole. Possible solutions to this well-known problem are examined with an emphasis on general dynamic constraints. Irradiation by the QSO is insufficient for stability even if the outer disc is strongly warped. Marginal local gravitational instability enhances viscous transport but extends the stable regions only modestly. Compton cooling in the observed QSO radiation field rules out hot thick discs unless the local accretion rate is vastly super-Eddington. The formation of stars or stellar-mass black holes, and the release of energy in these objects by fusion or accretion, may help to stabilize the remaining gas in an otherwise standard disc. But at fixed mass accretion rate, the energy inputs required for stability increase with radius; beyond a parsec, they approach the total QSO luminosity and are probably unsustainable by stars. Magnetic torques from a wind or corona, and gravitational torques from bars or global spirals, may increase the accretion speed and reduce the density of the disc. But dynamical arguments suggest that the accretion speed is at most sonic, so that instability still sets in beyond about a parsec. Alternatively, the QSO could be fed by stellar collisions in a very dense stellar cluster, but the velocity dispersion would have to be much higher than observed in nearby galactic nuclei containing quiescent black holes. In view of these difficulties, we suggest that QSO discs do not extend beyond a thousand Schwarzschild radii or so. Then they must be frequently replenished with gas of small specific angular momentum.
The outer parts of steady, geometrically thin, optically thick, viscously-driven accretion discs around quasi-stellar objects (QSOs) are predicted to be self-gravitating (e.g. Shlosman & Begelman 1987). This is due to the high mass accretion rate required to feed a QSO and the relatively shallow gravitational potential at large radius. Under standard assumptions, the Toomre stability parameter Q falls below unity at r ≳ 10−2 pc ∼103RS (Section 2). It seems unlikely that a strongly self-gravitating disc can persist much longer than an orbital time.
There are, however, several reasons to suspect that QSO discs are larger than 10−2 pc. The mass of the disc within this radius is much less than that of the black hole. Hence, in order to increase the mass of the hole substantially, or equivalently, to fuel QSO activity for a Salpeter (1964) timeSanders et al. 1989). Nearby Seyferts and other active galaxies, thought to be low-luminosity analogues of QSOs, sometimes have resolvable nuclear discs on scales ∼102 pc. Most spectacularly, very long baseline interferometry (VLBI) observations of maser emission in NGC 4258 and 1068 indicate discs on parsec scales and have been used to measure the mass of their central black holes (Nakai, Inoue & Miyoshi 1993; Greenhill et al. 1995; Greenhill & Gwinn 1997). It is usually difficult to know from observation whether the outer disc is as dense as would be required to support the central luminosity on the assumption of a radially constant accretion rate.
At least on larger galactic scales, it is believed that strong self-gravity leads to star formation (Martin & Kennicutt 2001). If the same is true of QSO discs, there is a danger that most of the gas would form stars, leaving little to fuel the QSO.
Several theoretical attempts have been made to modify the standard α-disc model so as to extend gaseous QSO discs self-consistently into the self-gravitating regime. The options include: driving accretion with global bars or density waves so as to increase the radial velocity and lower the surface density needed to sustain a given mass accretion rate (Shlosman & Begelman 1989); heating or mechanically stirring the disc with embedded stars or black holes so as to raise the temperature and lower the density of the gas (Collin & Zahn 1999a,b); or allowing the disc to fragment into colliding gas clumps whose epicyclic motions stabilize it against self-gravity, at least on scales larger than the clumps themselves (Kumar 1999).
Our purpose is to re-examine self-gravity in QSO accretion discs with an emphasis on dynamical and energetic constraints. These constraints are most severe for massive and luminous systems, so our interest is in black-hole masses M ≳ 108 M⊙ and luminosities close to the Eddington limit. Because many possibilities must be considered, we use simple methods and often sacrifice mathematical precision to simplicity and generality. Vertical disc structure is represented by algebraic rather than differential relations, permitting easy derivation of scaling laws. In general, however, we are disappointed in our original hope of finding an hospitable theoretical environment for extended discs.
The plan of the paper is as follows. In Section 2 we review self-gravity in steady α discs, including irradiation from the central source (Section 2.3), and enhancements to α by local gravitationally-driven turbulence in a moderately self-gravitating disc (Section 2.4). In Section 3 we consider discs that are stabilized by additional heating beyond that due to the dissipation of orbital energy, but in which angular momentum continues to be transported by an α viscosity. In contrast to the usual situation, the total energy released per unit mass accreted is strongly dependent on the outer radius of the disc (Section 3). In Section 4 we briefly explore various alternatives to viscous thin-disc accretion. Schemes that involve a thin disc but invoke faster-than-viscous angular-momentum transport include global spiral waves (Section 4.5) and magnetized disc winds (Section 4.3). In every such case, we argue that the accretion velocity is bounded by the sound speed. Other alternatives explored in Section 4 are less like thin discs: quasi-spherical flows (Section 4.1), collisional star clusters (Section 4.2) and clumpy discs (Section 4.6).
All of these alternatives to the standard α disc face severe theoretical difficulties, or else seem unlikely to permit a centrifugally supported accretion flow beyond ∼0.1 pc. The situation is summarized in Table 1. We conclude in Section 5 that the gas probably does not circularize beyond this radius and must be supplied to the galactic nucleus with low specific angular momentum.
2 Steady α Discs
To begin, we review some salient features of a standard theoretical accretion disc: one that is steady, optically thick, and geometrically thin (Pringle 1981). The mass accretion rate is constant with radius and time. Angular momentum is carried inward by the accreting gas at the same rate that it is transferred outward by viscosity, and the radiation of energy from the surfaces of the disc is balanced by viscous dissipation.
Some of the assumptions above are relaxed in later sections. The following assumptions, however, apply throughout this paper. The rotation curve is dominated by a Newtonian point mass M, as relativistic corrections are important only at small radii where self-gravity is negligible. Elsewhere, though self-gravity may be important for local stability, the disc is thin and its mass <M so that the rotation curve is not much affected. We do not consider radii so large that the potential of the stars becomes important (r ≫ 1 pc). Except where otherwise stated (Section 4.3), magnetic pressure contributes little to the thickness of the disc even if magnetic stresses dominate angular-momentum transport, as indicated by numerical simulations of magnetorotational instability (MRI; Balbus & Hawley 1998). Contrary suggestions, however, do exist (Pariev, Blackman & Boldyrev 2002).Shakura & Sunyaev 1973) viscosity parameter α. Σ is the surface mass density, and β≡pgas/p =pgas/(pgas+prad) is the ratio of gas pressure to total pressure at the midplane. The mechanism likely to be responsible for the viscosity of most discs is magnetorotationally-driven turbulence, for which simulations indicate α = 10−3–10−1 (Balbus & Hawley 1998). As discussed below, where the disc is self-gravitating, α may be as large as ∼0.3. Therefore, we often write α0.01≡ 102α or α0.3≡α/0.3. The parameter b is a switch that determines whether the viscosity is proportional to gas pressure (b = 1) or total pressure (b = 0). In the latter case, radiation-pressure-dominated regions of α discs are viscously unstable (Lightman & Eardley 1974), but it is possible that b = 0 in an average sense. Although the average surface density would be lower and the corresponding higher for b = 0 than for b = 1, the viscous instability is expected to produce overdense rings in which these trends may well be reversed. Miller & Stone 2000), so that the factor should be even smaller. Hence equation (5) probably overestimates the midplane temperature. Other things being equal, this would cause self-gravity to be underestimated, but we assume that the error is not significant.
2.1 Onset Of Self-Gravityequation (6) with equation (2), we obtain the simple relation
To see why self-gravity is inevitable, consider the radiation-pressure-dominated case, β≪ 1, so that the isothermal sound speed c2s = 4σT4/3cρ. From equation (3) and (5), . Eliminating T between these relations and using Σ = 2csρ/Ω leads toequation (4) for in equation (7) leads to Section 2.4).
In the appendix, we briefly derive the radial run of Σ and T that follow from the standard assumptions laid out in this section. From these we derive expressions for β(r) (equation A3) and Q(r) (equation A4). β increases monotonically with radius, and Q decreases. For typical parameters, pgas =prad at about the same radius where Q = 1, r∼ 103RS∼ 10−2 pc. This appears to be a coincidence, as the two radii rQ =1 and rβ =1 depend differently on the accretion parameters, especially α and . We generally presume . But if the disc were not in a steady state or if there were a massive wind, then the accretion rate at small radii, where most of the QSO luminosity presumably originates, could be very different from at large radii where self-gravity is problematic. Therefore, it is interesting to calculate the local value of that would yield Q = 1 at each radius. We begin by expressing β in terms of Q using equations (4), (6) and (7):equation (A3), this equation does not derive from equation (3) and (5); it does not assume a balance between radiative losses and viscous dissipation. The opacity κe.s. enters Equation (11) only because it is used in the definition of the Eddington ratio lE. Equation (11) remains true regardless of the heat source, whereas equation (A3) is specific to viscous dissipation.
Eliminating β between these two equations, we obtain loci of constant Q in the plane, as shown in Fig. 1. These have been drawn for constant opacity κ = κe.s. regardless of temperature, but, because of the very strong dependence of on r, we expect that a more realistic opacity law would yield similar results. The unstable region where Q < 1 is the interior of the wedge. Thus, below some radius depending on α and M, there are no gravitationally unstable solutions for any . Beyond this critical radius, there are two marginally stable solutions, which diverge very quickly from one another in . The upper branch is very strongly radiation-pressure dominated (β≪ 1) and very optically thick (κΣ≫ 1), so that it bears a resemblance to the supercritical solutions studied by Begelman & Meier (1982). Their solutions, however, besides assuming constant , are quasi-spherical and advection dominated out to the ‘trapping radius’Fig. 1 are calculated from the thin-disc equations assuming local balance between radiative losses and viscous heating. The relative thickness h/r increases with r on the upper branch. For b = 0 and M8 = 1, h/r = 0.1 at 1.2 pc and h/r→ 1 at 7.3 pc; for b = 1, h/r = 1 is already reached at 0.11 pc. Thus, the vertical limit of both plots has been set approximately where the thin-disc assumption breaks down.
2.2 Time Dependence
Steady-state accretion is mathematically convenient but the evidence for it is scant, because no QSO has been monitored for a significant fraction of its growth time (equation (1)). The observed luminosity constrains at small r, but at large r could be different. Suppose the black hole shines at L ∼LE with a duty cycle f ≪ 1, and that it scarcely accretes at other times. A steady inflow rate ∼fLE/εc2 at large radii is sufficient to support this behaviour, if ε is the radiative efficiency during the bright phase. Assuming that the black hole gains its mass by accretion, f≳tSal/t≳ 10−2, where t is the age of the universe when the QSO shone, estimated for zQSO≳ 2 and standard cosmological parameters. Thus, the accretion rate in the outer disc can typically be reduced by at most two orders of magnitude compared to the peak rate ((equation 4)). From Fig. 1, however, we see that a reduction of at least four orders of magnitude is needed to stabilize a disc at 0.1 pc if M8 = 1.
From this we conclude that a variable accretion rate will not solve the problem of self-gravity in the outer disc.
Although flared or warped outer regions can be warmed by irradiation from the inner parts (where most of the total disc luminosity originates), we now show that this effect is not enough to stabilize the disc at r ∼ rQ = 1.
To do so, irradiation must raise the midplane temperature T, not just the surface temperature Teff, substantially. As a result, the disc will be nearly isothermal from surface to midplane. The vertical pressure gradient, and therefore the disc thickness, will be dominated by the gas even if β≪ 1: that is, h2≈pgas/ρΩ2. If the effective viscosity derives from local magnetorotational turbulence, then it probably scales more directly with thickness than with sound speed because MRIs do not require compression. Then in all relevant respects the disc behaves as if β = 1 and cs2 → pgas/ρ = kBT/m. From equation (7), it then follows that the minimum temperature that the irradiation must provide to ensure gravitational stability is
2.4 Local Gravitational Turbulence
It has occasionally been proposed that a partially self-gravitating disc may transport angular momentum by spiral waves (e.g. Cameron 1978; Paczyński 1978; Larson 1984; Lin & Pringle 1987). Others have suggested that gravitational instabilities are intrinsically global and therefore not reducible to a local viscosity (Balbus & Papaloizou 1999). It may be that either opinion can be correct depending upon the ratio of disc thickness to radius, because the most unstable wavelength in a Q = 1 disc is ∼h. We consider local gravitational turbulence here, as it can be accommodated within the α model; a (much larger) upper bound on angular momentum transport by global spirals is given in Section 4.5.
By careful two-dimensional (2D) simulations, Gammie (2001) finds that gravitationally-driven turbulence can be local and can support α approaching unity; in fact, in the absence of any other viscosity mechanism, his discs self-regulate themselves so thatequation (14) differs from that of Gammie by a factor γ2D because we define α in terms of isothermal rather than adiabatic sound speed.] Gammie finds that self-regulation fails and the disc fragments if Ωtth≲ 3 for γ2D = 2. The latter result is not entirely surprising; in the absence of pressure support, the natural time-scale for collapse is tdyn = Ω−1, and we would not expect collapse to be prevented by thermal energy unless tth > tdyn. The converse statement, that fragmentation can be indefinitely postponed if Ωtth > 3, is not at all obvious but is supported by Gammie's simulations.
We are not aware of any direct numerical simulations of discs that are unstable to both gravitational and magnetorotational modes at the same radii, as is likely to be the case for QSO discs (Menou & Quataert 2001). It would be interesting to explore this, as there might be a synergy between the two instabilities.
QSO discs are expected to be very thin in the regions of interest to us, so that Gammie's results may be applicable. Equation (8) implies h =lERS/6ε in a radiation-pressure dominated disc, and hence h/r ≲ 0.003 at the innermost radius where Q ∼ 1 (equation (10)). We expect that, from this radius outward, α will rise smoothly from the value supported by magnetorotational turbulence (perhaps αm.h.d. ∼ 10−2) up to the maximum allowed by gravitational turbulence, αgrav,max ∼ 0.3, in such a way that Q ≈ constant. It follows from equation (9) that the ratio of outer to inner radii of this region is rather modest: ∼(αgrav,max/αm.h.d.)2/9∼ 2. At still larger radii, additional sources of energy are required in order to prevent catastrophic fragmentation.
3 Constant-Q Discs
It is not a new idea that accretion discs may regulate themselves at Q≈ 1. The heating required to offset radiative cooling and thereby maintain Q is usually supposed to come at the expense of the orbital energy of the gas (e.g. Paczyński 1978; Gammie 2001). In that case, gravitational instabilities and turbulence can be subsumed into the α parameter, at least to the extent that they act locally (Section 2.4). Bertin (1997) and Bertin & Lodato (1999) have suggested replacing the energy equation (3) with the constraint of constant Q (equation (6)), while retaining α to describe angular-momentum transport (equation (7)). In this section we adopt that approach and work out some consequences for the structure of the disc, with due attention to the effects of radiation pressure, which Bertin et al. seem to have neglected. It begs the question where the energy required to heat the disc may come from. Bertin & Lodato (1999) and especially Lodato & Bertin (2001) appear to suggest that the energy might ultimately derive from the gas orbits but yet be independent of local angular-momentum transport. We do not understand how this is possible. Instead, we attribute any surplus of the local radiative cooling above viscous dissipation to thermonuclear sources in stars, or perhaps to accretion on to small black holes that these stars evolve into. In Section 3.2 we calculate this surplus in terms of and the outer radius of the disc, and compare it to what might be available by converting a significant fraction of the disc gas into stars.
Although they may seem artificial, these assumptions are probably appropriate for the discs of spiral galaxies, and for the local interstellar medium (ISM) in particular. The local Galactic magnetic field is consistent with simulations of magnetorotational turbulence: namely, a magnetic energy density somewhat less than the thermal pressure, a predominantly toroidal orientation, and fluctuations comparable to the mean (Brandenburg et al. 1995). It is plausible therefore that there is a non-zero average magnetic stress −BrBϕ/4π = αpgas that systematically transfers angular momentum outward (Sellwood & Balbus 1999). Taking pgas/kB≈ 2000 cm−3 K−1, ρ≈ 0.3 mH cm−3 (Spitzer 1978), and circular velocity V0≈ 200 km s−1 (Binney & Merrifield 1998), the implied radial drift velocity vr≈−αpgas/ρV0≈−0.3 α km s−1 is small enough to have escaped detection. Perhaps coincidentally, equation (7) predicts , about half the Eddington rate for the Galaxy's 2.5 × 106 M⊙ central black hole. The implied viscous heating rate αpgasΩ0≈ 2 × 10−29α0.1 erg cm−3 is negligible compared to the inferred radiative cooling rate of the gas, ∼2 × 10−26 erg cm−3 (Spitzer 1978). Presumably, the temperature of the ISM is maintained by stars. An important difference between the local ISM and QSO discs is that the former is very optically thin, especially to absorption, which means that the energy input from stars is inefficiently radiated. We estimate below that constant-Q QSO discs remain optically thick out to at least 0.1 pc. A more careful treatment will be given in a later paper using realistic opacity tables (Sirko & Goodman 2003).
3.1 Density and Temperature
The midplane density in a constant-Q disc follows from equation (6)equation (10)) is ∼10−8M−4/38 g cm−3. The ratio β/(1 −β) = 3cρ/4σT3 is determined by equation (11), which shows that β ≪ 1 at r≫ 10RS if Q ∼ O(1). The temperature itself is equation (6),
Beyond 104–105RS≈ 0.1–1 M8 pc, the above formulae predict T≲ 5000 K, so that the opacity κ ≪ κe.s. and the disc becomes optically thin. (This assumes M8 =lE = ε0.1 = α0.3 = 1. Dust will raise opacity again at T ≲ 1700 K.) The disc must then be supported by gas pressure, notwithstanding equation (11) which presumes that the radiation is trapped. But in a gas-pressure dominated disc, the minimum temperature for gravitational stability is equation (12), which is about an order of magnitude larger than predicted by the formulae above. Note that this temperature is nearly independent of the disc rotation curve, whether dominated by the mass of the black hole or of the bulge stars, as Ω is subsumed into Q and . So, in a marginally gravitationally stable disc, there must be an extended region where the temperature adjusts itself within a limited range (5000 K ≲T≲ 104 K) so that the disc is marginally optically thin. At the low densities relevant here, the maximum opacity is κmax≈ 10κe.s. and is achieved at T≈ 104 K (Kurucz 1992; Keady & Kilcrease 2000). Hence the outer edge of the region in question should end at Σ ≈ κmax−1 ≈ 0.3 cm2 g−1, which occurs (assuming Q = 1 and T = 104 K) atequation (32) and equation (33)) rather than , because rthin lies outside the black hole's sphere of influence (equation (25)). Optically-thin constant-Q discs have been considered by Bertin & Lodato (2001) in a somewhat abstract context, but apparently without regard to radiation pressure, which remains important for the vertical hydrostatic balance even at small τ; see Sirko & Goodman (2003).
The parts of the disc beyond rthin would radiate predominantly in optical emission lines with velocity widths ∼σbulge, so that they might be identified with the QSO narrow-line region. However, we now see that the energy required to maintain the disc at constant and Q ≳ 1 all the way out to rthin would be prohibitive.
3.2 Energeticsequation (3)), which is largest at small radii; given a torque-free inner boundary at r =rmin, the global efficiency is ε = ε′(rmin). But the constant-Q discs require additional energy inputs, so that ε′ is generally larger than RS/4r and tends in fact to increase with radius. equation (23) or (24) are valid is the smallest among the radii rthin, RN and rout, where the first is the radius at which a Q = 1 disc would become optically thin (equation (21)), the second is the radius within which the black hole dominates the potential, equation (25) has used the recently reported M–σbulge relation (Gebhardt et al. 2000; Merritt & Ferrarese 2001, see Section 4.2 below). Actually, equation (25) probably underestimates the distance to which the black hole dominates the circular velocity, because the density profile in the cusp of bright bulges is considerably shallower than r−2.
At r = 10 pc ≈ 106M8−1RS, equation (23) and (24) imply ε′ ≈ 0.003 M8−1/6 and ε′ ≈ 0.1, respectively. The former efficiency is barely compatible with thermonuclear burning even if all of the disc gas is processed through high-mass stars. The latter is unsustainable by stars and in fact comparable to the efficiency of the central engine.
Even if the energy required to maintain Q∼ 1 can be found, these discs may have severe thermal instabilities. In a steady viscously heated disc, the thermal time (defined as the internal energy per unit area divided by the power radiated per unit area) is tth∼ (αΩ)−1, hence longer than the local dynamical time if α < 1 (Pringle 1981). While this does not guarantee thermal stability, it at least ensures a clear distinction between thermal and dynamical perturbations. For these constant-Q discs, however, the thermal time is shorter by the ratio of viscous heating to total heating. In view of equation (22) and the remarks following it, we can write αΩtth∼ε′(r)(4r/RS), and from equations (23)–(24), it follows that tth ≪ Ω−1 at r∼ 1 pc. In such circumstances, it is far from clear whether Q = 1 is the correct stability criterion. Unfortunately, we cannot say much more about this without understanding how the auxiliary heat sources respond to changes in local density and temperature, and with what time lag.
4 Non-Standard Discs and Other Alternatives
As shown in Section 2, a geometrically thin, optically thick accretion disc in a typical high-luminosity QSO would be self-gravitating at radii r≳ 102–103RS, or 10−3–10−2 pc, where RS = 2GM/c2 is the Schwarzschild radius of the central black hole. Large global radiative efficiency, ε ≳ 0.1, probably requires a thin disc near RS, but it does not much depend upon the nature of the flow at large radius. So, in this section, we consider whether the self-gravitating part of the disc can be replaced by some other form of accretion.
4.1 Hot, Quasi-Spherical Flows
These have been proposed as models for accretion at very low accretion rate and low radiative efficiency (Rees et al. 1982; Ichimaru 1977; Narayan & Yi 1994). At sufficiently low density, the gas cooling time and ion—electron thermal equilibration time are longer than the accretion time, so that the ion temperature (Ti) is approximately virial and the thickness of the disc is comparable to its radius. We suppose that angular-momentum transport is efficient (viscosity parameter α∼ 1), or that the angular momentum of the gas is very small to begin with, so that the accretion velocity (vr) is comparable to the free-fall velocity. The large vr and Ti combine with the low to make the density low, as required. It is unclear just how low must be for this mode of accretion to sustain itself, because angular-momentum transport is not well understood and collisionless processes may enhance the thermal coupling of ions and electrons.
A quasi-spherical flow is not viable when the luminosity is close to the Eddington limit, however, because of inverse-Compton cooling. For spherical free fall on to a source of luminosity L = lELEdd, the inverse-Compton cooling rate of free electrons isequation (26), tC < tff at r < 25lE2M8 pc. Hence the electrons assume the colour temperature of the central source (or even less, see below): TC≲ 105M−1/48. This is much less than the virial temperature, Tvir≈ 107M8rpc−1, at all radii of interest to the present paper.1 The electron density is Spitzer 1978) is much shorter than the flow time: Spitzer 1978). If ne obeys equation (28), then radiative cooling is actually faster than inverse-Compton cooling:
In short, both the ions and the electrons of a quasi-spherical flow on to a near-Eddington QSO cool in much less than a free-fall time. A thin disc will form unless the specific angular momentum of the flow is negligible.
4.2 Collisional Stellar Cluster
Dense stellar clusters have occasionally been nominated as precursors to QSO black holes, either by relativistic collapse (Zel'dovich & Podurets 1966; Shapiro & Teukolsky 1985; Ebisuzaki et al. 2001) or by collisions among non-degenerate stars (Spitzer & Saslaw 1966; Rees 1978). Our interest, however, is in stellar collisions as the main source of fuel for an already very massive black hole. This has been studied by McMillan, Lightman & Cohn (1981) and Illarionov & Romanova (1988), among others, who show that in order to supply a 108 M⊙ black hole at its Eddington rate, the velocity dispersion of such a cluster must be ≳ 103 km s−1. To demonstrate the robustness of this conclusion, we make some oversimplified but conservative estimates here.
In order to provide high radiative efficiency (ε), stars must be disrupted and their gaseous debris circularized before accretion. If Mb.h. ≲ 108 M⊙, main-sequence stars scattered on to loss-cone orbits are likely to be tidally disrupted rather than swallowed whole (Hills 1975). About half of the tidal debris is unbound and promptly escapes from the black hole (Lacy, Townes & Hollenbach 1982; Evans & Kochanek 1989), and much of the remainder is likely to be swallowed at low radiative efficiency before the gas circularizes (Cannizzo, Lee & Goodman 1990; Ayal, Livio & Piran 2000). We limit our discussion to Mb.h.≳ 108 M⊙, as required for the most luminous QSOs, and we assume that stars are disrupted by stellar collisions. We ignore loss-cone effects because stars swallowed whole do not contribute to the QSO luminosity.Spitzer & Saslaw 1966), but it will lead to a conservative estimate of . Gravitational focusing is unimportant at the high velocity dispersions relevant here. equation (29) and (30) leads to
The coefficient C depends upon the density profile ρ(r) of the stellar cluster. It can be arbitrarily large if ρ(r) is sufficiently steep, but then the collision rate is dominated by stars at small radius. These tightly-bound stars represent only a small fraction of the cluster mass and would be consumed in much less than the growth time of the black hole, unless their total mass (Mt.b.) is ≳ Mb.h., in which case the tightly-bound population is substantially self-gravitating and we redefine M*≡Mt.b.. We might suppose that Mt.b. ≪ Mb.h. and that the tightly-bound stars were continuously replenished by two-body relaxation from a reservoir of more weakly bound stars; but the larger cluster would then have to expand to conserve energy, reducing the collision rate and the fueling of the QSO. For these reasons, we consider clusters for which C is dominated by stars near the half-mass radius. For example, if the stars have isotropically distributed orbits with a common semimajor axis a in a Keplerian potential, then , , and C≈ 196. This leads toMerritt & Ferrarese (2001) use the central velocity dispersion, which is usually little different from σe: equation (31) is similar to the empirical equation (32) and (33), but the normalization is very different. Extrapolated to 700 km s−1, the empirical relations predict instead of 108 M⊙.
Therefore, if QSOs were fuelled by dense stellar clusters, these clusters must have been an order of magnitude more tightly bound than the surrounding bulge, and they would have been a dynamically distinct stellar component. There seems to be very little trace of this tightly-bound stellar population in present-day bulges. How would such a component form? A likely possibility is gaseous dissipation followed by star formation. As will be seen, accretion in a thin viscous disc may lead to just such a result. We have discussed fuelling the QSO by a disc and fuelling it by stellar collisions as though these were mutually exclusive possibilities, but perhaps the two occur in concert.
4.3 Wind-Driven Discs
In principle at least, a magnetized wind can remove angular momentum from a thin disc rather efficiently. Compared to a viscous disc of the same sound speed (cs) and accretion rate , a wind-driven disc might have an accretion velocity that is larger by a factor ∼rΩ/cs = r/h. The surface density would be correspondingly reduced, as would the tendency toward self-gravity.
The main point to notice about equation (34) is that the area of the top and bottom faces is much larger than that of the sides. This means that, for the same field strength, a magnetized wind can be much more effective at removing angular momentum from the accreting gas than a field confined to the gas layer. The accretion speed driven by a magnetic ‘viscosity’,equation (34) by a factor ∼(h/r)(Br/Bz).
Viscous stresses are normally presumed to be smaller than thermal pressure, i.e. α < 1. For a magnetic viscosity, , where α < 1 if the field strength is less than the equipartition value. A superequipartition field would shut off the crucial MRI (Balbus & Hawley 1998), and would probably escape the disc by interchange or Parker instabilities. We now give two arguments to show that, even in a wind-driven accretion disc, the field must be at or below equipartition.Blandford & Payne (1982), in order that centrifugal force should drive the wind outward, the poloidal field lines must make an angle of at most 60○ with the surface of the disc, so that . Therefore equation (34), we have
The angular momentum extracted from the face of the disc by the field must be carried off by the wind. It is problematic whether a strong wind can be launched from the disc (Ogilvie & Livio 2001), and whether rapid wind-driven accretion is stable (Cao & Spruit 2002). Apart from these difficulties, the consideration of the magnetic flux,
To summarize this subsection, accretion driven by magnetized winds is even less well understood than viscous accretion but might allow substantially higher accretion velocities and lower surface densities, perhaps by factors up to ∼r/αh.
4.4 Thin Discs With Strongly Magnetized Coronae
This is a variant of Section 4.3 in which most of the field lines are not open but re-attach to the disc at large distances Δr≫h (Galeev et al. 1979; Heyvaerts & Priest 1989). The vertical magnetic scale height is then H ∼ Δr ≫ h. The angular momentum flux carried through the corona,
The evidence for magnetized coronae that dominate angular-momentum transport is suggestive but inconclusive. Local simulations of MRI generally find that the scale height of the field exceeds that of the gas, but only by factors of order unity; they also find α∼ 10−2–10−1 rather than ∼r/h (Brandenburg et al. 1995; Stone et al. 1996; Miller & Stone 2000). Possibly, H/h is limited by the fact that the smallest dimension of the computational domain is ≲h. Global simulations of MRI have been performed for relatively thick discs only, so that it is difficult to distinguish scalings with h from scalings with r (Matsumoto & Shibata 1997; Hawley 2000). Global simulations of thin discs in three dimensions may not be available for some time because of the very large numbers of grid cells needed to resolve both the disc and the corona. Merloni & Fabian (2001) argue that X-ray observations of accreting black holes — active galactic nuclei (AGN) and X-ray binaries — demand a strongly magnetized corona, at least in the innermost part of the disc. On the other hand, observations of eclipsing cataclysmic variables indicate that X-rays are emitted from the disc—star boundary layer rather than an extended corona (Mukai et al. 1997; Ramsay et al. 2001).
4.5 Global Spiral WavesLynden-Bell & Kalnajs 1972). The above approximation is good for tightly-wrapped waves, μ ≫ m ≠ 0, which carry relatively little flux for a given density contrast Σm/Σ0. An exact formula for logarithmic spirals with Σm ∝ r−3/2 is Kalnajs 1971):
Suppose that the gravitational torque is balanced by the advection of angular momentum with the accreting gas. In other words, , so that there is no secular change in the angular momentum within radius r. The gravitationally-driven accretion speed is thenLee & Goodman (1999), but only for weak self-gravity (Q ≫ 1) and without dissipation or accretion.
In steady accretion on to a central mass that dominates the rotation curve, the advected angular-momentum flux , so Σ∝r−5/4 rather than r−3/2. Presumably the slight change in power-law index does not change the results (40) and (41) much.
4.6 Clumpy Discs
Can QSO discs persist at Q ≪ 1 without fragmenting entirely into stars?Lin & Pringle (1987) have suggested that complete fragmentation does not necessarily result from such a low Q. Numerical simulations suggest that it does (Monaghan & Lattanzio 1991; Gammie 2001).
The question is all the more urgent because of rather direct evidence for parsec-scale accretion discs in nearby AGN, if not QSOs, from VLBI observations of maser emission. If the nuclear disc of NGC 1068 is in a steady state, then the nuclear luminosity implies Q ∼ 10−3 at r ∼ 1 pc (Kumar 1999). NGC 4258 is much less luminous, and estimates of range from based on modeling the maser emission itself (Neufeld & Maloney 1995), to (Gammie, Narayan & Blandford 1999; Kumar 1999) for an assumed central advection dominated accretion flow (ADAF); at the former rate, Q∼ 1 at the outer edge of the masering region, ∼0.2 pc (Maoz 1995), while in the latter, Q∼ 10−2.
Kumar (1999) has suggested a clumpy rather than smooth disc, in which accretion occurs by gravitational scattering and physical collisions among clumps rather than an α viscosity. These clumps are supposed to be gas clouds rather than fully formed stars, in order to provide appropriate conditions for maser amplification. Although the model deals with the stability and accretion rate of the clumpy disc as a whole, it does not ask what prevents the individual clumps from collapsing. In the application to NGC 1068, the masses, radii and surface temperatures of the clumps are quoted as Mc ∼ 103 M⊙, Rc ∼ 0.1 pc and Teff,c ≈ 500 K; the virial temperature implied by this mass and radius is ∼2000 K. These clumps would be moderately optically thick at the wavelength corresponding to Teff,c. Their characteristic thermal time — the time required to radiate their binding energy — is
To answer the question raised above, we therefore believe that no QSO accretion disc, whether smooth or clumpy, can persist at Q ≪ 1.
Although the accretion velocity of a thin viscous disc is very strongly subsonic, with a Mach number , the discussion of Section 4 suggests that might approach ∼0.1 if accretion is driven by large-scale magnetic or gravitational fields rather than a local effective viscosity. Self-consistent solutions that achieve this bound would be very interesting to pursue.
Even at near-sonic accretion speeds, QSO discs become self-gravitating at radii less than a parsec if there is no important source of heating other than dissipation of orbital energy. Combining equation (3) and (5) with instead of the viscous relation (2), and assuming that gas pressure dominates (β ≈ 1), we find that
On the other hand, we have seen (Section 3.2) that, when angular-momentum transport is viscous, then it is unlikely that stars can supply enough additional heat to stabilize the disc beyond 1 pc, especially if viscosity is proportional to gas pressure rather than total pressure as viscous stability probably demands. In the latter case, even low-mass black holes embedded in the disc are probably inadequate. These statements assume that enough free gas remains in the disc to supply the central black hole at its Eddington rate. If all the gas converts to stars, stability may result, but the quasar is quenched. Perhaps a combination of stellar (or embedded-black-hole) heating and a superviscous accretion speed may allow an extended gravitationally stable disc; we hope to explore this possibility in a future paper.
Given the serious theoretical difficulties of all proposed mechanisms for keeping Q ≳ 1 at large radii, we are forced to take seriously the only remaining possibility: that QSO discs do not exist much beyond Rsg ∼ 10−2 pc (equation (10)) — at least not in a state of centrifugal support, vertical hydrostatic equilibrium, and steady accretion. Yet the mass within this radius is smaller than that of the central black hole by a factor ∼(h/r) ∼ 10−2 (since the midplane density ≈ M/2πr3 at Q = 1). Hence, in order to grow the black hole, the disc must be replenished, either steadily or intermittently, by infall of low-angular-momentum material. In order that the gas not circularize outside 10−2rpc, its specific angular momentum must be ≲ 70 M1/28 km s−1 pc, or some three orders of magnitude smaller than that of most stars in ellipticals and bulges (Binney & Merrifield 1998). Such a small ratio may seem unlikely, but on the other hand, Mbh/Mbulge ≈ 10−3 (McLure & Dunlop 2002). So perhaps the QSO is fuelled by the low-angular-momentum tail of gas that forms the bulge. This gas would arrive at the outer edge of the disc in a vertically broad infall, perhaps already carrying dust formed from metals injected by outflows from the bulge or QSO disc itself, and hence taking the place of a warped outer disc as the source of reprocessed infrared light. The picture this calls to mind is similar, except in scale, to the standard scenario for the formation of a protostar (Shu, Adams & Lizano 1987).
I thank Charles Gammie, Pawan Kumar, Kristen Menou, Ramesh Narayan, Bohdan Paczyński and Jonathan Tan for helpful discussions. This work was supported by the NASA Origins Programme under grant NAG 5-8385.
Appendix A: Viscously Heated Discs
For completeness, in this appendix we give formulae for the midplane temperature (T) and other characteristics of a steady optically thick disc heated by viscous dissipation only.equations (A1) and (A2) do not depend on β, which in any case is a known function of density and temperature: equation (A1) to eliminate cs from equation (7) yields the gravitational stability parameter