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.     

Waves and instabilities in a differentially rotating disc containing a poloidal magnetic field

The theory of waves and instabilities in a differentially rotating disc containing a poloidal magnetic field is developed within the framework of ideal magnetohydrodynamics. A continuous spectrum, for which the eigenfunctions are localized on individual magnetic surfaces, is identified but is found not to contain any instabilities associated with differential rotation. The normal modes of a weakly magnetized thin disc are studied by extending the asymptotic methods used previously to describe the equilibria. Waves propagate radially in the disc according to a dispersion relation which is determined by solving an eigenvalue problem at each radius. The dispersion relation for a hydrodynamic disc is re-examined and the modes are classified according to their behaviour in the limit of large wavenumber. The addition of a magnetic field introduces new, potentially unstable, modes and also breaks up the dispersion diagram by causing avoided crossings. The stability boundary to the magnetorotational instability in the parameter space of polytropic equilibria is located by solving directly for marginally stable equilibria. For a given vertical magnetic field in the disc, bending of the field lines has a stabilizing effect and it is shown that stable equilibria exist which are capable of launching a predominantly centrifugally driven wind.     

The Balbus–Hawley instability in weakly ionized discs

MHD in protostellar discs is modified by the Hall current when the ambipolar diffusion approximation breaks down. Here I examine the Balbus–Hawley (magnetorotational) instability of a weak, vertical magnetic field within a weakly ionized disc. Vertical stratification is neglected, and a linear analysis is undertaken for the case in which the wavevector of the perturbation is parallel to the magnetic field.
The growth rate depends on whether the initial magnetic field is parallel or antiparallel to the angular momentum of the disc. The parallel case is less (more) unstable than the antiparallel case if the Hall current is dominated by negative (positive) species. The less-unstable orientation is stable for χ ≲0.5, where χ is the ratio of a generalized neutral–ion collision frequency to the Keplerian frequency. The other orientation has a formal growth rate of the order of the Keplerian angular frequency even in the limit χ →0! In this limit the wavelength of the fastest-growing mode tends to infinity, so the minimum level of ionization for instability is determined by the requirement that a wavelength fit within a disc scaleheight. In the ambipolar diffusion case, this requires χ > v _A c _s; in the Hall case this imposes a potentially much weaker limit,         

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.     

Unstable modes of non-axisymmetric gaseous discs

We present a perturbation theory for studying the instabilities of non-axisymmetric gaseous discs. We perturb the dynamical equations of self-gravitating fluids in the vicinity of a non-axisymmetric equilibrium, and expand the perturbed physical quantities in terms of a complete basis set and a small non-axisymmetry parameter ε. We then derive a linear eigenvalue problem in matrix form, and determine the pattern speed, growth rate and mode shapes of the first three unstable modes. In non-axisymmetric discs, the amplitude and the phase angle of travelling waves are functions of both the radius R and the azimuthal angle φ. This is due to the interaction of different wave components in the response spectrum. We demonstrate that wave interaction in unstable discs, with small initial asymmetries, can develop dense clumps during the phase of exponential growth. Local clumps, which occur on the major spiral arms, can constitute seeds of gas giant planets in accretion 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.     

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.     

Thermal conduction between an accretion disc and a corona in active galactic nuclei: vertical structure of the transition layer

We study the vertical structure of the transition layer between an accretion disc and a corona in the context of the existence of a two-phase medium in thermally unstable regions. The disc is illuminated by hard X-ray radiation, and satisfies the condition of hydrostatic equilibrium. We take into account the energy exchange between the hot, Compton-heated corona (∼10⁸ K) and cool disc (∼10⁴ K) arising from both radiative processes and thermal conduction. In the case including thermal conduction, we perform a local stability analysis, and conclude that thermal conduction does not suppress thermal instability. In spite of the continuous temperature profile T ( τ ) there are regions of strong temperature gradient, in which spontaneous perturbations can lead to cloud condensation in the transition layer. We determine the minimum size λ _TC of such a perturbation.     

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.     

Detecting ionized bubbles in redshifted 21-cm maps

