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.

We present an extension to the special relativistic, ideal magnetohydrodynamics (MHD) equations, designed to capture effects due to resistivity. The extension takes the simple form of an additional source term which, when implemented numerically, is shown to emulate the behaviour produced by a fully resistive MHD description for a range of initial data. The extension is developed from first principle arguments, and thus requires no fine tuning of parameters, meaning it can be applied to a wide range of dynamical systems. Furthermore, our extension does not suffer from the same stiffness issues arising in resistive MHD, and thus can be evolved quickly using explicit methods, with performance benefits of roughly an order of magnitude compared to current methods.

Type I X-ray bursts are thermonuclear burning events which occur on the surfaces of accreting neutron stars. Burning begins in a localised spot in the star's ocean layer before propagating across the entire surface as a deflagration. On the scale of the entire star, the burning front can be thought of as discontinuity. To model this, we investigated the reactive Riemann problem for relativistic deflagrations and detonations and developed a numerical solver. Unlike for the Newtonian Riemann problem, where only the velocity perpendicular to the interface is relevant, in the relativistic case the tangential velocity becomes coupled through the Lorentz factor and can alter the waves present in the solution. We investigated whether a fast tangential velocity may be able to cause a deflagration wave to transition to a detonation. We found that such a transition is possible, but only for tangential velocities that are a significant fraction of the speed of light or for systems already on the verge of transitioning. Consequently, it is highly unlikely that this transition would occur for a burning front in a neutron star ocean without significant contributions from additional multidimensional effects.

pyro is a Python-based simulation framework designed for ease of implementation and exploration of hydrodynamics methods. It is built in a object-oriented fashion, allowing for the reuse of the core components and fast prototyping of new methods.

We discuss the challenges of modeling X-ray bursts in multi-dimensions, review the different calculations done to date, and discuss our new set of ongoing simulations. We also describe algorithmic improvements that may help in the future to offset some of the expense of these simulations, and describe what may be possible with exascale computing.

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.

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.

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

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

We perform 3+1 general relativistic simulations of rotating core collapse in the context of the collapsar model for long gamma-ray bursts. We employ a realistic progenitor, rotation based on results of stellar evolution calculations, and a simplified equation of state. Our simulations track self-consistently collapse, bounce, the postbounce phase, black hole formation, and the subsequent early hyperaccretion phase. We extract gravitational waves from the spacetime curvature and identify a unique gravitational wave signature associated with the early phase of collapsar formation.

We investigate the collapse of differentially rotating supermassive stars (SMSs) by means of 3+1 hydrodynamic simulations in general relativity. We particularly focus on the onset of collapse to understand the final outcome of collapsing SMSs. We find that the estimated ratio of the mass between the black hole (BH) and the surrounding disk from the equilibrium star is roughly the same as the results from numerical simulation. This suggests that the picture of axisymmetric collapse is adequate, in the absence of nonaxisymmetric instabilities, to illustrate the final state of the collapse. We also find that quasi-periodic gravitational waves continue to be emitted after the quasinormal mode frequency has decayed. We furthermore have found that when the newly formed BH is almost extreme Kerr, the amplitude of the quasi-periodic oscillation is enhanced during the late stages of the evolution. Geometrical features, shock waves, and instabilities of the fluid are suggested as a cause of this amplification behaviour. This alternative scenario for the collapse of differentially rotating SMSs might be observable by LISA.

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 results from the first 2D and 3D simulations of the collapse of rotating stellar iron cores in general relativity employing a finite-temperature equation of state and an approximate treatment of deleptonization during collapse. We compare fully nonlinear and conformally flat spacetime evolution methods and find that the conformally flat treatment is sufficiently accurate for the core-collapse supernova problem. We focus on the gravitational wave (GW) emission from rotating collapse, core bounce, and early postbounce phases. Our results indicate that the GW signature of these phases is much more generic than previously estimated. In addition, we track the growth of a nonaxisymmetric instability of dominant m = 1 character in one of our models that leads to prolonged narrow-band GW emission at ~930 Hz over several tens of milliseconds.

We investigate the role of rotational instabilities in the context of black hole formation in relativistic stars. In addition to the standard scenario - an axially symmetric dynamical instability forming a horizon at the star's center - the recently found low-$T/|W|$ instabilities are shown to lead to fragmentation and off-center horizon formation in differentially rotating stars. This process might be an alternative pathway to produce SMBHs from supermassive stars with inefficient angular momentum transport.

