As evidenced by the coincident detections of GW170817 and GRB 170817A, short gamma-ray bursts are likely associated with neutron star-neutron star merger events. Although rare, some bursts display episodes of early emission, with precursor flares being observed up to $\sim 10$ seconds prior to the main burst. As the stars inspiral due to gravitational wave emission, the exertion of mutual tidal forces leads to the excitation of stellar oscillation modes, which may come into resonance with the orbital motion. Mode amplitudes increase substantially during a period of resonance as tidal energy is deposited into the star. The neutron star crust experiences shear stress due to the oscillations and, if the resonant amplitudes are large enough, may become over-strained. This over-straining can lead to fractures or quakes which release energy, thereby fueling precursor activity prior to the merger. Using some simple Maclaurin spheroid models, we investigate the influence of magnetic fields and rapid rotation on tidally-forced $f$- and $r$- modes, and connect the associated eigenfrequencies with the orbital frequencies corresponding to precursor events seen in, for example, GRB 090510.

Many gamma-ray bursts are followed by periods of extended emission. At least in some cases, the burst afterglow may be powered by a rapidly rotating, highly-magnetised neutron star, which spins down due to electromagnetic and gravitational wave emission. Such a remnant is likely to strongly precess in the early stages of its life, which would lead to modulations in the X-ray luminosity as the triaxiality of the system evolves over time. Using a radiation profile appropriate for a precessing, oblique rotator, we find that Swift-XRT data of a long (080602) and a short (090510) burst matches the model with significantly higher accuracy (mean-square residuals dropping by $\gtrsim 200 \%$ in the early stages of the extended emission) than for an orthogonal rotator. We interpret this as evidence for precession in newborn magnetars.

Some pairs of neutron star models can intuitively be thought of as being `closer' together than others, in the sense that more precise observations might be required to distinguish between them than would be necessary for other pairs. In this paper, borrowing ideas from the study of geometrodynamics, we introduce a mathematical formalism to define a geometric distance between stellar models, to provide a quantitative meaning to this notion of closeness. In particular, it is known that the set of all metrics on a Riemannian manifold itself admits the structure of a Riemannian manifold (`configuration manifold'), which comes equipped with a canonical metric. By thinking of a stationary star as being a particular $3+1$ metric, the structure of which is determined through the Tolman-Oppenheimer-Volkoff relations and their generalisations, points on a suitably restricted configuration manifold can be thought of as representing different stars, and distances between these points can be computed. We develop the necessary mathematical machinery to build the configuration manifold of neutron star models, and provide some worked examples to illustrate how this space might be used in future studies of stellar structure.

