Formation of low-mass X-ray novae with black holes from triple systems

We apply a population synthesis technique to study the formation and evolution of low-mass X-ray binaries with black holes, observed as X-ray novae, from hierarchical triple systems. A scenario is suggested in which an inner close binary system evolves into an X-ray system with a large mass ratio. The high rate of accretion onto the neutron star leads to a common envelope stage, which may result in the formation of a Thorne-Zytkow (TZ) object. During its evolution, the envelope of the TZ object expands, encompassing the third star. The recurrent common-envelope stage decreases the size of the orbit of the third star, leading to the formation of a lowmass X-ray nova with a black hole. The dynamical stability of triple systems automatically ensures that only lowmass X-ray novae form. We also consider the possible formation of an X-ray nova from a binary in the case of asymmetrical core collapse during a supernova explosion.     

Close binaries containing Supermassive Black Holes

We consider the evolution of binary systems formed by a Supermassive Black Hole (SMBH) residing in the center of a galaxy or a globular cluster and a star in its immediate vicinity. The star is assumed to fill its Roche lobe, and the SMBH accretes primarily the matter of this star. The evolution of such a system is mainly determined by the same processes as for an ordinary binary. The main differences are that the donor star is irradiated by hard radiation emitted during accretion onto the SMBH; in a detached system, nearly all the donor wind is captured by the black hole, which strongly affects the evolution of the semi-major axis; it is not possible for companions of the most massive SMBHs to fill their Roche lobes, since the corresponding orbital separations are smaller than the radius of the last stable orbit in the gravitational field of the SMBH. Moreover, there may not be efficient exchange between the orbital angular momentum and the angular momentum of the overflowing matter in such systems. Our computations assumed that, if the characteristic timescale for mass transfer is smaller than the thermal timescale of the star, no momentum exchange occurs. Absorption of incident external radiation in the stellar envelope was treated using the same formalism that was used when computing the radiative transfer in the stellar atmosphere. Numerical simulations show that Roche-lobe overflow is possible for a broad range of initial system parameters. The evolution of semi-detached systems containing a star and a SMBH nearly always ends with the dynamical disruption of the star. Stars with masses close to the solar mass are destroyed immediately after they fill their Roche lobes. During the accretion of matter of disrupted stars, the SMBH can achieve quasar luminosities. If the SMBH accretes ambient gas as well as gas stripped from stars, the star is subject to additional radiation in the detached phase of its evolution, strengthening its stellar wind. This leads to an increase of the semi-major axis and subsequent decrease of the probability of Roche-lobe overflow during the subsequent evolution of the system.     

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.     

Signs of magnetic accretion in X-ray pulsars

The spin-down mechanism of accreting neutron stars is discussedwith an application to one of the best studied X-ray pulsars GX301-2. We show that the maximum possible spin-down torque applied to a neutron star from the accretion flow can be evaluated as K _sd ^(t) = ??²/(r _m r _cor)^3/2. The spin-down rate of the neutron star in GX301-2 can be explained provided the magnetospheric radius of the neutron star is smaller than its canonical value. We calculate the magnetospheric radius considering the mass-transfer in the binary system in the frame of the magnetic accretion scenario suggested by V.F. Shvartsman. The spin-down rate of the neutron star expected within this approach is in a good agreement with that derived from observations of GX301-2.     

Gamma-ray bursts—tracers of the history of star formation in the Universe

The rate of gamma-ray bursts (GRBs) in the Galaxy is estimated assuming that these events result from the formation of rapidly rotating Kerr black holes during the core collapse of massive, helium, Wolf-Rayet secondary components in very close binary systems. This process brings about rapid rotation of the cores of such Wolf-Rayet stars, inevitably resulting in the formation of Kerr black holes during type Ib,c supernovae. The current rate of formation of Kerr black holes (GRBs) in the Galaxy is about 3×10^?5/year. Collimation of the gamma-ray radiation into a small solid angle (about 0.1–0.01 sr) brings this rate into consistency with the observed rate of GRBs, estimated to be 10^?6–10^?7/year. Possible immediate progenitors of GRBs are massive X-ray binaries with X-ray luminosities of 10³⁸–10⁴⁰ erg/s. Due to the short lifetimes of the progenitors and the very high brightnesses of GRBs, the GRB rate can provide information about the history of star formation in the Universe on the Hubble time scale. A model in which the star-formation rate is determined by the conditions for ionization of the interstellar gas, whose density and volume are determined by supernovae, yields a Galactic star-formation history that can be viewed as representing the history of star formation in the Universe. The theoretical history of star formation is in satisfactory agreement with the history reconstructed from observations. The theoretical model for the history of star formation in the Galaxy can also be used to assess the influence of dust on optical observations of supernovae and GRBs in galaxies of various ages.     

Did the SN 1987A outburst leave a compact remnant?

We have computed the ejection of a massive envelope by a star during a type II supernova explosion in the presence of a compact remnant (a neutron star or black hole). This problem is of interest because of the possible presence of a compact remnant following the SN 1987A explosion. The computations demonstrate that a fairly large amount of matter is left in the neighborhood of the compact gravitating body. We present computations of the accretion rate onto the surface of the compact remnant. The estimated luminosity exceeds that observed for SN 1987A in various frequency ranges by several orders of magnitude.     

A prolonged flare in the blazar 3C 454.3

