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.     

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.     

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.     

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

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.     

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.     

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 evolution of GRS 1915+105—an X-ray binary with a black hole

A brief review of the observed parameters of binary systems with black holes is presented. We discuss in detail the evolutionary status of the X-ray binary GRS 1915+105, which contains a massive black hole. Numerical simulations of the evolution of GRS 1915+105 at the X-ray stage indicate that the most probable initial mass of the optical component (donor star) is (1.5–)M_⊙. Two possible scenarios are suggested for the evolution of the system prior to the formation of the black hole. If the initial mass of the optical component was (2.5–)M_⊙, the system underwent a common-envelope phase; in this case, the initial mass of the black hole progenitor did not exceed ～50M_⊙. If the initial mass of the donor was (1.5–2.5)M_⊙, a scenario without a common envelope is possible, with the initial mass of the black hole progenitor being smaller than ～50M_⊙. The lack of information about the initial mass-ratio distribution for binary components for small q and the uncertainty of the system parameters make it impossible to give preference to a particular scenario for the system's prior evolution.     

Masses of the visual components and black holes in X-ray novae: Effects of proximity of the components

It is shown that the approximation of the complex, tidally distorted shape of a star as a circular disc with local line profiles and a linear limb-darkening law, which is usually applied when deriving equatorial stellar rotation velocities from line profiles, leads to overestimation of the equatorial velocity V _rot sin i and underestimation of the component mass ratio q = M _x /M _v . A formula enabling correction of the effect of these simplifying assumptions on the shape of a star is used to re-determine the mass ratios q and the masses of the black holes M _x and visual components M _v in low-mass X-ray binary systems containing black holes. Taking into account the tidal–rotational distortion of the stellar shape can significantly increase the mass ratios q = M _x /M _v , reducing M _v , while M _x changes only slightly. The resulting distribution of M _v attains its maximum near M _v ? 0.35M _⊙, in disagreement with the results of population synthesis computations realizing standard models for Galactic X-ray novae with black holes. Possible ways to overcome this inconsistency are discussed. The derived distribution of M _x also differs strongly from the mass distribution for massive stars in the Galaxy.     

Neutron stars in close binaries with elliptical component orbits

We consider the evolution of close binaries in which the initial secondary component is a nondegenerate helium star with mass M_He = 0.4–60 M_⊙, while the initially more massive primary has evolved into a black hole, neutron star, or degenerate dwarf. The neutron star is assumed to originate as a result of the evolution of a helium star with a mass of 2.5 M_⊙ ≤ M_He ≤ 10 M_⊙ after the explosion of a type Ib,c supernova. If the axial rotation of the helium star before the explosion is rigid-body and synchronized with the orbital rotation, for P_orb ≤ 0.16 day, the rotational energy of the young neutron star will exceed the energy of an ordinary supernova. If the magnetic field of the neutron star is sufficiently strong, the necessary conditions for a magnetic-rotational supernova are provided. The initial rotational period of a young neutron star originating in a system with an orbital period shorter than ～50 days is shorter than ～4 s, which, according to observations, is required for the appearance of a radio pulsar. A helium star whose mass exceeds ～10 M_⊙ in a close binary with an orbital period shorter than one day and with the axial rotation of the helium presupernova synchronous with the orbital rotation evolves into a Kerr black hole, whose formation is likely to be accompanied by a gamma-ray burst with a duration longer than two seconds. In particular, we consider close binaries in which the second supernova results in the formation of a neutron star that remains in the binary. The theoretical distribution of orbital periods and eccentricities for such systems is consistent with that observed for radio pulsars in the Galactic disk in binaries with compact components and orbital eccentricities exceeding ～0.09, providing an explanation for the observed correlation between the orbital eccentricities and orbital periods for these systems.     

The probable quadruple systems FK Aql and FZ Del

Orbital-period variations of the eclipsing binaries FK Aql and FZ Del are analyzed. For each of the systems, a superposition of two cyclic variations of their orbital periods is found. FK Aql may be a quadruple system that contains two more bodies, besides the eclipsing binary, with masses M ₃ ? 1.75M _⊙ and M ₄ ? 1.47M _⊙, and the corresponding periods 15 and 82 yrs. This could also be a triple system with a third body of mass M ₃ ? 1.75M _⊙ and a period of the long-period orbit P ₃ = 15 yrs, or with a third body of mass M ₃ ? 1.30M _⊙ and a period of the long-period orbit P ₃ = 82 yrs. FZ Del may be a quadruple system with the additional componentmasses M ₃ ? 0.2M _⊙ and M ₄ ? 0.3M _⊙, with the periods 10.2 and 53.7 yrs. This could also be a triple system with a third-body mass M ₃ ? 0.2M _⊙ and a period of the long-period orbit P ₃ = 10.2 yrs. In both systems, the residual period variations could be due to magnetic cycles of the secondary. The period variations of the eclipsing binary FZ Del could also be due to apsidal motion, together with the influence of a third body or the effects of magnetic activity.     

