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.

Neutron stars may harbour the true ground state of matter in the form of strange quark matter. If present, this type of matter is expected to be a color superconductor, a consequence of quark pairing with respect to the color/flavor degrees of freedom. The stellar magnetic field threading the quark core becomes a color-magnetic admixture and, in the event that superconductivity is of type II, leads to the formation of color-magnetic vortices. In this Letter we show that the volume-averaged color-magnetic vortex tension force should naturally lead to a significant degree of non-axisymmetry in systems like radio pulsars. We show that gravitational radiation from such color-magnetic `mountains' in young pulsars like the Crab and Vela could be observable by the future Einstein Telescope, thus becoming a probe of paired quark matter in neutron stars. The detectability threshold can be pushed up toward the sensitivity level of Advanced LIGO if we invoke an interior magnetic field about a factor ten stronger than the surface polar field.

We present the first fully relativistic calculations of the crustal strain induced in a neutron star by a binary companion at the late stages of inspiral, employing realistic equations of state for the fluid core and the solid crust. We show that while the deep crust is likely to fail only shortly before coalescence, there is a large variation in elastic strain, with the outermost layers failing relatively early on in the inspiral. We discuss the significance of the results for both electromagnetic and gravitational-wave astronomy.

We discuss the response of neutron stars to the tidal interaction in a compact binary system, as encoded in the Love number associated with the induced deformation. This problem is of interest for gravitational-wave astronomy as there may be a detectable imprint on the signal from the late stages of binary coalescence. Previous work has focussed on simple barotropic neutron star models, providing an understanding of the role of the stellar compactness and overall density profile. We add realism to the discussion by developing the framework required to model stars with varying composition and an elastic crust. These effects are not expected to be significant for the next generation of detectors but it is nevertheless useful to be able to quantify them. Our results show that (perhaps surprisingly) internal stratification has no impact whatsoever on the Love number. We also show that crust elasticity provides a (predictably) small correction to existing models.

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.

In this paper we reconsider the problem of magnetic field diffusion in neutron star cores. We model the star as consisting of a mixture of neutrons, protons and electrons, and allow for particle reactions and binary collisions between species. Our analysis is in much the same spirit as that of Goldreich & Reisenegger (1992), and we content ourselves with rough estimates of magnetic diffusion timescales, rather than solving accurately for some particular field geometry. However, our work improves upon previous treatments in one crucial respect: we allow for superfluidity in the neutron star matter. We find that the consequent mutual friction force, coupling the neutrons and charged particles, together with the suppression of particles collisions and reactions, drastically affect the ambipolar magnetic field diffusion timescale. In particular, the addition of superfluidity means that it is unlikely that there is ambipolar diffusion in magnetar cores on the timescale of the lifetimes of these objects, contradicting an assumption often made in the modelling of the flaring activity commonly observed in magnetars. Our work suggests that if a decaying magnetic field is indeed the cause of magnetar activity, the field evolution is likely to take place outside of the core, and might represent Hall/Ohmic diffusion in the stellar crust, or else that a mechanism other than standard ambipolar diffusion is active, e.g. flux expulsion due to the interaction between neutron vortices and magnetic fluxtubes.

Mature neutron stars are cold enough to contain a number of superfluid and superconducting components. These systems are distinguished by the presence of additional dynamical degrees of freedom associated with superfluidity. In order to consider models with mixtures of condensates we need to develop a multifluid description that accounts for the presence of rotational neutron vortices and magnetic proton fluxtubes. We also need to model the forces that impede the motion of vortices and fluxtubes, and understand how these forces act on the condensates. This paper concerns the development of such a model for the outer core of a neutron star, where superfluid neutrons co-exist with a type II proton superconductor and an electron gas. We discuss the hydrodynamics of this system, focusing on the role of the entrainment effect, the magnetic field, the vortex/fluxtube tension and the dissipative mutual friction forces. Out final results can be directly applied to a number of interesting astrophysical scenarios, e.g. associated with neutron star oscillations or the evolution of the large scale magnetic field.

A new formulation of Carter's constant for geodesic motion in Kerr black holes is given. It is shown that Carter's constant corresponds to the total angular momentum plus a precisely defined part which is quadratic in the linear momenta. The characterization is exact in the weak field limit obtained by letting the gravitational constant go to zero. It is suggested that the new form can be useful in current studies of the dynamics of extreme mass ratio inspiral (EMRI) systems emitting gravitational radiation.

We present the results of the first study of global oscillations of relativistic stars with both elastic crusts and interpenetrating superfluid components. For simplicity, we focus on the axial quasi-normal modes. Our results demonstrate that the torsional crust modes are essentially unaffected by the coupling to the gravitational field. This is as expected since these oscillations are known to be weak gravitational-wave sources. In contrast, the presence of a loosely coupled superfluid neutron component in the crust can have a significant effect on the oscillation spectrum. We show that the entrainment between the superfluid and the crust nuclei is a key parameter in the problem. Our analysis highlights the need for a more detailed understanding of the coupled crust-superfluid at the microphysical level. Our numerical results have, even though we have not considered magnetised stars, some relevance for efforts to carry out seismology based on quasi-periodic oscillations observed in the tails of magnetar flares. In particular, we argue that the sensitive dependence on the entrainment may have to be accounted for in attempts to match theoretical models to observational data.

We investigate the role of neutron star superfluidity for magnetar oscillations. Using a plane-wave analysis we estimate the effects of a neutron superfluid in the elastic crust region. We demonstrate that the superfluid imprint is likely to be more significant than the effects of the crustal magnetic field. We also consider the region immediately beneath the crust, where superfluid neutrons are thought to coexist with a type II proton superconductor. Since the magnetic field in the latter is carried by an array of fluxtubes, the dynamics of this region differs from standard magnetohydrodynamics. We show that the presence of the neutron superfluid (again) leaves a clear imprint on the oscillations of the system. Taken together, our estimates show that the superfluid components cannot be ignored in efforts to carry out "magnetar seismology". This increases the level of complexity of the modelling problem, but also points to the possibility of using observations to probe the superfluid nature of supranuclear matter.

We estimate the maximal deformation that can be sustained by a rotating neutron star with a crystalline colour superconducting quark core. Our results suggest that current gravitational-wave data from LIGO have already reached the level where a detection would have been possible over a wide range of the poorly constrained QCD parameters. This leads to the non-trivial conclusion that compact objects \emph{do not} contain maximally strained colour crystalline cores drawn from this range of parameter space. We discuss the uncertainties associated with our simple model and how it can be improved in the future.

Rotating deformed neutron stars are important potential sources for groundbased gravitational-wave interferometers such as LIGO, GE0600 and VIRGO. One mechanism that may lead to significant non-asymmetries is the internal magnetic field. It is well known that a magnetic star will not be spherical and, if the magnetic axis is not aligned with the spin axis, the deformation will lead to the emission of gravitational waves. The aim of this paper is to develop a formalism that would allow us to model magnetically deformed stars, using both realistic equations of state and field configurations. As a first step, we consider a set of simplified model problems. Focusing on dipolar fields, we determine the internal magnetic field which is consistent with a given neutron star model. We then calculate the associated deformation. We conclude by discussing the relevance of our results for current gravitational-wave detectors and future prospects.

This is the fourth paper in a series that attempt to put forward a consistent framework for modelling solid regions in neutron stars. Here we turn our attention to axial perturbations of spherically symmetric spacetimes using a gauge invariant approach due to one of us. Using the formalism developed in the first paper in the series it turns out that the matter perturbations are neatly expressible in terms of a ``metric'' tensor field depending only on the speeds of shear wave propagation along the principal directions in the solid. The results are applicable to a wide class of elastic materials and does not assume material isotropy nor quasi-Hookean behaviour. The perturbation equations are then specialised to a static background and are given by two coupled wave equations. Our formalism is thus slightly simpler than the previously existing results of Shumaker & Thorne, where an additional initial value equation needs to be solved. The simplification is mainly due to the gauge invariance of our approach and shows up also in somewhat simpler boundary conditions. We also give a first order formulation suitable for numerical integration of the quasi-normal mode problem of a neutron star. The relations between the gauge independent variables and the, in general, gauge dependent perturbed metric and strain tensor are explicitly given.

Recent observations of quasi-periodic oscillations in the aftermath of giant flares in soft gamma-ray repeaters suggest a close coupling between the seismic motion of the crust after a major quake and the modes of oscillations in a magnetar. In this paper we consider the purely elastic modes of oscillation in the crust of a neutron star in the relativistic Cowling approximation (disregarding any magnetic field). We determine the axial crust modes for a large set of stellar models, using a state-of-the-art crust equation of state and a wide range of core masses and radii. We also devise useful approximate formulae for the mode-frequencies. We show that the relative crust thickness is well described by a function of the compactness of the star and a parameter describing the compressibility of the crust only. Considering the observational data for SGR 1900+14 and SGR 1806-20, we demonstrate how our results can be used to constrain the mass and radius of an oscillating neutron star.

We use a simple toy-model to discuss global MHD modes of a neutron star, taking into account the magnetic coupling between the elastic crust and the fluid core. Our results suggest that the notion of pure torsional crust modes is not useful for the coupled system, all modes excite Alfven waves in the core. However, we also show that the modes that are most likely to be excited by a fractured crust, eg. during a magnetar flare, are such that the crust and the core oscillate in concert. For our simple model, the frequencies of these modes are similar to the ``pure crustal'' frequencies. In addition, our model provides a natural explanation for the presence of lower frequency ($< 30$~Hz) QPOs seen in the December 2004 giant flare of SGR 1806-20.

At densities below the neutron drip threshold, a purely elastic solid model (including, if necessary, a frozen-in magnetic field) can provide an adequate description of a neutron star crust, but at higher densities it will be necessary to allow for the penetration of the solid lattice by an independently moving current of superfluid neutrons. In order to do this, the previously available category of relativistic elasticity models is combined here with a separately developed category of relativistic superfluidity models in a unified treatment based on the use of an appropriate Lagrangian master function. As well as models of the purely variational kind, in which the vortices flow freely with the fluid, such a master function also provides a corresponding category of non-dissipative models in which the vortices are pinned to the solid structure.

We discuss ultracompact stellar objects which have multiple necks in their optical geometry. There are in fact physically reasonable equations of state for which the number of necks can be arbitrarily large. The proofs of these statements rely on a recent regularized formulation of the field equations for static spherically symmetric models due to Nilsson and Uggla. We discuss in particular the equation of state p=\rho-\rho_s which plays a central role in this context.