Magnetohydrodynamical non-radiative accretion flows in two dimensions

We present the results of axisymmetric, time-dependent magnetohydrodynamic simulations of accretion flows around black holes. The calculations begin from a rotationally supported thick torus which contains a weak poloidal field. Accretion is produced by growth and saturation of the magnetorotational instability (MRI) provided that the wavelength of the fastest growing mode is less than the thickness of the torus. Using a computational grid that spans more than two decades in radius, we compare the time-averaged properties of the flow with previous hydrodynamical simulations. The net mass accretion rate is small compared with the mass inflow and outflow rates at large radii associated with turbulent eddies. Turbulence is driven by the MRI rather than convection. The two-dimensional structure of the time-averaged flow is significantly different compared with the hydrodynamical case. We discuss the limitations imposed on our results by the assumption of axisymmetry and the relatively small radial domain.     

A steady-state solution for warped accretion discs

We consider a thin accretion disc warped due to the Bardeen–Petterson effect, presenting both analytical and numerical solutions for the situation in which the two viscosity coefficients vary with radius as a power law, with the two power-law indices not necessarily equal. The analytical solutions are compared with numerical ones, showing that our new analytical solution is more accurate than the previous one, which overestimated the inclination change in the outer disc. Our new analytical solution is appropriate for moderately warped discs, while for extremely misaligned discs only a numerical solution is appropriate.     

Magnetized tori around Kerr black holes: analytic solutions with a toroidal magnetic field

The dynamics of accretion discs around galactic and extragalactic black holes may be influenced by their magnetic field. In this paper, we generalize the fully relativistic theory of stationary axisymmetric tori in Kerr metric of Abramowicz, Jaroszynski & Sikora by including strong toroidal magnetic field and construct analytic solutions for barotropic tori with constant angular momentum. This development is particularly important for the general relativistic computational magnetohydrodynamics that suffers from the lack of exact analytic solutions that are needed to test computer codes.     

The evolution of misaligned accretion discs and spinning black holes

In this paper, we consider the process of alignment of a spinning black hole and a surrounding misaligned accretion disc. We use a simplified set of equations, that describe the evolution of the system in the case where the propagation of warping disturbances in the accretion disc occurs diffusively, a situation likely to be common in the thin discs in active galactic nuclei (AGN). We also allow the direction of the hole spin to move under the action of the disc torques. In such a way, the evolution of the hole–disc system is computed self-consistently. We consider a number of different situations and we explore the relevant parameter range, by varying the location of the warp radius R _w and the propagation speed of the warp. We find that the dissipation associated with the twisting of the disc results in a large increase in the accretion rate through the disc, so that AGN accreting from a misaligned disc are likely to be significantly more luminous than those accreting from a flat disc. We compute explicitly the time-scales for the warping of the disc and for the alignment process and compare our results with earlier estimates based on simplified steady-state solutions. We also confirm earlier predictions that, under appropriate circumstances, accretion can proceed in a counter-aligned fashion, so that the accreted material will spin-down the hole, rather than spinning it up. Our results have implication in a number of different observational features of AGN such as the orientation and shape of jets, the shape of X-ray iron lines and the possibility of obscuration and absorption of X-ray by the outer disc as well as the general issue of the spin history of black holes during their growth.     

A comparative analysis of standard accretion discs spectra: an application to ultraluminous X-ray sources

We compare standard models of accretion discs around black holes (BHs) that include the appropriate zero-torque inner boundary condition and relativistic effects on the emission and propagation of radiation. The comparison is performed adopting the multicolour disc blackbody model (MCD) as reference and looking for the parameter space in which it is in statistical agreement with 'more physical' accretion disc models. We find simple 'recipes' that can be used for adjusting the estimates of the physical inner radius of the disc, the BH mass and the accretion rate inferred using the parameters of the MCD fits. We applied these results to four ultraluminous X-ray sources for which MCD spectral fits of their X-ray soft spectral components have been published and find that, in three cases (NGC 1313 X-1, X-2 and M 81 X-9), the BH masses inferred for a standard disc around a Schwarzschild BH are in the interval  ∼100–200 M_⊙  . Only if the BH is maximally rotating are the masses comparable to the much larger values previously derived in the literature.     

Self-similar solution of hot accretion flows with ordered magnetic field and outflow