Estimation of the accuracy of methods for determining component masses for low-mass X-ray binary systems

Modern modeling of the population of low-mass X-ray binary systems containing black holes applying standard assumptions leads to a lack of agreement between the modeled and observed mass distributions for the optical components, with the observed masses being lower. This makes the task of estimating the systematic errors in the derived component masses due to imperfect models relevant. To estimate the influence of systematic errors in the derived masses of stars in X-ray binary systems, we considered two approximations for the tidally deformed star in a Roche model. Approximating the star as a sphere with a volume equal to that of the Roche lobe leads to slight overestimation of the equatorial rotational velocity V _rot sin i, and hence to slight underestimation of the mass ratio q = M _x /M _v . Approximating the star as a flat, circular disk with constant local line profiles and a linear limb-darkening law (a classical rotational broadeningmodel) is an appreciably cruder approach, and leads to overestimation of V _rot sin i by about 20%. In the case of high values of q = M _x /M _v , this approximation leads to substantial underestimation of the mass ratio q, which can reach several tens of percent. The mass of the optical star is overestimated by a factor of 1.5 in this case, while the mass of the black hole is changed only slightly. Since most estimates of component mass ratios for X-ray binary systems are carried out using a classical rotational broadening model for the lines, this leads to the need for appreciable corrections to (reductions of) previously published masses for the optical stars, which enhances the contradiction with the standard evolutionary scenario for low-mass X-ray binaries containing black holes.     

The formation and evolution of the most massive stars

We consider the formation of massive stars under the assumption that a young star accretes material from the protostellar cloud through its accretion disk while losing gas in the polar directions via its stellar wind. The mass of the star reaches its maximum when the intensity of the gradually strengthening stellar wind of the young star becomes equal to the accretion rate. We show that the maximum mass of the forming stars increases with the temperature of gas in the protostellar cloud T ₀, since the rate at which the protostellar matter is accreted increases with T ₀. Numerical modeling indicates that the maximum mass of the forming stars increases to ～900 M _⊙ for T ₀ ～ 300 K. Such high temperatures of the protostellar gas can be reached either in dense star-formation regions or in the vicinity of bright active galactic nuclei. It is also shown that, the lower the abundance of heavy elements in the initial stellar material Z, the larger the maximum mass of the star, since the mass-loss rate due to the stellar wind decreases with decreasing Z. This suggests that supermassive stars with masses up to 10⁶ M _⊙ could be formed at early stages in the evolution of the Universe, in young galaxies that are almost devoid of heavy elements. Under the current conditions, for T ₀ = (30–100) K, the maximum mass of a star can reach ～100M _⊙, as is confirmed by observations. Another opportunity for the most massive stars to increase their masses emerges in connection with the formation and early stages of evolution of the most massive close binary systems: the most massive stars can be produced either by coalescence of the binary components or via mass transfer in such systems.     

The mass of the black hole in the X-ray binary M33 X-7 and the evolutionary status of M33 X-7 and IC 10 X-1

We have analyzed the observed radial-velocity curve for the X-ray binary M33 X-7 in a Roche model. We have analyzed the dependence between the component masses and the degree of filling of the optical star’s Roche lobe to obtain the ratio of the masses of the optical star and compact object. For the most probable mass of the optical star, m _v = 70 M⊙, the mass of the compact object is m _x = 15.55 ± 3.20 M⊙. It has been shown that black holes with masses of mx = 15 M⊙ and even higher can form in binaries. We present characteristic evolutionary tracks for binary systems passing through an evolutionary stage with properties similar to M33 X-7-type objects. According to population-synthesis analyses, such binaries should be present in galaxies with masses of at least 10¹¹ M⊙. The present number of such systems in M33 should be of the order of unity. We have also studied the evolutionary status of the X-ray binary IC 10 X-1 with a Wolf-Rayet component, which may contain a massive black hole. The final stages of the evolution of the M33 X-7 and IC 10 X-1 systems should be accompanied by the radiation of gravitational waves.     

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

