Accretion disc dynamics in extreme mass ratio compact binaries

An analysis is presented of a numerical investigation of the dynamics and geometry of accretion discs in binary systems with mass ratios   q = M ₂/ M ₁ < 0.1  , applicable to ultracompact X-ray binaries, AM CVn stars and very short period cataclysmic variables. The steady-state geometry of the disc in the binary reference frame is found to be quite different from that expected at higher mass ratios. For   q ∼ 0.1  , the disc takes on the usual elliptical shape, with the major axis aligned perpendicular to the line of centres of the two stars. However, at smaller mass ratios the elliptical gaseous orbits in the outer regions of the disc are rotated in the binary plane. The angle of rotation increases with gas temperature, but is found to vary inversely with q . At   q = 0.01  , the major axis of these orbits is aligned almost parallel to the line of centres of the two stars. These effects may be responsible for the similar disc structure inferred from Doppler tomography of the AM CVn star GP Com, which has   q = 0.02  . The steady-state geometry at low mass ratios is not predicted by an inviscid, restricted three-body model of gaseous orbits; it is related to the effects of tidal-viscous truncation of the disc near the Roche lobe boundary. Since the disc geometry can be inferred observationally for some systems, it is proposed that this may offer a useful diagnostic for the determination of mass ratios in ultracompact binaries.     

Atlas of tilted accretion discs and source to negative superhumps

Using smoothed particle hydrodynamics, we numerically simulate steady-state accretion discs for cataclysmic variable dwarf novae systems that have a secondary-to-primary mass ratio  0.35 ≤ q ≤ 0.55  . After these accretion discs have come to quasi-equilibrium, we rotate each disc out of the orbital plane by  δ= (1, 2, 3, 4, 5 or 20)°  to induce negative superhumps. For accretion discs tilted  5°  , we generate light curves and associated Fourier transforms for an atlas on negative superhumps and retrograde precession. Our simulation results suggest that accretion discs need to be tilted more than 3° for negative superhumps to be statistically significant. We also show that if the disc is tilted enough such that the gas stream strikes a disc face, then a dense cooling ring is generated near the radius of impact.
In addition to the atlas, we study these artificially tilted accretion discs to find the source to negative superhumps. Our results suggest that the source is additional light from innermost disc annuli, and this additional light waxes and wanes with the amount of gas stream overflow received as the secondary orbits. The nodes, where the gas stream transitions from flowing over to under the disc rim (and vice versa), precess in the retrograde direction.     

On the precession of accretion discs in X-ray binaries

In this Letter, recent results on the nodal precession of accretion discs in close binaries are applied to the discs in some X-ray binary systems. The ratio between the tidally forced precession period and the binary orbital period is given, as well as the condition required for the rigid precession of gaseous Keplerian discs. Hence the minimum precessional period that may be supported by a fluid Keplerian disc is determined. It is concluded that near-rigid body precession of tilted accretion discs can occur and generally reproduce observationally inferred precession periods, for reasonable system parameters. In particular, long periods in SS 433, Her X-1, LMC X-4 and SMC X-1 can be fitted by the tidal model. It is also found that the precession period that has been tentatively put forward for Cyg X-2 cannot be accommodated by a tidally precessing disc model for any realistic choice of system parameters.     

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.     

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.     

Disc precession in the M31 dipping X-ray binary Bo 158?

We present results from three XMM–Newton observations of the M31 low mass X-ray binary (LMXB) XMMU J004314.4+410726.3 (Bo 158), spaced over 3 d in 2004 July. Bo 158 was the first dipping LMXB to be discovered in M31. Periodic intensity dips were previously seen to occur on a 2.78-h period, due to absorption in material that is raised out of the plane of the accretion disc. The report of these observations stated that the dip depth was anticorrelated with source intensity. In light of the 2004 XMM–Newton observations of Bo 158, we suggest that the dip variation is due to precession of the accretion disc. This is to be expected in LMXBs with a mass ratio ≲0.3 (period ≲4 h), as the disc reaches the 3:1 resonance with the binary companion, causing elongation and precession of the disc. A smoothed particle hydrodynamics simulation of the disc in this system shows retrograde rotation of a disc warp on a period of  ∼11 P _orb  , and prograde disc precession on a period of  29 ± 1 P _orb  . This is consistent with the observed variation in the depth of the dips. We find that the dipping behaviour is most likely to be modified by the disc precession, hence we predict that the dipping behaviour repeats on an  81 ± 3 h  cycle.     

