Estimation of plasma parameters in an accretion column near the surface of accreting white dwarfs from their flux variability

We consider the behavior of matter in the accretion column that emerges under accretion in binary systems near the surface of a white dwarf. The plasma heated in a standing shock wave near the white dwarf surface efficiently radiates in the X-ray energy band. We suggest a method for estimating post-shock plasma parameters, such as the density, temperature, and height of the hot zone, from the power spectrum of its X-ray luminosity variability. The method is based on the fact that the flux variability amplitude for the hot region at various Fourier frequencies depends significantly on its cooling time, which is determined by the parameters of the hot zone in the accretion column. This allows the density and temperature of the hot matter to be estimated. We show that the characteristic cooling time can be efficiently determined from the break frequency in the power spectrum of the X-ray flux variability for accreting white dwarfs. The currently available X-ray instruments do not allow such measurements to be made because of an insufficient collecting area, but this will most likely become possible with new-generation large-area X-ray spectrometers.     

Diagnostics of the early explosion phase of a classical nova using its X-ray emission: A model for the X-ray outburst of CI Camelopardalis in 1998

We have computed a spherically symmetric model for the interaction of matter ejected during the outburst of a classical nova with the stellar wind from its optical component. This model is used to describe the intense X-ray outburst (the peak 3–20 keV flux was ～2 Crab) of the binary system CI Camelopardalis in 1998. According to our model, the stellar wind from the optical component heated by a strong shock wave produced when matter is ejected from the white dwarf as the result of a thermonuclear explosion on its surface is the emission source in the standard X-ray band. Comparison of the calculated and observed time dependences of the mean radiation temperature and luminosity of the binary system during its outburst has yielded very important characteristics of the explosion. We have been able to measure the velocity of the ejected matter immediately after the onset of the explosion for the first time: it follows from our model that the ejected matter had a velocity of ～2700 km s^?1 even on 0.1–0.5 day after the outburst onset and it flew with such a velocity for the first 1–1.5 day under an external force, possibly, the radiation pressure from the white dwarf. Subsequently, the matter probably became transparent and began to decelerate. The time dependence of the mean radiation temperature at late expansion phases has allowed us to estimate the mass of the ejected matter, ～10^?7–10^?6 M _⊙. The mass loss rate in the stellar wind required to explain the observed peak luminosity of the binary system during its outburst has been estimated to be \(\dot M\) ～ (1 ? 2) × 10^?6 M _⊙ yr^?1.     

Peculiarities of the accretion flow in the system HL CMa

The properties of the aperiodic luminosity variability for the dwarf novaHLCMa are considered. The variability of the system HL CMa is shown to be suppressed at frequencies above 0.7 × 10^?2 Hz. Different variability suppression mechanisms related to the radiation reprocessing time, partial disk evaporation, and characteristic variability formation time are proposed. It has been found that the variability suppression frequency does not change when the system passes from the quiescent state to the outburst one, suggesting that the accretion flow geometry is invariable. It is concluded from the optical and Xray luminosities of the system that the boundary layer on the white dwarf surface is optically thick in both quiescent and outburst states. The latter implies that the optically thick part of the accretion flow (disk) reaches the white dwarf surface. The accretion rate in the system and the accretion flow geometry and temperature have been estimated from the variability power spectra and spectral characteristics in a wide energy range, from the optical to X-ray ones.     

On the change of the inner boundary of an optically thick accretion disk around white dwarfs using the dwarf nova SS Cyg as an example

We present the results of our studies of the aperiodic optical flux variability for SS Cyg, an accreting binary systemwith a white dwarf. The main set of observational data presented here was obtained with the ANDOR/iXon DU-888 photometer mounted on the RTT-150 telescope, which allowed a record (for CCD photometers) time resolution up to 8 ms to be achieved. The power spectra of the source’s flux variability have revealed that the aperiodic variability contains information about the inner boundary of the optically thick flow in the binary system. We show that the inner boundary of the optically thick accretion disk comes close to the white dwarf surface at the maximum of the source’s bolometric light curve, i.e., at the peak of the instantaneous accretion rate onto the white dwarf, while the optically thick accretion disk is truncated at distances 8.5 × 10⁹ cm ∼10R _WD in the low state. We suggest that the location of the inner boundary of the accretion disk in the binary can be traced by studying the parameters of the power spectra for accreting white dwarfs. In particular, this allows the mass of the accreting object to be estimated.     

ASCA X-ray observations of EX Hya: spin-resolved spectroscopy

