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.     

The evolution of close binary systems with intermediate-mass black holes and ultra-luminous X-ray sources

The results of numerical studies of the evolution of a close binary system containing a black hole with a mass of ～3000M_⊙ are presented. Such a black hole could form in the center of a sufficiently rich and massive globular cluster. The secondary could be a main-sequence star, giant, or degenerate dwarf that fills or nearly fills its Roche lobe. The numerical simulations of the evolution of such a system take into account the magnetic wind of the donor together with the wind induced by X-ray irradiation from the primary, the radiation of gravitational waves by the system, and the nuclear evolution of the donor. Mass transfer between the components is possible when the donor fills its Roche lobe, and also via the black hole’s capture of some material from the induced stellar wind. The computations show that the evolution of systems with solar-mass donors depends only weakly on the mass of the accretor. We conclude that the observed ultra-luminous X-ray sources (L _X ? 10³⁸ erg/s) in nearby galaxies could include accreting black holes with masses of 10²?10⁴M_⊙. Three scenarios for the formation of black holes with such masses in the cores of globular clusters are considered: the collapse of superstars with the corresponding masses, the accretion of gas by a black hole with a stellar initial mass (<100M_⊙), and the tidal accumulation of stellar black holes. We conclude that the tidal accumulation of stellar-mass black holes is the main scenario for the formation of intermediate-mass black holes (10²?10⁴M_⊙) in the cores of globular clusters.     

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.     

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.     

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 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.     

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.     

The mass of the compact object in the low-mass X-ray binary 2S 0921-630

We interpret the observed radial-velocity curve of the optical star in the low-mass X-ray binary 2S 0921-630 using a Roche model, taking into account the X-ray heating of the optical star and screening of X-rays coming from the relativistic object by the accretion disk. Consequences of possible anisotropy of the X-ray radiation are considered. We obtain relations between the masses of the optical and compact (X-ray) components, m _v and m _x , for orbital inclinations i = 60°, 75°, and 90°. Including X-ray heating enabled us to reduce the compact object’s mass by ～0.5–1 M _⊙, compared to the case with no heating. Based on the K0III spectral type of the optical component (with a probable mass of m _v ? 2.9 M _⊙), we concluded that m _x ? 2.45?2.55 M _⊙ (for i = 75°?90°). If the K0III star has lost a substantial part of its mass as a result of mass exchange, as in the V404 Cyg and GRS 1905+105 systems, and its mass is m _v ? 0.65?0.75 M _⊙, the compact object’s mass is close to the standard mass of a neutron star, m _x ? 1.4 M _⊙ (for i = 75°?90°). Thus, it is probable that the X-ray source in the 2S 0921-630 binary is an accreting neutron star.     

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.     

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.     

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.     

The mass of the compact object in the X-ray binary 4U 1700-37

The results of a systematic analysis of master radial-velocity curves for the X-ray binary 4U 1700-37 are presented. The dependence of the mass of the X-ray component on the mass of the optical component is derived in a Roche model based on a fit of the master radial-velocity curve. The parameters of the optical star are used to estimate the mass of the compact object in three ways. The masses derived based on information about the surface gravity of the optical companion and various observational data are 2.25 _?0.24 ^+0.23 M_⊙ and 2.14 _?0.56 ^+0.50 M_⊙. The masses based on the radius of the optical star, 21.9R_⊙, are 1.76 _?0.21 ^+0.20 M_⊙ and 1.65 _?0.56 ^+0.78 M_⊙. The mass of the optical component derived from the mass-luminosity relation for X-ray binaries, 27.4M_⊙, yields masses for the compact object of 1.41 _?0.08 ⁺ M_⊙ and 1.35 _?0.18 ^+0.18 M_⊙.     

Evolution of Stars Paired with Intermediate-Mass Black Holes

Under certain conditions, stars close to intermediate-mass black holes (IMBHs) can form close binary systems with these objects, in which the Roche lobe can be filled by the star and intense accretion of the star’s matter onto the IMBH is possible. Recently, accreting IMBHs have been associated with hyperluminous X-ray sources (HLXs), whose X-ray luminosities can exceed 10⁴¹ erg/s. In this paper, the evolution of star—IMBH binary systems is investigated assuming that the IMBH mainly accretes the matter of its companion star, and that the presence of gas in the vicinity of the IMBH does not appreciably affect changes in the orbit of the star. The computations take into account all processes determining the evolution of ordinary binary systems, as well as the irradiation of a star by hard radiation during the accretion of its matter onto the IMBH. The absorption of external radiation in the stellar envelope was calculated applying the same formalism that is used to calculate the opacity of the stellar matter. The computations also assumed that, if the characteristic time for the mass transfer is less than the thermal time scale of the star, there is no exchange betwween the orbital angular momentum of the system and the angular momentum of the matter flowing onto the IMBH.

