Slim accretion discs: a model for ADAF-SLE transitions

We numerically construct slim, global, vertically integrated models of optically thin, transonic accretion discs around black holes, assuming a regularity condition at the sonic radius and boundary conditions at the outer radius of the disc and near the black hole. In agreement with several previous studies, we find two branches of shock-free solutions, in which the cooling is dominated either by advection, or by local radiation. We also confirm that the part of the accretion flow where advection dominates is in some circumstances limited in size: it does not extend beyond a certain outer limiting radius. New results found in our paper concern the location of the limiting radius and properties of the flow near to it. In particular, we find that beyond the limiting radius, the advective dominated solutions match on to Shapiro, Lightman & Eardley (SLE) discs through a smooth transition region. Therefore, the full global solutions are shock-free and unlimited in size. There is no need for postulating an extra physical effect (e.g. evaporation) for triggering the ADAF-SLE transition. It occurs due to standard accretion processes described by the classic slim disc equations.Comment: 12 pages, 7 figures, MNRAS accepte

Convection-Dominated Accretion Flows

Non-radiating, advection-dominated, accretion flows are convectively unstable. We calculate the two-dimensional (r-theta) structure of such flows assuming that (1) convection transports angular momentum inwards, opposite to normal viscosity and (2) viscous transport by other mechanisms (e.g., magnetic fields) is weak (alpha << 1). Under such conditions convection dominates the dynamics of the accretion flow and leads to a steady state structure that is marginally stable to convection. We show that the marginally stable flow has a constant temperature and rotational velocity on spherical shells, a net flux of energy from small to large radii, zero net accretion rate, and a radial density profile proportional to r^{-1/2}, flatter than the r^{-3/2} profile characteristic of spherical accretion flows. This solution accurately describes the full two-dimensional structure of recent axisymmetric numerical simulations of advection-dominated accretion flows.Comment: final version accepted by ApJ; discussion expanded, references adde

The Magnetohydrodynamics of Convection-Dominated Accretion Flows

