For a Schwarzchild black hole of mass $M$, we consider a test particle falling from rest at infinity and becoming trapped, at late time, on the unstable circular orbit of radius $r=4GM/c^2$. When the particle is endowed with a small mass, $\mu\ll M$, it experiences an effective gravitational self-force, whose conservative piece shifts the critical value of the angular momentum and the frequency of the asymptotic circular orbit away from their geodesic values. By directly integrating the self-force along the orbit (ignoring radiative dissipation), we numerically calculate these shifts to $O(\mu/M)$. Our numerical values are found to be in agreement with estimates first made within the Effective One Body formalism, and with predictions of the first law of black-hole-binary mechanics (as applied to the asymptotic circular orbit). Our calculation is based on a time-domain integration of the Lorenz-gauge perturbation equations, and it is a first such calculation for an unbound orbit. We tackle several technical difficulties specific to unbound orbits, illustrating how these may be handled in more general cases of unbound motion. Our method paves the way to calculations of the self-force along hyperbolic-type scattering orbits. Such orbits can probe the two-body potential down to the "light ring", and could thus supply strong-field calibration data for eccentricity-dependent terms in the Effective One Body model of merging binaries.

Black holes are unique among astrophysical sources: they are the simplest macroscopic objects in the Universe, and they are extraordinary in terms of their ability to convert energy into electromagnetic and gravitational radiation. Our capacity to probe their nature is limited by the sensitivity of our detectors. The LIGO/Virgo interferometers are the gravitational-wave equivalent of Galileo's telescope. The first few detections represent the beginning of a long journey of exploration. At the current pace of technological progress, it is reasonable to expect that the gravitational-wave detectors available in the 2035-2050s will be formidable tools to explore these fascinating objects in the cosmos, and space-based detectors with peak sensitivities in the mHz band represent one class of such tools. These detectors have a staggering discovery potential, and they will address fundamental open questions in physics and astronomy. Are astrophysical black holes adequately described by general relativity? Do we have empirical evidence for event horizons? Can black holes provide a glimpse into quantum gravity, or reveal a classical breakdown of Einstein's gravity? How and when did black holes form, and how do they grow? Are there new long-range interactions or fields in our universe, potentially related to dark matter and dark energy or a more fundamental description of gravitation? Precision tests of black hole spacetimes with mHz-band gravitational-wave detectors will probe general relativity and fundamental physics in previously inaccessible regimes, and allow us to address some of these fundamental issues in our current understanding of nature.

The grand challenges of contemporary fundamental physics---dark matter, dark energy, vacuum energy, inflation and early universe cosmology, singularities and the hierarchy problem---all involve gravity as a key component. And of all gravitational phenomena, black holes stand out in their elegant simplicity, while harbouring some of the most remarkable predictions of General Relativity: event horizons, singularities and ergoregions. The hitherto invisible landscape of the gravitational Universe is being unveiled before our eyes: the historical direct detection of gravitational waves by the LIGO-Virgo collaboration marks the dawn of a new era of scientific exploration. Gravitational-wave astronomy will allow us to test models of black hole formation, growth and evolution, as well as models of gravitational-wave generation and propagation. It will provide evidence for event horizons and ergoregions, test the theory of General Relativity itself, and may reveal the existence of new fundamental fields. The synthesis of these results has the potential to radically reshape our understanding of the cosmos and of the laws of Nature. The purpose of this work is to present a concise, yet comprehensive overview of the state of the art in the relevant fields of research, summarize important open problems, and lay out a roadmap for future progress.

[Abridged] This review surveys the theory of gravitational self-force in curved spacetime and its application to the gravitational two-body problem in the extreme-mass-ratio regime. We first lay the relevant formal foundation, describing the rigorous derivation of the equation of self-forced motion using matched asymptotic expansions and other ideas. We then review the progress that has been achieved in numerically calculating the self-force and its physical effects in the astrophysical scenario of a compact object inspiralling into a (rotating) massive black hole. We highlight the way in which, nowadays, self-force calculations make a fruitful contact with other approaches to the two-body problem and help inform an accurate universal model of binary black hole inspirals, valid across all mass ratios. We conclude with a summary of the state of the art, open problems and prospects. Our review is aimed at non-specialist readers and is for the most part self-contained and non-technical; only elementary-level acquaintance with General Relativity is assumed. Where useful, we draw on analogies with familiar concepts from Newtonian gravity or classical electrodynamics.