In general relativity, the Kerr metric uniquely represents the geometry surrounding an isolated, rotating black hole. An identification of significant non-Kerr features in some astrophysical source would then provide a `smoking-gun' for the break-down of general relativity in the strong-field regime. On the other hand, Kerr black holes are common to many other theories of gravity, and thus a validation of the Kerr metric does not necessarily favour general relativity amongst all possibilities. The nature of gravitational perturbations will however differ between different theories of gravity. Future precision tests involving gravitational waves from oscillating black holes, such as identifications of the quasi-normal mode spectrum from ring-down, will thus be able to probe the underlying theory, even if the object is Kerr. Here, we write down the equations governing metric perturbations of a Kerr black hole in $f(R)$ gravity in a form that is more conducive to numerical study.

Fast radio bursts are millisecond-duration radio pulses of extragalactic origin. A recent statistical analysis has found that the burst energetics of the repeating source FRB 121102 follow a power-law, with an exponent that is curiously consistent with the Gutenberg-Richter law for earthquakes. This hints that repeat-bursters may be compact objects undergoing violent tectonic activity. For young magnetars, possessing crustal magnetic fields which are both strong ($B \gtrsim 10^{15}$ G) and highly multipolar, Hall drift can instigate significant field rearrangements even on $\lesssim$ century long timescales. This reconfiguration generates zones of magnetic stress throughout the outer layers of the star, potentially strong enough to facilitate frequent crustal failures. In this paper, assuming a quake scenario, we show how the crustal field evolution, which determines the resulting fracture geometries, can be tied to burst properties. Highly anisotropic stresses are generated by the rapid evolution of multipolar fields, implying that small, localised cracks can occur sporadically throughout the crust during the Hall evolution. Each of these shallow fractures may release bursts of energy, consistent in magnitude with those seen in the repeating sources FRB 121102 and FRB 180814.J0422+73.

Modified theories of gravity are often built such that they contain general relativity as a limiting case. This inclusion property implies that the Kerr metric is common to many families of theories. For example, all analytic $f(R)$ theories with vanishing constant term admit the Kerr solution. In any given theory, however, the response of the gravitational field to astrophysical disturbances is tied to the structure of the field equations. As such, even if black holes are Kerr, the underlying theory can, in principle, be probed through gravitational distortions. In this paper, we study linear perturbations of a Kerr black hole in $f(R)$ gravity using the Newman-Penrose formalism. We show that, as in general relativity, the equations governing the perturbed metric, which depend on the quadratic term of the function $f$, completely decouple.

A magnetically confined mountain on the surface of an accreting neutron star simultaneously reduces the global magnetic dipole moment through magnetic burial and generates a mass quadrupole moment, which emits gravitational radiation. Previous mountain models have been calculated for idealized isothermal and adiabatic equations of state. Here these models are generalised to include non-zero, finite thermal conduction. Grad-Shafranov equilibria for three representative, polytropic equations of state are evolved over many conduction time-scales with the magnetohydrodynamic solver PLUTO. It is shown that conduction facilitates the flow of matter towards the pole. Consequently the buried magnetic field is partially resurrected starting from an initially polytropic Grad-Shafranov equilibrium. The poleward mass current makes the star more prolate, marginally increasing its detectability as a gravitational wave source, though to an extent which is likely to be subordinate to other mountain physics. Thermal currents also generate filamentary hot spots $(\gtrsim 10^{8} \text{ K})$ in the mountain, especially near the pole where the heat flux is largest, with implications for type I X-ray bursts.

Some modified theories of gravity are known to predict monopolar, in addition to the usual quadrupolar and beyond, gravitational radiation in the form of `breathing' modes. For the same reason that octupole and higher-multipole terms often contribute negligibly to the overall wave strain, monopole terms tend to dominate. We investigate both monopolar and quadrupolar continuous gravitational radiation from neutron stars deformed through internal magnetic stresses. We adopt the Parameterised-Post-Newtonian formalism to write down equations describing the leading-order stellar properties in a theory-independent way, and derive some exact solutions for stars with mixed poloidal-toroidal magnetic fields. We then turn to the specific case of scalar-tensor theories to demonstrate how observational upper limits on the gravitational-wave luminosity of certain neutron stars may be used to place constraints on modified gravity parameters, most notably the Eddington parameter $\gamma$. For conservative, purely poloidal models with characteristic field strength given by the spindown minimum, upper limits for the Vela pulsar yield $1-\gamma \lesssim 4.2 \times 10^{-3}$. For models containing a strong toroidal field housing $\sim 99\%$ of the internal magnetic energy, we obtain the bound $1- \gamma \lesssim 8.0 \times 10^{-7}$. This latter bound is an order of magnitude tighter than those obtained from current Solar system experiments, though applies to the strong-field regime.

The oscillation spectrum of a perturbed neutron star is intimately related to the physical properties of the star, such as the equation of state. Observing pulsating neutron stars therefore allows one to place constraints on these physical properties. However, it is not obvious exactly how much can be learnt from such measurements. If we observe for long enough, and precisely enough, is it possible to learn everything about the star? A classical result in the theory of spectral geometry states that one cannot uniquely `hear the shape of a drum'. More formally, it is known that an eigenfrequency spectrum may not uniquely correspond to a particular geometry; some `drums' may be indistinguishable from a normal-mode perspective. In contrast, we show that the drum result does not extend to perturbations of simple neutron stars within general relativity -- in the case of axial (toroidal) perturbations of static, perfect fluid stars, a quasi-normal mode spectrum uniquely corresponds to a stellar profile. We show in this paper that it is not possible for two neutron stars, with distinct fluid profiles, to oscillate in an identical manner. This result has the information-theoretic consequence that gravitational waves completely encode the properties of any given oscillating star: unique identifications are possible in the limit of perfect measurement.

