The rotational evolution of pulsars is governed by torque mechanisms whose mathematical structure encodes fundamental symmetries of the underlying physics. We demonstrate that the standard spin-down equation f˙=−sf−rf3−gf5 derives from a discrete antisymmetry requirement, namely invariance of the torque under reversal of rotation sense, which restricts the frequency dependence to odd integer powers. We show that physically motivated plasma processes systematically break this symmetry, introducing fractional frequency exponents: viscous Ekman pumping at the crust–superfluid boundary layer (f3/2), magnetohydrodynamic turbulent dissipation via Kolmogorov and Sweet–Parker cascades (f10/3, f11/3), non-linear superfluid vortex dynamics (f5/2), and saturated r-mode oscillations (f7−2β). The central result is an exact analytical resolution of the long-standing Crab pulsar braking index puzzle: the observed n=2.51±0.01, which has defied explanation for nearly four decades, emerges naturally from the superposition of magnetic dipole radiation (f˙∝f3) and boundary layer Ekman pumping (f˙∝f3/2), with analytically derived coefficients yielding a dipole-component surface field Bp=6.2×1012 G—higher than the standard PP˙ estimate of 3.8×1012 G, because that formula conflates dipole and non-dipole torques, but lower than applying the Larmor formula to the full spin-down rate (7.6×1012 G), since 32.7% of the total torque is non-radiative boundary-layer dissipation. We develop the Riemann–Liouville fractional calculus formalism for these equations, showing that fractional derivatives break time-translation symmetry through intrinsic memory effects, with solutions expressed in terms of Mittag-Leffler and Fox H-functions that interpolate continuously between exponential (fully symmetric) and power-law (scale-free symmetric) relaxation. Lambert–Tsallis Wq functions with non-extensive parameter q encoding broken statistical symmetry enable equation-of-state-independent inference of neutron star compactness and tidal deformability. Our framework establishes a unified symmetry-based classification of pulsar spin-down mechanisms and predicts frequency-dependent braking indices evolving at rate dn/dt∼2×10−4 yr−1, yielding Δn≈0.01 over 50 years—testable with current pulsar timing programmes. The formalism provides a coherent theoretical foundation connecting plasma microphysics at the neutron star interior to macroscopic observables in electromagnetic and gravitational wave channels.
Chishtie et al. (Mon,) studied this question.