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

Characterization of the Crab Pulsar's Timing Noise

We present a power spectral analysis of the Crab pulsar's timing noise, mainly using radio measurements from Jodrell Bank taken over the period 1982-1989. The power spectral analysis is complicated by nonuniform data sampling and the presence of a steep red power spectrum that can distort power spectra measurement by causing severe power ``leakage''. We develop a simple windowing method for computing red noise power spectra of uniformly sampled data sets and test it on Monte Carlo generated sample realizations of red power-law noise. We generalize time-domain methods of generating power-law red noise with even integer spectral indices to the case of noninteger spectral indices. The Jodrell Bank pulse phase residuals are dense and smooth enough that an interpolation onto a uniform time series is possible. A windowed power spectrum is computed revealing a periodic or nearly periodic component with a period of about 568 days and a 1/f^3 power-law noise component with a noise strength of 1.24 +/- 0.067 10^{-16} cycles^2/sec^2 over the analysis frequency range 0.003 - 0.1 cycles/day. This result deviates from past analyses which characterized the pulse phase timing residuals as either 1/f^4 power-law noise or a quasiperiodic process. The analysis was checked using the Deeter polynomial method of power spectrum estimation that was developed for the case of nonuniform sampling, but has lower spectral resolution. The timing noise is consistent with a torque noise spectrum rising with analysis frequency as f implying blue torque noise, a result not predicted by current models of pulsar timing noise. If the periodic or nearly periodic component is due to a binary companion, we find a companion mass > 3.2 Earth masses.Comment: 53 pages, 9 figures, submitted to MNRAS, abstract condense

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

Observations of 20 millisecond pulsars in 47 Tucanae at 20 cm

We have used a new observing system on the Parkes radio telescope to carry out a series of pulsar observations of the globular cluster 47 Tucanae at 20-cm wavelength. We detected all 11 previously known pulsars, and have discovered nine others, all of which are millisecond pulsars in binary systems. We have searched the data for relatively short orbital period systems, and found one pulsar with an orbital period of 96 min, the shortest of any known radio pulsar. The increased rate of detections with the new system resulted in improved estimates of the flux density of the previously known pulsars, determination of the orbital parameters of one of them, and a coherent timing solution for another one. Five of the pulsars now known in 47 Tucanae have orbital periods of a few hours and implied companion masses of only ~ 0.03 Msun. Two of these are eclipsed at some orbital phases, while three are seen at all phases at 20 cm but not always at lower frequencies. Four and possibly six of the other binary systems have longer orbital periods and companion masses ~ 0.2 Msun, with at least two of them having relatively large orbital eccentricities. All 20 pulsars have rotation periods in the range 2-8 ms.Comment: 15 pages, 6 embedded EPS figures, to be published in The Astrophysical Journa

Anomalous scattering of highly dispersed pulsars

We report multifrequency measurements of scatter broadening times for nine highly dispersed pulsars over a wide frequency range (0.6 -- 4.9 GHz). We find the scatter broadening times to be larger than expected and to scale with frequency with an average power-law index of $3.44\pm 0.13$, i.e. significantly less than that expected from standard theories. Such possible discrepancies have been predicted very recently by Cordes & Lazio.Comment: 7 pages, 4 figures, accepted for publication in ApJ Letter

Timing the millisecond pulsars in 47 Tucanae

In the last 10 years 20 millisecond pulsars have been discovered in the globular cluster 47 Tucanae. Hitherto, only 3 of these had published timing solutions. Here we improve upon these 3 and present 12 new solutions. These measurements can be used to determine a variety of physical properties of the pulsars and of the cluster. The 15 pulsars have positions determined with typical uncertianties of only a few milliarcsec and they are all located within 1.2 arcmin of the cluster centre. We have also measured the proper motions of 5 of the pulsars, which are consistent with the proper motion of 47 Tuc based on Hipparcos data. The period derivatives measured for many of the pulsars are dominated by the dynamical effects of the cluster gravitational field, and are used to constrain the surface mass density of the cluster. All pulsars have characteristic ages T > 170 Myr and magnetic fields B < 2.4e9 Gauss, and the average T > 1 Gyr. We have measured the rate of advance of periastron for the binary pulsar J0024-7204H, implying a total system mass 1.4+-0.8 solar masses.Comment: 17 pages, 11 included 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

Nucleon Superfluidity vs Observations of Cooling Neutron Stars

Cooling simulations of neutron stars (NSs) are performed assuming that stellar cores consist of neutrons, protons and electrons and using realistic density profiles of superfluid critical temperatures $T_{cn}(\rho)$ and $T_{cp}(\rho)$ of neutrons and protons. Taking a suitable profile of $T_{cp}(\rho)$ with maximum $\sim 5 \times 10^9$ K one can obtain smooth transition from slow to rapid cooling with increasing stellar mass. Adopting the same profile one can explain the majority of observations of thermal emission from isolated middle--aged NSs by cooling of NSs with different masses either with no neutron superfluidity in the cores or with a weak superfluidity, $T_{cn} < 10^8$ K. The required masses range from $\sim 1.2 M_\odot$ for (young and hot) RX J0822-43 and (old and warm) PSR 1055-52 and RX J1856-3754 to $\approx 1.45 M_\odot$ for the (colder) Geminga and Vela pulsars. Observations constrain the $T_{cn}(\rho)$ and $T_{cp}(\rho)$ profiles with respect to the threshold density of direct Urca process and maximum central density of NSs.Comment: 4 pages, 2 figures, AA Letters, accepte

Isolated pulsar spin evolution on the P-Pdot Diagram

We look at two contrasting spin-down models for isolated radio pulsars and, accounting for selection effects, synthesize observable populations. While our goal is to reproduce all of the observable characteristics, in this paper we pay particular attention to the form of the spin period vs. period derivative (P-Pdot) diagram and its dependence on various pulsar properties. We analyse the initial spin period, the braking index, the magnetic field, various beaming models, as well as the pulsar's luminosity. In addition to considering the standard magnetic dipole model for pulsar spin-down, we also consider the recent hybrid model proposed by Contopoulos & Spitkovsky. The magnetic dipole model, however, does a better job of reproducing the observed pulsar population. We conclude that random alignment angles and period dependent luminosity distributions are essential to reproduce the observed P-Pdot diagram. We also consider the time decay of alignment angles, and attempt to reconcile various models currently being studied. We conclude that, in order to account for recent evidence for the alignment found by Weltevrede & Johnston, the braking torque on a neutron star should not depend strongly on the inclination. Our simulation code is publically available and includes a web-based interface to examine the results and make predictions for yields of current and future surveys.Comment: 9 pages, 4 figure