by Warburton, Niels and Akcay, Sarp and Barack, Leor and Gair, Jonathan R. and Sago, Norichika
4.3 pages, 3 figures

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.

by Tiec, Alexandre Le and Barausse, Enrico and Buonanno, Alessandra
5 pages, 1 figure

Using the first law of binary black-hole mechanics, we compute the binding energy E and total angular momentum J of two non-spinning compact objects moving on circular orbits with frequency Omega, at leading order beyond the test-particle approximation. By minimizing E(Omega) we recover the exact frequency shift of the Schwarzschild innermost stable circular orbit induced by the conservative piece of the gravitational self-force. Comparing our results for the coordinate invariant relation E(J) to those recently obtained from numerical simulations of comparable-mass non-spinning black-hole binaries, we find a remarkably good agreement, even in the strong-field regime. Our findings confirm that the domain of validity of perturbative calculations may extend well beyond the extreme mass-ratio limit.

by Dolan, Sam R. and Wardell, Barry and Barack, Leor
30 pages, 5 figures, 3 tables

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 ($latex m$-)modes, taking advantage of the axial symmetry of the Kerr background, (iii) numerical evolution of the individual $latex 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).

by Gralla, Samuel E.

A difficulty with previous treatments of the gravitational self-force is that an explicit formula for the force is available only in a particular gauge (Lorenz gauge), where the force in other gauges must be found through a transformation law. For a class of gauges satisfying a “parity condition” ensuring that the Hamiltonian center of mass of the particle is well-defined, I show that the gravitational self-force is always given by the angle-average of the bare gravitational force. To derive this result I replace the computational strategy of previous work with a new approach, wherein the form of the force is first fixed up to a gauge-invariant piece by simple manipulations, and then that piece is determined by working in a gauge designed specifically to simplify the computation. This offers significant computational savings over the Lorenz gauge, since the Hadamard expansion is avoided entirely and the metric perturbation takes a very simple form. I also show that the rest mass of the particle does not evolve due to first-order self-force effects. Finally, I consider the “mode sum regularization” scheme for computing the self-force in black hole background spacetimes, and use the angle-average form of the force to show that the same mode-by-mode subtraction may be performed in all gauges satisfying the parity condition. This helps provide a practical foundation for the computation of self-forces in the Kerr background.

by Aoudia, Sofiane and Spallicci, Alessandro D. A. M.
12 pages, 13 figures (additional figures and text revised in v2 arXiv)

Perturbations of Schwarzschild-Droste black holes in the Regge-Wheeler gauge benefit from the availability of a wave equation and from the gauge invariance of the wave function, but lack smoothness. Nevertheless, the even perturbations belong to the C\textsuperscript{0} continuity class, if the wave function and its derivatives satisfy specific conditions on the discontinuities, known as jump conditions, at the particle position. These conditions suggest a new way for dealing with finite element integration in time domain. The forward time value in the upper node of the $latex (t, r^*$) grid cell is obtained by the linear combination of the three preceding node values and of analytic expressions based on the jump conditions. The numerical integration does not deal directly with the source term, the associated singularities and the potential. This amounts to an indirect integration of the wave equation. The known wave forms at infinity are recovered and the wave function at the particle position is shown. In this series of papers, the radial trajectory is dealt with first, being this method of integration applicable to generic orbits of EMRI (Extreme Mass Ratio Inspiral).

by Barack, Leor and Sago, Norichika
29 pages, 4 eps figures

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 $latex \mu$ is much smaller than the black hole mass $latex M$, and explore post-geodesic corrections of $latex 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 $latex O(\mu/M)$ conservative correction to the periastron advance of the orbit, as a function of the (gauge dependent) semi-latus 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 “red shift” invariant. We compute the $latex 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.

by Canizares, Priscilla and Sopuerta, Carlos F.
6 pages, 7 figures, submitted to proceedings of the 8th International LISA Symposium, Stanford, June 28 – July 2, 2010

