Using nonlinear, fully relativistic, simulations we investigate the dynamics and gravitational-wave signature associated with instabilities in neutron star post-merger remnants. For simplified models of the remnant, we establish the presence of instability in stars with moderate $T/|W|$, the ratio between the kinetic and the gravitational potential energies. Detailed analysis of the density oscillation pattern reveals a \emph{local} instability in the inner region of the more realistic differential rotation profile. We apply Rayleigh's inflection theorem and Fj{\o}rtoft's theorem to analyze the stability criteria concluding that this \emph{inner local} instability originates from a shear instability close to the peak of the angular velocity profile and that it later evolves into a fast-rotating $m=2$ oscillation pattern. We discuss the importance of the presence of a co-rotation point in the fluid, its connection with the shear instability, and comparisons to the Rossby Wave and Papaloizou Pringle Instabilities considered in the wider literature.

With the first detections of binary neutron star mergers by gravitational-wave detectors, it proves timely to consider how the internal structure of neutron stars affects the way in which they can be asymmetrically deformed. Such deformations may leave measurable imprints on gravitational-wave signals and can be sourced through tidal interactions or the formation of mountains. We detail the formalism that describes fully-relativistic neutron star models with elastic crusts undergoing static perturbations. This formalism primes the problem for studies into a variety of mechanisms that can deform a neutron star. We present results for a barotropic equation of state and a realistic model for the elastic crust, which enables us to compute relevant quantities such as the tidal deformability parameter. We find that the inclusion of an elastic crust provides a very small correction to the tidal deformability. The results allow us to demonstrate when and where the crust starts to fail during a binary inspiral and we find that the majority of the crust will remain intact up until merger.

Gravitational wave astronomy is expected to provide independent constraints on neutron star properties, such as their dense matter equation of state. This is possible with the measurements of binary components' tidal deformability, which alter the point-particle gravitational waveforms of the late inspiral phase of neutron-star binaries. Although current gravitational wave detectors are not sensitive enough for a precise determination of the individual tidal deformations of the components, a large number of combined observations with future detectors will decrease uncertainties in this quantity. Here we provide a first study of the tidal deformability effects due to the elasticity/solidity of the crust (hadronic phase) in a hybrid neutron star, as well as the influence of a quark-hadronic phase density jump on tidal deformations. We employ the framework of nonradial perturbations with zero frequency and study hadronic phases presenting elastic aspects when perturbed (with the shear modulus approximately $1\%$ of the pressure). We find that the relative tidal deformation change in a hybrid star with a perfect-fluid quark phase and a hadronic phase presenting an elastic part is never larger than about $2-4\%$ (with respect to a perfect-fluid counterpart). These maximum changes occur when the elastic region of a hybrid star is larger than approximately $60\%$ of the star's radius, which may happen when its quark phase is small and the density jump is large enough, or even when a hybrid star has an elastic mixed phase. For other cases, the relative tidal deformation changes due to an elastic crust are negligible ($10^{-5}-10^{-1}\%$), therefore unlikely to be measured even with third generation detectors. Thus, only when the size of the elastic hadronic region of a hybrid star is over half of its radius, the effects of elasticity could have a noticeable impact on tidal deformations.

During their most recent observing run, LIGO/Virgo reported the gravitational wave (GW) transient S191110af, a burst signal at a frequency of 1.78 kHz that lasted for 0.104 s. While this signal was later deemed non-astrophysical, genuine detections of uncertain origin will occur in the future. Here we study the potential for detecting GWs from neutron star fluid oscillations, which have mode frequency and duration matching those of S191110af and which can be used to constrain the equation of state of nuclear matter. Assuming that such transient oscillations can be excited to energies typical of a pulsar glitch, we use measured properties of known glitching pulsars to estimate the amplitude of GWs produced by such events. We find that current GW detectors may observe nearby pulsars undergoing large events with energy similar to Vela pulsar glitch energies, while next generation detectors could observe a significant number of events. Finally, we show that it is possible to distinguish between GWs produced by rapidly rotating and slowly rotating pulsars from the imprint of rotation on the f-mode frequency.