We present a new method for calculation of the gravitational self-force (GSF) in Kerr geometry, based on a time-domain reconstruction of the metric perturbation from curvature scalars. In this approach, the GSF is computed directly from a certain scalar-like self-potential that satisfies the time-domain Teukolsky equation on the Kerr background. The approach is computationally much cheaper than existing time-domain methods, which rely on a direct integration of the linearized Einstein's equations and are impaired by mode instabilities. At the same time, it retains the utility and flexibility of a time-domain treatment, allowing calculations for any type of orbit (including highly eccentric or unbound ones) and the possibility of self-consistently evolving the orbit under the effect of the GSF. Here we formulate our method, and present a first numerical application, for circular geodesic orbits in Schwarzschild geometry. We discuss further applications.

Vacuum perturbations of the Kerr metric can be reconstructed from the corresponding perturbation in either of the two Weyl scalars $\psi_0$ or $\psi_4$, using a procedure described by Chrzanowski and others in the 1970s. More recent work, motivated within the context of self-force physics, extends the procedure to metric perturbations sourced by a particle in a bound geodesic orbit. However, the existing procedure leaves undetermined a certain stationary, axially-symmetric piece of the metric perturbation. In the vacuum region away from the particle, this "completion" piece corresponds simply to mass and angular-momentum perturbations of the Kerr background, with amplitudes that are, however, a priori unknown. Here we present and implement a rigorous method for finding the completion piece. The key idea is to impose continuity, off the particle, of certain gauge-invariant fields constructed from the full (completed) perturbation, in order to determine the unknown amplitude parameters of the completion piece. We implement this method in full for bound (eccentric) geodesic orbits in the equatorial plane of the Kerr black hole. Our results provide a rigorous underpinning of recent results by Friedman {\it et al.}\ for circular orbits, and extend them to non-circular orbits.

We present a calculation of the scalar field self-force (SSF) acting on a scalar-charge particle in a strong-field orbit around a Kerr black hole. Our calculation specializes to circular and equatorial geodesic orbits. The analysis is an implementation of the standard mode-sum regularization scheme: We first calculate the multipole modes of the scalar-field perturbation using numerical integration in the frequency domain, and then apply a certain regularization procedure to each of the modes. The dissipative piece of the SSF is found to be consistent with the flux of energy and angular momentum carried by the scalar waves through the event horizon and out to infinity. The conservative (radial) component of the SSF is found to be attractive (inward pointing) for $r_0>r_{\rm c}(a)$ and repulsive (outward pointing) for $r_0<r_{\rm c}(a)$, where $a$ is the Kerr spin parameter, $r_0$ is the Boyer-Lindquist orbital radius, and $r_{\rm c}$ is a critical $a$-dependent radius at which the conservative SSF vanishes. When the motion is retrograde the conservative SSF is repulsive for all $r_0$ (as in the Schwarzschild case). The dominant conservative effect of the SSF in Schwarzschild spacetime is known to be of 3rd post-Newtonian (PN) order (with a logarithmic running). Our numerical results suggest that the leading-order PN correction due to the black hole's spin arises from spin-orbit coupling at 3PN, which dominates the overall SSF effect at large $r_0$. In PN language, the change-of-sign of the radial SSF is attributed to an interplay between the spin-orbit term ($\propto -ar_0^{-4.5}$) and the "Schwarzschild" term ($\propto r_0^{-5}\log r_0$).

