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.     

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 space velocities of radio pulsars

Known models proposed to explain the high space velocities of pulsars based on asymmetry of the transport coefficients of different sorts of neutrinos or electromagnetic radiation can be efficient only in the presence of high magnetic fields (to 10¹⁶ G) or short rotation periods for the neutron stars (of the order of 1 ms). This current study shows that the observed velocities are not correlated with either the pulsar periods or their surface magnetic fields. The initial rotation periods are estimated in a model for the magnetedipolar deceleration of their spin, aßsuming that the pulsar ages are equal to their kinematic ages. The initial period distribution is bimodal, with peaks at 5 ms and 0.5 s, and similar to the current distribution of periods. It is shown that asymmetry of the pulsar electromagnetic radiation is insufficient to give rise to additional acceleration of pulsars during their evolution after the supernova explosion that gave birth to them. The observations testify to deceleration of the motion, most likely due to the influence of the interstellar medium and interactions with nearby objects. The time scale for the exponential decrease in the magnetic field τ_D and in the angle between the rotation axis and magnetic moment τ_ß are estimated, yielding τ_β = 1.4 million years. The derived dependence of the transverse velocity of a pulsar on the angle between the line of sight and the rotation axis of the neutron star corresponds to the expected dependence for acceleration mechanisms associated with asymmetry of the radiation emitted by the two poles of the star.     

The formation of millisecond radio pulsars in close binary systems

An analysis of the basic parameters of a sample of radio and X-ray pulsars that are members of close binary systems is used to separate them into several families according to the nature of the pulsar companions and the previous evolution of the systems. To quantitatively describe the main parameters of close binaries containing neutron stars, we have performed numerical modeling of their evolution. The main driving forces of the evolution of these systems are the nuclear evolution of the donor, the magnetically coupled and radiation-induced stellar winds of the donor, and gravitational-wave radiation. We have considered donors that are low-mass stars in various stages of their evolution, nondegenerate helium stars, and degenerate stars. The systems studied are either the products of the normal evolution of close binaries with large initial component-mass ratios or result from inelastic collisions of old neutron stars with single and binary low-mass, main-sequence stars in the dense cores of globular clusters. The formation of single millisecond pulsars requires either the dynamical disruption of a low-mass (?0.1M_⊙) donor or its complete evaporation under the action of the X-ray radiation of the millisecond pulsar. The observed properties of binary radio pulsars with eccentric orbits combined with the bimodal spatial-velocity distribution of single radio pulsars suggest that it may be possible to explain the observed rotational and spatial motions of all radio pulsars as a result of their formation in close binaries. In this case, neutron stars formed from massive single stars or the components of massive wide binaries probably cannot acquire the high spatial velocities or rapid rotation rates that are required for the birth of a radio pulsar.     

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.     

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 distribution of space velocities of radio pulsars

The distribution of the directions of the space velocities of 67 radio pulsars is shown to be strongly anisotropic. This anisotropy cannot be explained by the structure of our Galaxy or by various types of solar motions. Pulsars with stronger surface magnetic fields B have higher velocities V. The mean value of V for B < 10¹⁰ G is 108 km/s, while 〈V〉 = 340 km/s for B > 10¹⁰ G. These results must be taken into account when identifying a mechanism to explain the observed pulsar velocities and their anisotropy.     

A dynamical study of 12 wide visual binaries

We have carried out the first dynamical studies of the relative motions of the components of 12 wide binaries: ADS 497, ADS 2427 (GL 130.1), ADS 3593, ADS 5436, ADS 6646, ADS 6783, ADS 9559 (δ Boo), ADS 10329, ADS 10759 (ψ Dra), ADS 14878, ADS 15229, and ADS 16558. The analysis is based on series of photographic observations made with the 26-inch refractor of Pulkovo Observatory, supplemented by data from the WDS Catalog, radial velocities and HIPPARCOS parallaxes. We used the parameters of the apparent motions, which can yield the orbit and mass of a binary from observations over a short arc of the orbit of the order of 5°–10°. The orbits and masses for six stars with orbital periods from 500 to 20 000 years have been reliably determined. For five stars with periods from 2000 to 100 000 years, we obtained probable orbits and mass estimates. For ADS 497 and ADS 10329, the mass excesses in their systems are 1 and 3 M _⊙, respectively, relative to the masses expected from the mass-luminosity dependence. In three cases, for ADS 3593, ADS 14878, and ADS 16558, the obtained estimates of the total mass of the components are 46, 25, and 7300 M _⊙, respectively, if the parallaxes are correct. A perturbation with amplitude 0.3″ and period 40 yrs has been detected in the orbital motion of the system ADS 10759, possibly indicating the presence of an invisible satellite with a mass of the order of 0.4 M _⊙. The orientations of the orbits in the Galactic reference frame have been determined for all 12 binaries. The planes of most of the obtained orbits are steeply inclined to the plane of the Galaxy, as was noted in previous studies.     

