Structure of gas flows in Z And in quiescence and during outbursts

A new mechanism is proposed to account for transitions between the quiescent and active states of symbiotic stars. A numerical study of the gas dynamics of the flows in the symbiotic star Z And shows that even small variations in the velocity of the wind from the cool giant can abruptly change the flow structure near the hot component. Such changes alter the accretion regime as the wind velocity increases: disk accretion makes a transition to accretion from the flow. Our calculations indicate that the accretion rate increases by a factor of several tens over a short time interval (～0.1 of the orbital period) during the rearrangement of the flow, when the accretion disk is destroyed.     

Doppler mapping of the SS Cyg system during outburst

Spectral observations of the SS Cyg system in its active state are used to construct Hβ and Hγ Doppler tomograms. These are compared with analogous tomograms for the quiescent state and synthetic tomograms derived from the results of three-dimensional gas-dynamical modeling. The parameters of the accretion disk during the outburst are estimated. An explanation for the observed flow pattern is proposed, based on a numerical model with an elliptical accretion disk.     

Formation of a radiative wind and accretion disk in microquasars. Generation of flare activity in the disk due to an increase in the supply of material to the Lagrange point L<Subscript>1</Subscript>. The system LMC X-3

Three-dimensional numerical hydrodynamical modeling of a radiative wind and accretion disk in a close binary system with a compact object is carried out, using the massive X-ray binary LMC X-3 as an example. This system contains a precessing disk, and may have relativistic jets. These computations show that an accretion disk with a radius of about 0.20 (in units of the component separation) forms from the radiative wind from the donor when the action of the wind on the central source is taken into account, when the accretion rate is equal to the observed value (about 3.0 × 10^?8 M _⊙/year, which corresponds to the case when the donor overflows its Roche lobe by nearly 1%). It is assumed that the speed of the donor wind at infinity is about 2200 km/s. The disk that forms is geometrically thick and nearly cylindrical in shape, with a low-density tunnel at its center extending from the accretor through the disk along the rotational axis. We have also modeled a flare in the disk due to short-term variations in the supply of material through the Lagrange point L₁, whose brightnesses and durations are able to explain flares in cataclysmic variables and X-ray binaries. The accretion disk is not formed when the donor underfills its Roche lobe by 0.5%, which corresponds to an accretion rate onto the compact object of 2.0 × 10^?9 M _⊙/year. In place of a disk, an accretion envelope with a radius of about 0.03 forms, within which gas moves along very steep spiral trajectories before falling onto the compact object. As in the accretion-disk case, a tunnel forms along the rotational axis of the accretion envelope; a shock forms behind the accretor, where flares occur in a compact region a small distance from the accretor at a rate of about six flares per orbital period, with amplitudes of about 10 ^m or more. The flare durations are two to four minutes, and the energies of individual particles at the flare maximum are about 100–150 keV. These flares appear to be analogous to the flares observed in gamma-ray and X-ray burst sources. We accordingly propose a model in which these phenomena are associated with massive, close X-ray binary systems with component-mass ratios exceeding unity, in which the donor does not fill its Roche lobe. Although no accretion disk forms around the compact object, an accretion region develops near the accretor, where the gamma-ray and X-ray flares occur.     

Hydrodynamical modeling of mass transfer in the close binary system <Emphasis Type="Italic">β</Emphasis> Lyr

