Vortex-mediated mutual friction governs the coupling between the superfluid and normal components in neutron star interiors. By, for example, comparing precise timing observations of pulsar glitches with theoretical predictions it is possible to constrain the physics in the interior of the star, but to do so an accurate model of the mutual friction coupling in general relativity is needed. We derive such a model directly from Carter's multifluid formalism, and study the vortex structure and coupling time-scale between the components in a relativistic star. We calculate how general relativity modifies the shape and the density of the quantized vortices and show that, in the quasi-Schwarzschild coordinates, they can be approximated as straight lines for realistic neutron star configurations. Finally, we present a simple universal formula (given as a function of the stellar compactness alone) for the relativistic correction to the glitch rise-time, which is valid under the assumption that the superfluid reservoir is in a thin shell in the crust or in the outer core. This universal relation can be easily employed to correct, a posteriori, any Newtonian estimate for the coupling time-scale, without any additional computational expense.

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.

We introduce a Newtonian model for the deformations of a compressible, autogravitating, and continuously stratified neutron star. The present framework can be applied to a number of astrophysical scenarios as it allows us to account for a great variety of loading forces. In this first analysis, the model is used to study the impact of a frozen adiabatic index in the estimate of rotation-induced deformations: we assume a polytropic equation of state for the matter at equilibrium but, since chemical reactions may be slow, the perturbations with respect to the unstressed configuration are modelled by using a different adiabatic index. We quantify the impact of a departure of the adiabatic index from its equilibrium value on the stressed stellar configuration and we find that a small perturbation can cause large variations both in displacements and strains. As a first practical application, we estimate the strain developed between two large glitches in the Vela pulsar showing that, starting from an initial unstressed configuration, it is not possible to reach the breaking threshold of the crust, namely to trigger a starquake. In this sense, the hypothesis that starquakes could trigger the unpinning of superfluid vortices is challenged and, for the quake to be a possible trigger, the solid crust must never fully relax after a glitch, making the sequence of starquakes in a neutron star an history-dependent process.

We study a class of Newtonian models for the deformations of non-magnetized neutron stars during their spin-down. The models have all an analytical solution, and thus allow to understand easily the dependence of the strain on the star's main physical quantities, such as radius, mass and crust thickness. In the first ``historical'' model the star is assumed to be comprised of a fluid core and an elastic crust with the same density. We compare the response of stars with different masses and equations of state to a decreasing centrifugal force, finding smaller deformations for heavier stars: the strain angle is peaked at the equator and turns out to be a decreasing function of the mass.We introduce a second, more refined, model in which the core and the crust have different densities and the gravitational potential of the deformed body is self-consistently accounted for. Also in this case the strain angle is a decreasing function of the stellar mass, but its maximum value is at the poles and is always larger than the corresponding one in the one-density model by a factor of two. Finally, within the present analytic approach, it is possible to estimate easily the impact of the Cowling approximation: neglecting the perturbations of the gravitational potential, the strain angle is 40\% of the one obtained with the complete model.

Pulsar glitches, sudden jumps in frequency in otherwise steadily spinning down radio pulsars, offer a unique glimpse into the superfluid interior of neutron stars. The exact trigger of these events remains, however, elusive and this has hampered attempts to use glitch observations to constrain fundamental physics. In this paper we propose a new method to measure the mass of glitching pulsars, using observations of the maximum glitch recorded in a star, together with state of the art microphysical models of the pinning interaction between superfluid vortices and ions in the crust. Studying systematically all the presently observed large glitchers, we find an inverse correlation between size of the maximum glitch and the pulsar mass. Our procedure will allow current and future observations of glitching pulsars to constrain not only the physics of glitch models but also the equation of state of dense matter in neutron star interiors.

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.

Large pulsar glitches (like the ones detected in the Vela) are though to be a consequence of the superfluid component present in the interior of mature neutron stars. However strong entrainment challenges this picture. We study the impact of entrainment on the angular momentum that can be exchanged between the normal component and the superfluid during a glitch by means of a consistent global model. This allows to estimate the maximum angular momentum reservoir stored into the superfluid component of the star: the essential ingredient are newly calculated mesoscopic pinning forces that block the superfluid vorticity in the crust of the neutron star. This method can also provide a quantitative test for global models of rotating neutron stars, as well as for microphysical inputs present in literature (like entrainment parameters and pinning forces).