Linear transient phenomena induced by flow non-normality in thin self-gravitating astrophysical discs are studied using the shearing sheet approximation. The considered system includes two modes of perturbations: vortex and (spiral density) wave. It is shown that self-gravity considerably alters the vortex mode dynamics; its transient (swing) growth may be several orders of magnitude stronger than in the non-self-gravitating case and two to three times larger than the transient growth of the wave mode. Based on this finding, we comment on the role of vortex mode perturbations in a gravitoturbulent state. We also describe the linear coupling of the perturbation modes, caused by the differential character of disc rotation. The coupling is asymmetric: vortex mode perturbations are able to excite wave mode perturbations, but not vice versa. This asymmetric coupling lends additional significance to the vortex mode as a participant in spiral density waves and shock manifestations in astrophysical discs.     

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.     

On the migration of a system of protoplanets

The evolution of a system consisting of a protoplanetary disc with two embedded Jupiter-sized planets is studied numerically. The disc is assumed to be flat and non-self-gravitating; this is modelled by the planar (two-dimensional) Navier–Stokes equations. The mutual gravitational interaction of the planets and the star, and the gravitational torques of the disc acting on the planets and the central star are included. The planets have an initial mass of one Jupiter mass M _Jup each, and the radial distances from the star are one and two semimajor axes of Jupiter, respectively.
During the evolution a joint wide annular gap is created by the planets. Both planets increase their mass owing to accretion of gas from the disc: after about 2500 orbital periods of the inner planet it has reached a mass of 2.3  M _Jup, while the outer planet has reached a mass of 3.2  M _Jup. The net gravitational torques exerted by the disc on the planets result in an inward migration of the outer planet on time-scales comparable to the viscous evolution time of the disc. The semimajor axis of the inner planet remains constant as there is very little gas left in its vicinity to induce any migration. When the distance of close approach eventually becomes smaller than the mutual Hill radius, the eccentricities increase strongly and the system may become unstable.
If disc depletion occurs rapidly enough before the planets come too close to each other, a stable system similar to our own Solar system may remain. Otherwise the orbits may become unstable and produce systems like υ And.     

Disc structure around strongly magnetic accretors: a full disc solution with turbulent diffusivity

The problem of the effect of a strongly magnetic star on a surrounding accretion disc is considered. For stellar rotation periods greater than a critical value, a numerical solution is found for a steady disc with turbulent magnetic diffusion, including electron scattering opacity and radiation pressure. Inside the corotation radius, the extraction of disc angular momentum by magnetic coupling to the star becomes strong and this leads to enhanced viscous stress and dissipation. The resulting elevated temperature causes electron scattering opacity and radiation pressure to become significant further from the star than in the absence of its magnetic field. The disc ends as its height increases rapidly due to the large central pressure, its density decreases and magnetically induced viscous instability occurs.     

Linearized perturbation on stationary inflow solutions in an inviscid and thin accretion disc

The influence of a linearized perturbation on stationary inflow solutions in an inviscid and thin accretion disc has been studied here, and it has been argued that a perturbative technique would indicate that all possible classes of inflow solutions would be stable. The choice of the driving potential, Newtonian or pseudo-Newtonian, would not particularly affect the arguments which establish the stability of solutions. It has then been surmised that in the matter of the selection of a particular solution, adoption of a non-perturbative technique, based on a more physical criterion, as in the case of the selection of the transonic solution in spherically symmetric accretion, would give a more conclusive indication concerning the choice of a particular branch of the flow.     

The photoevaporation of discs around young stars in massive clusters

We present models in which the photoevaporation of discs around young stars by an external ultraviolet source (as computed by Adams et al.) is coupled with the internal viscous evolution of the discs. These models are applied to the case of the Orion Nebula Cluster (ONC), where the presence of a strong ultraviolet field from the central OB stars, together with a detailed census of circumstellar discs and photoevaporative flows, is well established. In particular we investigate the constraints that are placed on the initial disc properties in the ONC by the twin requirement that most stars possess a disc on a scale of a few astronomical unit (au), but that only a minority (<20 per cent) are resolved by Hubble Space Telescope ( HST ) at a scale of 50 au. We find that these requirements place very weak constraints on the initial radius distribution of circumstellar discs: the resulting size distribution readily forgets the initial radius distribution, owing to the strong positive dependence of the photoevaporation rate on disc radius. Instead, the scarcity of large discs reflects the relative scarcity of initially massive discs (with mass  >0.1 M_⊙  ). The ubiquity of discs on a small scale, on the other hand, mainly constrains the time-span over which the discs have been exposed to the ultraviolet field (<2 Myr). We argue that the discs that are resolved by HST represent a population of discs in which self-gravity was important at the time that the dominant central OB star switched on, but that, according to our models, self-gravity is unlikely to be important in these discs at the present time. We discuss the implications of our results for the so-called proplyd lifetime problem.     

