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.     

Evolution of close binaries and gamma-ray bursts

We analyze the late stages of evolution of massive (M ₀ ? 8 M _⊙) close binaries, from the point of view of possible mechanisms for the generation of gamma-ray bursts. It is assumed that a gamma-ray burst requires the formation of a massive (～1 M _⊙), compact (R ? 10 km) accretion disk around a Kerr black hole or neutron star. Such Kerr black holes are produced by core collapses of Wolf-Rayet stars in very close binaries, as well as by mergers of neutron stars and black holes or two neutron stars in binaries. The required accretion disks can also form around neutron stars that were formed via the collapse of ONeMg white dwarfs. We estimate the Galactic rate of events resulting in the formation of rapidly rotating relativistic objects. The computations were carried out using the “Scenario Machine.”     

Acceleration of spatial motions of stars by close-binary supermassive black holes in galactic nuclei

The conditions for the acceleration of the spatial motions of stars by close-binary supermassive black holes (SMBHs) in galactic nuclei are analyzed in order to derive the velocity distribution for stars ejected from galaxies by such black holes. A close binary system consisting of two SMBHs in circular orbits was subject to a spherically symmetrical “barrage” of solar-mass stars with various initial velocities. The SMBHs were treated as point objects with Newtonian gravitational fields. Models with binary component-mass ratios of 1, 0.1, 0.01, and 0.001 were studied. The results demonstrate the possibility of accelerating neutron stars, stellar-mass black holes, and degenerate dwarfs to velocities comparable to the relative orbital velocities of the binary-SMBH components. In the stage when the binary components are merging due to the action of gravitational-wave radiation, this velocity can approach the speed of light. The most massive binary black-holes (M ? 10⁹M_⊙) can also accelerate main-sequence stars with solar or subsolar masses to such velocities.     

Formation and evolution of Ae and Be stars

Several scenarios for the formation of accretion and decretion disks in single and binary Ae and Be stars are proposed. It is shown that, in order for a rapidly rotating main-sequence Be star to lose mass via a disk, the star’s rotation must be quasi-rigid-body. Estimates show that such rotation can be maintained by the star’s magnetic field, which is probably a relict field. The evolution of single Be main-sequence stars is numerically simulated allowing for mass loss via the stellar wind and rotational mass loss assuming rigid-body rotation. The stellar wind is the factor that determines the maximum mass of Be stars, which is close to 30M _⊙. The evolution of Be stars in close binaries is analyzed in the approximation adopted in our scenario. Long gamma-ray bursts can be obtained as a result of the collapse of rapidly rotating oxygen—neon degenerate dwarfs—the accreting companions of Be stars—into neutron stars.     

Formation of close binary black holes merging due to gravitational-wave radiation

The conditions for the formation of close-binary black-hole systems merging over the Hubble time due to gravitational-wave radiation are considered in the framework of current ideas about the evolution of massive close-binary systems. The original systems whose mergers were detected by LIGO consisted of main-sequence stars with masses of 30–100M _⊙. The preservation of the compactness of a binary black hole during the evolution of its components requires either the formation of a common envelope, probably also with a low initial abundance of metals, or the presence of a “kick”—a velocity obtained during a supernova explosion accompanied by the formation of a black hole. In principle, such a kick can explain the relatively low frequency of mergers of the components of close-binary stellar black holes, if the characteristic speed of the kick exceeds the orbital velocities of the system components during the supernova explosion. Another opportunity for the components of close-binary systems to approach each other is related to their possible motion in a dense molecular cloud.     

The evolutionary status of ultraluminous X-ray sources

We analyze models for quasi-stationary, ultraluminous X-ray sources (ULXs) with luminosities 10³⁸–10⁴⁰ erg/s exceeding the Eddington limit for a ～1.4M_⊙ neutron star. With the exception of relatively rare stationary ULXs that are associated with supernova remnants or background quasars, most ULXs are close binary systems containing a massive stellar black hole (BH) that accretes matter donated by a stellar companion. To explain the observed luminosities of ～10⁴⁰ erg/s, the mass of the BH must be ～40M_⊙ if the accreted matter is helium and ～60M_⊙ if the accreted matter has the solar chemical composition. We consider donors in the form of main-sequence stars, red giants, red supergiants, degenerate helium dwarfs, heavy disks that are the remnants of disrupted degenerate dwarfs, helium nondegenerate stars, and Wolf-Rayet stars. The most common ULXs in galaxies with active star formation are BHs with Roche-lobe-filling main-sequence companions with masses ～7M_⊙ or close Wolf-Rayet companions, which support the required mass-exchange rate via their strong stellar winds. The most probable candidate ULXs in old galaxies are BHs surrounded by massive disks and close binaries containing a BH and degenerate helium-dwarf, red-giant, or red-supergiant donor.     

