When an accretion disk falls prey to the runaway instability, a large portion of its mass is devoured by the black hole within a few dynamical times. Despite decades of effort, it is still unclear under what conditions such an instability can occur. The technically most advanced relativistic simulations to date were unable to find a clear sign for the onset of the instability. In this work, we present three-dimensional relativistic hydrodynamics simulations of accretion disks around black holes in dynamical space-time. We focus on the configurations that are expected to be particularly prone to the development of this instability. We demonstrate, for the first time, that the fully self-consistent general relativistic evolution does indeed produce a runaway instability.

An indirect, hybrid Monte Carlo discretization of general relativistic kinetic theory suitable for the development of numerical schemes for radiation transport is presented. The discretization is based on surface flux estimators obtained from a local decomposition of the distribution function, and can handle optically thick regions by means of formal solutions within each cell. Furthermore, the scheme is designed for parallel implementation, and it admits the use of adaptive techniques by virtue of leaving all probability density functions unspecified. Some considerations for numerical uses of the scheme are discussed.

We perform nonlinear general relativistic ideal magnetohydrodynamic simulations of poloidal magnetic fields in rotating polytropic neutron stars. We have three primary goals: i) to understand the nature of magnetohydrodynamic instabilities inherent to poloidal magnetic fields in non-rotating and rotating neutron stars, ii) to explore the possible space of stable equilibrium configurations and iii) to understand gravitational wave emissions caused by the catastrophic reconfiguration of magnetic fields associated with giant magnetar flares. Our key physical contributions can be summarized as follows: i) gravitational waves from f-modes caused by magnetar flares are unlikely to be detected in the current or near-future generation of gravitational waves observatories, ii) gravitational waves from Alfven waves propagating inside the neutron star are more likely candidates, although this interpretation relies on the unknown damping time of these modes, iii) any magnetic field equilibria derived from our simulations are characterized as non-axisymmetric, with approximately 65% of their magnetic energy in the poloidal field, iv) rotation acts to separate the timescales of different instabilities in our system, with the varicose mode playing a more major role due to a delayed kink instability and v) despite the slowing growth rate of the kink mode, it is always present in our simulations, even for models where the rotational period is of the same order as the Alfven timescale.

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.

Are giant flares in magnetars viable sources of gravitational radiation? Few theoretical studies have been concerned with this problem, with the small number using either highly idealized models or assuming a magnetic field orders of magnitude beyond what is supported by observations. We perform nonlinear general-relativistic magnetohydrodynamics simulations of large-scale hydromagnetic instabilities in magnetar models. We utilise these models to find gravitational wave emissions over a wide range of energies, from 10^40 to 10^47 erg. This allows us to derive a systematic relationship between the surface field strength and the gravitational wave strain, which we find to be highly nonlinear. In particular, for typical magnetar fields of a few times 10^15 G, we conclude that a direct observation of f-modes excited by global magnetic field reconfigurations is unlikely with present or near-future gravitational wave observatories, though we also discuss the possibility that modes in a low-frequency band up to 100 Hz could be sufficiently excited to be relevant for observation.

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.

We present Horizon, a new graphics processing unit (GPU)-accelerated code to solve the equations of general relativistic magnetohydrodynamics in a given spacetime. We evaluate the code in several test cases, including magnetized Riemann problems and rapidly rotating neutron stars, and measure the performance benefits of the GPU acceleration in comparison to our CPU-based code Thor. We find substantial performance gains in comparison to a quad-core CPU both in single- and double-precision accuracy, and discuss these findings in the context of future numerical modeling efforts.

Self-gravitating relativistic disks around black holes can form as transient structures in a number of astrophysical scenarios such as binary neutron star and black hole-neutron star coalescences, as well as the core-collapse of massive stars. We explore the stability of such disks against runaway and non-axisymmetric instabilities using three-dimensional hydrodynamics simulations in full general relativity using the THOR code. We model the disk matter using the ideal fluid approximation with a $\Gamma$-law equation of state with $\Gamma=4/3$. We explore three disk models around non-rotating black holes with disk-to-black hole mass ratios of 0.24, 0.17 and 0.11. Due to metric blending in our initial data, all of our initial models contain an initial axisymmetric perturbation which induces radial disk oscillations. Despite these oscillations, our models do not develop the runaway instability during the first several orbital periods. Instead, all of the models develop unstable non-axisymmetric modes on a dynamical timescale. We observe two distinct types of instabilities: the Papaloizou-Pringle and the so-called intermediate type instabilities. The development of the non-axisymmetric mode with azimuthal number m = 1 is accompanied by an outspiraling motion of the black hole, which significantly amplifies the growth rate of the m = 1 mode in some cases. Overall, our simulations show that the properties of the unstable non-axisymmetric modes in our disk models are qualitatively similar to those in Newtonian theory.

Time-division SQUID multiplexers are used in many applications that require exquisite control of systematic error. One potential source of systematic error is the pickup of external magnetic fields in the multiplexer. We present measurements of the field sensitivity figure of merit, effective area, for both the first stage and second stage SQUID amplifiers in three NIST SQUID multiplexer designs. These designs include a new variety with improved gradiometry that significantly reduces the effective area of both the first and second stage SQUID amplifiers.

Rapidly rotating neutron stars can be unstable to the gravitational-wave-driven CFS mechanism if they have a neutral point in the spectrum of nonaxisymmetric f-modes. We investigate the frequencies of these modes in two sequences of uniformly rotating polytropes using nonlinear simulations in full general relativity, determine the approximate locations of the neutral points, and derive limits on the observable frequency band available to the instability in these sequences. We find that general relativity enhances the detectability of a CFS-unstable neutron star substantially, both by widening the instability window and enlarging the band into the optimal range for interferometric detectors like LIGO, VIRGO, and GEO-600.

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

Many systems of interest in general relativistic astrophysics, including neutron stars, accreting compact objects in X-ray binaries and active galactic nuclei, core collapse, and collapsars, are assumed to be approximately spherically symmetric or axisymmetric. In Newtonian or fixed-background relativistic approximations it is common practice to use spherical polar coordinates for computational grids; however, these coordinates have singularities and are difficult to use in fully relativistic models. We present, in this series of papers, a numerical technique which is able to use effectively spherical grids by employing multiple patches. We provide detailed instructions on how to implement such a scheme, and present a number of code tests for the fixed background case, including an accretion torus around a black hole.

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