High resolution, pulse to pulse observation of the 2016 Vela glitch and its relaxation provided an opportunity to probe the neutron star internal structure and dynamics with unprecedented detail. We use the observations of this glitch to infer superfluid characteristics in the framework of the vortex creep model. The glitch rise time constraint of 12.6 seconds put stringent limits on the angular momentum exchange between the crustal superfluid and the observed crust. Together with the observed excess acceleration in the rotation rate as compared to the post-glitch equilibrium value this discriminates crustal superfluid-crust lattice and core superfluid-crustal normal matter coupling time-scales. An evident decrease in the crustal rotation rate immediately before the glitch is consistent with the formation of a new vortex trap zone that initiates the large scale vortex unpinning avalanche. Formation of vortex trap by a crust breaking quake induces short-lived magnetospheric changes. The long term post-glitch spin-down rate evolution reveals the moments of inertia and recoupling time-scales of the superfluid layers participating in the glitch and leads to an estimation of the time to the next glitch which agrees with the time interval between the 2016 and 2019 glitches. Our results are consistent with theoretical estimates of effective neutron and proton masses in the superfluid. We also constrain the vortex line-flux tube pinning energy per intersection as 2 MeV.

The Crab pulsar displayed its largest glitch on 2017 November. An extended initial spin-up phase of this largest glitch was resolved, for the first time with high cadence of observations both in radio and X-rays on a time-scale of 2 days. A combination of crustquake and vortex unpinning models is invoked to account for the extended spin-up, magnitude and post-glitch relaxation characteristics of this glitch. We evaluate the extended spin-up followed by the familiar spin-down as due to the creep response to the initial induced inward motion of some vortex lines pinned to broken crustal plates moving inward towards the rotation axis, together with the common and familiar post-glitch creep response to the sudden outward motion of vortices unpinned at the glitch. Our analysis confirms that the number of unpinned vortices participating in glitches are similar in all Crab glitches, and within an order of magnitude in all glitches from all pulsars. This typical number of unpinned vortices is related to the broken plate size in quakes as triggers for vortex unpinning avalanches. The physical determinant of this universal broken plate size is in turn the critical strain angle in the neutron star crust. Occurrence of this largest Crab glitch after a relatively long inactive period is consistent with accumulation of the pinned vorticity to be tapped.

The two component model of neutron star dynamics describing the behaviour of the observed crust coupled to the superfluid interior has so far been applied to radio pulsars for which the external torques are constant on dynamical timescales. We recently solved this problem under arbitrary time dependent external torques. Our solutions pertain to internal torques that are linear in the rotation rates, as well as to the extremely non-linear internal torques of the vortex creep model. Two-component models with linear or nonlinear internal torques can now be applied to magnetars and to neutron stars in binary systems, with strong variability and timing noise. Time dependent external torques can be obtained from the observed spin-down (or spin-up) time series, $\dot{\Omega}(t)$.

Timing observations of rapidly rotating neutron stars revealed a great number of glitches, observed both from canonical radio pulsars and magnetars. Among them, 76 glitches have shown exponential relaxation(s) with characteristic decay times ranging from several days to a few months, followed by a more gradual recovery. Glitches displaying exponential relaxation with single or multiple decay time constants are analysed in terms of a model based on the interaction of the vortex lines with the toroidal arrangement of flux tubes in the outer core of the neutron star. Model results agree with the observed timescales in general. Thus, the glitch phenomenon can be used to deduce valuable information about neutron star structure, in particular on the interior magnetic field configuration which is unaccessible from surface observations. One immediate conclusion is that the magnetar glitch data are best explained with a much cooler core and therefore require that direct Urca type fast cooling mechanisms should be effective for magnetars.

The two component model describes neutron star dynamics incorporating the response of the superfluid interior. Conventional solutions and applications involve constant external torques, as appropriate for radio pulsars on dynamical timescales. We present the general solution of two component dynamics under arbitrary time dependent external torques, with internal torques that are linear in the rotation rates, or with the extremely non-linear internal torques due to vortex creep. The two-component model incorporating the response of linear or nonlinear internal torques can now be applied not only to radio pulsars but also to magnetars and to neutron stars in binary systems, with strong observed variability and noise in the spin-down or spin-up rates. Our results allow the extraction of the time dependent external torques from the observed spin-down (or spin-up) time series, $\dot{\Omega}(t)$. Applications are discussed.

Treatment of the vortex motion in the superfluids of the inner crust and the outer core of neutron stars is a key ingredient in modeling a number of pulsar phenomena, including glitches and magnetic field evolution. After recalculating the microscopic vortex velocity in the inner crust, we evaluate the velocity for the vortices in the outer core for the first time. The vortex motion between pinning sites is found to be substantially faster in the inner crust than in the outer core, $v_0^{\rm crust} \sim 10^{7}\mbox{\cms} \gg v_0^{\rm core} \sim 1\mbox{\cms}$. One immediate result is that vortex creep is always in the nonlinear regime in the outer core in contrast to the inner crust, where both nonlinear and linear regimes of vortex creep are possible. Other implications for pulsar glitches and magnetic field evolution are also presented.

The recently discovered rotationally powered pulsar PSR J1640-4631 is the first to have a braking index measured, with high enough precision, that is greater than three. An inclined magnetic rotator in vacuum or plasma would be subject not only to spin-down but also to an alignment torque. The vacuum model can address the braking index only for an almost orthogonal rotator that is incompatible with the single peaked pulse profile. The magnetic dipole model with the corotating plasma predicts braking indices between $3-3.25$. We find that the braking index of $3.15$ is consistent with two different inclination angles, $18.5\pm 3$ degrees and $56 \pm 4$ degrees. The smaller angle is preferred given the pulse profile has a single peak and the radio output of the source is weak. We infer the change in the inclination angle to be at the rate $-0.23$ degrees per century, three times smaller in absolute value than the rate recently observed for the Crab pulsar.

Glitches are sudden changes in rotation frequency and spin-down rate, observed from pulsars of all ages. Standard glitches are characterized by a positive step in angular velocity ($\Delta\Omega$ $ > $ $0$) and a negative step in the spin-down rate ($\Delta \dot \Omega$ $ < $ $0$) of the pulsar. There are no glitch-associated changes in the electromagnetic signature of rotation-powered pulsars in all cases so far. For the first time, in the last glitch of PSR J1119-6127, there is clear evidence for changing emission properties coincident with the glitch. This glitch is also unusual in its signature. Further, the absolute value of the spin-down rate actually decreases in the long term. This is in contrast to usual glitch behaviour. In this paper we extend the vortex creep model in order to take into account these peculiarities. We propose that a starquake with crustal plate movement towards the rotational poles of the star induces inward vortex motion which causes the unusual glitch signature. The component of the magnetic field perpendicular to the rotation axis will decrease, giving rise to a permanent change in the pulsar external torque.

A region of toroidally oriented quantized flux lines must exist in the proton superconductor in the core of the neutron star. This region will be a site of vortex pinning and creep. Entrainment of the neutron superfluid with the crustal lattice leads to a requirement of superfluid moment of inertia associated with vortex creep in excess of the available crustal moment of inertia. This will effect constraints on the equation of state. The toroidal flux region provides the moment of inertia necessary to complement the crust superfluid with postglitch relaxation behavior fitting the observations.