We present a numerical code for calculating the local gravitational self-force acting on a pointlike particle in a generic (bound) geodesic orbit around a Schwarzschild black hole. The calculation is carried out in the Lorenz gauge: For a given geodesic orbit, we decompose the Lorenz-gauge metric perturbation equations (sourced by the delta-function particle) into tensorial harmonics, and solve for each harmonic using numerical evolution in the time domain (in 1+1 dimensions). The physical self-force along the orbit is then obtained via mode-sum regularization. The total self-force contains a dissipative piece as well as a conservative piece, and we describe a simple method for disentangling these two pieces in a time-domain framework. The dissipative component is responsible for the loss of orbital energy and angular momentum through gravitational radiation; as a test of our code we demonstrate that the work done by the dissipative component of the computed force is precisely balanced by the asymptotic fluxes of energy and angular momentum, which we extract independently from the wave-zone numerical solutions. The conservative piece of the self force does not affect the time-averaged rate of energy and angular-momentum loss, but it influences the evolution of the orbital phases; this piece is calculated here for the first time in eccentric strong-field orbits. As a first concrete application of our code we recently reported the value of the shift in the location and frequency of the innermost stable circular orbit due to the conservative self-force [Phys. Rev. Lett.\ {\bf 102}, 191101 (2009)]. Here we provide full details of this analysis, and discuss future applications.

We present a numerical code for calculating the self force on a scalar charge moving in a bound (eccentric) geodesic in the equatorial plane of a Kerr black hole. We work in the frequency domain and make use of the method of extended homogeneous solutions [Phys.\ Rev.\ D {\bf 78}, 084021 (2008)], in conjunction with mode-sum regularization. Our work is part of a program to develop a computational architecture for fast and efficient self-force calculations, alternative to time-domain methods. We find that our frequency-domain method outperforms existing time-domain schemes for small eccentricities, and, remarkably, remains competitive up to eccentricities as high as $\sim 0.7$. As an application of our code we (i) compute the conservative scalar-field self-force correction to the innermost stable circular equatorial orbit, as a function of the Kerr spin parameter; and (ii) calculate the variation in the rest mass of the scalar particle along the orbit, caused by the component of the self force tangent to the four-velocity.

For a self-gravitating particle of mass \mu in orbit around a Kerr black hole of mass M >> \mu, we compute the O(\mu/M) shift in the frequency of the innermost stable circular equatorial orbit (ISCEO) due to the conservative piece of the gravitational self-force acting on the particle. Our treatment is based on a Hamiltonian formulation of the dynamics in terms of geodesic motion in a certain locally-defined effective smooth spacetime. We recover the same result using the so-called first law of binary black-hole mechanics. We give numerical results for the ISCEO frequency shift as a function of the black hole's spin amplitude, and compare with predictions based on the post-Newtonian approximation and the effective one-body model. Our results provide an accurate strong-field benchmark for spin effects in the general relativistic two-body problem.

Calculations of the gravitational self-force (GSF) in curved spacetime require as input the metric perturbation in a sufficiently regular gauge. A basic challenge in the program to compute the GSF for orbits around a Kerr black hole is that the standard procedure for reconstructing the perturbation is formulated in a class of radiation gauges, in which the particle singularity is non-isotropic and extends away from the particle's location. Here we present two practical schemes for calculating the GSF using a radiation-gauge reconstructed metric as input. The schemes are based on a detailed analysis of the local structure of the particle singularity in the radiation gauges. We identify 3 types of radiation gauges: two containing a radial string-like singularity emanating from the particle, either in one direction ("half-string" gauges) or both directions ("full-string" gauges); and a third type containing no strings but with a jump discontinuity across a surface intersecting the particle. Based on a flat-space example, we argue that the standard mode-by-mode reconstruction procedure yields the "regular half" of a half-string solution, or (equivalently) either of the regular halves of a no-string solution. For the half-string case, we formulate the GSF in a locally deformed radiation gauge that removes the string singularity near the particle. We derive a mode-sum formula for the GSF in this gauge, analogous to the standard Lorenz-gauge formula but with modified regularization parameters. For the no-string case, we formulate the GSF directly, without a local deformation, and we derive a mode-sum formula that requires no correction to the parameters but involves a certain averaging procedure. We explain the consistency of our results with Gralla's invariance theorem, and discuss the correspondence between our method and a related approach by Friedman et al.

