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.     

Dependence of the absorption-line profiles and radial-velocity curve of the optical star in an X-ray binary on the orbital inclination and component-mass ratio

Theoretical absorption-line profiles and radial-velocity curves for tidally deformed optical stars in X-ray binary systems are calculated assuming LTE. The variations in the profile shapes and radial-velocity curve of the optical star are analyzed as a function of the orbital inclination of the X-ray binary system. The dependence of the shape of the radial-velocity curve on the orbital inclination i increases with decreasing component-mass ratio q = m _x /m _v . The integrated line profiles and radial-velocity curves of the optical star are calculated for the Cyg X-1 binary, which are then used to estimate the orbital inclination and mass of the relativistic object: i < 43° andm _x = 8.2–12.8 M_⊙. These estimates are in good agreement with earlier results of fitting the radial-velocity curve of Cyg X-1 using a simpler model (i < 45°, m _x = 9.0–13.2 M_⊙).     

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

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 mass of the compact object in the X-ray binary her X-1/HZ her

We have obtained the first estimates of the masses of the components of the Her X-1/HZ Her X-ray binary system taking into account non-LTE effects in the formation of the H _γ absorption line: m _x = 1.8 M _⊙ and m _v = 2.5 M _⊙. These mass estimates were made in a Roche model based on the observed radial-velocity curve of the optical star, HZ Her. The masses for the X-ray pulsar and optical star obtained for an LTE model lie are m _x = 0.85 ± 0.15 M _⊙ and m _v = 1.87 ± 0.13 M _⊙. These mass estimates for the components of Her X-1/HZ Her derived from the radial-velocity curve should be considered tentative. Further mass estimates from high-precision observations of the orbital variability of the absorption profiles in a non-LTE model for the atmosphere of the optical component should be made.     

A “tsunami” in the protoplanetary disk of SV Cep

The results of ~15 years of photometric observations of the UX Ori star SV Cep in the near-infrared (JHKL) are presented. They demonstrate the presence of a cyclic component with a period of ~7 years in the variations of the IR fluxes. This is clearly seen in all four IR bands, but is absent in the optical. The variation amplitude is highest in the K band: ΔK ≈ 0.68 ^m . The shape of the variations differs slightly in the transition from J to L. However, it is reproduced with good accuracy during two cycles, suggesting a periodic process is observed. If the periodic perturbations in the circumstellar disk of SV Cep are due to a companion’s orbitalmotion, the orbital semi-major axis should be ~5AU, foramass of SVCep of 2.6M _⊙. The absence of a seven-year period in the optical light curve of SV Cep means that the observed period cannot be due to variations in the circumstellar extinction. The IR brightness variations could be due to the companion’s motion along an eccentric orbit, resulting in a periodic modulation of the rate of accretion onto the star.     

Parameters of V404 Cyg—A black-hole binary

We present the results of spectroscopic observations of the X-ray binary V404 Cyg obtained on the 6-m telescope of the Special Astrophysical Observatory in 2001–2002. We have used a statistical approach to interpret the radial-velocity curve of V404 Cyg. We derived the dependence of the mass of the X-ray emitting component m_x on the mass of the optical component m_v via an analysis of the radial-velocity curve based on profiles of the CaI 6439.075 Å absorption line synthesized in a Roche model. Using the orbital inclination estimated from the ellipticity of the optical component, i=54°–64°, and the component-mass ratio q=m_x/m_v=16.7 found from the rotational broadening of the spectral lines, we obtain m _s =10.65±1.95M_⊙ for the mass of the black hole.     

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.     

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

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

Evolution of the masses of neutron stars in binary systems

We study the growth of the masses of neutron stars in binary systems due to the accumulation of mass from the optical donors accreted onto the neutron-star surface. Possible scenarios for this accretion are considered. The masses and magnetic-field strengths of radio pulsars derived using population-synthesis methods are compared to the observational data. The population-synthesis analysis indicates that a neutron star can increase its mass from the standard value of m _x ? 1.35M_⊙ to the Oppenheimer-Volkoff limit, m _x ? 2.5M_⊙, via accretion from a companion.     

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.     

Radial pulsations of helium stars and a model for BX CIR

