Doppler tomography of Algols

The technique of Doppler tomography has been influential in the study of mass transfer in Algol‐type interacting binaries. The Algols contain a hot blue dwarf star with a magnetically‐active late‐type companion. In the close Algols, the gas stream flows directly into the photosphere of the blue mass‐gaining star because it does not have enough room to avoid impact with that star. Doppler tomograms of the Algols have been produced from over 2500 time‐resolved spectra at wavelengths corresponding to Hα, Hβ, He I (6678 Å), Si II (6371 Å) and Si IV (1394 ° A). These tomograms display images of accretion structures that include a gas stream, accretion annulus, accretion disk, stream‐star impact region, and occasionally a source of chromospheric emission associated with the cool, mass‐losing companion. Some Algol systems alternate between streamlike and disk‐like states, and provide direct evidence of active mass transfer within the Algols. This work produced the very first images of the gas stream for the entire class of interacting binaries, and demonstrated that the Algols are far more active than formerly believed, with variability on time scales of weeks to months. (© 2004 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Substantial stream–disc overflow found in three-dimensional SPH simulations of cataclysmic variables

We study numerically the interaction of the infalling gas stream and the rim of the accretion disc in cataclysmic variables. The simulations were performed with a smoothed particle hydrodynamics scheme with high spatial resolution. Parameters of the systems AM CVn, OY Car, DQ Her, U Gem and IP Peg were used for the simulations. The simulations cover a wide range of orbital periods, mass ratios and mass transfer rates, as well as different thermal states of the accretion disc. The main result of this study is that the accretion stream is not stopped at the impact region (the bright spot at the outer rim of the disc). In fact, after undergoing the shock interaction, most of the matter is deflected vertically and flows in a more or less diffuse stream to inner parts of the disc, hitting the disc surface close to the circularization radius at orbital phase 0.5. This is a common feature in all systems for all simulated parameters. This stream overflow can cause the X-ray absorption dips observed in cataclysmic variables (CVs) and low-mass X-ray binaries (LMXBs) around orbital phase 0.7, if the inclination is at least 65°. Under certain circumstances, namely a sudden increase of the mass transfer rate from the secondary or a rather small disc, parts of the overflowing stream bounce off the disc surface after hitting it at orbital phase ≈0.5. Another absorption region can be expected around orbital phase 0.2.
In our simulations most of the infalling matter reaches the inner disc very quickly. This must alter the evolution of the quiescent disc and the outburst behaviour considerably compared with purely viscous transport of the material through the disc from the outer rim, and therefore should be taken into account in dwarf nova outburst cycle calculations. To our knowledge, the consequences of such a massive stream overflow for the dwarf nova outburst cycle have not been considered yet.     

Fossil magnetic field of accretion disks of young stars

We elaborate the model of accretion disks of young stars with the fossil large-scale magnetic field in the frame of Shakura and Sunyaev approximation. Equations of the MHD model include Shakura and Sunyaev equations, induction equation and equations of ionization balance. Magnetic field is determined taking into account ohmic diffusion, magnetic ambipolar diffusion and buoyancy. Ionization fraction is calculated considering ionization by cosmic rays and X-rays, thermal ionization, radiative recombinations and recombinations on the dust grains. Analytical solution and numerical investigations show that the magnetic field is coupled to the gas in the case of radiative recombinations. Magnetic field is quasi-azimuthal close to accretion disk inner boundary and quasi-radial in the outer regions. Magnetic field is quasi-poloidal in the dusty “dead” zones with low ionization degree, where ohmic diffusion is efficient. Magnetic ambipolar diffusion reduces vertical magnetic field in 10 times comparing to the frozen-in field in this region. Magnetic field is quasi-azimuthal close to the outer boundary of accretion disks for standard ionization rates and dust grain size a _d=0.1 μm. In the case of large dust grains (a _d>0.1 μm) or enhanced ionization rates, the magnetic field is quasi-radial in the outer regions. It is shown that the inner boundary of dusty “dead” zone is placed at r=(0.1–0.6) AU for accretion disks of stars with M=(0.5–2)?M _⊙. Outer boundary of “dead” zone is placed at r=(3–21) AU and it is determined by magnetic ambipolar diffusion. Mass of solid material in the “dead” zone is more than 3?M _⊕ for stars with M≥1?M _⊙.     

The Sun and the transcendental world of binaries

The theory of gravity says that a binary with orbital frequency ν should be a source of gravitational waves at the double frequency and higher harmonics. This implies that long-term exposure of an ensemble of binaries to gravity waves with frequency ν _G can result in (a) a lack of binaries with frequencies near frequency ν _G /2 and its higher harmonics (the effect of unstable orbits) and/or (b) an excess of binaries whose orbital frequencies are “absolutely” incommensurable with ν _G /2 and its higher harmonics (the effect of stable orbits). It is assumed that the stable-orbit frequencies are almost equal to multiples of πν _G /2 and ν _G /2π, where π plays the role of a “perfect” factor ensuring the best antiresonance of binaries. The statistical analysis of frequencies of 5774 Galactic close binary systems (CBSs) with periods P less than 10 days is based on calculating the resonance spectrum that indicates the best common multiple for a given set of frequencies with allowance for the factor π. The CBS distribution turns out to be modulated by the frequency ν _G = 104.4(5) μHz, and this effect is the most pronounced for superfast and compact rotators, such as cataclysmic variables (CVs) and related objects. The frequency ν _G is, within the error, equal to the “enigmatic” frequency ν₀ = 104.160(1) μHz com discovered earlier in the power spectra of the Sun and brightness variations of some extragalactic sources. This confirms the existence of a “coherent cosmic oscillation” of the Universe with frequency ν₀(ν _G ). The new astrophysical phenomenon naturally explains an “CV period gap” at frequencies ≈πν _G /3 (P ≈ 153 min) and maxima at the neighboring frequencies ≈πν _G /2 and ≈πν _G /4 (P ≈ 102 and ≈204 min, respectively). The remarkable and “mysterious” role of the transcendental number π for the world of binaries is emphasized, and the mystery of physical nature of the “universal” oscillation ν₀(ν _G ) is highlighted.     

Type-Ia supernovae in semidetached binaries

We consider the evolutionary scenarios for close binaries that lead to the formation of semidetached systems in which a white dwarf can accumulate the Chandrasekhar mass through mass accretion from its companion, a main sequence star or a subgiant of mass M ～ 2M_⊙. Such dwarfs probably explode as type-Ia supernovae or collapse to form a neutron star. The population synthesis method is used to analyze the dependence of the model rate of these events in the Galaxy on the common envelope parameter, the mass transfer rate, and the response of a main-sequence star to helium accretion at an intermediate evolutionary stage. The rate of explosions in semidetached systems of this type in the Galaxy was found to be no higher than ?0.2×10^?3 yr^?1, which is less than 10% of the lower level for the empirically estimated SNe Ia rate.     

An analysis of active close binaries (CB) based on photometric measurements

For the purpose of analysing light curves of active CB with an accretion disc being at the evolutional phase of an intensive matter exchange between the components a model for light-curve synthesis has been realized where the attention is given to systems like W Ser not sufficiently examined yet with regard that in them the accretion disc is formed around an ordinary star. In the paper one uses the elements presented in the first paper of this series. The model can successfully describe the essential characteristics of the observed light curves due to existence of an accretion disc and a hot spot, as well as those originated in the temperature distribution along the disc radius. The system components are considered in the framework of the nonsynchronous Roche model and the accretion disc of a constant thickness lies in the orbital plane around the star capturing the matter of the neighbouring component.The primary surrounded by the disc is situated relatively well within the Roche oval and its rotation can be significantly nonsynchronous. Near the Lagrange equilibrium pointL ₁ flows from the secondary (which fills the Roche limit) the gas stream nourishing the disc. In the zone where the stream touches the lateral side of the disc a hot spot is formed.The proposed model enables estimating of the basic orbital and physical parameters of active type W Ser CB (Djuraevi, 1991) on the basis of photometric measurements by applying the inverse-problem method.     

慧眼-HXMT观测到的黑洞双星系统中的W-K 关系

通过几十年的观测研究, 黑洞X射线双星(X-Ray Binary, XRB)部分特征被揭示. 然而, 吸积盘结构尚不确定. 黑洞XRB功率密度谱的截断频率与准周期振荡(Quasi Periodic Oscillation, QPO)的相关性质(W-K关系)可以限制吸积盘结构. 利用慧眼-HXMT (Hard X-ray Modulation Telescope)观测到的5个黑洞XRB的数据, 对黑洞XRB的W-K关系进行了研究, 结果表明在慧眼-HXMT观测的3个探测器能段中W-K关系成立. 此外在MAXI J1535-571之中存在截断频率和吸积盘内半径的相关关系, 这和截断的吸积盘结构一致. 如果观测到的功率密度谱来自质量吸积率的扰动传播, 可以推测吸积盘内半径接近最内圆形稳定轨道, 此黑洞可能是高自旋系统.     

Small-scale inviscid accretion discs around black holes

Gas falling quasi-spherically on to a black hole forms an inner accretion disc if its specific angular momentum l exceeds l ∗∼ r _g c , where r _g is the Schwarzschild radius. The standard disc model assumes l ≫ l ∗. We argue that, in many black hole sources, accretion flows have angular momenta just above the threshold for disc formation, l ≳ l ∗, and assess the accretion mechanism in this regime. In a range l ∗< l < l _cr, a small-scale disc forms in which gas spirals fast into the black hole without any help from horizontal viscous stresses. Such an 'inviscid' disc, however, interacts inelastically with the feeding infall. The disc–infall interaction determines the dynamics and luminosity of the accretion flow. The inviscid disc radius can be as large as 14 r _g, and the energy release peaks at 2 r _g. The disc emits a Comptonized X-ray spectrum with a break at ∼100 keV. This accretion regime is likely to take place in wind-fed X-ray binaries and is also possible in active galactic nuclei.     

The outburst properties of AM CVn stars

We briefly summarise the observational properties of ultra‐compact binaries called AM CVn stars. We analyse their outbursts originating from the thermal‐viscous instability in helium accretion disc. We present our preliminary results in applying the model of Dwarf Novae outbursts to helium discs. We can calculate models of outbursts of reasonable amplitude of 2 mag with a constant α parameter throughout the calculation. Setting the mass transfer rate close to its upper critical value produces model lightcurves that resemble short superoutbursts (© 2010 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

A new model for QPOs in accreting black holes: application to the microquasar GRS 1915+105

In this paper we extend the idea suggested previously by Pétri (Astron. Astrophys. 439:L27, 2005a; 443:777, 2005b) (papers I and II) that the high frequency quasi-periodic oscillations (HF-QPOs) observed in low-mass X-ray binaries (LMXBs) may be explained as a resonant oscillation of the accretion disk with a rotating asymmetric background (gravitational or magnetic) field imposed by the compact object. Here, we apply this general idea to black hole binaries. It is assumed that a test particle experiences a similar parametric resonance mechanism such as the one described in paper I and II but now the resonance is induced by the interaction between a spiral density wave in the accretion disk, excited close to the innermost stable circular orbit, and vertical epicyclic oscillations. We use the Kerr spacetime geometry to deduce the characteristic frequencies of this test particle. The response of the test particle is maximal when the frequency ratio of the two strongest resonances is equal to 3:2 as observed in black hole candidates. Finally, applying our model to the microquasar GRS 1915+105, we reproduce the correct value of several HF-QPOs. Indeed the presence of the 168/113/56/42/28 Hz features in the power spectrum time analysis is predicted. Moreover, based only on the two HF-QPO frequencies, our model is able to constrain the mass M _BH and angular momentum a _BH of the accreting black hole. We show the relation between M _BH and a _BH for several black hole binaries. For instance, assuming a black hole weakly or mildly rotating, i.e. a _BH≤0.5?G? M _BH/c ², we find that for GRS 1915+105 its mass satisfies 13?M _⊙≤M _BH≤20?M _⊙. The same model applied to two other well-known BHCs gives for GRO J1655-40 a mass 5?M _⊙≤M _BH≤7?M _⊙ and for XTE J1550-564 a mass 8?M _⊙≤M _BH≤11?M _⊙. This is consistent with other independent estimations of the black hole mass. Finally for H1743-322, we found the following bounds, 9?M _⊙≤M _BH≤13?M _⊙.     

A stochastic Monte Carlo approach to model real star cluster evolution – II. Self-consistent models and primordial binaries

The new approach outlined in Paper I to follow the individual formation and evolution of binaries in an evolving, equal point-mass star cluster is extended for the self-consistent treatment of relaxation and close three- and four-body encounters for many binaries (typically a few per cent of the initial number of stars in the cluster mass). The distribution of single stars is treated as a conducting gas sphere with a standard anisotropic gaseous model. A Monte Carlo technique is used to model the motion of binaries, their formation and subsequent hardening by close encounters, and their relaxation (dynamical friction) with single stars and other binaries. The results are a further approach towards a realistic model of globular clusters with primordial binaries without using special hardware. We present, as our main result, the self-consistent evolution of a cluster consisting of 300 000 equal point-mass stars, plus 30 000 equal-mass binaries over several hundred half-mass relaxation times, well into the phase where most of the binaries have been dissolved and evacuated from the core. The cluster evolution is about three times slower than found by Gao et al. Other features are rather comparable. At every moment we are able to show the individual distribution of binaries in the cluster.     

Accretion in Strong Gravity: from Galactic to Supermassive Black Holes

The galactic black hole binary systems give an observational template showing how the accretion flow changes as a function of increasing mass accretion rate, or L/L_Edd. These data can be synthesised with theoretical models of the accretion flow to give a coherent picture of accretion in strong gravity, in which the major hard-soft spectral transition is triggered by a change in the nature and geometry of the inner accretion flow from a hot, optically thin plasma to a cool, optically thick accretion disc. However, a straightforward application of these models to AGN gives clear discrepancies in overall spectral shape. Either the underlying accretion model is wrong, despite its success in describing the Galactic systems and/or there is additional physics which breaks the simple scaling from stellar to supermassive black holes.     

X-ray binaries

Summary The various types and classes of X-ray binary are reviewed high-lighting recent results. The high mass X-ray binaries (HMXRBs) can be used to probe the nature of the mass loss from the OB star in these systems. Absorption measurements through one orbital cycle of the supergiant system X1700-37 are well modelled by a radiation driven wind and also require a gas stream trailing behind the X-ray source. In Cen X-3 the gas stream is accreted by the X-ray source via an accretion disk. Changes in the gas stream can cause the disk to thicken and the disk to obscure the X-ray source. How close the supergiant is to corotation seems to be as much a critical factor in these systems as how close it is to filling its Roche lobe. In the Be star X-ray binaries a strong correlation between the neutron stars rotation period and its orbital period has been explained as due to the neutron star being immersed in a dense, slow moving equatorial wind from the Be star. For the X-ray pulsars in the transient Be X-ray binaries a centrifugal barrier to accretion is important in determining the X-ray lightcurve and the spin evolution. The X-ray orbital modulations from the low mass X-ray binaries, LMXRBs, include eclipses by the companion and/or periodic dipping behaviour from structure at the edge of the disk. The corresponding optical modulations show a smooth sinusoidal like component and in some cases a sharp eclipse by the companion. The orbital period of the LMXRB XB1916-05 is 1% longer in the optical compared to that given by the X-ray dip period. The optical period has been interpreted as the orbital period, but this seems inconsistent with the well established view of the origin of the X-ray modulations in LMXRB. A new model is presented that assumes the X-ray dip period is the true orbital period. The 5.2 h eclipsing LMXRB XB2129+47 recently entered a low state and optical observations unexpectedly reveal an F star which is too big to fit into the binary. This is probably the first direct evidence that an X-ray binary is part of a hierarchical triple. Finally the class of X-ray binaries containing black hole candidates is reviewed focusing on the value of using X-ray signatures to identify new candidates.     

Supersonic steady accretion of a rotating cold “Iron Gas” onto a neutron star after gravitational collapse

We analytically generalize the well-known solution of steady supersonic spherically symmetric gas accretion onto a star (Bondi 1952) for an iron atmosphere with completely degenerate electrons with an arbitrary degree of relativity. This solution is used for typical physical conditions in the vicinity of protoneutron stars produced by gravitational collapse with masses M ₀=(1.4?1.8)M _⊙ and over a wide range of nonzero “iron gas” densities at infinity, ρ_∞=(10⁴?5×10⁶)g cm^?3. Under these conditions, we determine all accretion parameters, including the accretion rate, whose value is ～(10?50)M _⊙s^?1 at M ₀=1.8M _⊙ (it is a factor of 1.7 lower for M ₀=1.4M _⊙, because the accretion rate is exactly ∝M ₀ ² ). We take into account the effect of accreting-gas rotation in a quasi-one-dimensional approximation, which has generally proved to be marginal with respect to the accretion rate.     

A photometric and spectroscopic study of WW And – an Algol-type,long period binary system with an accretion disc

We have analysed the available spectra of WW And and for the first time obtained a reasonably well defined radial velocity curve of the primary star. Combined with the available radial velocity curve of the secondary component, these data led to the first determination of the spectroscopic mass ratio of the system at q_spec = 0.16 ± 0.03. We also determined the radius of the accretion disc from analysis of the double-peaked H_α emission lines. Our new, high-precision, Johnson VRI and the previously-available Strömgren vby light curves were modelled with stellar and accretion disc models. A consistent model for WW And – a semidetached system harbouring an accretion disc which is optically thick in its inner region, but optically thin in the outer parts – agrees well with both spectroscopic and photometric data.     

An analysis of the X and UV emissions of active galactic nuclei

Under the assumption of spherical accretion, we used the synchrotron radiation mechanism to derive the relation between the 2 keV X luminosity (1_X) and the 2500 A UV luminosity (1_op) of active galactic nuclei. We found that, when the accretion rate is low, 1_X is directly proportional to 1_op. When the accretion rate is high, because the X photons come from a dense region close to the central black hole, they are repeatedly Compton scattered by the thermal electrons and “dragged” into the black hole. In this case, by calculating the relation between the dragging region and the accretion rate, we found 1_X proportional to 1_opp to the power 0.30, thus explaining the tendency towards “saturation” in the observed 1_X − 1_op curve.     

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.     

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.     

Accretion disc dynamics in extreme mass ratio compact binaries

An analysis is presented of a numerical investigation of the dynamics and geometry of accretion discs in binary systems with mass ratios   q = M ₂/ M ₁ < 0.1  , applicable to ultracompact X-ray binaries, AM CVn stars and very short period cataclysmic variables. The steady-state geometry of the disc in the binary reference frame is found to be quite different from that expected at higher mass ratios. For   q ∼ 0.1  , the disc takes on the usual elliptical shape, with the major axis aligned perpendicular to the line of centres of the two stars. However, at smaller mass ratios the elliptical gaseous orbits in the outer regions of the disc are rotated in the binary plane. The angle of rotation increases with gas temperature, but is found to vary inversely with q . At   q = 0.01  , the major axis of these orbits is aligned almost parallel to the line of centres of the two stars. These effects may be responsible for the similar disc structure inferred from Doppler tomography of the AM CVn star GP Com, which has   q = 0.02  . The steady-state geometry at low mass ratios is not predicted by an inviscid, restricted three-body model of gaseous orbits; it is related to the effects of tidal-viscous truncation of the disc near the Roche lobe boundary. Since the disc geometry can be inferred observationally for some systems, it is proposed that this may offer a useful diagnostic for the determination of mass ratios in ultracompact binaries.