The light curves of soft X-ray transients

We show that the light curves of soft X-ray transients (SXTs) follow naturally from the disc instability picture, adapted to take account of irradiation by the central X-ray source during the outburst. Irradiation prevents the disc from returning to the cool state until central accretion is greatly reduced. This happens only after most of the disc mass has been accreted by the central object, on a viscous time-scale, accounting naturally for the exponential decay of the outburst on a far longer time-scale (τ20–40 d) than seen in dwarf novae, without any need to manipulate the viscosity parameter α. The accretion of most of the disc mass in outburst explains the much longer recurrence time of SXTs compared with dwarf novae. This picture also suggests an explanation of the secondary maximum seen in SXT light curves about 50–75 d after the start of each outburst, since central irradiation triggers the thermal instability of the outer disc, adding to the central accretion rate one viscous time later. The X-ray outburst decay constant τ should on average increase with orbital period, but saturate at a roughly constant value ∼40 d for orbital periods longer than about a day. The bolometric light curve should show a linear rather than an exponential decay at late times (a few times τ). Outbursts of long-period systems should be entirely in the linear decay regime, as is observed in GRO J1744−28. UV and optical light curves should resemble the X-rays but have decay time-scales up to 2–4 times longer.     

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.     

Magnetically moderated outbursts of WZ Sagittae

We argue that the quiescent value of the viscosity parameter of the accretion disc in WZ Sge may be  α_cold∼ 0.01  , in agreement with estimates of α_cold for other dwarf novae. Assuming the white dwarf in WZ Sge to be magnetic, we show that, in quiescence, material close to the white dwarf can be propelled to larger radii, depleting the inner accretion disc. The propeller therefore has the effect of stabilizing the inner disc and allowing the outer disc to accumulate mass. The outbursts of WZ Sge are then regulated by the (magnetically determined) evolution of the surface density of the outer disc at a radius close to the tidal limit. Numerical models confirm that the recurrence time can be significantly extended in this way. The outbursts are expected to be superoutbursts since the outer disc radius is forced to exceed the tidal (3:1 resonance) radius. The large, quiescent disc is expected to be massive, and to be able to supply the observed mass accretion rate during outburst. We predict that the long-term spin evolution of the white dwarf spin will involve a long cycle of spin-up and spin-down phases.     

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.     

Investigating the structure of the accretion disc in WZ Sge from multiwaveband time-resolved spectroscopic observations – I

We present the first of two papers describing an in-depth study of multiwaveband phase-resolved spectroscopy of the unusual dwarf nova WZ Sge. In this paper we present an extensive set of Doppler maps of WZ Sge covering optical and infrared emission lines, and describe a new technique for studying the accretion discs of cataclysmic variables using ratioed Doppler maps. Applying the ratioed Doppler map technique to our WZ Sge data shows that the radial temperature profile of the disc is unlike that predicted for a steady state α disc. Time-averaged spectra of the accretion disc line flux (with the bright spot contribution removed) show evidence in the shapes of the line profiles for the presence of shear broadening in a quiescent non-turbulent accretion disc. From the positions of the bright spots in the Doppler maps of different lines, we conclude that the bright spot region is elongated along the ballistic stream, and that the density of the outer disc is low. The velocity of the outer edge of the accretion disc measured from the H α line is found to be 723±23 km s⁻¹. Assuming that the accretion disc reaches to the 3:1 tidal resonance radius, we derive a value for the primary star mass of 0.82 M_⊙. We discuss the implications of our results on the present theories of WZ Sge type dwarf nova outbursts.     

Infrared spectroscopy of cataclysmic variables – III. Dwarf novae below the period gap and nova-like variables

We present K -band spectra of the short-period dwarf novae YZ Cnc, LY Hya, BK Lyn, T Leo, SW UMa and WZ Sge, the nova-like variables DW UMa, V1315 Aql, RW Tri, VY Scl, UU Aqr and GP Com, and a series of field dwarf stars with spectral types ranging from K2 to M6.
The spectra of the dwarf novae are dominated by emission lines of H  i and He  i . The large velocity and equivalent widths of these lines, in conjunction with the fact that the lines are double-peaked in the highest inclination systems, indicate an accretion disc origin. In the case of YZ Cnc and T Leo, for which we obtained time-resolved data covering a complete orbital cycle, the emission lines show modulations in their equivalent widths that are most probably associated with the bright spot (the region where the gas stream collides with the accretion disc). There are no clear detections of the secondary star in any of the dwarf novae below the period gap, yielding upper limits of 10–30 per cent for the contribution of the secondary star to the observed K -band flux. In conjunction with the K -band magnitudes of the dwarf novae, we use the derived secondary star contributions to calculate lower limits to the distances to these systems.
The spectra of the nova-like variables are dominated by broad, single-peaked emission lines of H  i and He  i – even the eclipsing systems we observed do not show the double-peaked profiles predicted by standard accretion disc theory. With the exception of RW Tri, which exhibits Na  i , Ca  i and ¹²CO absorption features consistent with a M0V secondary contributing 65 per cent of the observed K -band flux, we find no evidence for the secondary star in any of the nova-like variables. The implications of this result are discussed.     