The results of hydrodynamical calculations of radially pulsating helium stars with masses 0.5M_⊙≤M≤0.9M_⊙, bolometric luminosities 600L_⊙≤5×10³L_⊙, and effective temperatures 1.5×10⁴ K≤T_eff≤3.5×10⁴ K are presented. The pulsation instability of these stars is due to the effects of ionization of iron-group elements in layers with temperatures T～2×10⁵ K. The calculations were carried out using opacities for the relative mass abundances of hydrogen and heavy elements X=0 and Z=0.01, 0.015, and 0.02. Approximate formulas for the pulsation constant Q over the entire range of pulsation instability of the hot helium stars in terms of the mass M, radius R, effective temperature T_eff, and heavy-element abundance Z are derived. The instability of BX Cir to radial pulsations with the observed period Π=0.1066 d occurs only for a mass M≥0.55M_⊙, effective temperature T_eff≥23000 K, and heavy-element abundance Z≥0.015. The allowed mass of BX Cir is in the range 0.55M_⊙≤M≤0.8M_⊙, which corresponds to luminosities 800L_⊙≤M≤1400L_⊙ and mean radii 1.7R_⊙?R?2.1R_⊙.     

PN G068.1+11.0: A young pre-cataclysmic variable with an extremely hot primary

An analysis of spectroscopic and photometric data for the young pre-cataclysmic variable (PCV) PN G068.1+11.0, which passed through its common-envelope stage relatively recently, is presented. The spectroscopic and photometric data were obtained with the 6-m telescope and Zeiss-1000 telescope of the Special Astrophysical Observatory. The light curves show sinusoidal brightness variations with the orbital-period time scale and brightness-variation amplitudes of Δm = 1_.^m41, 1_.^m62, and 1_.^m57 in the B, V, and R bands, respectively. The system’s spectrum exhibits weak HI (Hβ–Hδ) andHeII λλ4541, 4686, 5411 Å absorption lines during the phases of minimum brightness, as well as HI, HeII, CIII, CIV, NIII, and OII emission lines whose intensity variations are synchronized with variations of the integrated brightness of the system. The emission-line formation in the spectra can be fully explained by the effects of fluorescence of the ultraviolet light from the primary at the surface of the cool star. All the characteristics of the optical light of PN G068.1+11.0 confirm that it is a young PCV containing sdO subdwarf. The radial velocities were measured from a blend of lines of moderately light elements, CIII+NIII λ4640 Å, which is formed at the surface of the secondary due to reflection effects. The ephemeris of the system has been improved through a joint analysis of the radial-velocity curves and light curves of pre-cataclysmic variable, using modelling of the reflection effects. The fundamental parameters of PN G068.1+11.0 have been determined using two evolutionary tracks for planetary-nebula nuclei of different masses (0.7 M_⊙and 0.78M_⊙). The model spectra for the system and a comparison with the observations demonstrate the possibility of refining the components’ effective temperatures if the quality of the spectra used is improved.     

IR photometry of the symbiotic star BF Cyg: Detection of the red giant’s ellipsoidal brightness variability

We present the results of our IR photometric observations of the classical symbiotic star BF Cyg acquired in 1978–2003. The variability range in the J and K bands was ～0.2^m. A periodic component in the cool star’s brightness variations is clearly visible, its period being half the orbital one and its J amplitude being ～0.15^m. This component is associated with the ellipsoidal shape of the red giant, which model calculations show fills its Roche lobe. This is required in order to reproduce ellipsoidal brightness variability with such a large amplitude: the calculated amplitude for a red giant filling 90% of its Roche lobe is half the observed value. At the same time, it was not possible to confidently chose the optimum component-mass ratio, q = M _giant /M_hot, and orbital inclination, i, from possible values in the ranges q = 2–4, i = 70°–90°. Including the contribution from the hot radiation sources (the hot component and ionized envelope), which vary with a period equal to the orbital period, has a considerable influence on the estimated parameters associated with the red giant’s ellipsoidal brightness variations, and this contribution cannot be neglected. The deviations of the observed from the calculated light curve are irregular, with the rms deviation being σ(O-C) ≈ 0.04^m.     

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.