Gravitational wave signals from compact astrophysical sources such as those observed by LIGO and Virgo require a high-accuracy, theory-based waveform model for the analysis of the recorded signal. Current inspiral-merger-ringdown models are calibrated only up to moderate mass ratios, thereby limiting their applicability to signals from high-mass ratio binary systems. We present EMRISur1dq1e4, a reduced-order surrogate model for gravitational waveforms of 13,500M in duration and including several harmonic modes for non-spinning black hole binary systems with mass-ratios varying from 3 to 10,000 thus vastly expanding the parameter range beyond the current models. This surrogate model is trained on waveform data generated by point-particle black hole perturbation theory (ppBHPT) both for large mass-ratio and comparable mass-ratio binaries. We observe that the gravitational waveforms generated through a simple application of ppBHPT to the comparable mass-ratio cases agree remarkably (and surprisingly) well with those from full numerical relativity after a rescaling of the ppBHPT's total mass parameter. This observation and the EMRISur1dq1e4 surrogate model will enable data analysis studies in the high-mass ratio regime, including potential intermediate mass-ratio signals from LIGO/Virgo and extreme-mass ratio events of interest to the future space-based observatory LISA.

Using inspiral and plunge trajectories we construct with a generalized Ori-Thorne algorithm, we use a time-domain black hole perturbation theory code to compute the corresponding gravitational waves. The last cycles of these waveforms are a superposition of Kerr quasinormal modes. In this paper, we examine how the modes' excitations vary as a function of source parameters, such as the larger black hole's spin and the geometry of the smaller body's inspiral and plunge. We find that the mixture of quasinormal modes that characterize the final gravitational waves from a coalescence is entirely determined by the spin $a$ of the larger black hole, an angle $I$ which characterizes the misalignment of the orbital plane from the black hole's spin axis, a second angle $\theta_{\rm fin}$ which describes the location at which the small body crosses the black hole's event horizon, and the direction sgn$(\dot\theta_{\rm fin})$ of the body's final motion. If these large-mass-ratio results hold at less extreme mass ratios, then measuring multiple ringdown modes of binary black hole coalescence gravitational waves may provide important information about the source's binary properties, such as the misalignment of the orbit's angular momentum with black hole spin. This may be particularly useful for large mass binaries, for which the early inspiral waves are out of the detectors' most sensitive band.

The coalescence of two black holes generates gravitational waves that carry detailed information about the properties of those black holes and their binary configuration. The final coalescence cycles are in the form of a {\it ringdown}: a superposition of quasi-normal modes of the merged remnant black hole. Each mode has an oscillation frequency and decay time that in general relativity is determined by the remnant's mass and spin. Measuring the frequency and decay time of multiple modes makes it possible to measure the remnant's mass and spin, and to test the waves against the predictions of gravity theories. In this {\it Letter}, we show that the relative amplitudes of these modes encodes information about a binary's {\it geometry}. Focusing on the large mass-ratio limit, which provides a simple-to-use tool for effectively exploring parameter space, we demonstrate how a binary's geometry is encoded in the relative amplitudes of these modes, and how to parameterize the modes in this limit. Although more work is needed to assess how well this carries over to less extreme mass ratios, our results indicate that measuring multiple ringdown modes from coalescence may aid in measuring important source properties, such as the misalignment of its members' spins and orbit.

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

The recent detection of gravitational waves has generated interest in alternatives to the black hole interpretation of sources. One set of such alternatives involves a prediction of gravitational wave "echoes". We consider two aspects of possible echoes: First, general features of echoes coming from spacetime reflecting conditions. We find that the detailed nature of such echoes does not bear any clear relationship to quasi-normal frequencies. Second, we point out the pitfalls in the analysis of local reflecting "walls" near the horizon of rapidly rotating black holes.

Modelling of gravitational waves from binary black hole inspiral has played an important role in the recent observations of such signals. The late-stage ringdown phase of the gravitational waveform is often associated with the null particle orbit ("light ring") of the black hole spacetime. With simple models we show that this link between the light ring and spacetime ringing is based more on the history of specific models than on an actual constraining relationship. We also show, in particular, that a better understanding of the dissociation of the two may be relevant to the astrophysically interesting case of rotating (Kerr) black holes.

