Abstract
Almost all pulsars with anomalous positive measurements (corresponding to anomalous braking indices in the range 5 < |n| < 100), including all the pulsars with observed large glitches (ΔΩ/Ω > 10−7) as well as post-glitch or interglitch measurements, obey the scaling between and glitch parameters originally noted in the Vela pulsar. Negative second derivative values can be understood in terms of glitches that were missed or remained unresolved. We discuss the glitch rates and a priori probabilities of positive and negative braking indices according to the model developed for the Vela pulsar. This behaviour supports the universal occurrence of a non-linear dynamical coupling between the neutron star crust and an interior superfluid component. The implied lower limit to dynamical energy dissipation in a neutron star with spindown rate is . Thermal luminosities and surface temperatures due to dynamical energy dissipation are estimated for old neutron stars which are spinning down as rotating magnetic dipoles beyond the pulsar death line.