Precessing warped accretion discs in X-ray binaries

We study the radiation-driven warping of accretion discs in the context of X-ray binaries. The latest evolutionary equations are adopted, which extend the classical alpha theory to time-dependent thin discs with non-linear warps. We also develop accurate, analytical expressions for the tidal torque and the radiation torque, including self-shadowing.
We investigate the possible non-linear dynamics of the system within the framework of bifurcation theory. First, we re-examine the stability of an initially flat disc to the Pringle instability. Then we compute directly the branches of non-linear solutions representing steadily precessing discs. Finally, we determine the stability of the non-linear solutions. Each problem involves only ordinary differential equations, allowing a rapid, accurate and well-resolved solution.
We find that radiation-driven warping is probably not a common occurrence in low-mass X-ray binaries. We also find that stable, steadily precessing discs exist for a narrow range of parameters close to the stability limit. This could explain why so few systems show clear, repeatable 'superorbital' variations. The best examples of such systems, Her X-1, SS 433 and LMC X-4, all lie close to the stability limit for a reasonable choice of parameters. Systems far from the stability limit, including Cyg X-2, Cen X-3 and SMC X-1, probably experience quasi-periodic or chaotic variability as first noticed recently by Wijers and Pringle. We show that radiation-driven warping provides a coherent and persuasive framework but that it does not provide a generic explanation for the long-term variabilities in all X-ray binaries.     

Warped accretion discs and the long periods in X-ray binaries

Precessing accretion discs have long been suggested as explanations for the long periods observed in a variety of X-ray binaries, most notably Her X-1/HZ Her. We show that an instability of the response of the disc to the radiation reaction force from the illumination by the central source can cause the disc to tilt out of the orbital plane and precess in something like the required manner. The rate of precession and disc tilt obtained for realistic values of system parameters compare favourably with the known body of data on X-ray binaries with long periods. We explore other possible types of behaviour than steadily precessing discs that might be observable in systems with somewhat different parameters. At high luminosities, the inner disc tilts through more than 90°, i.e., it rotates counter to the usual direction, which may explain the torque reversals in systems such as 4U 1626−67.     

Doppler tomography of accretion in binaries

