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.     

Parameters of radio pulsars in binary systems and globular clusters

The parameters of radio pulsars in binary systems and globular clusters are investigated. It is shown that such pulsars tend to have short periods (of the order of several milliseconds). Themagnetic fields of most of the pulsars considered are weak (surface fields of the order of 10⁸?10⁹ G). This corresponds to the generally accepted view that short-period neutron stars are spun up by angular momentum associated with the stellar wind from a companion. However, the fields at the light cylinders in these objects are two to three orders of magnitude higher than for the main population of single neutron stars. The dependence of the pulse width on the period does not differ from the corresponding dependences for single pulsars, assuming the emission is generated inside the polar cap, at moderate distances from the surface or near the light cylinder. The radio luminosities of pulsars in binary systems do not show the correlation with the rate of loss of rotational energy that is characteristic for single pulsars, probably due to the influence of accreting matter from a companion. Moreover, accretion apparently decreases the power of the emergent radiation, and can explain the observed systematic excess of the radio luminosity of single pulsars compared to pulsars in binary systems. The distributions and dependences presented in the article support generally accepted concepts concerning the processes occurring in binary systems containing neutron stars.     

Acceleration of close binaries in supernova explosions

We investigate the orientation of the orbital planes of X-ray binary systems relative to the direction of the additional velocity acquired by the binary in a supernova explosion or as a result of radiative acceleration. In the second case, the acceleration occurs due to X-ray radiation during a stage of intense accretion onto the neutron star, which has an asymmetric magnetic field. Observational consequences that could enable estimation of the role of each acceleration mechanism are discussed. The results are also applicable to binary millisecond radio pulsars, assuming that they have gone through an accretion stage.     

Magnetorotational supernova explosions and the formation of neutron stars in close binary systems

The formation of neutron stars in the closest binary systems (P _orb<12 h) gives the young neutron star/pulsar a high rotational velocity and energy. The presence of a magnetic field of 3×10¹¹–3×10¹³ G, as is observed for radio pulsars, enables the neutron star to transfer ～10⁵¹ erg of its rotational energy to the envelope over a time scale of less than an hour, leading to a magnetorotational supernova explosion. Estimates indicate that about 30% of all type-Ib,c supernovae may be the products of magnetorotational explosions. Young pulsars produced by such supernovae should exhibit comparatively slow rotation (P _rot>0.01 s), since a large fraction of their rotational angular momentum is lost during the explosion. The magnetorotational mechanism for the ejection of the envelope is also reflected by the shape of the envelope. It is possible that the Crab radio pulsar is an example of a product of a magnetorotational supernova. A possible scenario for the formation of the close binary radio pulsar discovered recently by Lyne et al. is considered.     

Mass exchange in close binaries on the dynamical time scale during the stage of contracting onto the main sequence

We consider mass exchange in close binaries containing low-mass, fully convective components on the dynamical time scale. We present the results of three-dimensional hydrodynamical simulations of this process in close binaries contracting onto the main sequence (binary pre-MS stars). Our results suggest that some systems with superplanets could be formed by this process, which could be considered a special type of binary star. We have determined the ranges of the relative donor masses that allow the merger of the binary components and the formation of systems with superplanets or the survival of the binary during the mass exchange. These process should result in a deficit of binaries with similar component masses.     

The nature hypervelocity stars

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

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.     

Period clustering of anomalous X-ray pulsars

The question of why the observed periods of anomalous X-ray pulsars (AXPs) and soft gamma-ray repeaters (SGRs) cluster in the range 2–12 s is discussed. The possibility that AXPs and SGRs are the descendants of high-mass X-ray binaries that have disintegrated in core-collapse supernova explosions is investigated. The spin periods of neutron stars in high-mass X-ray binaries evolve towards the equilibrium period, which is a few seconds, on average. After the explosion of its massive companion, the neutron star becomes embedded in a dense gaseous envelope, and accretion from this envelope leads to the formation of a residual magnetically levitating disk. It is shown that the expected mass of the disk in this case is 10^?7–10^?8 M_⊙, which is sufficient to support accretion at the rate 10¹⁴–10¹⁵ g/s over a few thousand years. During this period, the star manifests itself as an isolated X-ray pulsar with a number of parameters similar to those of AXPs and SGRs. The periods of such pulsars can cluster if the lifetime of the residual disk does not exceed the spin-down timescale of the neutron star.     

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.     

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

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

     

Wolf-Rayet stars 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 mass spectrum of black holes in close binary systems