In this lecture, we give a first introduction to neutron stars, based on fundamental physical principles. After outlining their outstanding macroscopic properties, as obtained from observations, we infer the extreme conditions of matter in their interiors. We then describe two crucial physical phenomena which characterize compact stars, namely the gravitational stability of strongly degenerate matter and the neutronization of nuclear matter with increasing density, and explain how the formation and properties of neutron stars are a direct consequence of the extreme compression of matter under strong gravity. Finally, we describe how multi-wavelength observations of different external macroscopic features (e.g. maximum mass, surface temperature, pulsar glitches) can give invaluable information about the exotic internal microscopic scenario: super-dense, isospin-asymmetric, superfluid, bulk hadronic matter (probably deconfined in the most central regions) which can be found nowhere else in the Universe. Indeed, neutrons stars represent a unique probe to study the properties of the low-temperature, high-density sector of the QCD phase diagram. Moreover, binary systems of compact stars allow to make extremely precise measurements of the properties of curved space-time in the strong field regime, as well as being efficient sources of gravitational waves.

In this article we present an analytical two-component model for pulsar rotational dynamics. Under the assumption of axial symmetry, implemented by a paraxial array of straight vortices that thread the entire neutron superfluid, we are able to project exactly the 3D hydrodynamical problem to a 1D cylindrical one. In the presence of density dependent entrainment the superfluid rotation is non-columnar: we circumvent this by using an auxiliary dynamical variable directly related to the areal density of vortices. The main result is a system of differential equations that take consistently into account the stratified spherical structure of the star, the dynamical effects of non-uniform entrainment, the differential rotation of the superfluid component and its coupling to the normal crust. These equations represent a mathematical framework in which to test quantitatively the macroscopic consequences of the presence of a stable vortex array, a working hypothesis widely used in glitch models. Even without solving the equations explicitly, we are able to draw some general quantitative conclusions; in particular, we show that the reservoir of angular momentum (corresponding to recent values of the pinning forces) is enough to reproduce the largest glitch observed in the Vela pulsar, provided its mass is not too large.

The crust of a neutron star is thought to be comprised of a lattice of nuclei immersed in a sea of free electrons and neutrons. As the neutrons are superfluid their angular momentum is carried by an array of quantized vortices. These vortices can pin to the nuclear lattice and prevent the neutron superfluid from spinning down, allowing it to store angular momentum which can then be released catastrophically, giving rise to a pulsar glitch. A crucial ingredient for this model is the maximum pinning force that the lattice can exert on the vortices, as this allows us to estimate the angular momentum that can be exchanged during a glitch. In this paper we perform, for the first time, a detailed and quantitative calculation of the pinning force \emph{per unit length} acting on a vortex immersed in the crust and resulting from the mesoscopic vortex-lattice interaction. We consider realistic vortex tensions, allow for displacement of the nuclei and average over all possible orientation of the crystal with respect to the vortex. We find that, as expected, the mesoscopic pinning force becomes weaker for longer vortices and is generally much smaller than previous estimates, based on vortices aligned with the crystal. Nevertheless the forces we obtain still have maximum values of order $f_{\rm{pin}}\approx 10^{15}$ dyn/cm, which would still allow for enough angular momentum to be stored in the crust to explain large Vela glitches, if part of the star is decoupled during the event.

Several glitches have been observed in young, isolated radio pulsars, while a clear detection in accretion-powered X-ray pulsars is still lacking. We use the Pizzochero snowplow model for pulsar glitches as well as starquake models to determine for the first time the expected properties of glitches in accreting pulsars and their observability. Since some accreting pulsars show accretion-induced long-term spin-up, we also investigate the possibility that anti-glitches occur in these stars. We find that glitches caused by quakes in a slow accreting neutron star are very rare and their detection extremely unlikely. On the contrary, glitches and anti-glitches caused by a transfer of angular momentum between the superfluid neutron vortices and the non-superfluid component may take place in accreting pulsars more often. We calculate the maximum jump in angular velocity of an anti-glitch and we find that it is expected to be about 1E-5 - 1E-4 rad/s. We also note that since accreting pulsars usually have rotational angular velocities lower than those of isolated glitching pulsars, both glitches and anti-glitches are expected to have long rise and recovery timescales compared to isolated glitching pulsars, with glitches and anti-glitches appearing as a simple step in angular velocity. Among accreting pulsars, we find that GX 1+4 is the best candidate for the detection of glitches with currently operating X-ray instruments and future missions such as the proposed Large Observatory for X-ray Timing (LOFT).

