The Gamma Ray Pulsar Population

We apply a likelihood analysis to pulsar detections, pulsar upper limits, and diffuse background measurements from the OSSE and EGRET instruments on the Compton Gamma Ray Observatory to constrain the luminosity law for gamma-ray pulsars and some properties of the gamma-ray pulsar population. We find that the dependence of luminosity on spin period and dipole magnetic field is much steeper at OSSE than at EGRET energies (50-200 keV and >100 MeV, respectively), suggesting that different emission mechanisms are responsible for low- and high-energy gamma-ray emission. Incorporating a spin-down model and assuming a pulsar spatial distribution, we estimate the fraction of the Galactic gamma-ray background due to unidentified pulsars and find that pulsars may be an important component of the OSSE diffuse flux, but are most likely not important at EGRET energies. Using measurements of the diffuse background flux from these instruments, we are able to place constraints on the braking index, initial spin period, and magnetic field of the Galactic pulsar population. We are also able to constrain the pulsar birthrate to be between 1/(25 yr) and 1/(500 yr). Our results are based on a large gamma-ray beam, but they do not scale in a simple way with beam size. With our assumed beam size, the implied gamma-ray efficiency for the EGRET detections is no more than 20%. We estimate that about 20 of the 169 unidentified EGRET sources are probably gamma-ray pulsars. We use our model to predict the pulsar population that will be seen by future gamma-ray instruments and estimate that GLAST will detect roughly 750 gamma-ray pulsars as steady sources, only 120 of which are currently known radio pulsars.Comment: 32 pages, including figures. submitted to Ap

Absolute Timing of the Crab Pulsar with RXTE

We have monitored the phase of the main X-ray pulse of the Crab pulsar with the Rossi X-ray Timing Explorer (RXTE) for almost eight years, since the start of the mission in January 1996. The absolute time of RXTE's clock is sufficiently accurate to allow this phase to be compared directly with the radio profile. Our monitoring observations of the pulsar took place bi-weekly (during the periods when it was at least 30 degrees from the Sun) and we correlated the data with radio timing ephemerides derived from observations made at Jodrell Bank. We have determined the phase of the X-ray main pulse for each observation with a typical error in the individual data points of 50 us. The total ensemble is consistent with a phase that is constant over the monitoring period, with the X-ray pulse leading the radio pulse by 0.0102+/-0.0012 period in phase, or 344+/-40 us in time. The error estimate is dominated by a systematic error of 40 us in the radio data, arising from uncertainties in the variable amount of pulse delay due to interstellar scattering and instrumental calibration. The statistical error is 0.00015 period, or 5 us. The separation of the main pulse and interpulse appears to be unchanging at time scales of a year or less, with an average value of 0.4001+/-0.0002 period. There is no apparent variation in these values with energy over the 2-30 keV range. The lag between the radio and X-ray pulses may be constant in phase (rotational) or constant in time (linear pathlength). We are not (yet) able to distinguish between these two interpretations.Comment: 11 pages, 2 figure

Neutron star magnetic field evolution, crust movement and glitches

Spinning superfluid neutrons in the core of a neutron star interact strongly with co-existing superconducting protons. One consequence is that the outward(inward) motion of core superfluid neutron vortices during spin-down(up) of a neutron star may alter the core's magnetic field. Such core field changes are expected to result in movements of the stellar crust and changes in the star's surface magnetic field which reflect those in the core below. Observed magnitudes and evolution of the spin-down indices of canonical pulsars are understood as a consequence of such surface field changes. If the growing crustal strains caused by the changing core magnetic field configuration in canonical spinning-down pulsars are relaxed by large scale crust-cracking events, special properties are predicted for the resulting changes in spin-period. These agree with various glitch observations, including glitch activity, permanent shifts in spin-down rates after glitches in young pulsars, the intervals between glitches, families of glitches with different magnitudes in the same pulsar, the sharp drop in glitch intervals and magnitudes as pulsar spin-periods approach 0.7s, and the general absence of glitching beyond this period.Comment: LaTex, 28 pages, 8 figs, accepted for publication in Ap

Are We Seeing Magnetic Axis Reorientation in the Crab and Vela Pulsars?

Variation in the angle $\alpha$ between a pulsar's rotational and magnetic axes would change the torque and spin-down rate. We show that sudden increases in $\alpha$, coincident with glitches, could be responsible for the persistent increases in spin-down rate that follow glitches in the Crab pulsar. Moreover, changes in $\alpha$ at a rate similar to that inferred for the Crab pulsar account naturally for the very low braking index of the Vela pulsar. If $\alpha$ increases with time, all pulsar ages obtained from the conventional braking model are underestimates. Decoupling of the neutron star liquid interior from the external torque cannot account for Vela's low braking index. Variations in the Crab's pulse profile due to changes in $\alpha$ might be measurable.Comment: 14 pages and one figure, Latex, uses aasms4.sty. Accepted to ApJ Letter

Temporal variations in scattering and dispersion measure in the Crab Pulsar and their effect on timing precision

We have measured variations in scattering time scales in the Crab Pulsar over a 30-year period, using observations made at 610 MHz with the 42-ft telescope at Jodrell Bank Observatory. Over more recent years, where regular Lovell Telescope observations at frequencies around 1400 MHz were available, we have also determined the dispersion measure variations, after disentangling the scattering delay from the dispersive delay. We demonstrate a relationship between scattering and dispersion measure variations, with a correlation coefficient of $0.56\pm0.01$. The short time scales over which these quantities vary, the size of the variations, and the close correlation between scattering and dispersion measure all suggest that the effects are due to discrete structures within the Crab Nebula, with size scales of $\sim6$ AU (corresponding to an angular size of $\sim2$ mas at an assumed distance of 2200 pc). We mitigate the effects of scattering on the observed pulse shape by using the measured scattering information to modify the template used for generating the pulse arrival times, thus improving the precision to which the pulsar can be timed. We test this on timing data taken during periods of high scattering, and obtain a factor of two improvement in the root mean square of the timing residuals.Comment: 10 pages, 7 figures. Accepted for publication in MNRA

Birth and Evolution of Isolated Radio Pulsars

We investigate the birth and evolution of Galactic isolated radio pulsars. We begin by estimating their birth space velocity distribution from proper motion measurements of Brisken et al. (2002, 2003). We find no evidence for multimodality of the distribution and favor one in which the absolute one-dimensional velocity components are exponentially distributed and with a three-dimensional mean velocity of 380^{+40}_{-60} km s^-1. We then proceed with a Monte Carlo-based population synthesis, modelling the birth properties of the pulsars, their time evolution, and their detection in the Parkes and Swinburne Multibeam surveys. We present a population model that appears generally consistent with the observations. Our results suggest that pulsars are born in the spiral arms, with a Galactocentric radial distribution that is well described by the functional form proposed by Yusifov & Kucuk (2004), in which the pulsar surface density peaks at radius ~3 kpc. The birth spin period distribution extends to several hundred milliseconds, with no evidence of multimodality. Models which assume the radio luminosities of pulsars to be independent of the spin periods and period derivatives are inadequate, as they lead to the detection of too many old simulated pulsars in our simulations. Dithered radio luminosities proportional to the square root of the spin-down luminosity accommodate the observations well and provide a natural mechanism for the pulsars to dim uniformly as they approach the death line, avoiding an observed pile-up on the latter. There is no evidence for significant torque decay (due to magnetic field decay or otherwise) over the lifetime of the pulsars as radio sources (~100 Myr). Finally, we estimate the pulsar birthrate and total number of pulsars in the Galaxy.Comment: 27 pages, including 15 figures, accepted by Ap