Ultraviolet spectroscopic study of EU UMa and ST LMi from HST and IUE observations

We present the first orbit-resolved ultraviolet spectroscopic observations of the two polar systems EU UMa and ST LMi obtained with the Hubble Space Telescope Faint Object Spectrograph (HST FOS) and International Ultraviolet Explorer (IUE) during their intermediate- and low-luminosity states in the period between 1982 and 2003. Different line profiles of the two systems showing variations of line fluxes at different orbital phases are presented. This paper focuses on the C IV emission line at 1550 Å produced in the accretion stream, presenting calculations of spectral line fluxes, ultraviolet luminosities, and accretion rates for the two systems. Our analysis of the spectroscopic data reveals changes with orbital phase of the emission line profiles that correspond to the light curve variations of both EU UMa and ST LMi in the optical and infrared bands. The variations of line fluxes are attributed to variations of both density and temperature as a result of a changing rate of mass transfer from the secondary star to the white dwarf. The ultraviolet luminosity and accretion rate of EU UMa are smaller than the ultraviolet luminosity and accretion rate of ST LMi.     

Ionization structure in accretion shocks with a composite cooling function

We have investigated the ionization structure of the post-shock regions of magnetic cataclysmic variables, using an analytic density and temperature structure model in which effects caused by bremsstrahlung and cyclotron cooling are considered. We find that in the majority of the shock-heated region where H- and He-like lines of the heavy elements are emitted, the collisional-ionization and corona-condition approximations are justified. We have calculated the line emissivity and ionization profiles for iron as a function of height within the post-shock flow. For low-mass white dwarfs, line emission takes place near the shock. For high-mass white dwarfs, most of the line emission takes place in regions well below the shock and hence it is less sensitive to the shock temperature. Thus, the line ratios are useful to determine the white dwarf masses for the low-mass white dwarfs, but the method is less reliable when the white dwarfs are massive. Line spectra can, however, be used to map the hydrodynamic structure of the post-shock accretion flow.     

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.     

Return mapping of phases and the analysis of the gravitational clustering hierarchy

We present spectroscopic and high-speed photometric data of the eclipsing polar V895 Cen. We find that the eclipsed component is consistent with it being the accretion regions on the white dwarf. This is in contrast to Stobie et al. who concluded that the eclipsed component was not the white dwarf. Further, we find no evidence for an accretion disc in our data. From our Doppler tomography results, we find that the white dwarf has   M ≳0.7 M_⊙  . Our indirect imaging of the accretion stream suggests that the stream is brightest close to the white dwarf. When we observed V895 Cen in its highest accretion state, emission was concentrated along field lines leading to the upper pole. There is no evidence for enhanced emission at the magnetic coupling region.     

The remarkable rapid X-ray, ultraviolet, optical and infrared variability in the black hole XTE J1118+480

The transient black-hole binary XTE J1118+480 exhibited dramatic rapid variability at all wavelengths which were suitably observed during its 2000 April–July outburst. We examine time-resolved X-ray, ultraviolet, optical and infrared data spanning the plateau phase of the outburst. We find that both X-ray and infrared bands show large amplitude variability. The ultraviolet and optical variability is more subdued, but clearly correlated with that seen in the X-rays. The ultraviolet, at least, appears to be dominated by the continuum, although the lines are also variable. Using the X-ray variations as a reference point, we find that the ultraviolet (UV) variability at long wavelengths occurs later than that at short wavelengths. Uncertainty in the Hubble Space Telescope timing prohibits a determination of the absolute lag with respect to the X-rays, however. The transfer function is clearly not a delta-function, exhibiting significant repeatable structure. For the main signal we can rule out an origin in reprocessing on the companion star – the lack of variation in the lags is not consistent with this, given a relatively high orbital inclination. Weak reprocessing from the disc and/or companion star may be present, but is not required, and another component must dominate the variability. This could be variable synchrotron emission correlated with X-ray variability, consistent with our earlier interpretation of the infrared (IR) flux as due to synchrotron emission rather than thermal disc emission. In fact, the broad-band energy distribution of the variability from IR to X-rays is consistent with expectations of optically thin synchrotron emission. We also follow the evolution of the low-frequency quasi-periodic oscillation in X-rays, UV, and optical. Its properties at all wavelengths are similar, indicating a common origin.     

White dwarf masses in magnetic cataclysmic variables: multi-temperature fits to Ginga data

