Search

Next-generation gravitational-wave instruments are expected to constrain the equation of state of dense nuclear matter by observing binaries involving neutron stars. We highlight a problematic systematic error in finite-temperature merger simulations, where shock heating associated with the neutron-star surface gives rise to elevated temperatures. We demonstrate the severe implications of this artificial heating by computing static and dynamical tidal parameters for neutron stars immersed in simulation temperature profiles. The induced systematic errors must be addressed if we want to build robust gravitational-wave signal models for neutron-star, or indeed neutron star-black hole, binaries., Comment: 10 pages, 5 figures, 1 table

We present a proof-of-principle implementation of the first fully covariant filtering scheme applied to relativistic fluid turbulence. The filtering is performed with respect to special observers, identified dynamically as moving with the "bulk of the flow". This means that filtering does not depend on foliations of spacetime but rather on the intrinsic fibration traced out by the observers. The covariance of the approach means that the results may be lifted into an arbitrary, curved spacetime. This practical step follows theoretical work showing that the residuals introduced by filtering a fine-scale ideal fluid can be represented by a non-ideal fluid prescription at the coarse scale. We interpret such non-ideal terms using a simple first-order gradient model, which allows us to extract effective turbulent viscosities and conductivity. A statistical regression on these terms shows that the majority of their variation may be explained by the thermodynamic properties of the filtered fluid and invariants of its flow, such as the shear and vorticity. This serves as a validation of the method and enables us to fit a functional, power-law form for the non-ideal coefficients -- an approach that may be used practically to give a sub-grid closure model in large-eddy simulations., Comment: 23 pages, 18 figures, as published in Phys. Rev. D

Mature neutron stars are expected exhibit gravity g-modes due to stratification caused by varying composition. These modes are affected by nuclear reactions, leading to complex (damped) mode frequencies and the suppression of high order g-modes, in contrast with the common non-dissipative analysis which leads to an infinite g-mode spectrum. Focusing on the transition between the fast and slow reaction regimes, we examine the impact of nuclear reactions on the g-mode spectrum. The general framework for the analysis is presented along with sample numerical results for the BSk21 equation of state and the standard Urca reactions., Comment: 11 pages, 3 figures

The magneto-rotational instability (MRI) - which is due to an interplay between a sheared background and the magnetic field - is commonly considered a key ingredient for developing and sustaining turbulence in the outer envelope of binary neutron star merger remnants. To assess whether (or not) the instability is active and resolved, criteria originally derived in the accretion disk literature - thus exploiting the symmetries of such systems - are often used. In this paper we discuss the magneto-shear instability as a truly local phenomenon, relaxing common symmetry assumptions on the background on top of which the instability grows. This makes the discussion well-suited for highly dynamical environments such as binary mergers. We find that - although this is somewhat hidden in the usual derivation of the MRI dispersion relation - the instability crucially depends on the assumed symmetries. Relaxing the symmetry assumptions on the background we find that the role of the magnetic field is significantly diminished, as it affects the modes' growth but does not drive it. This suggests that we should not expect the standard instability criteria to provide a faithful indication/diagnostic of what "is actually going on" in mergers. We conclude by making contact with a suitable filtering operation, as this is key to separating background and fluctuations in highly dynamical systems.

We re-visit the calculation of mode oscillations in the ocean of a rotating neutron star, which may be excited during thermonuclear X-ray bursts. Our present theoretical understanding of ocean modes relies heavily on the traditional approximation, commonly employed in geophysics. The approximation elegantly decouples the radial and angular sectors of the perturbation problem by neglecting the vertical contribution from the Coriolis force. However, as the implicit assumptions underlying it are not as well understood as they ought to be, we examine the traditional approximation and discuss the associated mode solutions. The results demonstrate that, while the approximation may be appropriate in certain contexts, it may not be accurate for rapidly rotating neutron stars. In addition, using the shallow-water approximation, we show analytically how the solutions that resemble r-modes change their nature in neutron-star oceans to behave like gravity waves. We also outline a simple prescription for lifting Newtonian results in a shallow ocean to general relativity, making the result more realistic., Comment: 16 pages, 2 figures. Missing text in Sec. III discovered in press

