Classification of radio pulsars using the principle-components method

The principle-components method is used as a basis to analyze the distributions of known radio pulsars in spaces of eigenvectors of correlation matrices for various samples of pulsars and classification parameters (from 4 to 11 parameters characterizing the physical and kinematic properties of the objects). Pulsars with periods P < 0.1 s form a separate cluster, far from the cluster formed by “normal” pulsars with P ～ 1 s, in all the studied spaces. These two groups also differ appreciably in their other parameters (period derivatives, magnetic fields, pulse widths). In particular, the spatial velocities of short-period pulsars (106 km/s) are appreciably lower than those displayed by long-period pulsars (334 km/s). The distributions of the pulsars at southern (Z < 0) and northern (Z > 0) Galactic latitudes do not differ; i.e., there is no anisotropy in the motions in these two directions perpendicular to the Galactic plane, or in the corresponding distributions of the pulsar parameters.     

Differences in the parameters of radio pulsars with short and long periods

A comparative analysis of various parameters of pulsars with short (P < 0.1 s) and long (P > 0.1 s) periods is carried out. There is no correlation between the radio and gamma-ray luminosities of the pulsars and their surfacemagnetic fields, but there is a correlation between the X-ray luminosity and the surfacemagnetic field. A dependence of the X-ray and gamma-ray luminosities on the magnetic field at the light cylinder is also found. This result provides evidence for the formation of hard, non-thermal emission at the periphery of the magnetosphere. An appreciable positive correlation between the luminosity and the rate of rotational energy loss by the neutron star is observed, supporting the idea that all radio pulsars have the same basic source of energy. The efficiency of the transformation of rotational energy into radiation is significantly higher in long-period pulsars. The dependence of the pulse width on the pulsar period is steeper for pulsars with short periods than for those with long periods. The results obtained support earlier assertions that there are differences in the processes generating the emission in pulsars with P < 0.1 s and those with P > 0.1 s.     

The geometry of radio pulsar magnetospheres

Data on the profiles and polarization of the 10- and 20-cm emission of radio pulsars are used to calculate the angle β between the rotational axis of the neutron star and its magnetic moment. It is shown that, for these calculations, it is sufficient to use catalog values of the pulse width at the 10% level W ₁₀, since the broadening of the observed pulses due to the transition to the full width W ₀ and narrowing of the pulses associated with the emission of radiation along tangents to the field lines approximately cancel each other out. The angles β ₁ are calculated for 283 pulsars at 20 cm and 132 pulsars at 10 cm, assuming that the line of sight passes through the center of the emission cone. The mean values of these angles are small and similar for the two wavelengths (〈β ₁〉 = 18° at λ = 10 cm and 〈β ₁〉 = 14° at λ = 20 cm). The angle β ₂ is estimated for several dozen pulsars for the case when the orientation of the angle to the line of sight is arbitrary. The mean value of β ₂ at 10 cm is found to be 〈β ₂〉 = 33.9° if the maximum derivative of the polarization position angle C is positive and 〈β ₂〉 = 52.1° ifC < 0. We find at 20 cm 〈β ₂〉 = 33.9° ifC > 0 and 〈β2〉 = 54.1° ifC < 0. The values at the two wavelengths are equal within the errors, and close to the β ₂ value obtained earlier at 30 cm (〈β ₂〉 = 36.4° if C >0 and 〈β2〉 = 49.1° if C < 0). The mean 〈β ₂〉 for the entire set of data can be taken to be 43.5°. The period dependence of the pulse width W(P) √ P ^−0.25 differs from the relation that is usually used in the polar-cap model, W(P) √ P ^−0.5. This difference could be associated with the rate of development of plasma instabilities near the surface of the neutron star (in the region where high-frequency radiation is generated). The role of the quadrupole component of the magnetic field is not important here. There is no dependence of the angle β on the pulsar age (z distance, luminosity L, or characteristic age τ = P/(2dP/dt)) for the studied sample.     

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.     

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.     

On Braking Mechanisms in Radio Pulsars

Astronomy Reports - The dependence dP/dt(P) is analyzed for three groups of pulsars with different period values: P&nbsp;&gt;&nbsp;2&nbsp;s, 0.1 s &lt; P &lt; 2 s, and P...     

Rotational periods and other parameters of magnetars