Photometry and spectroscopy of the system V4641 sagittarii in quiescence

We present the results of our CCD photometric and moderate-dispersion spectroscopic observations of the binary system V4641 Sgr, which contains a black hole of mass ≈9.5M_⊙ and a normal B9III star. The photometric light curve reveals an ellipticity effect with very high amplitudes in V and R, 0.40^m and 0.37^m, and the color curve shows that the surface temperature is nonuniform. All this testifies to tidal distortion of the normal star's surface due to the massive companion and to a high inclination of the orbit to the line of sight. In June and July 2002, during quiescence, we obtained data during three flares with amplitudes up to 0.26^m. In particular, spectroscopic observations were acquired near the time of the black hole's inferior conjunction. One hour before conjunction, a depression by EW=0.5 Å was observed in the red wing of the Hα absorption line, interpreted as absorption by gas flowing in the direction from the observer toward the normal star. This flow is apparently associated with a rarefied gas disk around the black hole, and the conjunction grazes the stellar surface if the orbital inclination is close to 70.7°. The maximum velocity along a circular Keplerian orbit is 650 km/s at a distance of R=0.15–0.20a from the black hole (where a is the component separation). Thus, we find the mass of the black hole to be M _BH =7.1–9.5M_⊙, confirming the model of Orosz et al. (2001).     

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.     

Nonlinear radial pulsations of Wolf-Rayet stars

The evolution of Population I stars with initial masses 60 M _⊙ ≤ M _ZAMS ≤ 120 M _⊙ is computed up to the Wolf-Rayet stage, when the central helium abundance decreases to Y _c ≈ 0.05. Several models from evolutionary sequences in the core helium-burning stage were used as initial conditions when solving the equations of radiative hydrodynamics for self-exciting stellar radial pulsations. The low-density envelope surrounding the compact core during the core helium burning is unstable against radial oscillations in a wide range of effective temperatures extending to T _eff ～ 10⁵ K. The e-folding time of the amplitude growth is comparable to the dynamical time scale of the star, and, when the instability ceases growing, the radial displacement of the outer layers is comparable to the stellar radius. Evolutionary changes of the stellar radius and luminosity are accompanied by a decrease in the amplitude of radial pulsations, but, at the effective temperature T _eff ≈ 10⁵ K, the stellar oscillations are still nonlinear, with a maximum expansion velocity of the outer layers of about one-third the local escape velocity. The period of the radial oscillations decreases from 9 hr to 4 min as stellar mass decreases from M = 28 M _⊙ to M = 6 M _⊙ in the course of evolution. The nonlinear oscillations lead to a substantial increase of the radii of the Lagrangian mass zones compared to their equilibrium radii throughout the instability region. The instability of Wolf-Rayet stars against radial oscillations is due to the action of the κ mechanism in the iron-group ionization zone, which has a temperature of T ～ 2 × 10⁵ K.     

Orbital period variations of the eclipsing binaries TU Cnc,VZ Leo,and OS Ori

Variations of the orbital periods of the eclipsing binaries TU Cnc, VZ Leo, and OS Ori are analyzed. Secular period decreases were earlier believed to occur in these systems. It is demonstrated that the period variations of TU Cnc can be represented using the light-time effect corresponding to the orbital motion of the eclipsing binary with a period of 78.6 years around the center ofmass of the triple system, with the mass of the third body being M ₃ > 0.82M _⊙. With the same accuracy, the period variations of VZ Leo and OS Ori can be represented either solely using the light-time effect, or a superposition of a secular period decrease and the light-time effect. For VZ Leo, the period of the long-term orbit is 63.8 years in the former case and 67.9 years in the latter case. Similar masses for the third body are indicated in both cases: M ₃ > 0.55M _⊙ and M ₃ > 0.61M _⊙. For OS Ori, the period of the long-term orbit is 46 years and M ₃ > 0.5M _⊙ in the former case, and the period is 36 years and M ₃ > 0.6M _⊙ in the latter case.     

The mass of the black hole in the X-ray binary LMC X-1

A dynamical estimate of the mass of the black hole in the LMC X-1 binary system is obtained in the framework of a Roche model for the optical star, based on fitting of the He I 4471 Å and He II 4200 Å absorption lines assuming LTE. The mass of the black hole derived from the radial-velocity curve for the He II 4200 Å line is m_x = 10.55 M_⊙, close to the value found earlier based on a model with two point bodies [1].