Light curve morphology study of UW CrB – evidence for a 5 d superorbital period

Since its discovery in 1990, UW CrB (also known as MS1603+2600) has remained a peculiar source without firm classification. Our current understanding is that it is an accretion disc corona (ADC) low-mass X-ray binary. In this paper, we present results from our photometric campaign dedicated to studying the changing morphology of the optical light curves. We find that the optical light curves show remarkable evidence for strongly evolving light curve shapes. In addition, we find that these changes show a modulation at a period of ∼5 d. We interpret these changes as either due to strong periodic accretion disc warping or due to other geometrical changes because of disc precession at a period of 5 d. Finally, we have detected 11 new optical bursts, the phase distribution of which supports the idea of a vertically extended asymmetric accretion disc.     

Launching of accretion disc winds along dynamo-generated magnetic fields

The problem of magnetic field generation and advection in accretion discs is considered, in the context of wind launching and angular momentum extraction. A dipole-symmetry solution of the dynamo equations is found, with force-free boundary conditions appropriate for matching to a wind solution. Consideration of the curved field geometry and diffusive nature of the disc enables the position of the sonic point to be calculated and related to the field inclination at the disc surface. A critical inclination of 20° to the horizontal results, for which the sonic point lies in the disc surface and there is no potential barrier to wind launching. Hence the wind mass-loss rate will only become excessive, leading to disc disruption, for large field bending. The compressional effect of the horizontal magnetic field enhances the wind mass flux.     

Non-linear fluid dynamics of eccentric discs

A new theory of eccentric accretion discs is presented. Starting from the basic fluid-dynamical equations in three dimensions, I derive the fundamental set of one-dimensional equations that describe how the mass, angular momentum and eccentricity vector of a thin disc evolve as a result of internal stresses and external forcing. The analysis is asymptotically exact in the limit of a thin disc, and allows for slowly varying eccentricities of arbitrary magnitude. The theory is worked out in detail for a Maxwellian viscoelastic model of the turbulent stress in an accretion disc. This generalizes the conventional alpha viscosity model to account for the non-zero relaxation time of the turbulence, and is physically motivated by a consideration of the nature of magnetohydrodynamic turbulence. It is confirmed that circular discs are typically viscously unstable to eccentric perturbations, as found by Lyubarskij, Postnov & Prokhorov, if the conventional alpha viscosity model is adopted. However, the instability can usually be suppressed by introducing a sufficient relaxation time and/or bulk viscosity. It is then shown that an initially uniformly eccentric disc does not retain its eccentricity as had been suggested by previous analyses. The evolutionary equations should be useful in many applications, including understanding the origin of planetary eccentricities and testing theories of quasi-periodic oscillations in X-ray binaries.     

The excitation, propagation and dissipation of waves in accretion discs: the non-linear axisymmetric case

We analyse the non-linear propagation and dissipation of axisymmetric waves in accretion discs using the ZEUS-2D hydrodynamics code. The waves are numerically resolved in the vertical and radial directions. Both vertically isothermal and thermally stratified accretion discs are considered. The waves are generated by means of resonant forcing, and several forms of forcing are considered. Compressional motions are taken to be locally adiabatic  ( γ =5/3)  . Prior to non-linear dissipation, the numerical results are in excellent agreement with the linear theory of wave channelling in predicting the types of modes that are excited, the energy flux by carried by each mode, and the vertical wave energy distribution as a function of radius. In all cases, waves are excited that propagate on both sides of the resonance (inwards and outwards). For vertically isothermal discs, non-linear dissipation occurs primarily through shocks that result from the classical steepening of acoustic waves. For discs that are substantially thermally stratified, wave channelling is the primary mechanism for shock generation. Wave channelling boosts the Mach number of the wave by vertically confining the wave to a small cool region at the base of the disc atmosphere. In general, outwardly propagating waves with Mach numbers near resonance  ℳ_r≳0.01  undergo shocks within a distance of order the resonance radius.