Strmömgren zones of type II supernovae

The emission measures EM in the directions of supernova remnants and pulsars are considered as functions of their ages t. The resulting plot has a well-defined lower boundary, which can be approximated by the expression EM_min∝1/t. The quantity EM_min increases with decreasing age t and does not level off or reach a maximum until t?500 yr. It is concluded that the bulk of the radiative energy that goes into ionizing and heating the interstellar gas is released at early stages of the supernova remnant’s evolution. We suggest that most of the kinetic energy of the supernova shell is converted into thermal energy and radiated at remnant ages t<100 yr, when the supernova shell, which is expanding at an enormous speed (about 10⁴ km/s), overtakes the shell produced by the presupernova in the supergiant stage. We have estimated the ionization energy E?10⁵¹ erg, diameter L?60 pc, and electron density N_e?7 cm^?3 of the HII regions around the supernovae (the supernova Strömgren zones). A list of objects that can be reliably identified as Strömgren zones of type II supernovae is presented. The plot of pulsar pulse broadening τ as a function of the pulsar age t also has a well-defined lower boundary, for which τ∝t^?2 when t≥1000 yr. This suggests that turbulence develops during the first thousand years after the supernova outburst. It is also concluded that turbulence plays an important role in the formation and evolution of the Strömgren zones of type II supernovae.     

Evolutionary connection between low-mass detached binaries,contact W UMa systems,and blue stragglers

We analyze the distribution of close binary stars in the orbital semimajor axis—primary mass plane. The reduced spatial density of stars with semimajor axes below 10R_⊙ is confirmed. We identify the area in this plane occupied by precursors of W UMa stars, assuming that the driving force causing the components to approach each other is their magnetic stellar wind. This scenario enables us to estimate the rate of formation (0.02/year) and lifetime (10⁸ yr) of W UMa stars. We derive a theoretical estimate of the ratio of the number of blue stragglers, N _BS , and of horizontal-branch stars, N _HB , in globular clusters based on the hypothesis that all blue stragglers are the result of component mergers in W UMa contact binaries. This ratio is N _BS /N _HB =0.4, close to the observed value for 62 Galactic globular clusters. We discuss possible reasons for the considerable dispersion of the observed estimates of this ratio for different clusters in our Galaxy.     

The H<Subscript>2</Subscript>O supermaser region in Orion KL: Epoch 1982.9

Radio interferometric observations of an H₂O maser flare in the Orion Nebula at epoch 1982.9 have been used to determine the flare’s spatial structure. Antennas in the Crimea, Effelsberg, and Onsala were used. The emission region consists of three groups of components. The angular sizes of the components are 0.2–0.9 mas, and the widths of the emitted lines are 0.2–0.7 km/s. The velocities of the components are correlated with their relative positions, which correspond to expanding concentric rings. Assuming a 1 M_⊙ protostar in a Keplerian approximation, the radius of the inner ring R is 15 AU, the velocity of its rotation V_rot is 8.98 km/s, and the radial component of the velocity V_rad is 1.79 km/s. For the outer ring, R=15.7 AU, V_rot=8.79 km/s, and V_rad=2.61 km/s.     

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.     

Variations in the orbital periods of the Algol-type eclipsing binaries RZ Cas and Z Dra

A detailed study of variations of the orbital periods of the Algol-type eclipsing binary systems RZ Cas and Z Dra is presented. The fairly complex variations of the periods of both systems can be represented as a superposition of a secular increase of the period, slow periodic fluctuations, and quasiperiodic oscillations with a small amplitude occurring on timescales of decades. The secular increase of the period can be explained by the steady mass transfer from the less massive to the more massive component with conservation of the total angular momentum. The mass-transfer rate is 5.7 × 10^?9M_⊙/yr for RZ Cas and 3.0×10^?8M_⊙/yr for Z Dra. To explain the long-period cyclic variations of the orbital periods of RZCas and Z Dra, it must be assumed that the eclipsing binaries move in long-period orbits. RZ Cas moves with a period of 133 yr around a third body withmass M₃ > 0.55M_⊙, while Z Dra moves with a period of 60 yr around a third body with mass M₃ > 0.7M_⊙. The residual fluctuations of the periods may be due to a superposition of variations due to magnetic cycles and non-stationary ejections of matter.     

The impact of reflection effects on the parameters of the old pre-cataclysmic variables MS Peg and LM Com