The X-ray and extreme-ultraviolet flux evolution of SS Cygni throughout outburst

We present the most complete multiwavelength coverage of any dwarf nova outburst: simultaneous optical, Extreme Ultraviolet Explorer and Rossi X-ray Timing Explorer observations of SS Cygni throughout a narrow asymmetric outburst. Our data show that the high-energy outburst begins in the X-ray waveband 0.9–1.4 d after the beginning of the optical rise and 0.6 d before the extreme-ultraviolet rise. The X-ray flux drops suddenly, immediately before the extreme-ultraviolet flux rise, supporting the view that both components arise in the boundary layer between the accretion disc and white dwarf surface. The early rise of the X-ray flux shows that the propagation time of the outburst heating wave may have been previously overestimated.
The transitions between X-ray and extreme-ultraviolet dominated emission are accompanied by intense variability in the X-ray flux, with time-scales of minutes. As detailed by Mauche & Robinson, dwarf nova oscillations are detected throughout the extreme-ultraviolet outburst, but we find they are absent from the X-ray light curve.
X-ray and extreme-ultraviolet luminosities imply accretion rates of  3 × 10¹⁵ g s⁻¹  in quiescence,  1 × 10¹⁶ g s⁻¹  when the boundary layer becomes optically thick, and  ∼10¹⁸ g s⁻¹  at the peak of the outburst. The quiescent accretion rate is two and a half orders of magnitude higher than predicted by the standard disc instability model, and we suggest this may be because the inner accretion disc in SS Cyg is in a permanent outburst state.     

A transient I-band excess in the optical spectrum of the accreting millisecond pulsar SAX J1808.4–3658

The optical counterpart of the transient, millisecond X-ray pulsar SAX J1808.4–3658 was observed in four colours ( BVRI ) for five weeks during the 2005 June–July outburst. The optical fluxes declined by ∼2 mag during the first 16d and then commenced quasi-periodic secondary outbursts, with time-scales of several days, similar to those seen in 2000 and 2002. The broad-band spectra derived from these measurements were generally consistent with emission from an X-ray heated accretion disc. During the first 16d decline in intensity the spectrum became redder. We suggest that the primary outburst was initiated by a viscosity change driven instability in the inner disc and note the contrast with another accreting millisecond pulsar, XTE J0929−314, for which the spectrum becomes bluer during the decline. On the night of 2005 June 5 (HJD 245 3527) the I -band flux was ∼0.45-mag brighter than on the preceding or following nights whereas the BV and R bands showed no obvious enhancement. A type I X-ray burst was detected by the Rossi X-ray Timing Explorer spacecraft during this I -band integration. It seems unlikely that reprocessed radiation from the burst was sufficient to explain the observed increase. We suggest that a major part of the I -band excess was due to synchrotron emission triggered by the X-ray burst. Several other significant short duration changes in V − I were detected. One occurred at about HJD 245 3546 in the early phase of the first secondary outburst and may be due to mass-transfer instability or to another synchrotron emission event.     

Progress on a model for β Lyrae

It is generally accepted that β Lyrae is a massive mass-transfer system in which the more massive gainer is embedded in a thick accretion disc. No detailed, quantitative accretion disc model of the system, which fits both the available photometric and spectroscopic data, has been published.
This paper considers a thick accretion disc model. An empirical two-temperature model for the accretion disc rim produces an accurate fit to the UBV and OAO2 photometry. The two-temperature rim model leads to synthetic spectra with a much larger Balmer jump than is present in spectrophotometric scans. The same model fails to produce IR light curves with the 0.5^p eclipse deeper than the 0.0^p eclipse, as shown in the observations of Jameson & Longmore.
The paper next reports initial results for a model in which an extra source of continuum radiation arises from Thomson scattering of radiation from the gainer, which is hidden from the observer by the thick accretion disc.
The adopted rate of mass transfer cannot supply the energy radiated by the accretion disc rim, nor can any credible larger rate. The most likely energy source is the mass gainer.     

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.     

Investigating the structure of the accretion disc in WZ Sge from multiwaveband time-resolved spectroscopic observations – II

