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.     

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.     

Transfer of energy and angular momentum in the magnetic coupling between a rotating black hole and the surrounding accretion disc

The transfer of energy and angular momentum in the magnetic coupling (MC) of a rotating black hole (BH) with its surrounding accretion disc is discussed based on a mapping relation derived by considering the conservation of magnetic flux with two basic assumptions: (i) the magnetic field on the horizon is constant, (ii) the magnetic field on the disc surface varies as a power law with the radial coordinate of the disc. The following results are obtained: (i) the transfer direction of energy and angular momentum between the BH and the disc depends on the position of a co-rotation radius relative to the MC region on the disc, which is eventually determined by the BH spin; (ii) the evolution characteristics of a rotating BH in the MC process without disc accretion are depicted in a parameter space, and a series of values of the BH spin are given to indicate the evolution characteristics; (iii) the efficiency of converting accreted mass into radiation energy of a BH–disc system is discussed by considering the coexistence of disc accretion and the MC process; (iv) the MC effects on disc radiation and the emissivity index are discussed and it is concluded that they are consistent with the recent XMM–Newton observation of the nearby bright Seyfert 1 galaxy MCG–6-30-15 with reference to a variety of parameters of the BH–disc system.     

Numerical simulations of astrophysical jets from Keplerian discs – III. The effects of mass loading

We present 2.5D time-dependent simulations of the non-linear evolution of non-relativistic outflows from the surface of Keplerian accretion discs. The gas is accelerated from the surface of the disc (which is a fixed platform in these simulations) into a cold corona in stable hydrostatic equilibrium. We explore the dependence of the resulting jet characteristics upon the mass loading of the winds. Two initial configurations of the threading disc magnetic field are studied: a potential field and a uniform vertical field configuration.
We show that the nature of the resulting highly collimated, jet-like outflows (steady or episodic) is determined by the mass load of the disc wind. The mass load controls the interplay between the collimating effects of the toroidal field and the kinetic energy density in the outflow. In this regard, we demonstrate that the onset of episodic behaviour of jets appears to be determined by the quantity     which compares the speed for a toroidal Alfvén wave to cross the diameter of the jet, with the flow speed v _p along the jet. This quantity decreases with increasing load. For sufficiently large N (small mass loads), disturbances appear to grow leading to instabilities and shocks. Knots are then generated and the outflow becomes episodic. These effects are qualitatively independent of the initial magnetic configuration that we employed and are probably generic to a wide variety of magnetized accretion disc models.     

Magnetic and vertical shear instabilities in accretion discs

The stability properties of magnetized discs rotating with angular velocity Ω = Ω( s ,  z ), dependent on both the radial and the vertical coordinates s and z , are considered. Such a rotation law is adequate for many astrophysical discs (e.g., galactic and protoplanetary discs, as well as accretion discs in binaries). In general, the angular velocity depends on height, even in thin accretion discs. A linear stability analysis is performed in the Boussinesq approximation, and the dispersion relation is obtained for short-wavelength perturbations. Any dependence of Ω on z can destabilize the flow. This concerns primarily small-scale perturbations for which the stabilizing effect of buoyancy is strongly suppressed due to the energy exchange with the surrounding plasma. For a weak magnetic field, instability of discs is mainly associated with vertical shear, whilst for an intermediate magnetic field the magnetic shear instability, first considered by Chandrasekhar and Velikhov, is more efficient. This instability is caused by the radial shear which is typically much stronger than the vertical shear. Therefore the growth time for the magnetic shear instability is much shorter than for the vertical shear instability. A relatively strong magnetic field can suppress both these instabilities. The vertical shear instability could be the source of turbulence in protoplanetary discs, where the conductivity is low.     

Effect of rotation and self-gravity on the propagation of MHD waves

Magnetohydrodynamics waves and instabilities in rotating, self-gravitating, anisotropic and collision-less plasma were investigated. The general dispersion relation was obtained using standard mode analysis by constructing the linearized set of equations. The wave mode solutions and stability properties of the dispersion relations are discussed in the propagations transverse and parallel to the magnetic field. These special cases are discussed considering the axis of rotation to be in transverse and along the magnetic field. In the case of propagation transverse to the magnetic field with axis of rotation parallel to the magnetic field, we derived the dispersion relation modified by rotation and self-gravitation. In the case of propagation parallel to the magnetic field with axis of rotation perpendicular to the magnetic field, we obtained two separate modes affected by rotation and self-gravitation. This indicates that the Slow mode and fire hose instability are not affected by rotation. Numerical analysis was performed for oblique propagation to show the effect of rotation and self-gravitation. It is found that rotation has an effect of reducing the value of the phase speeds on the fast and Alfven wave modes, but self-gravitation affect only on the Slow modes, thereby reducing the phase speed compare to the ideal magneto hydrodynamic (MHD) case.     

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.     