Modeling of gravitational waves (GWs) from binary black hole inspiral brings together early post-Newtonian waveforms and late quasinormal ringing waveforms. Attempts to bridge the two limits without recourse to numerical relativity involve predicting the time of the peak GW amplitude. This prediction will require solving the question of why the peak of the "source," i.e., the peak of the binary angular velocity, does not correspond to the peak of the GW amplitude. We show here that this offset can be understood as due to the existence two distinct components of the radiation: the "direct" radiation analogous to that in flat spacetime, and "scattered" radiation associated with curved spacetime. The time dependence of these two components, and of their relative phases determines the location of the peak amplitude. We use a highly simplified model to clarify the twocomponent nature of the source, then demonstrate that the explanation is valid also for an extreme mass ratio binary inspiral.

Using the effective-one-body (EOB) formalism and a time-domain Teukolsky code, we generate inspiral, merger, and ringdown waveforms in the small mass-ratio limit. We use EOB inspiral and plunge trajectories to build the Teukolsky equation source term, and compute full coalescence waveforms for a range of black hole spins. By comparing EOB waveforms that were recently developed for comparable mass binary black holes to these Teukolsky waveforms, we improve the EOB model for the (2,2), (2,1), (3,3), and (4,4) modes. Our results can be used to quickly and accurately extract useful information about merger waves for binaries with spin, and should be useful for improving analytic models of such binaries.

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 over-spin 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 over-spin the spacetime above the Kerr limit. Here we analyze back-reaction 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.

The thermonuclear explosion of a C/O white dwarf as a Type Ia supernova (SN Ia) generates a kinetic energy comparable to that released by a massive star during a SN II event. Current observations and theoretical models have established that SNe Ia are asymmetric, and therefore--like SNe II--potential sources of gravitational wave (GW) radiation. We perform the first detailed calculations of the GW emission for a SN Ia of any type within the single-degenerate channel. The gravitationally-confined detonation (GCD) mechanism predicts a strongly-polarized GW burst in the frequency band around 1 Hz. Third-generation spaceborne GW observatories currently in planning may be able to detect this predicted signal from SNe Ia at distances up to 1 Mpc. If observable, GWs may offer a direct probe into the first few seconds of the SNe Ia detonation.

Jacobson and Sotiriou showed that rotating black holes could be spun-up past the extremal limit by the capture of non-spinning test bodies, if one neglects radiative and self-force effects. This would represent a violation of the Cosmic Censorship Conjecture in four-dimensional, asymptotically flat spacetimes. We 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.

Spectacular breakthroughs in numerical relativity now make it possible to compute spacetime dynamics in almost complete generality, allowing us to model the coalescence and merger of binary black holes with essentially no approximations. The primary limitation of these calculations is now computational. In particular, it is difficult to model systems with large mass ratio and large spins, since one must accurately resolve the multiple lengthscales which play a role in such systems. Perturbation theory can play an important role in extending the reach of computational modeling for binary systems. In this paper, we present first results of a code which allows us to model the gravitational waves generated by the inspiral, merger, and ringdown of a binary system in which one member of the binary is much more massive than the other. This allows us to accurately calibrate binary dynamics in the large mass ratio regime. We focus in this analysis on the recoil imparted to the merged remnant by these waves. We closely examine the "antikick", an anti-phase cancellation of the recoil arising from the plunge and ringdown waves, described in detail by Schnittman et al. We find that, for orbits aligned with the black hole spin, the antikick grows as a function of spin. The total recoil is smallest for prograde coalescence into a rapidly rotating black hole, and largest for retrograde coalescence. Amusingly, this completely reverses the predicted trend for kick versus spin from analyses that only include inspiral information.

This is the second in a series of papers whose aim is to generate ``adiabatic'' gravitational waveforms from the inspiral of stellar-mass compact objects into massive black holes. In earlier work, we presented an accurate (2+1)D finite-difference time-domain code to solve the Teukolsky equation, which evolves curvature perturbations near rotating (Kerr) black holes. The key new ingredient there was a simple but accurate model of the singular source term based on a discrete representation of the Dirac-delta function and its derivatives. Our earlier work was intended as a proof of concept, using simple circular, equatorial geodesic orbits as a testbed. Such a source is effectively static, in that the smaller body remains at the same coordinate radius and orbital inclination over an orbit. (It of course moves through axial angle, but we separate that degree of freedom from the problem. Our numerical grid has only radial, polar, and time coordinates.) We now extend the time-domain code so that it can accommodate dynamic sources that move on a variety of physically interesting world lines. We validate the code with extensive comparison to frequency-domain waveforms for cases in which the source moves along generic (inclined and eccentric) bound geodesic orbits. We also demonstrate the ability of the time-domain code to accommodate sources moving on interesting non-geodesic worldlines. We do this by computing the waveform produced by a test mass following a ``kludged'' inspiral trajectory, made of bound geodesic segments driven toward merger by an approximate radiation loss formula.