We present the results of population syntheses obtained using our “scenario machine.” The mass spectra of black holes in X-ray binary systems before and after the stage of accretion from an optical companion are obtained for various evolutionary scenarios. The results of the model computations are compared to observational data. The observational data are used to estimate the fraction of a presupernova’s mass that collapses into a black hole. This model can explain the formation of low-mass (2–4M_⊙) black holes in binary systems with optical companions. We show that the number of low-mass black holes in the Galaxy is sufficiently high for them to be detected. The population-synthesis results suggest that the vast majority of low-mass black holes are formed via the accretion-induced collapse of neutron stars. The percentage of low-mass black holes in binary systems that form due to accretion-induced collapse is 2–15% of the total number of black holes in binaries, depending on the evolutionary scenario.     

The evolution of stars paired with supermassive black holes

A star located in the close vicinity of a supermassive black hole (SMBH) in a galactic nucleus or a globular-cluster core could form a close binary with the SMBH, with the star possibly filling its Roche lobe. The evolution of such binary systems is studied assuming that the SMBH mainly accretes matter from the companion star and that the presence of gas in the vicinity of the SMBH does not appreciably influence variations in the star’s orbit. The evolution of the star–SMBH system is mainly determined by the same processes as those determining the evolution of ordinary binaries. The main differences are that the star is subject to an incident flux of hard radiation arising during the accretion of matter by the SMBH, and, in detached systems, the SMBH captures virtually all the wind emitted by its stellar companion, which appreciably influences the evolution of the major axis of the orbit. Moreover, the exchange between the orbital angular momentum and the angular momentum of the overflowing matter may not be entirely standard in such systems. The computations assume that there will be no such exchange of angular momentum if the characteristic timescale for mass transfer is shorter than the thermal time scale of the star. The absorption of external radiation in the stellar envelope was computed using the same formalism applied when computing the opacity of the stellar matter. The numerical simulations show that, with the adopted assumptions, three types of evolution are possible for such a binary system, depending on the masses and the initial separation of the SMBH and star. Type I evolution leads to the complete destruction of the star. Only this type of evolution is realized for low-mass main-sequence (MS) stars, even those with large initial separations from their SMBHs. Massive MS stars will also be destroyed if the initial separation is sufficiently small. However, two other types of evolution are possible for massive stars, with a determining role in the time variations of the parameters of the star–SMBH system being played by the possible growth of the massive star into a red giant during the time it is located in the close vicinity of the SMBH. Type II evolution can be realized for massive MS stars that are initially farther from the SMBH than in the case of disruption. In this case, the massive star fills its Roche lobe during its expansion, but is not fully destroyed; the star retreats inside its Roche lobe after a period of intense mass loss. This type of evolution is characterized by an increase in the orbital period of the system with time. As a result, the remnant of the star (its former core) is preserved as a white dwarf, and can end up at a fairly large distance from the SMBH. Type III evolution can be realized formassiveMSstars that are initially located still farther from their SMBHs, and also for massive stars that are already evolved at the initial time. In these cases, the star moves away from the SMBH without filling its Roche lobe, due to its intense stellar wind. The remnants of such stars can also end up at a fairly large distances from their SMBHs.     

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.     

Radio pulsars in close binaries with neutron stars

The initial distribution of young radio pulsars, reconstructed from the observed distribution of their spatial velocities distorted by observational selection effects, taking into account the age and spatial distribution of radio pulsars with measured spatial velocities, appears to be bimodal. Most young pulsars are formed with velocities of ～100 km/s. Some fraction of young radio pulsars display an almost flat velocity distribution (dv/dv ≈ const) from 150 to 1000 km/s. Scenario modeling in the absence of an additional (kick) velocity acquired by the young neutron star during its formation in a supernova explosion can reproduce the initial velocity distribution of radio pulsars, but results in a higher fraction of radio pulsars in binaries than is observed. Assuming a random initial Maxwellian kick velocity of ～100 km/s makes it possible to reduce the fraction of radio pulsars in binaries to the observed value (<1%), while leaving the velocity distribution for radio pulsars close to the observed bimodal initial distribution. Such kick velocities are also able to explain the observed distribution of radio pulsars in close binaries in the eccentricity—orbital period plane.     

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

Evolution of stars near bright quasars

We consider the evolution of stars that are located in the near vicinity of bright quasars and illuminated by their hard radiation. The absorption of the external radiation flux in the stellar envelope was calculated applying the formalism used to determine the opacity of the stellar material. The simulations show that the external irradiation heats the outer layers of the star, altering their structure and reducing the thickness of the convective envelope in low-mass stars, but leaves the inner layers essentially unaffected. In addition, the irradiation substantially increases the mass loss in the stars. This is important for our understanding of the evolution of quasar masses, since this increased mass loss by nearby stars can supply the quasar with additional accreting gas. The results can also be applied to detached binaries in which a low-mass star is irradiated by a very massive companion.     

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.     

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.