We consider a resistive multi-fluid framework from the 3+1 space-time foliation point-of-view, paying particular attention to issues relating to the use of multi-parameter equations of state and the associated inversion from evolved to primitive variables. We highlight relevant numerical issues that arise for general systems with relative flows. As an application of the new formulation, we consider a three-component system relevant for hot neutron stars. In this case we let the baryons (neutrons and protons) move together, but allow heat and electrons to exhibit relative flow. This reduces the problem to three momentum equations; overall energy-momentum conservation, a generalised Ohm's law and a heat equation. Our results provide a hierarchy of increasingly complex models and prepare the ground for new state-of-the-art simulations of relevant scenarios in relativistic astrophysics.

We develop an action principle to construct the field equations for a multi-fluid system containing charge-neutral fluids, plasmas, and dissipation (via resistive interactions), by combining the standard, Maxwell action and minimal coupling of the electromagnetic field with a recently developed action for relativistic dissipative fluids. We use a pull-back formalism from spacetime to abstract matter spaces to build unconstrained variations for both the charge-neutral fluids and currents making up the plasmas. Using basic linear algebra techniques, we show that a general "relabeling" invariance exists for the abstract matter spaces. With the field equations in place, a phenomenological model for the resistivity is developed, using as constraints charge conservation and the Second Law of Thermodynamics. A minimal model for a system of electrons, protons, and heat is developed using the Onsager procedure for incorporating dissipation.

We develop a new framework for the modelling of charged fluid dynamics in general relativity. The model, which builds on a recently developed variational multi-fluid model for dissipative fluids, accounts for relevant effects like the inertia of both charge currents and heat and, for mature systems, the decoupling of superfluid components. We discuss how the model compares to standard relativistic magnetohydronamics and consider the connection between the fluid dynamics, the microphysics and the underlying equation of state. As illustrations of the formalism, we consider three distinct two-fluid models describing i) an Ohm's law for resistive charged flows, ii) a relativistic heat equation, and iii) an equation representing the momentum of a decoupled superfluid component. As a more complex example, we also formulate a three-fluid model which demonstrates the thermo-electric effect. This framework allows us to model neutron stars (and related systems) at a hierarchy of increasingly complex levels, and should enable us to make progress on a range of exciting problems in astrophysics and cosmology.

Building on a recently improved understanding of the problem of heat flow in general relativity, we develop a hydrodynamical model for coupled finite temperature superfluids. The formalism is designed with the dynamics of the outer core of a mature neutron star (where superfluid neutrons are coupled to a conglomerate of protons and electrons) in mind, but the main ingredients are relevant for a range of analogous problems. The entrainment between material fluid components (the condensates) and the entropy (the thermal excitations) plays a central role in the development. We compare and contrast the new model to previous results in the literature, and provide estimates for the relevant entrainment coefficients that should prove useful in future applications. Finally, we consider the sound-wave propagation in the system in two simple limits, demonstrating the presence of second sound if the temperature is sub-critical, but absence of this phenomenon above the critical temperature for superfluidity.

Two-stream instability requires, essentially, two things to operate: a relative flow between two fluids and some type of interaction between them. In this letter we provide the first demonstration that this mechanism may be active in a cosmological context. Building on a recently developed formalism for cosmological models with two, interpenetrating fluids with a relative flow between them, we show that two-stream instability may be triggered during the transition from one fluid domination to the other. We also demonstrate that the cosmological expansion eventually shuts down the instability by driving to zero the relativeflow and the coupling between the two fluids.

We present a Newtonian multi-fluid formalism for superfluid neutron star cores, focussing on the additional dissipative terms that arise when one takes into account the individual dynamical degrees of freedom associated with the coupled "fluids". The problem is of direct astrophysical interest as the nature of the dissipative terms can have significant impact on the damping of the various oscillation modes of the star and the associated gravitational-wave signatures. A particularly interesting application concerns the gravitational-wave driven instability of f- and r-modes. We apply the developed formalism to two specific three-fluid systems: (i) a hyperon core in which both Lambda and Sigma^- hyperons are present, and (ii) a core of deconfined quarks in the colour-flavour-locked phase in which a population of neutral K^0 kaons is present. The formalism is, however, general and can be applied to other problems in neutron-star dynamics (such as the effect of thermal excitations close to the superfluid transition temperature) as well as laboratory multi-fluid systems.

