Much of our physical intuition about black hole quasinormal modes in general relativity comes from the eikonal/geometric optics approximation. According to the well-established eikonal model, the fundamental quasinormal mode represents wavepackets orbiting in the vicinity of the black hole's geodesic photon ring, slowly peeling off towards the event horizon and infinity. Besides its strength as a "visualisation" tool, the eikonal approximation also provides a simple quantitative method for calculating the mode frequency, in close agreement with rigorous numerical results. In this paper we move away from Einstein's theory and its garden-variety black holes and go on to consider spherically symmetric black holes in modified theories of gravity through the lens of the eikonal approximation. The quasinormal modes of such black holes are typically described by a set of coupled wave equations for the various field degrees of freedom. Considering a general, theory-agnostic, system of two equations for two perturbed fields, we derive eikonal formulae for the complex fundamental quasinormal mode frequency. In addition we show that the eikonal modes can be related to the extremum of an effective potential and its associated "photon ring". As an application of our results we consider a specific example of a modified theory of gravity with known black hole quasinormal modes and find that these are well approximated by the eikonal formulae.

The grand challenges of contemporary fundamental physics---dark matter, dark energy, vacuum energy, inflation and early universe cosmology, singularities and the hierarchy problem---all involve gravity as a key component. And of all gravitational phenomena, black holes stand out in their elegant simplicity, while harbouring some of the most remarkable predictions of General Relativity: event horizons, singularities and ergoregions. The hitherto invisible landscape of the gravitational Universe is being unveiled before our eyes: the historical direct detection of gravitational waves by the LIGO-Virgo collaboration marks the dawn of a new era of scientific exploration. Gravitational-wave astronomy will allow us to test models of black hole formation, growth and evolution, as well as models of gravitational-wave generation and propagation. It will provide evidence for event horizons and ergoregions, test the theory of General Relativity itself, and may reveal the existence of new fundamental fields. The synthesis of these results has the potential to radically reshape our understanding of the cosmos and of the laws of Nature. The purpose of this work is to present a concise, yet comprehensive overview of the state of the art in the relevant fields of research, summarize important open problems, and lay out a roadmap for future progress.

With the doors beginning to swing open on the new gravitational wave astronomy, this review provides an up-to-date survey of the most important physical mechanisms that could lead to emission of potentially detectable gravitational radiation from isolated and accreting neutron stars. In particular we discuss the gravitational wave-driven instability and asteroseismology formalism of the f- and r-modes, the different ways that a neutron star could form and sustain a non-axisymmetric quadrupolar "mountain" deformation, the excitation of oscillations during magnetar flares and the possible gravitational wave signature of pulsar glitches. We focus on progress made in the recent years in each topic, make a fresh assessment of the gravitational wave detectability of each mechanism and, finally, highlight key problems and desiderata for future work.

We make the first measurement of the braking index $n$ of two putative millisecond magnetars born in short gamma-ray bursts. We measure $n=2.9\pm0.1$ and $n=2.6\pm0.1$ for millisecond magnetars born in GRB 130603B and GRB 140903A respectively. The neutron star born in GRB 130603B has the only known braking index consistent with the fiducial $n=3$ value. This value is ruled out with 99.95\% confidence for GRB 140903A. We discuss possible causes of $n<3$ braking indices in millisecond magnetars, showing that several models can account for the measurement of the braking index in GRB 140903A, while it is more difficult to account for a braking index consistent with $n=3$.

One of the central goals of the newborn field of gravitational wave astronomy is to test gravity in the highly nonlinear, strong field regime characterizing the spacetime of black holes. In particular, "black hole spectroscopy" (the observation and identification of black hole quasinormal mode frequencies in the gravitational wave signal) is expected to become one of the main tools for probing the structure and dynamics of Kerr black holes. In this paper we take a significant step towards that goal by constructing a "post-Kerr" quasinormal mode formalism. The formalism incorporates a parametrized but general perturbative deviation from the Kerr metric and exploits the well-established connection between the properties of the spacetime's circular null geodesics and the fundamental quasinormal mode to provide approximate, eikonal limit formulae for the modes' complex frequencies. The resulting algebraic toolkit can be used in waveform templates for ringing black holes with the purpose of measuring deviations from the Kerr metric. As a first illustrative application of our framework, we consider the Johannsen-Psaltis deformed Kerr metric and compute the resulting deviation in the quasinormal mode frequency relative to the known Kerr result.

