# 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. (N o 1—5 cn ■ 1 ' > 00 cn rn o X 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 assumptions. 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| . II. FORMALISM AND OUTLINE 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, respectively. 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 giO) = v + + a2or'^ + asor^ -f 0{r^) gW = aoir^^ + flu + 0-21?' + asir'^ -f 0{r^) g(2) = ao2r^^ + ai2r^^ + 022 + a32r -t- 0{r^) g(3) = aosr^^ + aur^'^ + a23r^^ + 033 + 0{r), (2) 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-^) ,(2) ,(3) + 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^^) + ... (3) 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. III. LOCAL METRIC IN P-GAUGE 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) (4) (5) (6) i(3) a2or^ + a2ir + 022 + Oif"^) aaor^ + agif^ -I- 032^ + 0{f°)) + i{d2of'^ + a2if -f 022 + 0{f-^)) (7) where A^(") is the 7i*^-order term of the Schwarzschild metric in Cartesian isotropic coordinates. (8) (9) (10) 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) g^ ('^-fMnV[%(- 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 (14) |,=o = e.,fen^S^,„n'n"(-r-3 -|- -Mf^ - -M^f) 9 + n,£kin''n''{~-r 7 -Mf^ -M^f) 9 1] \t-- + e,,kn'5B\n\-f'' + -Mr) + O(f0) (15) + n,n,{-Mf'^ --M'^f)] -,k fol ^m 1 -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. with 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' 2 (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} (18) (19) (20) (21) (22) n (23) 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 i^0jfe/|7 £ijk — V(fci?|o|i|o|j)l7 3 (24) (25) (26) (27) 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 B,. 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 vq 10 r^£,.,n'n-' (28) Vi = r^ {2£ijn^ - ni£jkn^n'') + r^ ( -£ijkn^n'' - -n^£kim.n''n^n'"'' - - i;etjkn''n^&i\ . \ / (29) 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^), (17) M'^^'>r-^+a2ir + a3ir^+0{r^) (30) with (aoo)2i = 2M£ijn'n^ k„l (010)21 = ^Me.^kB in (aij)2i = -2M{5ij£kin''n' + 2£ij) and 2 ("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) (31) (32) (33) (34) (35) (36) (0^)31 = -S£ijn^n^ 2 (0^^)31 - -e^jkn'SB\n^ (0*^)31 = -SijS£kin'^nK (38) (39) (40) 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) ,H (03ir^ + 0(r3). (42) 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) S£i 6B,. I^.\h"] "■QiQj 1 e'\Rl^]Uh"] 'Ojkl I (43) (44) ''' 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^) with (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 (45) (46) (47) (48) (49) (50) (51) 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- tion. 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). (52) 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]) , 9 = r] + + 020^^ + a,or^ + Oir^) hP ^X(i)j.-i + + 021?' + asir'^ + 0{r^) f = X(2)^-2 +0 + 022 + a32r + 0{r^) kP = 7W(3)r-3 + + 0(r-^) •) (53) 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. IV. GLOBAL METRIC IN P-SMOOTH GAUGES 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) Gi%\+' 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 by = h^ + h" -C(^g + 0{r^), (56) (57) 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^ . 8 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'^ (58) (59) Plugging equation ((58|) into the linearized Einstein equa- tion ([M)) gives G(i)[/i«] = _G(i)[/i' s^ (for r > 0). (60) 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), (61) where the G° source S^^' is given throughout the mani- fold by 5(i) = -G(i)[/i^]. (62) 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). (63) where £"^1/ is the Lorenz-gauge linearized Einstein tensor (a well-studied hyperbolic wave operator on g), E^,,[h] = VV^H^, - 2P" /iJ„ P- (64) 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 E^Ah""'] c(i) V(^V"ff,V, (65) 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 -H C^g + 0{r^), (66) 10 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 4'i.J/^ 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 have SSi. 6Br ^2 i [^Ojka'T' \ 7 ^^Rojkl (68) (69) 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 A,- ~ ?/^ 7 (70) B^u = = (v.Cm)!- y (71) we have ^^'= = 2 (- -d,d,h^^ + 25o9(./if)o - d^d^hf^ + ^QO^ij A^ and _B^^ that B, ifii^) 2 V''l7 u"V„A„ = B„„u' v„B^, = i?^,„;3u"A^ - M"r«,[/i^]| (77) (78) (79) 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 1 :h, 00 An Ai — Bio B,o = -£^,A^ - d^h^ ly- ^U'-OJ T o^'i^OO %,]=6,,,S',.A^ + %/ij|o, (80) (81) (82) (83) 2i?fcoe'(,a,)' -t- T^A^e ,(j6,)' + 2Bki,£_ i) ^O'S'i A. 5,-, (72) SB^, = {djdih, Ok dkdih^j) R S'^'BkiB,, + -CS, oo'-'u 2i?H,B,, '^ - Ao4 + TjAfefi,/- 5fc[ji ^ H'T'Ok^j 2e\£/B BkO^ijl^ kl J "feO 1 3 Ak€iii£ k'^ijlX (73) (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'), (74) 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. .(1) " 7Q/3 1 W = n i^ahpj + Vpha^ - V^/la/j) (76) 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 11 6 + 0{r^) 1 \ J.2 SfcOt jllJl kt jiOt 12 3 1 6 -^BjaeumBk" + —AjeumBk"' - -AQtij'^Bkira - i;Bjk£ii - Bij£ki - -A.fjfc; ) x'x'^x^ + Oir"), -Idkdihf^ + ^dAhfk - ^h^jSki + leurnhlljBkA x^x'^x' (84) (85) r 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 J" ] Cfh' (87) 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 m) J = ] CehP 2 (A-^es ^E.g) (86) 3 - ] (88) and plug in to the second-order Einstein equation ([35 get to 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). 12 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 with G(i)[j«] = 5(2) ^(2)^_G(i)[jS]„G'(2)[/i] (90) (91) where the effective source is bounded but potentially dis- continuous. 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 V. MOTION IN P-SMOOTH GAUGES 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 ^5 (92) 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 ) (94) (95) 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 ||51}. 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 5""(A)=ff A/i^« + X^f" 0(A^ f^^\{C^C^g-C^g), (96) (97) (98) 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 BG h" + h""- + 0{r^) = h"'' + 0{r'), (99) (100) 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 (101) 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 13 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 have 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^ and (i)i (103) - hlz('> + z^''^-'d,{-d,h^, + \d.h^^ (104) 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 VI. CHOICE OF P-GAUGE 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 gauges. 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. 14 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. VII. INSPIRAL 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- 15 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. VIII. SUMMARY OF PRESCRIPTION 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. ACKNOWLED GMENTS 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 Z^(A;t). 5M''(A;a;) ffM''(A;a;) Z^iX;r) Taylor expanding the metric components and worldline coordinate position Z'^{X; r) gives = 3l5(^) + Ag(l'(^) + A^ff^lH^) + OiX') (Al) = Z(")^(r) + AZ(i)'^(t) + A2z(2)A'(r) + 0{X^), (A2) 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) r'A" 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 Z(i)^(r) Z'(2)A'(r) ?^lz(")(r) Z^^^^{T)+^^\zmir) Z(0)(t) and ^e^ 0'(2) = „(2) v-^ .giO) ^ch^"}-c. (A4) (A5) (A6) (A7) 16 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'^Z^ = 0, Z(t) (A8) 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) _ (1) 7(1) _ f) „(1) , ^ft J^) c 7(1)3 (A9) and ^^^U-iOna<^n^-.aW^«'' o5oo 9ou + Z('^'^d,g^;i + ]-Z'^'^-d,dJ^ z^'\ = -dog'o (AlO) 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^\ = -^^gfl 1. (1) 1^ IJ-9aP n(0) 77(1) ^(0)q^(0)/J^ (All) 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) So -re + 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. 17 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 J 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 ^^f~C^h^ ■^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 12 4a 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 ^ "•oo -2 /^?. -2A, 1 r 6 (B4) 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 —r 9 - 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^ nin.jn'^n'-£u + 0{r^) (B5) J aS 1 u 2r2 3(5„ -f 4A'=(5„nfc 1 ' 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^tun.pn,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- 18 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 (2008) [4] S. Gralla, Phys. Rev. D 84 084050 (2011) [5] L. Barack and D. Golbourn, Phys. Rev. D 76 044020 (2007) [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 (1957) [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, |http://xact.es| [27] Wolfram Research, Inc., Mathematica, Version 8.0, Champaign, IL (2010). [28] S. Detweiler and B.F. Whiting, Phys. Rev. D, 67, 024025, (2003) [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 (2011) [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 (2009) [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, (1997) [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 (2008) [38] M. Bruni, S. Matarrese, S. MoUerach and S. Sonego Class. Quant. Grav. 14 2585-2606 (1997)