We present results from the first 2+1 and 3+1 simulations of the collapse of rotating stellar iron cores in general relativity employing a finite-temperature equation of state and an approximate treatment of deleptonization during collapse. We compare full 3+1 and conformally-flat spacetime evolution methods and find that the conformally-flat treatment is sufficiently accurate for the core-collapse supernova problem. We focus on the gravitational wave (GW) emission from rotating collapse, core bounce, and early postbounce phases. Our results indicate that the GW signature of these phases is much more generic than previously estimated. In addition, we track the growth of a nonaxisymmetric instability of dominant m = 1 character in two of our models that leads to prolonged narrow-band GW emission at ~ 920-930 Hz over several tens of milliseconds.

In a recent publication, we have demonstrated that differentially rotating stars admit new channels of black hole formation via fragmentation instabilities. Since a higher order instability of this kind could potentially transform a differentially rotating supermassive star into a multiple black hole system embedded in a massive accretion disk, we investigate the dependence of the instability on parameters of the equilibrium model. We find that many of the models constructed exhibit non-axisymmetric instabilities with corotation points, even for low values of T/|W|, which lead to a fission of the stars into one, two or three fragments, depending on the initial perturbation. At least in the models selected here, an m=1 mode becomes unstable at lower values of T/|W|, which would seem to favor a scenario where one black hole with a massive accretion disk forms. In this case, we have gained evidence that low values of compactness of the initial model can lead to a stabilization of the resulting fragment, thus preventing black hole formation in this scenario.

We present the first calculation of gravitational wave emission produced in the gravitational collapse of uniformly rotating neutron stars to black holes in fully three-dimensional simulations. The initial stellar models are relativistic polytropes which are dynamically unstable and with angular velocities ranging from slow rotation to the mass-shedding limit. An essential aspect of these simulations is the use of progressive mesh-refinement techniques which allow to move the outer boundaries of the computational domain to regions where gravitational radiation attains its asymptotic form. The waveforms have been extracted using a gauge-invariant approach in which the numerical spacetime is matched with the non-spherical perturbations of a Schwarzschild spacetime. Overall, the results indicate that the waveforms have features related to the properties of the initial stellar models (in terms of their w-mode oscillations) and of the newly produced rotating black holes (in terms of their quasi-normal modes). While our waveforms are in good qualitative agreement with those computed by Stark and Piran in two-dimensional simulations, our amplitudes are about one order of magnitude smaller and this difference is mostly likely due to our less severe pressure reduction. For a neutron star rotating uniformly near mass-shedding and collapsing at 10 kpc, the signal-to-noise ratio computed uniquely from the burst is S/N ~ 0.25, but this grows to be S/N <~ 4 in the case of LIGO II.

We investigate new paths to black hole formation by considering the general relativistic evolution of a differentially rotating polytrope with toroidal shape. We find that this polytrope is unstable to nonaxisymmetric modes, which leads to a fragmentation into self-gravitating, collapsing components. In the case of one such fragment, we apply a simplified adaptive mesh refinement technique to follow the evolution to the formation of an apparent horizon centered on the fragment. This is the first study of the one-armed instability in full general relativity.

We present a new three-dimensional fully general-relativistic hydrodynamics code using high-resolution shock-capturing techniques and a conformal traceless formulation of the Einstein equations. Besides presenting a thorough set of tests which the code passes with very high accuracy, we discuss its application to the study of the gravitational collapse of uniformly rotating neutron stars to Kerr black holes. The initial stellar models are modelled as relativistic polytropes which are either secularly or dynamically unstable and with angular velocities which range from slow rotation to the mass-shedding limit. We investigate the gravitational collapse by carefully studying not only the dynamics of the matter, but also that of the trapped surfaces, i.e. of both the apparent and event horizons formed during the collapse. The use of these surfaces, together with the dynamical horizon framework, allows for a precise measurement of the black-hole mass and spin. The ability to successfully perform these simulations for sufficiently long times relies on excising a region of the computational domain which includes the singularity and is within the apparent horizon. The dynamics of the collapsing matter is strongly influenced by the initial amount of angular momentum in the progenitor star and, for initial models with sufficiently high angular velocities, the collapse can lead to the formation of an unstable disc in differential rotation. All the simulations performed with uniformly rotating initial data and a polytropic or ideal-fluid equation of state show no evidence of shocks or of the presence of matter on stable orbits outside the black hole.