We present the first multi-fluid analysis of a dense neutron star core with a deconfined colour-flavour-locked superconducting quark component. Accounting only for the condensate and (finite temperature) phonons, we make progress by taking over results for superfluid $^4$He. The resultant two-fluid model accounts for a number of additional viscosity coefficients (compared to the Navier-Stokes equations) and we show how they enter the dissipation analysis for an oscillating star. We provide simple estimates for the gravitational-wave driven r-mode instability, demonstrating that the various phonon processes that we consider are not effective damping agents. Even though the results are likely of little direct astrophysical importance (since we consider an overly simplistic stellar model) our analysis represents significant technical progress, laying the foundation for more detailed numerical studies and preparing the ground for the inclusion of additional aspects (in particular associated with kaons) of the problem.

We develop the formalism for determining the quasinormal modes of general relativistic multi-fluid compact stars in such a way that the impact of superfluid gap data can be assessed. Our results represent the first attempt to study true multi-layer dynamics, an important step towards considering realistic superfluid/superconducting compact stars. We combine a relativistic model for entrainment with model equations of state that explicity incorporate the symmetry energy. Our analysis emphasises the many different parameters that are required for this kind of modelling, and the fact that standard tabulated equations of state are grossly incomplete in this respect. To make progress, future equations of state need to provide the energy density as a function of the various nucleon number densities, the temperature (i.e. entropy), and the entrainment among the various components.

We discuss wave propagation in rotating superfluid neutron star cores, taking into account the vortex mediated mutual friction force. For models where the two fluids co-rotate in the unperturbed state, our analysis clarifies the role of chemical coupling and entrainment for sound and inertial waves. We also investigate the mutual friction damping, providing results that demonstrate the well-known fact that sound waves propagating along a vortex array are undamped. We show that the same is not true for inertial waves, which are damped by the mutual friction regardless of the propagation direction. We then include the vortex tension, which arises due to local vortex curvature. Focussing on purely transverse inertial waves, we derive the small correction that the tension induces in the wave frequency. Finally, we allow for a relative linear flow in the background (along the rotation axis). In this case we show how the mutual friction coupling may induce a dynamical instability in the inertial waves. We discuss the critical flow required for the instability to be present, its physical interpretation and the possible relevance it may have for neutron star physics.

We analyse the implications of superfluid turbulence for neutron star physics. We begin by extending our previous results for the mutual friction force for a straight vortex array to account for the self-induced flow which arises when the vortices are curved. We then discuss Vinen's phenomenological model for isotropic turbulence, and derive the associated (Gorter-Mellink) form for the mutual friction. We compare this derivation to a more recent analysis of Schwarz, which sheds light on various involved issues. Having discussed isotropic turbulence, we argue that this case is unlikely to be relevant for neutron stars. Instead we expect a rotating neutron star to exhibit polarised turbulence, where relative flow drives the turbulence and rotation counteracts it. Based on recent results for superfluid Helium, we construct a phenomenological model that should have the key features of such a polarised turbulent system.

We discuss vortex-mediated mutual friction in the two-fluid model for superfluid neutron star cores. Our discussion is based on the general formalism developed by Carter and collaborators, which makes due distinction between transport velocity and momentum for each fluid. This is essential for an implementation of the so-called entrainment effect, whereby the flow of one fluid imparts momentum in the other and vice versa. The mutual friction follows by balancing the Magnus force that acts on the quantised neutron vortices with a resistive force due to the scattering of electrons off of the magnetic field with which each vortex core is endowed. We derive the form of the macroscopic mutual friction force which is relevant for a model based on smooth-averaging over a collection of vortices. We discuss the coefficients that enter the expression for this force, and the timescale on which the two interpenetrating fluids in a neutron star core are coupled. This discussion confirms that our new formulation accords well with previous work in this area.

