A magnetic scaleheight: the effect of toroidal magnetic fields on the thickness of accretion discs

We consider accreting systems in which the central object interacts, via the agency of its magnetic field, with the disc that surrounds it. The disc is turbulent and, so, has a finite effective conductivity. The field sweeps across the face of the disc, thereby forming a current that is directed radially within the disc. In turn, this disc current creates a toroidal field, where the interaction between the disc current and the toroidal field produces a Lorentz force that compresses the disc. We investigate this compression, which creates a magnetic scaleheight of the disc that can be much smaller than the conventional scaleheight. We derive an analytic expression for the magnetic scaleheight and apply it to fully ionized discs.     

A semi-analytic solution for the radial and vertical structure of accretion discs with a magnetic wind

A semi-analytic method is presented for solving for the radial and vertical structures of an accretion disc, with a magnetically channelled wind flowing from its surfaces. Both magnetic and turbulent viscous effects are taken into account, and the essential wind properties are related to the disc structure. The angular momentum removed by the wind plays a major part in driving the inflow through the disc, with photospheric temperatures being sufficient to generate the required wind mass flux. The magnetic field is generated by an αω-dynamo, but the method of solution should have application with other magnetic field sources. Self-consistent disc-wind solutions result, with rms turbulent Mach numbers which are in good agreement with those found in simulations of turbulence generated from magnetic shearing instabilities.     

Photon bubbles in accretion discs

We show that radiation-dominated accretion discs are likely to suffer from a 'photon bubble' instability similar to that described by Arons in the context of accretion on to neutron star polar caps. The instability requires a magnetic field for its existence. In an asymptotic regime appropriate to accretion discs, we find that the overstable modes obey the remarkably simple dispersion relation
ο²=−i gkF ( B , k ).
Here g is the vertical gravitational acceleration, B is the magnetic field, and F is a geometric factor of order unity that depends on the relative orientation of the magnetic field and the wavevector. In the non-linear outcome it seems likely that the instability will enhance vertical energy transport and thereby change the structure of the innermost parts of relativistic accretion discs.     

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.     

Turbulence and its parameterization in accretion discs

