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.     

Mass transfer in tidally unstable compact binaries

The 2001 outburst of WZ Sagittae has shown the most compelling evidence yet for an enhancement of the mass-transfer rate from the donor star during a dwarf nova outburst in the form of hotspot brightening. I show that, even in this extreme case, the brightening can be attributed to tidal heating near the interaction point of an accretion stream with the expanding edge of an eccentric accretion disc, with no need at all for an increase in the mass-transfer rate. Furthermore, I confirm previous suggestions that an increase in mass-transfer rate through the stream damps any eccentricity in an accretion disc and suppresses the appearance of superhumps, in contradiction to observations. Tidal heating is expected to be most significant in systems with small mass ratios. It follows that systems like WZ Sagittae – which has a tiny mass ratio – are those most likely to show a brightening in the hotspot region.     

Superhumps: confronting theory with observation

We review the theory and observations related to the 'superhump' precession of eccentric accretion discs in close binary systems. We agree with earlier work, although for different reasons, that the discrepancy between observation and dynamical theory implies that the effect of pressure in the disc cannot be neglected. We extend earlier work that investigates this effect to include the correct expression for the radius at which resonant orbits occur. Using analytic expressions for the accretion disc structure, we derive a relationship between the period excess and mass ratio with the pressure effects included. This is compared to the observed data, recently derived results for detailed integration of the disc equations and the equivalent empirically derived relations and used to predict values for the mass ratio based on measured values of the period excess for 88 systems.     

Tilted accretion discs in cataclysmic variables: tidal instabilities and superhumps

We investigate the growth of tidal instabilities in accretion discs in a binary star potential, using three-dimensional numerical simulations. As expected from analytic work, the disc is prone to an eccentric instability provided that it is large enough to extend to the 3:1 resonance. The eccentric disc leads to positive superhumps in the light curve. It has been proposed that negative superhumps might arise from a tilted disc, but we find no evidence that the companion gravitational tilt instability can grow fast enough in a fluid disc to create a measurable inclination. The origin of negative superhumps in the light curves of cataclysmic variables remains a puzzle.     

Viscous overstability and eccentricity evolution in three-dimensional gaseous discs

We investigate the growth or decay rate of the fundamental mode of even symmetry in a viscous accretion disc. This mode occurs in eccentric discs and is known to be potentially overstable. We determine the vertical structure of the disc and its modes, treating radiative energy transport in the diffusion approximation. In the limit of very long radial wavelength, an analytical criterion for viscous overstability is obtained, which involves the effective shear and bulk viscosity, the adiabatic exponent, and the opacity law of the disc. This differs from the prediction of a two-dimensional model. On shorter wavelengths (a few times the disc thickness), the criterion for overstability is more difficult to satisfy because of the different vertical structure of the mode. In a low-viscosity disc a third regime of intermediate wavelengths appears, in which the overstability is suppressed as the horizontal velocity perturbations develop significant vertical shear. We suggest that this effect determines the damping rate of eccentricity in protoplanetary discs, for which the long-wavelength analysis is inapplicable and overstability is unlikely to occur on any scale. In thinner accretion discs and in decretion discs around Be stars overstability may occur only on the longest wavelengths, leading to the preferential excitation of global eccentric modes.     

An alpha theory of time-dependent warped accretion discs

The non-linear fluid dynamics of a warped accretion disc was investigated in an earlier paper by developing a theory of fully non-linear bending waves in a thin, viscous disc. That analysis is extended here to take proper account of thermal and radiative effects by solving an energy equation that includes viscous dissipation and radiative transport. The problem is reduced to simple one-dimensional evolutionary equations for mass and angular momentum, expressed in physical units and suitable for direct application. This result constitutes a logical generalization of the alpha theory of Shakura & Sunyaev to the case of a time-dependent warped accretion disc. The local thermal–viscous stability of such a disc is also investigated.     

Kinematics of the ultracompact helium accretor AM Canum Venaticorum

We report on the results from a five-night campaign of high-speed spectroscopy of the 17-min binary AM Canum Venaticorum (AM CVn), obtained with the 4.2-m William Herschel Telescope on La Palma.
We detect a kinematic feature that appears to be entirely analogous to the 'central spike' known from the long-period, emission-line AM CVn stars GP Com, V396 Hya and SDSS J124058.03−015919.2, which has been attributed to the accreting white dwarf. Assuming that the feature indeed represents the projected velocity amplitude and phase of the accreting white dwarf, we derive a mass ratio   q = 0.18 ± 0.01  for AM CVn. This is significantly higher than the value found in previous, less direct measurements. We discuss the implications for AM CVn's evolutionary history and show that a helium star progenitor scenario is strongly favoured. We further discuss the implications for the interpretation of AM CVn's superhump behaviour, and for the detectability of its gravitational-wave signal with the Laser Interferometer Space Antenna ( LISA ).
In addition, we demonstrate a method for measuring the circularity or eccentricity of AM CVn's accretion disc, using stroboscopic Doppler tomography. We test the predictions of an eccentric, precessing disc that are based on AM CVn's observed superhump behaviour. We limit the effective eccentricity in the outermost part of the disc, where the resonances that drive the eccentricity are thought to occur, to   e = 0.04 ± 0.01  , which is smaller than previous models indicated.     

