Stability of accretion discs threaded by a strong magnetic field

We study the stability of poloidal magnetic fields anchored in a thin accretion disc. The two-dimensional hydrodynamics in the disc plane is followed by a grid-based numerical simulation including the vertically integrated magnetic forces. The three-dimensional magnetic field outside the disc is calculated in a potential field approximation from the magnetic flux density distribution in the disc. For uniformly rotating discs we confirm numerically the existence of the interchange instability as predicted by Spruit, Stehle & Papaloizou . In agreement with predictions from the shearing sheet model, discs with Keplerian rotation are found to be stabilized by the shear, as long as the contribution of magnetic forces to support against gravity is small. When this support becomes significant, we find a global instability which transports angular momentum outwardly and allows mass to accrete inwardly. The instability takes the form of a m =1 rotating 'crescent', reminiscent of the purely hydrodynamic non-linear instability previously found in pressure-supported discs. A model where the initial surface mass density Σ( r ) and B _z ( r ) decrease with radius as power laws shows transient mass accretion during about six orbital periods, and settles into a state with surface density and field strength decreasing approximately exponentially with radius. We argue that this instability is likely to be the main angular momentum transport mechanism in discs with a poloidal magnetic field sufficiently strong to suppress magnetic turbulence. It may be especially relevant in jet-producing discs.     

The stability of accretion discs with inflow driven purely by magnetic winds

A model is presented for an accretion disc in which the inflow is driven purely by the angular momentum removed in a centrifugally accelerated magnetic wind. Turbulent discs around compact stars are considered, with the required magnetic field being generated in the disc by a simple dynamo. The turbulent magnetic Prandtl number, N _p, measures the ratio of turbulent viscosity to turbulent magnetic diffusivity. Formally, the hypothetical limit   N _p→ 0  corresponds to the magnetic wind torque dominating the viscous torque, but in practice the inflow is magnetically controlled for   N _p≲ 0.1  .
The suggestion by previous authors that purely magnetic wind-driven discs may be unstable is investigated. A detailed steady solution is found which allows perturbations to the thermal balance and vertical equilibrium to be calculated, and hence the effect of perturbations to the magnetic diffusivity, η, to be assessed. For a standard parametrized form of η, the wind-driven angular momentum balance is found to be linearly unstable. An increase in the inflow rate leads to increased bending of the poloidal magnetic field and an enhanced wind mass loss rate. This increases the angular momentum loss rate which drives further inflow. There is a resultant increase in η, due to the temperature perturbation, but this does not relieve field bending sufficiently to prevent the instability.     

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.     

Global structure of self-excited magnetic fields arising from the magnetic shear instability

The global structure of a self-excited magnetic field arising from the magnetic shear instability has been simulated in spherical geometry by a 3D fully non-linear approach. In order to model the structure of an accretion disc we prescribe a rotation profile of the Brandt type which is Keplerian in the outer regions but yields rigid rotation at the inner core. We performed a whole series of runs at different dynamo numbers with an increasing number of modes in spectral space, thereby checking the influence of the numerical resolution in our simulations. Starting from arbitrary small perturbations, the magnetic and kinetic energies grow by several orders of magnitude as soon as a certain azimuthal resolution of at least m =15 was used at a dynamo number of order C =10⁵. Several phases of the transition to turbulence are realized and interpretations are given for the respective effects occurring at each stage. The resulting magnetic field is highly non-axisymmetric and possesses a pronounced inhomogeneous vortex structure of twisted flux tubes. The flow is almost axisymmetric but shows a Kolmogorov-like behaviour for small scales. The outer surface of the shell is penetrated by magnetic field lines in spot-like regions, which are located mainly in the equatorial plane. For very high dynamo numbers we find a cyclic behaviour of the averaged magnetic field amplitude. The problem of angular momentum transport is discussed in terms of the ShakuraSunyaev viscosity alpha , which depends exponentially on the radial distance and adopts values in the range 10³10⁵.     

Hall equilibrium of thin Keplerian discs embedded in mixed poloidal and toroidal magnetic fields

Axisymmetric steady-state weakly ionized Hall–magnetohydrodynamic (MHD) Keplerian thin discs are investigated by using asymptotic expansions in the small disc aspect ratio ε. The model incorporates the azimuthal and poloidal components of the magnetic fields in the leading order in ε. The disc structure is described by an appropriate Grad–Shafranov equation for the poloidal flux function ψ that involves two arbitrary functions of ψ for the toroidal and poloidal currents. The flux function is symmetric about the mid-plane and satisfies certain boundary conditions at the near-horizontal disc edges. The boundary conditions model the combined effect of the primordial as well as the dipole-like magnetic fields. An analytical solution for the Hall equilibrium is achieved by further expanding the relevant equations in an additional small parameter δ that is inversely proportional to the Hall parameter. It is thus found that the Hall equilibrium discs fall into two types: Keplerian discs with (i) small  ( R _d∼δ⁰)  and (ii) large  ( R _d≳δ^{− k} , k > 0)  radius of the disc. The numerical examples that are presented demonstrate the richness and great variety of magnetic and density configurations that may be achieved under the Hall–MHD equilibrium.     

