Observations of pulsar glitches have the potential to provide constraints on the dynamics of the high density interior of neutron stars. However, to do so, realistic glitch models must be constructed and compared to the data. We take a step towards this goal by testing non-linear models for the mutual friction force, which is responsible for the exchange of angular momentum between the neutron superfluid and the observable normal component in a glitch. In particular, we consider a non-linear dependence of the drag force on the relative velocity between superfluid vortices and the normal component, in which the contributions of both kelvon and phonon excitations are included. This non-linear model produces qualitatively new features, and is able to reproduce the observed bimodal distribution of glitch sizes in the pulsar population. The model also suggests that the differences in size distributions in individual pulsars may be due to the glitches being triggered in regions with different pinning strengths, as stronger pinning leads to higher vortex velocities and a qualitatively different mutual friction coupling with respect to the weak pinning case. Glitches in pulsars that appear to glitch quasi-periodically with similar sizes may thus be due to the same mechanisms as smaller events in pulsars that have no preferred glitch size, but simply originate in stronger pinning regions, possibly in the core of the star.

The interiors of mature neutron stars are expected to be superfluid. Superfluidity of matter on the microscopic scale can have a number of large scale, potentially observable consequences, as the superfluid component of the star can now flow relative to the `normal' component that is tracked by electromagnetic emission. The most spectacular of such phenomena are pulsar glitches, sudden spin-ups observed in many pulsars. A background flow of a normal-fluid component with respect to the superfluid is also known to lead to a number of instabilities in laboratory superfluids, possibly leading to turbulence and modifying the nature of the mutual friction coupling between the two fluids. In this paper we consider modes of oscillation in the crust and core of a superfluid neutron star, by conducting a plane-wave analysis. We explicitly account for a background flow between the two components (as would be expected in the presence of pinning) and the entrainment, and we consider both standard (Hall-Vinen) and isotropic (Gorter-Mellink) forms of the mutual friction. We find that for standard mutual friction there are families of unstable inertial and sound waves both in the case of a counter-flow along the superfluid vortex axis and for counterflow perpendicular to the vortex axis and find that entrainment leads to a quantitative difference between instabilities in the crust and core of the star. For isotropic mutual friction we find no unstable modes, and speculate that instabilities in a straight vortex array may be linked to glitching behaviour, which then ceases until the turbulence has decayed.

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.

The dynamics of quantised vorticity in neutron star interiors is at the heart of most pulsar glitch models. However, the large number of vortices (up to $\approx 10^{13}$) involved in a glitch and the huge disparity in scales between the femtometer scale of vortex cores and the kilometre scale of the star makes quantum dynamical simulations of the problem computationally intractable. In this paper we take a first step towards developing a mean field prescription to include the dynamics of vortices in large scale hydrodynamical simulations of superfluid neutron stars. We consider a one dimensional setup and show that vortex accumulation and differential rotation in the neutron superfluid lead to propagating waves, or `avalanches', as solutions for the equations of motion for the superfluid velocities. We introduce an additional variable, the fraction of free vortices, and test different prescriptions for its advection with the superfluid flow. We find that the new terms lead to solutions with a linear component in the rise of a glitch, and that, in specific setups, they can give rise to glitch precursors and even to decreases in frequency, or `anti-glitches'.