Radiatively inefficient accretion flows onto black holes are unstable due to both an outwardly decreasing entropy (`convection') and an outwardly decreasing rotation rate (the `magnetorotational instability'; MRI). Using a linear magnetohydrodynamic stability analysis, we show that long-wavelength modes are primarily destabilized by the entropy gradient and that such `convective' modes transport angular momentum inwards. Moreover, the stability criteria for the convective modes are the standard Hoiland criteria of hydrodynamics. By contrast, shorter wavelength modes are primarily destabilized by magnetic tension and differential rotation. These `MRI' modes transport angular momentum outwards. The convection-dominated accretion flow (CDAF) model, which has been proposed for radiatively inefficient accretion onto a black hole, posits that inward angular momentum transport and outward energy transport by long-wavelength convective fluctuations are crucial for determining the structure of the accretion flow. Our analysis suggests that the CDAF model is applicable to a magnetohydrodynamic accretion flow provided the magnetic field saturates at a sufficiently sub-equipartition value (plasma beta >> 1), so that long-wavelength convective fluctuations can fit inside the accretion disk. Numerical magnetohydrodynamic simulations are required to determine whether such a sub-equipartition field is in fact obtained.Comment: 17 pages including 3 figures. Accepted for publication in ApJ. New appendix and figure were added; some changes of the text were made in response to the referee

Statistical theory of thermal instability

A new statistical approach is presented to study the thermal instability process of optically thin unmagnetized plasma. In this approach the time evolution of mass distribution function over temperature is calculated. This function characterizes the statistical properties of the multiphase medium of arbitrary spaced three-dimensional structure of arbitrary temperature perturbations. We construct our theory under the isobarical condition (P=const over space), which is satisfied in the short wavelength limit. The developed theory is illustrated in the case of thermal instability of a slowly expanding interstellar cloud. Numerical solutions of equations of the statistical theory are constucted and compared with hydrodynamical solutions. The results of both approaches are identical in the short wavelength range when the isobarity condition is satisfied. Also the limits of applicability of the statistical theory are estimated. The possible evolution of initial spectrum of perturbations is discussed. The proposed theory and numerical models can be relevant to the formation of the two-phases medium in the ~1pc region around quasars. Then small warm (T~10000K) clouds are formed as the result of thermal instability in an expanded gas fragment, which is a product of either a star-star or star-accretion disk collision.Comment: 11 pages, 8 figures, submitted to MNRA

Two-dimensional models of hydrodynamical accretion flows into black holes

We present a systematic numerical study of two-dimensional axisymmetric accretion flows around black holes. The flows have no radiative cooling and are treated in the framework of the hydrodynamical approximation. The models calculated in this study cover the large range of the relevant parameter space. There are four types of flows, determined by the values of the viscosity parameter $\alpha$ and the adiabatic index $\gamma$: convective flows, large-scale circulations, pure inflows and bipolar outflows. Thermal conduction introduces significant changes to the solutions, but does not create a new flow type. Convective accretion flows and flows with large-scale circulations have significant outward-directed energy fluxes, which have important implications for the spectra and luminosities of accreting black holes.Comment: 43 pages, 23 figures, submitted to Ap

Synchrotron Radiation From Radiatively Inefficient Accretion Flow Simulations: Applications to Sgr A*

We calculate synchrotron radiation in three-dimensional pseudo-Newtonian magnetohydrodynamic simulations of radiatively inefficient accretion flows. We show that the emission is highly variable at optically thin frequencies, with order of magnitude variability on time-scales as short as the orbital period near the last stable orbit; this emission is linearly polarized at the 20-50 % level due to the coherent toroidal magnetic field in the flow. At optically thick frequencies, both the variability amplitude and polarization fraction decrease significantly with decreasing photon frequency. We argue that these results are broadly consistent with the observed properties of Sgr A* at the Galactic Center, including the rapid infrared flaring.Comment: Accepted for publication in Ap

Spectral Models of Convection-Dominated Accretion Flows

For small values of the dimensionless viscosity parameter, namely $\alpha\lesssim 0.1$, the dynamics of non-radiating accretion flows is dominated by convection; convection strongly suppresses the accretion of matter onto the central object and transports a luminosity $\sim 10^{-3}-10^{-2} \dot M c^2$ from small to large radii in the flow. A fraction of this convective luminosity is likely to be radiated at large radii via thermal bremsstrahlung emission. We show that this leads to a correlation between the frequency of maximal bremsstrahlung emission and the luminosity of the source, $\nu_{\rm peak} \propto L^{2/3}$. Accreting black holes with X-ray luminosities $10^{-4} L_{Edd}\gtrsim L_X(0.5-10{\rm keV}) \gtrsim 10^{-7}L_{Edd}$ are expected to have hard X-ray spectra, with photon indices $\Gamma\sim2$, and sources with $L_X\lesssim 10^{-9}L_{Edd}$ are expected to have soft spectra, with $\Gamma\sim3.5$. This is testable with {\it Chandra} and {\it XMM}.Comment: final version accepted by ApJ; significant modifications from previous versio

Three-dimensional MHD Simulations of Radiatively Inefficient Accretion Flows

We present three-dimensional MHD simulations of rotating radiatively inefficient accretion flows onto black holes. In the simulations, we continuously inject magnetized matter into the computational domain near the outer boundary, and we run the calculations long enough for the resulting accretion flow to reach a quasi-steady state. We have studied two limiting cases for the geometry of the injected magnetic field: pure toroidal field and pure poloidal field. In the case of toroidal field injection, the accreting matter forms a nearly axisymmetric, geometrically-thick, turbulent accretion disk. The disk resembles in many respects the convection-dominated accretion flows found in previous numerical and analytical investigations of viscous hydrodynamic flows. Models with poloidal field injection evolve through two distinct phases. In an initial transient phase, the flow forms a relatively flattened, quasi-Keplerian disk with a hot corona and a bipolar outflow. However, when the flow later achieves steady state, it changes in character completely. The magnetized accreting gas becomes two-phase, with most of the volume being dominated by a strong dipolar magnetic field from which a thermal low-density wind flows out. Accretion occurs mainly via narrow slowly-rotating radial streams which `diffuse' through the magnetic field with the help of magnetic reconnection events.Comment: 35 pages including 3 built-in plots and 14 separate jpg-plots; version accepted by Ap