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.

We critically revisit the definition of thermal equilibrium, in its operational formulation, provided by standard thermodynamics. We show that it refers to experimental conditions which break the covariance of the theory at a fundamental level and that, therefore, it cannot be applied to the case of moving bodies. We propose an extension of this definition which is manifestly covariant and can be applied to the study of isolated systems in special relativity. The zeroth law of thermodynamics is, then, proven to establish an equivalence relation among bodies which have not only the same temperature, but also the same center of mass four-velocity.

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.

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.