Glitches are sudden accelerations of the rotation rate $\nu$ of neutron stars, believed to be driven by the neutron superfluid inside the star's crust and core. They present a wide phenomenology and their amplitudes $\Delta\nu$ span over six orders of magnitude. The 21 known glitches of the Vela pulsar are amongst the largest observed (typically $\Delta\nu/\nu\sim10^{-6}$) and have very similar characteristics. We wish to explore the population of small-amplitude rotational changes in the Vela pulsar and determine the rate of occurrence and sizes of its smallest glitches. We use high cadence observations of the Vela pulsar taken between 1981 and 2005 by the Mount Pleasant Radio Observatory. An automated, systematic search was carried out that investigates whether a significant change of spin frequency $\nu$ and/or the spin-down rate $\dot{\nu}$ takes place at any given time. Our study reveals numerous events of all possible signatures (i.e. combinations of $\Delta\nu$ and $\Delta\dot{\nu}$ signs), usually small, with $|\Delta\nu|/\nu<10^{-9}$, which contribute to Vela's timing noise. We also find two glitches that have not been reported before, with respective sizes $\Delta\nu/\nu$ of $(5.55\pm0.03)\times10^{-9}$ and $(38\pm4)\times10^{-9}$. The latter glitch is followed by an exponential-like recovery with characteristic timescale of $\sim30$ d. The Vela pulsar presents an under-abundance of small glitches comparatively to many other glitching pulsars, which appears genuine and not a result of observational biases. Besides typical glitches, the smooth spin down of the pulsar is also affected by an almost continuous activity that can be partially characterised by step-like changes in $\nu$, $\dot{\nu}$ or both.

Pulsar glitches offer an insight into the dynamics of superfluids in the high density interior of a neutron star. To model these phenomena, however, one needs to have an understanding of the dynamics of a turbulent array of superfluid vortices moving through a pinning lattice. In this paper we develop a theoretical approach to describe vortex mediated mutual friction in a pinned, turbulent and rotating superfluid. Our model is then applied to the study of the post glitch rotational evolution in the Vela pulsar and in PSR J0537-6910. We show that in both cases a turbulent model fits the evolution of the spin frequency derivative better than a laminar one. We also predict that the second derivative of the frequency after a glitch should be correlated with the waiting time since the previous glitch, which we find to be consistent with observational data for these pulsars. The main conclusion of this paper is that in the post-glitch rotational evolution of these two pulsars we are most likely observing the response to the glitch of a pinned turbulent region of the star (possibly the crust) and not the laminar response of a regular straight vortex array.

Pulsar glitches are attributed to the sudden re-coupling of very weakly coupled large scale superfluid components in the neutron star interior. This process leads to rapid exchange of angular momentum and an increase in spin frequency. The transfer of angular momentum is regulated by a dissipative mutual friction, whose strength defines the spin-up timescale of a glitch. Hence, observations of glitch rises can be used to shed light on the dominant microphysical interactions at work in the high density interior of the star. We present a simple analytical model, complemented with more detailed numerical simulations, which produces a fast spin-up followed by a more gradual rise. Such features are observed in some large glitches of the Crab pulsar, including the largest recent glitch of 2017. We also use observations to constrain the mutual friction coefficient of the glitch-driving region for two possible locations: the inner crust and outer core of the star. We find that the features of Crab glitches require smaller values of the mutual friction coefficient than those needed to explain the much faster Vela spin-ups. This suggests a crustal origin for the former but an outer core contribution for the latter.

We discuss the unique spin evolution of the young X-ray pulsar PSR J0537-6910, a system in which the regular spin down is interrupted by glitches every few months. Drawing on the complete timing data from the Rossi X-ray Timing Explorer (RXTE, from 1999-2011), we argue that a trend in the inter-glitch behaviour points to an effective braking index close to $n=7$, much larger than expected. This value is interesting because it would accord with the neutron star spinning down due to gravitational waves from an unstable r-mode. We discuss to what extent this, admittedly speculative, scenario may be consistent and if the associated gravitational-wave signal would be within reach of ground based detectors. Our estimates suggest that one may, indeed, be able to use future observations to test the idea. Further precision timing would help enhance the achievable sensitivity and we advocate a joint observing campaign between the Neutron Star Interior Composition ExploreR (NICER) and the LIGO-Virgo network.