One method of obtaining the mass of the white dwarf in magnetic cataclysmic variables (mCVs) is through their hard X-ray spectra. However, previous mass estimates using this method give lower limits because the temperature of the plasma in the post-shock region (where the hard X-rays are emitted) is lower than the temperature of the shock itself. In AM Her systems, the additional cooling of the post-shock plasma by cyclotron emission will further lower the derived mass. Here we present estimates of the masses of the white dwarf in 13 mCVs derived using Ginga data and a model in which X-rays are emitted from a multi-temperature emission region with the appropriate temperature and density profile. We include in the model reflection from the surface of the white dwarf and a partially ionized absorber. We are able to achieve good fits to the data. We compare the derived masses with previous estimates and the masses for larger samples of isolated white dwarfs and those in CVs.     

Self-regulated gravitational accretion in protostellar discs

We present a numerical model for the evolution of a protostellar disc that has formed self-consistently from the collapse of a molecular cloud core. The global evolution of the disc is followed for several million years after its formation. The capture of a wide range of spatial and temporal scales is made possible by use of the thin-disc approximation. We focus on the role of gravitational torques in transporting mass inward and angular momentum outward during different evolutionary phases of a protostellar disc with disc-to-star mass ratio of order 0.1. In the early phase, when the infall of matter from the surrounding envelope is substantial, mass is transported inward by the gravitational torques from spiral arms that are a manifestation of the envelope-induced gravitational instability in the disc. In the late phase, when the gas reservoir of the envelope is depleted, the distinct spiral structure is replaced by ongoing irregular non-axisymmetric density perturbations. The amplitude of these density perturbations decreases with time, though this process is moderated by swing amplification aided by the existence of the disc's sharp outer edge. Our global modelling of the protostellar disc reveals that there is typically a residual non-zero gravitational torque from these density perturbations, i.e. their effects do not exactly cancel out in each region. In particular, the net gravitational torque in the inner disc tends to be negative during first several million years of the evolution, while the outer disc has a net positive gravitational torque. Our global model of a self-consistently formed disc shows that it is also self-regulated in the late phase, so that it is near the Toomre stability limit, with a near-uniform Toomre parameter Q ≈ 1.5–2.0. Since the disc also has near-Keplerian rotation, and comparatively weak temperature variation, it maintains a near-power-law surface density profile proportional to r ^−3/2.     

X-ray polarization signatures of Compton scattering in magnetic cataclysmic variables

Compton scattering within the accretion column of magnetic cataclysmic variables (mCVs) can induce a net polarization in the X-ray emission. We investigate this process using Monte Carlo simulations and find that significant polarization can arise as a result of the stratified flow structure in the shock-ionized column. We find that the degree of linear polarization can reach levels up to ∼8 per cent for systems with high accretion rates and low white dwarf masses, when viewed at large inclination angles with respect to the accretion column axis. These levels are substantially higher than previously predicted estimates using an accretion column model with uniform density and temperature. We also find that for systems with a relatively low-mass white dwarf accreting at a high accretion rate, the polarization properties may be insensitive to the magnetic field, since most of the scattering occurs at the base of the accretion column where the density structure is determined mainly by bremsstrahlung cooling instead of cyclotron cooling.     

An eclipse of the X-ray flux from the dwarf nova OY Carinae in quiescence

We present a phase-resolved ROSAT HRI X-ray light curve of the dwarf nova OY Car in quiescence. The X-ray flux is eclipsed at the same time as the optical eclipse of the primary, and the region of X-ray emission is comparable in size to the white dwarf. We use subsequent optical observations to update the orbital ephemeris of the system.     

Emission-line oscillations in the dwarf nova V2051 Ophiuchi

We have detected coherent oscillations, at multiple frequencies, in the line and continuum emission of the eclipsing dwarf nova V2051 Ophiuchi using the 10-m Keck II telescope. Our own novel data acquisition system allowed us to obtain very fast spectroscopy using a continuous readout of the CCD on the LRIS spectrograph. This is the first time that dwarf nova oscillations have been detected and resolved in the emission lines. The accretion disc is highly asymmetric with a stronger contribution from the blueshifted side of the disc during our observations. The disc extends from close to the white dwarf out to the outer regions of the primary Roche lobe.
Continuum oscillations at 56.12 s and its first harmonic at 28.06 s are most likely to originate on the surface of a spinning white dwarf with the fundamental period corresponding to the spin period. Balmer and helium emission lines oscillate with a period of 29.77 s at a mean amplitude of 1.9 per cent. The line kinematics and the eclipse constraints indicate an origin in the accretion disc at a radius of 12±2 R _wd. The amplitude of the emission-line oscillation modulates (0–4 per cent) at a period of 488 s, corresponding to the Kepler period at R =12 R _wd. This modulation is caused by the beating between the white dwarf spin and the orbital motion in the disc.
The observed emission-line oscillations cannot be explained by a truncated disc as in the intermediate polars. The observations suggest a non-axisymmetric bulge in the disc, orbiting at 12 R _wd, is required. The close correspondence between the location of the oscillations and the circularization radius of the system suggests that stream overflow effects may be of relevance.     

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