We present results of two-dimensional hydrodynamical simulations of mass transfer in the close binary system β Lyr for various radii of the accreting star and coefficients describing the interaction of the gaseous flow and the main component (primary). We take the stellar wind of the donor star into account and consider various assumptions about the radiative cooling of the gaseous flow. Our calculations show that the initial radius of the flow corresponding to our adopted mass-transfer rate through the inner Lagrange point (L₁) of (1–4) × 10^?5M_⊙/yr is large: 0.22–0.29 (in units of the orbital separation). In all the models, the secondary loses mass through both the inner and outer (L₁ and L₂) Lagrange points, which makes the mass transfer in the system nonconservative. Calculations for various values of the primary radius show a strong dependence on the coefficient fv that models the flow-primary interaction. When the radius of the primary is 0.5, there is a strong interaction between the gas flow from L1 and the flow reflected from the primary surface. For other values of the primary radius (0.1 and 0.2), the flow does not interact directly with the primary. The flow passes close to the primary and forms an accretion disk whose size is comparable to that of the Roche lobe and a dense circum-binary envelope surrounding both the disk and the binary components. The density in the disk varies from 10¹² to 10¹⁴ cm^?3, and is 10¹⁰–10¹² cm^?3 in the circum-binary envelope. The temperature in the accretion disk ranges from 30000 to 120000 K, while that in the circum-binary envelope is 4000–18000 K. When radiative cooling is taken into account explicitly, the calculations reveal the presence of a spiral shock in the accretion disk. The stellar wind blowing from the secondary strongly interacts with the accretion disk, circum-binary envelope, and flow from L₂. When radiative cooling is taken into account explicitly, this wind disrupts the accretion disk.     

Photometric parameters of the dwarf nova SS Cygni in the quiescent state

The mean 1983–1996 UBV light curves of the dwarf nova SS Cyg are used to derive the binary parameters in the quiescent state. Solutions are obtained for a classical hot-spot model and a model with an energy source lying outside the accretion disk. Photometric and spectroscopic data are combined to infer the masses and radii of the binary components. The white dwarf in SS Cyg is one and a half times as massive as the red dwarf, q=M _wd /M _rd ～1.45, M _rd ～0.46M _⊙ and M _wd ～0.66M _⊙. The orbital inclination of the system is i?51°–54°. The contribution of the accretion disk to the total flux in the quiescent state is estimated to be ～47–49% and ～54% in the VU and B filters, respectively. The hot spot contributes less than ～3% to the total optical flux. In the “non-classical” hot-spot model, the disk and bulge contributions are 27 and 2–8%, respectively, depending on the orbital phase. The shape of the mean light curves of SS Cyg suggests asymmetric heating of the red-dwarf surface in the quiescent state by high-temperature radiation generated in the hot-spot region.     

Three-dimensional hydrodynamical modeling of accretion-disk formation in microquasars. Cen X-3

We have carried out three-dimensional hydrodynamical modeling of the formation of an accretion disk around a compact object due to radiative wind of a massive donor in a close binary system. The massive X-ray binary Cen X-3, which has a precessing accretion disk and may possess relativistic jets, is considered as an example. The computations show that, when the action of the central compact object on the formation of the wind is taken into account, the radiative wind forms an accretion disk with a radius of 0.16 (in units of the orbital separation), which accretes at a rate close to 1 × 10^?8 M _⊙/yr. In this model, the disk is spherically symmetrical and geometrically thick, with a tunnel going from the accretor to the upper layers of the disk along the accretor’s rotational axis at the disk center. The number density of the gas in the tunnel is five orders of magnitude lower than in the disk. The wind-disk interaction at the outer boundary of the disk produces a strong shock (wind-disk shock) directed toward the donor. The black-body emission of the disk and tunnel is nonstationary, and resembles the outbursts observed in Cen X-3. An analysis of the location of the region of nonstationary emission suggests that the outbursts occur in the wind-disk shock.     

Three-dimensional hydrodynamical modeling of mass transfer in the close binary system β Lyr