Implementing a family of differential rotation laws inspired by binary neutron-star merger remnants, we consider the impact of the rotation profile on the low-T/W instability. We use time evolutions of the linearised dynamical equations, in Newtonian gravity, to study non-axisymmetric oscillations and identify the unstable modes. The presence and evolution of the low-T/W instability is monitored with the canonical energy and angular momentum, while the growth time is extracted from the evolved kinetic energy. The results for the new rotation laws highlight similarities with the commonly considered j-constant law. The instability sets in when an oscillation mode co-rotates with the star (i.e. whenever there is a point where the mode's pattern speed matches the bulk angular velocity) and grows faster deep inside the co-rotation region. However, the new profiles add features, like an additional co-rotation point to the problem, which affect the onset of instability. The rotation laws influence more drastically the oscillation frequencies of the l=m=2 f-mode in fast rotating models, but affect the instability growth time at any rotation rate. We also identify models where the low-T/W instability appears to be triggered by inertial modes. We discuss to what extent the inferred qualitative behaviour is likely to be of observational relevance.

The application of standard accretion theory to observations of X-ray binaries provides valuable insights into neutron star properties, such as their spin period and magnetic field. However, most studies concentrate on relatively old systems, where the neutron star is in its late propeller, accretor, or nearly spin equilibrium phase. Here we use an analytic model from standard accretion theory to illustrate the evolution of high-mass X-ray binaries early in their life. We show that a young neutron star is unlikely to be an accretor because of the long duration of ejector and propeller phases. We apply the model to the recently discovered ~4000 yr old high-mass X-ray binary XMMU J051342.6-672412 and find that the system's neutron star, with a tentative spin period of 4.4 s, cannot be in the accretor phase and has a magnetic field B > (a few)x10^13 G, which is comparable to the magnetic field of many older high-mass X-ray binaries and is much higher than the spin equilibrium inferred value of (a few)x10^11 G. The observed X-ray luminosity could be the result of thermal emission from a young cooling magnetic neutron star or a small amount of accretion that can occur in the propeller phase.

Two coupled, interpenetrating fluids suffer instabilities beyond certain critical counterflows. For ideal fluids, an energetic instability occurs at the point where a sound mode inverts its direction due to the counterflow, while dynamical instabilities only occur at larger relative velocities. Here we discuss two relativistic fluids, one of which is dissipative. Using linearized hydrodynamics, we show that in this case the energetic instability turns dynamical, i.e., there is an exponentially growing mode, and this exponential growth only occurs in the presence of dissipation. This result is general and does not rely on an underlying microscopic theory. It can be applied to various two-fluid systems for instance in the interior of neutron stars. We also point out that under certain circumstances the two-fluid system exhibits a mode analogous to the r-mode in neutron stars that can become unstable for arbitrarily small values of the counterflow.

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.

We develop a phenomenological, physically motivated, model for the effective tidal deformability of a neutron star binary, adding the frequency dependence (associated with the star's fundamental mode of oscillation) that comes into play during the late stages of an inspiral. Testing the model against alternative descriptions, we find that it provides an accurate representation of the tidal deformability. The simplicity of the model makes it an attractive alternative for a gravitational-wave data analysis implementation as it facilitates an inexpensive construction of a large number of templates covering the relevant parameter space.

The fastest-spinning neutron stars in low-mass X-ray binaries, despite having undergone millions of years of accretion, have been observed to spin well below the Keplerian break-up frequency. We simulate the spin evolution of synthetic populations of accreting neutron stars in order to assess whether gravitational waves can explain this behaviour and provide the distribution of spins that is observed. We model both persistent and transient accretion and consider two gravitational-wave-production mechanisms that could be present in these systems: thermal mountains and unstable $r$-modes. We consider the case of no gravitational-wave emission and observe that this does not match well with observation. We find evidence for gravitational waves being able to provide the observed spin distribution; the most promising mechanisms being a permanent quadrupole, thermal mountains and unstable $r$-modes. However, based on the resultant distributions alone it is difficult to distinguish between the competing mechanisms.

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.

In this paper we describe the potential of the enhanced X-ray Timing and Polarimetry (eXTP) mission for studies related to accretion flows in the strong field gravity regime around both stellar-mass and supermassive black-holes. eXTP has the unique capability of using advanced 'spectral-timing-polarimetry' techniques to analyze the rapid variations with three orthogonal diagnostics of the flow and its geometry, yielding unprecedented insight into the inner accreting regions, the effects of strong field gravity on the material within them and the powerful outflows which are driven by the accretion process.

In this White Paper we present the potential of the enhanced X-ray Timing and Polarimetry (eXTP) mission for studies related to Observatory Science targets. These include flaring stars, supernova remnants, accreting white dwarfs, low and high mass X-ray binaries, radio quiet and radio loud active galactic nuclei, tidal disruption events, and gamma-ray bursts. eXTP will be excellently suited to study one common aspect of these objects: their often transient nature. Developed by an international Consortium led by the Institute of High Energy Physics of the Chinese Academy of Science, the eXTP mission is expected to be launched in the mid 2020s.

Pulsar glitches provide a unique way to study neutron star microphysics because short post-glitch dynamics are directly linked to strong frictional processes on small scales. To illustrate this connection between macroscopic observables and microphysics, we review calculations of vortex interactions focusing on Kelvin wave excitations and determine the corresponding mutual friction strength for realistic microscopic parameters in the inner crust. These density-dependent crustal coupling profiles are combined with a simplified treatment of the core coupling and implemented in a three-component neutron star model to construct a predictive framework for glitch rises. As a result of the density-dependent dynamics, we find the superfluid to transfer angular momentum to different parts of the crust and the core on different timescales. This can cause the spin frequency change to become non-monotonic in time, allowing for a maximum value much larger than the measured glitch size, as well as a delay in the recovery. The exact shape of the calculated glitch rise is strongly dependent on the relative strength between the crust and core mutual friction, providing the means to probe not only the crustal superfluid but also the deeper neutron star interior. To demonstrate the potential of this approach, we compare our predictive model with the first pulse-to-pulse observations recorded during the December 2016 glitch of the Vela pulsar. Our analysis suggests that the glitch rise behavior is relatively insensitive to the crustal mutual friction strength as long as $\mathcal{B} \gtrsim 10^{-3}$, while being strongly dependent on the core coupling strength, which we find to be in the range $3 \times 10^{-5} \lesssim \mathcal{B}_{\rm core} \lesssim 10^{-4}$.

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.

A superconductor of paired protons is thought to form in the core of neutron stars soon after their birth. Minimum energy conditions suggest magnetic flux is expelled from the superconducting region due to the Meissner effect, such that the neutron star core is largely devoid of magnetic fields for some nuclear equation of state and proton pairing models. We show via neutron star cooling simulations that the superconducting region expands faster than flux is expected to be expelled because cooling timescales are much shorter than timescales of magnetic field diffusion. Thus magnetic fields remain in the bulk of the neutron star core for at least 10^6-10^7 yr. We estimate the size of flux free regions at 10^7 yr to be <~ 100 m for a magnetic field of 10^11 G and possibly smaller for stronger field strengths. For proton pairing models that are narrow, magnetic flux may be completely expelled from a thin shell of approximately the above size after 10^5 yr. This shell may insulate lower conductivity outer layers, where magnetic fields can diffuse and decay faster, from fields maintained in the highly conducting deep core.

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.

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.

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.

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.

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.

Pulsars are rotating neutron stars that are renowned for their timing precision, although glitches can interrupt the regular timing behavior when these stars are young. Glitches are thought to be caused by interactions between normal and superfluid matter in the star. We update our recent work on a new technique using pulsar glitch data to constrain superfluid and nuclear equation of state models, demonstrating how current and future astronomy telescopes can probe fundamental physics such as superfluidity near nuclear saturation and matter at supranuclear densities. Unlike traditional methods of measuring a star's mass by its gravitational effect on another object, our technique relies on nuclear physics knowledge and therefore allows measurement of the mass of pulsars which are in isolation.

Neutron stars are astrophysical laboratories of many extremes of physics. Their rich phenomenology provides insights into the state and composition of matter at densities which cannot be reached in terrestrial experiments. Since the core of a mature neutron star is expected to be dominated by superfluid and superconducting components, observations also probe the dynamics of large-scale quantum condensates. The testing and understanding of the relevant theory tends to focus on the interface between the astrophysics phenomenology and nuclear physics. The connections with low-temperature experiments tend to be ignored. However, there has been dramatic progress in understanding laboratory condensates (from the different phases of superfluid helium to the entire range of superconductors and cold atom condensates). In this review, we provide an overview of these developments, compare and contrast the mathematical descriptions of laboratory condensates and neutron stars and summarise the current experimental state-of-the-art. This discussion suggests novel ways that we may make progress in understanding neutron star physics using low-temperature laboratory experiments.

Carter and Lichnerowicz have established that barotropic fluid flows are conformally geodesic and obey Hamilton's principle. This variational approach can accommodate neutral, or charged and poorly conducting, fluids. We show that, unlike what has been previously thought, this approach can also accommodate perfectly conducting magnetofluids, via the Bekenstein-Oron description of ideal magnetohydrodynamics. When Noether symmetries associated with Killing vectors or tensors are present in geodesic flows, they lead to constants of motion polynomial in the momenta. We generalize these concepts to hydrodynamic flows. Moreover, the Hamiltonian descriptions of ideal magnetohydrodynamics allow one to cast the evolution equations into a hyperbolic form useful for evolving rotating or binary compact objects with magnetic fields in numerical general relativity. Conserved circulation laws, such as those of Kelvin, Alfv\'en and Bekenstein-Oron, emerge simply as special cases of the Poincar\'e-Cartan integral invariant of Hamiltonian systems. We use this approach to obtain an extension of Kelvin's theorem to baroclinic (non-isentropic) fluids, based on a temperature-dependent time parameter. We further extend this result to perfectly or poorly conducting baroclinic magnetoflows. Finally, in the barotropic case, such magnetoflows are shown to also be geodesic, albeit in a Finsler (rather than Riemann) space.

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.

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.

The X-ray source 1E 161348-5055 in the supernova remnant RCW 103 recently exhibited X-ray activity typical of magnetars, i.e. neutron stars with magnetic fields > 10^14-10^15 G. However, 1E 161348-5055 has an observed period of 6.67 hr, in contrast to magnetars which have a spin period of seconds. Here we describe a simple model which can explain the spin evolution of 1E 161348-5055, as well as other magnetars, from an initial period of milliseconds that would be required for dynamo generation of magnetar-strength magnetic fields. We propose that the key difference between 1E 161348-5055 and other magnetars is the persistence of a remnant disk of small total mass. This disk caused 1E 161348-5055 to undergo ejector and propeller phases in its life, during which strong torques caused a rapid increase of its spin period. By matching its observed spin period and ~1-3 kyr age, we find that 1E 161348-5055 has the (slightly) highest magnetic field of all known magnetars, with B~5x10^15 G, and that its disk had a mass of ~10^24 g, comparable to that of the asteroid Ceres.

One of the primary science goals of the next generation of hard X-ray timing instruments is to determine the equation of state of the matter at supranuclear densities inside neutron stars, by measuring the radius of neutron stars with different masses to accuracies of a few percent. Three main techniques can be used to achieve this goal. The first involves waveform modelling. The flux we observe from a hotspot on the neutron star surface offset from the rotational pole will be modulated by the star's rotation, giving rise to a pulsation. Information about mass and radius is encoded into the pulse profile via relativistic effects, and tight constraints on mass and radius can be obtained. The second technique involves characterising the spin distribution of accreting neutron stars. The most rapidly rotating stars provide a very clean constraint, since the mass-shedding limit is a function of mass and radius. However the overall spin distribution also provides a guide to the torque mechanisms in operation and the moment of inertia, both of which can depend sensitively on dense matter physics. The third technique is to search for quasi-periodic oscillations in X-ray flux associated with global seismic vibrations of magnetars (the most highly magnetized neutron stars), triggered by magnetic explosions. The vibrational frequencies depend on stellar parameters including the dense matter equation of state. We illustrate how these complementary X-ray timing techniques can be used to constrain the dense matter equation of state, and discuss the results that might be expected from a 10m$^2$ instrument. We also discuss how the results from such a facility would compare to other astronomical investigations of neutron star properties. [Modified for arXiv]

Pulsars are known for their superb timing precision, although glitches can interrupt the regular timing behavior when the stars are young. These glitches are thought to be caused by interactions between normal and superfluid matter in the crust of the star. However, glitching pulsars such as Vela have been shown to require a superfluid reservoir that greatly exceeds that available in the crust. We examine a model in which glitches tap the superfluid in the core. We test a variety of theoretical superfluid models against the most recent glitch data and find that only one model can successfully explain up to 45 years of observational data. We develop a new technique for combining radio and X-ray data to measure pulsar masses, thereby demonstrating how current and future telescopes can probe fundamental physics such as superfluidity near nuclear saturation.

The presence of superconducting and superfluid components in the core of mature neutron stars calls for the rethinking of a number of key magnetohydrodynamical notions like resistivity, the induction equation, magnetic energy and flux-freezing. Using a multi-fluid magnetohydrodynamics formalism, we investigate how the magnetic field evolution is modified when neutron star matter is composed of superfluid neutrons, type-II superconducting protons and relativistic electrons. As an application of this framework, we derive an induction equation where the resistive coupling originates from the mutual friction between the electrons and the vortex/fluxtube arrays of the neutron and proton condensates. The resulting induction equation allows the identification of two timescales that are significantly different from those of standard magnetohydrodynamics. The astrophysical implications of these results are briefly discussed.

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.

Crustquake events may be connected with both rapid spin-up `glitches' within the regular slowdown of neutron stars, and high-energy magnetar flares. We argue that magnetic field decay builds up stresses in a neutron star's crust, as the elastic shear force resists the Lorentz force's desire to rearrange the global magnetic-field equilibrium. We derive a criterion for crust-breaking induced by a changing magnetic-field configuration, and use this to investigate strain patterns in a neutron star's crust for a variety of different magnetic-field models. Universally, we find that the crust is most liable to break if the magnetic field has a strong toroidal component, in which case the epicentre of the crustquake is around the equator. We calculate the energy released in a crustquake as a function of the fracture depth, finding that it is independent of field strength. Crust-breaking is, however, associated with a characteristic local field strength of $2.4\times 10^{14}$ G for a breaking strain of $0.001$, or $2.4\times 10^{15}$ G at a breaking strain of $0.1$. We find that even the most luminous magnetar giant flare could have been powered by crustal energy release alone.

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.

The SKA will discover tens of thousands of pulsars and provide unprecedented data quality on these, as well as the currently known population, due to its unrivalled sensitivity. Here, we outline the state of the art of our understanding of magnetospheric radio emission from pulsars and how we will use the SKA to solve the open problems in pulsar magnetospheric physics.

With an average density higher than the nuclear density, neutron stars (NS) provide a unique test-ground for nuclear physics, quantum chromodynamics (QCD), and nuclear superfluidity. Determination of the fundamental interactions that govern matter under such extreme conditions is one of the major unsolved problems of modern physics, and -- since it is impossible to replicate these conditions on Earth -- a major scientific motivation for SKA. The most stringent observational constraints come from measurements of NS bulk properties: each model for the microscopic behaviour of matter predicts a specific density-pressure relation (its `Equation of state', EOS). This generates a unique mass-radius relation which predicts a characteristic radius for a large range of masses and a maximum mass above which NS collapse to black holes. It also uniquely predicts other bulk quantities, like maximum spin frequency and moment of inertia. The SKA, in Phase 1 and particularly in Phase 2 will, thanks to the exquisite timing precision enabled by its raw sensitivity, and surveys that dramatically increase the number of sources: 1) Provide many more precise NS mass measurements (high mass NS measurements are particularly important for ruling out EOS models); 2) Allow the measurement of the NS moment of inertia in highly relativistic binaries such as the Double Pulsar; 3) Greatly increase the number of fast-spinning NS, with the potential discovery of spin frequencies above those allowed by some EOS models; 4) Improve our knowledge of new classes of binary pulsars such as black widows and redbacks (which may be massive as a class) through sensitive broad-band radio observations; and 5) Improve our understanding of dense matter superfluidity and the state of matter in the interior through the study of rotational glitches, provided that an ad-hoc campaign is developed.

We study the low T/W instability associated with the f-mode of differentially rotating stars. Our stellar models are described by a polytropic equation of state and the rotation profile is given by the standard j-constant law. The properties of the relevant oscillation modes, including the instability growth time, are determined from time evolutions of the linearised dynamical equations in Newtonian gravity. In order to analyse the instability we monitor also the canonical energy and angular momentum. Our results demonstrate that the l=m=2 f-mode becomes unstable as soon as a co-rotation point develops inside the star (i.e. whenever there is a point where the mode's pattern speed matches the bulk angular velocity). Considering various degrees of differential rotation, we show that the instability grows faster deep inside the co-rotation region and deduce an empirical relation that correlates the mode frequency and the star's parameters, which captures the main features of the l=m=2 f-mode growth time. This function is proportional to the product of the kinetic to gravitational energy ratio and the gradient of the star's spin, strengthening further the relationship between the co-rotation point and the low T/W instability. We briefly consider also the l=m=2 r-mode and demonstrate that it never moves far inside the co-rotation region even for significant differential rotation.

The Large Observatory For x-ray Timing (LOFT) was studied within ESA M3 Cosmic Vision framework and participated in the final down-selection for a launch slot in 2022-2024. Thanks to the unprecedented combination of effective area and spectral resolution of its main instrument, LOFT will study the behaviour of matter under extreme conditions, such as the strong gravitational field in the innermost regions of accretion flows close to black holes and neutron stars, and the supra-nuclear densities in the interior of neutron stars. The science payload is based on a Large Area Detector (LAD, 10 m 2 effective area, 2-30 keV, 240 eV spectral resolution, 1 deg collimated field of view) and a WideField Monitor (WFM, 2-50 keV, 4 steradian field of view, 1 arcmin source location accuracy, 300 eV spectral resolution). The WFM is equipped with an on-board system for bright events (e.g. GRB) localization. The trigger time and position of these events are broadcast to the ground within 30 s from discovery. In this paper we present the status of the mission at the end of its Phase A study.

We study the oscillations of relativistic stars, incorporating key physics associated with internal composition, thermal gradients and crust elasticity. Our aim is to develop a formalism which is able to account for the state-of-the-art understanding of the complex physics associated with these systems. As a first step, we build models using a modern equation of state including composition gradients and density discontinuities associated with internal phase-transitions (like the crust-core transition and the point where muons first appear in the core). In order to understand the nature of the oscillation spectrum, we carry out cooling simulations to provide realistic snapshots of the temperature distribution in the interior as the star evolves through adolescence. The associated thermal pressure is incorporated in the perturbation analysis, and we discuss the presence of $g$-modes arising as a result of thermal effects. We also consider interface modes due to phase-transitions and the gradual formation of the star's crust and the emergence of a set of shear modes.

Rapidly rotating neutron stars in Low Mass X-ray Binaries have been proposed as an interesting source of gravitational waves. In this chapter we present estimates of the gravitational wave emission for various scenarios, given the (electromagnetically) observed characteristics of these systems. First of all we focus on the r-mode instability and show that a 'minimal' neutron star model (which does not incorporate exotica in the core, dynamically important magnetic fields or superfluid degrees of freedom), is not consistent with observations. We then present estimates of both thermally induced and magnetically sustained mountains in the crust. In general magnetic mountains are likely to be detectable only if the buried magnetic field of the star is of the order of $B\approx 10^{12}$ G. In the thermal mountain case we find that gravitational wave emission from persistent systems may be detected by ground based interferometers. Finally we re-asses the idea that gravitational wave emission may be balancing the accretion torque in these systems, and show that in most cases the disc/magnetosphere interaction can account for the observed spin periods.

Recently Strohmayer and Mahmoodifar presented evidence for a coherent oscillation in the X-ray light curve of the accreting millisecond pulsar XTE J1751-305, using data taken by RXTE during the 2002 outburst of this source. They noted that a possible explanation includes the excitation of a non-radial oscillation mode of the neutron star, either in the form of a g-mode or an r-mode. The r-mode interpretation has connections with proposed spin-evolution scenarios for systems such as XTE J1751-305. Here we examine in detail this interesting possible interpretation. Using the ratio of the observed oscillation frequency to the star's spin frequency, we derive an approximate neutron star mass-radius relation which yields reasonable values for the mass over the range of expected stellar radius (as constrained by observations of radius-expansion burst sources). However, we argue that the large mode amplitude suggested by the Strohmayer and Mahmoodifar analysis would inevitably lead to a large spin-down of the star, inconsistent with its observed spin evolution, regardless of whether the r-mode itself is in a stable or unstable regime. We therefore conclude that the r-mode interpretation of the observed oscillation is not consistent with our current understanding of neutron star dynamics and must be considered unlikely. Finally we note that, subject to the availability of a sufficiently accurate timing model, a direct gravitational-wave search may be able to confirm or reject an r-mode interpretation unambiguously, should such an event, with a similar inferred mode amplitude, recur during the Advanced detector era.

Many pulsars are formed with a binary companion from which they can accrete matter. Torque exerted by accreting matter can cause the pulsar spin to increase or decrease, and over long times, an equilibrium spin rate is achieved. Application of accretion theory to these systems provides a probe of the pulsar magnetic field. We compare the large number of recent torque measurements of accreting pulsars with a high-mass companion to the standard model for how accretion affects the pulsar spin period. We find that many long spin period (P > 100 s) pulsars must possess either extremely weak (B < 10^10 G) or extremely strong (B > 10^14 G) magnetic fields. We argue that the strong-field solution is more compelling, in which case these pulsars are near spin equilibrium. Our results provide evidence for a fundamental link between pulsars with the slowest spin periods and strong magnetic fields around high-mass companions and pulsars with the fastest spin periods and weak fields around low-mass companions. The strong magnetic fields also connect our pulsars to magnetars and strong-field isolated radio/X-ray pulsars. The strong field and old age of our sources suggests their magnetic field penetrates into the superconducting core of the neutron star.

We consider a new mechanism for damping the oscillations of a mature neutron star. The new dissipation channel arises if superfluid vortices are forced to cut through superconducting fluxtubes. This mechanism is interesting because the oscillation modes need to exceed a critical amplitude in order for it to operate. Once it acts the effect is very strong (and nonlinear) leading to efficient damping. The upshot of this is that modes are unlikely to ever evolve far beyond the critical amplitude. We consider the effect of this new dissipation channel on the r-modes, that may be driven unstable by the emission of gravitational waves. Our estimates show that the fluxtube cutting leads to a saturation threshold for the instability that can be smaller than that of other proposed mechanisms. This suggests that the idea may be of direct astrophysical relevance.

We construct hydromagnetic neutron star equilibria which allow for a non-zero electric current distribution in the exterior. The novelty of our models is that the neutron star's interior field is in equilibrium with its magnetosphere, thus bridging the gap between previous work in this area which either solves for the interior assuming a vacuum exterior or solves for the magnetosphere without modelling the star itself. We consider only non-rotating stars in this work, so our solutions are most immediately applicable to slowly-rotating systems such as magnetars. Nonetheless, we demonstrate that magnetospheres qualitatively resembling those expected for both magnetars and pulsars are possible within our framework. The "inside-out" approach taken in this paper should be more generally applicable to rotating neutron stars, where the interior and exterior regions are again not independent but evolve together.

Interferometric detectors will very soon give us an unprecedented view of the gravitational-wave sky, and in particular of the explosive and transient Universe. Now is the time to challenge our theoretical understanding of short-duration gravitational-wave signatures from cataclysmic events, their connection to more traditional electromagnetic and particle astrophysics, and the data analysis techniques that will make the observations a reality. This paper summarizes the state of the art, future science opportunities, and current challenges in understanding gravitational-wave transients.

Pulsars are rotating neutron stars that are seen to slow down, and the spin-down rate is thought to be due to magnetic dipole radiation. This leads to a prediction for the braking index n, which is a combination of spin period and its first and second time derivatives. However, all observed values of n are below the predicted value of 3. Here we provide a simple model that can explain the rotational evolution of young pulsars, including the n=2.51 of the 958-year-old pulsar in the Crab nebula. The model is based on a decrease in effective moment of inertia due to an increase in the fraction of the stellar core that becomes superfluid as the star cools via neutrino emission. The results suggest that future large radio monitoring campaigns of pulsars will yield measurements of the neutron star mass, nuclear equation of state, and superfluid properties.

We present recent work on using astronomical observations of neutron stars to reveal unique insights into nuclear matter that cannot be obtained from laboratories on Earth. First, we discuss our measurement of the rapid cooling of the youngest neutron star in the Galaxy; this provides the first direct evidence for superfluidity and superconductivity in the supra-nuclear core of neutron stars. We show that observations of thermonuclear X-ray bursts on neutron stars can be used to constrain properties of neutron superfluidity and neutrino emission. We describe the implications of rapid neutron star rotation rates on aspects of nuclear and superfluid physics. Finally, we show that entrainment coupling between the neutron superfluid and the nuclear lattice leads to a less mobile crust superfluid; this result puts into question the conventional picture of pulsar glitches as being solely due to the crust superfluid and suggests that the core superfluid also participates.

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.

The LOFT mission concept is one of four candidates selected by ESA for the M3 launch opportunity as Medium Size missions of the Cosmic Vision programme. The launch window is currently planned for between 2022 and 2024. LOFT is designed to exploit the diagnostics of rapid X-ray flux and spectral variability that directly probe the motion of matter down to distances very close to black holes and neutron stars, as well as the physical state of ultra-dense matter. These primary science goals will be addressed by a payload composed of a Large Area Detector (LAD) and a Wide Field Monitor (WFM). The LAD is a collimated (<1 degree field of view) experiment operating in the energy range 2-50 keV, with a 10 m^2 peak effective area and an energy resolution of 260 eV at 6 keV. The WFM will operate in the same energy range as the LAD, enabling simultaneous monitoring of a few-steradian wide field of view, with an angular resolution of <5 arcmin. The LAD and WFM experiments will allow us to investigate variability from submillisecond QPO's to year-long transient outbursts. In this paper we report the current status of the project.

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 examine to what extent the inferred surface temperature of magnetars in quiescence can constrain the presence of a superfluid in the neutron star core and the role of magnetic field decay in the core. By performing detailed simulations of neutron star cooling, we show that extremely strong heating from field decay in the core cannot produce the high observed surface temperatures nor delay the onset of neutron superfluidity in the core. We verify the results of Kaminker et al., namely that the high magnetar surface temperatures require heating in the neutron star crust, and crust heating is decoupled from cooling/heating in the core. Therefore, because crust heating masks core heating, it is not possible to conclude that magnetar cores are in a non-superfluid state purely from high surface temperatures. From our interior temperature evolutions and after accounting for proton superconductivity in the core, we find that neutron superfluidity in the core occurs less than a few hundred years after neutron star formation. This onset time is unaffected by heating due to core field decay at fields < 10^16 G. Thus all known neutron stars, including magnetars, without a core containing exotic particles, should have a core of superfluid neutrons and superconducting protons.

Superfluid hydrodynamics affects the spin-evolution of mature neutron stars, and may be key to explaining timing irregularities such as pulsar glitches. However, most models for this phenomenon exclude the global instability required to trigger the event. In this paper we discuss a mechanism that may fill this gap. We establish that small scale inertial r-modes become unstable in a superfluid neutron star that exhibits a rotational lag, expected to build up due to vortex pinning as the star spins down. Somewhat counterintuitively, this instability arises due to the (under normal circumstances dissipative) vortex-mediated mutual friction. We explore the nature of the superfluid instability for a simple incompressible model, allowing for entrainment coupling between the two fluid components. Our results recover a previously discussed dynamical instability in systems where the two components are strongly coupled. In addition, we demonstrate for the first time that the system is secularly unstable (with a growth time that scales with the mutual friction) throughout much of parameter space. Interestingly, large scale r-modes are also affected by this new aspect of the instability. We analyse the damping effect of shear viscosity, which should be particularly efficient at small scales, arguing that it will not be sufficient to completely suppress the instability in astrophysical systems.

Pulsar glitches are traditionally viewed as a manifestation of vortex dynamics associated with a neutron superfluid reservoir confined to the inner crust of the star. In this Letter we show that the non-dissipative entrainment coupling between the neutron superfluid and the nuclear lattice leads to a less mobile crust superfluid, effectively reducing the moment of inertia associated with the angular momentum reservoir. Combining the latest observational data for prolific glitching pulsars with theoretical results for the crust entrainment we find that the required superfluid reservoir exceeds that available in the crust. This challenges our understanding of the glitch phenomenon, and we discuss possible resolutions to the problem.

The advanced interferometer network will herald a new era in observational astronomy. There is a very strong science case to go beyond the advanced detector network and build detectors that operate in a frequency range from 1 Hz-10 kHz, with sensitivity a factor ten better in amplitude. Such detectors will be able to probe a range of topics in nuclear physics, astronomy, cosmology and fundamental physics, providing insights into many unsolved problems in these areas.

Einstein gravitational-wave Telescope (ET) is a design study funded by the European Commission to explore the technological challenges of and scientific benefits from building a third generation gravitational wave detector. The three-year study, which concluded earlier this year, has formulated the conceptual design of an observatory that can support the implementation of new technology for the next two to three decades. The goal of this talk is to introduce the audience to the overall aims and objectives of the project and to enumerate ET's potential to influence our understanding of fundamental physics, astrophysics and cosmology.

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 consider general relativistic magnetohydrodynamics from a charged multifluids point-of-view, taking a variational approach as our starting point. We develop the case of two charged components in detail, accounting for a phenomenological resistivity, providing specific examples for pair plasmas and proton-electron systems. We discuss both cold, low velocity, plasmas and hot systems where we account for a dynamical entropy component. The results for the cold case (which accord with recent work in the literature) provide a complete model for resistive relativistic magnetohydrodynamics, clarifying the assumptions that lead to various models that have been used in astrophysical applications. The analysis of the hot case is (as far as we are aware) novel, accounting for the relaxation times that are required to ensure causality and demonstrating the explicit coupling between fluxes of heat and charge.

We construct two-fluid equilibrium configurations for neutron stars with magnetic fields, using a self-consistent and nonlinear numerical approach. The two-fluid approach - likely to be valid for large regions of all but the youngest NSs - provides us with a straightforward way to introduce stratification and allows for more realistic models than the ubiquitous barotropic assumption. In all our models the neutrons are modelled as a superfluid, whilst for the protons we consider two cases: one where they are a normal fluid and another where they form a type-II superconductor. We consider a variety of field configurations in the normal-proton case and purely toroidal fields in the superconducting case. We find that stratification allows for a stronger toroidal component in mixed-field configurations, though the poloidal component remains the largest in all our models. We provide quantitative results for magnetic ellipticities of NSs, both in the normal- and superconducting-proton cases.

Traditionally, the subject of hydromagnetic equilibrium in neutron stars has been addressed in the context of standard magnetohydrodynamics, with matter obeying a barotropic equation of state. In this paper we take a step towards a more realistic treatment of the problem by considering neutron stars with interior superfluid components. In this multifluid model stratification associated with a varying matter composition (the relative proton to neutron density fraction) enters as a natural ingredient, leading to a non-barotropic system. After formulating the hydromagnetic equilibrium of superfluid/superconducting neutron stars as a perturbation problem, we focus on the particular case of a three-fluid system consisting of superfluid neutrons and normal protons and electrons. We determine the equilibrium structure of dipolar magnetic fields with a mixed poloidal-toroidal composition. We find that, with respect to barotropic models, stratification has the generic effect of leading to equilibria with a higher fraction of magnetic energy stored in the toroidal component. However, even in models with strong stratification the poloidal and toroidal components are comparable, with the former contributing the bulk of the magnetic energy.

We consider the pinning of superfluid (neutron) vortices to magnetic fluxtubes associated with a type II (proton) superconductor in neutron star cores. We demonstrate that core pinning affects the spin-down of the system significantly, and discuss implications for regular radio pulsars and magnetars. We find that magnetars are likely to be in the pinning regime, while most radio pulsars are not. This suggests that the currently inferred magnetic field for magnetars may be overestimated. We also obtain a new timescale for the magnetic field evolution which could be associated with the observed activity in magnetars, provided that the field has a strong toroidal component.

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 consider the astrophysical constraints on the gravitational-wave driven r-mode instability in accreting neutron stars in low-mass X-ray binaries. We use recent results on superfluid and superconducting properties to infer the core temperature in these neutron stars and show the diversity of the observed population. Simple theoretical models indicate that many of these systems reside inside the r-mode instability region. However, this is in clear disagreement with expectations, especially for the systems containing the most rapidly rotating neutron stars. The inconsistency highlights the need to re-evaluate our understanding of the many areas of physics relevant to the r-mode instability. We summarize the current status of our understanding, and we discuss directions for future research which could resolve this dilemma.

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.

This paper revisits the problem of heat conduction in relativistic fluids, associated with issues concerning both stability and causality. It has long been known that the problem requires information involving second order deviations from thermal equilibrium. Basically, any consistent first-order theory needs to remain cognizant of its higher-order origins. We demonstrate this by carrying out the required first-order reduction of a recent variational model. We provide an analysis of the dynamics of the system, obtaining the conditions that must be satisfied in order to avoid instabilities and acausal signal propagation. The results demonstrate, beyond any reasonable doubt, that the model has all the features one would expect of a real physical system. In particular, we highlight the presence of a second sound for heat in the appropriate limit. We also make contact with previous work on the problem by showing how the various constraints on our system agree with previously established results.

We study the effects of an elastic crust on the oscillation spectrum of superfluid neutron stars. Within the two fluid formalism, we consider Newtonian stellar models that include the relevant constituents of a mature neutron stars. The core is formed by a mixture of superfluid neutrons and a conglomerate of charged particles, while the inner crust is described by a lattice of nuclei permeated by superfluid neutrons. We linearise the Poisson and the conservation equations of nonrotating superfluid stars and study the effects of elasticity, entrainment and composition stratification on the shear and acoustic modes. In both the core and the crust, the entrainment is derived from recent results for the nucleon effective mass. Solving the perturbation equations as an eigenvalue problem, we find that the presence of superfluid neutrons in the crust and their large effective mass may have significant impact on the star's oscillation spectrum.

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.

We compare the rotation rate of neutron stars in low-mass X-ray binaries (LMXBs) with the orbital period of the binaries. We find that, while short orbital period LMXBs span a range of neutron star rotation rates, all the long period LMXBs have fast rotators. We also find that the rotation rates are highest for the systems with the highest mean mass accretion rates, as can be expected if the accretion rate correlates with the orbital period. We show that these properties can be understood by a balance between spin-up due to accretion and spin-down due to gravitational radiation. Our scenario indicates that the gravitational radiation emitted by these systems may be detectable by future ground-based gravitational wave detectors.

Gravitational-wave astronomy is an area of great promise, yet to be realized. While we are waiting for the first (undisputed!) direct detection of these elusive waves it is useful to take stock and consider the challenges that need to be met if we want this field to reach its full potential. This write-up provides a brief introduction to some of the key ideas and the current state-of-play, and lists a range of modelling questions that need to be considered in the future.

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.

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.

In this paper we consider a simple two-fluid model for pulsar glitches. We derive the basic equations that govern the spin evolution of the system from two-fluid hydrodynamics, accounting for the vortex mediated mutual friction force that determines the glitch rise. This leads to a simple "bulk" model that can be used to describe the main properties of a glitch event resulting from vortex unpinning. In order to model the long term relaxation following the glitch our model would require additional assumptions regarding the repinning of vortices, an issue that we only touch upon briefly. Instead, we focus on comparing the phenomenological model to results obtained from time-evolutions of the linearised two-fluid equations, i.e. a "hydrodynamic" model for glitches. This allows us to study, for the first time, dynamics that was "averaged" in the bulk model, i.e. consider the various neutron star oscillation modes that are excited during a glitch. The hydro-results are of some relevance for efforts to detect gravitational waves from glitching pulsars, although the conclusions drawn from our rather simple model are pessimistic as far as the detectability of these events is concerned.

We study time evolutions of superfluid neutron stars, focussing on the nature of the oscillation spectrum, the effect of mutual friction force on the oscillations and the hydrodynamical spin-up phase of pulsar glitches. We linearise the dynamical equations of a Newtonian two-fluid model for rapidly rotating backgrounds. In the axisymmetric equilibrium configurations, the two fluid components corotate and are in beta-equilibrium. We use analytical equations of state that generate stratified and non-stratified stellar models, which enable us to study the coupling between the dynamical degrees of freedom of the system. By means of time evolutions of the linearised dynamical equations, we determine the spectrum of axisymmetric and non-axisymmetric oscillation modes, accounting for the contribution of the gravitational potential perturbations, i.e. without adopting the Cowling approximation. We study the mutual friction damping of the superfluid oscillations and consider the effects of the non-dissipative part of the mutual friction force on the mode frequencies. We also provide technical details and relevant tests for the hydrodynamical model of pulsar glitches discussed by Sidery, Passamonti and Andersson (2010). In particular, we describe the method used to generate the initial data that mimic the pre-glitch state, and derive the equations that are used to extract the gravitational-wave signal.

In order to establish whether the unstable r-modes in a rotating neutron star provide a detectable source of gravitational waves, we need to understand the details of the many dissipative processes that tend to counteract the instability. It has been established that the bulk viscosity due to exotic particles, like hyperons, may be particularly important in this respect. However, the effects of hyperon superfluidity have so far not been fully accounted for. While the associated suppression of the reaction rates that give rise to the bulk viscosity has been estimated, superfluid aspects of the fluid dynamics have not been considered. In this paper we determine the r-mode instability window for a neutron star with a $\Sigma^{-}$ hyperon core, using the appropriate multifluid formalism including, for the first time, the effect of the "superfluid" bulk viscosity coefficients. We demonstrate that, even though the extra terms may increase the bulk viscosity damping somewhat, their presence does not affect the qualitative features of the r-mode instability window.

We discuss different ways that neutron stars can generate gravitational waves, describe recent improvements in modelling the relevant scenarios in the context of improving detector sensitivity, and show how observations are beginning to test our understanding of fundamental physics. The main purpose of the discussion is to establish promising science goals for third-generation ground-based detectors, like the Einstein Telescope, and identify the various challenges that need to be met if we want to use gravitational-wave data to probe neutron star physics.

The effects of nonlinear oscillations in compact stars are attracting considerable current interest. In order to study such phenomena in the framework of fully nonlinear general relativity, highly accurate numerical studies are required. We present a numerical scheme specifically tailored for studies, based on formulating the time evolution in terms of deviations from a stationary equilibrium configuration. Using this technique, we investigate nonlinear effects associated with radial oscillations of neutron stars for a wide range of amplitudes. In particular, we discuss mode coupling due to nonlinear interactions, the occurrence of resonance phenomena, shock formation near the stellar surface as well as the capacity of nonlinearities to stabilize perturbatively unstable neutron star models.

Using time evolutions of the relevant linearised equations we study non-axisymmetric oscillations of rapidly rotating and superfluid neutron stars. We consider perturbations of Newtonian axisymmetric background configurations and account for the presence of superfluid components via the standard two-fluid model. Within the Cowling approximation, we are able to carry out evolutions for uniformly rotating stars up to the mass-shedding limit. This leads to the first detailed analysis of superfluid neutron star oscillations in the fast rotation regime, where the star is significantly deformed by the centrifugal force. For simplicity, we focus on background models where the two fluids (superfluid neutrons and protons) co-rotate, are in beta-equilibrium and coexist throughout the volume of the star. We construct sequences of rotating stars for two analytical model equations of state. These models represent relatively simple generalisations of single fluid, polytropic stars. We study the effects of entrainment, rotation and symmetry energy on non-radial oscillations of these models. Our results show that entrainment and symmetry energy can have a significant effect on the rotational splitting of non-axisymmetric modes. In particular, the symmetry energy modifies the inertial mode frequencies considerably in the regime of fast rotation.

We use time-evolutions of the linear perturbation equations to study the oscillations of rapidly rotating neutrons stars. Our models account for the buoyancy due to composition gradients and we study, for the first time, the nature of the resultant g-modes in a fast spinning star. We provide detailed comparisons of non-stratified and stratified models. This leads to an improved understanding of the relationship between the inertial modes of a non-stratified star and the g-modes of a stratified system. In particular, we demonstrate that each g-mode becomes rotation-dominated, i.e. approaches a particular inertial mode, as the rotation rate of the star is increased. We also discuss issues relating to the gravitational-wave driven instability of the various classes of oscillation modes.

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 develop a new perturbative framework for studying the r-modes of rotating superfluid neutron stars. Our analysis accounts for the centrifugal deformation of the star, and considers the two-fluid dynamics at linear order in the perturbed velocities. Our main focus is on a simple model system where the total density profile is that of an $n=1$ polytrope. We derive a partially analytic solution for the superfluid analogue of the classical r-mode. This solution is used to analyse the relevance of the vortex mediated mutual friction damping, confirming that this dissipation mechanism is unlikely to suppress the gravitational-wave driven instability in rapidly spinning superfluid neutron stars. Our calculation of the superfluid r-modes is significantly simpler than previous approaches, because it decouples the r-mode from all other inertial modes of the system. This leads to the results being clearer, but it also means that we cannot comment on the relevance of potential avoided crossings (and associated "resonances") that may occur for particular parameter values. Our analysis of the mutual friction damping differs from previous studies in two important ways. Firstly, we incorporate realistic pairing gaps which means that the regions of superfluidity in the star's core vary with temperature. Secondly, we allow the mutual friction parameters to take the whole range of permissible values rather than focussing on a particular mechanism. Thus, we consider not only the weak drag regime, but also the strong drag regime where the fluid dynamics is significantly different.

We investigate the oscillations of slowly rotating superfluid stars, taking into account the vortex mediated mutual friction force that is expected to be the main damping mechanism in mature neutron star cores. Working to linear order in the rotation of the star, we consider both the fundamental f-modes and the inertial r-modes. In the case of the (polar) f-modes, we work out an analytic approximation of the mode which allows us to write down a closed expression for the mutual friction damping timescale. The analytic result is in good agreement with previous numerical results obtained using an energy integral argument. We extend previous work by considering the full range of permissible values for the vortex drag, e.g. the friction between each individual vortex and the electron fluid. This leads to the first ever results for the f-mode in the strong drag regime. Our estimates provide useful insight into the dependence on, and relevance of, various equation of state parameters. In the case of the (axial) r-modes, we confirm the existence of two classes of modes. However, we demonstrate that only one of these sets remains purely axial in more realistic neutron star models. Our analysis lays the foundation for companion studies of the mutual friction damping of the r-modes at second order in the slow-rotation approximation, the first time evolutions for superfluid neutron star perturbations and also the first detailed attempt at studying the dynamics of superfluid neutron stars with both a relative rotation between the components and mutual friction.

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.

The standard explanation for large pulsar glitches involves transfer of angular momentum from an internal superfluid component to the star's crust. This model requires an instability to trigger the sudden unpinning of the vortices by means of which the superfluid rotates. This Letter describes a new instability that may play this role. The instability, which is associated with the inertial r-modes of a superfluid neutron star, sets in once the rotational lag in the system reaches a critical level. We demonstrate that our model is in good agreement with observed glitch data, suggesting that the superfluid r-mode instability may indeed be the mechanism that triggers large pulsar glitches.

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.

We discuss short wavelength (inertial wave) instabilities present in the standard two-fluid neutron star model when there is sufficient relative flow along the superfluid neutron vortex array. We demonstrate that these instabilities may be triggered in precessing neutron stars, since the angular velocity vectors of the neutron and proton fluids are misaligned during precession. The presence of such an instability would render the standard, solid body rotation, model for free precession inconsistent. Our results suggest that the standard (Eulerian) slow precession that results for weak drag between the vortices and the charged fluid (protons and electrons) is not seriously constrained by the existence of the instability. In contrast, the fast precession, which results when vortices are strongly coupled to the charged component, is generally unstable. This implies that fast precession may not be realised in astrophysical systems

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 a new superfluid instability occuring in the interior of mature neutron stars with implications for freely precessing neutron stars. This short-wavelength instability is similar to the instability which is responsible for the formation of turbulence in superfluid Helium. Its existence raises serious questions about our understanding of neutron star precession and complicates attempts to constrain neutron star interiors using such observations.

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.

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.

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.

Motivated by the possibility of radiation driven instabilities in rotating magnetic stars, we study the stability properties of general linear perturbations of a stationary and axisymmetric, infinitely conducting perfect fluid configuration threaded by a magnetic field and surrounded by vacuum. We develop a Lagrangian perturbation framework which enables us to formulate a strict stability criterion based on the notion of a canonical energy (a functional of the fluid displacement $\xi $ and its first time derivative). For any given choice of $\{\xi,\partial_t \xi \} $, the sign of the canonical energy determines whether the configuration is stable or not at the linear level. Our analysis provides the first complete description of the stability problem for a magnetic star, allowing for both rotation and the presence of a magnetic field in the exterior vacuum region. A key feature of the Lagrangian formulation is the existence of so-called `trivial' fluid displacements, which do not represent true physical perturbations. In order for the stability criterion to make rigorous sense one has to isolate these trivials and consider only the physical `canonical' displacements. We discuss this problem and formulate a condition which must be satisfied by all canonical displacements. Having obtained a well-defined stability criterion we provide examples which indicate that the magnetic field has a stabilising effect on radiation driven instabilities.

We examine radial oscillations of superfluid neutron stars at finite internal temperatures. For this purpose we generalize the description of relativistic superfluid hydrodynamics to the case of superfluid mixtures. We show that in a neutron star at hydrostatic and beta-equilibrium the red-shifted temperature gradient is smoothed out by neutron superfluidity (but not by proton superfluidity). We calculate radial oscillation modes of neutron stars assuming "frozen" nuclear composition in the pulsating matter. The resulting pulsation frequencies show a strong temperature dependence in the temperature range (0.1-1) T_cn, where T_cn is the critical temperature of neutron superfluidity. Combining our results with thermal evolution, we obtain a significant evolution of the pulsation spectrum, associated with highly efficient Cooper pairing neutrino emission, for 20 years after superfluidity onset.

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.

In this brief summary I describe our recent work on superfluid neutron star dynamics. I review results on shear viscosity, hyperon bulk viscosity, vortex mediated mutual friction and the modelling of multifluid systems in general. For each problem I provide a set of questions that need to be addressed by future work.

The aim of this paper is to compare the two cases of an isolated neutron star, with a non-accreted crust, and that of an accreting neutron star, with an accreted crust, and try to estimate which one of the two would make a better source of gravitational waves. In order to do this we must evaluate the maximum ``mountain'' that the crust can sustain in these two cases. We first do this using the formalism of Ushomirsky, Cutler and Bildsten (2000) and find that the maximum quadrupole is very similar in the two cases, with the non-accreted crust sustaining a slightly larger mountain. We then develop a perturbation formalism for the problem, that allows us to drop the Cowling approximation and have more control over the boundaries. The use of this formalism confirms that there is not much difference between the two cases, but leads to results approximately one order of magnitude larger than those we obtain with the formalism of Ushomirsky, Cutler and Bildsten (2000).

Motivated by their gravitational-wave driven instability, we investigate the influence of the crust on r-mode oscillations in a neutron star. Using a simplistic model of an elastic neutron star crust with constant shear modulus, we carry out an analytic calculation with the main objective of deriving an expression for the ``slippage'' between the core and the crust. Our analytic estimates support previous numerical results, and provide useful insights into the details of the problem.

We investigate the damping of neutron star r-modes due to the presence of a viscous boundary (Ekman) layer at the interface between the crust and the core. Our study is motivated by the possibility that the gravitational-wave driven instability of the inertial r-modes may become active in rapidly spinning neutron stars, eg. in low-mass X-ray binaries, and the fact that a viscous Ekman layer at the core-crust interface provides an efficient damping mechanism for these oscillations. We review various approaches to the problem and carry out an analytic calculation of the effects due to the Ekman layer for a rigid crust. Our analytic estimates support previous numerical results, and provide further insight into the intricacies of the problem. We add to previous work by discussing the effect that compressibility and composition stratification have on the boundary layer damping. We show that, while stratification is unimportant for the r-mode problem, composition suppresses the damping rate by about a factor of two (depending on the detailed equation of state).

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.

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.

In this paper we discuss the spin-equilibrium of accreting neutron stars in LMXBs. We demonstrate that, when combined with a naive spin-up torque, the observed data leads to inferred magnetic fields which are at variance with those of galactic millisecond radiopulsars. This indicates the need for either additional spin-down torques (eg. gravitational radiation) or an improved accretion model. We show that a simple consistent accretion model can be arrived at by accounting for radiation pressure in rapidly accreting systems (above a few percent of the Eddington accretion rate). In our model the inner disk region is thick and significantly sub-Keplerian, and the estimated equilibrium periods are such that the LMXB neutron stars have properties that accord well with the galactic millisecond radiopulsar sample. The implications for future gravitational-wave observations are also discussed briefly.

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.

This paper has been withdrawn by the authors.

Recent numerical simulations indicate the presence of dynamical instabilities of the f-mode in differentially rotating stars even at very low values of $T/|W|$, the ratio of kinetic to potential energy. In this Letter we argue that these may be shear instabilities which occur when the degree of differential rotation exceeds a critical value and the f-mode develops a corotation point associated with the presence of a continuous spectrum. Our explanation, which is supported by detailed studies of a simple shell model, offers a straightforward way of understanding all of the key features of these instabilities.

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.

An understanding of the oscillations of differentially rotating systems is key to many areas of astrophysics. It is of particular relevance to the emission of gravitational waves from oscillating neutron stars, which are expected to possess significant differential rotation immediately after birth or binary merger. In a previous paper we analysed the normal modes of a simple system exhibiting differential rotation. In this complementary paper we address the initial value problem for the same simple model using both analytical methods and numerical time evolutions. We derive a necessary and sufficient condition for dynamical shear instability. We discuss the dynamical behaviour of the continuous spectrum in response to an initial perturbation, and show that certain singular solutions within the continuous spectrum appear physically indistinguishable from the discrete modes outside the continuous spectrum.

We discuss the generalisation of the so-called traditional approximation, well known in geophysics, to general relativity. We show that the approximation is applicable for rotating relativistic stars provided that one focuses on relatively thin radial shells. This means that the framework can be used to study waves in neutron star oceans. We demonstrate that, once the effects of the relativistic frame-dragging are accounted for, the angular problem reduces to Laplace's tidal equation. We derive the dispersion relation for various classes of waves in a neutron star ocean and show that the combined effects of the frame-dragging and the gravitational redshift typically lower the frequency of a mode by about 20%.

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.

An understanding of the dynamics of differentially rotating systems is key to many areas of astrophysics. We investigate the oscillations of a simple system exhibiting differential rotation, and discuss issues concerning the role of corotation points and the emergence of dynamical instabilities. This problem is of particular relevance to the emission of gravitational waves from oscillating neutron stars, which are expected to possess significant differential rotation immediately after birth or binary merger.

In this article we show that the, mainly spacetime, w-modes may become unstable in a rotating ultracompact ($R<3M$) relativistic star. We provide results for the axial modes of a rotating star in the slow rotation approximation, and demonstrate that an instability may become active in stars with an ergosphere. These results suggest that the so-called ``ergoregion instability'' discussed by other authors is, in fact, associated with the w-modes of the star. The unstable modes are found to grow on a timescale of seconds to minutes. We discuss briefly whether the instability is likely to have any astrophysical relevance.

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.

This article provides an overview of stellar instabilities as sources of gravitational waves. The aim is to put recent work on secular and dynamical instabilities in compact stars in context, and to summarize the current thinking about the detectability of gravitational waves from various scenarios. As a new generation of kilometer length interferometric detectors are now coming online this is a highly topical theme. The review is motivated by two key questions for future gravitational-wave astronomy: Are the gravitational waves from various instabilities detectable? If so, what can these gravitational-wave signals teach us about neutron star physics? Even though we may not have clear answers to these questions, recent studies of the dynamical bar-mode instability and the secular r-mode instability have provided new insights into many of the difficult issues involved in modelling unstable stars as gravitational-wave sources.

We study the inertial modes of slowly rotating, fully relativistic compact stars. The equations that govern perturbations of both barotropic and non-barotropic models are discussed, but we present numerical results only for the barotropic case. For barotropic stars all inertial modes are a hybrid mixture of axial and polar perturbations. We use a spectral method to solve for such modes of various polytropic models. Our main attention is on modes that can be driven unstable by the emission of gravitational waves. Hence, we calculate the gravitational-wave growth timescale for these unstable modes and compare the results to previous estimates obtained in Newtonian gravity (i.e. using post-Newtonian radiation formulas). We find that the inertial modes are slightly stabilized by relativistic effects, but that previous conclusions concerning eg. the unstable r-modes remain essentially unaltered when the problem is studied in full general relativity.

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.

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.

The loss of angular momentum due to unstable r-modes in hot young neutron stars has been proposed as a mechanism for achieving the spin rates inferred for young pulsars. One factor that could have a significant effect on the action of the r-mode instability is fallback of supernova remnant material. The associated accretion torque could potentially counteract any gravitational-wave induced spin-down, and accretion heating could affect the viscous damping rates and hence the instability. We discuss the effects of various external agents on the r-mode instability scenario within a simple model of supernova fallback on to a hot young magnetized neutron star. We find that the outcome depends strongly on the strength of the star's magnetic field. Our model is capable of generating spin rates for young neutron stars that accord well with initial spin rates inferred from pulsar observations. The combined action of r-mode instability and fallback appears to cause the spin rates of neutron stars born with very different spin rates to converge, on a timescale of about a year. The results suggest that stars with magnetic fields $\le10^{13}$ G could emit a detectable gravitational wave signal for perhaps several years after the supernova event. Stars with higher fields (magnetars) are unlikely to emit a detectable gravitational wave signal via the r-mode instability. The model also suggests that the r-mode instability could be extremely effective in preventing young neutron stars from going dynamically unstable to the bar-mode.

We present the results of numerical time evolutions of the linearised perturbations of rapidly and rigidly rotating Newtonian polytropes while making the Cowling approximation. The evolution code runs stably for hundreds of stellar rotations, allowing us to compare our results with previously published eigenmode calculations, for instance the f-mode calculations of Ipser & Lindblom, and the r-mode calculations of Karino et al. The mode frequencies were found to be in agreement within the expected accuracy. We have also examined the inertial modes recently computed by Lockitch & Friedman, and we were able to extend their slow-rotation results into the rapid rotation regime. In the longer term, this code will provide a platform for studying a number of poorly understood problems in stellar oscillation theory, such as the effect of differential rotation and gravitational radiation reaction on normal mode oscillations and, with suitable modifications, mode-mode coupling in the mildly non-linear regime.

We investigate the relevance of the gravitational-wave driven r-mode instability for strange stars. We find that the unstable r-modes affect strange stars in a way that is quite distinct from the neutron star case. For accreting strange stars we show that the onset of r-mode instability does not lead to the thermo-gravitational runaway that is likely to occur in neutron stars. Instead, the strange star evolves towards a quasi-equilibrium state on a timescale of about a year. This mechanism could thus explain the clustering of spin-frequencies inferred from kHz QPO data in Low-mass X-ray binaries. For young strange stars we show that the r-mode driven spin-evolution is also distinct from the neutron star case. In a young strange star the r-mode undergoes short cycles of instability during the first few months. This is followed by a quasi-adiabatic phase where the r-mode remains at a small, roughly constant, amplitude for thousands of years. Another distinguishing feature from the neutron star case is that the r-modes in a strange star never grow to large amplitudes. Our results suggest that the r-modes in strange star emit a persistent gravitational-wave signal that should be detectable with large-scale interferometers given an observation time of a few months. If detected, these signals would provide unique evidence for the existence of strange stars, which would put useful constraints on the parameters of QCD.

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.

The non-linear behaviour of oscillation modes in compact stars is a topic of considerable current interest. Accurate numerical studies of such phenomena are likely to require powerful new approaches to both fluid and spacetime computations. We propose that a key ingredient of such methods will be the non-linear evolution of deviations from the background stationary equilibrium star. We investigate the feasibility of this approach by applying it to non-linear radial oscillations of a neutron star, and explore numerically various non-linear features of this problem, for a large range of amplitudes. Quadratic and higher order mode coupling and non-linear transfer of energy is demonstrated and analysed in detail.

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.

In this paper we model the gravitational wave emission of a freely precessing neutron star. The aim is to estimate likely source strengths, as a guide for gravitational wave astronomers searching for such signals. We model the star as a partly elastic, partly fluid body with quadrupolar deformations of its moment of inertia tensor. The angular amplitude of the free precession is limited by the finite breaking strain of the star's crust. The effect of internal dissipation on the star is important, with the precession angle being rapidly damped in the case of a star with an oblate deformation. We then go on to study detailed scenarios where free precession is created and/or maintained by some astrophysical mechanism. We consider the effects of accretion torques, electromagnetic torques, glitches and stellar encounters. We find that the mechanisms considered are either too weak to lead to a signal detectable by an Advanced LIGO interferometer, or occur too infrequently to give a reasonable event rate. We therefore conclude that, using our stellar model at least, free precession is not a good candidate for detection by the forthcoming laser interferometers.

We present a model of a freely precessing neutron star which is then compared against pulsar observations. The aim is to draw conclusions regarding the structure of the star, and test theoretical ideas of crust-core coupling and superfluidity. We argue that, on theoretical grounds, it is likely that the core neutron superfluid does not participate in the free precession of the crust. We apply our model to the handful of proposed observations of free precession that have appeared in the literature. Assuming crust-only precession, we find that all but one of the observations are consistent with there not being any pinned crustal superfluid at all; the maximum amount of pinned superfluid consistent with the observations is about 10^-10 of the total stellar moment of inertia. However, the observations do not rule out the possibility that the crust and neutron superfluid core precess as a single unit. In this case the maximum amount of pinned superfluid consistent with the observations is about 10^-8 of the total stellar moment of inertia. Both of these values are many orders of magnitude less than the 10^-2 value predicted by many theories of pulsar glitches. We conclude that superfluid pinning, at least as it affects free precession, needs to be reconsidered.

In this review we summarize the current understanding of the gravitational-wave driven instability associated with the so-called r-modes in rotating neutron stars. We discuss the nature of the r-modes, the detailed mechanics of the instability and its potential astrophysical significance. In particular we discuss results regarding the spin-evolution of nascent neutron stars, the detectability of r-mode gravitational waves and mechanisms limiting the spin-rate of accreting neutron stars in binary systems.

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 study the r-modes and rotational ``hybrid'' modes of relativistic stars. As in Newtonian gravity, the spectrum of low-frequency rotational modes is highly sensitive to the stellar equation of state. If the star and its perturbations obey the same one-parameter equation of state (as with isentropic stars), there exist {\it no pure r-modes at all} - no modes whose limit, for a star with zero angular velocity, is an axial-parity perturbation. Rotating stars of this kind similarly have no pure g-modes, no modes whose spherical limit is a perturbation with polar parity and vanishing perturbed pressure and density. We compute the post-Newtonian corrections to the $l=m$ r-modes of isentropic and non-isentropic uniform density stars.

We present a simple spin evolution model that predicts that rapidly rotating accreting neutron stars will mainly be confined to a narrow range of spin-frequencies; P= 1.5-5 ms. This is in agreement with current observations of both neutron stars in the Low-Mass X-ray Binaries and millisecond radiopulsars. The main ingredients in the model are; i) the instability of r-modes above a critical spin rate, ii) thermal runaway due to heat released as viscous damping mechanisms counteract the r-mode growth, and iii) a revised estimate of the strength of dissipation due to the presence of a viscous boundary layer at the base of the crust in an old and relatively cold neutron star. We discuss the gravitational waves that are radiated during the brief r-mode driven spin-down phase. We also briefly touch on how the new estimates affects the predicted initial spin periods of hot young neutron stars.