We revisit the problem of the emission of gravitational waves from a test mass orbiting and thus perturbing a Kerr black hole. The source term of the Teukolsky perturbation equation contains a Dirac delta function which represents a point particle. We present a technique to effectively model the delta function and its derivatives using as few as four points on a numerical grid. The source term is then incorporated into a code that evolves the Teukolsky equation in the time domain as a (2+1) dimensional PDE. The waveforms and energy fluxes are extracted far from the black hole. Our comparisons with earlier work show an order of magnitude gain in performance (speed) and numerical errors less than 1% for a large fraction of parameter space. As a first application of this code, we analyze the effect of finite extraction radius on the energy fluxes. This paper is the first in a series whose goal is to develop adiabatic waveforms describing the inspiral of a small compact body into a massive Kerr black hole.

We study the Klein-Gordon equation for a massive scalar field in Kerr spacetime in the time-domain. We demonstrate that under conditions of super-radiance, the scalar field becomes unstable and its amplitude grows without bound. We also estimate the growth rate of this instability.

The late-time tails of a massive scalar field in the spacetime of black holes are studied numerically. Previous analytical results for a Schwarzschild black hole are confirmed: The late-time behavior of the field as recorded by a static observer is given by $\psi(t)\sim t^{-5/6}\sin [\omega (t)\times t]$, where $\omega(t)$ depends weakly on time. This result is carried over to the case of a Kerr black hole. In particular, it is found that the power-law index of -5/6 depends on neither the multipole mode $\ell$ nor on the spin rate of the black hole $a/M$. In all black hole spacetimes, massive scalar fields have the same late-time behavior irrespective of their initial data (i.e., angular distribution). Their late-time behavior is universal.

The presence of spiral structure in isolated galaxies is a problem that has only been partially explained by theoretical models. Because the rate and pattern of star formation in the disk must depend only on mechanisms internal to the disk, we may think of the spiral galaxy as a self-regulated system far from equilibrium. This paper uses this idea to look at a reaction-diffusion model for the formation of spiral structures in certain types of galaxies. In numerical runs of the model, spiral structure forms and persists over several revolutions of the disk, but eventually dies out.

Withdrawn; paper superceded by astro-ph/0105392

In this paper we consider the collision of spinning holes using first order perturbation theory of black holes (Teukolsky formalism). With these results (along with ones, we published in the past) one can predict the properties of the gravitational waves radiated from the late stage inspiral of two spinning, equal mass black holes. Also we note that the energy radiated by the head-on collision of two spinning holes with spins (that are equal and opposite) aligned along the common axis is more than the case in which the spins are perpendicular to the axis of the collision.

Using approximate techniques we study the final moments of the collision of two (individually non-spinnning) black holes which inspiral into each other. The approximation is based on treating the whole space-time as a single distorted black hole. We obtain estimates for the radiated energy, angular momentum and waveforms for the gravitational waves produced in such a collision. The results can be of interest for analyzing the data that will be forthcoming from gravitational wave interferometric detectors, like the LIGO, GEO, LISA, VIRGO and TAMA projects.

Using several approximations, we calculate an estimate of the gravitational radiation emitted when two equal mass black holes coalesce at the end of their binary inspiral. We find that about 1% of the mass energy of the pair will emerge as gravitational waves during the final ringdown and a negligible fraction of the angular momentum will be radiated.

We solve the problem of expressing the Weyl scalars $\psi $ that describe gravitational perturbations of a Kerr black hole in terms of Cauchy data. To do so we use geometrical identities (like the Gauss-Codazzi relations) as well as Einstein equations. We are able to explicitly express $\psi $ and $\partial _t\psi $ as functions only of the extrinsic curvature and the three-metric (and geometrical objects built out of it) of a generic spacelike slice of the spacetime. These results provide the link between initial data and $\psi $ to be evolved by the Teukolsky equation, and can be used to compute the gravitational radiation generated by two orbiting black holes in the close limit approximation. They can also be used to extract waveforms from spacetimes completely generated by numerical methods.