A spectroscopic study of V603 Aquilae: stellar parameters and continuumline variations

A time-resolved spectroscopic study of V603 Aql (Nova Aquilae 1918) is presented. An orbital period of P _orb=01385±00002, consistent with previous results, and a radial velocity semi-amplitude of K =20±3 km s¹ are obtained from the radial velocity variations of the H emission line. Similar K values are also found in H , H , and He  i emission lines. Using the measured FWHM of the H line and assuming that the derived semi-amplitude is that of the white dwarf, we deduce a most likely mass ratio of q =0.24±0.05 and stellar masses of M ₂=0.29±0.04 M and M ₁=1.2±0.2 M for the secondary and primary (the white dwarf) star, respectively. The dynamical solution also indicates a very low orbital inclination, i =13°±2°. We find that the continuum and line variations are modulated with both the positive and the negative superhump periods, indicating that they arise from similar regions of the accretion disc. Moreover, we find, for the first time from spectroscopy, evidence of negative superhumps in addition to the positive superhumps. Positive superhumps are explained within the disc instability model as caused by an eccentric disc surrounding the white dwarf, which is precessing (apsidal advance) because of tidal instabilities, causing the observed positive superhumps. A nodal precession in the accretion disc is currently believed to be the cause of the observed negative superhumps. The low value of q is consistent with the expected value for systems that show superhumps, in accordance with the eccentric disc model. We find no evidence of periodicity associated with the spin period.     

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.     

Superhumps in V348 Pup

The eclipsing nova-like cataclysmic variable star V348 Pup exhibits a persistent luminosity modulation with a period 6 per cent longer than its 2.44-h orbital period ( P _orb). This has been interpreted as a 'positive superhump' resulting from a slowly precessing non-axisymmetric accretion disc gravitationally interacting with the secondary. We find a clear modulation of mid-eclipse times on the superhump period, which agrees well with the predictions of a simple precessing eccentric disc model. Our modelling shows that the disc light centre is on the far side of the disc from the donor star when the superhump reaches maximum light. This phasing suggests a link between superhumps in V348 Pup and late superhumps in SU UMa systems. Modelling of the full light curve and maximum entropy eclipse mapping both show that the disc emission is concentrated closer to the white dwarf at superhump maximum than at superhump minimum. We detect additional signals consistent with the beat periods between the implied disc precession period and both and     

The dynamics of eccentric accretion discs in superhump systems

We have applied an eccentric accretion disc theory in simplified form to the case of an accretion disc in a binary system, where the disc contains the 3:1 Lindblad resonance. This is relevant to the case of superhumps in SU Ursae Majoris cataclysmic variables and other systems, where it is thought that this resonance leads to growth of eccentricity and a modulation in the light curve due to the interaction of a precessing eccentric disc with tidal stresses. A single differential equation is formulated which describes the propagation, resonant excitation and viscous damping of eccentricity. The theory is first worked out in the simple case of a narrow ring and leads to the conclusion that the eccentricity distribution is locally suppressed by the presence of the resonance, creating a dip in the eccentricity at the resonant radius. Application of this theory to the superhump case confirms this conclusion and produces a more accurate expression for the precession rate of the disc than has been previously accomplished with simple dynamical estimates.     

Eccentric discs in binaries with intermediate mass ratios: superhumps in the VY Sculptoris stars

We investigate the role of the eccentric disc resonance in systems with mass ratios q ≳1/4, and demonstrate the effects that changes in the mass flux from the secondary star have upon the disc radius and structure. The addition of material with low specific angular momentum to the outer edge of a disc restricts that disc radially. Should the mass flux from the secondary be reduced, it is possible for the disc in a system with mass ratio as large as 1/3 to expand to the 3:1 eccentric inner Lindblad resonance and for superhumps to be excited.     

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.     

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.     

Echoes from an irradiated disc in GRO J1655–40

Oscillations observed in the light curve of Nova V1974 Cygni 1992 since the summer of 1994 have been interpreted as permanent superhumps. From simple calculations based on the tidal disc instability model of Osaki, and assuming that the accretion disc is the dominant optical source in the binary system, we predict that the nova will evolve to become an SU UMa system as its brightness declines from its present value by another 2–3 mag. Linear extrapolation of its current rate of fading (in magnitude units) puts the time of this phase transition within the next 2–4 yr. Alternatively, the brightness decline will stop before the nova reaches that level, and the system will continue to show permanent superhumps in its light curve. It will then be similar to two other old novae, V603 Aql and CP Pup, which still display the permanent superhump phenomenon 80 and 56 yr, respectively, after their eruptions. We suggest that non-magnetic novae with short orbital periods could be progenitors of permanent superhump systems.     

A normal and superoutburst study of the eclipsing SU UMa star: DV Ursae Majoris