We develop a framework for constructing initial data sets for perturbations about spherically symmetric matter distributions. This framework facilitates setting initial data representing astrophysical sources of gravitational radiation involving relativistic stars. The procedure is based on Lichnerowicz-York's conformal approach to solve the constraints in Einstein's equations. The correspondence of these initial data sets in terms of the standard gauge perturbation variables in the Regge-Wheeler perturbation variables is established, and examples of initial data sets of merging neutron stars under the close-limit approximation are presented.

We develop a close-limit approximation to the head-on collision of two neutron stars similar to that used to treat the merger of black hole binaries . This approximation can serve as a useful benchmark test for future fully nonlinear studies. For neutron star binaries, the close-limit approximation involves assuming that the merged object can be approximated as a perturbed, stable neutron star during the ring-down phase of the coalescence. We introduce a prescription for the construction of initial data sets, discuss the physical plausibility of the various assumptions involved, and briefly investigate the character of the gravitational radiation produced during the merger. The numerical results show that several of the merged objects fluid pulsation modes are excited to a significant level.

We study the problem of detecting, and infering astrophysical information from, gravitational waves from a pulsating neutron star. We show that the fluid f and p-modes, as well as the gravitational-wave w-modes may be detectable from sources in our own galaxy, and investigate how accurately the frequencies and damping rates of these modes can be infered from a noisy gravitational-wave data stream. Based on the conclusions of this discussion we propose a strategy for revealing the supranuclear equation of state using the neutron star fingerprints: the observed frequencies of an f and a p-mode. We also discuss how well the source can be located in the sky using observations with several detectors.

