We compute adiabatic waveforms for extreme mass-ratio inspirals (EMRIs) by “stitching” together a long inspiral waveform from a sequence of waveform snapshots, each of which corresponds to a particular geodesic orbit. We show that the complicated total waveform can be regarded as a sum of “voices.” Each voice evolves in a simple way on long timescales, a property which can be exploited to efficiently produce waveform models that faithfully encode the properties of EMRI systems. We look at examples for a range of different orbital geometries: spherical orbits, equatorial eccentric orbits, and one example of generic (inclined and eccentric) orbits. To our knowledge, this is the first calculation of a generic EMRI waveform that uses strong-field radiation reaction. We examine waveforms in both the time and frequency domains. Although EMRIs evolve slowly enough that the stationary phase approximation (SPA) to the Fourier transform is valid, the SPA calculation must be done to higher order for some voices, since their instantaneous frequency can change from chirping forward ($\dot{f}>0$) to chirping backward ($\dot{f}<0$). The approach we develop can eventually be extended to more complete EMRI waveform models—for example, to include effects neglected by the adiabatic approximation, such as the conservative self-force and spin-curvature coupling.

We develop new strategies to build numerical relativity surrogate models for eccentric binary black hole systems, which are expected to play an increasingly important role in current and future gravitational-wave detectors. We introduce a new surrogate waveform model, NRSur2dq1Ecc, using 47 nonspinning, equal-mass waveforms with eccentricities up to 0.2 when measured at a reference time of $5500M$ before merger. This is the first waveform model that is directly trained on eccentric numerical relativity simulations and does not require that the binary circularizes before merger. The model includes the (2,2), (3,2), and (4,4) spin-weighted spherical harmonic modes. We also build a final black hole model, NRSur2dq1EccRemnant, which models the mass, and spin of the remnant black hole. We show that our waveform model can accurately predict numerical relativity waveforms with mismatches $\approx {10}^{-3}$, while the remnant model can recover the final mass and dimensionless spin with absolute errors smaller than $\approx 5\times {10}^{-4}M$ and $\approx 2\times {10}^{-3}$ respectively. We demonstrate that the waveform model can also recover subtle effects like mode mixing in the ringdown signal without any special ad hoc modeling steps. Finally, we show that despite being trained only on equal-mass binaries, NRSur2dq1Ecc can be reasonably extended up to mass ratio $q\approx 3$ with mismatches $\simeq {10}^{-2}$ for eccentricities smaller than $\sim 0.05$ as measured at a reference time of $2000M$ before merger. The methods developed here should prove useful in the building of future eccentric surrogate models over larger regions of the parameter space.

The relaxation of a distorted black hole to its final state provides important tests of general relativity within the reach of current and upcoming gravitational wave facilities. In black hole perturbation theory, this phase consists of a simple linear superposition of exponentially damped sinusoids (the quasinormal modes) and of a power-law tail. How many quasinormal modes are necessary to describe waveforms with a prescribed precision? What error do we incur by only including quasinormal modes, and not tails? What other systematic effects are present in current state-of-the-art numerical waveforms? These issues, which are basic to testing fundamental physics with distorted black holes, have hardly been addressed in the literature. We use numerical relativity waveforms and accurate evolutions within black hole perturbation theory to provide some answers. We show that (i) a determination of the fundamental $l=m=2$ quasinormal frequencies and damping times to within 1% or better requires the inclusion of at least the first overtone, and preferably of the first two or three overtones; (ii) a determination of the black hole mass and spin with precision better than 1% requires the inclusion of at least two quasinormal modes for any given angular harmonic mode ($\ell$, $m$). We also improve on previous estimates and fits for the ringdown energy radiated in the various multipoles. These results are important to quantify theoretical (as opposed to instrumental) limits in parameter estimation accuracy and tests of general relativity allowed by ringdown measurements with high signal-to-noise ratio gravitational wave detectors.

A general class of loop quantizations for anisotropic models is introduced and discussed, which enhances loop quantum cosmology by relevant features seen in inhomogeneous situations. The main new effect is an underlying lattice which is being refined during dynamical changes of the volume. In general, this leads to a new feature of dynamical difference equations which may not have constant step-size, posing new mathematical problems. It is discussed how such models can be evaluated and what lattice refinements imply for semiclassical behavior. Two detailed examples illustrate that stability conditions can put strong constraints on suitable refinement models, even in the absence of a fundamental Hamiltonian which defines changes of the underlying lattice. Thus, a large class of consistency tests of loop quantum gravity becomes available. In this context, it will also be seen that quantum corrections due to inverse powers of metric components in a constraint are much larger than they appeared recently in more special treatments of isotropic, free scalar models where they were artificially suppressed.

We describe the hyperboloidal compactification for Teukolsky equations in Kerr spacetime. We include null infinity on the numerical grid by attaching a hyperboloidal layer to a compact domain surrounding the rotating black hole and the orbit of an inspiralling point particle. This technique allows us to study, for the first time, gravitational waveforms from large- and extreme-mass-ratio inspirals in Kerr spacetime extracted at null infinity. Tests and comparisons of our results with previous calculations establish the accuracy and efficiency of the hyperboloidal layer method.

A classical thought-experiment to destroy black holes was envisaged by Wald in 1974: it consists of throwing particles with large angular momentum into an extremal black hole, checking whether their capture can overspin the black hole past the extremal limit and create a naked singularity. Wald showed that in the test-particle limit, particles that would be otherwise capable of producing naked singularities are simply scattered. Recently, Jacobson and Sotiriou showed that if one considers instead a black hole that is almost, but not exactly extremal, then in the absence of backreaction effects particle capture could indeed overspin the spacetime above the Kerr limit. Here we analyze backreaction effects and show that for some of the trajectories giving rise to naked singularities, radiative effects can be neglected. However, for these orbits the conservative self-force is important, and seems to have the right sign to prevent the formation of naked singularities.