The observation of the first pulse-to-pulse glitch in the Vela pulsar opens a new window on theoretical speculations about the internal dynamics of neutron stars, as it allows to test models for the description of the first moments in a glitch. We improve the analytical study of the minimal three-component model for pulsar glitches, by solving it with generic initial conditions for the two initial lags of the superfluid components. The purpose is to use this solution to fit the data of the 2016 Vela glitch - by employing a Bayesian approach - and to obtain a probability distribution for the physical parameters of the model and for observational parameters like the glitch rise time and the relaxation timescale. A physically reasonable, non-informative prior has been set on the different parameters of the model, so that the posterior distribution can be compared with state-of-the-art information obtained from microphysical calculations. By considering a model with a tightened prior on the moment of inertia fractions and by comparing it with the original model by means of Bayesian model selection, we study the possibility of a crust-limited superfluid reservoir. The posterior distribution for the moment of inertia fractions of the superfluid components, the coupling parameters and the initial velocity lags between the components has been obtained. Analysis of the inferred posterior also confirms the presence of an overshoot in that glitch, and sets an upper limit of 6 seconds on the glitch rise timescale. The comparison between the two models with different priors on the moment of inertia fractions seems to indicate the necessity of a core participation in the glitch phenomenon, regardless of the uncertain strength of the entrainment coupling.

During the spin-up phase of a large pulsar glitch - a sudden decrease of the rotational period of a neutron star - the angular velocity of the star may overshoot, namely reach values greater than that observed for the new post-glitch equilibrium. These transient phenomena are expected on the basis of theoretical models for pulsar internal dynamics and their observation has the potential to provide an important diagnostic for glitch modelling. We present a simple criterion to assess the presence of an overshoot, based on the minimal analytical model that is able to reproduce an overshooting spin-up. We employ it to fit the data of the 2016 glitch of the Vela pulsar, obtaining estimates of the fractional moments of inertia of the internal superfluid components involved in the glitch, of the rise and decay timescales of the overshoot, and of the mutual friction parameters between the superfluid components and the normal one. We study the cases with and without strong entrainment in the crust: in the former, we find indication of a large inner core strongly coupled to the observable component and of a reservoir of angular momentum extending into the core to densities below nuclear saturation, while in the latter a large reservoir extending above nuclear saturation and a standard normal component without inner core are suggested.

Observations of pulsar glitches may provide insights on the internal physics of neutron stars and recent studies show how it is in principle possible to constrain pulsar masses with timing observations. The reliability of these estimates depend on the current uncertainties about the structure of neutron stars and on our ability to model the dynamics of the superfluid neutrons in the internal layers. We assume a simplified model for the rotational dynamics of a neutron star and estimate an upper bound to the mass of 25 pulsars from their largest glitch and average activity: the aim is to understand to which extent the mass constraints are sensitive to the choice of the unknown structural properties of neutron stars, like the extension of the superfluid region and the equation of state. Reasonable values, within the range measured for neutron star masses, are obtained only if the superfluid domain extends for at least a small region inside the outer core, which is compatible with calculations of the neutron S-wave pairing gap. Moreover, the mass constraints stabilise when the superfluid domain extends to densities over nuclear saturation, irrespective of the equation of state tested.

Pinning of vortex lines in the inner crust of a spinning neutron star may be the mechanism that enhances the differential rotation of the internal neutron superfluid, making it possible to freeze some amount of angular momentum which eventually can be released, thus causing a pulsar glitch. We investigate the general relativistic corrections to pulsar glitch amplitudes in the slow-rotation approximation, consistently with the stratified structure of the star. We thus provide a relativistic generalization of a previous Newtonian model that was recently used to estimate upper bounds on the masses of glitching pulsars. We find that the effect of general relativity on the glitch amplitudes obtained by emptying the whole angular momentum reservoir is less than 30\%. Moreover we show that the Newtonian upper bounds on the masses of large glitchers obtained from observations of their largest recorded event differ by less than a few percent from those calculated within the relativistic framework. This work can also serve as a basis to construct more sophisticated models of angular momentum reservoir in a relativistic context. In particular, we present two alternative scenarios for "rigid" and "slack" vortex lines, and we generalize the Feynman-Onsager relation to the case when both entrainment coupling between the fluids and a strong gravitational field are present.