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.     

Three-dimensional hydrodynamical modeling of the two-component wind and accretion disk in the close binary β Lyrae

This paper continues a series of studies on three-dimensional hydrodynamical modeling of mass transfer in the binary system β Lyr. The model takes into account the stellar wind from the donor star, which outflows at a rate of , as demonstrated by radio observations. This stellar wind should appreciably influence the formation of the envelope in the binary. Computations have shown that the interaction of the matter flow from the Lagrangian point L₁ and the accretor wind leads to the formation of an optically and geometrically thick gaseous envelope around the accretor. The matter flow meets the accretor wind, spreads out, accumulates over the outer edge of the wind, and forms a geometrically thick envelope (disk). The wind flows freely at the center of the disk, over the accretor poles. Jet-like structures arise beyond the wind-propagation region, above the thick accretion disk. The matter flowing from the outer edge of the disk interacts with the donor wind, leading to the formation of a standing shock between L₁ and the outer edge of the disk, in the direction corresponding to orbital phase 0.25. This shock is able to explain the origin of the X-ray radiation from the binary β Lyr.     

The binary systems IC 10 X-1 and NGC 300 X-1: Accretion of matter from an intense Wolf–Rayet stellar wind onto a black hole

The current evolutionary stage of the binary systems IC 10 X-1 and NGC 300 X-1, which contain a massive black hole and a Wolf–Rayet star with a strong stellar wind that does not fill its Roche lobe, is considered. The high X-ray luminosity and X-ray properties testify to the presence of accretion disks in these systems. The consistency of the conditions for the existence of such a disk and the possibility of reproducing the observed X-ray luminosity in the framework of the Bondi–Hoyle–Littleton theory for a spherically symmetric stellar wind is analyzed. A brief review of information about the mass-loss rates of Wolf–Rayet stars and the speeds of their stellar winds is given. The evolution of these systems at the current stage is computed. Estimates made using the derived parameters show that it is not possible to achieve consistency, since the conditions for the existence of an accretion disk require that the speed of the Wolf–Rayetwind be appreciably lower than is required to reproduce the observedX-ray luminosity. Several explanations of this situation are possible: (1) the real pattern of the motion of the stellar-wind material in the binary is substantially more complex than is assumed in the Bondi–Hoyle–Littleton theory, changing the conditions for the formation of an accretion disk and influencing the accretion rate onto the black hole; (2) some of the accreting material leaves the accretor due to X-ray heating; (3) the accretion efficiency in these systems is nearly an order of magnitude lower than in the case of accretion through a thin disk onto a non-rotating black hole; (4) the intensity of the Wolf–Rayet wind is one to two orders of magnitude lower than has been suggested by modern studies.     

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.     

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.     

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.     

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.     

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 gas-dynamical modeling of changes in the flow structure during the transition from the quiescent to the active state in symbiotic stars

The results of three-dimensional modeling of the flow structure in the classical symbiotic system Z Andromedae are presented. Outbursts in systems of this type occur when the accretion rate exceeds the upper limit of the steady-burning range. Therefore, in order to realize the transition from a quiescent to an active state, it is necessary to find a mechanism capable of sufficiently increasing the accretion rate on the time scales typical for outburst development. Our calculations provide support for a mechanism for the transition from quiescence to outburst in classical symbiotic systems suggested earlier based on two-dimensional calculations (Bisikalo et al., 2002). Our results show that an accretion disk forms in the system for a wind velocity of 20 km s^?1. The accretion rate for the solution with the disk is ～22.5–25% of the mass-loss rate of the donor, which is ～4.5?5 × 10^?8M_⊙ yr^?1 for Z And. This value is in agreement with the steady-burning range for the white-dwarf masses usually accepted for this system. When the wind velocity increases from 20 to 30 km s^?1, the accretion disk is destroyed and the disk material falls onto the accretor surface. This process is followed by an approximately twofold jump in the accretion rate. The resulting growth in the accretion rate is sufficient so as to exceed the upper limit of the steady-burning range, thus bringing the system into an active state. The time during which the accretion rate is above the steady-burning value is in very good agreement with observations. Our analysis leads us to conclude that small variations in the donor wind velocity can lead to the transition from disk accretion to wind accretion and, as a consequence, to the transition from a quiescent to an active state in classical symbiotic stars.     

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.     