Negative Reynolds stress generation by accretion disc convection

The phenomenon of negative viscosity-alpha in convectively unstable Keplerian accretion discs is discussed. The convection is considered as a random flow with an axisymmetric mesoscale pattern. Its correlation tensor is computed with a time-averaging procedure using Kley's 2D hydrocode. There is a distinct anisotropy between the turbulence intensities in the radial and azimuthal directions, i.e. the radial velocity rms dominates the azimuthal one. As a consequence, an extra term in the expression for the turbulent transport of angular momentum appears which does not vanish for rigid rotation ('Λ-effect'). It is negative ('inwards transport') and even seems to dominate the positive contribution of the eddy viscosity representing outwards transport of angular momentum. For a turbulence model close to that of the mixing-length theory, the rotational influence on the anisotropy of the turbulence intensities,     , and the covariance  〈 u ' _R u ' _φ 〉  – representing the angular momentum transport – is computed and compared with the accretion disc simulations. Indeed, the negative angular momentum transport can be explained with the observed dominance of the radial turbulence intensity. If, on the other hand, in turbulence fields the azimuthal intensity would dominate or the turbulence is even isotropic, then we always find a positive transport of the angular momentum.     

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.     

The response of the Fe Kα line to changes in the X-ray illumination of accretion discs

X-ray reflection spectra from photoionized accretion discs in active galaxies are presented for a wide range of illumination conditions. The energy, equivalent width (EW) and flux of the Fe K α line are shown to depend strongly on the ratio of illuminating flux to disc flux,   F _x/ F _disc  , the photon index of the irradiating power law, Γ, and the incidence angle of the radiation, i . When   F _x/ F _disc≤2  a neutral Fe K α line is prominent for all but the largest values of Γ. At higher illuminating fluxes an He-like Fe K α line at 6.7 keV dominates the line complex. With a high-energy cut-off of 100 keV, the thermal ionization instability seems to suppress the ionized Fe K α line when  Γ≤1.6  . The Fe K α line flux correlates with   F _x/ F _disc  , but the dependence weakens as iron becomes fully ionized. The EW is roughly constant when   F _x/ F _disc  is low and a neutral line dominates, but then declines as the line progresses through higher ionization stages. There is a strong positive correlation between the Fe K α EW and Γ when the line energy is at 6.7 keV, and a slight negative one when it is at 6.4 keV. This is a potential observational diagnostic of the ionization state of the disc. Observations of the broad Fe K α line, which take into account any narrow component, would be able to test these predictions. Ionized Fe K α lines at 6.7 keV are predicted to be common in a simple magnetic flare geometry. A model that includes multiple ionization gradients on the disc is postulated to reconcile the results with observations.     

Gravitational wave emission from a companion black hole in the presence of an accretion disc around a super-massive Kerr black hole

Gravitational wave signal characteristics from a binary black hole system in which the companion moves through the accretion disc of the primary are studied. We chose the primary to be a super-massive  ( M = 10⁸ M_⊙)  Kerr black hole and the companion to be a massive black hole  ( M = 10⁵ M_⊙)  to clearly demonstrate the effects. We show that the drag exerted on the companion by the disc is sufficient to reduce the coalescence time of the binary. The drag is primarily due to the fact that the accretion disc on a black hole deviates from a Keplerian disc and becomes sub-Keplerian due to inner boundary condition on the black hole horizon. We consider two types of accretion rates on to the companion. The companion is deeply immersed inside the disc and it can accrete at the Bondi rate which depends on the instantaneous density of the disc. However, an accretion disc can also form around the smaller black hole and it can accrete at its Eddington rate. Thus, this case is also studied and the results are compared. We find that the effect of the disc will be significant in reducing the coalescence time and one needs to incorporate this while interpreting gravitational wave signals emitted from such a binary system.     

Magnetic field, dust and axisymmetrical mass loss on the asymptotic giant branch

