Skip to main content

Full text of "Second Order Gravitational Self-Force"

See other formats

Second Order Gravitational Self-Force 

Samuel E. Gralla 
Department of Physics 
University of Maryland 
College Park, MD 20742-411 

The second-order gravitational self-force on a small body is an important problem for 
gravitational-wave astronomy of extreme mass-ratio inspirals. We give a first- principles deriva- 
tion of a prescription for computing the first and second perturbed metric and motion of a small 
body moving through a vacuum background spacetime. The procedure involves solving for a "regu- 
lar field" with a specified (sufficiently smooth) "effective source" , and may be applied in any gauge 
that produces a sufficiently smooth regular field. 






1 ' 






With the promise of gravitational- wave astronomy, the 
self-field corrections to the motion of a small body have 
left the domain of pure theory to become a topic in as- 
trophysics. The principle system of interest is a compact 
object orbiting a supermassive black hole, or "extreme 
mass-ratio inspiral" . To obtain sufficiently accurate grav- 
itational waveforms for data analysis and parameter ex- 
traction, one must go beyond the geodesic approximation 
to include deviations caused by the body's finite mass 
(e.g., pj]). In fact, simple scaling arguments (e.g., 0) 
suggest that even the leading self-force correction is not 
enough, and to achieve the desired accuracy one must 
keep terms second order in the mass ratio. While much 
theoretical and computational progress has been made on 
the first order problem, by contrast very little is known 
about second-order gravitational self-force. 

Our previous work Q (hereafter paper I) established a 
rigorous and systematic approach to the motion of small 
bodies in general relativity, based on a one-parameter- 
family of solutions to Einstein's equation. Key elements 
are a far-zone limit (associated with viewing the body 
from far away) where the body shrinks to zero size and 
perturbs the external universe, and a near-zone limit (as- 
sociated with viewing the body up close) where the body 
remains at fixed size and is perturbed by the external 
universe. We developed the basic elements of the formal- 
ism to all orders in perturbation theory, but derived an 
equation of motion only to first order in the size/mass of 
the body. The basic approach was to first compute the 
far-zone metric perturbation in some gauge, and then to 
seek a smooth gauge transformation such that the near- 
zone background metric becomes mass-centered (in that 
its mass dipole moment vanishes) . The value of the gauge 
vector on the background worldline then defines the per- 
turbed position in the original gauge. In paper I we de- 
rived an equation for the Lorenz gauge motion, while in 
a later paper [4] (hereafter paper II) we derived an equa- 
tion of motion holding in a larger class of gauges. 

In the present work we identify a suitable notion of 
mass-centered at second order and define the second- 
order motion in an analogous way. However, our deriva- 
tion of an equation of motion proceeds in an entirely dif- 
ferent manner. One major change is that the approach 
is taken "in reverse" : instead of beginning with an ex- 

pression for a metric perturbation and seeking a gauge 
transformation to some mass-centered gauge, we instead 
begin with a series expansion for the general metric per- 
turbation in a particular mass-centered gauge and con- 
sider the class of all smoothly related gauges. This leads 
to a prescription (via an "effective source" method [5-7]) 
for computing the metric perturbation in such gauges, as 
well as a simple description of the motion in terms of the 
"regular field" employed in the effective source method. 
In the present paper we assume for convenience that the 
spin and higher moments of the body are negligible, but 
there should be no obstacle of principle to relaxing these 