Our ability to interpret and glean useful information from the large body of observations of strongly magnetised neutron stars rests largely on our theoretical understanding of magnetic field equilibria. We answer the following question: is one free to arbitrarily prescribe magnetic equilibria such that fluid degrees of freedom can balance the equilibrium equations? We examine this question for various models for neutron star matter; from the simplest single-fluid barotrope to more realistic non-barotropic multifluid models with superfluid/superconducting components, muons and entropy. We do this for both axi- and non-axisymmetric equilibria, and in Newtonian gravity and general relativity. We show that, in axisymmetry, the most realistic model allows complete freedom in choosing a magnetic field equilibrium whereas non-axisymmetric equilibria are never completely arbitrary.

The bulk properties of spherically symmetric stars in general relativity can be obtained by integrating the Tolman-Oppenheimer-Volkoff (TOV) equations. In previous work we developed a "post-TOV" formalism - inspired by parametrized post-Newtonian theory - which allows us to classify in a parametrized, phenomenological form all possible perturbative deviations from the structure of compact stars in general relativity that may be induced by modified gravity at second post-Newtonian order. In this paper we extend the formalism to deal with the stellar exterior, and we compute several potential astrophysical observables within the post-TOV formalism: the surface redshift $z_s$, the apparent radius $R_{\rm app}$, the Eddington luminosity at infinity $L_{\rm E}^\infty$ and the orbital frequencies. We show that, at leading order, all of these quantities depend on just two post-TOV parameters $\mu_1$ and $\chi$, and we discuss the possibility to measure (or set upper bounds on) these parameters.

Observations of short gamma-ray bursts indicate ongoing energy injection following the prompt emission, with the most likely candidate being the birth of a rapidly rotating, highly magnetised neutron star. We utilise X-ray observations of the burst remnant to constrain properties of the nascent neutron star, including its magnetic field-induced ellipticity and the saturation amplitude of various oscillation modes. Moreover, we derive strict upper limits on the gravitational wave emission from these objects by looking only at the X-ray light curve, showing the burst remnants are unlikely to be detected in the near future using ground-based gravitational wave interferometers such as Advanced LIGO.

The study of long-term evolution of neutron star (NS) magnetic fields is key to understanding the rich diversity of NS observations, and to unifying their nature despite the different emission mechanisms and observed properties. Such studies in principle permit a deeper understanding of the most important parameters driving their apparent variety, e.g. radio pulsars, magnetars, x-ray dim isolated neutron stars, gamma-ray pulsars. We describe, for the first time, the results from self-consistent magneto-thermal simulations considering not only the effects of the Hall-driven field dissipation in the crust, but adding a complete set of proposed driving forces in a superconducting core. We emphasize how each of these core-field processes drive magnetic evolution and affect observables, and show that when all forces are considered together in vectorial form, the net expulsion of core magnetic flux is negligible, and will have no observable effect in the crust (consequently in the observed surface emission) on megayear time-scales. Our new simulations suggest that strong magnetic fields in NS cores (and the signatures on the NS surface) will persist long after the crustal magnetic field has evolved and decayed, due to the weak combined effects of dissipation and expulsion in the stellar core.

Besides their astrophysical interest, compact stars also provide an arena for understanding the properties of theories of gravity that differ from Einstein's general relativity. Numerous studies have shown that different modified theories of gravity can modify the bulk properties (such as mass and radius) of neutron stars for given assumptions on the microphysics. What is not usually stressed though is the strong degeneracy in the predictions of these theories for the stellar mass and radius. Motivated by this observation, in this paper we take an alternative route and construct a stellar structure formalism which, without adhering to any particular theory of gravity, describes in a simple parametrized form the departure from compact stars in general relativity. This "post-Tolman-Oppenheimer-Volkoff (TOV)" formalism for spherical static stars is inspired by the well-known parametrized post-Newtonian theory, extended to second post-Newtonian order by adding suitable correction terms to the fully relativistic TOV equations. We show how neutron star properties are modified within our formalism, paying special attention to the effect of each correction term. We also show that the formalism is equivalent to general relativity with an "effective" (gravity-modified) equation of state.

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.