The rotational periods P, period derivatives dP/dt, and magnetic fields B in the region where the emission of anomalous X-ray pulsars (AXPs) and soft gamma-ray repeaters (SGRs) is generated are calculated using a model that associates the emission of these objects with the existence of drift waves at the periphery of the magnetosphere of a neutron star. The values obtained for these parameters are P = 11?737 ms, dP/dt = 3.7 × 10^?16?5.5 × 10^?12, and log B (G) = 2.63?6.25. We find a dependence between the X-ray luminosity of AXPs and SGRs, L _x, and the rate at which they lose rotational energy, dE/dt, which is similar to the L _x(dE/dt) dependence for radio pulsars with detected X-ray emission. Within the errors, AXPs/SGRs and radio pulsars with short periods (P < 0.1 s) display the same slopes for their log(dP/dt)-log P relations and for the dependence of the efficiency of their transformation of rotational energy into radiation on their periods. A dipole model is used to calculate the surface magnetic fields of the neutron stars in AXPs and SGRs, which turn out to be, on average, comparable to the surface fields of normal radio pulsars (〈log B _s (G)〉 = 11.90).     

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.     

Detection of regular variations in the intensity and pulse time of arrival of the anomalous pulsar PSR B0943+10

Timing of the anomalous pulsar PSR B0943+10 during 2007–2013 was carried out on the Large Phased Array radio telescope of the Pushchino Radio Astronomy Observatory at 112 MHz. The astrometric and rotational parameters for epoch MJD=56 500 have been determined. Considerable deviations of the pulse times of arrival from the precalculated values with a characteristic period of several years due to the presence of correlated low-frequency noise in the pulsar spin phase have been detected. These deviations can be explained in a planetary model by the presence of two companions of the pulsar, whose orbital parameters have been determined. A continuous increase in the longitude of the pulse maximum within the emission window, the pulse width, and the intensity have been detected after each switch to the burst mode. Together with the changes in pulse shape, degree of linear polarization of the pulse, and drift rate of individual pulses detected earlier, this indicates that all the main parameters of the radio emission in the B mode are unstable. This distinguishes PSR B0943+10 from all other modes-witching pulsars. The origin of the observed properties of this pulsar are probably associated with the interaction of its extended magnetosphere with the surrounding medium.     

Optimal filtration and a pulsar time scale

An algorithm is proposed for constructing a group (ensemble) pulsar time based on the application of optimal Wiener filters. This algorithm makes it possible to separate the contributions of variations of the atomic time scale and of the pulsar rotation to barycentric residual deviations of the pulse arrival times. The method is applied to observations of the pulsars PSR B1855+09 and PSR B1937+21, and is used to obtain corrections to UTC relative to the group pulsar time PT_ens. Direct comparison of the terrestial time TT(BIPM06) and the group pulsar time PT_ens shows that they disagree by no more than 0.4 ± 0.17 μs. Based on the fractional instability of the time difference TT(BIPM06)-PT_ens, σ _z = (0.5 ± 2) × 10⁻¹⁵, a new limit for the energy density of the gravitational-wave background is established at the level Ω _g h ² ∼ 10⁻⁹.     

The drift model for AXPs and SGRs in the light of new observtional data

The latest observational data are analyzed to investigate their consistency with two known models for anomalous X-ray pulsars (AXPs) and soft gamma-ray repeaters (SGRs): the magnetar and drift models. The results of spectral measurements disagree with the predictions of theories that assume the presence of super-strong magnetic fields on AXPs and SGRs and associated processes for the generation of electron-positron plasma in the upper layers of the neutron-star magnetosphere. We present arguments against the use of magnetic-dipole braking for these objects. The rotational periods P, period derivatives dP/dt, and magnetic fields B of known AXPs and SGRs are calculated for the drift model. The mean values of these parameters in the sample used are 〈P〉 = 108 ms and 〈log B _s [G]〉 = 12.08. Overall, the measured profiles, polarizations, and spectra can be brought into agreement with the drift model.     

Detection of gravitational waves through observations of a group of pulsars

We suggest a new approach to the detection of gravitational waves using observations of a group of millisecond pulsars. In contrast to the usual method, based on increasing the accuracy of the arrival times of pulses by excluding possible distorting factors, our method supposes that the additive phase noise that is inevitably present even in the most accurate observational data has various spectral components, which have characteristic amplitudes and begin to appear on different time scales. We use the “Caterpillar” (Singular Spectral Analysis, SSA) method to decompose the signal into its components. Our initial data are the residuals of the pulse arrival times for six millisecond pulsars. We constructed the angular correlation function for components of the decomposition of a given number, whose theoretical form for the case of an isotropic and homogeneous gravitational-wave background is known. The individual decomposition components show a statistically significant agreement with the theoretical expectations (correlation coefficient ρ = 0.92 ± 0.10).     