We present results from calculations of the orbital evolution in eccentric binaries of nonrotating black holes with extreme mass-ratios. Our inspiral model is based on the method of osculating geodesics, and is the first to incorporate the full gravitational self-force (GSF) effect, including conservative corrections. The GSF information is encapsulated in an analytic interpolation formula based on numerical GSF data for over a thousand sample geodesic orbits. We assess the importance of including conservative GSF corrections in waveform models for gravitational-wave searches.

Bound geodesic orbits around a Kerr black hole can be parametrized by three constants of the motion: the (specific) orbital energy, angular momentum and Carter constant. Generically, each orbit also has associated with it three frequencies, related to the radial, longitudinal and (mean) azimuthal motions. Here we note the curious fact that these two ways of characterizing bound geodesics are not in a one-to-one correspondence. While the former uniquely specifies an orbit up to initial conditions, the latter does not: there is a (strong-field) region of the parameter space in which pairs of physically distinct orbits can have the same three frequencies. In each such isofrequency pair the two orbits exhibit the same rate of periastron precession and the same rate of Lense-Thirring precession of the orbital plane, and (in a certain sense) they remain "synchronized" in phase.

We compute the conservative piece of the gravitational self-force (GSF) acting on a particle of mass m_1 as it moves along an (unstable) circular geodesic orbit between the innermost stable circular orbit (ISCO) and the light ring of a Schwarzschild black hole of mass m_2>> m_1. More precisely, we construct the function h_{uu}(x) = h_{\mu\nu} u^{\mu} u^{\nu} (related to Detweiler's gauge-invariant "redshift" variable), where h_{\mu\nu} is the regularized metric perturbation in the Lorenz gauge, u^{\mu} is the four-velocity of m_1, and x= [Gc^{-3}(m_1+m_2)\Omega]^{2/3} is an invariant coordinate constructed from the orbital frequency \Omega. In particular, we explore the behavior of h_{uu} just outside the "light ring" at x=1/3, where the circular orbit becomes null. Using the recently discovered link between h_{uu} and the piece a(u), linear in the symmetric mass ratio \nu, of the main radial potential A(u,\nu) of the Effective One Body (EOB) formalism, we compute a(u) over the entire domain 0<u<1/3. We find that a(u) diverges at the light-ring as ~0.25 (1-3u)^{-1/2}, explain the physical origin of this divergence, and discuss its consequences for the EOB formalism. We construct accurate global analytic fits for a(u), valid on the entire domain 0<u<1/3 (and possibly beyond), and give accurate numerical estimates of the values of a(u) and its first 3 derivatives at the ISCO, as well as the O(\nu) shift in the ISCO frequency. In previous work we used GSF data on slightly eccentric orbits to compute a certain linear combination of a(u) and its first two derivatives, involving also the O(\nu) piece \bar d(u) of a second EOB radial potential {\bar D}(u,\nu). Combining these results with our present global analytic representation of a(u), we numerically compute {\bar d}(u)$ on the interval 0<u\leq 1/6.

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.

This is the second in a series of papers aimed at developing a practical time-domain method for self-force calculations in Kerr spacetime. The key elements of the method are (i) removal of a singular part of the perturbation field with a suitable analytic "puncture" based on the Detweiler--Whiting decomposition, (ii) decomposition of the perturbation equations in azimuthal ($m$-)modes, taking advantage of the axial symmetry of the Kerr background, (iii) numerical evolution of the individual $m$-modes in 2+1-dimensions with a finite difference scheme, and (iv) reconstruction of the physical self-force from the mode sum. Here we report an implementation of the method to compute the scalar-field self-force along circular equatorial geodesic orbits around a Kerr black hole. This constitutes a first time-domain computation of the self force in Kerr geometry. Our time-domain code reproduces the results of a recent frequency-domain calculation by Warburton and Barack, but has the added advantage of being readily adaptable to include the back-reaction from the self force in a self-consistent manner. In a forthcoming paper---the third in the series---we apply our method to the gravitational self-force (in the Lorenz gauge).

The general relativistic (Mercury-type) periastron advance is calculated here for the first time with exquisite precision in full general relativity. We use accurate numerical relativity simulations of spinless black hole binaries with mass ratios 1/8 < m1/m2 < 1 and compare with the predictions of several analytic approximation schemes. We find the effective-one-body model to be remarkably accurate, and, surprisingly, so also the predictions of self-force theory [replacing m1/m2 --> m1m2/(m1+m2)^2]. Our results can inform a universal analytic model of the two-body dynamics, crucial for ongoing and future gravitational-wave searches.

We study conservative finite-mass corrections to the motion of a particle in a bound (eccentric) strong-field orbit around a Schwarzschild black hole. We assume the particle's mass $\mu$ is much smaller than the black hole mass $M$, and explore post-geodesic corrections of $O(\mu/M)$. Our analysis uses numerical data from a recently developed code that outputs the Lorenz-gauge gravitational self-force (GSF) acting on the particle along the eccentric geodesic. First, we calculate the $O(\mu/M)$ conservative correction to the periastron advance of the orbit, as a function of the (gauge-dependent) semilatus rectum and eccentricity. A gauge-invariant description of the GSF precession effect is made possible in the circular-orbit limit, where we express the correction to the periastron advance as a function of the invariant azimuthal frequency. We compare this relation with results from fully nonlinear numerical-relativistic simulations. In order to obtain a gauge-invariant measure of the GSF effect for fully eccentric orbits, we introduce a suitable generalization of Detweiler's circular-orbit "redshift" invariant. We compute the $O(\mu/M)$ conservative correction to this invariant, expressed as a function of the two invariant frequencies that parametrize the orbit. Our results are in good agreement with results from post-Newtonian calculations in the weak-field regime, as we shall report elsewhere. The results of our study can inform the development of analytical models for the dynamics of strongly gravitating binaries. They also provide an accurate benchmark for future numerical-relativistic simulations.

To model the radiative evolution of extreme mass-ratio binary inspirals (a key target of the LISA mission), the community needs efficient methods for computation of the gravitational self-force (SF) on the Kerr spacetime. Here we further develop a practical `$m$-mode regularization' scheme for SF calculations, and give details of a first implementation. The key steps in the method are (i) removal of a singular part of the perturbation field with a suitable `puncture' to leave a sufficiently regular residual within a finite worldtube surrounding the particle's worldline, (ii) decomposition in azimuthal ($m$-)modes, (iii) numerical evolution of the $m$-modes in 2+1D with a finite difference scheme, and (iv) reconstruction of the SF from the mode sum. The method relies on a judicious choice of puncture, based on the Detweiler--Whiting decomposition. We give a working definition for the `order' of the puncture, and show how it determines the convergence rate of the $m$-mode sum. The dissipative piece of the SF displays an exponentially convergent mode sum, while the $m$-mode sum for the conservative piece converges with a power law. In the latter case the individual modal contributions fall off at large $m$ as $m^{-n}$ for even $n$ and as $m^{-n+1}$ for odd $n$, where $n$ is the puncture order. We describe an $m$-mode implementation with a 4th-order puncture to compute the scalar-field SF along circular geodesics on Schwarzschild. In a forthcoming companion paper we extend the calculation to the Kerr spacetime.

We calculate the gravitational self force acting on a pointlike particle of mass $\mu$, set in a circular geodesic orbit around a Schwarzschild black hole. Our calculation is done in the Lorenz gauge: For given orbital radius, we first solve directly for the Lorenz-gauge metric perturbation using numerical evolution in the time domain; We then compute the (finite) back-reaction force from each of the multipole modes of the perturbation; Finally, we apply the ``mode sum'' method to obtain the total, physical self force. The {\em temporal} component of the self force (which is gauge invariant) describes the dissipation of orbital energy through gravitational radiation. Our results for this component are consistent, to within the computational accuracy, with the total flux of gravitational-wave energy radiated to infinity and through the event horizon. The {\em radial} component of the self force (which is gauge dependent) is calculated here for the first time. It describes a conservative shift in the orbital parameters away from their geodesic values. We thus obtain the $O(\mu)$ correction to the specific energy and angular momentum parameters (in the Lorenz gauge), as well as the $O(\mu)$ shift in the orbital frequency (which is gauge invariant).

This review is concerned with the gravitational self-force acting on a mass particle in orbit around a large black hole. Renewed interest in this old problem is driven by the prospects of detecting gravitational waves from strongly gravitating binaries with extreme mass ratios. We begin here with a summary of recent advances in the theory of gravitational self-interaction in curved spacetime, and proceed to survey some of the ideas and computational strategies devised for implementing this theory in the case of a particle orbiting a Kerr black hole. We review in detail two of these methods: (i) the standard mode-sum method, in which the metric perturbation is regularized mode-by-mode in a multipole decomposition, and (ii) $m$-mode regularization, whereby individual azimuthal modes of the metric perturbation are regularized in 2+1 dimensions. We discuss several practical issues that arise, including the choice of gauge, the numerical representation of the particle singularity, and how high-frequency contributions near the particle are dealt with in frequency-domain calculations. As an example of a full end-to-end implementation of the mode-sum method, we discuss the computation of the gravitational self-force for eccentric geodesic orbits in Schwarzschild, using a direct integration of the Lorenz-gauge perturbation equations in the time domain. With the computational framework now in place, researchers have recently turned to explore the physical consequences of the gravitational self force; we will describe some preliminary results in this area. An appendix to this review presents, for the first time, a detailed derivation of the regularization parameters necessary for implementing the mode-sum method in Kerr spacetime.

The innermost stable circular orbit (ISCO) of a test particle around a Schwarzschild black hole of mass $M$ has (areal) radius $r_{\rm isco}= 6M G/c^2$. If the particle is endowed with mass $\mu(\ll M)$, it experiences a gravitational self-force whose conservative piece alters the location of the ISCO. Here we calculate the resulting shifts $\Delta r_{\rm isco}$ and $\Delta\Omega_{\rm isco}$ in the ISCO's radius and frequency, at leading order in the mass ratio $\mu/M$. We obtain, in the Lorenz gauge, $\Delta r_{\rm isco}=-3.269 (\pm 0.003)\mu G/c^2$ and $\Delta\Omega_{\rm isco}/\Omega_{\rm isco}=0.4870 (\pm 0.0006) \mu/M$. We discuss the implications of our result within the context of the extreme-mass-ratio binary inspiral problem.

We present a new, simple method for calculating the scalar, electromagnetic, and gravitational self forces acting on particles in orbit around a Kerr black hole. The standard ``mode-sum regularization'' approach for self-force calculations relies on a decomposition of the full (retarded) perturbation field into multipole modes, followed by the application of a certain mode-by-mode regularization procedure. In recent years several groups have developed numerical codes for calculating black hole perturbations directly in 2+1 dimensions (i.e., decomposing the azimuthal dependence into $m$-modes, but refraining from a full multipole decomposition). Here we formulate a practical scheme for constructing the self force directly from the 2+1-dimensional $m$-modes. While the standard mode-sum method is serving well in calculations of the self force in Schwarzschild geometry, the new scheme should allow a more efficient treatment of the Kerr problem.

Inspirals of stellar-mass compact objects into $\sim 10^6 M_{\odot}$ black holes are especially interesting sources of gravitational waves for LISA. We investigate whether the emitted waveforms can be used to strongly constrain the geometry of the central massive object, and in essence check that it corresponds to a Kerr black hole (BH). For a Kerr BH, all multipole moments of the spacetime have a simple, unique relation to $M$ and $S$, the BH's mass and spin; in particular, the spacetime's mass quadrupole moment is given by $Q=- S^2/M$. Here we treat $Q$ as an additional parameter, independent of $M$ and $S$, and ask how well observation can constrain its difference from the Kerr value. This was already estimated by Ryan, but for simplified (circular, equatorial) orbits, and neglecting signal modulations due to the motion of the LISA satellites. Here we consider generic orbits and include these modulations. We use a family of approximate (post-Newtonian) waveforms, which represent the full parameter space of Inspiral sources, and exhibit the main qualitative features of true, general relativistic waveforms. We extend this parameter space to include (in an approximate manner) an arbitrary value of $Q$, and construct the Fisher information matrix for the extended parameter space. By inverting the Fisher matrix we estimate how accurately $Q$ could be extracted from LISA observations. For 1 year of coherent data from the inspiral of a $10 M_{\odot}$ BH into rotating BHs of masses $10^{5.5} M_{\odot}$, $10^6 M_{\odot}$, or $10^{6.5} M_{\odot}$, we find $\Delta (Q/M^3) \sim 10^{-4}$, $10^{-3}$, or $10^{-2}$, respectively (assuming total signal-to-noise ratio of 100, typical of the brightest detectable EMRIs). These results depend only weakly on the eccentricity of the orbit or the BH's spin.

Captures of compact objects (COs) by massive black holes (MBHs) in galactic nuclei will be an important source for LISA, the space-based gravitational-wave (GW) detector. However, a large fraction of captures will not be individually resolvable--either because they are too distant, have unfavorable orientation, or have too many years to go before final plunge--and so will constitute a source of ``confusion noise,'' obscuring other types of sources. Here we estimate the shape and overall magnitude of the spectrum of confusion noise from CO captures. The overall magnitude depends on the capture rates, which are rather uncertain, so we present results for a plausible range of rates. We show that the impact of capture confusion noise on the total LISA noise curve ranges from insignificant to modest, depending on these rates. Capture rates at the high end of estimated ranges would raise LISA's overall (effective) noise level by at most a factor \sim 2. While this slightly elevated noise level would somewhat decrease LISA's sensitivity to other classes of sources, overall, this would be a pleasant problem for LISA to have: It would also imply that detection rates for CO captures were at nearly their maximum possible levels (given LISA's baseline design). This paper includes several other results that should be useful in further studies of LISA capture sources, including (i) a calculation of the total GW energy output from generic inspirals into Kerr MBHs, and (ii) an approximate GW energy spectrum for a typical capture.

Captures of stellar-mass compact objects (COs) by massive ($\sim 10^6 M_\odot$) black holes (MBHs) are potentially an important source for LISA, the proposed space-based gravitational-wave (GW) detector. The orbits of the inspiraling COs are highly complicated; they can remain rather eccentric up until the final plunge, and display extreme versions of relativistic perihelion precession and Lense-Thirring precession of the orbital plane. The strongest capture signals will be ~10 times weaker than LISA's instrumental noise, but in principle (with sufficient computing power) they can be disentangled from the noise by matched filtering. The associated template waveforms are not yet in hand, but theorists will very likely be able to provide them before LISA launches. Here we introduce a family of approximate (post-Newtonian) capture waveforms, given in (nearly) analytic form, for use in advancing LISA studies until more accurate versions are available. Our model waveforms include most of the key qualitative features of true waveforms, and cover the full space of capture-event parameters (including orbital eccentricity and the MBH's spin). Here we use our approximate waveforms to (i) estimate the relative contributions of different harmonics (of the orbital frequency) to the total signal-to-noise ratio, and (ii) estimate the accuracy with which LISA will be able to extract the physical parameters of the capture event from the measured waveform. For a typical source (a $10 M_\odot$ CO captured by a $10^6 M_\odot$ MBH at a signal-to-noise ratio of 30), we find that LISA can determine the MBH and CO masses to within a fractional error of $\sim 10^{-4}$, measure $S/M^2$ (where $S$ and $M$ are the MBH's mass and spin) to within $\sim 10^{-4}$, and determine the sky location of the source to within $\sim 10^{-3}$ stradians.

One of the most exciting prospects for the LISA gravitational wave observatory is the detection of gravitational radiation from the inspiral of a compact object into a supermassive black hole. The large inspiral parameter space and low amplitude of the signal makes detection of these sources computationally challenging. We outline here a first cut data analysis scheme that assumes realistic computational resources. In the context of this scheme, we estimate the signal-to-noise ratio that a source requires to pass our thresholds and be detected. Combining this with an estimate of the population of sources in the Universe, we estimate the number of inspiral events that LISA could detect. The preliminary results are very encouraging -- with the baseline design, LISA can see inspirals out to a redshift z=1 and should detect over a thousand events during the mission lifetime.