We develop a flux-conservative formalism for a Newtonian multi-fluid system, including dissipation and entrainment (i.e. allowing the momentum of one fluid to be a linear combination of the velocities of all fluids). Maximum use is made of mass, energy, and linear and angular momentum conservation to specify the equations of motion. Also used extensively are insights gleaned from a convective variational action principle, key being the distinction between each velocity and its canonically conjugate momentum. Dissipation is incorporated to second order in the ``thermodynamic forces'' via the approach pioneered by Onsager. An immediate goal of the investigation is to understand better the number, and form, of independent dissipation terms required for a consistent set of equations of motion in the multi-fluid context. A significant, but seemingly innocuous detail, is that one must be careful to isolate ``forces'' that can be written as total gradients, otherwise errors can be made in relating the net internal force to the net externally applied force. Our long-range aim is to provide a formalism that can be used to model dynamical multi-fluid systems both perturbatively and via fully nonlinear 3D numerical evolutions. To elucidate the formalism we consider the standard model for a heat-conducting, superfluid neutron star, which is believed to be dominated by superfluid neutrons, superconducting protons, and a highly degenerate, ultra-relativistic gas of normal fluid electrons. We determine that in this case there are, in principle, 19 dissipation coefficients in the final set of equations.

We discuss the effects of superfluidity on the shear viscosity in a neutron star core. Our study combines existing theoretical results for the viscosity coefficients with data for the various superfluid energy gaps into a consistent description. In particular, we provide a simple model for the electron viscosity which is relevant both when the protons form a normal fluid and when they become superconducting. This model explains in a clear way why proton superconductivity leads to a significant strengthening of the shear viscosity. We present our results in a form which permits the use of data for any given modern equation of state (our final formulas are explicitly dependent on the proton fraction). We discuss a simple description of the relevant superfluid pairing gaps, and construct a number of models (spanning the range of current uncertainty) which are then used to discuss the superfluid suppression of shear viscosity. We conclude by a discussion of a number of challenges that must be met if we are to make further progress in this area of research.

We have developed a theoretical model and a numerical code for stationary rotating superfluid neutron stars in full general relativity. The underlying two-fluid model is based on Carter's covariant multi-fluid hydrodynamic formalism. The two fluids, representing the superfluid neutrons on one hand, and the protons and electrons on the other, are restricted to uniform rotation around a common axis, but are allowed to have different rotation rates. We have performed extensive tests of the numerical code, including quantitative comparisons to previous approximative results for these models. The results presented here are the first ``exact'' calculations of such models in the sense that no approximations (other than that inherent in a discretized numerical treatment) are used. Using this code we reconfirm the existence of prolate-oblate shaped configurations. We studied the dependency of the Kepler rotation limit and of the mass-density relation on the relative rotation rate. We further demonstrate how one can simulate a (albeit fluid) neutron-star ``crust'' by letting one fluid extend further outwards than the other, which results in interesting cases where the Kepler limit is actually determined by the outermost but slower fluid.

We develop a Lagrangian perturbation framework for rotating non-relativistic superfluid neutron stars. This leads to the first generalization of classic work on the stability properties of rotating stars to models which account for the presence of potentially weakly coupled superfluid components. Our analysis is based on the standard two-fluid model expected to be relevant for the conditions that prevail in the outer core of mature neutron stars. We discuss the implications of our results for dynamical and secular instabilities of a simple neutron star model in which the two fluids are allowed to assume different (uniform) rotation rates.

Neutron stars that are cold enough should have two or more superfluids/supercondutors in their inner crusts and cores. The implication of superfluidity/superconductivity for equilibrium and dynamical neutron star states is that each individual particle species that forms a condensate must have its own, independent number density current and equation of motion that determines that current. An important consequence of the quasiparticle nature of each condensate is the so-called entrainment effect, i.e. the momentum of a condensate is a linear combination of its own current and those of the other condensates. We present here the first fully relativistic modelling of slowly rotating superfluid neutron stars with entrainment that is accurate to the second-order in the rotation rates. The stars consist of superfluid neutrons, superconducting protons, and a highly degenerate, relativistic gas of electrons. We use a relativistic $\sigma$ - $\omega$ mean field model for the equation of state of the matter and the entrainment. We determine the effect of a relative rotation between the neutrons and protons on a star's total mass, shape, and Kepler, mass-shedding limit.

