Populating the Galaxy with pulsars – I. Stellar and binary evolution

The computation of theoretical pulsar populations has been a major component of pulsar studies since the 1970s. However, the majority of pulsar population synthesis has only regarded isolated pulsar evolution. Those that have examined pulsar evolution within binary systems tend to either treat binary evolution poorly or evolve the pulsar population in an ad hoc manner. Thus, no complete and direct comparison with observations of the pulsar population within the Galactic disc has been possible to date. Described here is the first component of what will be a complete synthetic pulsar population survey code. This component is used to evolve both isolated and binary pulsars. Synthetic observational surveys can then be performed on this population for a variety of radio telescopes. The final tool used for completing this work will be a code comprised of three components: stellar/binary evolution, Galactic kinematics and survey selection effects. Results provided here support the need for further (apparent) pulsar magnetic field decay during accretion, while they conversely suggest the need for a re-evaluation of the assumed typical millisecond pulsar formation process. Results also focus on reproducing the observed     diagram for Galactic pulsars and how this precludes short time-scales for standard pulsar exponential magnetic field decay. Finally, comparisons of bulk pulsar population characteristics are made to observations displaying the predictive power of this code, while we also show that under standard binary evolutionary assumption binary pulsars may accrete much mass.     

Power spectrum analysis of the timing noise in 18 southern pulsars

Power spectra of the timing noise observed in 18 southern pulsars have been derived using a novel technique, based on the CLEAN algorithm. Most of the spectra are well described by a single- or double-component power-law model. Some of these spectra can be interpreted in the context of one or more of the current timing noise models. The results combined with those obtained from the time-domain analyses of the timing activity in these pulsars are used to assess the viability of the various theoretical models of pulsar timing noise.     

Gamma-ray pulsars: the pulse profiles and phase-resolved spectra of Geminga

We present a calculation of a three-dimensional pulsar magnetosphere model to explain high-energy emission from the Geminga pulsar with a thick outer gap. High-energy γ -rays are produced by primary accelerated particles with a power-law energy distribution through curvature radiation inside the outer gap. We also calculate the emission pattern, pulse profile and phase-resolved spectra of high-energy γ -rays of the Geminga pulsar, and find that its pulse profile is consistent with the observed one if the magnetic inclination and viewing angle are ∼50° and ∼86° respectively. We describe the relative phases among soft (thermal) X-rays, hard (non-thermal) X-rays, and γ -rays. Our results indicate that X-ray and γ -ray emission from the Geminga pulsar may be explained by the single thick outer gap model. Finally, we discuss the implications of the radio and optical emission of the Geminga pulsar.     

The Parkes Southern Pulsar Survey — III. Timing of long-period pulsars

The Parkes survey of the entire southern sky for millisecond and other low-luminosity pulsars has now been completed. The survey system and initial results were described by Manchester et al. and the final results and population studies are described by Lyne et al. A total of 298 pulsars were detected, including 101 new discoveries of which 17 were millisecond pulsars. Here we report on timing observations at Parkes and Jodrell Bank of the 84 long-period (non-millisecond) pulsars discovered in the survey, including timing solutions for 78 of them. Pulse profiles and arrival times were obtained at several different frequencies over intervals of more than a year, yielding a position, period, period derivative and dispersion measure for each pulsar. Pulse profiles at frequencies near 400 MHz and 600 or 1400 MHz are presented for most of the observed pulsars. Significant timing noise was detected for five pulsars and a glitch was observed in the period of PSR J1123–6259.     

PSR J1453+1902 and the radio luminosities of solitary versus binary millisecond pulsars

