The first-order relativistic fluid theories of dissipation proposed by Eckart and Landau-Lifshitz have been proved to be unstable. They admit solutions which start in proximity of equilibrium and depart exponentially from it. We show that this behaviour is due to the fact that the total entropy of these fluids, restricted to the dynamically accessible states, has no upper bound. As a result, these systems have the tendency to constantly change according to the second law of thermodynamics and the unstable modes represent the directions of growth of the entropy in state space. We, then, verify that the conditions of stability of Israel and Stewart's theory are exactly the requirements for the entropy to have an absolute maximum. Hence, we explain how the instability of the first-order theories is a direct consequence of the truncation of the entropy current at the first order, which turns the maximum into a saddle point of the total entropy. Finally, we show that recently proposed first-order stable theories, constructed using more general frames, do not solve the instability problem by providing a maximum for the entropy, but, rather, are made stable by allowing for small violations of the second law.

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.

Bulk viscosity plays a key role in several astrophysical phenomena linked to neutron stars, from the damping of modes of oscillation in mature stars to the dynamics of the hot remnant of a binary merger. To quantitatively model these systems and constrain the imprint of the equation of state of dense matter it is, however, necessary to resort to a general relativistic description of bulk viscosity. Several prescriptions exist to modify the standard Navier-Stokes equations, which are parabolic and lead to non-causal solutions, to obtain hyperbolic equations that can be used in a relativistic setting. The connection between microphysics and the parameters introduced is not, however, always transparent. In this paper we present a relativistic formalism where any thermodynamic process which contributes to the bulk viscosity is modelled as a set of chemical reactions, whose reaction coordinates are abstract parameters describing the displacement from local thermodynamic equilibrium. The result is a non-equilibrium thermodynamic theory for bulk-viscous fluids which does not rely on any near-equilibrium assumption, and naturally gives rise to telegraph-type equations for the reaction processes. We present also two concrete applications of our theory, to the bulk viscosity in neutron stars and to a non interacting neutron gas and compare our formulation with that of Israel and Stewart (1979), which for a mature neutron star emerges as a perturbative expansion of our theory near equilibrium.

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.

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.

The internal layers of neutron stars are expected to contain several superfluid components that can significantly affect their dynamics. The description of such objects should rely on hydrodynamic models in which it is possible to unambiguously assign the value of the thermodynamic variables from microscopic calculations of the properties of matter. In this work we consider the phenomenological approach to multifluids modelling championed by Carter and, studying the relaxation of the system towards equilibrium, we assign a precise thermodynamic interpretation to its variables. We show that in thermodynamic equilibrium the equation of state contains less state variables than those needed in the phenomenological model, implying the existence of a gauge freedom of the theory that can be used to simplify the hydrodynamic formulation in the non-dissipative limit. Once this is understood, it becomes easy to translate the different multifluid formalisms that have been proposed in the literature into Carter's form. Finally, we show that the usual concepts of affinity and reaction coordinates, as they are introduced in chemistry, are not affected by the presence of superfluid currents. In an effort to make the concepts clear, the formalism is developed step-by-step from first principles, providing model examples and several applications of practical relevance for the study of superfluid neutron star interiors.

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.

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 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.

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 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.

From the study of a sample of about 62.3 million well reconstructed K0S decays recorded by the KLOE detector at the DAFNE accelerator in Frascati, the lifetimes of K0S mesons parallel and antiparallel to the direction of motion of the Earth with respect to the Cosmic Microwave Background reference frame have been studied. No difference has been found, and a limit on a possible asymmetry of the lifetime with respect to the CMB has been set at 95% C.L.: A < 0.98 x 10-3. This is presently the best experimental limit on such quantity, and it is smaller of the speed, expressed in natural units, of the Solar System with respect to the CMB. The present limit might constrain possible Lorentz-violating anisotropical theories.