The relation between the observed mass distribution for compact stars and the mechanism for supernova explosions

The observed mass distribution for the compact remnants of massive stars (neutron stars and black holes) and its relationship to the possible mechanism of ejection of the envelopes of type II and Ib/c supernovae are analyzed. The observed lack of compact remnants with masses 1.5–3 M _⊙ suggests a magneto-rotational mechanism for the supernovae, and a soft equation of state for neutron stars with limiting masses near 1.5 M _⊙. The observational consequences of this hypothesis are discussed.     

Type Ia supernovae: Non-standard candles of the universe

We analyze the influence of the evolution of light absorption by gray dust in the host galaxies of type Ia supernovae (SN Ia) and the evolution of the mean combined mass of close-binary carbon-oxygen white dwarfs merging due to gravitational waves (SN Ia precursors) on the interpretation of Hubble diagrams for SN Ia. A significant increase in the mean SN Ia energy due to the higher combined masses of merging dwarfs should be observable at redshifts z > 2. The observed relation between the distance moduli and redshifts of SN Ia can be interpreted not only as evidence for accelerated expansion of the Universe, but also as indicating time variations of the gray-dust absorption of light from these supernovae in various types of host galaxies, observational selection effects, and the decreasing mean combined masses of merging degenerate dwarfs.     

A model for the population of helium stars in the Galaxy: Low-mass stars

We model the Galactic ensemble of helium stars using population synthesis techniques, assuming that all helium stars are formed in binaries. In this picture, single helium stars are produced by mergers of helium remnants of the components of close binaries (mainly, the merging of helium white dwarfs) or in the disruption of binaries with helium components during supernova explosions. The estimated total birthrate of helium stars in the Galaxy is 0.043 yr^?1; the total number is 4 × 10⁶; and the binarity rate is 76%. We construct a subsample of low-mass (M_He ? 2M_⊙) helium stars defined by observational selection effects: the limiting magnitude (V_He ≤ 16), ratio of the magnitudes of the components in binaries (V_He ≤ V_comp), and lower limit for the semiamplitude of the radial velocity required for detecting binarity (K_min = 30 km s^?1). The parameters of this subsample are in satisfactory agreement with observations of helium subdwarfs. In particular, the binarity rate in the selection-limited sample is 58%. We analyze the relations between the orbital periods and masses of helium subdwarfs and their companions in systems with various combinations of components. We predict that the overwhelming majority (～97%) of unobserved companions to helium stars will be white dwarfs, predominantly, carbon-oxygen white dwarfs.     

The Role of Gravitational Radiation in the Evolution of Stars and Galaxies

The conditions for the formation of close binaries containing main-sequence stars, degenerate dwarfs of various types, neutron stars, and black holes of various masses are considered. The paper investigates the evolution of the closest binary systems under the influence of their gravitational-wave radiation. The conditions under which the binary components can merge on a time scale shorter than the Hubble time as a result of their emission of gravitational waves are estimated. A self-consistent scenario model is used to estimate the frequency of such events in the Galaxy, their observable manifestations, the nature of the merger products, and the role of these events in the evolution of stars and galaxies. The conditions for the formation and evolution of supermassive binary black holes during collisions andmergers of galaxies in their dense clusters are studied.

     

Wolf-Rayet stars,black holes,and gamma-ray bursters in close binaries