A source of high-velocity white dwarfs

We investigate whether the recently observed population of high-velocity white dwarfs can be derived from a population of binaries residing initially within the thin disc of the Galaxy. In particular, we consider binaries where the primary is sufficiently massive to explode as a Type II supernova. A large fraction of such binaries are broken up when the primary then explodes as a supernova, owing to the combined effects of the mass loss from the primary and the kick received by the neutron star on its formation. For binaries where the primary evolves to fill its Roche lobe, mass transfer from the primary leads to the onset of a common envelope phase during which the secondary and the core of the primary spiral together as the envelope is ejected. Such binaries are the progenitors of X-ray binaries if they are not broken up when the primary explodes. For those systems that are broken up, a large number of the secondaries receive kick velocities ∼100–200 km s⁻¹ and subsequently evolve into white dwarfs. We compute trajectories within the Galactic potential for this population of stars and relate the birth rate of these stars over the entire Galaxy to those seen locally with high velocities relative to the local standard of rest (LSR) . We show that for a reasonable set of assumptions concerning the Galactic supernova rate and the binary population, our model produces a local number density of high-velocity white dwarfs compatible with that inferred from observations. We therefore propose that a population of white dwarfs originating in the thin disc may make a significant contribution to the observed population of high-velocity white dwarfs.     

Optical spectroscopy and Doppler tomography of Cygnus X-2

We present phase resolved optical spectroscopy and Doppler tomography of V1341 Cygni, the optical counterpart to the neutron star low-mass X-ray binary (LMXB) Cygnus X-2 (Cyg X-2). We derive a radial velocity (RV) curve for the secondary star, finding a projected RV semi-amplitude of   K ₂= 79 ± 3 km s⁻¹  , leading to a mass function of  0.51 ± 0.06 M_⊙, ∼30  per cent lower than the previous estimate. We tentatively attribute the lower value of K ₂ (compared to that obtained by other authors) to variations in the X-ray irradiation of the secondary star at different epochs of observations. The limited phase coverage and/or longer timebase of previous observations may also contribute to the difference in K ₂. Our value for the mass function implies a primary mass of  1.5 ± 0.3 M_⊙  , somewhat lower than previous dynamical estimates, but consistent with the value found by analysis of type-I X-ray bursts from this system. Our Doppler tomography of the broad He  ii λ4686 line reveals that most of the emission from this line is produced on the irradiated face of the donor star, with little emission from the accretion disc. In contrast, the Doppler tomogram of the N  iii λ4640.64 Bowen blend line shows bright emission from near the gas stream/accretion disc impact region, with fainter emission from the gas stream and secondary star. This is the first LMXB for which the Bowen blend is dominated by emission from the gas stream/accretion disc impact region, without comparable emission from the secondary star. This has implications for the interpretation of Bowen blend Doppler tomograms of other LMXBs for which the ephemeris may not be accurately known.     

The lower boundary of the accretion column in magnetic cataclysmic variables

Using a parametrized function for the mass loss at the base of the post-shock region, we have constructed a formulation for magnetically confined accretion flows which avoids singularities, such as the infinity in density, at the base associated with all previous formulations. With the further inclusion of a term allowing for the heat input into the base from the accreting white dwarf, we are also able to obtain the hydrodynamic variables to match the conditions in the stellar atmosphere. (We do not, however, carry out a mutually consistent analysis for the match.) Changes to the emitted X-ray spectra are negligible unless the thickness of mass leakage region at the base approaches or exceeds one per cent of the height of the post-shock region. In this case the predicted spectra from higher-mass white dwarfs will be harder, and fits to X-ray data will predict lower white dwarf masses than previous formulations.     

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.     

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.     

High-speed energy-resolved STJ photometry of the eclipsing binary UZ For