We present results concerning adiabatic inertial-mode oscillations of non-stratified superfluid neutron stars in Newtonian gravity, using the anelastic and slow-rotation approximations. We consider a simple two-fluid model of a superfluid neutron star, where one fluid consists of the superfluid neutrons and the second fluid contains all the comoving constituents (protons, electrons). The two fluids are assumed to be ``free'' in the sense that vortex-mediated forces like mutual friction or pinning are absent, but they can be coupled by the equation of state, in particular by entrainment. The stationary background consists of the two fluids rotating uniformly around the same axis with potentially different rotation rates. We study the special cases of co-rotating backgrounds, vanishing entrainment, and the purely toroidal r-modes, analytically. We calculate numerically the eigenfunctions and frequencies of inertial modes in the general case of non co-rotating backgrounds, and study their dependence on the relative rotation rate and entrainment. In these non-stratified models we find avoided crossings only between associated mode-pairs, eg. an ``ordinary'' mode and its ``superfluid'' counterpart, while other mode-frequencies generally cross as the background parameters are varied. We confirm (for the first time in a mode calculation) the onset of a ``two-stream instability'' at a critical relative background rotation rate, and study some of the properties of this instability for the inertial modes.

General relativistic superfluid neutron stars have a significantly more intricate dynamics than their ordinary fluid counterparts. Superfluidity allows different superfluid (and superconducting) species of particles to have independent fluid flows, a consequence of which is that the fluid equations of motion contain as many fluid element velocities as superfluid species. Whenever the particles of one superfluid interact with those of another, the momentum of each superfluid will be a linear combination of both superfluid velocities. This leads to the so-called entrainment effect whereby the motion of one superfluid will induce a momentum in the other superfluid. We have constructed a fully relativistic model for entrainment between superfluid neutrons and superconducting protons using a relativistic $\sigma - \omega$ mean field model for the nucleons and their interactions. In this context there are two notions of ``relativistic'': relativistic motion of the individual nucleons with respect to a local region of the star (i.e. a fluid element containing, say, an Avogadro's number of particles), and the motion of fluid elements with respect to the rest of the star. While it is the case that the fluid elements will typically maintain average speeds at a fraction of that of light, the supranuclear densities in the core of a neutron star can make the nucleons themselves have quite high average speeds within each fluid element. The formalism is applied to the problem of slowly-rotating superfluid neutron star configurations, a distinguishing characteristic being that the neutrons can rotate at a rate different from that of the protons.

We describe recent progress in the numerical study of the structure of rapidly rotating superfluid neutron star models in full general relativity. The superfluid neutron star is described by a model of two interpenetrating and interacting fluids, one representing the superfluid neutrons and the second consisting of the remaining charged particles (protons, electrons, muons). We consider general stationary configurations where the two fluids can have different rotation rates around a common rotation axis. The previously discovered existence of configurations with one fluid in a prolate shape is confirmed.

This paper provides the first study of a new dynamical instability in superfluids. This instability is similar to the two-stream instability known to operate in plasmas. It is analogous to the Kelvin-Helmholtz instability, but has the distinguishing feature that the two fluids are interpenetrating. The instability sets in once the relative flow between the two components of the system reaches a critical level. Our analysis is based on the two-fluid equations that have been used to model the dynamics of the outer core of a neutron star, where superfluid neutrons are expected to coexist with superconducting protons and relativistic electrons. These equations are analogous to the standard Landau model for superfluid Helium. We study this instability for two different model problems. First we analyze a local dispersion relation for waves in a system where one fluid is at rest while the other flows at a constant rate. This provides a proof of principle of the existence of the two-stream instability for superfluids. Our second model problem concerns two rotating fluids confined within an infinitesimally thin spherical shell. The aim of this model is to assess whether the two-stream instability may be relevant (perhaps as a trigger mechanism) for pulsar glitches. Our results for this problem show that the entrainment effect could provide a sufficiently strong coupling for the instability to set in at a relative flow small enough to be astrophysically plausible.

Mature neutron stars are expected to have several superfluid components. Strong evidence for this is provided by the glitches that have been observed in dozens of pulsars. The underlying idea behind most glitch models is that, as the neutron star crust spins down due to the emission of electromagnetic radiation, the superfluid component lags behind until a critical point is reached and angular momentum is transferred from the superfluid to the crust, leading to the spin-up associated with the glitch. In this Letter we describe a superfluid analogue of the two-stream instability that is well known in plasma physics, and provide evidence that this instability is likely to be relevant for neutron stars. This is a new physical mechanism which may play a key role in explaining the glitch mechanism and which could also prove to be relevant in laboratory experiments on superfluid Helium.

Isolated neutron stars undergoing non-radial oscillations are expected to emit gravitational waves in the kilohertz frequency range. To date, radio astronomers have located about 1,300 pulsars, and can estimate that there are about 2 times 10^8 neutron stars in the galaxy. Many of these are surely old and cold enough that their interiors will contain matter in the superfluid or superconducting state. In fact, the so-called glitch phenomenon in pulsars (a sudden spin-up of the pulsar's crust) is best described by assuming the presence of superfluid neutrons and superconducting protons in the inner crusts and cores of the pulsars. Recently there has been much progress on modelling the dynamics of superfluid neutron stars in both the Newtonian and general relativistic regimes. We will discuss some of the main results of this recent work, perhaps the most important being that superfluidity should affect the gravitational waves from neutron stars (emitted, for instance, during a glitch) by modifying both the rotational properties of the background star and the modes of oscillation of the perturbed configuration. Finally, we present an analysis of the so-called zero-frequency subspace (i.e. the space of time-independent perturbations) and determine that it is spanned by two sets of polar (or spheroidal) and two sets of axial (or toroidal) degenerate perturbations for the general relativistic system. As in the Newtonian case, the polar perturbations are the g-modes which are missing from the pulsation spectrum of a non-rotating configuration, and the axial perturbations should lead to two sets of r-modes when the degeneracy of the frequencies is broken by having the background rotate.

We develop a general formalism to treat, in general relativity, the nonradial oscillations of a superfluid neutron star about static (non-rotating) configurations. The matter content of these stars can, as a first approximation, be described by a two-fluid model: one fluid is the neutron superfluid, which is believed to exist in the core and inner crust of mature neutron stars; the other fluid is a conglomerate of all charged constituents. We use a system of equations that governs the perturbations both of the metric and of the matter variables, whatever the equation of state for the two fluids. The entrainment effect is explicitly included. We also allow for an outer envelope composed of ordinary fluid. We derive and implement the junction conditions for the metric and matter variables at the core/envelope interface. We investigate how the quasinormal modes of a superfluid star are affected by changes in the entrainment parameter, and unveil a series of avoided crossings between the various modes. We provide a proof that all modes of a two-fluid star must radiate gravitationally. We also discuss the future detectability of pulsations in a superfluid star and argue that it may be possible to use gravitational-wave data to constrain the parameters of superfluid neutron stars.

We discuss the possibility that future gravitational-wave detectors may be able to detect various modes of oscillation of old, cold neutron stars. We argue that such detections would provide unique insights into the superfluid nature of neutron star cores, and could also lead to a much improved understanding of pulsar glitches. Our estimates are based on a detector configuration with several narrowbanded (cryogenic) interferometers operating as a "xylophone" which could lead to high sensitivity at high frequencies. We also draw on recent advances in our understanding of the dynamics of pulsating superfluid neutron star cores.

We discuss the nature of the various modes of pulsation of superfluid neutron stars using comparatively simple Newtonian models. The matter in these stars is described in terms of a two-fluid model, where one fluid is the neutron superfluid, which is believed to exist in the core and inner crust of mature neutron stars, and the other fluid represents a conglomerate of all other constituents (crust nuclei, protons, electrons, etc.). In our model, we incorporate the non-dissipative interaction due to the so-called entrainment effect, whereby the momentum of one constituent (e.g. the neutrons) carries along part of the mass of the other constituent. We show that there is no independent set of propagating g-modes in a superfluid neutron star core. We also approach the toroidal pulsation modes of a slowly rotating superfluid star. We show that the superfluid equations support a new class of r-modes. Finally, the role of the entrainment effect on the superfluid mode frequencies is shown explicitly via solutions to local dispersion relations.

We develop a formalism to model slowly rotating superfluid Newtonian neutron stars. A two-fluid model is used to describe the matter, where one fluid is the superfluid neutrons that are believed to exist in the inner crust and core of mature neutron stars, while the other fluid is a charge neutral conglomerate of the remaining constituents (crust nuclei, core superconducting protons, electrons, etc). We include the entrainment effect, which is a non-dissipative interaction between the two fluids whereby a momentum induced in one of the fluids will cause part of the mass of the other fluid to be carried along. The equations that describe rotational equilibria (i.e. axisymmetric and stationary configurations) are approximated using the slow-rotation approximation; an expansion in terms of the rotation rates of the two fluids where only terms up to second-order are kept. Our formalism allows the neutrons to rotate at a rate different from that of the charged constituents. For a particular equation of state that is quadratic in the two mass-densities and relative velocities of the fluids, we find an analytic solution to the slow-rotation equations. This result provides an elegant generalisation to the two-fluid problem of the standard expressions for the one-fluid polytrope ${\cal E} \propto \rho^2$. The model equation of state includes entrainment and is general enough to allow for realistic values for, say, mass and radius of the star. It also includes a term that can be related to ``symmetry energy'' terms that appear in more realistic equations of state. We use the analytic solution to explore how relative rotation between the two fluids, the ``symmetry energy'' term, and entrainment affect the neutron star's local distribution of particles, as well as global quantities as the Kepler limit, ellipticity, and moments of inertia.

We present a general formalism to treat slowly rotating general relativistic superfluid neutron stars. As a first approximation, their matter content can be described in terms of a two-fluid model, where one fluid is the neutron superfluid, which is believed to exist in the core and inner crust of mature neutron stars, and the other fluid represents a conglomerate of all other constituents (crust nuclei, protons, electrons, etc.). We obtain a system of equations, good to second-order in the rotational velocities, that determines the metric and the matter variables, irrespective of the equation of state for the two fluids. In particular, allowance is made for the so-called entrainment effect, whereby the momentum of one constituent (e.g. the neutrons) carries along part of the mass of the other constituent. As an illustration of the developed framework, we consider a simplified equation of state for which the two fluids are described by different polytropes. We determine numerically the effects of the two fluids on the rotational frame-dragging, the induced changes in the neutron and proton densities and the inertial mass, as well as the change in shape of the star. We further discuss issues regarding conservation of the two baryon numbers, the mass-shedding (Kepler) limit and chemical equilibrium.

We develop a general formalism to treat, in general relativity, the linear oscillations of a two-fluid star about static (non-rotating) configurations. Such a formalism is intended for neutron stars, whose matter content can be described, as a first approximation, by a two-fluid model: one fluid is the neutron superfluid, which is believed to exist in the core and inner crust of mature neutron stars; the other fluid is a conglomerate of all other constituents (crust nuclei, protons, electrons, etc...). We obtain a system of equations which govern the perturbations both of the metric and of the matter variables, whatever the equation of state for the two fluids. As a first application, we consider the simplified case of two non-interacting fluids, each with a polytropic equation of state. We compute numerically the quasi-normal modes (i.e. oscillations with purely outgoing gravitational radiation) of the corresponding system. When the adiabatic indices of the two fluids are different, we observe a splitting for each frequency of the analogous single fluid spectrum. The analysis also substantiates the claim that w-modes are largely due to spacetime oscillations.