View source for Accretion discs/Analytic models of accretion discs/3.2 Thick discs

<table width="1111" border="0" cellpadding="0" cellspacing="0"><tr><td valign="top" align="left">

== Thick discs: assumptions ==

<div align="justify">For "thick discs" models of <nowiki>accretion discs</nowiki> one assumes that:
* Matter distribution is stationary and axially symmetric, i.e. matter quantities such as density <math>\epsilon</math> or pressure <math>P</math> are independent on time <math>t</math> and the azimuthal angle <math>\phi\ .</math>
* Matter moves on circular trajectories, i.e. the four velocity has the form <math>u^i = [u^t, u^{\phi}, 0, 0]\ .</math> The angular velocity is defined as <math>\Omega = u^{\phi}/u^t\ ,</math> and the angular momentum as <math>\ell = - u_{\phi}/u_t\ ,</math>
* <math>t_{dyn} \ll t_{the} <t_{vis}\ ,</math> with <math>t_{dyn}</math> being the dynamical timescale in which pressure force adjusts to the balance of gravitational and centrifugal forces, <math>t_{the}</math> being the thermal timescale in which the <nowiki>entropy</nowiki> redistribution occurs due to dissipative heating and cooling processes, and <math>t_{vis}</math> being the viscous timescale in which angular momentum distribution changes due to torque caused by dissipative stresses. Mathematically, this is equivalent to assume the stress-energy tensor in the form, <math>T^i_{~\nu} = u^{\mu}\,u_{\nu}\,(P + \epsilon) - \delta^{\mu}_{~\nu}\,P\ .</math> 