It is commonly believed that the magnetic field threading a neutron star provides the ultimate mechanism (on top of fluid viscosity) for enforcing long-term corotation between the slowly spun down solid crust and the liquid core. We show that this argument fails for axisymmetric magnetic fields with closed field lines in the core, the commonly used `twisted torus' field being the most prominent example. The failure of such magnetic fields to enforce global crust-core corotation leads to the development of a persistent spin lag between the core region occupied by the closed field lines and the rest of the crust and core. We discuss the repercussions of this spin lag for the evolution of the magnetic field, suggesting that, in order for a neutron star to settle to a stable state of crust-core corotation, the bulk of the toroidal field component should be deposited into the crust soon after the neutron star's birth.

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.

Quasi-periodic oscillations have been seen in the light curves following several magnetar giant flares. These oscillations are of great interest as they probably provide our first ever view of the normal modes of oscillation of neutron stars. The state-of-the-art lies in the study of the oscillations of elastic-magnetic stellar models, mainly with a view to relating the observed frequencies to the structure and composition of the star itself. We advance this programme by considering several new physical mechanisms that are likely to be important for magnetar oscillations. These relate to the superfluid/superconducting nature of the stellar interior, and the damping of the modes, both through internal dissipation mechanisms and the launching of waves into the magnetosphere. We make simple order-of-magnitude estimates to show that both the frequencies and the damping time of magnetar oscillations can evolve in time, identifying three distinct `pathways' that can be followed, depending upon the initial magnitude of the mode excitation. These results are interesting as they show that the information buried in magnetar QPOs may be even richer than previously thought, and motivate more careful examination of magnetar light curves, to search for signatures of the different types of evolution that we have identified.

Since the work of Hartle in the 1970s, and the subsequent development of the the Membrane Paradigm approach to black hole physics it has been widely accepted that superradiant scattering of gravitational waves bears strong similarities with the phenomenon of ``tidal friction'' (well-known from Newtonian gravity) operating in binary systems of viscous material bodies. In this paper we revisit the superradiance-tidal friction analogy within the context of ultracompact relativistic bodies. We advocate that as long as these bodies have non-zero viscosity they should undergo tidal friction that can be construed as a kind of superradiant scattering from the point of view of the dynamics of an orbiting test-body. In addition we consider the presence of anisotropic matter, which is required for at least some ultracompact bodies, if they are to sustain a radius very close to the gravitational radius. We find that the tidal friction/superradiance output is enhanced with increasing anisotropy and that strongly anisotropic systems exhibit an unconventional response to tidal and centrifugal forces. Finally, we make contact with the artificial system comprising a black hole with its horizon replaced by a mirror (sometimes used as a proxy for ultracompact material bodies) and discuss superradiance and tidal friction in relation to it.

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.

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.

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.

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

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.

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.

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.

We study the damping of the gravitational radiation-driven f-mode instability in rotating neutron stars by nonlinear bulk viscosity in the so-called supra-thermal regime. In this regime the dissipative action of bulk viscosity is known to be enhanced as a result of nonlinear contributions with respect to the oscillation amplitude. Our analysis of the f-mode instability is based on a time-domain code that evolves linear perturbations of rapidly rotating polytropic neutron star models. The extracted mode frequency and eigenfunctions are subsequently used in standard energy integrals for the gravitational wave growth and viscous damping. We find that nonlinear bulk viscosity has a moderate impact on the size of the f-mode instability window, becoming an important factor and saturating the mode's growth at a relatively large oscillation amplitude. We show similarly that nonlinear bulk viscosity leads to a rather high saturation amplitude even for the r-mode instability. In addition, we show that the action of bulk viscosity can be significantly mitigated by the presence of superfluidity in neutron star matter. Apart from revising the f-mode's instability window we provide results on the mode's gravitational wave detectability. Considering an f-mode-unstable neutron star located in the Virgo cluster and assuming a mode amplitude at the level allowed by bulk viscosity, we find that the emitted gravitational wave signal could be detectable by advanced ground-based detectors such as Advanced LIGO/Virgo and the Einstein Telescope.

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.

Rapidly spinning neutron stars are known to harbour pulsation modes that may become unstable and grow in amplitude by emitting gravitational radiation. Among the various stellar modes, the f-mode is the one typically considered as a promising source of gravitational radiation for ground-based detectors such as LIGO and VIRGO. Improving the existing work in Newtonian stellar models, we present the first calculation of the basic properties of the f-mode instability in rapidly rotating relativistic neutron stars, adopting the Cowling approximation. Using a relativistic polytropic stellar model, we obtain a minimum gravitational growth timescale (for the dominant l=m=4 mode) of the order of 10^3-10^4 s near the Kepler spin frequency Omega_K, which is substantially shorter than the Newtonian value. By accounting for dissipation in neutron star matter, i.e. shear/bulk viscosity and superfluid mutual friction, we calculate the associated f-mode instability window. For our specific stellar model, the instability is active above 0.92 \times Omega_K and for temperatures \sim (10^9 - 2 \times 10^{10}) K, characteristic of newborn neutron stars.