We present high-time-resolution optical photometry of the eclipsing binary UZ For using a superconducting tunnel junction (STJ) device, a photon-counting array detector with intrinsic energy resolution. Three eclipses of the ∼18-mag 126.5-min orbital binary were observed using a 6×6 array of tantalum STJs at the 4.2-m William Herschel Telescope on La Palma. The detector presently provides individual photon arrival-time accuracy to about 5 μs, and a wavelength resolution of about 60 nm at 500 nm, with each array element capable of counting up to ∼5000 photon s⁻¹. The data allow us to place accurate constraints on the accretion geometry from our time- and spectrally resolved monitoring, especially of the eclipse ingress and egress. We find that there are two small accretion regions, located close to the poles of the white dwarf. The positions of these are accurately constrained, and show little movement from eclipse to eclipse, even over a number of years. The colour of the emission from the two regions appears similar, although their X-ray properties are known to be significantly different: we argue that the usual accretion shock may be absent at the non-X-ray-emitting region, and instead the flow here interacts directly with the white dwarf surface; alternatively, a special grazing occultation of this region is required. There is no evidence for any quasi-periodic oscillations on time-scales of the order of seconds, consistent with relatively stable cyclotron cooling in each accretion region.     

The mass of the white dwarf in the recurrent nova U Scorpii

We present spectroscopy of the eclipsing recurrent nova U Sco. The radial velocity semi-amplitude of the primary star was found to be     from the motion of the wings of the He  ii λ 4686-Å emission line. By detecting weak absorption features from the secondary star, we find its radial velocity semi-amplitude to be     . From these parameters, we obtain a mass of     for the white dwarf primary star and a mass of     for the secondary star. The radius of the secondary is calculated to be     , confirming that it is evolved. The inclination of the system is calculated to be     , consistent with the deep eclipse seen in the light-curves. The helium emission lines are double-peaked, with the blueshifted regions of the disc being eclipsed prior to the redshifted regions, clearly indicating the presence of an accretion disc. The high mass of the white dwarf is consistent with the thermonuclear runaway model of recurrent nova outbursts, and confirms that U Sco is the best Type Ia supernova progenitor currently known. We predict that U Sco is likely to explode within ∼700 000 yr.     

Indirect imaging of the accretion stream in eclipsing polars – II. HU Aquarii

We apply our technique for indirect imaging of the accretion stream to the polar HU Aquarii, using eclipse profiles observed when the system was in a high accretion state. The accretion stream is relatively luminous, contributing as much as the accretion region on the white dwarf, or more, to the overall system brightness. We model the eclipse profiles using a model stream consisting of a ballistic trajectory from the L1 point followed by a magnetically channelled trajectory that follows a dipole field line out of the orbital plane. We perform model fits using two geometries: a stream that accretes on to both footpoints of the field line, and a stream that accretes only on to the footpoint of the field line above the orbital plane. The stream images indicate that the distribution of emission along the stream is not a simple function of the radial distance from the white dwarf. The stream is redirected by the magnetic field of the white dwarf at a distance 1.0–1.3×10¹⁰ cm from the white dwarf; this implies a mass transfer rate in the range 8–76×10¹⁶ g s⁻¹. The absorption dips in the light curve indicate that the magnetically entrained part of the stream moves from 42° to 48° from the line of centres over the three nights of observation. This is in close agreement with the results of the one-footpoint models, suggesting that this is the more appropriate geometry for these data. The stream images show that, in almost all sections of the stream, the flux peaks in B and is successively fainter in U , V and R .     

Polarimetry of the eclipsing polar RX J0929.1−2404

We report polarimetric, spectropolarimetric and photometric observations of the eclipsing ROSAT cataclysmic variable RX J0929.1−2404, which confirm that the system is a new polar (AM Herculis system). This brings the number of eclipsing polars to nine, with RX J0929.1−2404 being only the third such system above the period gap. Circular polarization variations from ∼−20 to 10 per cent are seen over the 3.39-h orbital period, with a minimum around the time of eclipse. The photopolarimetric data were modelled using arc-shaped cyclotron emission regions in a centred dipole geometry. Results imply that RX J0929.1−2404 is a 'two-pole' system, with one emission region partially visible at all orbital phases. Spectropolarimetry observations show some evidence for the presence of cyclotron humps in the continuum, with spacings consistent with a magnetic field strength of ∼20 MG. Photometry of the eclipses provides information on the size of the emission region, which is consistent with a hotspot on the surface of the white dwarf. The eclipse duration implies an inclination in the range 70°≲ i ≲78°.