I propose a mechanism for axisymmetrical mass loss on the asymptotic giant branch (AGB) that may account for the axially symmetric structure of elliptical planetary nebulae. The proposed model operates for slowly rotating AGB stars, having angular velocities in the range of 10⁻⁴ω _Kep  ω  10⁻² ω_Kep, where ω_Kep is the equatorial Keplerian angular velocity. Such angular velocities could be gained from a planet companion of mass  0.1  M _Jupiter, which deposits its orbital angular momentum to the envelope at late stages, or even from single stars that are fast rotators on the main sequence. The model assumes that dynamo magnetic activity results in the formation of cool spots, above which dust forms much more easily. The enhanced magnetic activity towards the equator results in a higher dust formation rate there, and hence higher mass-loss rate. As the star ascends the AGB, both the mass-loss rate and magnetic activity increase rapidly, and hence the mass loss becomes more asymmetrical, with higher mass-loss rate closer to the equatorial plane.     

Magnetorotational instability in stratified, weakly ionized accretion discs

We present a linear analysis of the vertical structure and growth of the magnetorotational instability in stratified, weakly ionized accretion discs, such as protostellar and quiescent dwarf novae systems. The method includes the effects of the magnetic coupling, the conductivity regime of the fluid and the strength of the magnetic field, which is initially vertical. The conductivity is treated as a tensor and is assumed to be constant with height.
We obtained solutions for the structure and growth rate of global unstable modes for different conductivity regimes, strengths of the initial magnetic field and coupling between ionized and neutral components of the fluid. The envelopes of short-wavelength perturbations are determined by the action of competing local growth rates at different heights, driven by the vertical stratification of the disc. Ambipolar diffusion perturbations peak consistently higher above the midplane than modes including Hall conductivity. For weak coupling, perturbations including the Hall effect grow faster and act over a more extended cross-section of the disc than those obtained using the ambipolar diffusion approximation.
Finally, we derived an approximate criterion for when Hall diffusion determines the growth of the magnetorotational instability. This is satisfied over a wide range of radii in protostellar discs, reducing the extent of the magnetic 'dead zone'. Even if the magnetic coupling is weak, significant accretion may occur close to the midplane, rather than in the surface regions of weakly ionized discs.     

Simulations of superhumps and superoutbursts

We numerically study the tidal instability of accretion discs in close binary systems using a two-dimensional SPH code. We find that the precession rate of tidally unstable, eccentric discs does not only depend upon the binary mass ratio q . Although the (prograde) disc precession rate increases with the strength of the tidal potential, we find that increasing the shear viscosity ν also has a significant prograde effect. Increasing the disc temperature has a retrograde impact upon the precession rate.   We find that motion relative to the binary potential results in superhump-like, periodic luminosity variations in the outer reaches of an eccentric disc. The nature and location of the luminosity modulation are functions of ν. Light curves most similar to observations are obtained for ν values appropriate for a dwarf nova in outburst.   We investigate the thermal–tidal instability model for superoutburst. A dwarf nova outburst is simulated by instantaneously increasing ν, which causes a rapid radial expansion of the disc. Should the disc encounter the 3: 1 eccentric inner Lindblad resonance and become tidally unstable, then tidal torques become much more efficient at removing angular momentum from the disc. The disc then shrinks and M _d increases. The resulting increase in disc luminosity is found to be consistent with the excess luminosity of a superoutburst.     

Small-scale inviscid accretion discs around black holes

Gas falling quasi-spherically on to a black hole forms an inner accretion disc if its specific angular momentum l exceeds l ∗∼ r _g c , where r _g is the Schwarzschild radius. The standard disc model assumes l ≫ l ∗. We argue that, in many black hole sources, accretion flows have angular momenta just above the threshold for disc formation, l ≳ l ∗, and assess the accretion mechanism in this regime. In a range l ∗< l < l _cr, a small-scale disc forms in which gas spirals fast into the black hole without any help from horizontal viscous stresses. Such an 'inviscid' disc, however, interacts inelastically with the feeding infall. The disc–infall interaction determines the dynamics and luminosity of the accretion flow. The inviscid disc radius can be as large as 14 r _g, and the energy release peaks at 2 r _g. The disc emits a Comptonized X-ray spectrum with a break at ∼100 keV. This accretion regime is likely to take place in wind-fed X-ray binaries and is also possible in active galactic nuclei.     

A model of magnetically induced disc–corona for black hole binaries