We model the non-linear ideal magnetohydrodynamics of poloidal magnetic fields in neutron stars in general relativity assuming a polytropic equation of state. We identify familiar hydromagnetic modes, in particular the 'sausage/varicose' mode and 'kink' instability inherent to poloidal magnetic fields. The evolution is dominated by the kink instability, which causes a cataclysmic reconfiguration of the magnetic field. The system subsequently evolves to new, non-axisymmetric, quasi-equilibrium end-states. The existence of this branch of stable quasi-equilibria may have consequences for magnetar physics, including flare generation mechanisms and interpretations of quasi-periodic oscillations.

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

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

The objects known as anomalous X-ray pulsars and soft gamma repeaters are commonly identified with magnetars, neutron stars with ultrastrong magnetic fields. The rotational history of these objects has, so far, revealed no evidence of free precession. At the same time these objects do not generally appear to have magnetic axes nearly parallel or orthogonal to their spin axes. In this paper we show that the combination of these two observations, together with simple rigid-body dynamics, leads to non-trivial predictions about the interior properties of magnetars: either (i) elastic stresses in magnetar crusts are close to the theoretical upper limit above which the crustal matter yields or (ii) there is a "pinned" superfluid component in the magnetar interior. As a potentially observable consequence of these ideas we point out that, in the case of no pinned superfluidity, magnetars of stronger magnetic field strength than those currently observed would have to be nearly aligned/orthogonal rotators.

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 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

One of the most exciting potential sources of gravitational waves for low-frequency, space-based gravitational wave (GW) detectors such as the proposed Laser Interferometer Space Antenna (LISA) is the inspiral of compact objects into massive black holes in the centers of galaxies. The detection of waves from such "extreme mass ratio inspiral" systems (EMRIs) and extraction of information from those waves require template waveforms. The systems' extreme mass ratio means that their waveforms can be determined accurately using black hole perturbation theory. Such calculations are computationally very expensive. There is a pressing need for families of approximate waveforms that may be generated cheaply and quickly but which still capture the main features of true waveforms. In this paper, we introduce a family of such "kludge" waveforms and describe ways to generate them. We assess performance of the introduced approximations by comparing "kludge" waveforms to accurate waveforms obtained by solving the Teukolsky equation in the adiabatic limit (neglecting GW backreaction). We find that the kludge waveforms do extremely well at approximating the true gravitational waveform, having overlaps with the Teukolsky waveforms of 95% or higher over most of the parameter space for which comparisons can currently be made. Indeed, we find these kludges to be of such high quality (despite their ease of calculation) that it is possible they may play some role in the final search of LISA data for EMRIs.

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.

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.

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.

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.

The future LISA detector will constitute the prime instrument for high-precision gravitational wave observations.LISA is expected to provide information for the properties of spacetime in the vicinity of massive black holes which reside in galactic nuclei.Such black holes can capture stellar-mass compact objects, which afterwards slowly inspiral,radiating gravitational waves.The body's orbital motion and the associated waveform carry information about the spacetime metric of the massive black hole,and it is possible to extract this information and experimentally identify (or not!) a Kerr black hole.In this paper we lay the foundations for a practical `spacetime-mapping' framework. Our work is based on the assumption that the massive body is not necessarily a Kerr black hole, and that the vacuum exterior spacetime is stationary axisymmetric,described by a metric which deviates slightly from the Kerr metric. We first provide a simple recipe for building such a `quasi-Kerr' metric by adding to the Kerr metric the deviation in the value of the quadrupole moment. We then study geodesic motion in this metric,focusing on equatorial orbits. We proceed by computing `kludge' waveforms which we compare with their Kerr counterparts. We find that a modest deviation from the Kerr metric is sufficient for producing a significant mismatch between the waveforms, provided we fix the orbital parameters. This result suggests that an attempt to use Kerr waveform templates for studying EMRIs around a non-Kerr object might result in serious loss of signal-to-noise ratio and total number of detected events. The waveform comparisons also unveil a `confusion' problem, that is the possibility of matching a true non-Kerr waveform with a Kerr template of different orbital parameters.