Viscous Damping of Non-Adiabatic MHD Waves in an Unbounded Solar Coronal Plasma

The damping of MHD waves in solar coronal magnetic field is studied taking into account thermal conduction and compressive viscosity as dissipative mechanisms. We consider viscous homogeneous unbounded solar coronal plasma permeated by a uniform magnetic field. A general fifth-order dispersion relation for MHD waves has been derived and solved numerically for different solar coronal regimes. The dispersion relation results three wave modes: slow, fast, and thermal modes. Damping time and damping per periods for slow- and fast-mode waves determined from dispersion relation show that the slow-mode waves are heavily damped in comparison with fast-mode waves in prominences, prominence–corona transition regions (PCTR), and corona. In PCTRs and coronal active regions, wave instabilities appear for considered heating mechanisms. For same heating mechanisms in different prominences the behavior of damping time and damping per period changes significantly from small to large wavenumbers. In all PCTRs and corona, damping time always decreases linearly with increase in wavenumber indicate sharp damping of slow- and fast-mode waves.     

修正的粘滞律及有磁吸积盘的稳定性研究

本文采用修正的粘滞定律及磁流体力学研究了薄吸积盘内区及外区的稳定性问题。运用微扰方法导出了色散方程，分析了四种情况下吸积盘的不稳定性，结果表明：在同时考虑磁场和修正的粘滞律时，吸积盘中存在着三种振荡模式，其中粘滞模式总是稳定的，磁声速模式（包括向里、向外传播两种模式）通常是不稳定的。这些结果为解释ＢＬ Ｌａｃ天体、Ｓｅｙｆｅｒｔ星系、类星体等活动星系核的光变现象提供了理论依据。     

Non-linear bending waves in Keplerian accretion discs

The non-linear dynamics of a warped accretion disc is investigated in the important case of a thin Keplerian disc with negligible viscosity and self-gravity. A one-dimensional evolutionary equation is formally derived that describes the primary non-linear and dispersive effects on propagating bending waves other than parametric instabilities. It has the form of a derivative non-linear Schrödinger (DNLS) equation with coefficients that are obtained explicitly for a particular model of a disc. The properties of this equation are analysed in some detail and illustrative numerical solutions are presented. The non-linear and dispersive effects both depend on the compressibility of the gas through its adiabatic index Γ. In the physically realistic case Γ < 3, non-linearity does not lead to the steepening of bending waves but instead enhances their linear dispersion. In the opposite case Γ > 3, non-linearity leads to wave steepening and solitary waves are supported. The effects of a small effective viscosity, which may suppress parametric instabilities, are also considered. This analysis may provide a useful point of comparison between theory and numerical simulations of warped accretion discs.     

A comparison of DA white dwarf temperatures and gravities from Lyman and Balmer line studies

The stability of turbulent accretion discs is considered, in which a magnetically influenced wind plays a major role in driving the inflow. The magnetic field is generated by a dynamo operating in the disc, involving radial shear and turbulence. The steady angular momentum balance is found to be linearly stable for a range of radial boundary conditions, and an expression is derived for the adjustment time-scale as a function of the equilibrium ratio of the magnetic and viscous disc torques.     

Drag-driven instability of a dust layer in a magnetized protoplanetary disc

We study drag-driven instability in a protoplanetary disc consisting of a layer of single-sized dust particles which are coupled to the magnetized gas aerodynamically and the particle-to-gas feedback is included. We find a dispersion relation for axisymmetric linear disturbances and growth rate of the unstable modes are calculated numerically. While the secular gravitational instability in the absence of particle-togas feedback predicts the dust layer is unstable, magnetic fields significantly amplify the instability if the Toomre parameter for the gas component is fixed. We also show that even a weak magnetic field is able to amplify the instability more or less irrespective of the dust-gas coupling.     

Waves in cosmic magnetic structures taking into account the anisotropic plasma pressure

This paper provides an analysis of magneto-sonic eigenwaves travelling in magnetic plasma structures based on the Chew-Goldberger-Low approximation, for which the plasma kinetic pressure is different along and across the magnetic field. The anisotropy does not lead to the emergence of new modes. The dependence of phase velocities of the waves, trapped by a single magnetic surface, on the pressure anisotropy is investigated. For a magnetic slab with field-free surroundings, the dispersion relations for the eigenwaves are obtained. The pressure anisotropy may change dispersion relations of such modes significantly. In particular, backward waves are possible in the case of strong anisotropy. The dependences of the thresholds for the mirror and hose instabilities on the system parameters are obtained. In particular, hose and mirror instabilities of such waves are absent for some wave number regions. The results are used to obtain the eigenwave characteristics in coronal loops and chromospheric flux tubes.     

Magnetohydrodynamic instabilities and turbulence in accretion disks