We present a case study for the relevance of the r-mode instability for accreting compact stars. Our estimates are based on approximations that facilitate back-of-the-envelope calculations. We discuss two different cases: 1) For recycled millisecond pulsars we argue that the r-mode instability may be active at rotation periods longer than the Kepler period (which provides the dynamical limit on rotation) as long as the core temperature is larger than about $2\times10^5$ K. Our estimates suggest that the instability may have played a role in the evolution of the fastest spinning pulsars, and that it may be presently active in the recently discovered 2.49 ms X-ray pulsar SAX J1808.4-3658 as well as the rapidly spinning neutron stars observed in low mass X-ray binaries. This provides a new explanation for the remarkably similar rotation periods inferred from kHz quasi-periodic oscillations in the LMXBs. The possibility that the rotation of recycled pulsars may be gravitational-radiation limited is interesting because the gravitational waves from a neutron star rotating at the instability limit may well be detectable with the new generation of interferometric detectors. 2) We also consider white dwarfs, and find that the r-mode instability may possibly be active in short period white dwarfs. Our order of magnitude estimates (for a white dwarf of $M=M_\odot$ and $R=0.01 R_\odot$ composed of $C^{12}$) show that the instability could be operating for rotational periods shorter than $P\approx 27-33$ s. This number is in interesting agreement with the observed periods ($>28$ s) of the rapidly spinning DQ Herculis stars. However, we find that the instability grows too slowly to affect the rotation of these stars significantly .