We present 3 yr of timing observations for PSR J1453+1902, a 5.79-ms pulsar discovered during a 430-MHz drift-scan survey with the Arecibo telescope. Our observations show that PSR J1453+1902 is solitary and has a proper motion of  8 ±  2  mas yr⁻¹. At the nominal distance of 1.2 kpc estimated from the pulsar's dispersion measure, this corresponds to a transverse speed of  46 ± 11   km s⁻¹  , typical of the millisecond pulsar population. We analyse the current sample of 55 millisecond pulsars in the Galactic disc and revisit the question of whether the luminosities of isolated millisecond pulsars are different from their binary counterparts. We demonstrate that the apparent differences in the luminosity distributions seen in samples selected from 430-MHz surveys can be explained by small-number statistics and observational selection biases. An examination of the sample from 1400-MHz surveys shows no differences in the distributions. The simplest conclusion from the current data is that the spin, kinematic, spatial and luminosity distributions of isolated and binary millisecond pulsars are consistent with a single homogeneous population.     

Low-frequency flux-density measurements and spectral features of millisecond pulsars

Radio flux-density measurements for a large sample of millisecond pulsars at a low frequency of 102 MHz are presented. Using higher frequency measurements, we construct their spectra in the frequency range from 102 MHz to 4.8 GHz, the widest one studied to date. The spectra of millisecond and normal pulsars have been found to differ. The spectra of millisecond pulsars have no low-frequency turnover typical of normal pulsars. The absence of a low-frequency turnover in the spectrum suggests that the emitting regions of millisecond and normal pulsars differ in geometry, which we interpret by deviation of the magnetic field from a dipole one or by compactness of the emitting region.     

Models for the formation of binary and millisecond radio pulsars

The peculiar combination of a relatively short pulse period and a relatively weak surface dipole magnetic field strength of binary radio pulsars finds a consistent explanation in terms of (i) decay of the surface dipole component of neutron-star magnetic fields on a timescale of (2–5) × 10⁶ yr, in combination with (ii) spin-up of the rotation of the neutron star during a subsequent mass-transfer phase. The four known binary radio pulsars appear to fall into two different categories. Two of them, PSR 0655 + 64 and PSR 1913 + 16, have short orbital periods (<25 h) and high mass functions, indicating companion masses 0.7M⊙ (∼1 (± 0.3) M⊙ and 1.4 M⊙, respectively). The other two, PSR 0820 + 02 and PSR 1953 + 29, have long orbital periods (117d), nearly circular orbits, and low, almost identical mass functions of about 3×10^-3 M⊙, suggesting companion masses of about 0.3M⊙. It is pointed out that these two classes of systems are expected to be formed by the later evolution of binaries consisting of a neutron star and a normal companion star, in which the companion was (considerably) more massive than the neutron star, or less massive than the neutron star, respectively. In the first case the companion of the neutron star in the final system will be a massive white dwarf, in a circular orbit, or a neutron star in an eccentric orbit. In the second case the final companion to the neutron star will be a low-mass (∼ 0.3 M⊙) helium white dwarf in a wide and nearly circular orbit. In systems of the second type the neutron star was most probably formed by the accretion-induced collapse of a white dwarf. This explains in a natural way why PSR 1953 + 29 has a millisecond rotation period and PSR 0820 + 02 has not. Among the binary models proposed for the formation of the 1.5-millisecond pulsar, the only ones that appear to be viable are those in which the companion disappeared by coalescence with the neutron star. In such models the companion may have been a red dwarf of mass 0.03M⊙, a neutron star, or a massive (>0.7M⊙) white dwarf. Only in the last-mentioned case is a position of the pulsar close to the galactic plane a natural consequence. In the first-mentioned case the progenitor system most probably was a cataclysmic-variable binary in which the white dwarf collapsed by accretion.     

The population of pulsars with interpulses and the implications for beam evolution

The observed fraction of pulsars with interpulses, their period distribution and the observed pulse width versus pulse period correlation are shown to be inconsistent with a model in which the angle α between the magnetic axis and the rotation axis is random. This conclusion appears to be unavoidable, even when non-circular beams are considered. Allowing the magnetic axis to align from a random distribution at birth with a time-scale of  ∼7 × 10⁷ yr  can, however, explain those observations well. The time-scale derived is consistent with that obtained via independent methods. The probability that a pulsar beam intersects the line of sight is a function of the angle α and therefore beam evolution has important consequences for evolutionary models and for estimations of the total number of neutron stars. The validity of the standard formula for the spin-down rate, which is independent of α, appears to be questionable.