The gravitational-wave signals emitted by Extreme-Mass-Ratio Inspirals will be hidden in the instrumental LISA noise and the foreground noise produced by galactic binaries in the LISA band. Then, we need accurate gravitational-wave templates to extract these signals from the noise and obtain the relevant physical parameters. This means that in the modeling of these systems we have to take into account how the orbit of the stellar-mass compact object is modified by the action of its own gravitational field. This effect can be described as the action of a local force, the self-force. We present a time-domain technique to compute the self-force for geodesic eccentric orbits around a non-rotating massive black hole. To illustrate the method we have applied it to a testbed model consisting of scalar charged particle orbiting a non-dynamical black hole. A key feature of our method is that it does not introduce a small scale associated with the stellar-mass compact object. This is achieved by using a multidomain framework where the particle is located at the interface between two subdomains. In this way, we just have to evolve homogeneous wave-like equations with smooth solutions that have to be communicated across the subdomain boundaries using appropriate junction conditions. The numerical technique that we use to implement this scheme is the pseudospectral collocation method. We show the suitability of this technique for the modeling of Extreme-Mass-Ratio Inspirals and show that it can provide accurate results for the self-force.

by Huerta, E. A. and Gair, Jonathan R
6 pages, 1 figure, submitted to proceedings of the 8th International LISA Symposium, Stanford, June 28 – July 2, 2010

It is not currently clear how important it will be to include conservative self-force (SF) corrections in the models for extreme-mass-ratio inspiral (EMRI) waveforms that will be used to detect such signals in LISA (Laser Interferometer Space Antenna) data. These proceedings will address this issue for circular-equatorial inspirals using an approximate EMRI model that includes conservative corrections at leading post-Newtonian order. We will present estimates of the magnitude of the parameter estimation errors that would result from omitting conservative corrections, and compare these to the errors that will arise from noise fluctuations in the detector. We will also use this model to explore the relative importance of the second-order radiative piece of the SF, which is not presently known.

by Flanagan, Eanna E. and Hinderer, Tanja
5 pages, 1 figure

We show that transient resonances occur in the two body problem in general relativity, in the highly relativistic, extreme mass-ratio regime for spinning black holes. These resonances occur when the ratio of polar and radial orbital frequencies, which is slowly evolving under the influence of gravitational radiation reaction, passes through a low order rational number. At such points, the adiabatic approximation to the orbital evolution breaks down, and there is a brief but order unity correction to the inspiral rate. Corrections to the gravitational wave signal’s phase due to resonance effects scale as the square root of the inverse of mass of the small body, and thus become large in the extreme-mass-ratio limit, dominating over all other post-adiabatic effects. The resonances make orbits more sensitive to changes in initial data (though not quite chaotic), and are genuine non-perturbative effects that are not seen at any order in a standard post-Newtonian expansion. Our results apply to an important potential source of gravitational waves, the gradual inspiral of white dwarfs, neutron stars, or black holes into much more massive black holes. It is hoped to exploit observations of these sources to map the spacetime geometry of black holes. However, such mapping will require accurate models of binary dynamics, which is a computational challenge whose difficulty is significantly increased by resonance effects. We estimate that the resonance phase shifts will be of order a few tens of cycles for mass ratios $latex \sim 10^{-6}$, by numerically evolving fully relativistic orbital dynamics supplemented with an approximate, post-Newtonian self-force.

by Shah, Abhay and Keidl, Tobias and Friedman, John and Kim, Dong-Hoon and Price, Larry
21 pages, 2 figures