We present a three-dimensional hydrodynamical modeling of mass transfer in the close binary system β Lyr taking radiative cooling into account explicitly. The assumed mass-transfer rate through the first Lagrangian point L₁ is 3.0 × 10^?5 M _⊙/yr. A flow with a radius of 0.14–0.16 (in units of orbital separation) is formed in the vicinity of L₁. This flow forms an accretion disk with a radius close to 23 R _⊙ and a thickness of about 10 R _⊙. The accretion disk is surrounded by an outer envelope that extends beyond the computational domain. A spiral shock forms at the outer boundary of the disk at orbital phase 0.25. Geometrically, the disk is toruslike, while the outer envelope is cylinder-like. In this model, which has low temperatures inside the computational domain, no jetlike structures form in the disk. It is possible that the jetlike structure in β Lyr arises due to the interaction of radiative wind from the accretor with the flow from L₁. In the model considered, a hot region exists over the poles of the accretor at a height of about 0.21. The amount of matter lost by the system is close to 10% of the mass flowing through L₁; i.e., the mass transfer in the system is almost conservative. For a mass-transfer rate of 3.0 × 10^?5 M _⊙/yr, the orbital period varies by 40.4 s/yr. This means that the observed variation of the orbital period of 19 s/yr should correspond to a mass-transfer rate close to 1.0 × 10^?5 M _⊙/yr.     

Three-dimensional hydrodynamical modeling of mass transfer in the close binary system β Lyr with an accretor wind

We present three-dimensional hydrodynamical modeling of mass transfer in the close binary system β Lyr taking into account explicitly radiative cooling and the stellar wind of the accretor. Our computations show that flow forces wind out from the orbital plane, where an accretion disk with a radius of 0.4–0.5 and a height of about 0.15–0.17 (in units of orbital separation) is formed. Gas motions directed upward from the orbital plane are initiated in the region of interaction of the flow from L₁ and the accretor wind (x = 0.91, y = ?0.17); i.e., a jetlike structure forms. This structure has the shape of a gas pillar above the orbital plane, where gas moves with the velocity of stellar wind. The number density of the gas in this structure is about 10¹⁴ cm^?3, and its temperature is 20 000–45 000 K. At heights of about 0.15–0.20 above the orbital plane, in the region between the jetlike structure and the disk, two spiral shocks form. It is possible that the emission lines observed in the spectrum of β Lyr binary originate in this region.     

Interpretation of light curves of the cataclysmic variable OY Car in a model with shockless interaction between a gaseous stream and the disk

To determine the parameters of the accretion disk and shock-wave region responsible for the formation of the orbital peak in the light curve of the binary system OY Car (an SU UMa-type variable), we have analyzed its U BV R and JK light curves using two gas-dynamical models with different regions of shock interaction: one with a hot line along the stream from the Lagrange point L1 and one with a hot spot on the accretion disk. The hot-line model can better describe the quiescent state of the system: the maximum X² for the optical light curves does not exceed 207, whereas the minimum residual for the hot-spot model is X²>290. The shape of the eclipse is almost identical in both models; the main differences are in interpreting out-of-eclipse portions of the light curves, whose shape can varyin the transition from one orbital cycle to another. The hot-spot model is not able to describe variations of the system’s brightness at orbital phases ?～0.1–0.6. The rather complex behavior of the observed flux in this phase interval can be explained in the hot-line model as being due to variations of the temperature and size of the system. Based on the analysis of a sequence of 20 B curves of OY Car, we conclude that the flux variations in the primary minimum are due to variations of the luminosity of the accretion disk, whereas the flux variability in the vicinity of the orbital peak is due to the combined effect of the radiation of the disk and hot line. The JK light curves of OY Car in the quiescent state and during a small flare also indicate preference for the hot-line model, since the primaryminimum and the flux near quadratures calculated using the hot-spot model are not consistent with the observations.     

Spectroscopy and doppler mapping of the binary SS Cyg during outburst

We carried out spectroscopy of the binary SSCyg in the Hα, Hβ, and Hγ lines in its active state in August and December 2006. We have estimated the parameters of the main flow elements contributing to the spectra. Profile variations during the orbital period are analyzed, and a Doppler tomogram computed for the Hα line. We consider the evolution of the line profiles with the development of the outburst. A phenomenological model explaining the observed outburst features is suggested. In this model, the main elements of the flow determining the shape of the spectral lines are the accretion disk, a toroidal shell formed in the inner parts of the disk, an expanding spherical shell around the accreting star, a region in front of the bow shock that forms due to the orbital motion of the disk in the circumbinary envelope, and the surface of the donor star near the inner Lagrange point, L1, which is heated by radiation from the accretor.     

