One of the key features observed in accretion discs is the strong and often chaotic time variability in the X-ray spectrum, known as
Quasi-periodic oscillations (QPOs). They were first noticed in dwarf novae, i.e., erupting (
cataclysmic variables) as incoherent pulses with time scales of 30-170s, along with coherent periodic oscillations with a period of 20s, the so called dwarf nova oscillations (DNOs).
These flux variations occur at all sorts of time scales and some of them are supposedly caused by an ensemble of waves and oscillation modes in the innermost region of the accretion disc. The field of studying the temporal behaviour of discs by means of oscillations is called
discoseismology
In general, oscillations are the result of restoring forces acting on perturbations. For instance, if one perturbs a fluid element radially inwards, it conserves its own angular momentum and will be rotating too slow for its new location. Centripetal forces consequently push it outwards again. These kind of
inertial oscillations are called
epicyclic oscillations.
One can then express the Eulerian perturbations of all physical quantities through a single function <math>\delta W \propto \delta p/\rho</math> which satisfies a second-order partial differential equation. Since the accretion disc is considered to be stationary and axisymmetric, the angular and time dependences are factored out as <math>\delta W = W(r,z)e^{i(m\phi - \sigma t)}\ ,</math> where the eigenfrequency <math>\sigma(r,z) = \omega - m\Omega</math> and <math>m</math> is the azimuthal wave number. It is assumed that the variation of oscillation modes in radial direction is much stronger than in vertical direction. The resulting in two partial differential equations for the functional amplitude <math>W(r,z) = W_r(r) W_y(r,y)\ ,</math> where the radial eigenfunction, <math>W_r\ ,</math> varies fast with <math>r</math> and the vertical eigenfunction, <math>W_z\ ,</math> varies slowly with <math>r\ ,</math> are give by
<math>\frac{d^2W_r}{dr^2}-\frac{1}{(\omega^2-\omega_r^2)}[\frac{d}{dr}(\omega^2-\omega_r^2)]\frac{dW_r}{dr} +\alpha^2(\omega^2-\omega_r^2)(1-\frac{\Psi}{\tilde{\omega}^2})W_r = 0</math> and
<math>(1-y^2)\frac{d^2W_y}{dy^2}-2gy\frac{dW_y}{dy}+2g\tilde{\omega}^2[1-(1-\frac{\Psi}{\tilde{\omega}^2})(1-y^2)]W_y = 0</math>
where <math>\omega_r(r)</math> and <math>\omega_{\theta}(r)</math> are the radial and vertical epicylic frequecy, respectively, <math>\tilde{\omega}(r) = \omega(r)/\omega_{\theta}(r)</math> and <math>\alpha(r) = \frac{dt}{d\tau}\frac{\sqrt{g_{rr}}}{c_s(r,0)}\ ,</math> using coefficients of the Kerr metric in Boyer-Lindquist coordinates. <math>\Psi</math> is the eigenvalue of the (WKB) separation function. The radial boundary conditions, e.g., depend on the type of mode and its capture zone (see below). Oscillations in accretion discs are studied by means of <math>\Psi(r,\sigma)</math> with the angular mode number <math>m\ ,</math> the vertical and radial mode numbers (number of nodes in the corresponding eigenfunction) <math>j</math> and <math>n\ ,</math> respectively.
</td>
</td>
|
500px
|
|
500px
|
| Figure 1: Schematic picture showing trapped axisymmetric g-modes with n=1 in the region between <math>r_1</math> and <math>r_2\ ,</math> below the radial epicyclic mode <math>\kappa(r)=\omega_r(r)\ .</math> p-modes only exist above the Keplerian frequency <math>\Omega_K\ .</math> (Figure credits: Kato et al., 1998) | |
Figure 2: Schematic picture showing the dependence of three characteristic disc frequencies on the spin <math>a</math> of the black hole. (Figure credits: Wagoner, 1999) |
Classification:
A mode oscillates in the radial range outside the inner disc, <math>r > r_i\ ,</math> where the quantity <math> (\omega^2 - \omega_r^2) (1 - \frac{\Psi}{\tilde{\omega}^2}) </math> is positive. The two points <math>r=r_{\pm}(m, a, \sigma)</math> refer to the location where the
Lindblad resonances occur.
-
p-modes are inertial acoustic modes defined by <math>\Psi < \tilde{\omega}^2</math> and are trapped where <math>\omega^2 > \omega_r^2</math> in two zones between the inner (<math>r_i</math>) and outer (<math>r_o</math>) radius of the disc, <math>r_i < r < r_- </math> or <math>r_+ < r < r_o\ .</math> The latter, shown in Figure 1, produce the stronger luminosity modulation. In the corotating frame these modes appear at frequencies slightly higher that the radial epicyclic frequency. Pressure is their main restoring force.
-
g-modes are inertial gravity modes defined by <math>\Psi > \tilde{\omega}^2</math> and are trapped where <math>\omega^2 < \omega_r^2</math> in the zone <math>r_- < r < r_+\ .</math> In the corotating frame these modes appear at low frequencies. Gravity is their main restoring force. - In Figure 1 <math>r_1=r_-</math> and <math>r_2=r_+\ ,</math> in figure 2 the spin dependency is plotted for <math>m=0\ .</math>
-
c-modes are corrugation modes defined by <math>\Psi = \tilde{\omega}^2\ .</math> They are non-radial (<math>m=1</math>) and vertically incompressible modes that appear near in inner disc edge and precess slowly around the rotational axis. In the corotating frame these modes appear at highest frequencies.
All modes have frequencies <math>\propto 1/M\ .</math> Since g-modes are trapped in the region of the temperature maximum of the disc, they are expected to be observed best.
- Wagoner, R., 1999, "Relativistic Diskoseismology", PhR, 311, 259
- Kato S., 2001, "Basic Properties of Thin-Disk Oscillations", PASJ, 53, 1
Geometrically thick accretion discs, i.e. accretion tori, always allow axisymmetric, incompressible modes corresponding to global oscillations of the entire torus at radial (<math>\sigma=\omega_r</math>) and vertical (<math>\sigma=\omega_{\theta}</math>) epicyclic frequencies. Other possible modes, provided <math>m=0\ ,</math> are basically acoustic (p-), surface gravity (g-) and internal inertial (c-) modes and can be found by solving the relativisic Papaloizou-Pringle equation
<math>
\frac{1}{(-g)^{1/2}}\left\{\partial_{\mu}\left[(-g)^{1/2}g^{\mu\nu}f^n\partial_{\nu}W\right]\right\} - \left(m^2g^{\phi\phi} - 2m\omega g^{t\phi} + \omega^2 g^{tt}\right)f^n W
= -\frac{2nA(\bar{\omega}-m\bar{\Omega})^2}{\beta^2 r^2_0}f^{n-1}W
</math>
together with the boundary condition that the Lagrangian perturbation in pressure at the unperturbed surface (<math>f=0</math>) vanishes:
<math>
\Delta p = (\delta p + \xi^{\alpha}\nabla_{\alpha}p) = 0
</math>
Classification:
-
<math>\times</math>-modes are surface gravity modes (k=2) derived from an eigenfunction <math>W = a xy\ ,</math> for some constant <math>a\ ,</math> which is odd in <math>x</math> and <math>y</math> and results in two modes.
<math>
\bar{\sigma}_0^2 = \frac{1}{2}\{\omega_r^2+\omega_{\theta}^2 \pm [(\omega_r^2+\omega_{\theta}^2)^2+4\kappa_0^2\omega_{\theta}^2]^{1/2}\},
</math>
where <math>\bar{\sigma}_0^2=\sigma_0/\Omega_0</math> taken at the location <math>r_0</math> of the pressure maximum in the torus centre (hence the index 0). <math>\kappa_0^2 = \frac{\mathcal{E}_0^2}{l_0 A_0^2}(\frac{g^{tt}_{,r}-l_0 g^{t\phi}_{,r}}{g_{rr}}\frac{dl}{dr})</math> is the squared frequency of the inertial oscillation in the fluid due to an angular momentum gradient. For constant angular momentum distribution this term vanishes. The positive square root of which gives the x-mode that is a surface gravity mode (Figure 8). The negative square root gives a purely incompressible inertial (c-) mode whose poloidal velocity field represents a circulation around the pressure maximum.
-
breathing- and <math>+</math>-modes are derived from an eigenfunction <math>W = a + bx + cy\ ,</math> for some constants <math>a,b,c\ .</math> The resulting eigenfrequencies are (for <math>b=c=0</math>) the zero corotation frequency mode as well as
<math>
\bar{\sigma}_0^2 = \frac{1}{2n}\{(2n+1)(\omega_r^2+\omega_{\theta}^2)-(n+1)\kappa_0^2 \pm [((2n+1)(\omega_{\theta}^2-\omega_r^2)^2]^{1/2} \},
</math>
-
breathing - modes have frequencies corresponding to the upper sign in the above equation. The torus cross section contracts and expands (Figure 6). Breathing modes are comparable to acoustic modes (k=0, j=1) in the incompressible Newtonian limit for <math>l=const.\ ,</math> while in the Keplerian limit the mode frequency becomes that of a vertical acoustic wave.
-
<math>+</math>-modes have frequencies corresponding to the lower sign. In the incompressible <math>n\rightarrow 0</math> limit they are comparable to (k=2) gravity modes.
333px
Figure 3: Non-oscillating torus |
333px
Figure 4: Radial epicyclic oscillations |
333px
Figure 5: Vertical epicyclic oscillations |
333px
Figure 6: <math>\times</math>-mode oscillations |
333px
Figure 7: breathing mode oscillations |
333px
Figure 8: <math>+</math>-mode oscillations |
To
non-axisymmetric oscillation modes, however, accretion tori are dynamically unstable. The instability, discovered by Papaloizou & Pringle (1984), affects all torus configurations, and most violently tori with a constant angular momentum distribution.
Whether or not hydrodynamical oscillation modes may survive such global instabilities or the presence of a weak magnetic field (MRI turbulence), is subject of current, numerical investigations.
- Papaloizou, J., Pringle, J., 1984, "The dynamical stability of differentially rotating discs with constant specific angular momentum", MNRAS, 208,721
- Blaes O. et al, 2006, "Oscillation modes of relativistic slender tori", MNRAS, 369, 1235