A non-linear theory of vertical resonances in accretion discs

An important and widely neglected aspect of the interaction between an accretion disc and a massive companion with a coplanar orbit is the vertical component of the tidal force. As shown by Lubow, the response of the disc to vertical forcing is resonant at certain radii, at which a localized torque is exerted, and from which a compressive wave (p mode) may be emitted. Although these vertical resonances are weaker than the corresponding Lindblad resonances, the   m =2  inner vertical resonance in a binary star is typically located within the tidal truncation radius of a circumstellar disc.
In this paper I develop a general theory of vertical resonances, allowing for non-linearity of the response, and dissipation by radiative damping and turbulent viscosity. The problem is reduced to a universal, non-linear ordinary differential equation with two real parameters. Solutions of the complex non-linear Airy equation are presented to illustrate the non-linear saturation of the resonance and the effects of dissipation. It is argued that the   m =2  inner vertical resonance is unlikely to truncate the disc in cataclysmic variable stars, but contributes to angular momentum transport and produces a potentially observable non-axisymmetric structure.     

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.     

The effect of radiation stress on gas motions in AGN accretion discs

In this paper we look at one of the effects of irradiation on a warped accretion disc in the context of active galactic nuclei (AGN). A warp will catch a substantial amount of the radiation emitted by the central object. We consider the fluid motions that may arise inside a warped disc when the surface is subject to a radiation stress, and also the net mass flows that result. We find that, to first order, we have a balance of the viscous and Coriolis-type forces. The radial radiation stress causes outward motion of the surface layer, but only the azimuthal Poynting–Robertson drag leads to an increase in the net accretion rate. We investigate the distribution of the velocity perturbations and find them to be significant in determining the local structure of the disc.
An unexpected result is that the picture changes significantly when we take into account the periodic illumination of the warped disc. A type of resonance at the local Keplerian rotation frequency causes a flow that penetrates the whole thickness of the disc; these flows are faster than the flows due to unchanging illumination. They totally dominate the induced flows in terms of sheer mass, but significant impact on disc structure still occurs only near the surface, where velocity perturbations typically go up to some kilometres per second.     

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.     

Iron line profiles and self-shadowing from relativistic thick accretion discs

We present Fe Kα line profiles from and images of relativistic discs with finite thickness around a rotating black hole using a novel code. The line is thought to be produced by iron fluorescence of a relatively cold X-ray-illuminated material in the innermost parts of the accretion disc and provides an excellent diagnostic of accretion flows in the vicinity of black holes. Previous studies have concentrated on the case of a thin, Keplerian accretion disc. This disc must become thicker and sub-Keplerian with increasing accretion rates. These can affect the line profiles and in turn can influence the estimation of the accretion disc and black hole parameters from the observed line profiles. We here embark on, for the first time, a fully relativistic computation which offers key insights into the effects of geometrical thickness and the sub-Keplerian orbital velocity on the line profiles. We include all relativistic effects such as frame-dragging, Doppler boost, time dilation, gravitational redshift and light bending. We find that the separation and the relative height between the blue and red peaks of the line profile diminish as the thickness of the disc increases. This code is also well suited to produce accretion disc images. We calculate the redshift and flux images of the accretion disc and find that the observed image of the disc strongly depends on the inclination angle. The self-shadowing effect appears remarkable for a high inclination angle, and leads to the black hole shadow being in this case, completely hidden by the disc itself.     

On an excitation mechanism for trapped inertial waves in discs around black holes

According to one model, high-frequency quasi-periodic oscillations (QPOs) can be identified with inertial waves, trapped in the inner regions of accretion discs around black holes due to relativistic effects. In order to be detected, their amplitudes need to reach large enough values via some excitation mechanism. We work out in detail a non-linear coupling mechanism suggested by Kato, in which a global warping or eccentricity of the disc has a fundamental role. These large-scale deformations combine with trapped modes to generate 'intermediate' waves of negative energy that are damped as they approach either their corotation resonance or the inner edge of the disc, resulting in amplification of the trapped waves. We determine the growth rates of the inertial modes, as well as their dependence on the spin of the black hole and the properties of the disc. Our results indicate that this coupling mechanism can provide an efficient excitation of trapped inertial waves, provided the global deformations reach the inner part of the disc with non-negligible amplitude.     

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.     

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.     

On the evolution of giant protoplanets forming in circumbinary discs