Observations and numerical magnetohydrodynamic (MHD) simulations indicate the existence of outflows and ordered large-scale magnetic fields in the inner region of hot accretion flows. In this paper, we present the self-similar solutions for advection-dominated accretion flows (ADAFs) with outflows and ordered magnetic fields. Stimulated by numerical simulations, we assume that the magnetic field has a strong toroidal component and a vertical component in addition to a stochastic component. We obtain the self-similar solutions to the equations describing the magnetized ADAFs, taking into account the dynamical effects of the outflow. We compare the results with the canonical ADAFs and find that the dynamical properties of ADAFs such as radial velocity, angular velocity and temperature can be significantly changed in the presence of ordered magnetic fields and outflows. The stronger the magnetic field is, the lower the temperature of the accretion flow will be and the faster the flow rotates. The relevance to observations is briefly discussed.     

Magnetized, mass-loaded, rotating accretion flows

We present a semi-analytical investigation of a simple one-dimensional, steady-state model for a mass-loaded, rotating, magnetized, hydrodynamical flow. Our approach is analogous to one used in early studies of magnetized winds. The model represents the infall towards a central point mass of the gas generated in a cluster of stars surrounding it, as is likely to occur in some active nuclei and starburst galaxies. We describe the properties of the different classes of infall solutions. We find that the flow becomes faster than the fast-mode speed, and hence decoupled from the centre, only for a limited range of parameter values, and when magnetic stresses are ineffective. Such flow is slowed as it approaches a centrifugal barrier, implying the existence of an accretion disc. When the flow does not become super-fast and the magnetic torque is insufficient, no steady solution extending inward to the centre exists. Finally, with a larger magnetic torque, solutions representing steady sub-Alfvénic flows are found, which can resemble spherical hydrodynamical infall. Such solutions, if applicable, would imply that rotation is not important and that any accretion disc formed would be of very limited size.     

Luminous hot accretion discs

We find a new two-temperature hot branch of equilibrium solutions for stationary accretion discs around black holes. In units of Eddington accretion rate defined as 10 L _Edd c ², the accretion rates to which these solutions correspond are within the range ̇ ₁≲ ̇ ≲1, where ̇ ₁ is the critical rate of advection-dominated accretion flow (ADAF). In these solutions, the energy loss rate of the ions by Coulomb energy transfer between the ions and electrons is larger than the viscously heating rate and it is the advective heating together with the viscous dissipation that balances the Coulomb cooling of ions. When ̇ ₁≲ ̇ ≲ ̇ ₂, where ̇ ₂∼5 ̇ ₁<1, the accretion flow remains hot throughout the disc. When ̇ ₂≲ ̇ ≲1, Coulomb interaction will cool the inner region of the disc within a certain radius ( r _tr∼several tens of Schwarzschild radii or larger depending on the accretion rate and the outer boundary condition) and the disc will collapse on to the equatorial plane and form an optically thick cold annulus. Compared with ADAF, these hot solutions are much more luminous because of the high accretion rate and efficiency; therefore, we call them luminous hot accretion discs.     

Hydrodynamical non-radiative accretion flows in two dimensions

Two-dimensional (axially symmetric) numerical hydrodynamical calculations of accretion flows that cannot cool through emission of radiation are presented. The calculations begin from an equilibrium configuration consisting of a thick torus with constant specific angular momentum. Accretion is induced by the addition of a small anomalous azimuthal shear stress which is characterized by a function ν . We study the flows generated as the amplitude and form of ν are varied. A spherical polar grid which spans more than two orders of magnitude in radius is used to resolve the flow over a wide range of spatial scales. We find that convection in the inner regions produces significant outward mass motions that carry away both the energy liberated by and a large fraction of the mass participating in the accretion flow. Although the instantaneous structure of the flow is complex and dominated by convective eddies, long-time averages of the dynamical variables show remarkable correspondence to certain steady-state solutions. The two-dimensional structure of the time-averaged flow is marginally stable to the Høiland criterion, indicating that convection is efficient. Near the equatorial plane, the radial profiles of the time-averaged variables are power laws with an index that depends on the radial scaling of the shear stress. A stress in which ν ∝ r ^1/2 recovers the widely studied self-similar solution corresponding to an ' α -disc'. We find that, regardless of the adiabatic index of the gas, or the form or magnitude of the shear stress, the mass inflow rate is a strongly increasing function of radius, and is everywhere nearly exactly balanced by mass outflow. The net mass accretion rate through the disc is only a fraction of the rate at which mass is supplied to the inflow at large radii, and is given by the local, viscous accretion rate associated with the flow properties near the central object.