We report on time-resolved photometry carried out during the 1995 short outburst and the 1997 long outburst in the eclipsing dwarf nova DV UMa. The revised orbital period is 0.0858526172 (67) d. We detected gigantic superhumps with an amplitude of ∼0.6 mag in the mid-phase of the 1997 outburst, revealing the SU UMa nature of DV UMa. The superhump period is 0.0887 (4) d. The superhumps became less clear during the late phase of the superoutburst, and we found two possible periods of 0.0885 (15) and 0.0764 (15). During both outbursts, the eclipse was wide and shallow near the maximum, and then became narrower and deeper, which is qualitatively well explained by the current disc instability theory.     

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.     

Eccentric instability of a two-dimensional fluid disc with a central massive object

In this paper we investigate, by linear modal analysis, the one-armed dynamical instability of a two-dimensional fluid disc that has a massive object at its centre. The model of the disc is chosen to avoid the artificial instabilities that originate from the unrealistic disc configurations that have been adopted in previous studies. We find a one-armed instability for which the central massive object is displaced from the centre, which is generally called the 'eccentric instability'. However, to excite the eccentric instability, the mass of the central object should be appreciably smaller than that of the disc, and this mass ratio is far smaller than what was originally proposed. The instability shown in this paper is likely to be excited in a stellar system with a central massive object, e.g. a galactic nucleus harbouring a massive black hole, and further studies are desirable via techniques such as numerical simulations.     

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).     

Detection of superhumps in the VY Scl-type nova-like variable KR Aur

We report on the detection of negative superhumps in KR Aur, a typical VY Scl star. The observations were obtained with a multi-channel photometer over 107 h. The analysis of the data clearly revealed brightness variations with a period of 3.771 (±0.005) h. This is 3.5 per cent shorter than P _orb, suggesting that the observed oscillation is a negative superhump. Negative superhumps in VY Scl stars are widespread. The discovery of powerful soft X-rays from V751 Cyg suggests that VY Scl stars may contain white dwarfs, on to which nuclear burning of the accreted material occurs. If this suspicion is correct, it is possible that the powerful radiation emerging from the white dwarf may cause a tilt of the accretion disc to the orbital plane, and its retrograde precession may produce the negative superhumps seen in VY Scl 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.     

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.     

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.     

Balmer absorption variability and the disc wind of V592 Cassiopeiae

We present results from time-series optical spectroscopy of the low-inclination, nova-like cataclysmic variable, V592 Cas. The data span the wavelength range from ∼4000 to 5000  Å, and include Balmer lines (Hβ to Hε) as well as He i and He ii . The Balmer lines are generally characterized by shallow absorption troughs with superimposed narrower central emission components. The absorption troughs are variable on time-scales of at least ∼20 min, but the fluctuations are asymmetric such that they are more dominant on the blueward side of the profile out to ∼−2000 km s⁻¹. Fourier analysis reveals modulation periods for emission radial velocities and absorption changes of ∼0.114 and 0.103 d; we support earlier suggestions that the latter is likely to be 1-day alias of the former, which is essentially the system orbital period. However, there is a phase lag of ∼0.3 between absorption and emission changes. The Balmer absorption changes are compared to the predicted behaviour caused by the (stream overflowing) SW Sex phenomenon. There are several discrepancies in this scenario, including the highly asymmetric nature of blueward changes, the velocity amplitudes, and the phase relation to the emission changes. The role of a disc wind in V592 Cas is also discussed; in particular, with this data set, we cannot rule in or rule out the possibility that axisymmetry of the outflow is broken because it is seated on a warped/tilted inner disc, which is implied by previous evidence for negative superhumps in this system.     

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.     

Thermal stability and nova cycles in permanent superhump systems

Archival data on permanent superhump systems are compiled to test the thermal stability of their accretion discs. We find that their discs are almost certainly thermally stable as expected. This result confirms Osaki's suggestion that permanent superhump systems form a new subclass of cataclysmic variables (CVs), with relatively short orbital periods and high mass-transfer rates. We note that if the high accretion rates estimated in permanent superhump systems represent their mean secular values, then their mass-transfer rates cannot be explained by gravitational radiation, therefore, either magnetic braking should be extrapolated to systems below the period gap or they must have mass-transfer cycles. Alternatively, a new mechanism that removes angular momentum from CVs below the gap should be invoked.
We suggest applying the nova cycle scenarios offered for systems above the period gap to the short orbital period CVs. Permanent superhumps have been observed in the two non-magnetic ex-novae with binary periods below the gap. Their post-nova magnitudes are brighter than their pre-outburst values. In one case (V1974 Cyg) it has been demonstrated that the pre-nova should have been a regular SU UMa system. Thus, it is the first nova whose accretion disc was observed to change its thermal stability. If the superhumps in this system indicate persistent high mass-transfer rates rather than a temporary change induced by irradiation from the hot post-nova white dwarf, it is the first direct evidence for mass-transfer cycles in CVs. The proposed cycles are driven by the nova eruption.     

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.