The ubiquity of variable radio emission and spin-down rates in pulsars

Pulsars are often lauded for their (relative) rotational and radio emission stability over long time scales. However, long-term observing programmes are identifying an increasing number of pulsars that deviate from this preconceived notion. Using Gaussian process regression and Bayesian inference techniques, we investigated the emission and rotational stability of 259 isolated radio pulsars that have been monitored using Murriyang, the Parkes 64 m radio telescope, over the past three decades. We found that 238 pulsars display significant variability in their spin-down rates, 52 of which also exhibit changes in profile shape. Including 23 known state-switching pulsars, this represents the largest catalogue of variable pulsars identified to date and indicates that these behaviours are ubiquitous among the wider population. The intensity of spin-down fluctuations positively scales with increasing pulsar spin-down rate, with only a marginal dependence on spin-frequency. This may have substantial implications for ongoing searches for gravitational waves in the ensemble timing of millisecond pulsars. We also discuss challenges in explaining the physical origins of quasi-periodic and transient profile/spin-down variations detected among a subset of our pulsars.