We analyse the spectral changes over the spin modulation in the intermediate polar EX Hya using archival ASCA data. We find that the modulation can be modelled as either (1) the effect of occultation of the accretion poles by the limb of the white dwarf, or (2) the effect of phase-dependent photoelectric absorption. We argue, on the basis of the partial X-ray eclipse, that the accretion columns in the system are tall, with shock height ∼ R _wd, and hence that the spin modulation is caused mainly by occultation. We find that the temperature distribution along the accretion shocks is incompatible with the calculations of Aizu, except for a restricted parameter regime with a high M _wd. Hence the material in the shock must cool faster than predicted by theory.     

Relation between the X-ray and optical luminosities in binary systems with accreting nonmagnetic white dwarfs

We investigate the relation between the optical (g-band) and X-ray (0.5–10 keV) luminosities of accreting nonmagnetic white dwarfs. According to the present-day counts of the populations of star systems in our Galaxy, these systems have the highest space density among the close binary systems with white dwarfs. We show that the dependence of the optical luminosity of accreting white dwarfs on their X-ray luminosity forms a fairly narrow one-parameter curve. The typical half-width of this curve does not exceed 0.2–0.3 dex in optical and X-ray luminosities, which is essentially consistent with the amplitude of the aperiodic flux variability for these objects. At X-ray luminosities L _x ～ 10³² erg s^?1 or lower, the optical g-band luminosity of the accretion flow is shown to be related to its X-ray luminosity by a factor ～2–3. At even lower X-ray luminosities (L _x ? 10³⁰ erg s^?1), the contribution from the photosphere of the white dwarf begins to dominate in the optical spectrum of the binary system and its optical brightness does not drop below M _g ～ 13–14. Using the latter fact, we show that in current and planned X-ray sky surveys, the family of accreting nonmagnetic white dwarfs can be completely identified to the distance determined by the sensitivity of an optical sky survey in this region. For the Sloan Digital Sky Survey (SDSS) with a limiting sensitivity m _g ～ 22.5, this distance is ～400–600 pc.     

Multiple X-ray bursts and the model of a spreading layer of accreting matter over the neutron star surface

We report the detection of series of close type I X-ray bursts consisting of two or three events with a recurrence time much shorter than the characteristic (at the observed mean accretion rate) time of matter accumulation needed for a thermonuclear explosion to be initiated on the neutron star surface during the JEM-X/INTEGRAL observations of several X-ray bursters. We show that such series of bursts are naturally explained in the model of a spreading layer of accreting matter over the neutron star surface in the case of a sufficiently high (? ? 1 × 10^?9 M _⊙ yr^?1) accretion rate (corresponding to a mean luminosity L _tot ? 1 × 10³⁷erg s^?1). The existence of triple bursts requires some refinement of the model—the importance of a central ring zone is shown. In the standard model of a spreading layer no infall of matter in this zone is believed to occur.     

X-ray and ultraviolet observations of the dwarf nova VW Hyi in quiescence

We present an analysis of X-ray and ultraviolet (UV) data of the dwarf nova VW Hyi that were obtained with XMM–Newton during the quiescent state. The X-ray spectrum indicates the presence of an optically thin plasma in the boundary layer that cools as it settles on to the white dwarf. The plasma has a continuous temperature distribution that is well described by a power law or a cooling flow model with a maximum temperature of 6–8 keV. We estimate from the X-ray spectrum a boundary layer luminosity of  8 × 10³⁰ erg s^-1  , which is only 20 per cent of the disc luminosity. The rate of accretion on to the white dwarf is  5 × 10⁻¹² M_⊙ yr⁻¹  , about half of the rate in the disc. From the high-resolution X-ray spectra, we estimate that the X-ray emitting part of the boundary layer is rotating with a velocity of 540 km s⁻¹, which is close to the rotation velocity of the white dwarf but is significantly smaller than the Keplerian velocity. We detect a 60-s quasi-periodic oscillation of the X-ray flux, which is likely to be due to the rotation of the boundary layer. The X-ray and the UV flux show strong variability on a time-scale of ∼1500 s. We find that the variability in the two bands is correlated and that the X-ray fluctuations are delayed by ∼100 s. The correlation indicates that the variable UV flux is emitted near the transition region between the disc and the boundary layer and that accretion rate fluctuations in this region are propagated to the X-ray emitting part of the boundary layer within ∼100 s. An orbital modulation of the X-ray flux suggests that the inner accretion disc is tilted with respect to the orbital plane. The elemental abundances in the boundary layer are close to their solar values.     

Capture radius and synchronization of the white dwarf in the unique magnetic cataclysmic system V1432 Aql