The young, fast-spinning, X-ray pulsar J0537-6910 displays an extreme glitch activity, with large spin-ups interrupting its decelerating rotation every ~100 days. We present nearly 13 years of timing data from this pulsar, obtained with the {\it Rossi X-ray Timing Explorer}. We discovered 22 new glitches and performed a consistent analysis of all 45 glitches detected in the complete data span. Our results corroborate the previously reported strong correlation between glitch spin-up size and the time to the next glitch, a relation that has not been observed so far in any other pulsar. The spin evolution is dominated by the glitches, which occur at a rate ~3.5 per year, and the post-glitch recoveries, which prevail the entire inter-glitch intervals. This distinctive behaviour provides invaluable insights into the physics of glitches. The observations can be explained with a multi-component model which accounts for the dynamics of the neutron superfluid present in the crust and core of neutron stars. We place limits on the moment of inertia of the component responsible for the spin-up and, ignoring differential rotation, the velocity difference it can sustain with the crust. Contrary to its rapid decrease between glitches, the spin-down rate increased over the 13 years, and we find the long-term braking index $n_{\rm l}=-1.22(4)$, the only negative braking index seen in a young pulsar. We briefly discuss the plausible interpretations of this result, which is in stark contrast to the predictions of standard models of pulsar spin-down.

Pulsars are rotating neutron stars that are renowned for their timing precision, although glitches can interrupt the regular timing behavior when these stars are young. Glitches are thought to be caused by interactions between normal and superfluid matter in the star. We update our recent work on a new technique using pulsar glitch data to constrain superfluid and nuclear equation of state models, demonstrating how current and future astronomy telescopes can probe fundamental physics such as superfluidity near nuclear saturation and matter at supranuclear densities. Unlike traditional methods of measuring a star's mass by its gravitational effect on another object, our technique relies on nuclear physics knowledge and therefore allows measurement of the mass of pulsars which are in isolation.

Pulsars are known for their superb timing precision, although glitches can interrupt the regular timing behavior when the stars are young. These glitches are thought to be caused by interactions between normal and superfluid matter in the crust of the star. However, glitching pulsars such as Vela have been shown to require a superfluid reservoir that greatly exceeds that available in the crust. We examine a model in which glitches tap the superfluid in the core. We test a variety of theoretical superfluid models against the most recent glitch data and find that only one model can successfully explain up to 45 years of observational data. We develop a new technique for combining radio and X-ray data to measure pulsar masses, thereby demonstrating how current and future telescopes can probe fundamental physics such as superfluidity near nuclear saturation.

Crustquake events may be connected with both rapid spin-up `glitches' within the regular slowdown of neutron stars, and high-energy magnetar flares. We argue that magnetic field decay builds up stresses in a neutron star's crust, as the elastic shear force resists the Lorentz force's desire to rearrange the global magnetic-field equilibrium. We derive a criterion for crust-breaking induced by a changing magnetic-field configuration, and use this to investigate strain patterns in a neutron star's crust for a variety of different magnetic-field models. Universally, we find that the crust is most liable to break if the magnetic field has a strong toroidal component, in which case the epicentre of the crustquake is around the equator. We calculate the energy released in a crustquake as a function of the fracture depth, finding that it is independent of field strength. Crust-breaking is, however, associated with a characteristic local field strength of $2.4\times 10^{14}$ G for a breaking strain of $0.001$, or $2.4\times 10^{15}$ G at a breaking strain of $0.1$. We find that even the most luminous magnetar giant flare could have been powered by crustal energy release alone.

With an average density higher than the nuclear density, neutron stars (NS) provide a unique test-ground for nuclear physics, quantum chromodynamics (QCD), and nuclear superfluidity. Determination of the fundamental interactions that govern matter under such extreme conditions is one of the major unsolved problems of modern physics, and -- since it is impossible to replicate these conditions on Earth -- a major scientific motivation for SKA. The most stringent observational constraints come from measurements of NS bulk properties: each model for the microscopic behaviour of matter predicts a specific density-pressure relation (its `Equation of state', EOS). This generates a unique mass-radius relation which predicts a characteristic radius for a large range of masses and a maximum mass above which NS collapse to black holes. It also uniquely predicts other bulk quantities, like maximum spin frequency and moment of inertia. The SKA, in Phase 1 and particularly in Phase 2 will, thanks to the exquisite timing precision enabled by its raw sensitivity, and surveys that dramatically increase the number of sources: 1) Provide many more precise NS mass measurements (high mass NS measurements are particularly important for ruling out EOS models); 2) Allow the measurement of the NS moment of inertia in highly relativistic binaries such as the Double Pulsar; 3) Greatly increase the number of fast-spinning NS, with the potential discovery of spin frequencies above those allowed by some EOS models; 4) Improve our knowledge of new classes of binary pulsars such as black widows and redbacks (which may be massive as a class) through sensitive broad-band radio observations; and 5) Improve our understanding of dense matter superfluidity and the state of matter in the interior through the study of rotational glitches, provided that an ad-hoc campaign is developed.