Wolf-Rayet stars with relativistic companions

The evolution of close binary systems containing Wolf-Rayet (WR) stars and black holes (BHs) is analyzed numerically. Both the stellar wind from the donor star itself and the induced stellar wind due to irradiation of the donor with hard radiation arising during accretion onto the relativistic component are considered. The mass and angular momentum losses due to the stellar wind are also taken into account at phases when the WR star fills its Roche lobe. It is shown that, if a WR star with a mass higher than ～10M _⊙ fills its Roche lobe in an initial evolutionary phase, the donor star will eventually lose contact with the Roche lobe as the binary loses mass and angular momentum via the stellar wind, suggesting that the semi-detached binary will become detached. The star will remain a bright X-ray source, since the stellar wind that is captured by the black hole ensures a near-Eddington accretion rate. If the initial mass of the helium donor is below ～5M _⊙, the donor may only temporarily detach from its Roche lobe. Induced stellar wind plays a significant role in the evolution of binaries containing helium donors with initial masses of ～2M _⊙. We compute the evolution of three observed WR-BH binaries: Cyg X-3, IC 10 X-1, and NGC 300 X-1, as well as the evolution of the SS 433 binary system, which is a progenitor of such systems, under the assumption that this binary will avoid a common-envelope stage in its further evolution, as it does in its current evolutionary phase.     

Three-dimensional hydrodynamical modeling of the formation of the accretion disk in the SS 433 binary system

Three-dimensional hydrodynamical modeling of the formation of the accretion disk in the SS 433 binary system is carried out with various types of cooling and numerical grids. These computations show that a thick accretion disk with a height of 0.25–0.30 (in units of the component separation) is formed around the compact object, from a flow with a large radius (0.2–0.3 in the same units) that forms in the vicinity of the inner Lagrangian point. This disk has the form of a flattened torus. The number of orbits of a particle of gas in the disk is 100–150, testifying to a minimal influence of numerical viscosity in these computations. The computations also show that the stream flowing from L₁ is nearly conservative, and spirals in the disk are not formed due to the influence of the donor gravitation.     

Features of the Flow Structure in the Vicinity of the Inner Lagrangian Point in Polars

The structures of plasma flows in close binary systems whose accretors have strong intrinsic magnetic fields are studied. A close binary system with the parameters of a typical polar is considered. The results of three-dimensional numerical simulations of the matter flow from the donor into the accretor Roche lobe are presented. Special attention is given to the flow structure in the vicinity of the inner Lagrangian point, where the accretion flow is formed. The interaction of the accretion-flow material from the donor’s envelope with the magnetic field of the accretor results in the formation of a hierarchical structure of the magnetosphere, because less dense areas of the accretion flow are stopped by the magnetic field of the white dwarf earlier than more dense regions. Taking into account this kind of magnetosphere structure can affect analysis results and interpretation of the observations.     

Formation of accretion disks in close-binary systems with magnetic fields

We have developed a three-dimensional numerical model and applied it to simulate plasma flows in semi-detached binary systems whose accretor possesses a strong intrinsic magnetic field. The model is based on the assumption that the plasma dynamics are determined by the slow mean flow, which forms a backdrop for the rapid propagation of MHD waves. The equations describing the slow motion of matter were obtained by averaging over rapidly propagating pulsations. The numerical model includes the diffusion of magnetic field by current dissipation in turbulent vortices, magnetic buoyancy, and wave MHD turbulence. A modified three-dimensional, parallel, numerical code was used to simulate the flow structure in close binary systems with various accretor magnetic fields, from 10⁵ to 10⁸ G. The conditions for the formation of the accretion disk and the criteria distinguishing the two types of flow corresponding to intermediate polars and polars are discussed.     