In this review article I attempt to summarise past and present-ongoing-work on the problem of the inspiral of a small body in the gravitational field of a much more massive Kerr black hole. Such extreme mass ratio systems, expected to occur in galactic nuclei, will constitute prime sources of gravitational radiation for the future LISA gravitational radiation detector. The article's main goal is to provide a survey of basic celestial mechanics in Kerr spacetime and calculations of gravitational waveforms and backreaction on the small body's orbital motion, based on the traditional `flux-balance' method and the Teukolsky black hole perturbation formalism.

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.

We present a new approximate method for constructing gravitational radiation driven inspirals of test-bodies orbiting Kerr black holes. Such orbits can be fully described by a semi-latus rectum $p$, an eccentricity $e$, and an inclination angle $\iota$; or, by an energy $E$, an angular momentum component $L_z$, and a third constant $Q$. Our scheme uses expressions that are exact (within an adiabatic approximation) for the rates of change ($\dot{p}$, $\dot{e}$, $\dot{\iota}$) as linear combinations of the fluxes ($\dot{E}$, $\dot{L_z}$, $\dot{Q}$), but uses quadrupole-order formulae for these fluxes. This scheme thus encodes the exact orbital dynamics, augmenting it with approximate radiation reaction. Comparing inspiral trajectories, we find that this approximation agrees well with numerical results for the special cases of eccentric equatorial and circular inclined orbits, far more accurate than corresponding weak-field formulae for ($\dot{p}$, $\dot{e}$, $\dot{\iota}$). We use this technique to study the inspiral of a test-body in inclined, eccentric Kerr orbits. Our results should be useful tools for constructing approximate waveforms that can be used to study data analysis problems for the future LISA gravitational-wave observatory, in lieu of waveforms from more rigorous techniques that are currently under development.

We study eccentric equatorial orbits of a test-body around a Kerr black hole under the influence of gravitational radiation reaction. We have adopted a well established two-step approach: assuming that the particle is moving along a geodesic (justifiable as long as the orbital evolution is adiabatic) we calculate numerically the fluxes of energy and angular momentum radiated to infinity and to the black hole horizon, via the Teukolsky-Sasaki-Nakamura formalism. We can then infer the rate of change of orbital energy and angular momentum and thus the evolution of the orbit. The orbits are fully described by a semi-latus rectum $p$ and an eccentricity $e$. We find that while, during the inspiral, $e$ decreases until shortly before the orbit reaches the separatrix of stable bound orbits (which is defined by $p_{s}(e)$), in many astrophysically relevant cases the eccentricity will still be significant in the last stages of the inspiral. In addition, when a critical value $p_{crit}(e)$ is reached, the eccentricity begins to increase as a result of continued radiation induced inspiral. The two values $p_{s}$, $p_{crit}$ (for given $e$) move closer to each other, in coordinate terms, as the black hole spin is increased, as they do also for fixed spin and increasing eccentricity. Of particular interest are moderate and high eccentricity orbits around rapidly spinning black holes, with $p(e) \approx p_{s}(e)$. We call these ``zoom-whirl'' orbits, because of their characteristic behaviour involving several revolutions around the central body near periastron. Gravitational waveforms produced by such orbits are calculated and shown to have a very particular signature. Such signals may well prove of considerable astrophysical importance for the future LISA detector.

We study the dynamical evolution of perturbations in the gravitational field of a collapsing fluid star. Specifically, we consider the initial value problem for a massless scalar field in a spacetime similar to the Oppenheimer-Snyder collapse model, and numerically evolve in time the relevant wave equation. Our main objective is to examine whether the phenomenon of parametric amplification, known to be responsible for the strong amplification of primordial perturbations in the expanding Universe, can efficiently operate during gravitational collapse. Although the time-varying gravitational field inside the star can, in principle, support such a process, we nevertheless find that the perturbing field escapes from the star too early for amplification to become significant. To put an upper limit in the efficiency of the amplification mechanism (for a scalar field) we furthermore consider the case of perturbations trapped inside the star for the entire duration of the collapse. In this extreme case, the field energy is typically amplified at the level ~ 1% when the star is about to cross its Schwarszchild radius. Significant amplification is observed at later stages when the star has even smaller radius. Therefore, the conclusion emerging from our simple model is that parametric amplification is unlikely to be of significance during gravitational collapse. Further work, based on more realistic collapse models, is required in order to fully assess the astrophysical importance of parametric amplification.