Some exact, nonlinear, vacuum gravitational wave solutions are derived for certain polynomial $f(R)$ gravities. We show that the boundaries of the gravitational domain of dependence, associated with events in polynomial $f(R)$ gravity, are not null as they are in general relativity. The implication is that electromagnetic and gravitational causality separate into distinct notions in modified gravity, which may have observable astrophysical consequences. The linear theory predicts that tachyonic instabilities occur, when the quadratic coefficient $a_{2}$ of the Taylor expansion of $f(R)$ is negative, while the exact, nonlinear, cylindrical wave solutions presented here can be superluminal for all values of $a_{2}$. Anisotropic solutions are found, whose wave-fronts trace out time- or space-like hypersurfaces with complicated geometric properties. We show that the solutions exist in $f(R)$ theories that are consistent with Solar System and pulsar timing experiments.

The Ernst formulation of the Einstein equations is generalised to accommodate $f(R)$ theories of gravity. It is shown that, as in general relativity, the axisymmetric $f(R)$ field equations for a vacuum spacetime that is either stationary or cylindrically symmetric reduce to a single, non-linear differential equation for a complex-valued scalar function. As a worked example, we apply the generalised Ernst equations to derive a $f(R)$ generalisation of the Zipoy-Voorhees metric, which may be used to describe the gravitational field outside of an ellipsoidal neutron star. We also apply the theory to investigate the phase speed of large-amplitude gravitational waves in $f(R)$ gravity in the context of soliton-like solutions that display shock-wave behaviour across the causal boundary.

