Two-temperature accretion flows in magnetic cataclysmic variables: structures of post-shock emission regions and X-ray spectroscopy

We use a two-temperature hydrodynamical formulation to determine the temperature and density structures of the post-shock accretion flows in magnetic cataclysmic variables (mCVs) and calculate the corresponding X-ray spectra. The effects of two-temperature flows are significant for systems with a massive white dwarf and a strong white-dwarf magnetic field. Our calculations show that two-temperature flows predict harder keV spectra than one-temperature flows for the same white-dwarf mass and magnetic field. This result is insensitive to whether the electrons and ions have equal temperature at the shock, but depends on the electron–ion exchange rate, relative to the rate of radiative loss along the flow. White-dwarf masses obtained by fitting the X-ray spectra of mCVs using hydrodynamic models including the two-temperature effects will be lower than those obtained using single-temperature models. The bias is more severe for systems with a massive white dwarf.     

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.     

Perturbative analysis of two-temperature radiative shocks with multiple cooling processes

The structure of the hot downstream region below a radiative accretion shock, such as that of an accreting compact object, may oscillate because of a global thermal instability. The oscillatory behaviour depends on the functional forms of the cooling processes, the energy exchanges of electrons and ions in the shock-heated matter, and the boundary conditions. We analyse the stability of a shock with unequal electron and ion temperatures, where the cooling consists of thermal bremsstrahlung radiation which promotes instability, plus a competing process which tends to stabilize the shock. The effect of transverse perturbations is considered also. As an illustration, we study the special case in which the stabilizing cooling process is of order 3/20 in density and 5/2 in temperature, which is an approximation for the effects of cyclotron cooling in magnetic cataclysmic variables. We vary the efficiency of the second cooling process, the strength of the electron–ion exchange and the ratio of electron and ion pressures at the shock, to examine particular effects on the stability properties and frequencies of oscillation modes.     

Stability analyses of two-temperature radiative shocks: formulation, eigenfunctions, luminosity response and boundary conditions

We present a general formulation for stability analyses of radiative shocks with multiple cooling processes, longitudinal and transverse perturbations, and unequal electron and ion temperatures. Using the accretion shocks of magnetic cataclysmic variables as an illustrative application, we investigate the shock instabilities by examining the eigenfunctions of the perturbed hydrodynamic variables. We also investigate the effects of varying the condition at the lower boundary of the post-shock flow from a zero-velocity fixed wall to several alternative types of boundaries involving the perturbed hydrodynamic variables, and the variations of the emission from the post-shock flow under different modes of oscillations. We found that the stability properties for flow with a stationary-wall lower boundary are not significantly affected by perturbing the lower boundary condition, and they are determined mainly by the energy-transport processes. Moreover, there is no obvious correlation between the amplitude or phase of the luminosity response and the stability properties of the system. Stability of the system can, however, be modified in the presence of transverse perturbation. The luminosity responses are also altered by transverse perturbation.     

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.     

Two- and three-dimensional numerical simulations of accretion discs in a close binary system

We perform 2D and 3D numerical simulations of an accretion disc in a close binary system using the simplified flux vector splitting (SFS) finite volume method. In our calculations, the gas is assumed to be ideal with γ =1.01, 1.05, 1.1 and 1.2 . The mass ratio of the mass-losing star to the mass-accreting star is unity. Our results show that spiral shocks are formed on the accretion disc in all cases. In 2D calculations we find that the smaller γ is, the more tightly the spiral winds. We observe this trend in 3D calculations as well in a somewhat weaker sense. Mach numbers in our discs are less than 10. These values are lower than the values in observed accretion discs in close binary systems.
Recently, Steeghs, Harlaftis & Horne found the first convincing evidence for spiral structure in the accretion disc of the eclipsing dwarf nova binary IP Pegasi, using the technique known as Doppler tomography. Although the Mach numbers in present calculations are rather low, we may claim that the spiral structure that we discovered in earlier numerical simulations is now found observationally.     

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.     

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.     