Results from two-color VR photometry of the unique cataclysmic magnetic variable star V1432 Aql and a theoretical model of these data are presented. The accuracy is improved by using the “mean-weighted comparison star” method. The derivative of the rotational period is dP/dt = −1.11(±0.016)·10⁻⁸. The characteristic synchronization time for the rotational and orbital motions of the white dwarf is 96.7±1.5 years, in good agreement with theory for the acceleration of an asynchronous propeller owing to the angular momentum of accreting matter. A third type of minimum detected in the light curve is interpreted in terms of the presence of an arc, or ring, rather than an accretion disk. A theoretical model is developed for determining the capture radius of accreted matter by the magnetic field of the white dwarf using the phase difference between the two types of minima associated with the axial rotation. This parameter is estimated to be 16–28 times the radius of the white dwarf for an inclined column model. A dependence of the main characteristics of the system on the mass of the white dwarf is derived which yields better values for the range of this quantity than those determined by indirect methods. For the assumed masses (M₁ = 0.9 M_⊙ and M₂ = 0.3 M_⊙) the estimated accretion rate is ∼7×10⁻¹⁰ M_⊙. It is shown that in a synchronizing polar the contribution to the change in the period by the variation in the angular momentum of the white dwarf is negligible compared to the accretion torque. In the future multicolor monitoring is needed for studying the spin-orbital synchronization and periodic changes in the accretion structure caused by “spinning” of the white dwarf. __________ Translated from Astrofizika, Vol. 50, No. 1, pp. 135–159 (February 2007).     

Broadband Spectrum of the X-ray Binary M33 X-6 from NuSTAR and Swift–XRT Data: An Extragalactic Z-Source?

We present the results of our study of the X-ray spectrum for the source X-6 in the nearby galaxy M33 obtained for the first time at energies above 10 keV from the data of the NuSTAR orbital telescope. The archival Swift–XRT data for energy coverage below 3 keV have been used, which has allowed the spectrum of M33 X-6 to be constructed in the wide energy range 0.3–20 keV. The spectrum of the source is well described by the model of an optically and geometrically thick accretion disk with a maximum temperature of ~2 keV and an inner radius of ~5 cos^?1/2θ km (where >θ is the unknown disk inclination angle with respect to the observer). There is also evidence for the presence of an additional hard component in the spectrum. The X-ray luminosity ofM33 X-6 measured for the first time in the wide energy range 0.3–20 keV is ~2 × 10³⁸ erg s^?1, with the luminosity in the hard 10–20 keV X-ray band being ~10% of the source’s total luminosity. The results obtained suggest that X-6 may be a Z-source, i.e., an X-ray binary with subcritical accretion onto a weakly magnetized neutron star.     

Accretion effects on compact members of binary stars

We explore the change in the period of axial rotation and in the radius of a magnetized compact star in a binary system, induced by the accretion on it of mass with angular momentum from the surface of its non-compact companion. No specific assumption is made concerning the accretion model, and the primary’s interior is described by the Fermi-Dirac statistics for degenerate matter. The rate of change with time of the period and radius is expressed in terms of the compact primary’s physical parameters and total absolute luminosity. The conditions are fully derived under which the above changes can be positive, negative or even vanish. In the case of the millisecond pulsars in binary X-ray sources the predicted values of the period time derivative, depending on the values of the accretion rate and. the absolute luminosity, can be positive or negative—if not vanishing—and they fall absolutely in the range 10⁻²¹ −10⁻¹⁷ ss⁻¹, in good agreement with current observational data. The corresponding rate of change of radius, either positive or negative, fall in the range of 10⁻³ −10⁻¹ cm y⁻¹. Finally, it is proved that the well-known bursters can be explained by thermonuclear flash due to gravitational instability in the accreted matter, but their explanation as a result of direct contraction could be possible only for quite high accretion rates (>10⁻⁷ M ⊙ y⁻¹). This last result indicates that, in contrast to the accretion-induced change in period, which can be of either sign irrespective of the primary’s age, the accretion-induced non-catastrophic contraction is impossible, while according to repent results the contraction in general is possible for young compact objects.     

4U 0115+63 from RXTE and INTEGRAL data: Pulse profile and cyclotron line energy

