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.     

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.     

Evolution of the masses of neutron stars in binary systems

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

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

Magnetic field models for HD 116458 and HD 126515

We have modeled the magnetic fields of the slowly rotating stars HD 116458 and HD 126515 using the “magnetic charge” technique. HD 116458 has a small angle between its rotation axis and dipole axis (β = 12°), whereas this angle is large for HD 126515 (β = 86°). Both stars can be described with a decentered-dipole model, with the respective displacements being r = 0.07 and r = 0.24 in units of the stellar radius. The decentered-dipole model is able to satisfactorily explain the phase relations for the effective field, B_e(P), and the mean surface field, B_s(P), for both stars, along with the fact that the B_e(P) phase relation for HD 126515 is anharmonic. We discuss the role of systematic measurement errors possibly resulting from instrumental or methodical effects in one or both of the phase relations. The displacement of the dipole probably reflects real asymmetry of the stellar field structure, and is not due to measurement errors. Using both phase relations, B_e(P) and B_s(P), in the modeling considerably reduces the influence of the nonuniform distribution of chemical elements on the stellar surface.     

A spindown mechanism for short-period radio pulsars

It is shown that a model with accretion in a “quasi-propeller” mode can explain the observed spindown of pulsars with periods P<0.1 s. The mean accretion rate for 39 selected objects is \(\dot M = 5.6 \times 10^{ - 11} M_ \odot /year\). If \(\dot M\) is constant during the pulsar’s lifetime, the neutron star will stop rotating after 10⁷ years. The mean magnetic field at the neutron-star surface calculated in this model, \(\bar H_0 = 6.8 \times 10^8 G\), is consistent to an order of magnitude with the values of H₀ for millisecond pulsars from known catalogs. However, the actual value of H₀ for particular objects can differ from the catalog values by appreciable factors, and these quantities must be recalculated using more adequate models. The accretion disk around the neutron star should not impede the escape of the pulsar’s radiation, since this radiation is generated near the light cylinder in pulsars with P<0.1 s. Pulsars such as PSR 0531+21 and PSR 0833-45 have probably spun down due to the effect of magnetic-dipole radiation. If the difference in the braking indices for these objects from n=3 is due to the effect of accretion, the accretion rate must be of the order of 10¹⁸ g/s.     

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

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

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.     

Gravitational darkening coefficients of the stars in the semidetached binary V Puppis

We have analyzed high-precision vby light curvesfor the semi-detached binary V Pup in a Roche model. They are consistent with the standard gravitational darkening coefficient for hot stars, β = 0.25, rather than the value β = 1.36 ± 0.04 derived by Kitamura and Nakamura [1] using a simpler model. We rigorously estimate the confidence intervals for the allowed gravitational darkening coefficients for a star filling its Roche lobe to be β = (?0.24, +1.29) for the 99% confidence level and β = (?0.21, +1.26) for the 67% confidence level.     

Modeling the wind of the Be star SS 2883

Observations of eclipses of the radio pulsar B1259-63 by the disk of its Be-star companion SS 2883 provide an excellent opportunity to study the winds of stars of this type. The eclipses lead to variations in the radio flux (due to variations in the free-free absorption), dispersion measure, rotation measure, and linear polarization of the pulsar. We have carried out numerical modeling of the parameters of the Be-star wind and compared the results with observations. The analysis assumes that the Be-star wind has two components: a disk wind in the equatorial plane of the Be star with a power-law fall-off in the electron density n _e with distance from the center of the star \(\rho (n_e \sim \rho ^{ - \beta _o } )\), and a spherical wind above the poles. The parameters for a disk model of the wind are estimated. The disk is thin (opening angle 7.5°) and dense (electron density at the stellar surface n_0e ～ 10¹² cm^?3, β₀ = 2.55). The spherical wind is weak (n_0e ? 10⁹ cm^?3, β₀ = 2). This is the first comparison of calculated and observed fluxes of the pulsating radio emission.     