We propose a model of magnetic connection (MC) of a black hole with its surrounding accretion disc based on large-scale magnetic field. The MC gives rise to transport of energy and angular momentum between the black hole and the disc, and the closed field lines pipe the hot matter evaporated from the disc, and shape it in the corona above the disc to form a magnetically induced disc–corona system, in which the corona has the same configuration as the large-scale magnetic field. We numerically solve the dynamic equations in the context of the Kerr metric, in which the large-scale magnetic field is determined by dynamo process and equipartition between magnetic pressure and gas pressure. Thus we can obtain a global solution rather than assuming the distribution of large-scale magnetic field beforehand. The main MC effects lie in three aspects. (1) The rotational energy of a fast-spinning black hole can be extracted, enhancing the dissipation in the accretion disc, (2) the closed field lines provide a natural channel for corona matter escaping from disc and finally falling into black hole and (3) the scope of the corona can be bounded by the conservation of magnetic flux. We simulate the high-energy spectra of this system by using Monte Carlo method, and find that the relative hardness of the spectra decreases as accretion rate or black hole spin a _* increases. We fit the typical X-ray spectra of three black hole binaries  (GRO J1655−40, XTE 1118+480 and GX 339−4)  in the low/hard or very high state.     

A possible energy mechanism for cosmological γ-ray bursts

We suggest that an extreme Kerr black hole with a mass ∼10⁶ M_⊙, a dimensionless angular momentum     and a marginally stable orbital radius     located in a normal galaxy may produce a γ -ray burst (GRB) by capturing and disrupting a star. During the capture period, a transient accretion disc is formed and a strong transient magnetic field ∼     lasting for     may be produced at the inner boundary of the accretion disc. A large amount of rotational energy of the black hole is extracted and released in an ultrarelativistic jet with a bulk Lorentz factor Γ larger than 10³ via the Blandford–Znajek process. The relativistic jet energy can be converted into γ -radiation via an internal shock mechanism. The GRB duration should be the same as the lifetime of the strong transient magnetic field. The maximum number of sub-bursts is estimated to be     because the disc material is likely to break into pieces with a size about the thickness of the disc h at the cusp     The shortest risetime of the burst estimated from this model is ∼     The model GRB density rate is also estimated.     

The runaway instability of thick discs around black holes I. The constant angular momentum case

We present results from a numerical study of the runaway instability of thick discs around black holes. This instability is an important issue for most models of cosmic gamma-ray bursts, where the central engine responsible for the initial energy release is such a system consisting of a thick disc surrounding a black hole. We have carried out a comprehensive number of time-dependent simulations aimed at exploring the appearance of the instability. Our study has been performed using a fully relativistic hydrodynamics code. The general relativistic hydrodynamic equations are formulated as a hyperbolic flux-conservative system and solved using a suitable Godunov-type scheme. We build a series of constant angular momentum discs around a Schwarzschild black hole. Furthermore, the self-gravity of the disc is neglected and the evolution of the central black hole is assumed to be that of a sequence of exact Schwarzschild black holes of varying mass. The black hole mass increase is thus determined by the mass accretion rate across the event horizon. In agreement with previous studies based on stationary models, we find that by allowing the mass of the black hole to grow the disc becomes unstable. Our hydrodynamical simulations show that for all disc-to-hole mass ratios considered (between 1 and 0.05), the runaway instability appears very fast on a dynamical time-scale of a few orbital periods, typically a few 10 ms and never exceeding 1 s for our particular choice of the mass of the black hole (2.5 M) and a large range of mass fluxes  ( m 10^-3 M s^-1)  . The implications of our results in the context of gamma-ray bursts are briefly discussed.     

The excitation of spiral density waves through turbulent fluctuations in accretion discs – I. WKBJ theory

We study and elucidate the mechanism of spiral density wave excitation in a differentially rotating flow with turbulence which could result from the magneto-rotational instability. We formulate a set of wave equations with sources that are only non-zero in the presence of turbulent fluctuations. We solve these in a shearing box domain, subject to the boundary conditions of periodicity in shearing coordinates, using a WKBJ method. It is found that, for a particular azimuthal wavelength, the wave excitation occurs through a sequence of regularly spaced swings during which the wave changes from leading to trailing form. This is a generic process that is expected to occur in shearing discs with turbulence. Trailing waves of equal amplitude propagating in opposite directions are produced, both of which produce an outward angular momentum flux that we give expressions for as functions of the disc parameters and azimuthal wavelength.
By solving the wave amplitude equations numerically, we justify the WKBJ approach for a Keplerian rotation law for all parameter regimes of interest. In order to quantify the wave excitation completely, the important wave source terms need to be specified. Assuming conditions of weak non-linearity, these can be identified and are associated with a quantity related to the potential vorticity, being the only survivors in the linear regime. Under the additional assumption that the source has a flat power spectrum at long azimuthal wavelengths, the optimal azimuthal wavelength produced is found to be determined solely by the WKBJ response and is estimated to be  2π H   , with H being the nominal disc scaleheight. In a following paper by Heinemann & Papaloizou, we perform direct three-dimensional simulations and compare results manifesting the wave excitation process and its source with the assumptions made and the theory developed here in detail, finding excellent agreement.