In order to extract the precise physical information encoded in the gravitational and electromagnetic signals from powerful neutron-star merger events, we need to include as much of the relevant physics as possible in our numerical simulations. This presents a severe challenge, given that many of the involved parameters are poorly constrained. In this paper we focus on the role of nuclear reactions. Combining a theoretical discussion with an analysis connecting to state-of-the-art simulations, we outline multiple arguments that lead to a reactive system being described in terms of a bulk viscosity. The results demonstrate that in order to properly account for nuclear reactions, future simulations must be able to handle different regimes where rather different assumptions/approximations are appropriate. We also touch upon the link to models based on the large-eddy-strategy required to capture turbulence., Comment: final version published in Phys. Rev. D

The first full angular analysis and an updated measurement of the decay-rate $CP$ asymmetry of the $D^0 \to \pi^+\pi^-\mu^+\mu^-$ and $D^0 \to K^+K^-\mu^+\mu^-$ decays are reported. The analysis uses proton-proton collision data collected with the LHCb detector at centre-of-mass energies of 7, 8 and 13 TeV. The data set corresponds to an integrated luminosity of 9 fb$^{-1}$. The full set of $CP$-averaged angular observables and their $CP$ asymmetries are measured as a function of the dimuon invariant mass. The results are consistent with expectations from the standard model and with $CP$ symmetry., Comment: All figures and tables, along with any supplementary material and additional information, are available at https://cern.ch/lhcbproject/Publications/p/LHCb-PAPER-2021-035.html (LHCb public pages)

Tests of lepton universality in $B^0\to K^0_S \ell^+ \ell^-$ and $B^+\to K^{*+} \ell^+ \ell^-$ decays where $\ell$ is either an electron or a muon are presented. The differential branching fractions of $B^0\to K^0_S e^+ e^-$ and $B^+\to K^{*+} e^+ e^-$ decays are measured in intervals of the dilepton invariant mass squared. The measurements are performed using proton-proton collision data recorded by the LHCb experiment, corresponding to an integrated luminosity of $9\, \mathrm{fb}^{-1}$. The results are consistent with the Standard Model and previous tests of lepton universality in related decay modes. The first observation of $B^0 \to K^0_S e^+ e^-$ and $B^+ \to K^{*+} e^+ e^-$ decays is reported., Comment: All figures and tables, along with any supplementary material and additional information, are available at https://cern.ch/lhcbproject/Publications/p/LHCb-PAPER-2021-038.html. v2 - Added two more leptoquark references, and fixed formatting bugs in arXiv title and abstract. v3 - Version of paper published in PRL

We consider a framework for non-ideal magnetohydrodynamics in general relativity, paying particular attention to the physics involved. The discussion highlights the connection between the microphysics (associated with a given equation of state) and the global dynamics (from the point of view of numerical simulations), and includes a careful consideration of the assumptions that lead to ideal and resistive magnetohydrodynamics. We pay particular attention to the issue of local charge neutrality, which tends to be assumed but appears to be more involved than is generally appreciated. While we do not resolve all the involved issues, we highlight how some of the assumptions and simplifications may be tested by simulations. The final formulation is consistent, both logically and physically, preparing the ground for a new generation of models of relevant astrophysical scenarios., Comment: 15 pages, 3 figures, MNRAS format

Rapidly spinning, deformed neutron stars have long been considered potential gravitational-wave emitters. However, so far only upper limits on the size of the involved quadrupole deformations have been obtained. For this reason, it is pertinent to ask how large a mountain can be before the neutron star crust fractures. This is the question we consider in this paper, which describes how mountains can be calculated in relativistic gravity. Formally, this is a perturbative calculation that requires a fiducial force to source the mountain. Therefore, we consider three simple examples and increase their deforming amplitudes until the crust yields. We demonstrate how the derived mountains depend on the equation of state by considering a range of models obtained from chiral effective field theory. We find that the largest mountains depend sensitively on both the mechanism that sources them and the nuclear-matter equation of state., Comment: 13 pages, 3 figures. Accepted for publication in MNRAS with minor revision

We provide a bird's-eye view of neutron-star seismology, which aims to probe the extreme physics associated with these objects, in the context of gravitational-wave astronomy. Focussing on the fundamental mode of oscillation, which is an efficient gravitational-wave emitter, we consider the seismology aspects of a number of astrophysically relevant scenarios, ranging from transients (like pulsar glitches and magnetar flares), to the dynamics of tides in inspiralling compact binaries and the eventual merged object and instabilities acting in isolated, rapidly rotating, neutron stars. The aim is not to provide a thorough review, but rather to introduce (some of) the key ideas and highlight issues that need further attention., Comment: 20 pages, 2 figures, submitted for special issue in Universe