Discovery of long-term superorbital periodicities in the pseudo-transient LMXBs: IGR J17098−3628 and EXO 0748−676

Long-term monitoring of the recently discovered X-ray transient, IGR J17098−3628, by the All-Sky Monitor on-board the Rossi X-Ray Timing Explorer , has shown that it displays a long-term (≈163 d) quasi-periodic modulation in the data spanning its 'active' state (i.e. approximately MJD 53450–54200). Furthermore, this light curve is not typical of 'classical' soft X-ray transients, in that J17098−3628 has remained active since its initial discovery, and may be more akin to the pseudo-transient EXO 0748−676, which is now classified as a persistent low-mass X-ray binary (LMXB). However, EXO 0748−676 recently entered a more active phase (since approximately MJD 53050), and since then we find that it too displays a quasi-periodic modulation (≈181 d) in its light curve. This must be a 'superorbital' modulation, as the orbital period of EXO 0748−676 is well established (3.8 h), and hence we interpret both objects' long periods as representing some intrinsic properties of the accretion disc (such as coupled precessional and warping effects). By analogy, we therefore suggest that IGR J17098−3628 is another member of this class of pseudo-transient LMXBs and is likely to have a <1 d orbital period.     

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.     

SPH simulations of negative (nodal) superhumps: a parametric study

Negative superhumps in cataclysmic variable systems result when the accretion disc is tilted with respect to the orbital plane. The line of nodes of the tilted disc precesses slowly in the retrograde direction, resulting in a photometric signal with a period slightly less than the orbital period. We use the method of smoothed particle hydrodynamics to simulate a series of models of differing mass ratio and effective viscosity to determine the retrograde precession period and superhump period deficit  ɛ₋  as a function of system mass ratio q . We tabulate our results and present fits to both  ɛ₋  and  ɛ₊  versus q , as well as compare the numerical results with those compiled from the literature of negative superhump observations. One surprising result is that while we find negative superhumps most clearly in simulations with an accretion stream present, we also find evidence for negative superhumps in simulations in which we shut off the mass transfer stream completely, indicating that the origin of the photometric signal is more complicated than previously believed.     

The standstill luminosity in Z Cam systems

We consider accretion discs in close binary systems. We show that the heating of a disc at the impact point of the accretion stream contributes significantly to the local energy budget at its outer edge. As a result, the thermal balance relation between local accretion rate and surface density (the 'S-curve') changes; the critical mass transfer rate above which no dwarf nova outbursts occur can be up to 40 per cent smaller than without impact heating. Standstills in Z Cam systems thus occur at smaller mass transfer rates than otherwise expected, and are fainter than the peak luminosity during the dwarf nova phase as a result.     

Substantial stream–disc overflow found in three-dimensional SPH simulations of cataclysmic variables