Oscillatory instability of radiative shocks with multiple cooling processes

The stand-off shock formed in the accretion flow on to a stationary wall, such as the surface of a white dwarf, may be thermally unstable, depending on the cooling processes which dominate the post-shock flow. Some processes lead to instability, while others tend to stabilize the shock. We consider competition between the destabilizing influence of thermal bremsstrahlung cooling, and a stabilizing process which is a power law in density and temperature. Cyclotron cooling and processes which are of order 1, 3/2 and 2 in density are considered. The relative efficiency and power-law indices of the second mechanism are varied, and particular effects on the stability properties and frequencies of oscillation modes are examined.     

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

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

A non-linear theory of vertical resonances in accretion discs

An important and widely neglected aspect of the interaction between an accretion disc and a massive companion with a coplanar orbit is the vertical component of the tidal force. As shown by Lubow, the response of the disc to vertical forcing is resonant at certain radii, at which a localized torque is exerted, and from which a compressive wave (p mode) may be emitted. Although these vertical resonances are weaker than the corresponding Lindblad resonances, the   m =2  inner vertical resonance in a binary star is typically located within the tidal truncation radius of a circumstellar disc.
In this paper I develop a general theory of vertical resonances, allowing for non-linearity of the response, and dissipation by radiative damping and turbulent viscosity. The problem is reduced to a universal, non-linear ordinary differential equation with two real parameters. Solutions of the complex non-linear Airy equation are presented to illustrate the non-linear saturation of the resonance and the effects of dissipation. It is argued that the   m =2  inner vertical resonance is unlikely to truncate the disc in cataclysmic variable stars, but contributes to angular momentum transport and produces a potentially observable non-axisymmetric structure.     

Corotational damping of discoseismic c modes in black hole accretion discs

Discoseismic c modes in accretion discs have been invoked to explain low-frequency variabilities observed in black hole X-ray binaries. These modes are trapped in the innermost region of the disc and have frequencies much lower than the rotation frequency at the disc inner radius. We show that because the trapped waves can tunnel through the evanescent barrier to the corotational wave zone, the c modes are damped due to wave absorption at the corotation resonance. We calculate the corotational damping rates of various c modes using the Wentzel-Kramers-Brillouin (WKB) approximation. The damping rate varies widely depending on the mode frequency, the black hole spin parameter and the disc sound speed, and is generally much less than 10 per cent of the mode frequency. A sufficiently strong excitation mechanism is needed to overcome this corotational damping and make the mode observable.     

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.     

GRB060602B = Swift J1749.4−2807: an unusual transiently accreting neutron-star X-ray binary

We present an analysis of the Swift Burst Alert Telescope (BAT) and X-ray telescope (XRT) data of GRB060602B, which is most likely an accreting neutron star in a binary system and not a gamma-ray burst. Our analysis shows that the BAT burst spectrum is consistent with a thermonuclear flash (type I X-ray burst) from the surface of an accreting neutron star in a binary system. The X-ray binary nature is further confirmed by the report of a detection of a faint point source at the position of the XRT counterpart of the burst in archival XMM–Newton data approximately six year before the burst and in more recent XMM–Newton data obtained at the end of 2006 September (nearly four months after the burst). Since the source is very likely not a gamma-ray burst, we rename the source Swift J1749.4−2807, based on the Swift /BAT discovery coordinates. Using the BAT data of the type I X-ray burst, we determined that the source is at most at a distance of  6.7 ± 1.3 kpc  . For a transiently accreting X-ray binary, its soft X-ray behaviour is atypical: its 2–10 keV X-ray luminosity (as measured using the Swift /XRT data) decreased by nearly three orders of magnitude in about 1 day, much faster than what is usually seen for X-ray transients. If the earlier phases of the outburst also evolved this rapidly, then many similar systems might remain undiscovered because the X-rays are difficult to detect and the type I X-ray bursts might be missed by all the sky surveying instruments. This source might be part of a class of very fast transient low-mass X-ray binary systems of which there may be a significant population in our Galaxy.