Numerical simulations have shown that, under these assumptions, three types of evolution are possible for such a binary system, depending on the mass of the IMBH and the star, as well as on the star’s initial distance from the IMBH. The first type ends with the destruction of the star. For low-mass main sequence (MS) stars, only this option is realized, even in the case of large initial distances from IMBH. For massive MS stars, the star is also destroyed if the mass of the IMBH is high and the initial distance of the star from the IMBH is sufficiently small.

The second type of evolution can occur for massive MS stars, which are initially located farther from the IMBH than in the first type of evolution. In this case, the massive star fills its Roche lobe during its evolutionary expansion, after which a stage of intense mass transfer begins. It is in this phase of the evolution that the star- IMBH system can manifest itself as a HLX, when its X-ray luminosity L_X exceeds 10⁴¹ erg/s for a fairly long time. Numerical simulations show that the initial mass of the donor star in systems with M_BH = (10³?10⁵)M_⊙ must be close to ～10 M_⊙ in this case. The characteristic duration of the HLX stage is 30 000–70 000 years. For smaller initial donor masses close to ～5M_⊙, L_X does not reach 10⁴¹ erg/s in the stage of intense mass transfer, but can exceed 10⁴⁰ erg/s. The duration of this stage of evolution is 300 000–800 000 years. A characteristic feature of this second type of evolution is an increase in the orbital period of the system over time. As a result, after a period of intense mass loss, the star “retreats” inside the Roche lobe. A remnant of the star in the form of a white dwarf is left behind, and can end up fairly far from the IMBH.

The third type of evolution can occur for massive MS stars that are initially even farther from the IMBH, as well as for massive stars that are already evolved at the initial time. In this case, conservative mass exchange in the presence of intense stellar wind leads to the star moving away from the IMBH, without filling its Roche lobe at all. For massive stars with sufficiently strong stellar winds (for example, stars with masses ～50M_⊙), the accretion rate of matter onto the IMBH in this case can reach values that are characteristic of HLXs. As in the case of the second type of evolution, the stellar remnant can remain at a fairly large distance from the IMBH.

     

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.     

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.     

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.     

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.     

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.     

Star formation region in Orion KL. epoch 1985.8

We have derived the fine structure of the region of the H₂O supermaserflare in the Orion Nebula at epoch 1985.8. This structure includes a chain of compact components that extends to 25 AU and has a width of 0.4 AU. The velocities of the components vary along the chain. The structure corresponds to an accretion disk separated into protoplanetary rings, viewed edge-on. The velocities of the components correspond to Keplerian motion around an object with a mass of M=0.3±0.2M_⊙. The velocity of the central object relative to the Local Standard of Rest is V_LSR=4.0±0.7 km/s. The radius of the inner part of the disk is 9±4 AU, while the radius of the outer disk is 35±6 AU. The rotational velocities of the inner and outer rings are 5±1 km/s and 2.5±0.5 km/s, respectively. The emission of the structure is amplified in the ambient medium—an envelope with velocities of 7.6±0.3 km/s. The rate at which the envelope is accreting onto the central object is 3.6±0.7 km/s. The gradient of the infall velocity is 1.1 km/s.     

Late stages of the evolution of close compact binaries: Type I supernovae,gamma-ray bursts,and supersoft X-ray sources

We consider the evolution of close binaries resulting in the most intensive explosive phenomena in the stellar Universe—Type Ia supernovae and gamma-ray bursts. For Type Ia supernovae, which represent thermonuclear explosions of carbon-oxygen dwarfs whose masses reach the Chandrasekhar limit during the accretion of matter from the donor star, we derive the conditions for the accumulation of the limiting mass by the degenerate dwarf in the close binary. Accretion onto the degenerate dwarf can be accompanied by supersoft X-ray radiation with luminosity 1–10⁴ L _⊙. Gamma-ray bursts are believe to accompany the formation and rapid evolution of compact accretion-decretion disks during the formation of relativistic objects—black holes and neutron stars. The rapid (～1 M _⊙/s) accretion of matter from these disks onto the central compact relativistic star results in an energy release of ～0.1 M _⊙ c ² ～ 10⁵³ erg in the form of gamma-rays and neutrinos over a time of 0.1–1000 s. Such disks can form via the collapse of the rapidly rotating cores of Type Ib, Ic supernovae, which are components in extremely close binaries, or alternately due to the collapse of accreting oxygen-neon degenerate dwarfs with the Chandrasekhar mass into neutron stars, or the merging of neutron stars with neutron stars or black holes in close binaries. We present numerical models of the evolution of some close binaries that result in Type Ia supernovae, and also estimate the rates of these supernovae (～0.003/year) and of gamma-ray bursts (～10^?4/year) in our Galaxy for various evolutionary scenarios. The collimation of the gamma-ray burst radiation within an opening angle of several degrees “matches” the latter estimate with the observed rate of these events, ～10^?7–10^?8/year calculated for a galaxy with the mass of our Galaxy.