This brief review introduces the notion of a relativistic multifluid system -- a multi-component system with identifiable relative flows -- and outlines a set of models for scenarios relevant for different astronomical observation channels. The specific problems used to illustrate the key principles include superfluid hydrodynamics (with relevance for radio and x-ray pulsar timing and gravitational-wave searches), heat flow (connecting to the problem of neutron star cooling and associated x-ray observations) and the coupling between matter and electromagnetism (linking to explosive phenomena like gamma-ray bursts and more subtle issues like the long-term evolution of a neutron star's magnetic field). We also comment on the coupling between matter and radiation, for which the multifluid approach would seem less appropriate. The main motivation of the survey is to illustrate less familiar aspects that come into play in multifluid problems, establish the relevant "language" and provide a platform for more detailed work on these issues., Comment: 19 pages, no figures, to appear in special issue for Frontiers Astronomy and Space Sciences

As mature neutron stars are cold (on the relevant temperature scale), one has to carefully consider the state of matter in their interior. The outer kilometer or so is expected to freeze to form an elastic crust of increasingly neutron-rich nuclei, coexisting with a superfluid neutron component, while the star's fluid core contains a mixed superfluid/superconductor. The dynamics of the star depend heavily on the parameters associated with the different phases. The presence of superfluidity brings new degrees of freedom -- in essence we are dealing with a complex multi-fluid system -- and additional features: Bulk rotation is supported by a dense array of quantised vortices, which introduce dissipation via mutual friction, and the motion of the superfluid is affected by the so-called entrainment effect. This brief survey provides an introduction to -- along with a commentary on our current understanding of -- these dynamical aspects, paying particular attention to the role of entrainment, and outlines the impact of superfluidity on neutron-star seismology., Comment: 18 pages, 4 figures, contribution to special issue

As the era of gravitational-wave astronomy has well and truly begun, gravitational radiation from rotating neutron stars remains elusive. Rapidly spinning neutron stars are the main targets for continuous-wave searches since, according to general relativity, provided they are asymmetrically deformed, they will emit gravitational waves. It is believed that detecting such radiation will unlock the answer to why no pulsars have been observed to spin close to the break-up frequency. We review existing studies on the maximum mountain that a neutron star crust can support, critique the key assumptions and identify issues relating to boundary conditions that need to be resolved. In light of this discussion, we present a new scheme for modelling neutron star mountains. The crucial ingredient for this scheme is a description of the fiducial force which takes the star away from sphericity. We consider three examples: a source potential which is a solution to Laplace's equation, another solution which does not act in the core of the star and a thermal pressure perturbation. For all the cases, we find that the largest quadrupoles are between a factor of a few to two orders of magnitude below previous estimates of the maximum mountain size., Comment: 13 pages, 8 figures. Accepted for publication in MNRAS

The relativistic fluid is a highly successful model used to describe the dynamics of many-particle systems moving at high velocities and/or in strong gravity. It takes as input physics from microscopic scales and yields as output predictions of bulk, macroscopic motion. By inverting the process-e.g., drawing on astrophysical observations-an understanding of relativistic features can lead to insight into physics on the microscopic scale. Relativistic fluids have been used to model systems as "small" as colliding heavy ions in laboratory experiments, and as large as the Universe itself, with "intermediate" sized objects like neutron stars being considered along the way. The purpose of this review is to discuss the mathematical and theoretical physics underpinnings of the relativistic (multi-) fluid model. We focus on the variational principle approach championed by Brandon Carter and collaborators, in which a crucial element is to distinguish the momenta that are conjugate to the particle number density currents. This approach differs from the "standard" text-book derivation of the equations of motion from the divergence of the stress-energy tensor in that one explicitly obtains the relativistic Euler equation as an "integrability" condition on the relativistic vorticity. We discuss the conservation laws and the equations of motion in detail, and provide a number of (in our opinion) interesting and relevant applications of the general theory. The formalism provides a foundation for complex models, e.g., including electromagnetism, superfluidity and elasticity-all of which are relevant for state of the art neutron-star modelling., Comment: Update of invited review for Living Reviews in Relativity, original version from 2007. Substantial revision just under 250 pages with about 450 references. Comments welcome

It is well known that a superfluid rotates by forming an array of quantized vortices. A relativistic formulation for superfluid vortex dynamics is required for a range of problems in astrophysics and cosmology, from neutron star interiors and radio pulsar glitches to possible dark matter condensates on galactic scales. This paper develops a formalism for such systems, extending the well-established variational approach to relativistic fluids to account for the presence of a collection of quantized vortices. The model is firmly anchored in the geometry of the problem (drawing on aspects from basic string dynamics) and accounts for elastic aspects associated with a vortex array, providing a precise foundation for applications which have so far been based on somewhat ad hoc phenomenology., Comment: 29 pages, RevTeX preprint format, 1 figure

Finite size effects come into play during the late stages of neutron star binary inspiral, with the tidal deformability of the supranuclear density matter leaving an imprint on the gravitational-wave signal. As demonstrated in the case of GW170817, this leads to a constraint on the neutron star radius (and hence the equation of state). A deeper understanding of the tidal response requires an analysis of both the state and composition of matter. While these aspects may not have dramatic impact, they could lead to systematic effects that need to be kept in mind as the observational data become more precise. As a step in this direction we explore the role of the composition of matter, which is likely to remain "frozen" during the late stages of binary inspiral. We provide the first in-depth analysis of the problem, including estimates of how composition impacts on the effective tidal deformability. The results provide improved insight into how aspects of physics that tend to be "ignored" impact on binary neutron star gravitational-wave signals., Comment: 27 pages (RevTex preprint format), a few figures

We explore the correlation between the pre-merger tidal deformability and the post-merger remnant oscillations seen in numerical simulation of neutron star binaries, with the aim of understanding to what extent the physics support the existence of such a relation. We consider the impact of thermal and rotational effects and argue that the proposed relation does, indeed, make sense and provide simple arguments that help explain the result., Comment: 6 pages, 5 figures, MNRAS style

We discuss the impact of nuclear reactions on the spectrum of gravity g-modes of a mature neutron star, demonstrating the anticipated disappearance of these modes when the timescale associated with the oscillations is longer than that of nuclear reactions. This is the expected result, but different aspects of the demonstration may be relevant for related problems in neutron star astrophysics. In particular, we develop the framework required for an explicit implementation of finite-time nuclear reactions in neutron star oscillation problems and demonstrate how this formulation connects with the usual bulk viscosity prescription. We also discuss implications of the absence of very high order g-modes for problems of astrophysical relevance., Comment: 7 pages, 2 figures

We introduce a phenomenological, physically motivated, model for the effective tidal deformability of a neutron star, adding the frequency dependence (associated with the star's fundamental mode of oscillation) that comes into play during the late stages of the binary inspiral. Testing the model against alternative descriptions, we demonstrate that it provides an accurate representation of the dynamical tide up to close to merger. The simplicity of the prescription makes it an attractive alternative for a gravitational-wave data analysis implementation, facilitating an inexpensive construction of a large number of templates covering the relevant parameter space., Comment: Revised version with added details. To appear in MNRAS

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 fixed

In the standard form of the relativistic heat equation used in astrophysics, information propagates instantaneously, rather than being limited by the speed of light as demanded by relativity. We show how this equation nonetheless follows from a more general, causal theory of heat propagation in which the entropy plays the role of a fluid. In deriving this result, however, we see that it is necessary to make some assumptions which are not universally valid: the dynamical timescales of the process must be long compared with the explicitly causal physics of the theory, the heat flow must be sufficiently steady, and the spacetime static. Generalising the heat equation (e.g. restoring causality) would thus entail retaining some of the terms we neglected. As a first extension, we derive the heat equation for the spacetime associated with a slowly-rotating star or black hole, showing that it only differs from the static result by an additional advection term due to the rotation, and as a consequence demonstrate that a hotspot on a neutron star will be seen to be modulated at the rotation frequency by a distant observer., Comment: 13 pages, 0 figures

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 format

We discuss the role of dynamical tidal effects for inspiralling neutron star binaries, focussing on features that may be considered "unmodelled" in gravitational-wave searches. In order to cover the range of possibilities, we consider i) individual oscillation modes becoming resonant with the tide, ii) the elliptical instability, where a pair of inertial modes exhibit a nonlinear resonance with the tide, and iii) the non-resonant p-g instability which may arise as high order p- and g-modes in the star couple nonlinearly to the tide. In each case, we estimate the amount of additional energy loss that needs to be associated with the dynamical tide in order for the effect to impact on an observed gravitational-wave signal. We explore to what extent the involved neutron star physics may be considered known and how one may be able to use observational data to constrain theory., Comment: 17 pages, 4 figures

