Time evolution of pulsar obliquity angle from 3D simulations of
  magnetospheres
  The rotational period of isolated pulsars increases over time due to the
extraction of angular momentum by electromagnetic torques. These torques also
change the obliquity angle $\alpha$ between the magnetic and rotational axes.
Although actual pulsar magnetospheres are plasma-filled, the time evolution of
$\alpha$ has mostly been studied for vacuum pulsar magnetospheres. In this
work, we self-consistently account for the plasma effects for the first time by
analysing the results of time-dependent 3D force-free and magnetohydrodynamic
simulations of pulsar magnetospheres. We show that if a neutron star is
spherically symmetric and is embedded with a dipolar magnetic moment, the
pulsar evolves so as to minimise its spin-down luminosity: both vacuum and
plasma-filled pulsars evolve toward the aligned configuration ($\alpha=0$).
However, they approach the alignment in qualitatively different ways. Vacuum
pulsars come into alignment exponentially fast, with $\alpha \propto
\exp(-t/\tau)$ and $\tau \sim$ spindown timescale. In contrast, we find that
plasma-filled pulsars align much more slowly, with $\alpha \propto
(t/\tau)^{-1/2}$. We argue that the slow time evolution of obliquity of
plasma-filled pulsars can potentially resolve several observational puzzles,
including the origin of normal pulsars with periods of $\sim1$ second, the
evidence that oblique pulsars come into alignment over a timescale of $\sim
10^7$ years, and the observed deficit, relative to an isotropic obliquity
distribution, of pulsars showing interpulse emission.