This is the second of two companion papers on computing the self-force in a radiation gauge; more precisely, the method uses a radiation gauge for the radiative part of the metric perturbation, together with an arbitrarily chosen gauge for the parts of the perturbation associated with changes in black-hole mass and spin and with a shift in the center of mass. We compute the conservative part of the self-force for a particle in circular orbit around a Schwarzschild black hole. The gauge vector relating our radiation gauge to a Lorenz gauge is helically symmetric, implying that the quantity h_{\alpha\beta} u^\alpha u^\beta (= h_{uu}) must have the same value for our radiation gauge as for a Lorenz gauge; and we confirm this numerically to one part in 10^{13}. As outlined in the first paper, the perturbed metric is constructed from a Hertz potential that is in term obtained algebraically from the the retarded perturbed spin-2 Weyl scalar, \psi_0 . We use a mode-sum renormalization and find the renormalization coefficients by matching a series in L = \ell + 1/2 to the large-L behavior of the expression for the self-force in terms of the retarded field h_{\alpha\beta}^{ret}; we similarly find the leading renormalization coefficients of h_{uu} and the related change in the angular velocity of the particle due to its self-force. We show numerically that the singular part of the self-force has the form f_{\alpha} \propto , the part of \nabla_\alpha \rho^{-1} that is axisymmetric about a radial line through the particle. This differs only by a constant from its form for a Lorenz gauge. It is because we do not use a radiation gauge to describe the change in black-hole mass that the singular part of the self-force has no singularity along a radial line through the particle and, at least in this example, is spherically symmetric to subleading order in \rho.

by Blanchet, Luc and Detweiler, Steven and Tiec, Alexandre Le and Whiting, Bernard F.
29 pages, 3 figures; to appear in the book “Mass and Motion in General Relativity”, proceedings of the C.N.R.S. School in Orleans, France, eds. L. Blanchet, A. Spallicci and B. F. Whiting

The relativistic motion of a compact binary system moving in circular orbit is investigated using the post-Newtonian (PN) approximation and the perturbative self-force (SF) formalism. A particular gauge-invariant observable quantity is computed as a function of the binary’s orbital frequency. The conservative effect induced by the gravitational SF is obtained numerically with high precision, and compared to the PN prediction developed to high order. The PN calculation involves the computation of the 3PN regularized metric at the location of the particle. Its divergent self-field is regularized by means of dimensional regularization. The poles proportional to 1/(d-3) which occur within dimensional regularization at the 3PN order disappear from the final gauge-invariant result. The leading 4PN and next-to-leading 5PN conservative logarithmic contributions originating from gravitational-wave tails are also obtained. Making use of these exact PN results, some previously unknown PN coefficients are measured up to the very high 7PN order by fitting to the numerical self-force data. Using just the 2PN and new logarithmic terms, the value of the 3PN coefficient is also confirmed numerically with very high precision. The consistency of this cross-cultural comparison provides a crucial test of the very different regularization methods used in both SF and PN formalisms, and illustrates the complementarity of these approximation schemes when modelling compact binary systems.

by Thornburg, Jonathan
27 pages, 12 eps figures (10 of them color, but all are viewable ok in black-and-white), uses RevTeX 4.1

If a small “particle” of mass $latex \mu M$ (with $latex \mu \ll 1$) orbits a Schwarzschild or Kerr black hole of mass $latex M$, the particle is subject to an $latex \O(\mu)$ radiation-reaction “self-force”. Here I argue that it’s valuable to compute this self-force highly accurately (relative error of $latex \ltsim 10^{-6}$) and efficiently, and I describe techniques for doing this and for obtaining and validating error estimates for the computation. I use an adaptive-mesh-refinement (AMR) time-domain numerical integration of the perturbation equations in the Barack-Ori mode-sum regularization formalism; this is efficient, yet allows easy generalization to arbitrary particle orbits. I focus on the model problem of a scalar particle in a circular geodesic orbit in Schwarzschild spacetime.

The mode-sum formalism gives the self-force as an infinite sum of regularized spherical-harmonic modes $latex \sum_{\ell=0}^\infty F_{\ell,\reg}$, with $latex F_{\ell,\reg}$ (and an “internal” error estimate) computed numerically for $latex \ell \ltsim 30$ and estimated for larger~$latex \ell$ by fitting an asymptotic “tail” series. Here I validate the internal error estimates for the individual $latex F_{\ell,\reg}$ using a large set of numerical self-force computations of widely-varying accuracies. I present numerical evidence that the actual numerical errors in $latex F_{\ell,\reg}$ for different~$latex \ell$ are at most weakly correlated, so the usual statistical error estimates are valid for computing the self-force. I show that the tail fit is numerically ill-conditioned, but this can be mostly alleviated by renormalizing the basis functions to have similar magnitudes.

