Extreme mass ratio inspirals, in which a stellar-mass object merges with a supermassive black hole, are prime sources for space-based gravitational wave detectors because they will facilitate tests of strong gravity and probe the spacetime around rotating compact objects. In the last few years of such inspirals, the total phase is in the millions of radians and details of the waveforms are sensitive to small perturbations. We show that one potentially detectable perturbation is the presence of a second supermassive black hole within a few tenths of a parsec. The acceleration produced by the perturber on the extreme mass-ratio system produces a steady drift that causes the waveform to deviate systematically from that of an isolated system. If the perturber is a few tenths of a parsec from the extreme-mass ratio system (plausible in as many as a few percent of cases) higher derivatives of motion might also be detectable. In that case, the mass and distance of the perturber can be derived independently, which would allow a new probe of merger dynamics.

GW Notes was born from the need for a journal where the distinct communities involved in gravitation wave research might gather. While these three communities - astrophysics, general relativity and data analysis - have made significant collaborative progress over recent years, we believe that it is indispensable to future advancement that they draw closer, and that they speak a common idiom. In this 6th GW Note (since we started numbering with #0), we present the work of Jonathan Thornburg, who has been fully-refereed, on the Capra research programme for capture of small compact objects by massive black holes.

Suppose a small compact object (black hole or neutron star) of mass $m$ orbits a large black hole of mass $M \gg m$. This system emits gravitational waves (GWs) that have a radiation-reaction effect on the particle's motion. EMRIs (extreme--mass-ratio inspirals) of this type will be important GW sources for LISA; LISA's data analysis will require highly accurate EMRI GW templates. In this article I outline the "Capra" research program to try to model EMRIs and calculate their GWs \textit{ab initio}, assuming only that $m \ll M$ and that the Einstein equations hold. Here we treat the EMRI spacetime as a perturbation of the large black hole's "background" (Schwarzschild or Kerr) spacetime and use the methods of black-hole perturbation theory, expanding in the small parameter $m/M$. The small body's motion can be described either as the result of a radiation-reaction "self-force" acting in the background spacetime or as geodesic motion in a perturbed spacetime. Several different lines of reasoning lead to the (same) basic $\O(m/M)$ "MiSaTaQuWa" equations of motion for the particle. Surprisingly, for a nonlinear field theory such as general relativity, modelling the small body as a point particle works well. The particle's own field is singular along the particle worldline so it's difficult to formulate a meaningful "perturbation" theory or equations of motion there. I discuss "mode-sum" and "puncture-function" regularization schemes that resolve this difficulty and allow practical self-force calculations, and I outline an important recent calculation of this type. Most Capra research to date has used 1st order perturbation theory. To obtain the very high accuracies needed to fully exploit LISA's observations of the strongest EMRIs, 2nd order perturbation theory will probably be needed.

Gravitational waveforms from the inspiral and ring-down stages of the binary black hole coalescences can be modelled accurately by approximation/perturbation techniques in general relativity. Recent progress in numerical relativity has enabled us to model also the non-perturbative merger phase of the binary black-hole coalescence problem. This enables us to \emph{coherently} search for all three stages of the coalescence of non-spinning binary black holes using a single template bank. Taking our motivation from these results, we propose a family of template waveforms which can model the inspiral, merger, and ring-down stages of the coalescence of non-spinning binary black holes that follow quasi-circular inspiral. This two-dimensional template family is explicitly parametrized by the physical parameters of the binary. We show that the template family is not only \emph{effectual} in detecting the signals from black hole coalescences, but also \emph{faithful} in estimating the parameters of the binary. We compare the sensitivity of a search (in the context of different ground-based interferometers) using all three stages of the black hole coalescence with other template-based searches which look for individual stages separately. We find that the proposed search is significantly more sensitive than other template-based searches for a substantial mass-range, potentially bringing about remarkable improvement in the event-rate of ground-based interferometers. As part of this work, we also prescribe a general procedure to construct interpolated template banks using non-spinning black hole waveforms produced by numerical relativity.

Recent progress in numerical relativity has enabled us to model the non-perturbative merger phase of the binary black-hole coalescence problem. Based on these results, we propose a phenomenological family of waveforms which can model the inspiral, merger, and ring-down stages of black hole coalescence. We also construct a template bank using this family of waveforms and discuss its implementation in the search for signatures of gravitational waves produced by black-hole coalescences in the data of ground-based interferometers. This template bank might enable us to extend the present inspiral searches to higher-mass binary black-hole systems, i.e., systems with total mass greater than about 80 solar masses, thereby increasing the reach of the current generation of ground-based detectors.

This paper presents a detailed discussion of the ``Newton's method'' algorithm for finding apparent horizons in 3+1 numerical relativity. We describe a method for computing the Jacobian matrix of the finite differenced $H(h)$ function by symbolically differentiating the finite difference equations, giving the Jacobian elements directly in terms of the finite difference molecule coefficients used in computing $H(h)$. Assuming the finite differencing scheme commutes with linearization, we show how the Jacobian elements may be computed by first linearizing the continuum $H(h)$ equations, then finite differencing the linearized (continuum) equations. We find this symbolic differentiation method of computing the $H(h)$ Jacobian to be {\em much} more efficient than the usual numerical perturbation method, and also much easier to implement than is commonly thought. When solving the discrete $H(h) = 0$ equations, we find that Newton's method generally converges very rapidly. However, if the initial guess for the horizon position contains significant high-spatial-frequency error components, Newton's method has a small (poor) radius of convergence. This is {\em not} an artifact of insufficient resolution in the finite difference grid; rather, it appears to be caused by a strong nonlinearity in the continuum $H(h)$ function for high-spatial-frequency error components in $h$. Robust variants of Newton's method can boost the radius of convergence by O(1) factors, but the underlying nonlinearity remains, and appears to worsen rapidly with increasing initial-guess-error spatial frequency. Using 4th~order finite differencing, we find typical accuracies for computed horizon positions in the $10^{-5}$ range for $\Delta\theta = \frac{\pi/2}{50}$.