Since its conception, Doppler tomography has matured into a versatile and widely used tool. It exploits the information contained in the highly‐structured spectral line‐profiles typically observed in mass‐transferring binaries. Using inversion techniques akin to medical imaging, it permits the reconstruction of Doppler maps that image the accretion flow on micro‐arcsecond scales. I summarise the basic concepts behind the technique and highlight two recent results; the use of donor star emission as a means to system parameter determination, and the real‐time movies of the evolving accretion flow in the cataclysmic variable WZ Sge during its 2001 outburst. I conclude with future opportunities in Doppler tomography by exploiting the combination of superior data sets, second generation reconstruction codes and simulated theoretical tomograms to delve deeper into the physics of accretion flows. (© 2004 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Impact of reverberation in flared accretion discs on temporal characteristics of X-ray binaries

Observations suggest that accretion discs in many X-ray binaries are likely flared. An outer edge of the disc intercepts radiation from the central X-ray source. Part of that radiation is absorbed and re-emitted in the optical/UV spectral ranges. However, a large fraction of that radiation is reflected and appears in the broad-band X-ray spectrum as a Compton reflection bump. This radiation is delayed and variability is somewhat smeared compared with the intrinsic X-ray radiation. We compute response functions for flat and flared accretion discs and for isotropic and anisotropic X-ray sources. A simple approximation for the response function which is valid in the broad range of the disc shapes and inclinations, inner and outer radii, and the plasma bulk velocity is proposed. We also study the impact of the X-ray reprocessing on temporal characteristics of X-ray binaries such as the power spectral density, auto- and cross-correlation functions, and time/phase lags. We propose a reprocessing model which explains the secondary peaks in the phase lag Fourier spectra observed in Cyg X-1 and other Galactic black hole sources. The position of the peaks could be used to determine the size of the accretion disc.     

Three-dimensional smoothed particle hydrodynamics simulations of radiation-driven warped accretion discs

We present three-dimensional smoothed particle hydrodynamics calculations of warped accretion discs in X-ray binary systems. Geometrically thin, optically thick accretion discs are illuminated by a central radiation source. This illumination exerts a non-axisymmetric radiation pressure on the surface of the disc, resulting in a torque that acts on the disc to induce a twist or warp. Initially planar discs are unstable to warping driven by the radiation torque and, in general, the warps also precess in a retrograde direction relative to the orbital flow. We simulate a number of X-ray binary systems which have different mass ratios, using a number of different luminosities for each. Radiation-driven warping occurs for all systems simulated. For mass ratios   q ∼ 0.1  a moderate warp occurs in the inner disc while the outer disc remains in the orbital plane (cf. X 1916−053). For less extreme mass ratios, the entire disc tilts out of the orbital plane (cf. Her X–1). For discs that are tilted out of the orbital plane in which the outer edge material of the disc is precessing in a prograde direction, we obtain both positive and negative superhumps simultaneously in the dissipation light curve (cf. V603 Aql).     

On the observability of spiral structures in cataclysmic variable accretion discs

We use the grid of hydrodynamic accretion disc calculations of Stehle to construct orbital phase‐dependent emission‐line profiles of thin discs carrying spiral density waves. The observational signatures of spiral waves are explored to establish the feasibility of detecting spiral waves in cataclysmic variable discs using prominent emission lines in the visible range of the spectrum. For high Mach number accretion discs ( M v _φ c _s≃ 15 – 30), we find that the spiral shock arms are so tightly wound that they leave few obvious fingerprints in the emission lines. Only a minor variation of the double peak separation in the line profile at a level of ∼8 per cent is produced. For accretion discs in outburst ( M ≃ 5 – 20) however, the lines are dominated by the emission from an m =2 spiral pattern in the disc. We show that reliable Doppler tomograms of spiral shock patterns can be reconstructed provided that a signal‐to‐noise ratio of at least 15, a wavelength resolution of ∼80 km s⁻¹ and a time resolution of ∼50 spectra per binary orbit are achieved. We confirm that the observed spiral pattern in the disc of IP Pegasi can be reproduced by tidal density waves in the accretion disc and demands the presence of a large, hot disc, at least in the early outburst stages.     

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.     

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.     

Instability of spherical accretion — I. Shock-free Bondi accretion

We examine the spatial stability of spherical adiabatic Bondi accretion on to a point gravitating mass against external perturbations. Both transonic critical and subsonic subcritical accretion are shown to be stable against purely radial acoustic, vortex or entropy perturbations. In the case of non-radial perturbations the amplitude of the perturbations grows without limit with smaller radii. Instability manifests itself only if the size of the accreting body is much less than the Bondi radius so that the inflow is highly supersonic or highly subsonic at the surface of the accretor in the case of critical or subcritical accretion respectively. These asymptotics hold and consequently the instability may develop for adiabatic index of accreting gas γ < 5/3. We suggest that this instability may lead to an essential thermalization of accreting flow thus, particularly, solving the problem of otherwise inefficient energy release in spherical accretion on to a black hole.     

The effective temperature distribution of steady-state, mass-losing accretion discs

Mass loss appears to be a common phenomenon among astrophysical accretion disc systems. An outflow emanating from an accretion disc can act as a sink for mass, angular momentum and energy, and can therefore alter the dissipation rates and effective temperatures across the disc. Here, the radial distributions of dissipation rate and effective temperature across a Keplerian, steady-state, mass-losing accretion disc are derived, using a simple, parametric approach that is sufficiently general to be applicable to many types of dynamical disc–wind models.
Effective temperature distributions for mass-losing accretion discs in cataclysmic variables are shown explicitly, with parameters chosen to describe both radiation-driven and centrifugally driven outflows. For realistic wind mass-loss rates of a few per cent, only centrifugally driven outflows – particularly those in which mass loss is concentrated in the inner disc – are likely to alter the effective temperature distribution of the disc significantly. Accretion discs that drive such outflows could produce spectra and eclipse light curves that are noticeably different from those produced by standard, conservative discs.     

The correlations between the twin kHz quasi-periodic oscillation frequencies of low-mass X-ray binaries

We analyzed the recently published kHz quasi-period oscillaiton (QPO) data in the neutron star low-mass X-ray binaries (LMXBs), in order to investigate the different correlations of the twin-peak kHz QPOs in bright Z sources and in the less luminous Atoll sources. We find a power-law relation  ν₁∼ν ^b ₂  between the upper and the lower kHz QPOs with different indices: b ≃ 1.5 for the Atoll source 4U 1728-34 and b ≃ 1.9 for the Z source Sco X-1. The implications of our results for the theoretical models for kHz QPOs are discussed.