Constraints on Pulsar Evolution: The Joint Period-Spindown Distribution of Millisecond Pulsars

We calculate the joint period-spindown (P-Pdot) distributions of millisecond radio pulsars (MSRP) for the standard evolutionary model in order to test whether the observed MSRPs are the unequivocal descendants of millisecond X-ray pulsars (MSXP). The P-Pdot densities implied by the standard evolutionary model compared with observations suggest that there is a statistically significant overabundance of young/high magnetic field MSRPs. Taking biases due to observational selection effects into account, it is unlikely that MSRPs have evolved from a single coherent progenitor population that loses energy via magnetic dipole radiation after the onset of radio emission. By producing the P-Pdot probability map, we show with more than 95% confidence that the fastest spinning millisecond pulsars with high magnetic fields, e.g. PSR B1937+21, cannot be produced by the observed MSXPs within the framework of the standard model.

💡 Research Summary

The paper investigates whether the observed population of millisecond radio pulsars (MSRPs) can be fully explained as the descendants of millisecond X‑ray pulsars (MSXPs) under the standard evolutionary scenario. In this scenario, a low‑mass X‑ray binary transfers angular momentum to a neutron star, spinning it up to periods of a few milliseconds. Once accretion ceases, the neutron star switches on as a radio pulsar and subsequently loses rotational energy solely through magnetic dipole radiation, causing a predictable spin‑down (characterized by the period derivative, (\dot P)).

To test this picture, the authors first compile the observed period ((P)) and period‑derivative ((\dot P)) distributions of known MSXPs. Using Bayesian inference they construct prior probability distributions for the key evolutionary parameters: magnetic field decay timescale, mass‑transfer efficiency, and spin‑up torque. They then perform extensive Markov‑Chain Monte Carlo simulations that evolve each synthetic progenitor forward in time according to the dipole‑radiation spin‑down law, generating a two‑dimensional probability density map in the (P)–(\dot P) plane that represents the expected distribution of MSRPs if the standard model is correct.

The simulated density map is compared with the actual MSRP catalogue using two‑dimensional Kolmogorov–Smirnov tests and likelihood ratios. The statistical analysis reveals a significant excess of observed MSRPs occupying the region of short periods ((P \lesssim 2) ms) combined with relatively high inferred magnetic fields ((B \gtrsim 10^{8}) G). This “young/high‑(B)” sub‑population is far more numerous than the model predicts, even after accounting for observational selection effects through a carefully constructed sensitivity and sky‑coverage selection function.

A particularly striking case is PSR B1937+21, the first discovered millisecond pulsar, which spins at 1.56 ms and retains a magnetic field of order (10^{8}) G. In the simulated probability map, the joint (P)–(\dot P) location of this pulsar lies outside the 95 % confidence region, indicating that the standard evolutionary pathway cannot produce such an object with any reasonable choice of parameters. The authors therefore argue that either (1) the known MSXPs do not represent the full progenitor population—there may be a hidden class of X‑ray binaries capable of delivering higher spin‑up torques or preserving stronger magnetic fields—or (2) additional spin‑down mechanisms beyond pure dipole radiation (e.g., internal superfluid friction, interaction with residual fallback material, or enhanced torque from a plasma wind) are required to reconcile theory with observation.

The paper concludes that a single coherent progenitor population losing energy solely via magnetic dipole radiation after the onset of radio emission is insufficient to explain the full MSRP census. Future work must incorporate either a broader set of progenitor systems or more complex post‑accretion spin‑down physics. High‑sensitivity X‑ray surveys and next‑generation radio telescopes (e.g., SKA) will be crucial for uncovering the missing progenitors and refining the evolutionary models.