We analyze the observations of the transient X-ray pulsar 4U 0115+63 with the RXTE and INTEGRAL observatories in a wide X-ray (3–100 keV) energy band during its intense outbursts in 1999 and 2004. The energy of the fundamental harmonic of the cyclotron resonance absorption line near the maximum of the X-ray flux from the source (luminosity range 5 × 10³⁷–2 × 10³⁸ erg s^?1) is ～11 keV. When the pulsar luminosity falls below ～5 × 10³⁷ erg s^?1, the energy of the fundamental harmonic is displaced sharply toward the high energies, up to ～16 keV. Under the assumption of a dipole magnetic field configuration, this change in cyclotron harmonic energy corresponds to a decrease in the height of the emitting region by ～2 km, while other spectral parameters, in particular, the cutoff energy, remain essentially constant. At a luminosity ～7 × 10³⁷ erg s^?1, four almost equidistant cyclotron line harmonics are recorded in the spectrum. This suggests that either the region where the emission originates is compact or the emergent spectrum from different (in height) segments of the accretion column is uniform. We have found significant pulse profile variations with energy, luminosity, and time. In particular, we show that the profile variations from pulse to pulse are not reduced to a simple modulation of the accretion rate specified by external conditions.     

Constraints on accretion column parameters for the intermediate polar LS Pegasi from the power spectrum of its optical light curve

The properties of the hot zone in the accretion flow near the surface of a magnetized white dwarf have been studied. For this purpose, the aperiodic optical variability of LS Peg, one of the brightest intermediate polars in the northern sky, has been investigated. The main radiation of the hot zone, which is then reemitted in the optical band, results from the radiation of an optically thin plasma heated during the passage of the accretion flow of a standing shock. Recently, Semena and Revnivtsev (2012) have shown that the aperiodic variability (flickering) of accreting magnetized white dwarfs should have a characteristic feature in the range of Fourier frequencies corresponding to the plasma cooling time in this hot region. The photometric brightness measurements for LS Peg made with the RTT-150 telescope using a high-speed ANDOR iXon CCD array have allowed the break frequency in the power spectrum of the source’s variability to be constrained. Constraints on the geometry of the accretion column for the white dwarf in LS Peg and on the plasma parameters in it have been obtained.     

Preliminary results of the spectral analysis of Suzaku data of SW Ursae Majoris and BZ Ursae Majoris

Herein, we present the preliminary results of the spectral analysis of the Suzaku data of SW Ursae Majoris (SW UMa) and BZ Ursae Majoris (BZ UMa) which were obtained in their quiescent states. The aim of this study was to analyse 0.5–6.0 keV energy range X-ray photons of the Suzaku satellite and to determine their emission mechanisms and the areas where they were radiated from these two dwarf novae. No study on this energy range for BZ UMa has been reported in the literature. We analysed the Suzaku data in the (0.5–6.0) keV and in (0.5–2.5) keV energy bands of the systems separatelly to search for possible differences between the emission mechanisms in the energy ranges; however, none were found. We obtained the best-fitted spectral models, flux, luminosity, temperature and mass accretion rate values for the systems. The white dwarf radius of BZ UMa, which has not been reported in the literature, was found to be 6.88 × 10⁸ cm. From this study, we concluded that the mass accretion rate values for both systems were lower than the critical mass accretion rate value in the boundary layers of cataclysmic variables, which indicates that the boundary layers are composed of optically thin thermal X-ray emitting gas in their quiescent state; these results support those of previous studies regarding SW UMa and some studies regarding BZ UMa. Another result was that there can be a coronal structure above the boundary layer which act as soft X-ray emitter for SW UMa.     

Modeling the luminosity function of galactic low-mass X-ray binaries

The evolution of the family of binaries with a low-mass star and a compact neutron star companion (low-mass X-ray binaries (LMXBs) with neutron stars) ismodeled by the method of population synthesis. Continuous Roche-lobe filling by the optical star in LMXBs is assumed to be maintained by the removal of orbital angular momentum from the binary by a magnetic stellar wind from the optical star and the radiation of gravitational waves by the binary. The developed model of LMXB evolution has the following significant distinctions: (1) allowance for the effect of the rotational evolution of a magnetized compact remnant on themass transfer scenario in the binary, (2) amore accurate allowance for the response of the donor star to mass loss at the Roche-lobe filling stage. The results of theoretical calculations are shown to be in good agreement with the observed orbital period-X-ray luminosity diagrams for persistent Galactic LMXBs and their X-ray luminosity function. This suggests that the main elements of binary evolution, on the whole, are correctly reflected in the developed code. It is shown that most of the Galactic bulge LMXBs at luminosities L _x > 10³⁷ erg s^?1 should have a post-main-sequence Roche-lobe-filling secondary component (low-mass giants). Almost all of the models considered predict a deficit of LMXBs at X-ray luminosities near ～10^36.5 erg s^?1 due to the transition of the binary from the regime of angular momentum removal by a magnetic stellar wind to the regime of gravitational waves (analogous to the widely known period gap in cataclysmic variables, accreting white dwarfs). At low luminosities, the shape of the model luminosity function for LMXBs is affected significantly by their transient behavior-the accretion rate onto the compact companion is not always equal to the mass transfer rate due to instabilities in the accretion disk around the compact object. The best agreement with observed binaries is achieved in the models suggesting that heavy neutron stars with masses 1.4–1.9M _⊙ can be born.     

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.     