Using AMR, fixed mesh refinement, and extended-precision floating-point arithmetic, I obtain the (contravariant) radial component of the self-force for a particle in a circular geodesic orbit of areal radius $latex r = 10M$ to within $latex 1$~ppm relative error.

by Barack, Leor and Sago, Norichika
42 pages

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.

by Blanchet, Luc and Detweiler, Steven and Tiec, Alexandre Le and Whiting, Bernard F.
32 pages, 2 figures

We continue a previous work on the comparison between the post-Newtonian (PN) approximation and the gravitational self-force (SF) analysis of circular orbits in a Schwarzschild background. We show that the numerical SF data contain physical information corresponding to extremely high PN approximations. We find that knowing analytically determined appropriate PN parameters helps tremendously in allowing the numerical data to be used to obtain higher order PN coefficients. Using standard PN theory we compute analytically the leading 4PN and the next-to-leading 5PN logarithmic terms in the conservative part of the dynamics of a compact binary system. The numerical perturbative SF results support well the analytic PN calculations through first order in the mass ratio, and are used to accurately measure the 4PN and 5PN non-logarithmic coefficients in a particular gauge invariant observable. Furthermore we are able to give estimates of higher order contributions up to the 7PN level. In our best fit we also confirm with high precision the value of the 3PN coefficient. This interplay between PN and SF efforts is important for the synthesis of template waveforms of extreme mass ratio inspirals to be analysed by the space-based gravitational wave instrument LISA. Our work will also have an impact on efforts that combine numerical results in a quantitative analytical framework so as to generate complete inspiral waveforms for the ground-based detection of gravitational waves by instruments such as LIGO and Virgo.

by Canizares, Priscilla and Sopuerta, Carlos F.
3 pages. To appear in Proceedings of the Twelfth Marcel Grossmann Meeting on General Relativity, edited by Thibault Damour, Robert T Jantzen and Remo Ruffini, World Scientific, Singapore, 2010

We introduce a new time-domain method for computing the self-force acting on a scalar particle in a Schwarzschild geometry. The principal feature of our method consists in the division of the spatial domain into several subdomains and locating the particle at the interface betweem two them. In this way, we avoid the need of resolving a small length scale associated with the presence of a particle in the computational domain and, at the same time, we avoid numerical problems due to the low differentiability of solutions of equations with point-like singular behaviour.

by Field, Scott E. and Hesthaven, Jan S. and Lau, Stephen R.
Uses revtex4, 23 pages, 9 figures, 3 tables

In the context of metric perturbation theory for non-spinning black holes, extreme mass ratio binary (EMRB) systems are described by distributionally forced master wave equations. Numerical solution of a master wave equation as an initial boundary value problem requires initial data. However, because the correct initial data for generic-orbit systems is unknown, specification of trivial initial data is a common choice, despite being inconsistent and resulting in a solution which is initially discontinuous in time. As is well known, this choice leads to a “burst” of junk radiation which eventually propagates off the computational domain. We observe another unintended consequence of trivial initial data: development of a persistent spurious solution, here referred to as the Jost junk solution, which contaminates the physical solution for long times. This work studies the influence of both types of junk on metric perturbations, waveforms, and self-force measurements, and it demonstrates that smooth modified source terms mollify the Jost solution and reduce junk radiation.

by Blanchet, Luc and Detweiler, Steven and Tiec, Alexandre Le and Whiting, Bernard F.
36 pages, 3 figures