The magnetospheric structure of radio pulsars with interpulses

Pulsars with interpulses—pulse components located between the main pulses—are studied. About 50 such objects are currently known. Methods developed earlier to determine the angle β between the rotation axis and the magnetic moment of the neutron star are used to investigate the geometry of the magnetospheres in these objects. In a number of pulsars, β < 20°, so that not only interpulses, but also radiation between pulses and a correlation between the behaviors of the interpulses and main pulses, is expected. In other pulses, this angle is greater than 60°, and interpulses can appear if the radiation cone is sufficiently broad and there is a favorable orientation of the line of sight of the observer. Thus, the earlier prediction that there should be two types of pulsars with interpulses—aligned and orthogonal—is supported. Estimates of the ages of the pulsars in these two groups indicate that aligned rotators are appreciably older than orthogonal rotators.     

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.     

Angles between the rotational axis and magnetic moment in 80 pulsars from observations near 1 GHz

Data on the pulse structure and variations of the linear polarization angle at frequencies near 1 GHz have been used to estimate the angles β between the rotational axis and magnetic moment of the neutron stars assocaited with 80 pulsars. The calculations applied several methods. The minimum values of β were estimated from the observed pulse width W ₁₀ at the 10% level for the entire sample. Maximum estimates of β were obtained for six sources with small polarization position angle derivatives. Equations for the angle β were derived for various forms of the observed profile, and solutions obtained for 34 pulsars. The β values calculated using different methods are compared. For three pulsars with known interpulses, the obtained values of β demonstrate that two (PSR B1055-52 and PSR 1822-09) are aligned rotators, whereas the other (PSR B1702-19) is an orthogonal rotator. A search for interpulses and interpulse emission in PSRB1641-45, PSR1642-03, and PSR 1944+17 is necessary, and a search for an interpulse at 180° from the main pulse is required in PSR B2321-61.     

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.     

Pulsed optical emission by radio pulsars

The luminosity L of radio pulsars due to synchrotron radiation by the primary beam at the magnetosphere periphery is derived. There is a strong correlation between the observed optical luminosities of radio pulsars and the parameter $\dot P/P^4$ (where P is the pulsar period). This correlation predicts appreciable optical emission from several dozen pulsars, in particular, from all those with P<0.1 s. Agreement with optical observations can be achieved for Lorentz factors of the secondary plasma γ_p=2–13. Plasma with such energies can be produced only when the magnetic-field structure near the neutron-star surface deviates substantially from a dipolar field. The peak frequency of the synchrotron spectrum should shift toward higher values as the pulsar period P decreases; this is, in agreement with observational data for 27 radio pulsars for which emission has been detected outside the radio band.     

Search for and detection of pulsars inmonitoring observations at 111 MHz

In the course of monitoring interplanetary scintillations of a large number of sources using the Big Scanning Antenna of the Lebedev Physical Institute, a search for pulsars with periods ≥0.4 s at declinations ?9? < δ < 42? and right ascensions 0h < α < 24h was simultaneously carried out. The search was conducted using four years of observations carried out at 110.25MHz in six frequency channels making up a 2.5 MHz band and having a time resolution of 100 ms. The initial identification of pulsar candidates was done using Fourier power spectra averaged over the entire observational period; the pulsar candidates were then verified using observations with higher frequency and time resolution: 32 frequency channels and a time resolution of 12.5 ms. Eighteen new pulsars were discovered in the studied area, whose main characteristics are presented.     

The radio emission of anomalous pulsars

Previously developed methods for estimating the angle β between the spin axis of a neutron star and its magnetic moment together with observational data for anomalous X-ray pulsars (AXPs) indicate that these objects are nearly aligned rotators, and that the drift model can be applied to them. The magnetospheres of aligned rotators are appreciably more extended than in pulsars with large values of β. With such extents for the magnetosphere, the conditions for the generation of transverse waves via the cyclotron instability are satisfied. The expected spectrum of the resulting radiation is very steep (its spectral index is α > 3), consistent with the observed radio spectra of known AXPs (α > 2). A large magnetosphere also favors the appearance of appreciable pitch angles for relativistic electrons, and therefore the generation of synchrotron emission. The maximum of this emission falls in the microwave range. This mechanism provides appreciable fluxes at frequencies of tens of gigahertz and can explain the observed enhanced AXP radiation in this range.