Providing a link between magnetars and radio pulsars, high magnetic field neutron stars are ideal targets to investigate how bursting/magnetospheric activity and braking torque variations are connected to rotational glitches. The last spin-up glitch of the highly magnetised pulsar J1119$-$6127 back in 2007 was the first glitch in a rotationally powered radio pulsar to be accompanied by radiative changes. Moreover, it was followed by an uncommon glitch relaxation that resulted in a smaller spin-down rate relative to the prediction of the pre-glitch timing model. Here, we present 4 years of new radio timing observations and analyse the total of 16 years of timing data for this source. The new data uncover an ongoing evolution of the spin-down rate, thereby allowing us to exclude permanent changes in the external or internal torque as a stand-alone cause of the peculiar features of the glitch recovery. Furthermore, no additional variations of the radio pulse profile are detected, strengthening the association of the previously observed transient emission features with the glitching activity. A self-consistent measurement of the braking index yields a value $n\simeq2.7$, indicating a trajectory in the $P-\dot{P}$ plane inclined towards the magnetars. Such a potential evolutionary link might be strengthened by a, possibly permanent, reduction of $\sim15\%$ in $n$ at the epoch of the 2007 glitch.

Glitches are sudden spin-up events that punctuate the steady spin down of pulsars and are thought to be due to the presence of a superfluid component within neutron stars. The precise glitch mechanism and its trigger, however, remain unknown. The size of glitches is a key diagnostic for models of the underlying physics. While the largest glitches have long been taken into account by theoretical models, it has always been assumed that the minimum size lay below the detectability limit of the measurements. In this paper we define general glitch detectability limits and use them on 29 years of daily observations of the Crab pulsar, carried out at Jodrell Bank Observatory. We find that all glitches lie well above the detectability limits and by using an automated method to search for small events we are able to uncover the full glitch size distribution, with no biases. Contrary to the prediction of most models, the distribution presents a rapid decrease of the number of glitches below ~0.05 $\mu$Hz. This substantial minimum size indicates that a glitch must involve the motion of at least several billion superfluid vortices and provides an extra observable which can greatly help the identification of the trigger mechanism. Our study also shows that glitches are clearly separated from all the other rotation irregularities. This supports the idea that the origin of glitches is different to that of timing noise, which comprises the unmodelled random fluctuations in the rotation rates of pulsars.

Pulsar glitches are sudden increases in the spin frequency of an otherwise steadily spinning down neutron star. These events are thought to represent a direct probe of the dynamics of the superfluid interior of the star. However glitches can differ significantly from one another, not only in size and frequency, but also in the post-glitch response of the star. Some appear as simple steps in frequency, while others also display an increase in spin-down rate after the glitch. Others still show several exponentially relaxing components in the post-glitch recovery. We show that if glitches are indeed due to large scale unpinning of superfluid vortices, the different regions in which this occurs and respective timescales on which they recouple can lead to the various observed signatures. Furthermore we show that this framework naturally accounts for the peculiar relaxations of glitches in Anomalous X-ray Pulsars.