In this paper we review the possibilities for magnetohydrodynamic processes to handle the angular momentum transport in accretion disks. Traditionally the angular momentum transport has been considered to be the result of turbulent viscosity in the disk, although the Keplerian flow in accretion disks is linearly stable towards hydrodynamic perturbations. It is on the other hand linearly unstable to some magnetohydrodynamic (MHD) instabilities. The most important instabilities are the Parker and Balbus-Hawley instabilities that are related to the magnetic buoyancy and the shear flow, respectively. We discuss these instabilities not only in the traditional MHD framework, but also in the context of slender flux tubes, that reduce the complexity of the problem while keeping most of the stability properties of the complete problem. In the non-linear regime the instabilities produce turbulence. Recent numerical simulations describe the generation of magnetic fields by a dynamo in the resulting turbulent flow. Eventually such a dynamo may generate a global magnetic field in the disk. The relation of the MHD-turbulence to observations of accretion disks is still obscure. It is commonly believed that magnetic fields can be highly efficient in transporting the angular momentum, but emission lines, short-time scale variability and non-thermal radiation, which a stellar astronomer would take as signs of magnetic variability, are more commonly observed during periods of low accretion rates. Received October 12, 1995 / Accepted November 16, 1995     

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 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.     

Burial of the polar magnetic field of an accreting neutron star – II. Hydromagnetic stability of axisymmetric equilibria

The theory of polar magnetic burial in accreting neutron stars predicts that a mountain of accreted material accumulates at the magnetic poles of the star, and that, as the mountain spreads equatorward, it is confined by, and compresses, the equatorial magnetic field. Here, we extend previous, axisymmetric, Grad–Shafranov calculations of the hydromagnetic structure of a magnetic mountain up to accreted masses as high as   M _a= 6 × 10⁻⁴ M_⊙  , by importing the output from previous calculations (which were limited by numerical problems and the formation of closed bubbles to   M _a < 10⁻⁴ M_⊙  ) into the time-dependent, ideal-magnetohydrodynamic code zeus-3d and loading additional mass on to the star dynamically. The rise of buoyant magnetic bubbles through the accreted layer is observed in these experiments. We also investigate the stability of the resulting hydromagnetic equilibria by perturbing them in zeus-3d . Surprisingly, it is observed that the equilibria are marginally stable for all   M _a≤ 6 × 10⁻⁴ M_⊙  ; the mountain oscillates persistently when perturbed, in a combination of Alfvén and acoustic modes, without appreciable damping or growth, and is therefore not disrupted (apart from a transient Parker instability initially, which expels <1 per cent of the mass and magnetic flux).     

The radial disc structure around a magnetic neutron star: analytic and semi-analytic solutions

The radial structure of a thin accretion disc is calculated in the presence of a central dipole magnetic field aligned with the rotation axis. The problem is treated using a modified expression for the turbulent magnetic diffusion, which allows the angular momentum equation to be integrated analytically. The governing algebraic equations are solved iteratively between 1 and 10⁴ stellar radii. An analytic approximation is provided that is valid near the disruption radius at about 100 stellar radii. At that point, which is approximately 60 per cent of the Alfvén radius and typically about 30 per cent of the corotation radius, the disc becomes viscously unstable. This instability results from the fact that both radiation pressure and opacity caused by electron scattering become important. This in turn is a consequence of the magnetic field which leads to an enhanced temperature in the inner parts. This is because the magnetic field gives rise to a strongly enhanced vertically integrated viscosity, so that the viscous torque can balance the magnetic torque.     

Interface modes and their instabilities in accretion disc boundary layers

We study global non-axisymmetric oscillation modes trapped near the inner boundary of an accretion disc. Observations indicate that some of the quasi-periodic oscillations (QPOs) observed in the luminosities of accreting compact objects (neutron stars, black holes and white dwarfs) are produced in the innermost regions of accretion discs or boundary layers. Two simple models are considered in this paper. The magnetosphere–disc model consists of a thin Keplerian disc in contact with a uniformly rotating magnetosphere with and low plasma density, while the star–disc model involves a Keplerian disc terminated at the stellar atmosphere with high density and small density scaleheight. We find that the interface modes at the magnetosphere–disc boundary are generally unstable due to Rayleigh–Taylor and/or Kelvin–Helmholtz instabilities. However, differential rotation of the disc tends to suppress Rayleigh–Taylor instability, and a sufficiently high disc sound speed (or temperature) is needed to overcome this suppression and to attain net mode growth. On the other hand, Kelvin–Helmholtz instability may be active at low disc sound speeds. We also find that the interface modes trapped at the boundary between a thin disc and an unmagnetized star do not suffer Rayleigh–Taylor or Kelvin–Helmholtz instability, but can become unstable due to wave leakage to large disc radii and, for sufficiently steep disc density distributions, due to wave absorption at the corotation resonance in the disc. The non-axisymmetric interface modes studied in this paper may be relevant to the high-frequency QPOs observed in some X-ray binaries and in cataclysmic variables.