Changes in the parameters of the accretion disk around the dwarf nova SDSS J090350.73+330036.1 during an outburst

R-band photometric light curves of the eruptive eclipsing binary SDSS J090350.73+330036.1 obtained during a superoutburst in May 2010 (JD 2455341-2455347) are analyzed. Observations covering an interval near the outburst maximum and the post-maximum decrease by 0.7 ^m are presented. Oscillations (superhumps) whose period differs from the orbital period by several percent are observed in the light curve together with eclipses, suggesting that the studied system is a SU UMa dwarf nova. A ??spiral arm?? model is used to fit the light curves and determine the parameters of the accretion disk and other components of the binary system. Together with a hot line, this model takes into account, geometrical inhomogeneities on the surface of the accretion disk, namely, two thickenings at its outer edge that decrease exponentially in the vertical direction with approach toward the white dwarf. The increase in the R-band flux from the system during the superoutburst mainly results from the enhanced luminosity of the accretion disk due to the increase in its radius by up to ??0.44a ₀ at the outburst maximum (a ₀ is the component separation), as well as a shallower radial temperature decrease law than in the canonical case. As the superoutburst faded, the disk radius decreased smoothly at the end of our observation (to ??0.33a ₀), the thickness of its outer edge and temperature of its boundary layer decreased, and the parameter ?? _g approached its canonical value. Deviations from the mean brightness of the system as a function of the superhump period P _sh are analyzed for each out-of-eclipse set of observations. Various factors affecting the appearance and amplitudes of superhumps in the orbital light curves are considered.     

Formation of the accretion disk in the SS 433 system with explicit radiative cooling,convective heat conduction,and radiation pressure

The paper continues three-dimensional hydrodynamical computations of the formation of an accretion disk in the SS 433 system, taking into account radiative cooling explicitly, convective thermal conductivity, and radiation pressure. The computational results show that the powerful, broad flow forms an optically thick accretion disk with a gas density of 10¹²–10¹⁴ cm^?3, a temperature of 15000–35 000 K, a radius of about 0.3, and a height of 0.2–0.3 (in units of the component separation). Spiral shocks form in the disk, and a narrow conelike cavity (tunnel) forms at the center. In this tunnel, gas is accelerated to relativistic speeds, leaving the system in the form of narrow jets.     

Non-steady-state accretion disks in X-ray novae: Outburst models for Nova Monocerotis 1975 and Nova Muscae 1991

We have fit outbursts of two X-ray novae (Nova Monocerotis 1975=A0620-00 and Nova Muscae GS 1991=1124-683) using a non-steady-state accretion-disk model. The model is based on a new solution for a diffusion-type equation for non-steady-state accretion and describes the evolution of a viscous α disk in a binary system after the peak of the outburst, when the matter in the disk is totally ionized. The accretion rate in the disk decreases according to a power law. We derive formulas for the accretion rate and effective temperature of the disk. The model has three free input parameters: the mass of the central object M, the turbulence parameter α, and the normalization parameter δt. The results of the modeling are compared with the observed X-ray and optical B and V light curves. The estimates for the turbulence parameter α are similar: 0.2–0.4 for A 0620-00 and 0.45–0.65 for GS 1124-683, suggesting a similar nature for the viscosity in the accretion disks around the compact objects in these sources. We have also derived the distances to these systems as functions of the masses of their compact objects.     

Morphology of the interaction between the stream and cool accretion disk in a semidetached binary system

