大质量黑洞吸积盘的径向自引力

本文研究了大质量黑洞吸积盘的自引力,用薄盘位形上积分的方法计算了吸积盘自引力的径向与垂向分量,着重讨论了径向自引力。主要结果为:对于大质量黑洞(M～10~8—10~(10)M_⊙)吸积盘,在(R/R_g)～10~5—10~4的距离上,径向自引力会超过中心天体引力。在这个距离上,吸积盘的动力学结构完全不同于开普勒盘。提出了径向自引力不稳定扰动作为一种能源机制。本文还得到吸积盘自引力与中心天体引力量级比较的两个判据,并由此得到大质量黑洞吸积盘外半径的近似解析估计。本文结果可用于类星体和星系核吸积盘。     

Accretion in Strong Gravity: from Galactic to Supermassive Black Holes

The galactic black hole binary systems give an observational template showing how the accretion flow changes as a function of increasing mass accretion rate, or L/L_Edd. These data can be synthesised with theoretical models of the accretion flow to give a coherent picture of accretion in strong gravity, in which the major hard-soft spectral transition is triggered by a change in the nature and geometry of the inner accretion flow from a hot, optically thin plasma to a cool, optically thick accretion disc. However, a straightforward application of these models to AGN gives clear discrepancies in overall spectral shape. Either the underlying accretion model is wrong, despite its success in describing the Galactic systems and/or there is additional physics which breaks the simple scaling from stellar to supermassive black holes.     

含平流双温吸积盘的振荡不稳定性研究

从流体动力学方程出发,用微扰法得出含平流双温吸积盘的径向、环向不稳定性的色散方程．并对平流和径向粘滞力对双温吸积盘的影响进行了较详细的讨论．结果表明：平流和径向粘滞力对声模有较大的影响,且不改变粘滞模和热模的稳定性质．而环向扰动对吸积盘的各种模有着较明显的作用．这一模型有利于解释活动天体的周期和准周期光变现象．     

黑洞热吸积流的理论研究

<正>活动星系核(AGN)、低光度活动星系核(LLAGN)以及X射线双星(XRB)系统中都存在黑洞吸积过程.黑洞吸积是人们理解相关天体的辐射、光变等现象的关键.本学位论文主要侧重于对目前非常流行的热吸积流(径移主导吸积流ADAF以及明亮热吸积流LHAF)的研究.该类吸积流是理解LLAGN以及处于硬态的XRB的基本理论模型.论文第1章详细介绍了相关的背景.首先讨论了天体环境中的黑洞及其所处的系统—AGN、     

Accretion disc outbursts: a new version of an old model

We have developed 1D time-dependent numerical models of accretion discs, using an adaptive grid technique and an implicit numerical scheme, in which the disc size is allowed to vary with time. The code fully resolves the cooling and heating fronts propagating in the disc. We show that models in which the radius of the outer edge of the disc is fixed produce incorrect results, from which probably incorrect conclusions about the viscosity law have been inferred. In particular we show that outside-in outbursts are possible when a standard bimodal behaviour of the Shakura–Sunyaev viscosity parameter α is used. We also discuss to what extent insufficient grid resolution has limited the predictive power of previous models. We find that the global properties (magnitudes, etc.) of transient discs can be addressed by codes using a high, but reasonable, number of fixed grid points. However, the study of the detailed physical properties of the transition fronts generally requires resolutions which are out of reach of fixed grid codes. It appears that most time-dependent models of accretion discs published in the literature have been limited by resolution effects, improper outer boundary conditions, or both.     

Accretion disc dynamics in extreme mass ratio compact binaries

An analysis is presented of a numerical investigation of the dynamics and geometry of accretion discs in binary systems with mass ratios   q = M ₂/ M ₁ < 0.1  , applicable to ultracompact X-ray binaries, AM CVn stars and very short period cataclysmic variables. The steady-state geometry of the disc in the binary reference frame is found to be quite different from that expected at higher mass ratios. For   q ∼ 0.1  , the disc takes on the usual elliptical shape, with the major axis aligned perpendicular to the line of centres of the two stars. However, at smaller mass ratios the elliptical gaseous orbits in the outer regions of the disc are rotated in the binary plane. The angle of rotation increases with gas temperature, but is found to vary inversely with q . At   q = 0.01  , the major axis of these orbits is aligned almost parallel to the line of centres of the two stars. These effects may be responsible for the similar disc structure inferred from Doppler tomography of the AM CVn star GP Com, which has   q = 0.02  . The steady-state geometry at low mass ratios is not predicted by an inviscid, restricted three-body model of gaseous orbits; it is related to the effects of tidal-viscous truncation of the disc near the Roche lobe boundary. Since the disc geometry can be inferred observationally for some systems, it is proposed that this may offer a useful diagnostic for the determination of mass ratios in ultracompact binaries.     