We study numerically the interaction of the infalling gas stream and the rim of the accretion disc in cataclysmic variables. The simulations were performed with a smoothed particle hydrodynamics scheme with high spatial resolution. Parameters of the systems AM CVn, OY Car, DQ Her, U Gem and IP Peg were used for the simulations. The simulations cover a wide range of orbital periods, mass ratios and mass transfer rates, as well as different thermal states of the accretion disc. The main result of this study is that the accretion stream is not stopped at the impact region (the bright spot at the outer rim of the disc). In fact, after undergoing the shock interaction, most of the matter is deflected vertically and flows in a more or less diffuse stream to inner parts of the disc, hitting the disc surface close to the circularization radius at orbital phase 0.5. This is a common feature in all systems for all simulated parameters. This stream overflow can cause the X-ray absorption dips observed in cataclysmic variables (CVs) and low-mass X-ray binaries (LMXBs) around orbital phase 0.7, if the inclination is at least 65°. Under certain circumstances, namely a sudden increase of the mass transfer rate from the secondary or a rather small disc, parts of the overflowing stream bounce off the disc surface after hitting it at orbital phase ≈0.5. Another absorption region can be expected around orbital phase 0.2.
In our simulations most of the infalling matter reaches the inner disc very quickly. This must alter the evolution of the quiescent disc and the outburst behaviour considerably compared with purely viscous transport of the material through the disc from the outer rim, and therefore should be taken into account in dwarf nova outburst cycle calculations. To our knowledge, the consequences of such a massive stream overflow for the dwarf nova outburst cycle have not been considered yet.     

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.     

A comparative analysis of standard accretion discs spectra: an application to ultraluminous X-ray sources

We compare standard models of accretion discs around black holes (BHs) that include the appropriate zero-torque inner boundary condition and relativistic effects on the emission and propagation of radiation. The comparison is performed adopting the multicolour disc blackbody model (MCD) as reference and looking for the parameter space in which it is in statistical agreement with 'more physical' accretion disc models. We find simple 'recipes' that can be used for adjusting the estimates of the physical inner radius of the disc, the BH mass and the accretion rate inferred using the parameters of the MCD fits. We applied these results to four ultraluminous X-ray sources for which MCD spectral fits of their X-ray soft spectral components have been published and find that, in three cases (NGC 1313 X-1, X-2 and M 81 X-9), the BH masses inferred for a standard disc around a Schwarzschild BH are in the interval  ∼100–200 M_⊙  . Only if the BH is maximally rotating are the masses comparable to the much larger values previously derived in the literature.     

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.     

Optical photometry and spectroscopy of the accretion-powered millisecond pulsar HETE J1900.1 − 2455

We present phase resolved optical photometry and spectroscopy of the accreting millisecond pulsar HETE J1900.1−2455. Our R -band light curves exhibit a sinusoidal modulation, at close to the orbital period, which we initially attributed to X-ray heating of the irradiated face of the secondary star. However, further analysis reveals that the source of the modulation is more likely due to superhumps caused by a precessing accretion disc. Doppler tomography of a broad Hα emission line reveals an emission ring, consistent with that expected from an accretion disc. Using the velocity of the emission ring as an estimate for the projected outer disc velocity, we constrain the maximum projected velocity of the secondary to be 200 km s⁻¹, placing a lower limit of  0.05 M_⊙  on the secondary mass. For a  1.4 M_⊙  primary, this implies that the orbital inclination is low, ≲20°. Utilizing the observed relationship between the secondary mass and the orbital period in short-period cataclysmic variables, we estimate the secondary mass to be ∼0.085  M_⊙  , which implies an upper limit of  ∼2.4 M_⊙  for the primary mass.     

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.     

XMM–Newton observations of UW CrB: detection of X-ray bursts and evidence for accretion disc evolution

UW CrB (MS 1603+2600) is a peculiar short-period X-ray binary that exhibits extraordinary optical behaviour. The shape of the optical light curve of the system changes drastically from night to night, without any changes in overall brightness. Here we report X-ray observations of UW CrB obtained with XMM–Newton . We find evidence for several X-ray bursts, confirming a neutron star primary. This considerably strengthens the case that UW CrB is an accretion disc corona system located at a distance of at least 5–7 kpc (3–5 kpc above the Galactic plane). The X-ray and Optical Monitor (ultraviolet–optical) light curves show remarkable shape variation from one observing run to another, which we suggest are due to large-scale variations in the accretion disc shape resulting from a warp that periodically obscures the optical and soft X-ray emission. This is also supported by the changes in phase-resolved X-ray spectra.     

The evolution of misaligned accretion discs and spinning black holes

