2008MNRAS.390..175W

A cellular automaton model of pulsar glitches is described, based on the superfluid vortex unpinning paradigm. Recent analyses of pulsar glitch data suggest that glitches result from scale-invariant avalanches, which are consistent with a self-organized critical system (SOCS). A cellular automaton provides a computationally efficient means of modelling the collective behaviour of up to 10¹⁶ vortices in the pulsar interior, whilst ensuring that the dominant aspects of the microphysics are not lost. The automaton generates avalanche distributions that are qualitatively consistent with a SOCS and with glitch data. The probability density functions of glitch sizes and durations are power laws, and the probability density function of waiting times between successive glitches is Poissonian, consistent with statistically independent events. The output of the model depends on the physical and computational parameters used. The fitted power-law exponents a and b (the size and duration distributions, respectively) decrease as the strength of the vortex pinning increases. Similarly, the exponents increase as the fraction of vortices that are pinned decreases. For the physical and computational parameters considered, one finds -4.3 ≤ a ≤ -2.0, -5.5 ≤ b ≤ -2.2, and mean glitching rates in the range 0.0023 ≤ λ ≤ 0.13 in units of inverse time.