We present an analysis of data from multi-frequency monitoring of the blazar 3C 454.3 in 2010–2012, when the source experienced an unusually prolonged flare with a duration of about two years. This corresponds to the orbital period of the companion in a scenario in which two supermassive black holes are present in the nucleus of 3C 454.3. The flare’s shape, duration, and amplitude can be explained as a result of precession, if the plane of the accretion disk and the orbital plane of the binary are coincident. We detected small-scale structure of the flare, on time scales of no more than a month. These features probably correspond to inhomogeneities in the accretion disk and surrounding regions, with sizes of the order of 10¹⁵ cm. We estimated the size of the accretion disk based on the dynamical and geometrical parameters of this binary system: its diameter is comparable to the size of the orbit of the supermassive binary black hole, and its thickness does not exceed the gravitational radius of the central black hole. The presence of characteristic small-scale features during the flare makes it possible to estimate the relative time delays of variations in different spectral ranges: from gamma-ray to millimeter wavelengths.     

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

Searching for Observational Evidence for Binary Star Systems in Gamma-ray Bursts

Astronomy Reports - We search for observational evidence that is in agreement with the model of a binary star system in gamma-ray bursts (GRBs), as presented in [1–3]. The model predicts...     

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.

     

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.     

Formation of planets during the evolution of single and binary stars

Current views of the origin and evolution of single and binary stars suggest that the planets can form aroundmain-sequence single and binary stars, degenerate dwarfs, neutron stars, and stellarmass black holes according to several scenarios. Planets can arise during the formation of a star mainly due to excess angular momentum leading to the formation of an accretion-decretion disk of gas and dust around a single star or the components of a binary. It is the evolution of such disks that gives rise to planetary systems. A disk can arise around a star during its evolution due to the accretion of matter from dense interstellar clouds of gas and dust onto the star, the accretion of mass froma companion in a binary system, and the loss of matter during the contraction of a rapidly rotating star, in particular, if the star rotates as a rigid body and the rotation accelerates with its evolution along the main sequence. The fraction of stars with planetary systems is theoretically estimated as 30–40%, which is close to the current observational estimate of ∼34%.     

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.     

Light-curve syntheses for close binary systems: Modeling spiral waves in an elliptical disk around a white dwarf

We present an algorithm for synthesizing the light curve of a close binary consisting of a normal star (a red dwarf that fills its Roche lobe) and a spherical star (a white dwarf). The spherical component is surrounded by an elliptical accretion disk with a complex shape: it is geometrically thin near the spherical star and geometrically thick at the edge of the disk. An additional complication is presented by the presence of a one-or two-armed spiral pattern at the inner surface of the disk. The maximum height of the spiral arm above the disk surface is located at ～9 R _d , and the height decreases exponentially as the arm approaches the inner regions of the disk. Shielding of the inner hot parts of the disk by the crests of the spirals results in the formation of “steps” in out-of-eclipse parts of the orbital light curves. The algorithm takes into account the presence of a “hot line” by the lateral surface of the disk, making it possible to model binary systems in both quiescence and outburst. In the latter case, the hot line degenerates into a small bulge at the outer lateral surface of the disk, which can be considered an analog of a hot spot. The algorithm was applied to the orbital light curve of the cataclysmic binary IP Peg during its October 30, 2000, outburst. To explain the variations of the out-of-eclipse brightness of the system during the outburst, it is necessary to include the presence of a one-armed spiral wave at the inner surface of the disk, close to the periastron of the elliptical disk. We have obtained the parameters of IP Peg during the outburst for various models of the system.     

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.     

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

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.     

VLBI of supernovae and gamma-ray bursts

Supernovae and gamma-ray bursts (GRBs) are among the brightest events in the universe. Excluding Type Ia supernovae and short GRBs, they are the result of the core collapse of a massive star with material being ejectedwith speeds of several 1000 km/s to nearly the speed of light, and with a neutron star or a black hole left over as the compact remnant of the explosion. Synchrotron radiation in the radio is generated in a shell when the ejecta interact with the surrounding medium and possibly also in the central region near the compact remnant itself. VLBI has allowed resolving some of these sources and monitoring their expansion in detail, thereby revealing characteristics of the dying star, the explosion, the expanding shock front, and the expected compact remnant. We report on updates of some of the most interesting results that have been obtained with VLBI so far. Movies of supernovae are available from our website. They show the evolution from shortly after the explosion to decades thereafter, in one case revealing an emerging compact central source, which may be associated with shock interaction near the explosion center or with the stellar corpse itself, a neutron star or a black hole.     

Synthetic light curves of close binaries in a cool-disk model. The “hot line” and “hot spot” as visual indicators of interaction between the flow and disk

We present a “combined” model taking into account visual manifestations of the interaction between the gas flow and the accretion disk in a close binary system in the form of a “hot line” and a “hot spot.” The binary consists of a red dwarf that fills its Roche lobe and a compact spherical star (a white dwarf or neutron star) surrounded with a thick ellipsoidal accretion disk of a complex shape. The disk thickness is not large near the compact star but increases according to a parabolic law towards its outer edge. The oblique collision of the gaseous flow with matter of the cool, rotating disk, whose outer edge has a temperature <10 000 K, creates an extended region of enhanced energy release. In the combined model, this region is represented with a hot line that coincides with the optically opaque part of the flow and is located outside the disk, together with a hot spot at the outer surface of the disk, on the leeward side of the flow. The synthetic light curves for the combinedmodel and a hot-line model demonstrate that both models are able to fairly accurately reproduce the shapes of both classical and atypical light curves of cataclysmic variables in quiescence. Our determination of the parameters of the cataclysmic variable OY Car from an analysis of its light curves using the two models shows that the basic characteristics of the close binary, such as the component mass ratio q = M ₁/M ₂, orbital inclination i, effective temperatures of the red dwarf (T ₂) and white dwarf (T ₁), and orientation of the disk α _e , remain the same within the errors. The parameters describing the size of the slightly elliptical disk and the radiation flux from the disk differ by several percent (∼ 2–8%). A more significant difference is detected in the parameters of the hot line, due to the different shape and alignment of the flow and the presence of an additional radiation source—the hot spot—on the disk.