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

Merging of components in close binaries: Type Ia supernovae,massive white dwarfs,and Ap stars

The “Scenario Machine” (a computer code designed for studies of the evolution of close binaries) was used to carry out a population synthesis for a wide range of merging astrophysical objects: main-sequence stars with main-sequence stars; white dwarfs with white dwarfs, neutron stars, and black holes; neutron stars with neutron stars and black holes; and black holes with black holes. We calculate the rates of such events, and plot the mass distributions for merging white dwarfs and main-sequence stars. It is shown that Type Ia supernovae can be used as standard candles only after approximately one billion years of evolution of galaxies. In the course of this evolution, the average energy of Type Ia supernovae should decrease by roughly 10%; the maximum and minimum energies of Type Ia supernovae may differ by no less than by a factor of 1.5. This circumstance must be taken into account at estimating the parameters of the Universe expansion acceleration. According to theoretical estimates, the most massive—as a rule, magnetic—white dwarfs probably originate from mergers of white dwarfs of lower mass. At least some magnetic Ap and Bp stars may form in mergers of low-mass main-sequence stars (M ? 1.5 M _⊙) with convective envelopes.     

Axially symmetrical models of the products of a conservative merger of a binary white dwarf with similar-mass components

The possibility of a conservative merger of a binary white dwarf whose components have similar masses is studied. Axially symmetrical models for single, rapidly rotating white dwarfs that are possible products of such mergers are constructed and their physical characteristics investigated. The merger products must be turbulent, and the viscosity of the electron gas is not sufficient to support the observed luminosities of massive, bright white dwarfs. The amount of dissipative energy and the timescale for its release are estimated.     

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.     

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.     

Equilibrium Configurations of Rotating White Dwarfs at Finite Temperatures

In this work, cold and hot, static and rotating white dwarf stars are investigated within the framework of classical physics, employing the Chandrasekhar equation of state. The main parameters of white dwarfs such as the central density, pressure, total mass and radius are calculated fulfilling the stability criteria for hot rotating stars. To construct rotating configurations the Hartle approach is involved. It is shown that the effects of finite temperatures become crucial in low-mass white dwarfs, whereas rotation is relevant in all mass range. The simultaneous accounting for temperature and rotation is critical in the calculation of the radii of white dwarfs. The results obtained in this work can be applied to explain a variety of observational data for white dwarfs from the Sloan Digital Sky Survey Data Releases.     

Population synthesis of rapidly rotating main-sequence stars with companions

Usingthe “Scenario Machine” (a specialized numerical code formodeling the evolution of large ensembles of binary systems), we have studied the physical properties of rapidly rotating main-sequence binary stars (Be stars) with white-dwarf companions and their abundance in the Galaxy. The calculations are the first to take into account the cooling of the compact object and the effect of synchronization of the rotation on the evolution of Be stars in close binaries. The synchronization time scale can be shorter than the main-sequence lifetime of a Be star formed during the first mass transfer. This strongly influences the distribution of orbital periods for binary Be stars. In particular, it can explain the observed deficit of short-period Be binaries. According to our computations, the number of binary systems in the Galaxy containing a Be star and white dwarf is large: 70–80% of all Be stars in binaries should have degenerate dwarf companions. Based on our calculations, we conclude that the compact components in these systems have high surface temperatures. Despite their high surface temperatures, the detection of white dwarfs in such systems is hampered by the fact that the entire orbit of the white dwarf is embedded in the dense circumstellar envelope of the primary, and all the extreme-UV and soft X-ray emission of the compact object is absorbed by the Be star’s envelope. It may be possible to detect the white dwarfs via observations of helium emission lines of Be stars of not very early spectral types. The ultraviolet continuum energies of these stars are not sufficient to produce helium line emission. We also discuss numerical results for Be stars with other evolved companions, such as helium stars and neutron stars, and suggest an explanation for the absence of Be-black-hole binaries.     

Formation of supermassive black holes due to the tidal deceleration of stellar black holes, globular clusters, and galaxies

The formation and evolution of supermassive (10²?10¹⁰ M _⊙) black holes (SMBHs) in the dense cores of globular clusters and galaxies is investigated. The raw material for the construction of the SMBHs is stellar black holes produced during the evolution of massive (25?150M _⊙) stars. The first SMBHs, with masses of ～1000M _⊙, arise in the centers of the densest and most massive globular clusters. Current scenarios for the formation of SMBHs in the cores of globular clusters are analyzed. The dynamical deceleration of the most massive and slowly moving stellar-mass (< 100M _⊙) black holes, accompanied by the radiation of gravitational waves in late stages, is a probable scenario for the formation of SMBHs in the most massive and densest globular clusters. The dynamical friction of the most massive globular clusters close to the dense cores of their galaxies, with the formation of close binary black holes due to the radiation of gravitational waves, leads to the formation of SMBHs with masses ? 10³ M _⊙ in these regions. The stars of these galaxies form galactic bulges, providing a possible explanation for the correlation between the masses of the bulge and of the central SMBHs. The deceleration of the most massive galaxies in the central regions of the most massive and dense clusters of galaxies could lead to the appearance of the most massive (to 10¹⁰ M _⊙) SMBHs in the cores of cD galaxies. A side product of this cascade scenario for the formation of massive galaxies with SMBHs in their cores is the appearance of stars with high spatial velocities (> 300 km/s). The velocities of neutron stars and stellar-mass black holes can reach ～10⁵ km/s.     

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.     

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