Single X-ray Bursts and the Model of a Spreading Layer of Accreting Matter over the Neutron Star Surface

The excess of the rate of type I X-ray bursts over that expected when the matter fallen between bursts completely burns out in a thermonuclear explosion which is observed in bursters with a high persistent luminosity (4 × 10³⁶ ? L_X ? 2 × 10³⁷ erg s^?1) is explained in terms of the model of a spreading layer of matter coming from the accretion disk over the neutron star surface. In this model the accreting matter settles to the stellar surface mainly in two high-latitude ring zones. Despite the subsequent spreading of matter over the entire star, its surface density in these zones turns out to be higher than the average one by 2–3 orders of magnitude, which determines the predominant ignition probability. The multiple events whereby the flame after the thermonuclear explosion in one ring zone (initial burst) propagates through less densematter to another zone and initiates a second explosion in it (recurrent burst) make a certain contribution to the observed excess of the burst rate. However, the localized explosions of matter in these zones, after which the burning in the zone rapidly dies out without affecting other zones, make a noticeably larger contribution to the excess of the burst rate over the expected one.     

An electrically powered binary star?

We propose a model for stellar binary systems consisting of a magnetic and a non-magnetic white dwarf pair which is powered principally by electrical energy. In our model the luminosity is caused by resistive heating of the stellar atmospheres arising from induced currents driven within the binary. This process is reminiscent of the Jupiter–Io system, but greatly increased in power because of the larger companion and stronger magnetic field of the primary. Electrical power is an alternative stellar luminosity source, following on from nuclear fusion and accretion. We find that this source of heating is sufficient to account for the observed X-ray luminosity of the 9.5-min binary RX J1914+24, and provides an explanation for its puzzling characteristics.     

X-ray confirmation of the new intermediate polar RX J1238–38

By finding a 2147-s X-ray pulsation in the recently identified ROSAT   source RX  1238–38 we confirm that it is a member of the intermediate polar class of cataclysmic variable. We analyse the spectral changes over the white dwarf spin cycle, but are unable to distinguish between competing mechanisms for the cause of the pulsation. RX  1238–38 has an anomalous ratio of spin period to orbital period, similar to that of EX Hya.     

Models for the formation of binary and millisecond radio pulsars

The peculiar combination of a relatively short pulse period and a relatively weak surface dipole magnetic field strength of binary radio pulsars finds a consistent explanation in terms of (i) decay of the surface dipole component of neutron-star magnetic fields on a timescale of (2–5) × 10⁶ yr, in combination with (ii) spin-up of the rotation of the neutron star during a subsequent mass-transfer phase. The four known binary radio pulsars appear to fall into two different categories. Two of them, PSR 0655 + 64 and PSR 1913 + 16, have short orbital periods (<25 h) and high mass functions, indicating companion masses 0.7M⊙ (∼1 (± 0.3) M⊙ and 1.4 M⊙, respectively). The other two, PSR 0820 + 02 and PSR 1953 + 29, have long orbital periods (117d), nearly circular orbits, and low, almost identical mass functions of about 3×10^-3 M⊙, suggesting companion masses of about 0.3M⊙. It is pointed out that these two classes of systems are expected to be formed by the later evolution of binaries consisting of a neutron star and a normal companion star, in which the companion was (considerably) more massive than the neutron star, or less massive than the neutron star, respectively. In the first case the companion of the neutron star in the final system will be a massive white dwarf, in a circular orbit, or a neutron star in an eccentric orbit. In the second case the final companion to the neutron star will be a low-mass (∼ 0.3 M⊙) helium white dwarf in a wide and nearly circular orbit. In systems of the second type the neutron star was most probably formed by the accretion-induced collapse of a white dwarf. This explains in a natural way why PSR 1953 + 29 has a millisecond rotation period and PSR 0820 + 02 has not. Among the binary models proposed for the formation of the 1.5-millisecond pulsar, the only ones that appear to be viable are those in which the companion disappeared by coalescence with the neutron star. In such models the companion may have been a red dwarf of mass 0.03M⊙, a neutron star, or a massive (>0.7M⊙) white dwarf. Only in the last-mentioned case is a position of the pulsar close to the galactic plane a natural consequence. In the first-mentioned case the progenitor system most probably was a cataclysmic-variable binary in which the white dwarf collapsed by accretion.