The problem of a compact binary system whose components move on circular orbits is addressed using two different approximation techniques in general relativity. The post-Newtonian (PN) approximation involves an expansion in powers of v/c<<1, and is most appropriate for small orbital velocities v. The perturbative self-force (SF) analysis requires an extreme mass ratio m1/m2<<1 for the components of the binary. A particular coordinate-invariant observable is determined as a function of the orbital frequency of the system using these two different approximations. The post-Newtonian calculation is pushed up to the third post-Newtonian (3PN) order. It involves the metric generated by two point particles and evaluated at the location of one of the particles. We regularize the divergent self-field of the particle by means of dimensional regularization. We show that the poles proportional to 1/(d-3) appearing in dimensional regularization at the 3PN order cancel out from the final gauge invariant observable. The 3PN analytical result, through first order in the mass ratio, and the numerical SF calculation are found to agree well. The consistency of this cross cultural comparison confirms the soundness of both approximations in describing compact binary systems. In particular, it provides an independent test of the very different regularization procedures invoked in the two approximation schemes.

by Spallicci, A. and Aoudia, S.
To be published on 21 Rencontres de Blois: Windows on the Universe, http://confs.obspm.fr/Blois2009/, 4 pages 1 figure

Through detection by low gravitational wave space interferometers, the capture of stars by supermassive black holes will constitute a giant step forward in the understanding of gravitation in strong field. The impact of the perturbations on the motion of the star is computed via the tail, the back-scattered part of the perturbations, or via a radiative Green function. In the former approach, the self-force acts upon the background geodesic, while in the latter, the geodesic is conceived in the total (background plus perturbations) field. Regularisations (mode-sum and Riemann-Hurwitz $latex \zeta$ function) intervene to cancel divergencies coming from the infinitesimal size of the particle. The non-adiabatic trajectories require the most sophisticated techniques for studying the evolution of the motion, like the self-consistent approach.

by Poisson, Eric
18 pages, 5 figures, lecture given at the CNRS School on Mass (Orleans, June 2008)

I present an overview of the methods involved in the computation of the scalar, electromagnetic, and gravitational self-forces acting on a point particle moving in a curved spacetime. For simplicity, the focus here will be on the scalar self-force. The lecture follows closely my review article on this subject published in Living Reviews in Relativity. I begin with a review of geometrical elements (Synge’s world function, the parallel propagator). Next I introduce useful coordinate systems (Fermi normal coordinates and retarded light-cone coordinates) in a neighborhood of the particle’s world line. I then present the wave equation for a scalar field in curved spacetime and the equations of motion for a particle endowed with a scalar charge. The wave equation is solved by means of a Green’s function, and the self-force is constructed from the field gradient. Because the retarded field is singular on the world line, the self-force must involve a regularized version of the field gradient, and I describe how the regular piece of the self-field can be identified. In the penultimate section of the lecture I put the construction of the self-force on a sophisticated axiomatic basis, and in the concluding section I explain how one can do better by abandoning the dangerous fiction of a point particle

by Detweiler, Steven
38 pages, 4 figures, Lecture given at the “School on Mass” (Orleans, France, June 2008), uses Springer’s “svmult.cls”

The gravitational field of a particle of small mass $latex \mu$ moving through curved spacetime, with metric $latex g_{ab}$, is naturally and easily decomposed into two parts each of which satisfies the perturbed Einstein equations through $latex O(\mu)$. One part is an inhomogeneous field $latex h^S_{ab}$ which, near the particle, looks like the Coulomb $latex \mu/r$ field with tidal distortion from the local Riemann tensor. This singular field is defined in a neighborhood of the small particle and does not depend upon boundary conditions or upon the behavior of the source in either the past or the future. The other part is a homogeneous field $latex h^R_{ab}$. In a perturbative analysis, the motion of the particle is then best described as being a geodesic in the metric $latex g_{ab}+h^R_{ab}$. This geodesic motion includes all of the effects which might be called radiation reaction and conservative effects as well.

by Vega, Ian and Diener, Peter and Tichy, Wolfgang and Detweiler, Steven
23 pages, 13 figures