In this paper, we consider the process of alignment of a spinning black hole and a surrounding misaligned accretion disc. We use a simplified set of equations, that describe the evolution of the system in the case where the propagation of warping disturbances in the accretion disc occurs diffusively, a situation likely to be common in the thin discs in active galactic nuclei (AGN). We also allow the direction of the hole spin to move under the action of the disc torques. In such a way, the evolution of the hole–disc system is computed self-consistently. We consider a number of different situations and we explore the relevant parameter range, by varying the location of the warp radius R _w and the propagation speed of the warp. We find that the dissipation associated with the twisting of the disc results in a large increase in the accretion rate through the disc, so that AGN accreting from a misaligned disc are likely to be significantly more luminous than those accreting from a flat disc. We compute explicitly the time-scales for the warping of the disc and for the alignment process and compare our results with earlier estimates based on simplified steady-state solutions. We also confirm earlier predictions that, under appropriate circumstances, accretion can proceed in a counter-aligned fashion, so that the accreted material will spin-down the hole, rather than spinning it up. Our results have implication in a number of different observational features of AGN such as the orientation and shape of jets, the shape of X-ray iron lines and the possibility of obscuration and absorption of X-ray by the outer disc as well as the general issue of the spin history of black holes during their growth.     

Black hole mergers: can gas discs solve the 'final parsec' problem?

We compute the effect of an orbiting gas disc in promoting the coalescence of a central supermassive black hole binary. Unlike earlier studies, we consider a finite mass of gas with explicit time dependence: we do not assume that the gas necessarily adopts a steady state or a spatially constant accretion rate, i.e. that the merging black hole was somehow inserted into a pre-existing accretion disc. We consider the tidal torque of the binary on the disc, and the binary's gravitational radiation. We study the effects of star formation in the gas disc in a simple energy feedback framework.
The disc spectrum differs in detail from that found before. In particular, tidal torques from the secondary black hole heat the edges of the gap, creating bright rims around the secondary. These rims do not in practice have uniform brightness either in azimuth or time, but can on average account for as much as 50 per cent of the integrated light from the disc. This may lead to detectable high-photon-energy variability on the relatively long orbital time-scale of the secondary black hole, and thus offer a prospective signature of a coalescing black hole binary.
We also find that the disc can drive the binary to merger on a reasonable time-scale only if its mass is at least comparable with that of the secondary black hole, and if the initial binary separation is relatively small, i.e.   a ₀≲ 0.05  pc. Star formation complicates the merger further by removing mass from the disc. In the feedback model we consider, this sets an effective limit to the disc mass. As a result, binary merging is unlikely unless the black hole mass ratio is ≲0.001. Gas discs thus appear not to be an effective solution to the 'last parsec' problem for a significant class of mergers.     

Superhumps in low-mass X-ray binaries

We propose a mechanism for the superhump modulations observed in optical photometry of at least two black-hole X-ray transients (SXTs). As in extreme mass-ratio cataclysmic variables (CVs), superhumps are assumed to result from the presence of the 3:1 orbital resonance in the accretion disc. This causes the disc to become non-axisymmetric and precess. However, the mechanism for superhump luminosity variations in low-mass X-ray binaries (LMXBs) must differ from that in CVs, where it is attributed to a tidally-driven modulation of the disc's viscous dissipation, varying on the beat between the orbital and disc precession period. By contrast in LMXBs, tidal dissipation in the outer accretion disc is negligible: the optical emission is overwhelmingly dominated by reprocessing of intercepted central X-rays. Thus a different origin for the superhump modulation is required. Recent observations and numerical simulations indicate that in an extreme mass-ratio system the disc area changes on the superhump period. We deduce that the superhumps observed in SXTs arise from a modulation of the reprocessed flux by the changing area. Therefore, unlike the situation in CVs, where the superhump amplitude is inclination-independent, superhumps should be best seen in low-inclination LMXBs, whereas an orbital modulation from the heated face of the secondary star should be more prominent at high inclinations. Modulation at the disc precession period (10 s of days) may indicate disc asymmetries such as warping. We comment on the orbital period determinations of LMXBs, and the possibility and significance of possible permanent superhump LMXBs.