A newly discovered instability in rotating neutron stars, driven by gravitational radiation reaction acting on the stars' $r$-modes, is shown here to set an upper limit on the spin rate of young neutron stars. We calculate the timescales for growth of linear perturbations due to gravitational radiation reaction, and for dissipation by shear and bulk viscosity, working to second order in a slow-rotation expansion within a Newtonian polytropic stellar model. The results are very temperature-sensitive: in hot neutron stars ($T>10^9$ K), the lowest-order $r$-modes are unstable, while in colder stars they are damped by viscosity. These calculations have a number of interesting astrophysical implications. First, the $r$-mode instability will spin down a newly born neutron star to a period close to the initial period inferred for the Crab pulsar, probably between 10 and 20 ms. Second, as an initially rapidly rotating star star spins down, an energy equivalent to roughly 1% of a solar mass is radiated as gravitational waves, which makes the process an interesting source for detectable gravitational waves. Third, the $r$-mode instability rules out the scenario whereby millisecond pulsars are formed by accretion-induced collapse of a white dwarf: the new star would be hot enough to spin down to much slower rates.

We present new results for pulsating neutron stars. We have calculated the eigenfrequencies of the modes that one would expect to be the most important gravitational-wave sources: the fundamental fluid f-mode, the first pressure p-mode and the first gravitational-wave w-mode, for twelve realistic equations of state. From this numerical data we have inferred a set of ``empirical relations'' between the mode-frequencies and the parameters of the star (the radius R and the mass M). Some of these relation prove to be surprisingly robust, and we show how they can be used to extract the details of the star (radius, mass, eos) from observed modes with errors no larger than a few percent.