Neutron stars involve extreme physics which is difficult (perhaps impossible) to explore in laboratory experiments. We have to turn to astrophysical observations, and try to extract information from the entire range of the electromagnetic spectrum. In addition, neutron stars may radiate gravitational waves through a range of scenarios. In this brief summary I outline some of the main ideas, focussing on what we do and do not know, and describe the challenges involved in trying to catch these faint whispers from the very edge of physics., Comment: Has appeared in Journal of Astrophysics and Astronomy special issue on 'Physics of Neutron Stars and related objects', celebrating the 75th birth-year of G. Srinivasan

The young, fast-spinning, X-ray pulsar J0537-6910 displays an extreme glitch activity, with large spin-ups interrupting its decelerating rotation every ~100 days. We present nearly 13 years of timing data from this pulsar, obtained with the {\it Rossi X-ray Timing Explorer}. We discovered 22 new glitches and performed a consistent analysis of all 45 glitches detected in the complete data span. Our results corroborate the previously reported strong correlation between glitch spin-up size and the time to the next glitch, a relation that has not been observed so far in any other pulsar. The spin evolution is dominated by the glitches, which occur at a rate ~3.5 per year, and the post-glitch recoveries, which prevail the entire inter-glitch intervals. This distinctive behaviour provides invaluable insights into the physics of glitches. The observations can be explained with a multi-component model which accounts for the dynamics of the neutron superfluid present in the crust and core of neutron stars. We place limits on the moment of inertia of the component responsible for the spin-up and, ignoring differential rotation, the velocity difference it can sustain with the crust. Contrary to its rapid decrease between glitches, the spin-down rate increased over the 13 years, and we find the long-term braking index $n_{\rm l}=-1.22(4)$, the only negative braking index seen in a young pulsar. We briefly discuss the plausible interpretations of this result, which is in stark contrast to the predictions of standard models of pulsar spin-down., Comment: Minor changes to match the MNRAS accepted version

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 ApJ

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., Comment: 27 pages, 1 figure

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., Comment: 29 pages, no figures

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., Comment: 35 pages, no figures