Distribution of inhomogeneities in the interstellar plasma in the directions of three distant pulsars from observations with the RadioAstron ground–space interferometer

The RadioAstron ground–space interferometer has been used to measure the angular sizes of the scattering disks of the three distant pulsars B1641–45, B1749–28, and B1933+16. The observations were carried out with the participation of the Westerbork Synthesis Radio Telescope; two 32-m telescopes at Torun, Poland and Svetloe, Russia (the latter being one antenna of the KVAZAR network); the Saint Croix VLBA antenna; the Arecibo radio telescope; the Parkes, Narrabri (ATCA), Mopra, Hobart, and Ceduna Australian radio telescopes; and the Hartebeesthoek radio telescope in South Africa. The full widths at half maximum of the scattering disks were 27 mas at 1668 MHz for B1641–45, 0.5 mas at 1668 MHz for B1749–28, and 12.3 at 316 MHz and 0.84 mas at 1668 MHz for B1933+16. The characteristic time scales for scatter-broadening of the pulses on inhomogeneities in the interstellar plasma τ_sc were also measured for these pulsars using various methods. Joint knowledge of the size of the scattering disk and the scatter-broadening time scale enables estimation of the distance to the effective scattering screen d. For B1641–45, d = 3.0 kpc for a distance to the pulsar D = 4.9 kpc, and for B1749–28, d = 0.95 kpc for D = 1.3 kpc. Observations of B1933+16 were carried out simultaneously at 316 and 1668 MHz. The positions of the screen derived using the measurements at the two frequencies agree: d ₁ = 2.6 and d ₂ = 2.7 kpc, for a distance to the pulsar of 3.7 kpc. Two screens were detected for this pulsar from an analysis of parabolic arcs in the secondary dynamic spectrum at 1668 MHz, at 1.3 and 3.1 kpc. The scattering screens for two of the pulsars are identified with real physical objects located along the lines of sight toward the pulsars: G339.1–04 (B1641–45) and G0.55–0.85 (B1749–28).     

Peculiarities of radio pulsars with emission outside the radio range

Parameters of 100 radio pulsars detected outside the radio range (he pulsars) are compared with those of pulsars radiating only in the radio (n pulsars). The periods of he pulsars are, on average, appreciably shorter than those of n pulsars: 〈P〉 = 0.10 and 0.56 s, respectively. The distribution of the magnetic field at the light cylinder is shifted toward higher magnetic fields for the pulsars with high-energy radiation, compared to the distribution for pulsars radiating only in the radio. The magnetic fields at the light cylinder are 〈B _lc〉 = 9×10³ G for he radio pulsars, and 〈B_lc〉 = 56 G formost purely radio pulsars. This suggests the generation of high-energy nonthermal radiation in radio pulsars at the peripheries of their magnetospheres. The distribution of the spin-energy loss rate dE/dt is uniform for he pulsars, and is characterized by a higher average value \(\left( {\left\langle {\log \frac{{dE}} {{dt}}} \right\rangle = 35.53} \right) \) , compared to n pulsars, \(\left( {\left\langle {\log \frac{{dE}} {{dt}}} \right\rangle = 32.60} \right) \) . The spatial distribution of he pulsars is nonuniform: they form two well separated clouds.     

A new model of a magnetar

A new model is put forward to explain the observed features of anomalous X-ray pulsars (AXPs) and soft gamma-ray repeaters (SGRs). It is shown that drift waves can be excited in the magnetosphere of a neutron star with a rotational period of P～0.1 s, surface magnetic field B_s～10¹² G, and angle between the rotational axis and magnetic moment β<10°. These waves lead to the formation of radiation pulses with a period of P_dr～10 s. The rate of loss of rotational energy by such a star (～10³⁷ erg/s) is sufficient to produce the observed increase in the period \((\dot P \sim 10^{ - 10} )\), the X-ray luminosities of AXPs and SGRs (～10³⁴–10³⁶ erg/s), and an injection of relativistic particles into the surrounding supernova remnant. A modulation of the constant component of the radiation with a period of P～0.1 s is predicted. In order for SGRs to produce gamma-ray bursts, an additional source of energy must be invoked. Radio pulsars with periods of P_obs>5 s can be described by the proposed model; in this case, their rotational periods are considerably less than P_obs and the observed pulses are due to the drift waves.     