Unstable modes of non-axisymmetric gaseous discs

We present a perturbation theory for studying the instabilities of non-axisymmetric gaseous discs. We perturb the dynamical equations of self-gravitating fluids in the vicinity of a non-axisymmetric equilibrium, and expand the perturbed physical quantities in terms of a complete basis set and a small non-axisymmetry parameter ε. We then derive a linear eigenvalue problem in matrix form, and determine the pattern speed, growth rate and mode shapes of the first three unstable modes. In non-axisymmetric discs, the amplitude and the phase angle of travelling waves are functions of both the radius R and the azimuthal angle φ. This is due to the interaction of different wave components in the response spectrum. We demonstrate that wave interaction in unstable discs, with small initial asymmetries, can develop dense clumps during the phase of exponential growth. Local clumps, which occur on the major spiral arms, can constitute seeds of gas giant planets in accretion discs.     

Accretion to magnetized stars through the Rayleigh–Taylor instability: global 3D simulations

We present results of 3D simulations of magnetohydrodynamics (MHD) instabilities at the accretion disc–magnetosphere boundary. The instability is Rayleigh–Taylor, and develops for a fairly broad range of accretion rates and stellar rotation rates and magnetic fields. It manifests itself in the form of tall, thin tongues of plasma that penetrate the magnetosphere in the equatorial plane. The shape and number of the tongues changes with time on the inner disc dynamical time-scale. In contrast with funnel flows, which deposit matter mainly in the polar region, the tongues deposit matter much closer to the stellar equator. The instability appears for relatively small misalignment angles, Θ≲ 30°, between the star's rotation and magnetic axes, and is associated with higher accretion rates. The hotspots and light curves during accretion through instability are generally much more chaotic than during stable accretion. The unstable state of accretion has possible implications for quasi-periodic oscillations and intermittent pulsations from accreting systems, as well as planet migration.     

The dependence on system parameters of strange-mode instabilities in accretion discs

A systematic study of the dependence on disc parameters and input physics, such as opacity and the treatment of convection, of strange-mode instabilities in thin accretion discs, which have been discovered recently, is presented. The instabilities are found to exist for a wide range of parameters, are partly very robust, and their growth rates can reach the dynamical range. Even discs on galactic scales around massive black holes are affected by them. Two groups of instabilities can be distinguished, the first of which is related to the radiation-pressure-dominated part of the disc, and the second to helium/hydrogen ionization. By application of the NAR approximation, both of them can be shown to be of mechanical origin, and the classical κ -mechanism can be excluded as the instability mechanism. A heuristic model for strange-mode instabilities proposed in the context of stellar strange-mode instabilities in luminous stars seems to be applicable to the group associated with dominant radiation pressure.     

Accretion disc viscosity: how big is alpha?

We consider observational and theoretical estimates of the accretion disc viscosity parameter α. We find that in thin, fully ionized discs, the best observational evidence suggests a typical range α∼ 0.1–0.4, whereas the relevant numerical simulations tend to derive estimates for α which are an order of magnitude smaller. We discuss possible reasons for this apparent discrepancy.     

The radial–azimuthal instability of a hot two-temperature accretion disc with advection

The radial–azimuthal instability of a hot two-temperature accretion disc with advection is examined in this paper. We find that the inclusion of very little advection has significant effects on two acoustic modes for a geometrically thin, cooling-dominated two-temperature disc, but has no effect on acoustic modes for a geometrically slim, cooling-dominated two-temperature disc. We also find that, when azimuthal perturbations are considered, the stability properties of the disc are different from those in the pure radial perturbation case. An increase of the azimuthal wavenumber will stabilize the acoustic modes but make the viscous and thermal modes more unstable for a geometrically thin, cooling-dominated two-temperature disc. It makes the thermal mode more unstable and the acoustic mode more stable, but only affects the instability of the viscous mode for short-wavelength perturbations for a geometrically slim, cooling-dominated two-temperature disc. For a geometrically slim, advection-dominated two-temperature disc, the increase of the azimuthal perturbation makes the I- and O-modes more stable and the thermal mode more unstable, but has no effect on the viscous mode.     