We analyze heating and cooling processes in accretion disks in binaries. For realistic parameters of the accretion disks in close binaries (\(\dot M \simeq 10^{ - 12} - 10^7 M_ \odot /yr\) and α?10^?1–10^?2), the gas temperature in the outer parts of the disk is from ～10⁴ to ～10⁶ K. Our previous gas-dynamical studies of mass transfer in close binaries indicate that, for hot disks (with temperatures for the outer parts of the disk of several hundred thousand K), the interaction between the stream from the inner Lagrange point and the disk is shockless. To study the morphology of the interaction between the stream and a cool accretion disk, we carried out three-dimensional modeling of the flow structure in a binary for the case when the gas temperature in the outer parts of the forming disk does not exceed 13 600 K. The flow pattern indicates that the interaction is again shockless. The computations provide evidence that, as is the case for hot disks, the zone of enhanced energy release (the “hot line”) is located beyond the disk and originates due to the interaction between the circumdisk halo and the stream.     

Structure of the accretion disk in the cataclysmic variable UX UMa in the transition to its quiescent state

Our analysis of BV RI light curves for the cataclysmic variable UX UMa obtained in intermediate activity states, in the transition between the active and quiescent states of the system on March 12, 1997 and May 3, 2000, shows that the shapes of these light curves cannot be adequately described using the standard hot-spot model. A model with a “hot line” near the edge of the disk and a two-armed spiral structure on the disk surface reproduces much better out-of-eclipse variations in the light curves. The presence of an extended hot line can explain the anomalous shape of the I light curve on March 12, 1997. The decrease in the observed luminosity of the system between March 12, 1997 and May 3, 2000 could be due to a decrease in the disk luminosity by a factor of 2–2.5; the higher disk luminosity on the earlier date is associated with appreciable deviations of the radial temperature distribution of the disk material from that for the standard model. The phases and depths of dips in the out-of-eclipse sections of the UX UMa light curves are due primarily to the parameters of the complex shape of the accretion disk, which has a spiral structure located mainly near its outer edge. The contribution of the spiral arms in the V filter reaches 20–50% of the total disk radiation. The crest of the first spiral wave in our model maintains its approximate position in azimuth; this structure could represent a bulge in a halo at the outer edge of the disk near orbital phases φ ～ 0.7, in the direction of the continuation of the extended shock in the disk itself. The position of the crest of the second spiral arm changes with time. This structure may represent a one-armed spiral wave near the apastron of the weakly elliptical disk. Finally, the observations testify to the presence of another spiral arm that is les clearly manifest in terms of both its luminosity and its height above the unperturbed disk surface. Thus, in an intermediate activity state of UX UMa, the surface of the accretion disk is distorted by the action of a two-armed spiral structure in the outer regions of the disk, which is asymmetric in both its luminosity and dimensions, and a bulge at the disk edge in the region of its interaction with the inflow to the disk.     

The Hydrogen and Helium Lines of the Symbiotic Binary Z And during Its Brightening at the End of 2002

High resolution observations in the region of the Hα, HeII λ 4686, and Hγ lines in the spectrum of the symbiotic binary Z And were performed during a small-amplitude flare at the end of 2002. The profiles of the hydrogen lines were double-peaked, and suggest that the lines may be emitted mainly by an optically thin accretion disk. Since the Hα line is strongly contaminated by emission from the envelope, the Hγ line is used to investigate the properties of the accretion disk. The Hα line has broad wings, believed to be determined mostly by radiation damping, although the high-velocity stellar wind from the compact object in the system may also contribute. The Hγ line has a broad emission component, assumed to be emitted mainly from the inner part of the accretion disk. The HeIIλ 4686 line also has a broad emission component, but is believed to arise in a region of high-velocity stellar wind. The outer radius of the accretion disk can be calculated from the shift between the peaks. Assuming that the orbital inclination can range from 47° to 76°, we estimate the outer radius to be 20–50 R_⨀. The behavior of the observed lines can be interpreted in the model proposed for the line spectrum during the first large 2000–2002 flare of this binary.     

The matter-flow structure in the SS Cyg system in its quiescent state from comparisons of observational and synthetic Doppler tomograms

