The dynamics of dissipative multi-fluid neutron star cores

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.Comment: RevTex, no figure

The Spin Distribution of Fast Spinning Neutron Stars in Low Mass X-Ray Binaries: Evidence for Two Sub-Populations

We study the current sample of rapidly rotating neutron stars in both accreting and non-accreting binaries in order to determine whether the spin distribution of accreting neutron stars in low-mass X-ray binaries can be reconciled with current accretion torque models. We perform a statistical analysis of the spin distributions and show that there is evidence for two sub-populations among low-mass X-ray binaries, one at relatively low spin frequency, with an average of ~300 Hz and a broad spread, and a peaked population at higher frequency with average spin frequency of ~575 Hz. We show that the two sub-populations are separated by a cut-point at a frequency of ~540 Hz. We also show that the spin frequency of radio millisecond pulsars does not follow a log-normal distribution and shows no evidence for the existence of distinct sub-populations. We discuss the uncertainties of different accretion models and speculate that either the accreting neutron star cut-point marks the onset of gravitational waves as an efficient mechanism to remove angular momentum or some of the neutron stars in the fast sub-population do not evolve into radio millisecond pulsars.Comment: Submitted to Ap

Relativistic Two-stream Instability

We study the (local) propagation of plane waves in a relativistic, non-dissipative, two-fluid system, allowing for a relative velocity in the "background" configuration. The main aim is to analyze relativistic two-stream instability. This instability requires a relative flow -- either across an interface or when two or more fluids interpenetrate -- and can be triggered, for example, when one-dimensional plane-waves appear to be left-moving with respect to one fluid, but right-moving with respect to another. The dispersion relation of the two-fluid system is studied for different two-fluid equations of state: (i) the "free" (where there is no direct coupling between the fluid densities), (ii) coupled, and (iii) entrained (where the fluid momenta are linear combinations of the velocities) cases are considered in a frame-independent fashion (eg. no restriction to the rest-frame of either fluid). As a by-product of our analysis we determine the necessary conditions for a two-fluid system to be causal and absolutely stable and establish a new constraint on the entrainment.Comment: 15 pages, 2 eps-figure

A Relativistic Mean Field Model for Entrainment in General Relativistic Superfluid Neutron Stars

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.Comment: 16 pages, 5 figures, submitted to PR

The dynamics of neutron star crusts: Lagrangian perturbation theory for a relativistic superfluid-elastic system

The inner crust of a mature neutron star is composed of an elastic lattice of neutron-rich nuclei penetrated by free neutrons. These neutrons can flow relative to the crust once the star cools below the superfluid transition temperature. In order to model the dynamics of this system, which is relevant for a range of problems from pulsar glitches to magnetar seismology and continuous gravitational-wave emission from rotating deformed neutron stars, we need to understand general relativistic Lagrangian perturbation theory for elastic matter coupled to a superfluid component. This paper develops the relevant formalism to the level required for astrophysical applications.Comment: 31 pages, double spacing, minor typos fixe

Lagrangian perturbation theory for a superfluid immersed in an elastic neutron star crust

The inner crust of mature neutron stars, where an elastic lattice of neutron-rich nuclei coexists with a neutron superfluid, impacts on a range of astrophysical phenomena. The presence of the superfluid is key to our understanding of pulsar glitches, and is expected to affect the thermal conductivity and hence the evolution of the surface temperature. The coupling between crust and superfluid must also be accounted for in studies of neutron star dynamics, discussions of global oscillations and associated instabilities. In this paper we develop Lagrangian perturbation theory for this problem, paying attention to key issues like superfluid entrainment, potential vortex pinning, dissipative mutual friction and the star's magnetic field. We also discuss the nature of the core-crust interface. The results provide a theoretical foundation for a range of interesting astrophysical applications.Comment: 13 pages, no figures, to appear in MNRA

R-mode oscillations and rocket effect in rotating superfluid neutron stars. I. Formalism

We derive the hydrodynamical equations of r-mode oscillations in neutron stars in presence of a novel damping mechanism related to particle number changing processes. The change in the number densities of the various species leads to new dissipative terms in the equations which are responsible of the {\it rocket effect}. We employ a two-fluid model, with one fluid consisting of the charged components, while the second fluid consists of superfluid neutrons. We consider two different kind of r-mode oscillations, one associated with comoving displacements, and the second one associated with countermoving, out of phase, displacements.Comment: 10 page

A detailed study of quasinormal frequencies of the Kerr black hole

We compute the quasinormal frequencies of the Kerr black hole using a continued fraction method. The continued fraction method first proposed by Leaver is still the only known method stable and accurate for the numerical determination of the Kerr quasinormal frequencies. We numerically obtain not only the slowly but also the rapidly damped quasinormal frequencies and analyze the peculiar behavior of these frequencies at the Kerr limit. We also calculate the algebraically special frequency first identified by Chandrasekhar and confirm that it coincide with the $n=8$ quasinormal frequency only at the Schwarzschild limit.Comment: REVTEX, 15 pages, 7 eps figure

The enigmatic spin evolution of PSR J0537-6910: r-modes, gravitational waves and the case for continued timing

We discuss the unique spin evolution of the young X-ray pulsar PSR J0537-6910, a system in which the regular spin down is interrupted by glitches every few months. Drawing on the complete timing data from the Rossi X-ray Timing Explorer (RXTE, from 1999-2011), we argue that a trend in the inter-glitch behaviour points to an effective braking index close to $n=7$, much larger than expected. This value is interesting because it would accord with the neutron star spinning down due to gravitational waves from an unstable r-mode. We discuss to what extent this, admittedly speculative, scenario may be consistent and if the associated gravitational-wave signal would be within reach of ground based detectors. Our estimates suggest that one may, indeed, be able to use future observations to test the idea. Further precision timing would help enhance the achievable sensitivity and we advocate a joint observing campaign between the Neutron Star Interior Composition ExploreR (NICER) and the LIGO-Virgo network.Comment: 10 pages, 4 figures, emulate ApJ forma

Slowly Rotating General Relativistic Superfluid Neutron Stars with Relativistic Entrainment

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.Comment: 30 pages, 10 figures, uses ReVTeX