A precondition for the radio emission of pulsars is the existence of strong, small-scale magnetic field structures (`magnetic spots') in the polar cap region. Their creation can proceed via crustal Hall drift out of two qualitatively and quantitatively different initial magnetic field configurations: a field confined completely to the crust and another which penetrates the whole star. The aim of this study is to explore whether these magnetic structures in the crust can deform the star sufficiently to make it an observable source of gravitational waves. We model the evolution of these field configurations, which can develop, within $\sim 10^4$ -- $10^5$ yr, magnetic spots with local surface field strengths $\sim 10^{14}$ G maintained over $\gtrsim 10^6$ yr. Deformations caused by the magnetic forces are calculated. We show that, under favourable initial conditions, a star undergoing crustal Hall drift can have ellipticity $\epsilon\sim 10^{-6}$, even with sub-magnetar polar field strengths, after $\sim 10^5$ yr. A pulsar rotating at $\sim 10^2$ Hz with such $\epsilon$ is a promising gravitational-wave source candidate. Since such large deformations can be caused only by a particular magnetic field configuration that penetrates the whole star and whose maximum magnetic energy is concentrated in the outer core region, gravitational wave emission observed from radio pulsars can thus inform us about the internal field structures of young neutron stars.

We construct multipole moments for stationary, asymptotically flat, spacetime solutions to higher-order curvature theories of gravity. The moments are defined using $3+1$ techniques involving timelike Killing vector constructions as in the classic papers by Geroch and Hansen. Using the fact that the Kerr-Newman metric is a vacuum solution to a particular class of $f(R)$ theories of gravity, we compute all its moments, and find that they admit recurrence relations similar to those for the Kerr solution in general relativity. It has been proposed previously that modelling the measured frequencies of quasi-periodic oscillations from galactic microquasars enables experimental tests of the no-hair theorem. We explore the possibility that, even if the no-hair relation is found to break down in the context of general relativity, there may be an $f(R)$ counterpart that is preserved. We apply the results to the microquasars GRS $1915$+$105$ and GRO J$1655$-$40$ using the diskoseismology and kinematic resonance models, and constrain the spins and `charges' [which are not really electric charges in the $f(R)$ context] of their black holes.

Quantum mechanical corrections to the hydromagnetic force balance equation, derived from the microscopic Schr\"{o}dinger-Pauli theory of quantum plasmas, modify the equilibrium structure and hence the mass quadrupole moment of a neutron star. It is shown here that the dominant effect --- spin paramagnetism --- is most significant in a magnetar, where one typically has $\mu_{B}|\boldsymbol{B}|\gtrsim k_B T_e$, where $\mu_{B}$ is the Bohr magneton, $\boldsymbol{B}$ is the magnetic field, and $T_e$ is the electron temperature. The spin paramagnetic deformation of a nonbarotropic magnetar with a linked poloidal-toroidal magnetic field is calculated to be up to ${{\sim 10}}$ times greater than the deformation caused solely by the Lorentz force. It depends on the degree of Pauli blocking by conduction electrons and the propensity to form magnetic domains, processes which are incompletely modelled at magnetar field strengths. The star becomes more oblate, as the toroidal field component strengthens. The result implies that existing classical predictions underestimate the maximum strength of the gravitational wave signal from rapidly spinning magnetars at birth. Turning the argument around, future gravitational-wave upper limits of increasing sensitivity will place ever-stricter constraints on the physics of Pauli blocking and magnetic domain formation under magnetar conditions.

The sudden spin-down event ('anti-glitch') observed in AXP 1E 2259+586 on 2012 April 21 was arguably caused by a decay of its internal toroidal magnetic field component, which turns a stable prolate configuration into an unstable one. We refine previous models of this process by modelling the star's magnetic field self-consistently as a 'twisted torus' configuration in non-barotropic equilibrium (which allows us to explore a greater range of equilibrium configurations). It is shown that, if the star's magnetic field is purely dipolar, the change in the toroidal field strength required to produce an anti-glitch of the observed size can be ~ 10 times larger than previously calculated. If the star has a quadrupolar magnetic field component, then an anti-glitch of similar magnitude can be produced via a decay of the quadrupole component, in addition to a decay of the toroidal component. We show that, if the quadrupole component decays, the minimum initial toroidal field strength and the change in toroidal field strength needed to produce the observed anti-glitch are lower than in the pure dipole twisted torus. In addition, we predict the maximum anti-glitch sizes, assuming that they are caused by a change in ellipticity, in four glitching magnetars and discuss the implications for energetics of accompanying X-ray bursts.

A recipe is presented to construct an analytic, self-consistent model of a non-barotropic neutron star with a poloidal-toroidal field of arbitrary multipole order, whose toroidal component is confined in a torus around the neutral curve inside the star, as in numerical simulations of twisted tori. The recipe takes advantage of magnetic-field-aligned coordinates to ensure continuity of the mass density at the surface of the torus. The density perturbation and ellipticity of such a star are calculated in general and for the special case of a mixed dipole-quadrupole field as a worked example. The calculation generalises previous work restricted to dipolar, poloidal-toroidal and multipolar, poloidal-only configurations. The results are applied, as an example, to magnetars whose observations (e.g., spectral features and pulse modulation) indicate that the internal magnetic fields may be at least one order of magnitude stronger than the external fields, as inferred from their spin downs, and are not purely dipolar.

We present some results regarding metric perturbations in general relativity and other metric theories of gravity. In particular, using the Newman Penrose variables, we write down and discuss the equations which govern tensor field perturbations of ranks $0, \pm 1$ and $\pm 2$ (scalar,vector,tensor bosonic perturbations) over certain space-times that admit specific background metrics expressible in isotropic coordinates. Armed with these equations, we are able to develop the Hadamard series which can be associated with the fundamental solution of the equations, wherein we introduce an inhomogeneous singularity at the physical space-time point of the perturbing particle. The Hadamard series is then used to generate closed form solutions by making choices for an appropriate ansatz solution. In particular, we solve for the spin-weighted electrostatic potential for the Reissner-Nordstrom black hole and for the fully dynamical potential for the Friedmann-Robertson-Walker cosmological solution.