Prescriptions for numerical self-force calculations have traditionally been designed for frequency-domain or (1+1) time-domain codes which employ a mode decomposition to facilitate in carrying out a delicate regularization scheme. This has prevented self-force analyses from benefiting from the powerful suite of tools developed and used by numerical relativists for simulations of the evolution of comparable-mass black hole binaries. In this work, we revisit a previously-introduced (3+1) method for self-force calculations, and demonstrate its viability by applying it to the test case of a scalar charge moving in a circular orbit around a Schwarzschild black hole. Two (3+1) codes originally developed for numerical relativity applications were independently employed, and in each we were able to compute the two independent components of the self-force and the energy flux correctly to within $latex < 1%$. We also demonstrate consistency between $latex t$-component of the self-force and the scalar energy flux. Our results constitute the first successful calculation of a self-force in a (3+1) framework, and thus open opportunities for the numerical relativity community in self-force analyses and the perturbative modeling of extreme-mass-ratio inspirals.

by Barack, Leor
Invited topical review for CQG; 61 pages, 4 eps figures; uses iopart.cls, iopart10.clo

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) $latex 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.

by Galley, Chad R. and Hu, Bei-Lok
17 pages, 3 figures, Invited contribution to the International Conference on Classical and Quantum Relativistic Dynamics of Particles and Fields (IARD) held at the Aristotle University, Thessaloniki, Greece, 22-26 June 2008. Proceedings to appear in Foundations of Physics

We present a new analytical framework for describing the dynamics of a gravitational binary system with unequal masses moving with arbitrary relative velocity, taking into account the backreaction from both compact objects in the form of tidal deformation, gravitational waves and self forces. Allowing all dynamical variables to interact with each other in a self-consistent manner this formalism ensures that all the dynamical quantities involved are conserved on the background spacetime and obey the gauge invariance under general coordinate transformations that preserve the background geometry. Because it is based on a generalized perturbation theory and the important new emphasis is on the self-consistency of all the dynamical variables involved we call it a gravitational perturbation theory with self-consistent backreaction (GP-SCB).

As an illustration of how this formalism is implemented we construct perturbatively a self-consistent set of equations of motion for an inspiraling gravitational binary, which does not require extra assumptions such as slow motion, weak-field or small mass ratio for its formulation. This case should encompass the inspiral and possibly the plunge and merger phases of binaries with otherwise general parameters (e.g., mass ratio and relative velocity) though more investigation is needed to substantiate it.

In the second part, we discuss how the mass ratio can be treated as a perturbation parameter in the post-Newtonian effective field theory (PN-EFT) approach, thus extending the work of Goldberger and Rothstein for equal mass binaries to variable mass ratios. We provide rough estimates for the higher post-Newtonian orders needed to determine the number of gravitational wave cycles, with a specified precision, that fall into a detector’s bandwidth.

by Barack, Leor and Sago, Norichika
4 pages

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

by Canizares, Priscilla and Sopuerta, Carlos F.
15 pages, 9 figures, Revtex 4. Minor changes to match published version

The description of the inspiral of a stellar-mass compact object into a massive black hole sitting at a galactic centre is a problem of major relevance for the future space-based gravitational-wave observatory LISA (Laser Interferometer Space Antenna), as the signals from these systems will be buried in the data stream and accurate gravitational-wave templates will be needed to extract them. The main difficulty in describing these systems lies in the estimation of the gravitational effects of the stellar-mass compact object on his own trajectory around the massive black hole, which can be modeled as the action of a local force, the self-force. In this paper, we present a new time-domain numerical method for the computation of the self-force in a simplified model consisting of a charged scalar particle orbiting a nonrotating black hole. We use a multi-domain framework in such a way that the particle is located at the interface between two domains so that the presence of the particle and its physical effects appear only through appropriate boundary conditions. In this way we eliminate completely the presence of a small length scale associated with the need of resolving the particle. This technique also avoids the problems associated with the impact of a low differentiability of the solution in the accuracy of the numerical computations. The spatial discretization of the field equations is done by using the pseudospectral collocation method and the time evolution, based on the method of lines, uses a Runge-Kutta solver. We show how this special framework can provide very efficient and accurate computations in the time domain, which makes the technique amenable for the intensive computations required in the astrophysically-relevant scenarios for LISA.