We have determined the main parameters of the old precataclysmic variable stars MS Peg and LM Com. The radial velocities of the components, reflection effects in the spectra, and light curves of the systems are studied based on model stellar atmospheres subject to external irradiation. Forty-seven moderate-resolution spectra for MS Peg and 57 for LM Com obtained with the 6-m telescope of the Special Astrophysical Observatory are used to derive the refined orbital periods of 0.1736660 days and 0.2586873 days, respectively; the orbital eccentricities do not exceed e=0.04. The mass (M _w =0.49e_⊙) and radius (e _w =0.015R_⊙) of the MS Peg primary calculated using the gravitational redshift correspond to those for a cooling carbon white dwarf with a thin hydrogen envelope. The parameters of the red dwarf (M _r =0.19M_⊙, T_eff=3560 K, R _r =0.18R_⊙) are close to those derived from evolutionary tracks for main-sequence M stars with solar chemical composition. The radius (R _r =0.22R_⊙) and temperature (T_eff=3650 K) of the LM Com secondary exceed theoretical estimates for main-sequence stars with masses of M _r =0.17M_⊙. The luminosity excess of the red dwarf in LM Com can be explained by a prolonged (T>5×10⁶ yrs) relaxation of the M star to its normal state after the binary leaves the common-envelope stage. For both systems, theoretical U, B, V, and R light curves and spectra calculated using the adopted sets of parameters are generally consistent with the observations. This confirms the radiative origin of the hot spots, the unimportance of horizontal radiative transport, and the absence of large-scale velocity fields with high values (V_trans>50 km/s) at the surfaces of the secondaries. Most of the emission lines in the spectra of these objects are formed under conditions close to thermalization, enabling modeling of their pro files in an LTE approximation. A strong λ3905 Å emission line has been identified as the 3s²3p4s 1P⁰-3s²3p² 1S SiI λ3905.52 Å line formed in the atmosphere of the hot spot. The observed intensity can be explained by non-LTE “superionization” of SiI atoms by soft UV radiation from the white dwarf. We suggest a technique for identifying binaries whose cool components are subject to UV irradiation based on observations of λ3905 Å emission in their spectra.     

Dynamical evolution of multiple stars

We have modeled the dynamical evolution of small groups of N=3–18 stars in the framework of the gravitational N-body problem, taking into account possible coalescences of stars and the ejection of single and binary stars from the system. The distribution of states is analyzed for a time equal to 300 initial crossing times of the system. The parameters of the binaries and stable triple systems formed, as well as those of ejected single stars, are studied. In most cases, the evolution of the group results in the formation of a binary or stable triple system. The orbital eccentricities of the binaries formed are distributed according to the law f(e)=2e. As a rule, stable triple systems display pronounced hierarchy (the mean ratio of the semimajor axes of the outer and inner binaries is about 20:1). Stars are ejected with velocities from several km/s to several tens of km/s. The results of the modeling are compared with the parameters of observed wide binaries and triple systems.     

The probable quadruple systems FK Aql and FZ Del

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

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.     

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.     

Physical and dynamical parameters of the multiple system HD 222326

We present the results of our study of the physical and dynamical parameters of the multiple system HD 222326. A new method for determining the individual radial velocities of components in wide binary and multiple systems in the case of small radial-velocity differences (δV _r ≤ the FWHMfor the line profiles) is suggested and tested for both model systems and the binary HD 10009. This testing yielded the component radial velocities V _{r 1,2} for HD 10009, enabling us to derive the center-of mass velocity, V _γ, for the first time. We determined the radial velocities of the components of HD 222326 from high-resolution spectra, and refined the orbital parameters of the subsystems using speckle-interferometric observations. A combined spectroscopic and speckle interferometric analysis enabled us to find the positions of the components in the spectral type—luminosity diagram and to estimate their masses. It is likely that the components are all in various evolutionary stages after leaving the main sequence. We analyzed the dynamical evolution of the system using numerical modeling in the gravitational three-body problem and the known stability criteria for triple systems. The system is probably stable on time scales of at least 10⁶ years. The presence of a fourth component in the system is also suggested.     

Magnetic fields and the $$\dot P - P$$ diagram for radio pulsars

Various mechanisms for the loss of angular momentum of neutron stars are analyzed. Theoretical predictions about the evolution of the period are compared with the observed distribution of pulsars on the log\(\dot P\)log(P) diagram. Pulsars with short periods (P≤0.1 s) cannot be fit well by any of the models considered. Their braking index is n=?1, which requires the development of a new braking mechanism. The evolution of pulsars with P>1.25 s is described by the law \(\dot P \propto P^2\), probably due to processes internal to the neutron stars. The observational data for pulsars with 0.1<P≤1.25 s can be fit with a hybrid model incorporating internal processes and magnetic-dipole losses. The magnetic fields in pulsar catalogs should be recomputed in accordance with the results obtained. For example, the magnetic fields obtained for two magnetars with P=5.16 s and P=7.47 s are B _s =1.7×10¹³ and 2.9×10¹³ G, which are lower than the critical field B_cr=4.4×10¹³ G. For a substantial fraction of pulsars, their characteristic ages \(\tau = P/2\dot P\) cannot serve as measures of their real ages.