The first numerical study of axial (toroidal) pulsation modes of a slowly rotating relativistic star is presented. The calculation includes terms of first order in $\epsilon \equiv \Omega \sqrt{R^3/M}<<1$ ($R$ is the radius, $M$ is the mass and $\Omega$ is the rotation frequency of the star), and accounts for effects due to the coriolis force. Effects due to the centrifugal flattening of the star enter at order $\epsilon^2$ and are not included in the analysis. It is shown that increased rotation tends to decrease the damping times for prograde modes, while retrograde become longer lived. Specifically, we show that rotation affects the axial gravitational-wave $w$-modes in this way. We also present the first relativistic calculation of the so-called $r$-modes (analogous to Rossby waves in the Earth's oceans). These have frequencies of the same order of magnitude as the rotation frequency of the star. The presented results indicate that the $r$-modes are unstable due to the emission of gravitational radiation for \underline{all} rotating perfect fluid stars. This is interesting since the previously considered gravitational-wave instability associated with (for example) the $f$-mode of the star sets in at a critical rotation rate. Because they are unstable also for the slowest rotating stars the $r$-modes may well be of considerable astrophysical importance.

An extensive survey of gravitational-wave modes for uniform density stars is presented. The study covers stars ranging in compactness from $R/M=100$ to the limit of stability in general relativity: $R/M = 9/4$. We establish that polar and axial gravitational-wave modes exist for all these stellar models. Moreover, there are two distinct families of both axial and polar modes. We discuss the physics of these modes and argue that one family is primarily associated with the interior of the star, while the second family is mainly associated with the stellar surface. We also show that the problem of axial perturbations has all the essential features of the polar problem as far as gravitational waves are concerned. This means that the axial problem is much more important than has previously been assumed. We also find some surprising features, such as avoided crossings between the polar gravitational-wave modes and the Kelvin f-mode as the star becomes very compact. This seems to suggest that the f-mode should be considered on equal footing with the polar w-modes for ultracompact stars. All modes may have the main character of trapped modes inside the curvature potential barrier for $R/M < 3$.