We consider the evolutionary status of observed close binary systems containing black holes and Wolf-Rayet (WR) stars. When the component masses and the orbital period of a system are known, the reason for the formation of a WR star in an initial massive system of two main-sequence stars can be established. Such WR stars can form due to the action of the stellar wind from a massive OB star (M_OB≥50M_⊙), conservative mass transfer between components with close initial masses, or the loss of the common envelope in a system with a large (up to ～25) initial component mass ratio. The strong impact of observational selection effects on the creation of samples of close binaries with black holes and WR stars is demonstrated. We estimate theoretical mass-loss rates for WR stars, which are essential for our understanding the observed ratio of the numbers of carbon and nitrogen WR stars in the Galaxy \(\dot M_{WR} (M_ \odot yr^{ - 1} ) = 5 \times 10^{ - 7} (M_{WR} /M_ \odot )^{1.3} \). We also estimate the minimum initial masses of the components in close binaries producing black holes and WR stars to be ～25M_⊙. The spatial velocities of systems with black holes indicate that, during the formation of a black hole from a WR star, the mass loss reaches at least several solar masses. The rate of formation of rapidly rotating Kerr black holes in close binaries in the Galaxy is ～3×10^?6 yr^?1. Their formation may be accompanied by a burst of gamma radiation, possibly providing clues to the nature of gamma-ray bursts. The initial distribution of the component mass ratios for close binaries is dN～dq=dM₂/M₁ in the interval 0.04?q₀≤1, suggesting a single mechanism for their formation.     

Evolution of close stellar binaries with black holes under the action of gravitational radiation and magnetic and induced stellar winds of the donor

We show that semi-detached close binary systems with massive (4–25M_⊙) black holes are formed in the evolution of massive stellar binaries in which the initial mass of the primary exceeds ～25M_⊙. The mass exchange in such systems is maintained by the nuclear evolution of the donor and by its magnetic and induced stellar winds. The donor in such systems can be a main-sequence star, subgiant, non-degenerate helium star, or white dwarf. The evolution of corresponding systems with black-hole masses of 10M_⊙ is investigated.     

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.     

The nature hypervelocity stars

We list and analyze the main currently known mechanisms for accelerating the space motions of stars. A high space velocity of a star can be a consequence of its formation in the early stages of the evolution of a massive galaxy, when it was spheroidal and non-stationary, so that stars were born with velocities close to the escape velocity for the galaxy. Another possibility is that the star arrived from another galaxy with a velocity that is high for our Galaxy. The decay of unstable close multiple stars or supernova explosions in close binaries can also provide velocities of up to several hundreds of km/s to main-sequence stars and velocities of up to ∼1000 km/s to degenerate stars, neutron stars, and stellar-mass black holes. The merger of components of a binary system containing two neutron stars or a neutron star and a black hole due to gravitational-wave radiation can accelerate the nascent black hole to a velocity∼1000 km/s. Hypervelocity relativistic stars can be born due to asymmetric neutrino ejection during a supernova explosion. Stars can be efficiently accelerated by single and binary supermassive black holes (with masses from several millions to several billions of solar masses) in the nuclei of galaxies. Thanks to their gravitational field and fast orbital motion (in the case of binary objects), supermassive black holes are able to accelerate even main-sequence stars to relativistic velocities.     

The origin of the rotational and spatial velocities of radio pulsars

We analyze possible origins of the observed high rotational and spatial velocities of radio pulsars. In particular, these can be understood if all radio pulsars originate in close binary systems with orbital periods of 0.1–100 days, with the neutron star being formed by a type Ib,c supernova. The high spatial velocities of pulsars (v _p up to 1000 km/s) reflect the high Keplerian velocities of the components of these binaries, while their short periods of rotation (P _p < 4 s) are due to the rapid rotation of the presupernova helium-star components with masses of 2.5–10 M_⊙, which is synchronous with their orbital rotation. Single massive stars or components in wide binaries are likely to produce only slowly rotating (P _p > 4 s) neutron stars or black holes, which cannot be radio pulsars. As a result, the rate of formation of radio pulsars should be a factor of a few lower than the rate of type II and type Ib,c supernovae estimated from observations. This scenario for the formation of radio pulsars is supported by (i) the bimodal spatial velocity distribution of radio pulsars; (ii) the coincidence of the observed spatial velocities of radio pulsars with the orbital velocities of the components of close binaries with nondegenerate helium presupernovae; (iii) the correlation between the orbital and rotational periods for 22 observed radio pulsars in binaries with elliptical orbits; and (iv) the similarity of the observed rate of formation of radio pulsars and the rate of type Ib,c supernovae.     

Models for highly flattened, rapidly rotating cool stars in a Newtonian approximation