Accretion disc turbulence is investigated in the framework of the shearing box approximation. The turbulence is either driven by the magneto‐rotational instability or, in the non‐magnetic case, by an explicit and artificial forcing term in the momentum equation. Unlike the magnetic case, where most of the dissipation occurs in the disc corona, in the forced hydrodynamic case most of the dissipation occurs near the midplane. In the hydrodynamic case evidence is presented for the stochastic excitation of epicycles. When the vertical and radial epicyclic frequencies are different (modeling the properties around rotating black holes), the beat frequency between these two frequencies appear to show up as a peak in the temporal power spectrum in some cases. Finally, the full turbulent resistivity tensor is determined and it is found that, if the turbulence is driven by a forcing term, the signs of its off‐diagonal components are such that this effect would not be capable of dynamo action by the shear–current effect. (© 2005 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

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.     

The dynamics of eccentric accretion discs in superhump systems

We have applied an eccentric accretion disc theory in simplified form to the case of an accretion disc in a binary system, where the disc contains the 3:1 Lindblad resonance. This is relevant to the case of superhumps in SU Ursae Majoris cataclysmic variables and other systems, where it is thought that this resonance leads to growth of eccentricity and a modulation in the light curve due to the interaction of a precessing eccentric disc with tidal stresses. A single differential equation is formulated which describes the propagation, resonant excitation and viscous damping of eccentricity. The theory is first worked out in the simple case of a narrow ring and leads to the conclusion that the eccentricity distribution is locally suppressed by the presence of the resonance, creating a dip in the eccentricity at the resonant radius. Application of this theory to the superhump case confirms this conclusion and produces a more accurate expression for the precession rate of the disc than has been previously accomplished with simple dynamical estimates.     

Magnetic field dragging and the vertical structure of thin accretion discs

The presence of an imposed vertical magnetic field may drastically influence the structure of thin accretion discs. If the field is sufficiently strong, the rotation law can depart from the Keplerian one. We consider the structure of a disc for a given eddy magnetic diffusivity but neglect details of the energy transport. The magnetic field is assumed to be in balance with the internal energy of the accretion flow. The thickness of the disc as well as the turbulent magnetic Prandtl number and the viscosity, α , are the key parameters of our model. The calculations show that the radial velocity can reach the sound speed for a magnetic disc if the thickness is comparable to that of a non-magnetic one. This leads to a strong amplification of the accretion rate for a given surface density. The inclination angle of the magnetic field lines can exceed the critical value 30° (required to launch cold jets) even for a relatively small magnetic Prandtl number of order unity. The toroidal magnetic fields induced at the disc surface are smaller than predicted in previous studies.     

Mass loss from the inner regions of accretion discs due to centrifugally driven magnetic wind flows

Wind flows and collimated jets are believed to be a feature of a range of disc accreting systems. These include active galactic nuclei, T Tauri stars, X-ray binaries and cataclysmic variables. The observed collimation implies large-scale magnetic fields and it is known that dipole-symmetry fields of sufficient strength can channel wind flows emanating from the surfaces of a disc. The disc inflow leads to the bending of the poloidal magnetic field lines, and centrifugally driven magnetic winds can be launched when the bending exceeds a critical value. Such winds can result in angular momentum transport at least as effective as turbulent viscosity, and hence they can play a major part in driving the disc inflow.
It is shown here that if the standard boundary condition of vanishing viscous stress close to the stellar surface is applied, together with the standard connection between viscosity and magnetic diffusivity, then poloidal magnetic field bending increases as the star is approached with a corresponding increase in the wind mass loss rate. A significant amount of material can be lost from the system via the enhanced wind from a narrow region close to the stellar surface. This occurs for a Keplerian angular velocity distribution and for a modified form of angular velocity, which allows for matching of the disc and stellar rotation rates through a boundary layer above the stellar surface. The enhanced mass loss is significantly affected by the behaviour of the disc angular velocity as the stellar surface is approached, and hence by the stellar rotation rate. Such a mechanism may be related to the production of jets from the inner regions of disc accreting systems.     

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.     

Hydrodynamic simulations of propagating warps and bending waves in accretion discs

We present the results of a study of propagating warp or bending waves in accretion discs. Three-dimensional hydrodynamic simulations were performed using smoothed particle hydrodynamics (SPH), and the results are compared with calculations based on the linear theory of warped discs.
We examine the response of a gaseous disc to an initially imposed warping disturbance under a variety of physical conditions. We consider primarily the physical regime in which the dimensionless viscosity parameter α < H r , where H r is the disc aspect ratio, so that bending waves are expected to propagate. We also performed calculations for disc models in which α > H r , where the warps are expected to evolve diffusively. Small-amplitude (linear) perturbations are studied in both Keplerian and slightly non-Keplerian discs, and we find that the results of the SPH calculations can be reasonably well fitted by those of the linear theory. The main results of these calculations are: (i) the warp in Keplerian discs when α < H r propagates with little dispersion, and damps at a rate expected from estimates of the code viscosity; (ii) warps evolve diffusively when α > H r ; (iii) the slightly non-Keplerian discs lead to a substantially more dispersive behaviour of the warps, which damp at a similar rate to the Keplerian case, when α < H r .
Initially imposed higher amplitude, non-linear warping disturbances were studied in Keplerian discs. The results indicate that non-linear warps can lead to the formation of shocks, and that the evolution of the warp becomes less wave-like and more diffusive in character.
This work is relevant to the study of the warped accretion discs that may occur around Kerr black holes or in misaligned binary systems, and is mainly concerned with discs in which α < H r . The results indicate that SPH can model the hydrodynamics of warped discs, even when using rather modest numbers of particles.     

Tidally distorted accretion discs in binary stars

The non-axisymmetric features observed in the discs of dwarf novae in outburst are usually considered to be spiral shocks, which are the non-linear relatives of tidally excited waves. This interpretation suffers from a number of problems. For example, the natural site of wave excitation lies outside the Roche lobe, the disc must be especially hot, and most treatments of wave propagation do not take into account the vertical structure of the disc.
In this paper I construct a detailed semi-analytical model of the non-linear tidal distortion of a thin, three-dimensional accretion disc by a binary companion on a circular orbit. The analysis presented here allows for vertical motion and radiative energy transport, and introduces a simple model for the turbulent magnetic stress. The   m =2  inner vertical resonance has an important influence on the amplitude and phase of the tidal distortion. I show that the observed patterns find a natural explanation if the emission is associated with the tidally thickened sectors of the outer disc, which may be irradiated from the centre. According to this hypothesis, it may be possible to constrain the physical parameters of the disc through future observations.     

Episodic accretion in magnetically layered protoplanetary discs

We study protoplanetary disc evolution assuming that angular momentum transport is driven by gravitational instability at large radii, and magnetohydrodynamic (MHD) turbulence in the hot inner regions. At radii of the order of 1 au such discs develop a magnetically layered structure, with accretion occurring in an ionized surface layer overlying quiescent gas that is too cool to sustain MHD turbulence. We show that layered discs are subject to a limit cycle instability, in which accretion on to the protostar occurs in ∼10⁴-yr bursts with ̇ ∼10⁻⁵ M_⊙ yr⁻¹, separated by quiescent intervals lasting ∼10⁵ yr where ̇ ≈10⁻⁸ M_⊙ yr⁻¹. Such bursts could lead to repeated episodes of strong mass outflow in young stellar objects. The transition to this episodic mode of accretion occurs at an early epoch ( t ≪1 Myr), and the model therefore predicts that many young pre-main-sequence stars should have low rates of accretion through the inner disc. At ages of a few Myr, the discs are up to an order of magnitude more massive than the minimum-mass solar nebula, with most of the mass locked up in the quiescent layer of the disc at r ∼1 au. The predicted rate of low-mass planetary migration is reduced at the outer edge of the layered disc, which could lead to an enhanced probability of giant planet formation at radii of 1–3 au.