Using this form of the stress-energy tensor, [http://adsabs.harvard.edu/cgi-bin/nph-bib_query?bibcode=1978A%26A....63..221A Abramowicz et al. 1978] have derived from the <nowiki>equilibrium</nowiki> condition <math>\nabla_{\mu}\,T^{\mu}_{~\nu} = 0</math> the relativistic "Euler" equation,</div>
<!-- ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ -->
<table width="100%"><tr><td width="80%" align="left" valign="middle">
:<math> 
\frac{\nabla_{\mu} P}{\epsilon + P} = \frac{\nabla_{\mu}\,g_{tt} + 2\Omega\,\nabla_{\mu}\,g_{t\phi} + \Omega^2\,\nabla_{\mu}\,g_{\phi\phi}}{g_{tt} + 2\Omega\,g_{t\phi} + \Omega^2\,g_{\phi\phi}} = \nabla_{\mu}\ln A + \frac{\ell\,\nabla_{\mu}\Omega}{1 - \ell\,\Omega},
~~~{\rm with}~~~A^2(r, \theta) = \frac{1}{g_{tt}(r,\theta) + 2\Omega\,g_{t\phi}(r,\theta) + \Omega^2\,g_{\phi\phi}(r,\theta)}.
</math>
</td><td width="20%" align="right" valign="middle"><math>(3.1.1)</math></td></tr></table>
<!-- ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ -->

== Equipressure surfaces: analytic solution in a general case ==

<div align="justify">The "equipressure" surfaces are defined by an implicit condition <math>P(r, \theta) = const\ ,</math> which may be solved to get the explicit form <math>\theta = \theta(r)\ .</math> Then, (3.2.1) implies that the function <math>\theta(r)</math> obeys,</div>
<!-- ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ -->
<table width="100%"><tr><td width="80%" align="left" valign="middle">
:<math>
\frac{d\theta}{dr} = 
- \frac{\partial_{r}\,g_{tt} + 2\Omega\,\partial_{r}\,g_{t\phi} + \Omega^2\,\partial_{r}\,g_{\phi\phi}}
{\partial_{\theta}\,g_{tt} + 2\Omega\,\partial_{\theta}\,g_{t\phi} + \Omega^2\,\partial_{\theta}\,g_{\phi\phi}}
=
- \frac{\partial_{r}\,g^{tt} - 2\ell\,\partial_{r}\,g^{t\phi} + \ell^2\,\partial_{r}\,g^{\phi\phi}}
{\partial_{\theta}\,g^{tt} - 2\ell\,\partial_{\theta}\,g^{t\phi} + \ell^2\,\partial_{\theta}\,g^{\phi\phi}}.
</math>
</td><td width="20%" align="right" valign="middle"><math>(3.2.2)</math></td></tr></table>
<!-- ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ -->
<div align="justify">If <math>\ell = \ell(r, \theta)</math> or, which is equivalent, <math>\Omega = \Omega(r, \theta)</math> are known functions, the equation (3.2.2) for the equipressure surfaces takes the form of a standard ordinary differential equation <math>d\theta/dr = f(r, \theta)</math> with known rhs, and it may be directly integrated. This has been done first by [http://adsabs.harvard.edu/abs/1980AcA....30....1J Jaroszynski et al., 1980] and then by several other authors. In the Figure below, we show an example from a recent paper by [http://adsabs.harvard.edu/abs/2008arXiv0812.2467L Lei et al., 2008] who assumed the angular momentum distribution in the form that depends on three constant parameters (<math>\eta, \beta, \gamma</math>),</div>
<!-- ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ -->
<table width="100%"><tr><td width="80%" align="left" valign="middle">
:<math>
{\ell}(r, \theta) = {\ell}_0\left\{ \frac{{\ell}_K(r)}{{\ell}_0}\right\}^{\beta}(\sin\theta)^{2\gamma} ~\left({\rm for}~~ r \geq r_{ms}\right)
,~~
{\ell}(r, \theta) = {\ell}_{ms}(\sin\theta)^{2\gamma} ~\left({\rm for}~~ r < r_{ms}\right).
~~{\rm Here}~~ {\ell}_0 \equiv \eta\,{\ell}_K(r_{ms}) ~~{\rm and}~~ {\ell}_{ms} = {\ell}_0\left\{\frac{{\ell}_K(r_{ms})}{{\ell}_0}\right\}^{\beta}.
</math>
</td><td width="20%" align="right" valign="middle"><math>(3.2.3)</math></td></tr></table>
<!-- ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ -->
<!-- ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ -->
<table width="1111" cellspacing="2" cellpadding="0" border="0" style="background: #cccccc"><tr>
<td width="551" align="center" valign="top" style="background: #ffffff">[[Image:thick-disc-01.jpg]]
<table width="100%" cellspacing="8" cellpadding="0" border="0"><tr><td align="center" valign="middle">
<div align="justify"><b>Fig. 1</b> Equipressure surfaces for a non rotating black hole, and with angular momentum distribution (3.2.3) with <math>\gamma = 0.00\ ,</math> <math>\eta = 1.085\ ,</math> <math>\beta = 0.00\ .</math> Angular momentum is constant everywhere in space.</div>
</td></tr></table><br>
</td>
<td width="550" align="center" valign="top" style="background: #ffffff">[[Image:thick-disc-02.jpg]]
<table width="100%" cellspacing="8" cellpadding="0" border="0"><tr><td align="center" valign="middle">
<div align="justify"><b>Fig. 2</b> Equipressure surfaces for a non rotating black hole, and with angular momentum distribution (3.2.3) with <math>\gamma = 0.50\ ,</math> <math>\eta = 1.085\ ,</math> <math>\beta = 0.00\ .</math> Angular momentum is constant radially, but it changes along 
the polar angle.</div></td></tr></table><br>
</td>
</tr><tr>
<td width="551" align="center" valign="top" style="background: #ffffff">[[Image:thick-disc-03.jpg]]
<table width="100%" cellspacing="8" cellpadding="0" border="0"><tr><td align="center" valign="middle">
<div align="justify"><b>Fig. 3</b> Equipressure surfaces for a non rotating black hole, and with angular momentum distribution (3.2.3) with <math>\gamma = 0.18\ ,</math> <math>\eta = 1.085\ ,</math> <math>\beta = 0.9\ ,</math> compared with global MHD numerical simulations by 
[http://adsabs.harvard.edu/abs/2007ApJ...668..417F Fragile et al., 2007], time averaged over a period corresponding to the orbital period at <math>25r_G</math> (color-coded).</div></td></tr></table>
</td> 
<td width="550" align="center" valign="top" style="background: #ffffff">[[Image:thick-disc-044.jpg]]
<table width="100%" cellspacing="8" cellpadding="0" border="0"><tr><td align="center" valign="middle">
<div align="justify"><b>Fig. 4</b> Equipressure surfaces for a non rotating black hole, and with angular momentum distribution (3.2.3) with <math>\gamma = 0.00\ ,</math> <math>\eta = 1.085\ ,</math> <math>\beta = 0.99\ .</math> Angular momentum does not depend on the polar angle, its radial distribution is almost Keplerian.</div></td></tr></table>
</td>
</tr></table>
<!-- ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ -->

== Barytropic thick discs and the von Zeipel theorem ==

<div align="justify"> From equation (3.1.1) one proves (see e.g. [http://adsabs.harvard.edu/cgi-bin/nph-bib_query?bibcode=1978A%26A....63..221A Abramowicz et al. 1978]) that for barytropic fluids <math>\epsilon = \epsilon(P)\ ,</math> the surfaces of constant angular velocity and of constant angular momentum coincide, i.e. <math>\ell = \ell(\Omega)\ .</math> This is often called "the von Zeipel theorem". The equation of state <math>\epsilon = \epsilon(P)</math> and the rotation law <math>\ell = \ell(\Omega)</math> are independent and may be separately assumed. When they are known, the analytic solution is given by,</div>
<!-- +++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ -->
<table width="100%"><tr><td width="80%" align="left" valign="middle">
:<math> 
W(P) \equiv \int \frac{dP}{\epsilon(P) + P} = \ln A + \int \frac{d\Omega}{1 - \Omega\,\ell(\Omega)} \equiv \ln A(r, \theta) + F(\Omega),
~~~{\rm and}~~~\ell = \ell(\Omega) = 
- \frac{g_{t\phi}(r, \theta) + \Omega\,g_{\phi\phi}(r, \theta)}{g_{tt}(r, \theta) + \Omega\,g_{t\phi}(r, \theta)}.
</math>
</td><td width="20%" align="right" valign="middle"><math>(3.2.4)</math></td></tr></table>
<!-- +++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ -->
<div align="justify">The functions <math>W = W(P)\ ,</math> <math>F = F(\Omega)</math> and <math>\Omega = \Omega(r, \theta)</math> are explicitly known, and therefore one knows explicitly location of the equipressure surfaces <math>P = P(r, \theta)\ .</math> In the special (but important) case <math>\ell = \ell_0 = const\ ,</math> the function <math>\theta = \theta(r)</math> that gives the location of equipressure surfaces <math>P(r, \theta) = const</math> is given explicitly (for a non-rotating <nowiki>black hole</nowiki>)</div>
<!-- +++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ -->
<table width="100%"><tr><td width="80%" align="left" valign="middle">
:<math> 
\sin^2\theta = \frac{\ell_0}{C_0\,r^2 + 2\,G\,M\,\left(1 - r_G/r\right)^{-1}}. ~~~{\rm Here}~~C_0=const ~~{\rm numerates~the~equipressure~surfaces}.
</math>
</td><td width="20%" align="right" valign="middle"><math>(3.2.5)</math></td></tr></table>
<!-- +++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ -->
<div align="justify">Putting <math>r_G = 0</math> in (3.2.5), one recovers the Newtonian formula for <math>\ell = \ell_0 = const</math> discs in the <math>-G\,M/r</math> potential.</div>

== The Roche lobe overflow ==
<!-- +++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ -->
<table width="100%" cellpadding="0" cellspacing="0" border="0"><tr>
<td align="center" valign="middle" width="100%">
[[Image:roche-lobe-acl.jpg|center]]
</td></tr><tr><td align="left" valign="top" width="100%">
<div align="justify"><b>Fig. 5</b> Taken from [http://adsabs.harvard.edu/abs/1980ApJ...242..772A Abramowicz et al., 1980] <b>(a)</b> At the location <math>r = r_{in}\ ,</math> called the "cusp", angular momentum in the disc equals the Keplerian one, <math>\ell_{disc}(r_{in}) = \ell_K(r_{in})\ .</math> For <math>r > r_{in}</math> it is <math>\ell_{disc}(r_{in}) > \ell_K(r_{in})</math> and <math>d\ell/dr > 0\ ,</math> and for <math>r < r_{in}</math> it is <math>\ell_{disc}(r_{in}) < \ell_K(r_{in})</math> and <math>d\ell/dr \approx 0\ .</math> <b>(b)</b> The particular equipotential surface <math>W = W_{in}\ ,</math> called the "Roche lobe", crosses itself at the cusp. For <math>r \gg r_{in}</math> the surface of the disc (<math>P = 0</math>) coincides with the equipotential <math>W = W_{S}\ .</math> <b>(c)</b> The non-zero potential difference <math>\Delta W = W_S - W_{in}</math> implies that no <nowiki>equilibrium</nowiki> is possible at radii around and smaller than the cusp. Instead, there will be dynamical mass loss from the disk with the accretion rate that for a polytropic fluid <math>P = K\rho^{1+1/n}</math> equals ([http://adsabs.harvard.edu/abs/1978A%26A....63..209K Kozlowski et al. 1978], [http://adsabs.harvard.edu/abs/1985PASJ...37..727A Abramowicz 1985])</div>
<!-- ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ -->
<table width="100%"><tr><td width="80%" align="left" valign="middle">
:<math> 
{\dot M}_{in} = (2\pi)^{3/2}\frac{\Gamma(n + 3/2)}{(1 + 1/n)^n (1 + 1/2)^{n + 3/2}\Gamma(n + 3)}K^{-n}\frac{r_{in}}{\Omega_K(r_{in})}\Delta W^{n + 1},
~~{\rm with}~~ \Gamma(m) ~~{\rm being~the~Euler~gamma~function}.
</math>
</td><td width="20%" align="right" valign="middle"><math>(3.2.6)</math></td></tr></table>
<!-- ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ -->
<div align="justify">The mass loss (3.2.5) induced by the Roche lobe overflow self-regulates the accretion rate in the innermost part of all types of <nowiki>accretion disks</nowiki> (thin, slim, adafs, thick) around black holes and sufficiently compact neutron stars. This self-regulated overflow has several important consequences:
* It locally stabilizes <nowiki>accretion discs</nowiki> against thermal and viscous <nowiki>instabilities</nowiki> [http://adsabs.harvard.edu/abs/1981Natur.294..235A (Abramowicz, 1981)] and globally against the Papaloizou and Pringle <nowiki>instability</nowiki> [http://adsabs.harvard.edu/abs/1987MNRAS.227..975B (Blaes, 1987)]
* The amount of overflow, and therefore <math>{\dot M}_{in}\ ,</math> may be modulated by global discs <nowiki>oscillations</nowiki>. In the case of neutron stars, this leads to a modulation of the luminosity of the boundary layer at the neutron star surface. Although <nowiki>oscillations</nowiki> originate in the disk, the are observed in radiation that comes from the boundary layer. This is relevant for the observed neutron star quasi periodic <nowiki>oscillation</nowiki> (QPO) [http://adsabs.harvard.edu/abs/2007AcA....57....1A (Horak et al., 2007)].
* The "runaway" <nowiki>instability</nowiki> occurs when the mass exchange between <nowiki>accretion disc</nowiki> and <nowiki>black hole</nowiki> causes the cusp to move deeper into the disc, increasing the mass lost rate <math>{\dot M}_{in}\ .</math> This effect was suggested by [http://adsabs.harvard.edu/abs/1983Natur.302..597A Abramowicz et al. (1983)], and studied by several other authors (e.g. [http://adsabs.harvard.edu/abs/2004MNRAS.349..841D Daigne and Font, 2004] or [http://adsabs.harvard.edu/abs/2008PhRvD..78f4037M Montero et al., 2008]). It may determine the life-time of a massive torus around a black hole, which is relevant to some models of gamma ray <nowiki>bursts</nowiki>.</div>
</td></tr><tr><td align="center" valign="middle" width="100%">
[[Image:roche-lobe-overflow.jpg|center]]
</td></tr><tr><td align="left" valign="top" width="100%">
<div align="justify"><b>Fig. 6</b> Taken from [http://adsabs.harvard.edu/abs/1997MNRAS.284..767I Igumenshchev and Beloborodov, 1997]. The analytic formula (3.2.5) is very accurate, as a comparison with the numerical simulations shows. <b>(a)</b> The numerically calculated equipressure structure close to the central black hole is remarkably similar to that calculated analytically. In particular, there is obviously a cusp there. <b>(b)</b> In a qualitative agreement with the analytic model, the radial fluid velocity (described by arrows) is small far away from the cusp <math>r \gg r_{in}\ ,</math> but large (no <nowiki>equilibrium</nowiki>) at <math>r \approx r_{in}\ .</math> <b>(c)</b> Even more impressive is the excellent quantitative agreement of the predictions of the analytic formula (3.2.5), represented by lines, with results of numerical simulations, represented by points (circles, squares and triangles).</div>  
</td></tr></table>
<!-- +++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ -->

== Super-Eddington luminosities of radiation pressure supported thick discs ("Polish doughnuts") ==

<div align="justify">
The maximal energy available from an object with the mass <math>M</math> (and gravitational radius <math>R_G = GM/c^2</math>) is <math>E_{max} = M\,c^2\ .</math> The minimal time in which this energy may be liberated is <math>t_{min} = R_G/c\ .</math> Thus, the maximal power <math>L_{max} = E_{max}/t_{min} = c^5/G \equiv L_{Planck} = 10^{58}\,[{\rm erg/sec}] = 10^{52}\,[{\rm Watts}]\ .</math> An object with mass <math>M</math> has the "gravitational cross-section" <math>\Sigma_{grav} = 4\,\pi\,R^2_G\ .</math> If radiation interacts with matter by electron scattering (with the Thomson cross section <math>\sigma_T = 8\pi\,e^4/3\,m^2_e</math>), its "radiation cross section" is <math>\Sigma_{rad} = (M/m_P)\sigma_T\ .</math> M.&nbsp;Sikora noticed (unpublished) that the upper limit for the radiative power of an object in which gravity and radiation pressure are in <nowiki>equilibrium</nowiki> is given by the Planck power and the object gravitational and radiative cross sections,</div>
<!-- ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ -->
<table width="100%"><tr><td width="80%" align="left" valign="middle">
:<math> 
(\,\rm{Eddington~luminosity}\,) \equiv L_{Edd} = L_{Planck}\frac{\Sigma_{grav}}{\Sigma_{rad}} = \frac{4\pi\,G\,M\,m_P\,c}{\sigma_T} 
= 1.4 \times 10^{38}\left( \frac{M}{M_{sun}}\right) [{\rm erg}/{\rm sec}].
</math>
</td><td width="20%" align="right" valign="middle"><math>(3.2.7)</math></td></tr></table>
<!-- ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ -->
<!-- ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ -->
<table width="1111" cellpadding="0" cellspacing="0" border="0"><tr>
<td width="250" valign="top" align="right">[[Image:Madau-super-eddington.jpg]]</td><td width="11">&nbsp;</td><td width="850" align="left" valign="top"><div align="justify"> A "Polish doughnut" is a radiation pressure supported thick <nowiki>accretion discs</nowiki> around a central black hole. Polish doughnuts have toroidal shapes, resembling a large sphere (<math>R \gg R_G</math>) with a pair of empty narrow funnels along the rotation axis. The total luminosity of a Polish doughnut may be approximated as <math>L/L_{Edd} \equiv \lambda \approx \log (R/R_G)\ .</math> The logarithm here is of a crucial importance. It prevents astrophysically realistic doughnuts (i.e. with <math>R < 10^6\,R_G\ ,</math> say) to have highly super-Eddington luminosities. The theory predicts for such "realistic" fat tori slightly super-Eddington total (isotropic) luminosities <math>\lambda \le 7\ .</math> However, because the funnels have solid angles <math>\Theta^2 \sim R_G/R\ ,</math> radiation in the funnels may be, in principle, collimated to highly super-Eddington values <math>\lambda \sim R/R_G \gg 1\ .</math> This simple estimate agrees with a more detailed modeling of the Polish doughnuts radiation field by [http://adsabs.harvard.edu/abs/1981MNRAS.196..257S Sikora (1981)] and [http://adsabs.harvard.edu/abs/1988ApJ...327..116M Madau (1988)] who obtained <math>\lambda \ge 10^2</math> for discs with <math>R/R_G \sim 10^2\ .</math> A typical value that follows from observational estimates for the non-blazar active galactic nuclei, e.g. by [http://adsabs.harvard.edu/abs/1987ApJ...321..305C Czerny and Elvis (1987)] is <math>\lambda \sim 10\ ,</math> but of course for blazars and other similar sources (e.g. for ULXs, if they are powered by stellar mass black holes, as argued e.g. by [http://adsabs.harvard.edu/abs/2008NewAR..51..775K King (2008)], it must be <math>\lambda \gg 10\ .</math> <b>Fig. 7</b> (Credit: [http://adsabs.harvard.edu/abs/1988ApJ...327..116M Madau (1988)]) shows radiation from a Polish doughnut seen at different inclination angles, inclination <math>0^{\circ}</math> corresponds to line of sight along the funnel axis.</div></td>
</tr></table>
<!-- ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ -->
<!-- ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ -->
<table width="100%"><tr><td width="80%" align="left" valign="middle">
:<math> 
(\,\rm{Eddington~accretion~rate}\,) \equiv {\dot M}_{Edd} = \frac{L_{Edd}}{c^2} 
= 1.5 \times 10^{17}\left( \frac{M}{M_{sun}}\right) [{\rm g}/{\rm sec}], ~~~{\dot m} = \frac{\dot M}{{\dot M}_{Edd}}.
</math>
</td><td width="20%" align="right" valign="middle"><math>(3.2.8)</math></td></tr></table>
<!-- ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ -->
<div align="justify">Theoretical predictions about super-Eddington accretion rates:
* Radiation pressure supported black hole thick <nowiki>accretion discs</nowiki> ("Polish doughnuts") have typically super-Eddington luminosities <math>\lambda > 1\ .</math> 
* These discs have very small accretion efficiency and therefore must have highly super-Eddington accretion rates <math>{\dot m} \gg 1\ .</math>
* Super-Eddington accretion does not necessarily imply strong outflows.

The assumption that <math>{\dot m} = 1</math> is the upper limit for the growth rate of the seed black holes, often adopted in the context of the cosmic structure formation, is false.</div>
</td></tr></table>

<hr>