The pulsar J1852+0040 in Kes 79 and the nature of central compact objects in supernova remnants

A possible model for the pulsar PSR J1852+0040 associated with the supernova remnant Kes 79 and detected in place of a central compact object in this remnant is discussed. The main observational properties of the pulsar can be understood as consequences of its weak surface magnetic field (B _s < 3 × 10¹¹ G) and short rotational period (P ～ 0.1 s). Its X-ray emission is thermal, and is generated in a small region near the surface of the neutron star due to cooling of the surface as the surface accretes matter from a relict disk surrounding the pulsar. The radio emission is generated in the outer layers of the pulsar magnetosphere by the synchrotron (cyclotron) mechanism. The optical luminosity of J1852+0040 is estimated to be L _opt < 10²⁸ erg/s. If the spectral features in another central compact object, 1E 1207.4+5209, are interpreted as electron cyclotron lines, this provides evidence for a weak surface magnetic field for this neutron star as well (B < 6 × 10¹⁰ G). The hypothesis that all central compact objects have weak surface fields makes it possible to explain the number of detected central compact objects, the absence of pulsar-wind nebulae associated with these objects, and the fact that no pulsar has yet been detected at the position of SN 1987a. We suggest that, after the supernova remnant has dissipated, the central compact object becomes a weak X-ray source (XDINS), whose weak emission is also due to the weakness of its magnetic field.     

The pulsar PSR B1931+24 as an orthogonal rotator

The angles of the magnetic moment μ and the line of sight L to the rotation axis Ω are estimated for the pulsar PSR B1921+24, which displays “on” and “off” periods in its radio emission. It is shown that this object is an orthogonal rotator, i.e., the angle β between μ and Ω is equal to 88°.2 and the angle between L and Ω is ζ = 98.7°, and that its rotation period should be twice the usually adopted value (P = 1.626 s). One possible reason for the peculiarities of this pulsar could be the precession of a relic disk in the equatorial region of the object. Further observations (in particular, in the infrared) are required to confirm the existence of such a disk. Polarization data for other pulsars whose radiation switches on and off (transients) are also required, to determine if they are likewise orthogonal rotators. Calculations for PSR B0656+14 show that β ∼ 20°, and the sharp increase of its pulse intensities is due to intrinsic reasons, and is not associated with a relic disk. Original Russian Text ? I.F. Malov, 2007, published in Astronomicheskiĭ Zhurnal, 2007, Vol. 84, No. 6, pp. 531–535.     

Spots and the Activity of Stars in the Hyades Cluster from Observations with the Kepler Space Telescope (K2)

Observations of the K2 mission (continuing the program of the Kepler Space Telescope) are used to estimate the spot coverage S (the fractional area of spots on the surface of an active star) for stars of the Hyades cluster. The analysis is based on data on the photometric variations of 47 confirmed single cluster members, together with their atmospheric parameters, masses, and rotation periods. The resulting values of S for these Hyades objects are lower than those stars of the Pleiades cluster (on average, by ΔS ~ 0.05?0.06). A comparison of the results of studies of cool, low-mass dwarfs in the Hyades and Pleiades clusters, as well as the results of a study of 1570 M stars from the main field observed in the Kepler SpaceMission, indicates that the Hyades stars are more evolved than the Pleiades stars, and demonstrate lower activity. The activity of seven solar-type Hyades stars (S = 0.013 ± 0.006) almost approaches the activity level of the present-day Sun, and is lower than the activity of solar-mass stars in the Pleiades (S = 0.031 ± 0.003). Solar-type stars in the Hyades rotate faster than the Sun (〈P〉 = 8.6 ^d ), but slower than similar Pleiades stars.     

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.