Accretion disc coronae as magnetic reservoirs

Most astrophysical sources powered by accretion on to a black hole, either of stellar mass or supermassive, when observed with hard X-rays show signs of a hot Comptonizing component in the flow, the so-called corona , with observed temperatures and optical depths lying in a narrow range (0.1≲ τ ≲1 and 1×10⁹ K≲ T ≲3×10⁹ K). Here we argue that these facts constitute strong supporting evidence for a magnetically dominated corona. We show that the inferred thermal energy content of the corona, in all black hole systems, is far too low to explain their observed hard X-ray luminosities, unless either the size of the corona is at least of the order of 10³ Schwarzschild radii, or the corona itself is in fact a reservoir , where the energy is mainly stored in the form of a magnetic field generated by a sheared rotator (probably the accretion disc). We briefly outline the main reasons why the former possibility is to be discarded, and the latter preferred.     

The response of a turbulent accretion disc to an imposed epicyclic shearing motion

We excite an epicyclic motion, the amplitude of which depends on the vertical position, z , in a simulation of a turbulent accretion disc. An epicyclic motion of this kind may be caused by a warping of the disc. By studying how the epicyclic motion decays, we can obtain information about the interaction between the warp and the disc turbulence. A high-amplitude epicyclic motion decays first by exciting inertial waves through a parametric instability, but its subsequent exponential damping may be reproduced by a turbulent viscosity. We estimate the effective viscosity parameter, α _v, pertaining to such a vertical shear. We also gain new information on the properties of the disc turbulence in general, and measure the usual viscosity parameter, α _h, pertaining to a horizontal (Keplerian) shear. We find that, as is often assumed in theoretical studies, α _v is approximately equal to α _h and both are much less than unity, for the field strengths achieved in our local box calculations of turbulence. In view of the smallness (∼0.01) of α _v and α _h we conclude that for β p _gas p _mag∼10 the time-scale for diffusion or damping of a warp is much shorter than the usual viscous time-scale. Finally, we review the astrophysical implications.     

Evolution of initially localized perturbations in stratified ionized discs

A detailed solution of an initial value problem of a vertically localized initial perturbation in rotating magnetized vertically stratified disc is presented. The appropriate linearized magnetohydrodynamics equations are solved by employing the Wentzel–Kramers–Brillouin (WKB) approximation and the results are verified numerically. The eigenfrequencies as well as eigenfunctions are explicitly obtained. It is demonstrated that the initial perturbation remains confined within the disc. It is further shown that thin enough discs are stable but as their thickness grows increasing number of unstable modes participate in the solution of the initial value problem. However, it is demonstrated that due to the localization of the initial perturbation, the growth time of the instability is significantly longer than the calculated inverse growth rate of the individual unstable eigenfunctions.     

Magnetic shear-flow instability in thin accretion discs

The possibility that the magnetic shear-flow instability (also known as the 'Balbus–Hawley' instability) might give rise to turbulence in a thin accretion disc is investigated through numerical simulations. The study is linear and the fluid disc is supposed to be incompressible and differentially rotating with a simple velocity profile with Ω∝ R ^{− q} . The simplicity of the model is counterbalanced by the fact that the study is fully global in all three spatial directions with boundaries on each side; finite diffusivities are also allowed. The investigation is also carried out for several values of the azimuthal wavenumber of the perturbations in order to analyse whether non-axisymmetric modes might be preferred, which may produce, in a non-linear extension of the study, a self-sustained magnetic field.
  We find the final pattern steady, with similar kinetic and magnetic energies and the angular momentum always transported outwards. Despite the differential rotation, there are only small differences for the eigenvalues for various non-axisymmetric eigensolutions. Axisymmetric instabilities are by no means preferred; in fact for Prandtl numbers between 0.1 and 1, the azimuthal wavenumbers m =0,1,2(10¹⁶ g s^-1). All three quantities appear to be equally readily excited. The equatorial symmetry is quadrupolar for the magnetic field and dipolar for the flow field system. The maximal magnetic field strength required to cause the instability is almost independent of the magnetic Prandtl number. With typical white dwarf values, a magnetic amplitude of 10⁵ G is estimated.