The high density interior of a neutron star is expected to contain superconducting protons and superfluid neutrons. Theoretical estimates suggest that the protons will form a type II superconductor in which the stellar magnetic field is carried by flux tubes. The strong interaction between the flux tubes and the neutron rotational vortices could lead to strong 'pinning', i.e. vortex motion could be impeded. This has important implications especially for pulsar glitch models as it would lead to a large part of the vorticity of the star being decoupled from the 'normal' component, to which the electromagnetic emission is locked. In this paper we explore the consequences of strong pinning in the core on the 'snowplow' model for pulsar glitches (Pizzochero 2011), making use of realistic equations of state and relativistic background models for the neutron star. We find that in general a large fraction of pinned vorticity in the core is not compatible with observations of giant glitches in the Vela pulsar. The conclusion is thus that either most of the core is in a type I superconducting state or that the interaction between vortices and flux tubes is weaker than previously assumed.

Many pulsars are observed to "glitch", i.e. show sudden jumps in their rotational frequency $\nu$, some of which can be as large as $\Delta \nu/\nu\approx 10^{-6}-10^{-5}$ in a subset of pulsars known as giant glitchers. Recently Pizzochero (2011) has shown that an analytic model based on realistic values for the pinning forces in the crust and for the angular momentum transfer in the star can describe the average properties of giant glitches, such as the inter-glitch waiting time, the step in frequency and that in frequency derivative. In this paper we extend the model (originally developed in Newtonian gravity and for a polytropic equation of state) to realistic backgrounds obtained by integrating the relativistic equations of stellar structure and using physically motivated equations of state to describe matter in the neutron star. We find that this more detailed treatment still reproduces the main features of giant glitches in the Vela pulsar and allows us to set constraints on the equation of state. In particular we find that stiffer equations of state are favoured and that it is unlikely that the Vela pulsar has a high mass (larger than $M\approx 1.5 M_\odot$).

The angular momentum transfer associated to Vela-like glitches has never been calculated {\em directly} within a realistic scenario for the storage and release of superfluid vorticity; therefore, the explanation of giant glitches in terms of vortices has not yet been tested against observations. We present the first physically reasonable model, both at the microscopic and macroscopic level (spherical geometry, n=1 polytropic density profile, density-dependent pinning forces compatible with vortex rigidity), to determine where in the star the vorticity is pinned, how much of it, and for how long. For standard neutron star parameters ($M=1.4 M_{\odot}, R_s=10$ km, $\dot{\Omega}=\dot{\Omega}_{\rm Vela}=-10^{-10}$ Hz s$^{-1}$), we find that maximum pinning forces of order $f_m\approx10^{15}$ dyn cm$^{-1}$ can accumulate $\Delta L_{\rm gl}\approx10^{40}$ erg s of superfluid angular momentum, and release it to the crust at intervals $\Delta t_{\rm gl}\approx3$ years. This estimate of $\Delta L_{\rm gl}$ is one order of magnitude smaller than what implied indirectly by current models for post-glitch recovery, where the core and inner-crust vortices are taken as physically disconnected; yet, it successfully yields the magnitudes observed in recent Vela glitches for {\em both} jump parameters, $\Delta\Omega_{\rm gl}$ and $\Delta\dot{\Omega}_{\rm gl}$, provided one assumes that only a small fraction ($<10%$) of the total star vorticity is coupled to the crust on the short timescale of a glitch. This is reasonable in our approach, where no layer of normal matter exists between the core and the inner-crust, as indicated by existing microscopic calculation. The new scenario presented here is nonetheless compatible with current post-glitch models.