We consider the detailed dynamics of an array of quantised superfluid vortices in the framework of general relativity, as required for quantitative modelling of realistic neutron star cores. Our model builds on the variational approach to relativistic (multi-) fluid dynamics, where the vorticity plays a central role. The description provides a natural extension of, and as it happens a better insight into, existing Newtonian models. In particular, we account for the mutual friction associated with scattering of a second "normal" component in the mixture off of the superfluid vortices., Comment: 9 pages, RevTeX

We consider the local dynamics of a realistic neutron star core, including composition gradients, superfluidity and thermal effects. The main focus is on the gravity g-modes, which are supported by composition stratification and thermal gradients. We derive the equations that govern this problem in full detail, paying particular attention to the input that needs to be provided through the equation of state and distinguishing between normal and superfluid regions. The analysis highlights a number of key issues that should be kept in mind whenever equation of state data is compiled from nuclear physics for use in neutron star calculations. We provide explicit results for a particular stellar model and a specific nucleonic equation of state, making use of cooling simulations to show how the local wave spectrum evolves as the star ages. Our results show that the composition gradient is effectively dominated by the muons whenever they are present. When the star cools below the superfluid transition, the support for g-modes at lower densities (where there are no muons) is entirely thermal. We confirm the recent suggestion that the g-modes in this region may be unstable, but our results indicate that this instability will be weak and would only be present for a brief period of the star's life. Our analysis accounts for the presence of thermal excitations encoded in entrainment between the entropy and the superfluid component. Finally, we discuss the complete spectrum, including the normal sound waves and, in superfluid regions, the second sound., Comment: 29 pages, 9 figures, submitted to MNRAS

This is a White Paper in support of the mission concept of the Large Observatory for X-ray Timing (LOFT), proposed as a medium-sized ESA mission. We discuss the potential of LOFT for the study of very faint X-ray binaries, orbital period distribution of black hole X-ray binaries and neutron star spin up. For a summary, we refer to the paper., Comment: White Paper in Support of the Mission Concept of the Large Observatory for X-ray Timing. (v2 few typos corrected)