We investigate the physical characteristics of single, rapidly rotating white dwarfs, which could form as a result of a merger of two white dwarfs with different masses and filled Roche lobes, due to the radiation of gravitational waves. When the merging of the binary components occurs without loss of mass and angular momentum, the merger products are subject to secular instability, and the density in their cores does not exceed ～10⁸ g/cm³. Models are constructed for rapidly rotating neutron stars, which could form after the collapse of rotating iron cores of evolved massive stars. Dynamically unstable neutron-star models are characterized by a shift of the maximum density from the rotational axis. The total momentum of such neutron stars is about half the maximum possible momentum for the evolved cores of massive stars.     

Rotational mass loss by Be stars

We have carried out a numerical study of rotational mass loss by rapidly rotating Be stars assuming preservation of rigid-body rotation during their main-sequence evolution. Evolutionary models are computed for stars with solar chemical composition and initial masses of 3, 10 and M_⊙. As a result of their rapid initial rotation, these stars can lose one to four percent of their initial mass during the main-sequence stage. The amount of mass lost increases with the initial mass of the star. The matter lost by Be stars can form gas-dust disks with masses comparable to the masses of planets, which, in principle, makes possible the formation of planetary systems around such stars.     

The evolutionary status of the most massive WNh stars in close binary systems. NGC 3603-A1

The evolution of the components of the unique, massive, close binary system NGC 3603-A1, which consists of stars of spectral types WN6ha and WN6h, is analyzed. The component masses are estimated to be 116 and 89M _⊙, close to the highest measured stellar masses. Numerical modeling of the evolution of the components has been carried out, taking into account mass loss via the stellar winds of the two massive stars. It is shown that the maximum possible initial component masses are close to 140 and 125M _⊙. The components are currently slightly evolved main-sequence stars, with a comparative low degree of helium enrichment at their surfaces. Further evolution of the system will lead to filling of the Roche lobe of the primary and subsequent evolution in a common envelope. This may lead to the merger of the components, with the evolution of the system ending in the formation of a singlemassive black hole after the second supernova explosion. Otherwise, depending on the masses of the resulting black holes, either a binary system of two black holes or two unbound black holes may form, accompanied by gamma-ray bursts.     

Evolution of the components of the close binary WR 20a

We have analyzed the evolution of the components of the unique massive binary system WR 20a, which consists of a Wolf-Rayet nitrogen star and an Of star with an extremely small separation. The estimated masses of the components are 83 and 82 M _⊙, which are among the highest stellar mass inferred. We have carried out numerical modeling of the evolution of the components, taking into account the mass loss due to the stellar wind inherent to massive stars. In a scenario in which the systemis detached from the time the components reach the main sequence until its present state, the initial component masses are inferred to be close to 110 M _⊙, if the initial masses of the stars were equal, or 120 and 100 M _⊙, if they were different. Currently, the components are evolved main-sequence stars, whose surfaces are relatively little enriched by helium. The further evolution of the system will result in one of the components filling its Roche lobe and evolution within a common envelope. As a result, the components may coalesce, leading to the formation of a single massive black hole the supernova explosion. Otherwise, depending on the masses of the resulting black holes, either a binary system with two black holes or two free black holes will be formed. In the latter case, gamma-ray bursts will be observed.     

Wolf-Rayet stars and relativistic objects: Distinctions between the mass distributions in close binary systems

The observed properties of Wolf-Rayet stars and relativistic objects in close binary systems are analyzed. The final masses M _CO ^f for the carbon-oxygen cores of WR stars in WR + O binaries are calculated taking into account the radial loss of matter via stellar wind, which depends on the mass of the star. The analysis includes new data on the clumpy structure of WR winds, which appreciably decreases the required mass-loss rates $\dot M_{WR}$ for the WR stars. The masses M _CO ^f lie in the range (1–2)M _⊙–(20–44)M _⊙ and have a continuous distribution. The masses of the relativistic objects M _x are 1–20M _⊙ and have a bimodal distribution: the mean masses for neutron stars and black holes are 1.35 ± 0.15M _⊙ and 8–10M _⊙, respectively, with a gap from 2–4M _⊙ in which no neutron stars or black holes are observed in close binaries. The mean final CO-core mass is $\overline M _{CO}^f = 7.4 - 10.3M_ \odot$ , close to the mean mass for the black holes. This suggests that it is not only the mass of the progenitor that determines the nature of the relativistic object, but other parameters as well-rotation, magnetic field, etc. One SB1R Wolf-Rayet binary and 11 suspected WR + C binaries that may have low-mass companions (main-sequence or subgiant M-A stars) are identified; these could be the progenitors of low-mass X-ray binaries with neutron stars and black holes.