Theoretical Modeling of Rotational Dynamics and Spin Evolution in Neutron Stars

How does the loss of rotational energy shape the spin evolution of a neutron star? Combining conservation of angular momentum and energy with electromagnetic and gravitational-wave torques, we derive analytic models of the rotation period in terms of the magnetic field \( B_{p} \), radius \( R \), moment of inertia \( I \), inclination \( α \) and ellipticity \( ε \). For a force-free dipole, \( P^{2}=P_{0}^{2}+2{\pi }^{2}B_{p}^{2}R^{6}\left(1+\sin^{2}\alpha \right)t/\left(Ic^{3}\right) \). For any torque coefficient \( K\left(t\right) \) the braking index obeys \( n=n_{0}-2{\tau }_{c}\,dlnK/dt \), so field decay gives \( n=3+4{\tau }_{c}/{\tau }_{D} \) and a finite asymptotic period. Combined dipole and gravitational-wave braking is solved exactly; a measured n between 3 and 5 gives the torque ratio \( \eta =\left(n-3\right)/\left(5-n\right) \), and the lifetime gravitational-wave energy fraction is \( 1-\ln\left(1+{\eta }_{0}\right)/{\eta }_{0} \). The coupled spin–inclination evolution conserves \( Ωcosα \) in vacuum and \( \varOmega \cos^{2}\alpha /\sin\alpha \) in force-free magnetospheres, where alignment raises the braking index by at most 1/4. Applied to observed pulsars, the models give an initial period of 19 ms for the Crab and show that braking indices below 3 require torques that grow with time.