We present the results of hydrodynamic simulations of Jovian mass protoplanets that form in circumbinary discs. The simulations follow the orbital evolution of the binary plus protoplanet system acting under their mutual gravitational forces, and forces exerted by the viscous circumbinary disc. The evolution involves the clearing of the inner circumbinary disc initially, so that the binary plus protoplanet system orbits within a low density cavity. Continued interaction between disc and protoplanet causes inward migration of the planet towards the inner binary. Subsequent evolution can take three distinct paths: (i) the protoplanet enters the 4 : 1 mean motion resonance with the binary, but is gravitationally scattered through a close encounter with the secondary star; (ii) the protoplanet enters the 4 : 1 mean motion resonance, the resonance breaks, and the planet remains in a stable orbit just outside the resonance; (iii) when the binary has initial eccentricity   e _bin≥ 0.2  , the disc becomes eccentric, leading to a stalling of the planet migration, and the formation of a stable circumbinary planet.
These results have implications for a number of issues in the study of extrasolar planets. The ejection of protoplanets in close binary systems provides a source of 'free-floating planets', which have been discovered recently. The formation of a large, tidally truncated cavity may provide an observational signature of circumbinary planets during formation. The existence of protoplanets orbiting stably just outside a mean motion resonance (4 : 1) in the simulations indicate that such sites may harbour planets in binary star systems, and these could potentially be observed. Finally, the formation of stable circumbinary planets in eccentric binary systems indicates that circumbinary planets may not be uncommon.     

Discovery of drifting high-frequency quasi-periodic oscillations in global simulations of magnetic boundary layers

We report on the numerical discovery of quasi-periodic oscillations (QPOs) associated with accretion through a non-axisymmetric magnetic boundary layer in the unstable regime, when two ordered equatorial streams form and rotate synchronously at approximately the angular velocity of the inner disc. The streams hit the star's surface producing hotspots. Rotation of the spots leads to high-frequency QPOs. We performed a number of simulation runs for different magnetospheric sizes from small to tiny, and observed a definite correlation between the inner disc radius and the QPO frequency: the frequency is higher when the magnetosphere is smaller. In the stable regime, a small magnetosphere forms and accretion through the usual funnel streams is observed, and the frequency of the star is expected to dominate the light curve. We performed exploratory investigations of the case in which the magnetosphere becomes negligibly small and the disc interacts with the star through an equatorial belt. We also performed investigation of somewhat larger magnetospheres where one or two ordered tongues may dominate over other chaotic tongues. In application to millisecond pulsars, we obtain QPO frequencies in the range of 350–990 Hz for one spot. The frequency associated with rotation of one spot may dominate if spots are not identical or antipodal. If the spots are similar and antipodal, then the frequencies are twice as high. We show that variation of the accretion rate leads to drift of the QPO peak.     

Implications of mean field accretion disc theory for vorticity and magnetic field growth

In addition to the scalar Shakura–Sunyaev α _ss turbulent viscosity transport term used in simple analytic accretion disc modelling, a pseudo-scalar transport term also arises. The essence of this term can be captured even in simple models for which vertical averaging is interpreted as integration over a half-thickness and each hemisphere is separately studied. The additional term highlights a complementarity between mean field magnetic dynamo theory and accretion disc theory treated as a mean field theory. Such pseudo-scalar terms have been studied, and can lead to large-scale magnetic field and vorticity growth. Here it is shown that vorticity can grow even in the simplest azimuthal and half-height integrated disc model, for which mean quantities depend only on radius. The simplest vorticity growth solutions seem to have scales and vortex survival times consistent with those required for facilitating planet formation. In addition, it is shown that, when the magnetic back-reaction is included to lowest order, the pseudo-scalar driving the magnetic field growth and that driving the vorticity growth will behave differently with respect to shearing and non-shearing flows: the former pseudo‐scalar can more easily reverse sign in the two cases.     

Comparisons and connections between mean field dynamo theory and accretion disc theory

The origin of large scale magnetic fields in astrophysical rotators, and the conversion of gravitational energy into radiation near stars and compact objects via accretion have been subjects of active research for a half century. Magnetohydrodynamic turbulence makes both problems highly nonlinear, so both subjects have benefitted from numerical simulations.However, understanding the key principles and practical modeling of observations warrants testable semi‐analytic mean field theories that distill the essential physics. Mean field dynamo (MFD) theory and alpha‐viscosity accretion disc theory exemplify this pursuit. That the latter is a mean field theory is not always made explicit but the combination of turbulence and global symmetry imply such. The more commonly explicit presentation of assumptions in 20th century textbook MFDT has exposed it to arguably more widespread criticism than incurred by 20th century alpha‐accretion theory despite complementary weaknesses. In the 21st century however, MFDT has experienced a breakthrough with a dynamical saturation theory that consistently agrees with simulations. Such has not yet occurred in accretion disc theory, though progress is emerging. Ironically however, for accretion engines, MFDT and accretion theory are presently two artificially uncoupled pieces of what should be a single coupled theory. Large scale fields and accretion flows are dynamically intertwined because large scale fields likely play a key role in angular momentum transport. I discuss and synthesize aspects of recent progress in MFDT and accretion disc theory to suggest why the two likely conspire in a unified theory (© 2010 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

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.