Doppler tomograms are constructed for the quiescent state of the SS Cyg system based on Hβ and Hγ spectral-line observations carried out in August 2006 with the 2-m telescope at Terskol Peak. Gasdynamical simulations combined with the Doppler tomograms enable identification of the main features of the flow. Comparisons of synthetic tomograms with observations indicate that an accretion disk is present in the quiescent system. In the tomograms, the luminosity is maximum at the arms of the spiral tidal shock at the shock front due to the interaction between the gas of the circum-binary envelope and material in the stream issuing from the Lagrangian point L₁ (the “hot line”), and in the region behind the bow shock due to the motion of the accretor and disk in the gas of the envelope. The contribution of this last element results in appreciable asymmetry of the tomograms.     

Ring structures in Orion KL

We present the results of studies of the superfine structure of H₂O maser sources in the Orion Nebula. Powerful, low-velocity, compact maser sources are distributed in eight active zones. Highly organized structures in the form of chains of compact components were revealed in two of these, in the molecular cloud OMC-1. The component sizes are ～0.1 AU and their brightness temperatures are T _b =10¹²?10¹⁶ K. The structures correspond to tangential sections of concentric rings viewed edge-on. The ring emission is concentrated in the azimuthal plane, decreasing the probability of their discovery. The formation of protostars is accompanied by the development of accretion disks and bipolar flows, with associated H₂O maser emission. The accretion disks are in the stage of fragmentation into protoplanetary rings. In a Keplerian approximation, the protostars have low masses, possibly evidence for instability of the systems. Supermaser emission of the rings is probably triggered by precession of the accretion disk. The molecular cloud’s radial velocity is V _LSR=7.74 km/s and its optical depth is τ≈5. The emission from components with velocities within the maser window is additionally amplified. The components’ emission is linearly polarized via anisotropic pumping.     

Period clustering of anomalous X-ray pulsars

The question of why the observed periods of anomalous X-ray pulsars (AXPs) and soft gamma-ray repeaters (SGRs) cluster in the range 2–12 s is discussed. The possibility that AXPs and SGRs are the descendants of high-mass X-ray binaries that have disintegrated in core-collapse supernova explosions is investigated. The spin periods of neutron stars in high-mass X-ray binaries evolve towards the equilibrium period, which is a few seconds, on average. After the explosion of its massive companion, the neutron star becomes embedded in a dense gaseous envelope, and accretion from this envelope leads to the formation of a residual magnetically levitating disk. It is shown that the expected mass of the disk in this case is 10^?7–10^?8 M_⊙, which is sufficient to support accretion at the rate 10¹⁴–10¹⁵ g/s over a few thousand years. During this period, the star manifests itself as an isolated X-ray pulsar with a number of parameters similar to those of AXPs and SGRs. The periods of such pulsars can cluster if the lifetime of the residual disk does not exceed the spin-down timescale of the neutron star.     

Light-curve syntheses for close binary systems: Modeling spiral waves in an elliptical disk around a white dwarf

We present an algorithm for synthesizing the light curve of a close binary consisting of a normal star (a red dwarf that fills its Roche lobe) and a spherical star (a white dwarf). The spherical component is surrounded by an elliptical accretion disk with a complex shape: it is geometrically thin near the spherical star and geometrically thick at the edge of the disk. An additional complication is presented by the presence of a one-or two-armed spiral pattern at the inner surface of the disk. The maximum height of the spiral arm above the disk surface is located at ～9 R _d , and the height decreases exponentially as the arm approaches the inner regions of the disk. Shielding of the inner hot parts of the disk by the crests of the spirals results in the formation of “steps” in out-of-eclipse parts of the orbital light curves. The algorithm takes into account the presence of a “hot line” by the lateral surface of the disk, making it possible to model binary systems in both quiescence and outburst. In the latter case, the hot line degenerates into a small bulge at the outer lateral surface of the disk, which can be considered an analog of a hot spot. The algorithm was applied to the orbital light curve of the cataclysmic binary IP Peg during its October 30, 2000, outburst. To explain the variations of the out-of-eclipse brightness of the system during the outburst, it is necessary to include the presence of a one-armed spiral wave at the inner surface of the disk, close to the periastron of the elliptical disk. We have obtained the parameters of IP Peg during the outburst for various models of the system.