Although pulsars are one of the most stable clocks in the universe, many of them are observed to 'glitch', i.e. to suddenly increase their spin frequency (\nu) with fractional increases that range from \Delta\nu/\nu \approx 10^{-11} to 10^{-5}. In this paper we focus on the 'giant' glitches, i.e. glitches with fractional increases in the spin rate of the order of \Delta\nu/{\nu} \approx 10^{-6}, that are observed in a sub class of pulsars including the Vela. We show that giant glitches can be modelled with a two-fluid hydrodynamical approach. The model is based on the formalism for superfluid neutron stars of Andersson and Comer (2006) and on the realistic pinning forces of Grill and Pizzochero (2011). We show that all stages of Vela glitches, from the rise to the post-glitch relaxation, can be reproduced with a set of physically reasonable parameters and that the sizes and waiting times between giant glitches in other pulsars are also consistent with our model.

In this lecture, we give a first introduction to neutron stars, based on fundamental physical principles. After outlining their amazing macroscopic properties, as obtained from observations, we infer the extreme conditions of matter in their interiors. We then describe two crucial physical phenomena which characterize compact stars, gravitational stability of strongly degenerate matter and neutronization of nuclear matter with increasing density, and explain how the formation and properties of neutron stars are a consequence of the extreme compression of matter under gravity. Finally, we describe how astronomical observations of various external macroscopic features can give invaluable information about the exotic microscopic scenario inside: neutrons stars represent a unique probe to study super-dense, isospin-asymmetric, superfluid, bulk hadronic matter.

We investigate to what extent the temperature dependence of the nuclear symmetry energy can affect the neutronization of the stellar core prior to neutrino trapping during gravitational collapse. To this end, we implement a one-zone simulation to follow the collapse until beta equilibrium is reached and the lepton fraction remains constant. Since the strength of electron capture on the neutron-rich nuclei associated to the supernova scenario is still an open issue, we keep it as a free parameter. We find that the temperature dependence of the symmetry energy consistently yields a small reduction of deleptonization, which corresponds to a systematic effect on the shock wave energetics: the gain in dissociation energy of the shock has a small yet non-negligible value of about 0.4 foe (1 foe = 10^51 erg) and this result is almost independent from the strength of nuclear electron capture. The presence of such a systematic effect and its robustness under changes of the parameters of the one-zone model are significative enough to justify further investigations with detailed numerical simulations of supernova explosions.

We performed a detailed study of the 2007 outburst of the 352s pulsar SAXJ2103.5+4545, a Be/X-ray transient observed by INTEGRAL, to study its spectral and temporal properties during the evolution of the outburst. SAXJ2103.5+4545 was observed with IBIS/ISGRI from 25 to 27 April 2007 and from 6 to 8 May 2007. The 20-100keV spectrum is well described by a bremsstrahlung model with a temperature kT = 24keV. The pulse profiles are variable with time and energy. A pulse period derivative of pdot = -3.4E-7 s/s has been observed during the outburst. Instead, a spin-down of pdot = 5.5E-9 s/s is observed between the 2007 outburst reported here and the previous one occurred in December 2004. This is the largest spin-down measured for SAXJ2103.5+4545 since its discovery. We estimate a neutron star magnetic field in the range (1.6-3)E13 G using the Ghosh & Lamb torque model.

We investigate when the energy that pins a superfluid vortex to the lattice of nuclei in the inner crust of neutron stars can be approximated by the energy that binds the vortex to a single nucleus. Indeed, although the pinning energy is the quantity relevant to the theory of pulsar glitches, so far full quantum calculations have been possible only for the binding energy. Physically, the presence of nearby nuclei can be neglected if the lattice is dilute, namely with nuclei sufficiently distant from each other. We find that the dilute limit is reached only for quite large Wigner-Seitz cells, with radii > 55 fm; these are found only in the outermost low-density regions of the inner crust. We conclude that present quantum calculations do not correspond to the pinning energies in almost the entire inner crust and thus their results are not predictive for the theory of glitches.

We calculate the force that pins vortices in the neutron superfluid to nuclei in the inner crust of rotating neutron stars, relying on a detailed microscopic description of both the vortex radial profile and the inner crust nuclear structure. The contribution to the pinning energy from pair condensation is estimated in the local density approximation with realistic nucleon-nucleon interactions. The kinetic contribution, not consistently included in previous approaches, is evaluated in the same approximation, and found to be relevant. The vortex-nucleus interaction turns out to be attractive for stellar desities greater than $\sim 10^13 g/cm^3$. In this region, we find values for the pinning force which are almost one order of magnitude lower than the ones obtained so far. This has direct consequences on the critical velocity differences for depinning.