A treatment of second-order gravitational self-force 
was given previously in [^. This approach is essentially 
axiomatic in that a number of properties that the force 
"should" have are assumed (principally, a list of ingre- 
dients from which it may be buili|j), and based on these 
assumptions a force expression is obtained. By contrast, 
our approach is fundamental in that we begin with Ein- 
stein's equation for extended bodies in the limit of small 
size, and proceed by defining perturbed position and 
computing an equation it satisfies. The approach of [2] 
also contains a serious practical drawback in that it re- 
quires the first-order metric perturbation to be expressed 
in a gauge where the (first-order) self-force is zero. Since 
a body will deviate secularly from its background mo- 
tion as it loses energy to gravitational- wave emission, this 
gauge can only remain useful for a limited time and is in- 
appropriate for calculations of inspiral (see discussion in 
section lVlI)) . By contrast, our equation of motion holds in 
a class of gauges encornpassing all possible motions. Fi- 
nally, the approach of [3| encounters divergences of both 
the "infrared" (at spatial infinity) and "ultraviolet" (at 
the particle) varieties; while the ultraviolet divergences 
are regularized, the infrared divergences are left infinite. 
By contrast, our derivation and result are well-defined. A 
recent paper on second-order perturbation theory [8| also 
involves regularizations, and concludes with an equation 

^ It is worth noting that one ingredient disallowed in |3| , the angle- 
average, does appear in an expression for the force that holds in 
a larger class of gauges [^ . 

whose mathematical legitimacy is unclearQ 

Very recently, Pound ^U^ has given an outline of a 
method to derive a prescription for computing the second 
order motion and metric perturbation of a small body. 
The approach appears to contain many features similar 
to our own, although insufficient detail was given in [lO[ 
to enable a proper comparison of his approach to ours. 

Our conventions are as follows. We forgo the abstract 
index notation and work with coordinate components of 
tensors throughout. We find this more convenient for 
discussing the perturbed position of the particle as well as 
for considering non-smooth coordinate transformations. 
Greek indices label spacetime tensor components, while 
time and space components are denoted by and mid- 
alphabet Latin indices i, 7, . . . , respectively. Our sign 
conventions are those of Wald [11| . 


The basic approach of paper I is to consider a one- 
parameter-family of spacetimes containing a body that 
shrinks to zero size and mass with the perturbation pa- 
rameter, A. We build an appropriately shrinking body 
into the family by assuming the existence of a second 
limit that is designed to maintain any such body at fixed 
size (effectively "zooming in" on it). More specifically, 
given a family of metrics g^iy(A) in coordinates (t,x^), 
we introduce a scaled metric g^^, = X^^g^iv and, for 
some time t^, scaled coordinates t = A~^(t — io) and 
ff = X^^x^ . Denoting the scaled metric in scaled coordi- 
nates by gjiv, we then have the simple relationship 

gp.p{>^;to;t,x'') = g^,.(A;t = to + At,^' = Ax*), (1) 

which relates components of the scaled metric in scaled 
coordinates to corresponding components of the original 
metric in the original coordinates. One can construct 
perturbation series either off of the original metric or the 
scaled metric; we define glU' = (1/?i!)(9a)"5^i/|a=o and 
gp]^' = (l/n!)(9A)"ffMpU=o, where derivatives are taken at 
fixed original and scaled coordinates, respectively. These 
series are referred to as the far-zone and near-zone series, 

While we will always work with coordinate compo- 
nents in the original cartesian-like coordinates (i,x'), it 
is convenient to introduce spherical-like variables r = 
^Sijx'^x^ and n' — x^/r (denoted n when representing 

^ Equation (26) of [8| contains both delta functions and a term 
written as the second-order Einstein tensor acting on the distri- 
bution /i^"^"^'. Since products of distributions are not defined in 
general, it would require further analysis to give meaning to this 
term. Since /i^'^'^* is sourced by a point particle delta function 
(within linearized theory), while point particle delta functions 
are not allowed in the full theory [^ , it would be surprising if the 
second-order Einstein tensor of /i^'®* were a valid distribution. 

a direction on the sphere). The assumptions of paper 
I (adopted identically here) give the curve r = A = 
(denoted 7) the interpretation of the lowest-order mo- 
tion of the particle, and imply that 7 (assumed timelike) 
is in fact a geodesic. This allows us to choose a coor- 
dinate system (such as Fermi normal coordinates) such 
that g*^°^ — V + 0{r^). After making such a choice, our 
assumptions give the form of the far-zone series (defined 
only for r > 0) to be 


= v 









= aoir^^ 


flu + 







= ao2r^^ 


ai2r^^ + 







= aosr^^ 


aur^'^ + 






where the (a^i^)„m (tensor component indices spuprcsscd 
above) are smooth functions of time and the two-sphere, 
anm — anm(t,n). Using equation (jT]), one may obtain 
an expression for the near-zone series in terms of the 

\(^fii^ )nm , 

g(°) = 77 + aoir-i + ao2r"^ + 003^^"^ + 0{f-^) 
g(i) = flu + ai2f-^ + ai3f"2 j^ 0{f-^) 



+ t{aoif '+ao2f ^+ao3r ^ + 0{f ")) 
= a2or^ + a,2ir + 022 + 023^^^ + 0{f^^) 
+ i(aii + ai2f-'^ + aiy.f-'^ -f 0(r"^)) 

+ -P(aoif-i -I- ao2r"^ + Sosr-^ + 0{f'^)) 

= 030^^ + aaif^ -f a32f -f 033 -I- 0{f^^) + ... 


where the i-dependence of g^"^-* (which goes up to P) is left 
unexpressed. In equation ([3]), tensor indices (scaled on 
the LHS and unsealed on the RHS) have been suppressed, 
an overdot refers to a derivative with respect to the time 
argument of the a„m, and the a„m are evaluated at t = tg. 
For example, in full notation the first line would read 

gfiv = Viiy + (a^i/)oi(io,«)?'~^ + .... 

Notice that the "columns" of equation ^ correspond 
to near-zone perturbations in equation ([S]) . For example, 
the near-zone background g^^' corresponds to the first 
column in equation ([2]), and the first near-zone pertur- 
bation g^^> is specified by the second column (stationary 
part) and the time derivative of the first column (part 
linear in t). The alignment adopted in equation ^ is 
a helpful visualization tool for the relationship between 
the near-zone and far-zone perturbation series. 

Equation ([S]) shows that near-zone background met- 
ric g^^' is stationary and asymptotically flat. Fur- 
thermore, the metric is written in adapted coordinates 
(the components are i-independent and asymptotically 
Minkowskian) , so that the mass dipole momenl[j pro- 
vides a measure of how "off center" the coordinates are. 

'^ By mass dipole moment of g'") we mean the 1/r^, (. = 1 part 
of (l/2)gQQ . Note that the mass dipole moment is contained in 
ao2, which is located at second order in the far-zone. 

In paper I we showed that a smooth first-order far-zone 
gauge transformation can always be made to eliminate 
the mass dipole moment of the near-zone background. 
Since the new coordinates are then mass-centered, the 
new perturbed motion should vanish, suggesting that the 
old perturbed motion be defined to be the value of the 
gauge vector on the worldline. This defines the first- 
order motion (in any allowed gaugqj) in terms of a far- 
zone gauge transformation to make the near-zone metric 
mass-centered at zeroth order. 

We would similarly like to define the second order mo- 
tion in terms of a far-zone gauge transformation that 
makes the near-zone metric mass-centered through first 
order. However, our lowest-order notion of mass-centered 
(vanishing mass dipole) was sensible only because the 
metric components of g^^' are always f-independent and 
asymptotically Minkowskian. It is clear from equation 
([3]) that the perturbed metric g^^'' + Xg^^'' will not nec- 
essarily satisfy these criteria. However, if a gauge can 
be found where g*^°^ -I- Ag*^^^ is in fact i-independent 
and asymptotically Minkowskian and furthermore has no 
mass dipole, then we may regard this gauge as mass- 
centered. We show below that such gauges do in fact 
exist, which defines the motion in these and smoothly 
related gauges. However, unlike in the lower order case, 
we do not show that all allowed gauges are smoothly re- 
lated to a mass-centered gauge. Instead, we simply find a 
mass-centered gauge and develop a prescription for work- 
ing in any of the (large) class of smoothly related gauges. 

The remainder of this paper is organized as follows. 
In section lllli we explicitly solve Einstein's equation in 
series in r and A to determine the general solution com- 
patible with our assumptions (to the relevant orders in 
r and A) up to coordinate freedom. In particular, this 
establishes a convenient local inertial coordinate system 
for the far-zone background metric (named RWZ coor- 
dinates), and a convenient mass-centered gauge choicqj 
(named P gauge) for the metric perturbations. We use 
our P gauge solution in two important ways. First, in 
section IIVI we use the explicit singular form of the P 
gauge solution to identify appropriate "singular fields" 
for use in an effective source prescription for comput- 
ing the global metric perturbation in smoothly related 
gauges (P-smooth gauges). Second, in section |V] we use 
the mass-centered property of the P gauge to deduce 
a simple prescription for determining the motion in P- 
smooth gauges. We do not ascribe any fundamental sta- 
tus to our particular class of P-smooth gauges, and in 
section lVTl we discuss how the paper could have proceeded 
(identically) if an alternative class of gauges were used. 

In section rvm we discuss incorporating our results into a 
formalism for long-term waveform generation. Finally in 
section rVIIII we summarize our prescription for comput- 
ing the second-order motion and metric perturbation of 
a small body. 


We now explicitly construct a gauge that is mass cen- 
tered in the sense discussed in section [Hi i.e., a gauge in 
which the near-zone metric is f-independent and asymp- 
totically Minkowskiian and through first order in A. We 
will call this gauge the P gauge, where the P stands 
for "particular" , in order to emphasize that other mass- 
centered gauges could have been chosen. (We discuss 
this freedom in more detail in section IVll ) We perform 
our computations using the near-zone perturbation se- 
ries. While it is necessary to proceed only to first order 
in A to establish that a gauge is mass-centered, the mass- 
centered coordinate choice influences the form of terms 
at higher order in near-zone perturbation theory, many 
of which will be needed for the later analysis of the pa- 
per. In performing our near-zone calculations, we will in 
fact have to proceed through third order in A. 

We begin our computations with the background near- 
zone metric. Since this metric is stationary and asymp- 
totically flat, it is characterized by multipole moments 
[13 | . We treat a body with negligible spin and higher mo- 
ments, and therefore take the spin and higher moments 
of this metric to vanish. Thus the near-zone background 
metric is simply the Schwarzschild (exterior) metriqfl for 
all time tg. While in principle the mass may depend on 
time to, in paper I it was shown to be constant. We label 
the mass by M and choose Cartesian isotropic coordi- 
nates for g^^^ (at all time to)- 

Since the metric components of g^'^' are then indepen- 
dent of to, by (121) the near-zone perturbation g^^^ must be 
independent of time t. Furthermore, equation ^ shows 
that the perturbation is asymptotically constant. Stan- 
dard Schwarzschild perturbation results [13, 14] then im- 
ply that its only physical effect can be to perturb the 
multipole moments of g^^' . In line with our choice to 
consider a body with no spin and higher moments, we set 
the perturbed spin and higher moments to zero. While 
in principle there could be a perturbation to the mass, it 

* Since the metric perturbations are singular, non-smooth gauge 
transformations are allowed, the definition of motion in one gauge 
does not automatically define the motion in all other gauges. 

^ We refer to finite-A coordinate transformations that preserve the 
metric components of the background metric as "gauge transfor- 
mations" . 

® While for convenience we will make statements about an "entire" 
metric, it should be borne in mind that we only require that such 
statements hold to the orders explicitly displayed in the paper. 
(These orders have been carefully chosen for consistency with all 
statements made.) For example, in this case we say that the 
metric is Schwarzschild, but in fact we only require it to match 
Schwarzschild to 0{f~^) (see equation (|3|). Thus we in fact 
only assume that the spin and quadrupole moments vanish — 
the effects of higher moments are automatically negligible at the 
present level of approximation (that is, these moments do not 
appear at the orders in A pursued in this paper). 

was shown in paper II that the perturbed mass does not 
evolve with time. Therefore this quantity may as weh 
be "renormalized" into the background mass, or equahy 
acceptably simply set to zero. We will set the perturbed 
mass (and higher moments) to zero. With these physi- 
cal choices the first near-zone perturbation is pure gauge, 
and we choose it to vanish. 

With the above choices the near-zone metric agrees 
with the ordinary (mass-centered) Schwarzschild metric 
in Cartesian isotropic coordinates through first order in 
A, and therefore fits our notion of being mass-centered. 
Explicitly, the perturbation series ([3]) is now given by 

g(°) - ry + A^tDf-^ + A4(2)--2 ^ _^(3)^-3 ^ ^(--4) 




a2or^ + a2ir + 022 + Oif"^) 

aaor^ + agif^ -I- 032^ + 0{f°)) 
+ i{d2of'^ + a2if -f 022 + 0{f-^)) 


where A^(") is the 7i*^-order term of the Schwarzschild 

metric in Cartesian isotropic coordinates. 



with t^ = (-1,0,0,0). 

We now consider the second-order near-zone metric 
perturbation, g'^'^\ Since the first near-zone perturba- 
tion vanishes, the second perturbation satisfies the lin- 
earized Einstein equation off of the Schwarzschild back- 
ground. Equation ([B]) shows that our perturbation is i- 
independent to the relevant order, so that we may re- 
strict consideration to stationary solutions. We use the 
Regge- Wheeler- Zerilli (RWZ) formalism [ij, [lj|, where 
one decomposes the perturbation into a sum of (tensor) 
spherical harmonic modes labeled by azimuthal number 
£. For I — Q and I — 1 modes, the general stationary 
solution (up to gauge) has r — J' 00 asymptotic behavior 
of f~'^^+^\ while for ^ > 1 the general stationary solution 
(up to gauge) is given by a linear combination of functions 
behaving as f~(^+i) and r^. From these considerations, 
comparison with equation ((6]) shows that the general so- 
lution (up to gauge) for g^'^> of our form is pure i — 2 
to the displayed orders. This solution is characterized 
by ten constants (one for each m-mode of each parity), 
which are conveniently represented as two constant sym- 
metric trace-free (STF) tensors in three-dimensional Eu- 
clidean space. (Excellent reviews of the STF approach to 
spherical harmonic decompositions are found in ^15|,[l6J.) 

In our computations we employ the RWZ formalism as 
presented in [l7l - [ig |. translating the results into Carte- 
sian isotropic coordinates and STF language. We use the 
closed-form expressions for the stationary master func- 
tions given in [20,]. We find that the general solution 
(in Regge- Wheeler gauge) for our 5^^' may be written in 
terms of two arbitrary STF tensors £ij and Bij by 

5^0^ = £,,7i''n^(-r2 + 2Mf + ^M^) + 0{f-^) (11) 

-g"^ = e,,,n=B\n\'^f^ + '^Mf - \m^) -f 0{f-^) (12) 



f^ - AMf - 5M^) 

+ ninj{2Mf - AM^)] + 0{f-^), (13) 

where eijk is the Cartesian Levi-Civita symbol. The STF 
tensors Sij and Bij may depend on the time to at which 
the near-zone limit is taken, but are independent of the 
near-zone background coordinates i, x* . Equations (jlip - 
(1131) determine the unknown functions 020,021,022 in a 
particular gauge. 

We now turn to the third order near-zone perturbation, 
g^'^\ This perturbation also satisfies the linearized Ein- 
stein equation off of the Schwarzschild background (on 
account of the vanishing of the first perturbation). We 
use the RWZ formalism to find the general solution for 
g^^' consistent with equation ([T]). From general consid- 
erations of the sort discussed for the second-order per- 
turbation, above, this solution contains only £ = 2 and 
£ = 3 modes (up to gauge). The i-dependence of the 
perturbation is fixed entirely by g^'^' (see equation (O). 
For the part independent of i, the RWZ formalism yields 

ffoo \t= 

^0 = £,jkn'n^n''{--f^ + -Mf^ + —M^f) 

9io It 

1-3 2,^_2 7 
^„fcn'ri-'n"'( r^ + -Mr H . 

+ SSiju'n^ {-f'^ + 2Mf) + 0{7^) 

,2 , 4,, o 2 


|,=o = e.,fen^S^,„n'n"(-r-3 -|- -Mf^ - -M^f) 


+ n,£kin''n''{~-r 




9 1] \t-- 

+ e,,kn'5B\n\-f'' + -Mr) + O(f0) (15) 

+ n,n,{-Mf'^ --M'^f)] 

-,k fol ^m 


-fn(,e,)fein"6'„n"(-f^ + 2Mf 

2 -M^f) 
12 ' 

+ SSkin'^n^ [S,j (-f^ - 4Mf) + 2Mfninj\ 

+ 0{f^), (16) 

where we have introduced STF tensors Eijk-^ Bijk, S£ij, 
and SBij , and the overdot indicates a derivative with re- 
spect to to. Equations ((T4l) - (|T6l) determine the unknown 
functions 030,031,032 in a mass-centered gauge. We have 
therefore constructed the perturbation series in a partic- 
ular gauge to the relevant orders. 

We note that the strategy of using the RWZ formahsm 
in the near-zone has been employed before, most notably 
in [2l|, |22|- The main difference (besides the different 
coordinate choices) is that [211, ^^ impose boudary con- 
ditions of regularity at the horizon of a black hole. By 
contrast, we treat an arbitrary body and impose bound- 
ary conditions only at infinity. Indeed, it is a result of 
our computations that no interior boundary conditions 
are necessary to fix the form at the relevant orders, which 
corresponds (after the analysis of the paper) to the result 
that the motion of the body is independent of its detailed 
composition at the perturbative orders considered. 


A. Far Zone Expressions 

Our formulae for the near-zone metric determine the 
displayed coefficients a„m in equation ([3]). These coeffi- 
cients may then be used to reconstruct the far-zone se- 
ries via equation ([2]). This gives the general far-zone 
solution to the Einstein equation to the relevant orders 
in r and A in a particular gauge. However, rather than 
using this gauge as our P gauge, we instead first make 
a particular first-order (far-zone) gauge transformation, 
which is designed to make the first-order metric pertur- 
bation satisfy the Lorenz condition, while preserving the 
mass-centered property. Using this gauge as our P gauge 
ensures that the P-smooth class includes Lorenz gauges, 
which facilitates comparison with previous work using 
Lorenz gauge, as well as allows Lorenz-gauge numerical 
results to be incorporated into second-order calculations 
based on our prescription. However, we emphasize that 
(unlike in some previous work) the Lorenz gauge plays 
no fundamental role in our analysis. Our choice to per- 
form an additional Lorenz-motivated gauge transforma- 
tion (i.e., our Lorenz-motivated choice of P gauge) affects 
the details of many complicated formulae throughout the 
paper, but has otherwise no influence on our prescription. 
We specifically discuss alternative choices of P gauge in 
section I VII 

0. Notation In the remainder of the paper we will re- 
fer to the far-zone background, first perturbation, second 
perturbation, and third perturbation as g, h, j, and fc, 
respectively. (That is, g = g'^^\ h = g^^^ and j = 5^^\ 
k = 5^'^-'.) This facilitates the introduction of many 
necessary new superscripts and other modifiers. In or- 
der to avoid any potential ambiguity, the one-parameter- 
family of metrics will always be referred to with the A- 
dependence indicated, g{X). 

1. Background. Reconstructing the background metric 
g from the near-zone solution gives 

(000)20 == -EijU^n' 

(aio)20 = -€ijkn^B\n^ 

{aij)20 = -SijSkin'^n^ 

(000)30 = --Eijk-n'^n'n^ 

2 2 • 

(ajo)30 = -e^jkn^B^^n^n"" - -n,£kin^-n} 





We may now interpret the STF tensors f^j, Bij, £ijk and 
^ijk by computing the Riemann tensor of g. It is then 
straightforward to see that 

£ij — i?OJOj I7 

BiA — — -e ,i?i 


£ijk — V(fci?|o|i|o|j)l7 



Thus our coordinate system for the background metric 
expresses an arbitrary vacuum metric in terms of the cur- 
vature on an arbitrary timelike geodesic x^ — 0. We refer 
to these coordinates as RWZ coordinates, after the use of 
the Regge- Wheeler- Zerilli gauge in solving the near-zone 
equations. Our metri c ag rees with the Fermi normal co- 
ordinate metric (e.g., 1231 ) to 0{r) and with the Thorne- 
Hartle-Zhang metric [23 - l24| to 0{r'^). The form of the 
metric at 0{r^) appears to be new. 

2. First Perturbation. Reconstructing the first far- 
zone perturbation gives a specific expression for h in 

terms of £,;.,, Bin, £. 

y ' 

^ij ^ ^ijkj i-^ijkj 


SSij, and SBij. 

Instead of 

adopting this expression as our P gauge, however, we 
first make a particular gauge transformation generated 
by the gauge vector 





Vi = r^ {2£ijn^ - ni£jkn^n'') 

+ r^ ( -£ijkn^n'' - -n^£kim.n''n^n'"'' - 

- i;etjkn''n^&i\ . 




This gauge transformation is designed to make h^ sat- 
isfy the Lorenz condition, as may be checked by direct 
computation using the formulae below. As will be dis- 
played explicitly in equation (I55|) . below, the transfor- 
mation does not affect the mass-centered property of the 
coordinates. After performing the transformation we de- 
note the resulting perturbation by h^ , which is given by 

g^V + a2Q{ty + asoity + 0{r^), 


M'^^'>r-^+a2ir + a3ir^+0{r^) (30) 


(aoo)2i = 2M£ijn'n^ 


(010)21 = ^Me.^kB in 

(aij)2i = -2M{5ij£kin''n' + 2£ij) 



("00)31 = -M£,,kn'n^n^ 

(«fo)3i = \Me,,kn^B\yn^ 

2 • 
- -M{£ijn^ - n.^£kin^in}) 

(ag)3i - M T-^^ffei^n'^n'n" - 2£,,fen'= 

+ ^Af (7i(,ej)fe,7i'=^'„7i™ - 2B(/e,)fcin'=) (37) 




(0^)31 = -S£ijn^n^ 

(0^^)31 - -e^jkn'SB\n^ 

(0*^)31 = -SijS£kin'^nK 


We have spht 031 into S and H pieces in order to make a 
similar spht h^ = h^ + h^ + 0{r^), 

h' = A^ilV-i + [a^^Ur + {a^^^h^ + 0{r') (41) 


(03ir^ + 0(r3). 


The reason for this spht wih become clear when the "sin- 
gular field" h^ is employed in the following section as part 
of a prescription for computing the metric perturbation. 
The guiding principle is that h^ be determined by the 
background metric (containing only £ij,Bij,£ijk,Bijk, 
and not the unknown 6£ij and 6Bij) and that the re- 
mainder h^ be C^ . There are many other choices besides 
ours that satisfy these properties, and we could equally 
well have made these choices. Our choices also have the 
additional properties that h^ and h^ separately solve 
the field equations to the displayed orders jj For future 
use, we relate S£ij and SBij to h^ by computing the lin- 
earized Riemann tensor of h^ , finding (cf. equations (p4)) 
and (EU) 






'Ojkl I 


''' This is most easily seen by checking that h^ is a solution, a 
computation that requires only the leading order term <; = 77 of 
the background. Since the sum h^ = h^ + h^ is by construction 
a solution for r > 0, it follows that h^ is also a solution for r > 0. 

Here RliJpcr[h] is defined as a function of a symmetric 

rank 2 tensor hf_,^ by R}}Jpa[dxg\\=o] = dxR^„pa{X)\\=o 
for a smooth one-parameter-family g{X). 

3. Second Perturbation. We reconstruct j from the 
near-zone solution and take into account the effects on j 
of the first-order gauge transformation, equations ([28])- 
([29]) . The second-order perturbation is given by 

f = M^^^r-^ + 022 + aj.2r + 0{r^) 


(000)22 - -iM^£^Jn'n^ 
(ajo)22 = —M'^€,jkn'n'-B\ 

{a,j)22 = SM^£^J - M^£kin''n^ (6% + n,n.^ 

5 3 

(000)32 = --£ijfen'n%'' -f -M8£ijn'-ri? 

2 2 • 

(0,0)32 = -e^j],niB\„{a^n^ - -n^£kln^n^ 

- ^Me^Jkn'SBW 

iazj)32 = ~-;;5ij£kimn''n^n"^ + -n(iej)kin''B\j^vJ 

- M5^j5£kin^nK 




Note that we have not made a second-order gauge 
transformation, analogous to the transformation ((28)) - 
([29| made at first order. The first-order gauge trans- 
formation was designed to make h satisfy the Lorenz 
condition, which was desirable because of the long his- 
tory of use of Lorenz gauge in both theoretical and com- 
putational work at first order. For second-order pertur- 
bation theory, the relevant previous work is [13, [2J, |23] , 
where the Lorenz condition was imposed on the second- 
order metric perturbationlj In the interest of compar- 
ison, we have investigated whether this condition may 
be imposed within our formalism. We have found that 
it appears necessary to introduce r log r terms in to the 
metric perturbation in order to impose this condition. 
This directly violates the metric form required by our 
assumptions (equation ^), and, if allowed, would lead 
to A log A terms in the near-zone series by equation ([l} . 
Since a smooth near-zone perturbation series is an essen- 
tial ingredient in our justification (see paper I) of the rel- 
evance of our perturbation series to small (but extended) 
bodies, we take the viewpoint that such a far-zone gauge 
is too singular to sensibly describe a small body, at least 
within our current approach|3 

Note that when specialized to a fiat background spacetime, 
this differs from the harmonic gauge condition used in Post- 
Newtonian theory by terms involving the first-order perturba- 

The appearance of log terms at second order in the Lorenz gauge 
was also found in the gauge-relaxed formalism of |25| . This has 
an analogous singular effect on the near-zone metric; this effect 
is not discussed. 

4- Third Perturbation. Reconstructing the third-order 
metric perturbation from the near-zone yields 

fc^=A^(3)r-3 + 0(r-i). 


While it would have been straightforward to compute the 
terms proportional to r~^ and r° (i.e., 023 and 033) from 
our near-zone expression (plus the effects of the first order 
gauge transformation, equations (l28)) - (l29l) ). these terms 
are not relevant for our analysis. Note that the first-order 
gauge transformation, equations (P5| - (P5|) . has had no 
effect on the displayed orders. We now collect the results 
of this section in the form of equation ([2]) , 


= r] + + 020^^ 

+ a,or^ + Oir^) 


^X(i)j.-i + + 021?' 

+ asir'^ + 0{r^) 


= X(2)^-2 +0 + 022 

+ a32r + 0{r^) 


= 7W(3)r-3 + + 0(r-^) 


where the a„m are now given explicitly by the formulae 
in this section. In this form it is easily seen that the near- 
zone background (first column) is Schwarzschild and the 
first near-zone perturbation (second column) vanishes, so 
that the P gauge is indeed mass-centered. 

5. Summary of Results. Equation (|53p , together with 
the preceding expressions for the a„m, is the main re- 
sult of this section. This expression provides a series ex- 
pansion in r for general zeroth, first, second, and third- 
order metric perturbation subject to our assumptions, 
expressed in a particular mass-centered gauge, known as 
P-gauge. For use in the following section, we have also 
isolated off a particular singular portion of h^ , denoted 
h . We have used the tensor analysis package xTensor 
|2a | for the software package Mathematica [27| to per- 
form many of the computations in this section. We have 
verified by direct computation (taking several hours on a 
personal computer) that the metric g+Xh^+X^j^'+X'^kP 
satisfies Einstein's equation to the relevant orders in A 
and r. 


In the previous section the general solution for the met- 
ric g{X) was determined in series in r and A, subject to 
particular coordinate choices. Since the motion is also 
known in these coordinates (it is given by the coordinates 
of the background geodesic 7), we have at some level de- 
termined the general solution to our problem. Of course, 
this general solution is of no use in practice, since it con- 
tains undetermined parameters (with no physical inter- 
pretation) and gives the metric only locally near r — 0. 
Nevertheless, this analysis has revealed the structure of 
the general solution near r = 0, which will allow us to 
develop a prescription for obtaining the global metric per- 
turbation in a P-smooth gauge in situations of physical 
interest, as described below. 

Given our assumptions on the one-parameter-family. 

Einstein's equation implies in the far-zone that 

GW[/i]=0 (forr>0) (54) 

Gi'}[h]^0 (forr>0), (55) 


where G^^-* and G^^-* are the first and second order Ein- 
stein operators, respectively. When combined with the 
assumed form of the metric perturbations near r — 
(equation ([5])), these equations provide the complete de- 
scription required to compute h and j in a given situation 
of interest (i.e., once suitable initial and/or boundary 
conditions have been prescribed). In practice, however, 
it may be difficult to ensure that a numerical solution 
have the correct divergent behavior near r = 0. Further- 
more, it is far from obvious how to ensure that the metric 
perturbation will be determined in a gauge for which we 
define the motion. 

A solution to both of these problems is to use our 
knowledge of the general P-gauge series solution near 
r = to "regularize" the differential equation. One sim- 
ply subtracts off the known singular behavior and evolves 
the regular remainder. This type of numerical technique 
was introduced into the field of self-force computation by 
[5I, la] , and is now generally known as the "effective source 
approach" Q. At first order, our approach is equivalent 
to the standard approach, except that we are not re- 
stricted to the Lorenz gauge, and instead allow the use 
of any gauge condition that gives rise to a sufficiently 
regular "regular field" . Our presentation of the method 
differs in that we do not make use of 5-function sources, 
instead working directly with our assumed form of the 
metric perturbation for r > Orn 

A. First Order 

In the previous section we constructed the general so- 
lution for the first-order metric perturbation in series in 
r in a particular gauge. We refer to this gauge as the P 
gauge and denote the perturbation by h^ . In a general 
smoothly related gauge, the metric perturbation is given 

= h^ + h" -C(^g + 0{r^), 



where the split of h^ into h^ and h^ was introduced 
in equation (|4l])EB Recall that h^ = 0{l/r) is a sin- 
gular approximate solution to the linearized Einstein 

^^ In paper I, we proved that our assumptions in fact imply a delta- 
function source for h (regarded as a distribution) at first or- 
der. From the point of view of developing an effective source 
description from our assumptions, such a delta- function descrip- 
tion would appear only as an unnecessary intermediary. 

^^ The error terms in equations I I56II and I I57II are redundant with 
those in the definitions of h^ , h° and h^ , but we include the 
error terms as a reminder of the local nature of h^ , h° and h^ . 


equation specified by the background curvature tensors 
£ij,Bij,£ijk,Bijk, while h^ = 0{r^) is a C^ approxi- 
mate solution given in terms of undetermined parame- 
ters SSij and 5Bij, which encode its Riemann curvature 
via equations (|4S| and (|44l) . The sum, h^ , represents 
the general solution with a particular gauge choice, up 
to 0{r^) errors. We emphasize that our /i^, h^ , and h^ 
are given only as approximate solutions near r — Q] none 
of these quantities has a finite-r definition from which 
the series expansions emerge. This is in contrast to the 
singular field of |2^| , which is defined in a normal neigh- 
borhood through the use of Hadamard Green's function 
techniques in the Lorenz gauge. We not checked if our h^ 
agrees with a series expansion of the Detweiler- Whiting 
singular fieldFn 

Implementing the effective source approach requires 
choosing an (arbitrary) extension of h^ to the entire man- 
ifold (minus r = 0). We will distinguish extended quan- 
tities with a "hat" : Let h^ denote an arbitrary function 
on the manifold (minus r — Q) such that h^ agrees with 
h^ (equation (J4T|) ') to all displayed orders in r (i.e., to 
0(r^)). We then define a global "regular field" h^ in 
terms of the metric perturbation h by 

h" = h- h'^ 


Plugging equation ((58|) into the linearized Einstein equa- 
tion ([M)) gives 

G(i)[/i«] = _G(i)[/i' 


(for r > 0). 


By construction we have that G'^-'[/i'^] is 0(r), so that 
the right-hand side is in fact 0{r). Thus the "source" 
— G^^) \Ji^] is G°, and we may in fact drop the requirement 
that r > 0. We may then write the first-order equation 
as simply 

G(i)[/i^i = 5(i), 


where the G° source S^^' is given throughout the mani- 
fold by 

5(i) = -G(i)[/i^]. 


The logic of the above argument has been that if one 
has a metric perturbation h satisfying equation (|54p and 
in a P-smooth gauge (equation (1551) ). then the effective 
source equation ((6T|) holds. In practice, we want to pro- 
ceed in the reverse direction: we wish to solve equation 

^^ In [231 a singular field was constructed in a manner similar to 
ours (but using different coordinate ciioices and notation). It was 
then claimed that this singular field agrees with the Detweiler- 
Whiting singular field up to errors of 0{r ). It seems likely that 
our singular field agrees with that of 22] (and therefore with the 
Detweiler- Whiting singular field) at this order. 

(1611) and thereby obtain an h satisfying equation (1541) and 
in a P-smooth gauge. Retracing the steps of the argu- 
ment in reverse, it is clear this will hold provided the 
solution h'^ of equation (|6ip is G^ at r = 0[3 Obtaining 
such an h^ will depend on the initial and/or boundary 
conditions chosen, as well on as the choice of gauge. 

We first discuss the choice of initial and/or boundary 
conditions for h^. We view the specification of a "phys- 
ical situation of interest" as a choice of initial and/or 
boundary conditions for the metric peturbation h. In 
principle, one would first determine such conditions in a 
P-smooth gauge and then infer the relevant conditions on 
h^ — h — h . In practice, determining appropriate initial 
conditions for h is likely to prove difficult, even without 
the added requirement of using a P-smooth gauge. Faced 
with difficulty determining appropriate initial data, the 
usual strategy is simply to choose inappropriate initial 
data and evolve in the hopes that at a later time (af- 
ter "spurious radiation" has left the system) the solution 
will nevertheless resemble the desired physical situation. 
We suggest that one employ this strategy at the level of 
the regular field h^, where one could simply choose triv- 
ial initial data (or a suitable generalization should triv- 
ial initial data confiict with any gauge conditions used). 
Effective source calculations made with the scalar wave 
equation [3] suggest that this strategy will prove effective 
in the gravitational case as well. 

We next discuss the choice of gauge. Since h^ is re- 
lated to the metric perturbation h by addition of a fixed 
quantity h , the usual arguments that h and h + C^g 
represent the same physical configuration imply that h^ 
and h^ -f C^g represent the same physical configuration. 
In particular, any gauge condition that is "allowed" for 
h will remain "allowed" for h^. For example, it is well 
known that one may impose the Lorenz condition on a 
smooth perturbation h, V^i/^y = with capitalization 
denoting trace-reversal, H^y ~ h^^^ — {l/2)g^j^^h. Sinn- 
larly, one may argue identically that it is always possi- 
ble to impose the Lorenz condition on the regular field, 
V^H^jj = 0, where capitalization denotes trace-reversal. 
In this case equation ([6T|) becomes 

i?^,[/i«] = ^(i). 


where £"^1/ is the Lorenz-gauge linearized Einstein tensor 
(a well-studied hyperbolic wave operator on g), 

E^,,[h] = VV^H^, - 2P" /iJ„ 



where the capitalization of the arbitrary perturbation h 
represents trace-reversal. Since _E^,y is a hyperbolic wave 
operator, it is expected that the G° source S*^^-* will give 

^^ It is clear that C^ solutions exist by the existence of P-smooth 
gauges, proved by construction in the previous section. 

rise to a C^ solution ^^13 and therefore that the gauge 
condition V^HJ^^, does in fact provide a metric perturba- 
tion h — h^ + h^ in a P-smooth gauge. Note, however, 
that this gauge differs from "the Lorenz gauge", which 
refers to the Lorenz condition on the full pertubation /i, 
V'^i/^i/ = 0. If one desires this condition to be satis- 
fied, one instead needs to enforce \/^H^^^ 
which case equation (I61[) becomes 

ViJ-^ , in 





where capitalization denotes trace-reversal. We see that 
the failure of h to satisfy the Lorenz condition appears 
as an extra effective source for h^. We have chosen our 
h^ to satisfy the Lorenz condition, V^iJ^^ — 0{r^), so 
that the failure comes only from the choice of extension, 
and the the right-hand-side of ([65|) remains C°. In par- 
ticular, the solution h^ should be C^, so that the Lorenz 
gauge is P-smooth. (Indeed, we chose the P gauge and 
hence the P-smooth class precisely so that Lorenz gauges 
would be included.) However, we emphasize that while it 
may be useful to use the Lorenz gauge to compare with 
previous work or to determine a first-order- accurate long- 
term evolution via a particular proposed prescription (see 
discussion in section rVII|) . for the purposes of determin- 
ing h there is no fundamental reason to prefer one gauge 
over another. 

In implementing the effective source method, above, 
we have made a convenient choice of "singular field" h^ . 
However, we emphasize that many other choices could 
have been made, with equivalent results. In particular, 
one may modify h by the addition of any given smooth 
function /. In this case the effective source —G^-^^[h^] 
will remain C'^ (though it will no longer be 0{r), since 
h^ is no longer a solution to all orders considered), and 
the full metric pertubation h = h^ + h will remain the 
same, provided the appropriate initial/boundary/gauge 
conditions for the corresponding new h^ (modified by 
— /) are chosen. However, while we delay a systematic 
discussion of the motion until section |Vl we note here 
that an advantage of our particular choice of h^ is that 
the first-order motion may be described as geodesic in 
the peturbation h^ (as in the original treatment of |28|). 
The basic point is that, as may be seen from equation 
dnH) with ^ = and equation (gl]), we have h^ = 0{r'^) 
in the P gauge, i.e., the regular field and its first deriva- 
tive vanish on the worldline. The statement that the 
perturbed motion vanishes (together with the statement 
that the background motion is geodesic) may then be ex- 
pressed equivalently as the statement of geodesic motion 

^* While general theorems on wave operators (e.g., prop. 7.4.7 of 
[29ll ) would guarantee only weaker regularity of the solution, ex- 
perience with the effective source method for scalar wave opera- 
tors [tII shows that sources of our type do in practice give rise to 
sufficiently regular solutions. 

in g + A/i^, which, as a covariant statement, will hold 
in any smoothly related gauge. This argument is given 
more formally and explicitly in section |Vl below. 

B. Determination of 5£ij, SBij, and ^^. 

Our next task is the identification of an appropriate 
singular field at second order. At first order, the singu- 
lar field was found by noting that the unknown tensors 
{6£ij , SBij } appeared in the P-gauge perturbation only in 
a smooth way, so that a singular part depending only on 
the known tensors {£ij,Bij,£ijk,Bijk} could be chosen. 
Furthermore, since smooth gauge transformations affect 
the metric perturbation only by addition of a smooth 
term, this choice of singular field guarantees that h — h^ 
is regular in all P-smooth gauges. At second order, how- 
ever, the unknown tensors {SSij, SBij} do appear as part 
of singular terms (see equation (j45p . where the ra32 terms 
are not differentiable) . Furthermore, an identification of 
a singular part, j^ , oi j^ does not guarantee that j — j^ 
is regular in all P-smooth gauges, since smooth gauge 
transformations change the second-order metric pertur- 
bation by a singular term, C^h^ (see equation (|A7|) and 
recall that h^ is singular). To correctly identify a sin- 
gular part of j will therefore require expressions for all 
of the unknown quantities {SSij, SBij, S,fj.} that appear in 
the expression for the general P-smooth metric pertur- 
bation, h = h^ — C^g. 

The relevant question is the following: given a pertur- 
bation h va a. P-smooth gauge (imagined, e.g., to have 
been numerically computed by the prescription given 
in the previous section), how can we express this per- 
turbation as h ^ h^ — C^g for some {S£ij,5Bij,S,^_i}l 
(We remind the reader that h^ is constructed from S£ij 
and SBij.) Below we find that there is precisely a ten- 
parameter freedom in the choice of {S£ij,SBij,£^^} that 
specifies a decomposition of the form h = h^ — C^g, 
and give a prescription for computing these quantities 
in terms of an integration of transport equations along 
7. As shown therein, knowledge of ^ determines S£ij and 
SBij, so that we may view the ten-parameter freedom 
in the decomposition as a ten-parameter freedom in the 
choice of ^. Since different such choices lead to different 
second-order metric perturbations but (by construction) 
preserve the first-order perturbation h, this freedom cor- 
responds to the infiuence of first-order gauge freedom on 
the second-order metric perturbation. The freedom in 
choice of SSij, SBij and ^ (at fixed h) may be viewed as 
first-order gauge freedom that manifests only at second 
order in the metric components. 

Since h is a specified function of known quantities, 
we may without loss of generality consider the smooth 
vacuum perturbation h^ = h — h . (Since we will work 
with the regular field only locally near r = in this 
section, we drop the hat in its notation.) From equation 
(l58l) we have 


C^g + 0{r^), 



where we remind the reader that h^ is a simple function 
of 5£ij and SBij , given by equation P^ . To show that 
S£ij and SBij (and hence h^) are determined by ^, we 
compute the Hnearized Riemann tensor of h^, yielding 


Equation (|75|) shows that solutions to equation ((74)) are 
determined everywhere by a choice of ^ and V^ at a sin- 
gle point. Equation (|74|) restricts this choice to a ten- 
parameter-family (such as "Killing data" ^^ and Vj^f^]). 
We now show constructively that all such choices lead to 
'^aPiSl^^'iyi ^ "^V ; solutions to equation ([71)) . 

-Ra/37J^^] - ^^RafSjS + 0{r). (67) Since r(i)[/i^] is 0{r), equations ^ and dTSj) give for 

''l-4'j,J/»'']+4'j,J-A5] + 0(r) 

(The definition of Rfj^upa [h] is given below equation (04]).) 
In the second line we have used the covariance property 
F'^ ''[C^g] — C^F^^\ holding for any covariant function of 
the metric F[g]. Using equations (03]) and ()i^ we then 



^2 i [^Ojka'T' \ 



where our Lie derivative expressions refer to components 
of Lie derivatives of the (background) Riemann tensor 
(rather than some kind of derivative of a component). 
Thus S£ij and 6Bij (and hence h^) are determined by 
^ and its first derivative on the worldline. In particular, 
introducing tensors A^ and B^^, defined along 7 by 


~ ?/^ 7 


B^u = 

= (v.Cm)!- 



we have 

^^'= = 2 


-d,d,h^^ + 25o9(./if)o - 


+ ^QO^ij 

A^ and _B^^ that 



2 V''l7 

u"V„A„ = B„„u' 

v„B^, = i?^,„;3u"A^ - M"r«,[/i^]| 


Equations ()75)) - (|75)) give transport rules for A and B 
along 7, while equation (|77l) gives a constraint (which 
is preserved by the transport.) In the RWZ coordinates 
these may be written 




Ai — Bio 
B,o = -£^,A^ - d^h^ 

ly- ^U'-OJ T o^'i^OO 

%,]=6,,,S',.A^ + %/ij|o, 


2i?fcoe'(,a,)' -t- T^A^e ,(j6,)' + 2Bki,£_ 



A. 5,-, 


SB^, = 





S'^'BkiB,, + -CS, 


2i?H,B,, '^ - Ao4 + TjAfefi,/- 


^ H'T'Ok^j 




J "feO 






(Note that while SBij is not given above in manifestly 
symmetric form, one may easily confirm its symmetry 
using the fact that h^ is a vacuum perturbation.) 
We now regard equation ()66p as an equation for ^, 

2V(^C.)=C-C + 0(r'), 


where h^ is constructed from ^ via equations ()68I) - (|M)) 
(equivalently dTDjl-dZS))) and (j^. Taking a derivative 
and employing manipulations normally used for Killing's 
equation (e.g., appendix C of [III), we have 

VaVpij + RpJ^s^~r['^ 

,c.p[h''-h"]+0{r'), (75) 
where r'^^^[ft,] is the perturbed Christoffel symbol with a 
lowered index. 


" 7Q/3 


W = n i^ahpj + Vpha^ - V^/la/j) 


where the overdot denotes a t derivative and /i^^ and its 
derivatives are evaluated at a;* = (i.e., on 7). Equa- 
tions (|80))-(l83)) (together with ([77))) determine A^ and 
Bfii, given a choice of initial data for {A^, Sio, 5[ij]}- 
The reader may recognize the last two terms of equa- 
tion (|5^ as the self-force on the particle, here taking 
a "perturbed geodesic equation" form. As discussed in 
more detail in our systematic treatment of the motion 
in section |V] below, our definition of motion implies that 
^■(i)^ = ^^1^ = A^, so that the transport equation for A^ 
is in fact the first-order equation of motion. However, in 
the present section we confine ourselves to the derivation 
of a prescription for computing the first and second order 
metric perturbation, for which the interpretation of f^|-y 
as giving the motion is entirely irrelevant. 

Given a choice of {A^^, Bio, Buji} at some point along 
7, equations (|5n)) - ([55)) determine these quantities every- 
where on 7. We now imagine that a choice has been 
made, so that A and B are known along 7. This deter- 
mines S£ij and SBij via equations ()72p - ()73p (equivalently 
dSHI-dni)) and hence h^ by equation (02). The right 
hand side of equation ([75)) is then "known" in terms of 
the value and derivative of ^ on the worldline (i.e., in 
terms of A and B) , so that we may determine ^ to higher 
order in r by expanding the left hand side in r and equat- 
ing orders in r. After some effort we obtain 



+ 0{r^) 

1 \ J.2 

SfcOt jllJl 

kt jiOt 




-^BjaeumBk" + —AjeumBk"' - -AQtij'^Bkira - i;Bjk£ii - Bij£ki - -A.fjfc; ) x'x'^x^ + Oir"), 

-Idkdihf^ + ^dAhfk - ^h^jSki + leurnhlljBkA x^x'^x' 




where h^ and its derivatives are evaluated on 7. We 
then check by direct computation that the above for- 
mula does give a solution to equation ([74]) (and not just 
([75]) ). provided that A and B satisfy the transport equa- 
tions (|80 |) - (j83| and that h^ is a vacuum perturbation. 
Thus equations ((5(1)) - ([55]) provide a ten-parameter-family 
of solutions for ^ to equation (|74l) and hence ((56|) . Since it 
was already shown that the general solution is at most a 
ten-parameter-family, the general solution is in fact a ten- 
parameter family, and all solutions may be constructed 
this way. 

The main results of this subsection are equations ([72]) , 
([73| . and (|84|) - (|85l) . which give expressions for 5£ij, SBij 
and ^^ in terms of an integration of the transport equa- 
tions (l80)) - (l83t for A^ and B^. We have used the tensor 
analysis package xTensor [26| for the software package 
Mathematica [27| to perform the extensive computations 
of this subsection. 

terms are regular, an appropriate singular field is thus 





One may now straightforwardly combine equations (j45p , 
((30l) . dill) -dial), and dUll-dlSl) to produce an expression 
for j'^ in terms of £ij, B^j, £ijk, Bijk, A^, and S^^,. This 
expression is given in equations (jB4P - (|B6P of appendix 
IB] We remind the reader that the choice of initial data 
for Afj_ and i?^^ constitutes a choice of first-order gauge 
freedom that manifests only at second (and higher) or- 
der. In particular, A^ represents the perturbed position 
of the particle, and in this sense the second-order singu- 
lar field — and hence effective source — "knows" about the 
first-order deviation from geodesic motion. 

Following the same logic as in the first-order case, one 
should compute j^ to 0{r) and then choose an arbitrary 
extension, j , to the entire manifold (minus 7). We then 
introduce a regular field j^ by 

C. Second Order Effective Source 

Equation (|45|) gives the general second-order metric 
perturbation in series in r in a particular gauge (the "P 
gauge"). In a smoothly related gauge, the second-order 
metric perturbation is given by (see equations (jASp and 


J = ] 


2 (A-^es 



3 - ] 


and plug in to the second-order Einstein equation ([35 


G^'-^lf] = -G(i)[j^] - G^''\h] (for r > 0). (89) 

While each term on the right-hand-side of equation (1891) 
blows up at r = 0, by construction the sum is 0(1)1^^1 

Since a prescription for computing 5£ij, SBij, and S, has 
now been given, the first two terms on the right-hand- 
side may be considered "known" . Since the remaining 

^^ To see this explicitly, note that equation II89I I holds to 0{r ^) if 
the hats are removed: G'-'^'^U'^ + j^ + 0{r^)] = -G'-^'>[h]. It then 
follows that GW[jSj _|_ G(2)[h] = -G^yR] + o(l) = 0(1). 


Thus the right-hand-side is in fact bounded (but not nec- 
essarily continuous) at r = 0. We may nevertheless drop 
the requirement that r > by interpreting ([89]) in a 

Sobalev (or distributional 1) sense. We therefore write 


G(i)[j«] = 5(2) 




where the effective source is bounded but potentially dis- 

As at first order, one may determine the perturbation 
j in a P-smooth gauge by solving equation (PH)) with ini- 
tial, boundary, and/or gauge conditions such that j^ is with 
sufficiently regular (in this case C^), and it appears that 
the Lorenz condition on the regular field, V^J^ = 0, 
(where capitalization denotes trace-reversal) should be 
an appropriate gauge choiceQ 


We have now developed a prescription for comput- 
ing the global metric perturbation in P-smooth gauges, 
where (by definition) the metric may be written 



J = f - C^hP + - [C^C^g - C^g) , (93) 

for smooth ^ and S. Since the motion is defined to vanish 
in P gauge, the motion in the above P-smooth gauge is 
given by (see equation (lASp ) 


Recall that we previously notated ^^|-y by A^ (see equa- 
tion ([70]) ). Thus our A^ in fact gives the first-order mo- 
tion, and the analysis of section lTVBI has in fact produced 
the first-order equation of motion in equations (I5n|) - (I5^ . 
From the point of view of the systematic calculations per- 
formed there, it comes as some surprise that the form of 
perturbed geodesic equation emerges. We now use a sim- 
ple (and trivial) argument to show why the form of the 

^^ To give a distributional interpretation we promote the entire 
right-hand-side of I I89I I to a distribution. We give no distribu- 
tional interpretation to each term separately. 

^"^ Unlike at first order, however, it is not possible to impose the 
Lorenz condition on the full metric perturbation j by our effective 
source method, since our second-order singular field violates the 
Lorenz condition by a singular amount. (The analog of equation 
II65I I would then contain a singular source term.) More discussion 
of this gauge condition can be found in the text below equation 

perturbed geodesic equation must in fact occur. This ar- 
gument also derives the second-order equation of motion 
in terms of a second-perturbed geodesic form. 

The argument proceeds as follows. In the P gauge, 
the description of motion is geodesic in the background 
metric g (since the perturbed motions vanish). To deter- 
mine the description in smoothly related gauges, use the 
P gauge to promote the background metric to a finite-A 
tensor, g^^iX) = g^iJ , where this equation holds only in 
the P gauge. Within the class of P-smooth gauges, one 
now has the invariant description of motion that Z'^{X) 
is geodesic in g^^{X) + 0{\^). Perturbatively, we have 


A/i^« + X^f" 





and it follows that the first and second perturbed po- 
sitions Z^^) and Z'^) must satisfy the first and second 
perturbed geodesic equation in first and second pertur- 
bations h^'^ and j^'^ . At first order, we have already 
found that the motion is given by the perturbed geodesic 
equation in our regular field /i^. But from (e.g.) equation 
(|551) we have 


h" + h""- + 0{r^) 

= h"'' + 0{r'), 


where the second line follows from the fact that h^ = 
0{r^). Since the perturbed geodesic equation (equation 
(jA9P ) includes only first spatial derivatives of the per- 
turbation, equation (jlOOl) shows that the statement of 
geodesic motion in h^ is equivalent to the statement of 
geodesic motion in h^'^ . This "explains" the appearance 
of the geodesic form in equation ([75)) , and suggests that 
the motion is more naturally regarded as geodesic in h^'^ 
(which happens to coincide with our choice of h^ to the 
relevant order). This viewpoint has fundamental appeal 
in that the motion, which is pure gauge, is given in terms 
of a pure gauge metric perturbation. 

To determine the second-order equation of motion, we 
could similarly proceed to directly "solve" equation ([M)) 
for S, as we did in section [TVBI to solve equation ([M)) for 
^ (though our goal there was the formulation of a second- 
order effective source). However, we may avoid this task 
by appealing to the above argument, which shows that 
the second-perturbed description of motion is the second- 
perturbed geodesic equation in perturbations h^'^ and 
jBG^ From equations ([M|), ([57]) and ([MD, we have for 
our particular choice of j^ that 


Since the second-perturbed geodesic equation (equation 
()A10|) ') contains only first spatial derivatives of the sec- 
ond perturbation, we may equally well use j^ instead of 


jBG jjj^ determining the motion. However, the second- 
perturbed geodesic equation also contains a term involv- 
ing the second spatial derivative of the first perturbation, 
and for this term the difference between h^ and h^'~^ is 
relevant, since these quantities agree only to 0{r) (see 
equation (j99| ). To solve for the second-order motion one 
must first determine h^^ . This may be accomplished by 
subtracting h^ from h^ (see equation ([M]) ). where h^ 
may be determined from equation (l42t with (l72|) - d73t . 
The motion is then given by solving equation (|A10|) with 

^(1) ^ h^G and g(2) ^ jBG^ 

Note, however, that the term relevant for the differ- 
ence between h^ and h^^ is simply {l/2)Z'^^'>^djd^h^„^ 

(appearing in the second line of the expression for Z ■ 
in equation (lAlOp ). and so it is in fact only necessary 
to consider the 00 component of h^ . In particular, we 

hio"" = C - Co = h^o + 6£^3^^ + 0{r^), (102) 

where 5£ij is given by equation (|72|) . 

Although we view the interpretation of geodesic mo- 
tion in the BG fields as having more fundamental sta- 
tus, it is the regular fields that will arise in practice, 
and we now explicitly present the final equations of mo- 
tion in terms of the regular fields. Using equations (j99| . 

dTUTj) and (Uni to relate {h^°,j^°} to {h",f}, the 
final equations of motion (equations (jA9p - (jA10[) with 
5(1) -^ h^'^ and .g(^) -^ j^^) become 

7(1) - fi h^ 




- hlz('> + z^''^-'d,{-d,h^, + \d.h^^ 


where 5Eij is given by equation ([7^ . In this form, the 
second-order equation of motion is seen to be geodesic 
in the regular fields, up to a correction term (the term 
proportional to 5£ij) that accounts for the fact that the 
motion is in fact only geodesic in the BG fields Fn 


The content of this paper has been the identification 
of a class of gauges for which the motion may be sensibly 
defined and the development of a prescription for com- 
puting the metric and motion in such gauges. This class 
was chosen by constructing a particular mass-centered 
gauge (called P gauge) and considering the class of all 
gauges related by smooth first and second order gauge 
vectors (P-smooth gauges). In constructing the P gauge 
many particular choices were made, and the reader may 
wonder the effect of making different choices, leading to a 
P' gauge and possibly distinct class of P'-smooth gauges. 

Suppose that the content of section IIIII were repeated, 
except that a different mass-centered gauge, called P' 
gauge, were chosen. For concreteness, the reader may 
imagine that we chose Cartesian Schwarzschild coordi- 
nates rather than Cartesian isotropic coordinates for the 
near-zone background metric, and did not make the ad- 
ditional first-order gauge transformation, equations ([28])- 
(P^)) . This would produce a P' gauge that is related to 
our P gauge by a first-order far-zone gauge vector of the 
form V^ — n^ + 0{r) (as well as by analogous second and 
third order gauge vectors) , which modifies the metric per- 
turbation by a singular amount (changing the structure 
of A^(-^) from isotropic-type to Schwarzschild- type) . Af- 
ter identifying an appropriate singular field (one option 
would be transforming the old singular field by V^), one 
could develop an effective source method to determine 
the metric perturbations in P'-smooth gauges. Since the 
P' gauge is mass-centered, the analysis of the motion 
will then proceed identically, leading to a prescription 
for determining the motion, {Z''^^\ Z'^^^}, in P'-smooth 

It is clear that the perturbations {h,j} and {h',j'} in 
P-smooth and (respectively) P'-smooth gauges thus con- 
structed will differ by a (possibly singular) gauge trans- 
formation (provided that the initial data differ by a gauge 
transformation), and thus represent the "same physics". 
The reader may further wonder whether {Z^^\ Z^2)| ^^d 
{Z'(-^\ Z'(2)| tj^ng constructed also represent the "same 
physics". However, since the gauge transformation law 
for a curve, equations (|A4p and (|A5p . does not make sense 
in the presence of singular gauge vectors, we have no a 
priori criterion with respect to which to check this type 
of gauge covariance property. Instead, we may view our 
definition of motion as (in principle) providing a gener- 
alized gauge transformation law for the motion that en- 
sures that {h, Z^^\j, Z(2)} and {h', Z'^^\j', Z'^^)} repre- 
sent the "same physics" . For smooth gauge vectors, the 
law trivially agrees with equations (jA4p and (jA5[) . For 
(non-smooth) gauge vectors that link a P-smooth gauge 
to a P'-smooth gauge (for particular known choices of 

Note that S£ij does not represent the perturbed Riemann tensor 

of h^ but rather that of h^ , which is related to that of h^ 
by equations I I67I I- II69II (see also II102II ). This accounts for the 
positive sign in front of the S£ij term. 


P and P'), it should be possible to derive such a law 
by writing the gauge transformation as the composition 
of a smooth transformation with the singular (but mass- 
centered-preserving) transformation that relates P gauge 
to P' gauge; the law is then simply be equations (IA4p 
and (|A5|) using the smooth transformation. For general 
gauges, the situation is clouded by the fact that the al- 
lowed form of the gauge transformation is conjectured 
but not known [30| , and further that (even restricting to 
the conjectured class) the class of gauges smoothly con- 
nected to a mass-centered gauge (i.e., those for which we 
can define the motion) is not known at second order. In 
the appendix of paper I (see also [30|]) we obtained some 
results at first orderjll we have not obtained analogous 
results at second order, where the situation is far more 
complicated. However, while such results would certainly 
be of some theoretical interest, we see no practical draw- 
back to simply working in a particular class of gauges 
(such as our P-smooth class) for which the motion can 
be sensibly defined and computed. 


Our perturbation expansion describes asymptotically 
small departures from a fixed background metric g and 
background worldline 7. This should allow one to in- 
vestigate local-in-time effects, such as second-order cor- 
rections to quantities already investigated at first order, 
including gravitational redshift |3ll| , stability of circular 
orbits [32] , periastron advance |33| , loss of energy and an- 
gular momentum, and "snapshot" waveforms [SJ]. How- 
ever, if the goal is to produce waveforms reflecting an 
entire inspiral, it is clear that our expansion off of a fixed 
background geodesic will eventually produce inaccurate 
results. In order to produce the waveform templates 
needed for gravitational- wave data analysis, therefore, it 
will be necessary to go beyond a perturbation expansion 
off of a fixed background geodesic. 

In principle, it seems clear that one should simply 
"patch together" a sequence of perturbation expansions 
off of a sequence of background geodesies. However, the 
details of implementing such a procedure appear to be 
quite problematic. For example, while it seems clear that 
the new background geodesic should be chosen tangent to 
the old perturbed motion and that initial data for the new 
perturbed motion should be trivial, it is far from obvious 
how to choose the initial data for the new metric pertur- 
bation, which satisfies a different field equation (with a 
different effective source). The whole procedure is fur- 
ther complicated by the choice of gauge: both the metric 

^^ More precisely, we showed that for first-order gauge vectors of 
the form ^'^ = F'^{t,n) + 0{r) for smooth F^*, the first-order 
motion changes by SZ'^^'' = (3/47r)(?i''_Fjn*), where the angle 
brackets denote an average over the sphere. 

perturbation and the position perturbation are gauge- 
dependent, and one would require a way of ensuring that 
the new choices are in the "same" gauge as the old. It 
is easy to see how carelessness in this matter can lead 
to unphysical results: Since the choice of the next back- 
ground geodesic depends on the choice of gauge, a naive 
proposal wherein one simply chooses "no incoming radia- 
tion" with some gauge choice at each step would produce 
a final waveform that depends on the gauge choices made. 

These difficulties are well-known, and a number of 
approaches have been developed. In paper I, we used 
the Hadamard form (e.g., [23|) of the Lorenz gauge re- 
tarded metric perturbation together with a point particle 
description to argue that the "MiSaTaQuWa equation" 
[35| — a modified linearized Einstein equation sourced by 
a point particle on a non-geodesic trajectory determined 
by an integrodifferential equation — should provide an ac- 
curate long-term description. Unfortunately this argu- 
ment has no natural generalization as it stands, since 
we have given no Hadamard or point particle descrip- 
tion at second order. A derivation of Pound [2^ directly 
obtains the MiSaTaQuWa equation by expanding in the 
acceleration of an unspecified worldline, and is a promis- 
ing route toward obtaining a second-order generalization. 
However, both of these approaches depend on the Lorenz 
gauge in an essential way (through its "relaxation"), 
and it has not been investigated whether analogous pre- 
scriptions based on relaxing alternative gauge conditions 
would produce the same physical waveform. Neverthe- 
less, it seems likely that MiSaTaQuWa equation provides 
a reliable — if computationally challenging — prescription 
for first-order- accurate long-term evolution. 

An alternative, "adiabatic" approach to long-term evo- 
lution has been pursued by Mino [36| and Hinderer and 
Flanagan [37|. Here, one considers bound orbits of a 
Kerr black hole and assumes adiabaticity in the sense 
that the radiation reaction timescale is much longer than 
the orbital timescale. This assumption allows one to use 
self-force results (such as would be provided at second- 
order by applying the prescription of this paper) to de- 
termine an adiabatic evolution of the orbital parameters 
of the background geodesic. As in the non-adiabatic ap- 
proaches, above, the gauge dependence of the prescrip- 
tion has not yet been carefully analyzed. However, it 
has been suggested (in both the Mino and the Hinderer- 
Flanagan approaches) that simple conditions reflecting 
"no secular growth over short timescales" should lead 
to a gauge-invariant waveform. If the relevant condition 
on the gauge can be precisely identified, it should be 
straightforward to choose such a gauge within our for- 
malism, since we allow a wide class of smoothly-related 
gauges. In particular, the Lorenz condition applied to 
the regular field is a locally defined gauge condition and 
therefore should lead to perturbations that do not ex- 
hibit secular growtho Thus the combination of our re- 

By contrast, the approach of [3 requires one to work in a mass- 


suits with the work of [36|, l37| appears to be a promising 
approach to producing second-order-accurate waveform 
templates for gravitational-wave astronomy of extreme 
mass-ratio inspirals. 


We conclude by summarizing the prescription for com- 
puting the first and second-order motion and metric. 
First, choose a vacuum background spacetime g, such as 
Schwarzschild or Kerr. Next choose a timelike geodesic, 
7, of that spacetime (representing the lowest-order mo- 
tion of the body) , and choose and a point 70 at which the 
perturbed motion is taken to be coincident. Determine 
a coordinate transformation between a global coodinate 
system for g and a local RWZ coordinate system about 
the geodesic, equation (fT7|). which in particular deter- 
mines STF curvature tensors {Eij,Bij,£ijk,Bijk\- Now 
compute h^ to 0{r^) in terms of these STF tensors from 
equations (jBll) - (IB3p and choose an arbitrary extension, 
h^ . Then compute the effective source, equation (|62|) . 
and solve equation (1611) for /i^, imposing a convenient 
gauge condition on h^ such that h^ is C^ . The first- 
order metric perturbation h is then given in a P-smooth 
gauge hy h — h^ + h^ . If one is stopping at first or- 
der, one may now determine the first-order motion Z^^' 
by integrating equation (|103p with trivial initial data at 

70EI ^ ^ 

If one is proceeding to second order, one should in- 
stead integrate equations (l80l) - ([83|) for A^ and Bf^^. (The 
integration for A'' is redundant with an integration for 
^■(i)^ ^ ^M via equation ^U^.) The initial data for 

Af^ = Z p a-iid BiQ = Ai = Z \ should be trivial (con- 
sistent with the interpretation of the particle being ini- 
tially coincident with the background worldline), while 
the initial data for Buj-\ is arbitrary (trivial being one 
allowed choice). Next compute j^ to 0(r) from equa- 
tions (IB4[) - (jB6p and choose an arbitrary extension, j^ . 
Then compute the second-order effective source, equa- 
tion (|M|). and solve equation (|M| for j^ , imposing a 
convenient gauge condition on j^ (such as the Lorenz 
condition) such that j^ is C^ . Finally, the second-order 
motion Z^ "> is given by integrating equation (jl04p with 
trivial initial data at 70, and the second order metric 
perturbation j is given by j = j -{- j^. The first-order 

motion was previously calculated as Z ^ — A^, and the 
first-order metric perturbation was previously calculated 
as ft = h^ + h^ . Second-order observables may be con- 
structed from the combination {h, Z'^^\j, Z^^^^}. 

centered gauge at first order, in wiiicii case the metric perturba- 
tion should exhibit secular growth. 

If one is stopping at first order, one only requires that h be C^ 
instead of C^ . Correspondingly, one may choose to compute h^ 
only to 0{r) when constructing the effective source. 


Support for this research was provided by NSF grant 
PHY08-54807 to the University of Chicago and by NASA 
through the Einstein Fellowship Program, grant PFl- 
120082. Some of the key ideas for this work arose in con- 
versation with Robert Wald. The author also acknowl- 
edges Abraham Harte for helpful comments. 

Appendix A: The Perturbed Geodesic Equation 

In this appendix we consider the perturbation expan- 
sion of a metric and worldline and discuss the gauge free- 
dom through second order. We then suppose that the 
worldline is geodesic and derive the perturbative descrip- 
tion through second order, expressing the results in the 
RWZ coordinate system (equation (fTT)) ) used throughout 
the paper. 

Fix a coordinate system x^ on a manifold M. Consider 
a smooth one-parameter- family of metrics ^^^(A) along 
with a smooth one-parameter-family of timelike curves 



Taylor expanding the metric components 
and worldline coordinate position Z'^{X; r) gives 

= 3l5(^) + Ag(l'(^) + A^ff^lH^) + OiX') 


= Z(")^(r) + AZ(i)'^(t) + A2z(2)A'(r) + 0{X^), 


where we have defined metric perturbations g/," = 
{l/nl)d\g^^\\=o and coordinate position peturbations 
^(")m ^ (l/n!)9A^''|A=o in the usual way These quan- 
tities depend on the choice of coordinates x^. Under a 
change x'^{X; x"), the coordinate position Z^{X; r) trans- 
forms by Z'^{X;t) = x'^{X; z^{X;t)), while the metric 
components transform via the tensor transformation law. 
Restricting to coordinate transformations that reduce to 
the identity at A = gives the "gauge freedom" within 
perturbation theory. Following [38| we write the coordi- 
nate transformation as 

^X^^ + lx^E'^+^d.^n + OiX''), (A3) 


SO that the smooth vector fields ^ and S are the first 
and second-order generators of the diffeomorphism cor- 
responding to the coordinate transformation. The trans- 
formation laws for Z'^{X) and g^j,,j{X) now give 







0'(2) = „(2) 









where primed perturbations are defined via Taylor expan- 
sion (in A) of components in the primed coordinate sys- 
tem. Under a change of coordinates for the background 
spacetime (i.e., a A- independent change of coordinates 
for M), the metric perturbations transform as tensors 
on M, while the coordinate position perturbations trans- 
form as vectors on Z^^'. Thus if we work exclusively 
within perturbation theory, we may remove A from the 
description and view the perturbations as tensor fields 
on the background spacetime that obey additional gauge 
transformation laws. 

Now suppose that each curve Z'^{X;t) satisfies the 
afhnely-parameterized geodesic equation in (7^1, (A; x). 


Z'^Z^ = 0, 



where an overdot denotes an ordinary derivative with 
respect to t (at fixed A). In perturbing equation (jA8l) it 
is convenient to choose RWZ coordinates (equation (fT7|) ) 
for the background metric, where Z^°)* = and Z^^^^^ = 
t = T. Plugging in the expansions (|A1[) and (jA2l) and 
collecting powers of A yields 

^1) _ 


7(1) _ f) „(1) , ^ft J^) c 7(1)3 




o5oo 9ou 

+ Z('^'^d,g^;i + ]-Z'^'^-d,dJ^ 

z^'\ = -dog'o 


where all quantities are evaluated on the background 
worldline 7 and the background metric (equal to 77^^ on 
7) is used to raise and lower indices. 

We have left equations (IA9|) and (jA10[) in coordinate 
form, which is sufficient for practical purposes, since ap- 
plying the prescription of this paper will require con- 
structing RWZ coordinates in any case. However, it 
is straightforward in principle to convert these expres- 
sions into covariant language using the formulae for the 
STF curvature tensors, equations ([M)) - (P7)) . as well as 
the fact that the background four-velocity is given by 

Z(^)°' — (1,0,0,0). For use in comparing to previous 
work, we give the covariant version of (jA9p . 

Z^\ = 




1^ IJ-9aP 

n(0) 77(1) 


Equation (jAlip differs from other equations sometimes 
referred to as the "perturbed geodesic equation" in two 
ways. The first difference is that we have no projec- 
tion orthogonal to the background worldline. This cor- 
responds to our choice of an affine parameter in the per- 
turbed spacetime, equation (jASp . as opposed to a pa- 
rameter such that the perturbed tangent vector Z^^^^^ is 
normalized in the background metric o The second dif- 
ference is that a "geodesic deviation" term (involving the 
Riemann tensor of the background) appears in our equa- 
tion. If one assumes a A-independent metric family, the 
definitions and calculations of this appendix reproduce 
standard derivations of the geodesic deviation equation 
(and provide a second-order generalization). A version 
of the perturbed geodesic equation without the geodesic 
deviation term would have to refer to a definition of the 
motion perturbation and/or metric peturbation that dif- 
fers from our straightforward Taylor expansion. 

Appendix B: Expressions for the Singular fields 

Here we display expressions for the first and second or- 
der singular fields in the RWZ coordinates for the back- 
ground metric (equation (|17p ). For convenience in dis- 
playing the results, we have let M ^' 1 (corresponding 
to a choice of units adapted to the small body), so that 
explicit factors of M do not appear. Factors of M may 
be restored on dimensional grounds, and explicit instruc- 
tions are given below. 

The first order singular field h is given by equation 
dH]) with equations ©, dHll-dSl) and ^-^. Collect- 
ing those equations together yields 

h^Q = - + 2r£,jn'n^ + ^r^S.jkn'n^n^ + 0{r^) (Bl) 



+ n^Sij ^ n^n^n^Sjk) + 0{r^) (B2) 

/ifj = -5.^ - 2r {2£,j -t- 5.,,£kin''n') 

f2 l^k 

-,k^l K? m 

- 6?i'^fyfc - 2,5,,£fcz™n'=n'n'") + 0{r^). (B3) 

^^ Note that in paper I we eflFectively used the latter parameteri- 
zation by demanding that our deviation vector be orthogonal to 
the background worldline. 


In units where M ^ 1, an explicit factor of M would 
multiply the entire right-hand-sides of the expressions 
given above. The second order singular field is given by 

f = f 


equations (gH), ©, (|46l)-(l5T]) and ^-^, ^ is given by 
equations ([Ml) - ([55|) . and h^ is given by equations ([50]) . 

C^h^ (equation (|87l) ). where j^ is given by gives 

dSU-dlOl) and dTSJ-dlSl). Computing j 


■^00 - ^2 

— tT" 

1 - A'm 

2/i^o + hf^n'ni] +\\ - ^d'h^^n^ - dkhf^n'n^n'' + 2A'n^£^j + n'n'£^J + AA'n^n^n''£Jk 



hd^hgn.nj - lOdidjh^^n'n^ + 2d^dahf^n'n' - 2dkdih^^n'n^n^n^ - 40A'e,kin^n''B/ + AA'eikin^n''B, 

Ah^o^'n^Sij - IQhf^n'n'n^n^Ski + n'n^n''Sijk + ^A'nin'n^n^Sjki] + 0{r^) 

JiO ^ 




1 r 


I2d,h^,n^ + Udohfy + iA^e.kin'^B,' + 8Aie,kin''B/ - e.kin^n'^B/ ~ &A^e,,in''Bk' 

+ 4A^eam.njn''n'^Bk 

k^l r? m 

1 r 


- 6e^"'B[u]n^n''Bj^ - Se^'^h^yn'^B,^ - Se^'^^ S'"' hjln' n^B^p + 6e/"'Byi]n^n''Bkm + ie'rhfiU^n'^Bkm 
hf^n^n''n'n"'BiP + 8A^ejk"'n''n'Bum + 8A^ e.k^n'^n'Bji^ - 8A^ e,,"'n''n'Bkhn - 2e,,"'n^n''n'Bkir 

■ impif'j 

- 76n^£ij + 10A^njn''£ik - 2A> nin''Ejk + 2A,n^n''£jk - I7n,n^n''£jk - 2A^ 


+ 0{r^) 





u 2r2 

3(5„ -f 4A'=(5„nfc 


' 2afe/i,^n^ - -,5,,9„/if;n'=n'n" ~ 8£,j + 4A'=nfc£,y + AH,,n'£ki 

+ 8dod jh^f.-run'' + Sdod^hQ^jUu -f SdodihQ^.Ujv}' - SdjOkh^Qriin'' - Sdodohfj^Uin'' - Sdidkh^QjijU^ 

- Sdodohf^UjU^ - 726ijdod^h^knkni + 48dod^h^^n,njnkni + iUijdkdih^^n^n^ + 36dijdodoh^in''n^ 

+ 169fc5i/i^n'=n' + 8djdihgn''n' + ^d^dikf^n'^n^ - 89,9j/if,n'^n' - 249^ 9, /i^gnin^n'^n' - 2Adodoh^ininjn''n^ 

- 46,jdmdph^in''n'n"'nP - IGA'^e.kiB/ - l%Ahumn,n'B,"' -f KSAh.imUkn'Br + We.i^h^^n'^n'Br 

- mA'^e^kiB/ - wA'^Eki^n.n'Bj"' + wA'^eumUkn'Sj"' + 16e,,„/i^,n'=n'ej'" + wA'^e.k^n.n'Bi"' 
+ 16 A'' e^kmUju' Br + lUA''ekmpS^Jn'n'''BlP ~ A'^ekmpU.n.n'n^BiP + YlA'^t.juBl - 32e,Mn'=fi/ 

- ^A^tumUjn'Br + 2QA^t,i^nun'Br + 12^'=e,fe,S/ - 32e,fe,n''^6/ - ^A^'^umn^n'Bf' -f 20^'=e,,„nfen'^j" 

- 4Aje,/„n''n'6fe" - 4Aiej7mn''n'^fe'" + ihtjimn^n^n^BiT + 35ei,„njn'=?i'^fe'" + 4A'=ejfe„ 71^71'^/™ 
+ 4v4'=e,fe„n,n'^r -I- 48v4'=efc„p%n'n"^,P - ■i2A^,n,n'n^Bi^ + ^A^t.^pn^nun'TTBi^ 

+ ^A^t^^pn^nkU^n'^Bi^ - 16/ito^u " 64/if,n'=7i'f,y - 96S[jfc]£:,'= - 48/ifj;£:,'= + \6h^Qn.^n''£,k + IG/if^n'^n'^:,, 

- 96B[ifc]£/ - A8hg£/ + lQh^Qn,n^£jk + l&hf^n^ n^ £ ji ~ 128^h^Qn^n^£ki - \6hf^^n^£ki + l^h^^n^n^n^n^Eki 
-h 96(5,jS[fe™]n'=n'£:," -h A8Sijh^„,n''n^£i"' - 96B[fe„]n,:njn'=n'£:,™ - 48/if„,ninjn'=n'£:,™ 

-f lQ5^jh'^in''n^rr''nF£rap - 96Ao'^»j + 48Ao(5ijn'=n'4/ - 48^ow^?^J«''"'4^ + 48A''£:ijfc - 96n'=fyfe + 48A'=7ifc?i'£y7 

- 48^'=%n'n™£Mm + 48^'=n,n,n'n"£:fe,,„ + 6,5,jn'=n'n"ffe,„ - l&AH,jnkn'n^nP£i„J + ©(r^) (B6) 

In units where M ^ 1, an explicit factor of M would 
multiply the terms involving /i^, A and -B, while an ex- 
plicit factor of M^ would multiply the remaining terms. 
In writing the above results we have used equations (j77p 
and (|8T|) to express B^^y entirely in terms of /i^^, A^ and 
Buj^ . Thus the second-order singular field depends on 

the quantities {A^ , i?[y j , h^^ , £ij , Bij , f^ fc , Bijk } ■ We re- 
mind the reader that A^^ represents the first-order mo- 
tion, A^ = Z'^^^'', so that the second-order singular field 
may be viewed as depending on the background met- 
ric (through {£ij, Bij, £ijk, Bijk}), the first-order metric 
(through the regular field h^ = h — h^), as well as choices 
of initial data for the first-order motion, Z'^'^^^, and rela- 


tive spatial coordinate alignment, B^^j] . 

Equations (|Bip - (|B6p for the first and second-order sin- 
gular fields are the main computational results of this 
paper. We have performed and verified the calculations 

leading to these expressions using the tensor analysis 
package xTensor [26| for the software package Mathe- 
matica [271. 

[1] L. Barack, Class. Quant. Grav. 26 213001 (2009) 
[2] E. Rosenthal, Phys. Rev. D 73 044034 (2006) 
E. Rosenthal, Phys. Rev. D 74 084018 (2006) 
[3] S. Gralla and R. Wald, Class. Quantum Crav. 25 205009 

[4] S. Gralla, Phys. Rev. D 84 084050 (2011) 
[5] L. Barack and D. Golbourn, Phys. Rev. D 76 044020 

[6] I. Vega and S. Detweiler, Phys. Rev. D 77 084008 (2008) 
[7] I. Vega, B. Wardell and P. Diener, Class. Quantum Grav. 

28 134010 (2011) 
[8] S. Detweiler, Phys. Rev. D 85 044048 (2012) 
[9] R. Geroch and J. Traschen Phys. Rev. D 36 1017 (1987) 
[10] A. Pound, arXiv: 120 1.5089 (2012) 
[11] R. Wald 1984 General Relativity (Chicago, IL: University 

of Chicago Press) 
[12] R. Geroch, J. Math. Phys. 11 2580 (1970) 
R. O. Hansen, J. Math. Phys. 15 46 (1974) 
R. Beig and W. Simon, Gen. Rel. Grav. 12 1003 (1980) 
P. Kundu, J. Math. Phys., 22, 1236 (1981) 
[13] T. Regge and J. Wheeler, Phys. Rev. 108 10631069 

[14] F. Zerilli, Phys. Rev. D 2 2141 (1970) 
[15] K. Thorne, Rev. Mod. Phys. 52 299339 (1980) 
[16] L. Blanchet and T. Damour, Phil. Trans. R. Soc. Lond. 

A 320 379-430 (1986) 
[17] K Martel, Phys. Rev. D 69 044025 (2004) 
[18] K. Martel and E. Poisson, Phys. Rev. D 71 104003 (2005) 
[19] S. Hopper and C. Evans, Phys. Rev. D 82 084010 (2010) 
[20] S. Field, J. Hesthaven and S. Lau, Phys. Rev. D 81 

124030 (2010) 
[21] S. Detweiler, Phys. Rev. Lett. 86 1931-1934 (2001) 
[22] S. Detweiler, Class. Quant. Grav. 22 S681-S716 (2005) 
[23] E. Poisson, A. Pound and I. Vega, Living Rev. Relativity 
14 (2011) 

[24] K. Thorne and J. Hartle Phys. Rev. D 31 1815 (1985) 

X. Zhang Phys. Rev. D 34 991 (1986) 
[25] A. Pound, Phys. Rev. D 81 024023 (2010) 
[26] J. Martn-Garcia, xAct: Efficient Tensor Computer Alge- 

bra, || 
[27] Wolfram Research, Inc., Mathematica, Version 8.0, 

Champaign, IL (2010). 
[28] S. Detweiler and B.F. Whiting, Phys. Rev. D, 67, 024025, 

[29] S. Hawking and G. Ellis 1973 The Large Scale Structure 

of Space- Time (Cambridge: Cambridge Univ. Press) 
[30] S. Gralla and R. Wald, Class. Quant. Grav. 28 177001 

[31] S. Detweiler, Phys. Rev. D 77 124026 (2008) N. Sago, 

L. Barack, S. Detweiler, Phys. Rev. D 78 124024 (2008) 

A. Shah, T. Keidl, J. Friedman, D. Kim, L. Price Phys. 

Rev. D 83 064018 (2011) 
[32] L. Barack and N. Sago, Phys. Rev. Lett. 102 191101 

[33] L. Barack and N. Sago, Phys. Rev. D 81 084021 (2010) 

L. Barack and N. Sago, Phys. Rev. D 83 084023 (2011) 
[34] S. Hughes, S. Drasco, E. Flanagan, J. Franklin Phys. Rev. 

Lett. 94 221101 (2005) 

S. Drasco and S. Hughes, Phys. Rev. D 73 024027 (2006) 
[35] Y. Mino, M. Sasaki, and T. Tanaka, Phys. Rev. D, 55, 

3457-3476, (1997) 

T.C. Quinn and R.M. Wald, Phys. Rev. D, 56, 3381-3394, 

[36] Y. Mino, Phys. Rev. D 67 084027 (2003) 

Y. Mino, Prog. Theor. Phys. 113 733-761 (2005) 

Y. Mino, Prog. Theor. Phys. 115 43-61 (2006) 

Y. Mino, Phys. Rev. D 77 044008 (2008) 
[37] T. Hinderer and E. Flanagan, Phys. Rev. D78 064028 

[38] M. Bruni, S. Matarrese, S. MoUerach and S. Sonego 

Class. Quant. Grav. 14 2585-2606 (1997)