Variability of the photospheric radiation of active K-M dwarfs and their X-ray luminosities

Variability of the photospheric radiation of 40 (dKe-dMe) dwarfs in the solar neighborhood due to variations in the spottedness of their surfaces is analyzed based on the behavior of their mean annual brightnesses over long time intervals. The amplitudes and characteristic time scales of the variations of the mean annual brightness are taken to be indicators of photospheric activity and were used to infer the levels of photospheric activity in the stars studied. The influence of axial rotation on the development of cyclic activity in young red dwarfs and F-M main-sequence stars is analyzed. The durations and amplitudes of the photospheric variability of rapidly rotating (dK0e-dK5e) stars testifies to a higher level of photospheric activity among red dwarfs and solar-type stars. The X-ray luminosities of these stars grow with the amplitude of the variations of the mean annual brightness. However, this is not typical of rapidly rotating M dwarfs, for which the X-ray emission varies by more than two orders of magnitude, although their degrees of spottedness are all virtually the same. A linear relationship between the X-ray and bolometric luminosities is observed for young (dKe-dMe) stars, with their ratios log(L _x/L _bol) being about ?3. These properties can be used to determine whether a red dwarf is a young star or is already on the main sequence.     

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.     

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.

     

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.     

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 quiescent X-ray emission of X-ray novae and the coronal emission of late-type stars

The X-ray luminosities and spectra of F-M stars of luminosity classes IV–V are analyzed. In dwarfs with rotational velocities of about 100 km/s, such as the optical components of low-mass X-ray novae with black holes, hot plasma can be confined in coronal loops even in the presence of fairly weak magnetic fields. Thus, the soft X-ray emission of such systems in their quiescent state (to 10³¹ erg/s) could be associated with the coronal emission of the optical component/dwarf. Two systems studied with subgiants (V1033 Sco and V404 Cyg) have X-ray luminosities 2×10³²–2×10³³ erg/s. The X-ray emission of a solar-type corona cannot provide such luminosities. However, a transition to a non-solar corona is possible in rapidly rotating subgiants—a dynamical corona whose X-ray emission can be one to two orders of magnitude higher than observed for more slowly rotating late-type subgiants in the solar neighborhood. This suggests that the quiescent X-ray emission of these two systems is provided by emission from the corona of the subgiant optical component.     

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.     

Formation of planetary systems and brown dwarfs around single stars

The conditions for the formation of planets and brown dwarfs around single main-sequence stars are considered in two scenarios. The formation of planets and brown dwarfs requires that the initial specific angular momentum of a solar-mass protostar be (0.32)×10¹⁸ cm²/s. The accreted matter of the protostar envelope forms a compact gas ring (disk) around the young star. If the viscosity of the matter in this ring (disk) is small, increasing its mass above a certain limit results in gravitational instability and the formation of a brown dwarf. If the viscosity of the gas is sufficiently large, the bulk of the protostar envelope material will be accreted by the young star, and the gas disk will grow considerably to the size of a protoplanetary dust disk due to the conservation of angular momentum. The formation of dust in the cool part of the extended disk and its subsequent collisional coalescence ultimately results in the formation of solar-type planetary systems.     

Spots and activity of Pleiades stars from observations with the Kepler Space Telescope (K2)

Observations of the K2 continuation of Kepler Space Telescope program are used to estimate the spot coverage S (the fractional spotted area on the surface of an active star) for stars of the Pleiades cluster. The analysis is based on data on photometric variations of 759 confirmed clustermembers, together with their atmospheric parameters, masses, and rotation periods. The relationship between the activity (S) of these Pleiades stars and their effective temperatures shows considerable change in S for stars with temperatures T _eff less than 6100 K (this can be considered the limiting value for which spot formation activity begins) and a monotonic increase in S for cooler objects (a change in the slope for stars with Teff ~ 3700 K). The scatter in this parameter ΔS about its mean dependence on the (V ?Ks)0 color index remains approximately the same over the entire (V?K _s )0 range, including cool, fully convective dwarfs. The computated S values do not indicate differences between slowly rotating and rapidly rotating stars with color indices 1.1 < (V?K _s )0 < 3.7. The main results of this study include measurements of the activity of a large number of stars having the same age (759 members of the Pleiades cluster), resulting in the first determination of the relationship between the spot-forming activity and masses of stars. For 27 stars with masses differing from the solarmass by nomore than 0.1M⊙, themean spot coverage is S = 0.031±0.003, suggesting that the activity of candidate young Suns is more pronounced than that of the present-day Sun. These stars rotate considerably faster than the Sun, with an average rotation period of 4.3^d. The results of this study of cool, low-mass dwarfs of the Pleiades cluster are compared to results from an earlier study of 1570 M stars.