Magnetic-field structure in the accretion disks of semi-detached binary systems

The results of three-dimensional MHD numerical simulations are used to investigate the characteristic properties of the magnetic-field structures in the accretion disks of semi-detached binary systems. It is assumed that the intrinsic magnetic field of the accretor star is dipolar. Turbulent diffusion of the magnetic field in the disk is taken into account. The SS Cyg system is considered as an example. The results of the numerical simulations show the intense generation of a predominantly toroidal magnetic field in the accretion disk. Magnetic zones with well defined structures for the toroidal magnetic field form in the disk, which are separated by current sheets in which there ismagnetic reconnection and current dissipation. Possible observational manifestations of such structures are discussed. It is shown that the interaction of a spiral precessional wave with the accretor’s magnetosphere could lead to quasi-periodic oscillations of the accretion rate.     

Models of astrophysical decretion-accretion diffusional disks

The formation of gaseous diffusional accretion-decretion disks is an important stage in the evolution of numerous astronomical objects. Matter is accreted onto the object in the accretion part of these disks, while the angular momentum of the accreted matter is transported from the central region to the periphery in the decretion part. Here, we consider general questions connected with the formation and evolution of diffusive accretion-decretion disks in various astrophysical objects. Such disks can be described using nonstationary diffusion models. The phenomenological parameters of these models are the coefficients in the relations for the characteristic turbulent velocity and mean free path of diffusion elements in the disk. We have developed a numerical technique to compute the disk evolution for a number of models (a massive disk, a disk with continuous accretion, a purely decretion disk). Analytical expressions estimating the basic parameters of accretion-decretion disks are presented. We discuss the relationship between the models considered and the classical α model of an accretion disk.     

Two-dimensional hydrodynamical modeling of mass transfer in semidetached binaries with asynchronously rotating components

We have modeled the mass transfer in the three semidetached binaries U Cep, RZ Sct, and V373 Cas taking into account radiative cooling both implicitly and explicitly. The systems have asynchronously rotating components and high mass-transfer rates of the order of 10^?6M_⊙/yr; they are undergoing various stages of their evolution. An accreting star rotates asynchronously if added angular momentum is redistributed over the entire star over a time that exceeds the synchronization time. Calculations have indicated that, in the model considered, mass transfer through the point L₁ is unable to desynchronize the donor star. The formation of an accretion disk and outer envelope depends on the component-mass ratio of the binary. If this ratio is of the order of unity, the flow makes a direct impact with the atmosphere of the accreting star, resulting in the formation of a small accretion disk and a relatively dense outer envelope. This is true of the disks in U Cep and V373 Cas. When the component-mass ratio substantially exceeds unity (the case in RZ Sct), the flow forms a large, dense accretion disk and less dense outer envelope. Taking into account radiative cooling both implicitly and explicitly, we show that a series of shocks forms in the envelopes of these systems.     

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.     

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.     

The influence of the inclination of the accretor’s magnetic axis on the accretion-disk structure in intermediate polars

Typical changes of the accretion-disk structures in intermediate polars are studied as a function of the inclination of the accretor’s magnetic field. Thre-dimensional numerical modeling was performed for seven differentmagnetic-axis inclinations. The results showthat the system forms a magnetosphere region, and that column accretion occurs. The action of the magnetic field tilts the inner parts of the disk along the magnetic axis of the accretor. The magnetic-field inclination appreciably influences matter transfer in the disk and accretion processes. Generation of toroidal magnetic field, magnetic braking, and alignment of the dipole magnetic field result in oscillations of the accretion rate. A direct relationship between the field inclination and the oscillation amplitude is found, as well as an inverse relationship between the field inclination and the oscillation period.