We present our second paper describing multiwaveband time-resolved spectroscopy of WZ Sge. We analyse the evolution of both optical and IR emission lines throughout the orbital period and find evidence, in the Balmer lines, for an optically thin accretion disc and an optically thick hotspot. Optical and IR emission lines are used to compute radial velocity curves. Fits to our radial velocity measurements give an internally inconsistent set of values for K ₁, γ and the phase of red-to-blue crossing. We present a probable explanation for these discrepancies, and provide evidence for similar behaviour in other short orbital period dwarf novae. Selected optical and IR spectra are measured to determine the accretion disc radii. Values for the disc radii are found to be strongly dependent on the assumed WD mass and binary orbital inclination. However, the separation of the peaks in the optical emission line (i.e., an indication of the outer disc radius) has been found to be constant during all phases of the supercycle period over the last 40 years.     

Stream overflow and dwarf nova eruptions

We consider the effects of accretion stream overflow on the viscous dynamics of accretion discs in dwarf novae. If the stream from the secondary star is geometrically thick enough, some fraction of its material can flow over and under the disc. The mass and specific angular momentum of the stream are then deposited not only at the point of collision with the outer disc, but also at those radii in the inner disc with geometric heights that are large enough to intercept the residual stream, or near the radius where the disc has the same specific angular momentum as the stream. The overflowing stream can alter the behaviour of heating fronts and cooling fronts in the disc. If the mass fraction of the overflowing stream is of order tens of per cent, the deposition of mass in the inner parts of the disc is sufficient to change the character of the eruption light curves significantly.     

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)     

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.     

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.     

Physics of accretion in the millisecond pulsar XTE J1751 − 305

We have undertaken an extensive study of X-ray data from the accreting millisecond pulsar XTE J1751 − 305 observed by RXTE and XMM–Newton during its 2002 outburst. In all aspects this source is similar to the prototypical millisecond pulsar SAX J1808.4 − 3658, except for the higher peak luminosity of 13 per cent of Eddington, and the optical depth of the hard X-ray source, which is larger by a factor ∼2. Its broad-band X-ray spectrum can be modelled by three components. We interpret the two soft components as thermal emission from a colder  ( kT ∼ 0.6 keV)  accretion disc and a hotter (∼1 keV) spot on the neutron star surface. We interpret the hard component as thermal Comptonization in plasma of temperature ∼40 keV and optical depth ∼1.5 in a slab geometry. The plasma is heated by the accretion shock as the material collimated by the magnetic field impacts on to the surface. The seed photons for Comptonization are provided by the hotspot, not by the disc. The Compton reflection is weak and the disc is probably truncated into an optically thin flow above the magnetospheric radius. Rotation of the emission region with the star creates an almost sinusoidal pulse profile with an rms amplitude of 3.3 per cent. The energy-dependent soft phase lags can be modelled by two pulsating components shifted in phase, which is naturally explained by a different character of emission of the optically thick spot and optically thin shock combined with the action of the Doppler boosting. The observed variability amplitude constrains the hotspot to lie within 3°–4° of the rotational pole. We estimate the inner radius of the optically thick accreting disc to be about 40 km. In that case, the absence of emission from the antipodal spot, which can be blocked by the accretion disc, gives the inclination of the system as ≳70°.     

Parallel tracks in infrared versus X-ray emission in black hole X-ray transient outbursts: a hysteresis effect?

We report the discovery of a new hysteresis effect in black hole X-ray binary state transitions, that of the near-infrared (NIR) flux (which most likely originates in the jets) versus X-ray flux. We find, looking at existing data sets, that the IR emission of black hole X-ray transients appears to be weaker in the low/hard state rise of an outburst than the low/hard state decline of an outburst at a given X-ray luminosity. We discuss how this effect may be caused by a shift in the radiative efficiency of the inflowing or outflowing matter, or variations in the disc viscosity or the spectrum/power of the jet. In addition we show that there is a correlation (in slope but not in normalization) between IR and X-ray luminosities on the rise and decline, for all three low-mass black hole X-ray binaries with well-sampled IR and X-ray coverage:   L _NIR∝ L ^0.5–0.7_X  . In the high/soft state this slope is much shallower;   L _NIR∝ L ^0.1–0.2_X  , and we find that the NIR emission in this state is most likely dominated by the viscously heated (as opposed to X-ray heated) accretion disc in all three sources.     

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.     

Indirect imaging of the accretion stream in eclipsing polars – III. HU Aquarii low state

We apply our technique for indirect imaging of the accretion stream to the polar HU Aqr, using eclipse profiles observed when the system was in a low-accretion state. The eclipse profile is different from that in the high state, and more variable from cycle to cycle. We find that the stream maps are brightest near the white dwarf and there is no significant brightening in the threading region. In the low state the stream threads on to the magnetic field closer to the L₁ point than in the high state, with a footpoint of the accreting field line at high latitude. We then produce maps of the accretion region from polarimetry using Stokes imaging. These show that the majority of the accretion occurs near the equator. The difference between the maps may be explained if most of the stream material is not emitting significantly in the low state. If so, neither the stream eclipse mapping nor Doppler tomography techniques will trace the bulk of the accretion flow between the two stars.     

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.