arxiv:gr-qc/ v1 17 Dec 1997

Transcription

1 3rd post-newtonin higher order ADM Hmilton dynmics for two-body point-mss systems rxiv:gr-qc/97075v 7 Dec 997 Piotr Jrnowski Institute of Physics, Bi lystok University Lipow 4, 5-44 Bi lystok, Polnd Gerhrd Schäfer Theoretisch-Physiklisches Institut, Friedrich-Schiller-Universität Mx-Wien-Pltz, Jen, Germny Abstrct The pper presents the conservtive dynmics of two-body point-mss systems up to the third post-newtonin order /c 6. The two-body dynmics is given in terms of higher order ADM Hmilton function which results from third post-newtonin Routh functionl for the totl field-plus-mtter system. The pplied regulriztion procedures, together with mking use of distributionl differentition of homogeneous functions, give unique results for the terms in the Hmilton function prt from the coefficient of the term νp i i r. The result suggests n invlidtion of the binry point-mss model t the third post-newtonin order. PACS numbers: 04.5.Nx, Db, Jd, Lf Introduction nd summry The clcultion of generl reltivistic equtions of motion for compct binry systems, in the recent pst, ws exclusively devoted to the obtention of higher order post-newtonin grvittionl rdition rection contributions. With the works by Iyer nd Will, Blnchet, nd the uthors 3 the rdition rection levels hve been completed up to 3.5 post-newtonin 3.5PN order, i.e. to the order /c 7/ beyond the Newtonin dynmics for the.5pn order see, e.g., the review 4. Quite recently, Gopkumr et l. 5 succeeded in giving first results for the grvittionl rdition rection in compct binry systems t 4.5PN order, pplying blnce equtions between fr-zone fluxes nd ner-zone losses of energy nd ngulr momentum. The conservtive 3PN nd 4PN orders of pproximtion, however, hve not been tckled so fr up to the PN level of pproximtion see, e.g., the review 6 nd the pper 7. It is the im of the present pper to develop the two-body point-mss dynmics up to the 3PN order of pproximtion using the cnonicl formlism of Arnowitt-Deser-Misner ADM 8. There re severl spects which mke such clcultion rther interesting: i the well-known problem of pplicbility of Dirc delt distributions for the source of the generl reltivistic grvittionl field; ii the need of the 3PN dynmics for better undertnding of the innermost stble circulr orbit for binry systems 9, 0, ; iii the importnce of the 3PN dynmics E-mil: pio@lph.fuwb.edu.pl E-mil: gos@tpi.uni-jen.de

2 to control the vibility of the PN filters for gvittionl wve mesurements from inspirlling compct binries. Within our cnonicl pproch, together with pplying Hdmrd s prtie finie technic, some specific nlytic regulriztion formul nd generliztion of it chieved by us the non-generlized formul hs been successfully used by Dmour 3, Dmour nd Schäfer 4, nd Kopeikin 5 t the PN level nd by us 3 t the 3.5PN level, nd the distributionl differentition of homogeneous functions, we were ble to explicitly clculte nd uniquely regulrize ll terms which occur t the 3PN order of pproximtion, but one notice: this is the only mbiguous term up to the 3.5PN order of pproximtion, inclusively. The coefficient of the following term turns out to be finite but mbiguous: νp i i r, or, in non-reduced-vrible form prt from fctor 4: m v i vj c rs i j Gm r m v i vj c rs i j Gm r where m,v i,r s nd m,v i,r s denote the msses, velocities, nd Schwrzschild rdii r s = Gm /c, in isotropic coordintes of the bodies nd, respectively, nd r their reltive distnce. This term describes the qudrupole-qudrupole interction of the kinetic energy tensor of ech body with the grvittionl tidl field of the other body scled to the Schwrzschild rdius of the former one. The bodies, in this interction term, re obviously not tidlly deformed this is expected from Ref. 6 where it is shown tht the tidl deformtion comes in t 5PN only but they seem to hve been ttributed extensions of Schwrzschild rdius size. Those extensions re beyond the pplicbility of the 3PN pproximtion so the obtined result my suggest tht fully consistent 3PN compct binry model needs extended bodies s source for the grvittionl field. Only within the complete theory the bodies Schwrzschild rdii re expected to enter the equtions for the binry system in consistent nd unique mnner. The pper is orgnized s follows. In Section we introduce the point-mss model nd the post-newtonin pproximtion scheme nd we develop the constrint equtions to the 3PN order of pproximtion. To obtin n utonomous Hmilton function for the bodies t 3PN order, dropping the dissiptive.5pn level, we introduce in Section 3 the Routh functionl for the totl field-plus-mtter system nd eliminte the field degrees of freedom for the bodies vribles. The utonomous Hmilton function comes out of higher order. The Section 4 is devoted to the explicit clcultion of the higher order utonomous Hmilton function pplying the regulriztion procedures of the Appendices B. nd B.. In the Section 5 thourough investigtion of the obtined mbiguity is undertken. In the Section 6 we compre our results with limiting expressions known from the literture. The Appendix A presents severl explicit metric coefficients nd some useful formule for inverse Lplcins. The Appendix B. is devoted to the Hdmrd s prtie finie regulriztion. The Appendix B. presents powerful nlytic regulriztion procedure bsed on n nlytic formul derived by Riesz nd generlized by us in this pper. In the Appendix B.3 the nlytic regulriztion procedure devised by Riesz is given. The Appendix B.4 shows the distributionl differentition of homogeneous functions. We use units in which 6πG = c =, where G is the Newtonin grvittionl constnt nd c the velocity of light. We employ the following nottion: x = x i i =,,3 denotes point in the 3-dimensionl Eucliden spce R 3 endowed with stndrd Eucliden metric nd sclr product denoted by dot. Letters nd b re body lbels, so x R 3 denotes the position of the th point mss. We lso define r := x x, r := r, n := r /r ; nd for b, r b := x x b, r b := r b, n b := r b /r b ; stnds here for the length of vector. The liner momentum vector of the th body is denoted by p = p i, nd m denotes its mss prmeter. Indices with round brckets, like in φ, give the order of the object in inverse powers of the,

3 3 velocity of light, in this cse, /c. We bbrevite δ x x by δ. An overdot, like in ẋ, mens the totl time derivtive. The prtil differentition with respect to x i is denoted by i or by comm, i.e., i φ φ,i ; the prtil differentition with respect to x i we denote by i. Throughout this pper we extensively used the computer lgebr system Mthemtic 6. The constrint equtions up to 3PN order We consider mny-body point-mss system which intercts with the grvittionl field ccording to the theory of generl reltivity. For such system the constrint equtions in the cnonicl formlism of ADM red see, e.g., Eqs..8 in 7 g / gr g ij π ij πij π ij = g ij p i p j m / δ, π ij j = g ij p j δ. Here g ij = 4 g ij re the field vribles the prefix 4 denotes four-dimensionl quntity, ll unmrked quntities re understood s three-dimensionl, g ij is the inverse of g ij g ij g jk = δ i k, g := det{g ij }, indictes the covrint derivtive with respect to g ij, R is the curvture sclr formed from the metric g ij ; π ij is the cnonicl conjugte to the field g ij. Sptil indices re rised nd lowered using g ij nd g ij, respectively. We use the following coordinte conditions see, e.g., Eqs..9 in 7 g ij = 8 φ 4 δ ij ij, 3 π ii = 0, 4 where ij is the trnsverse trceless prt of g ij δ ij. The trce-free field momentum π ij cn be split into two prts, longitudinl π ij nd trnsversl π ijtt one: π ij = π ij π ijtt, 5 where π ij cn be expressed in terms of single vector π i s follows: π ij = π i,j π j,i δ ijπ k,k π k,ijk. 6 If both the constrint equtions nd the coordinte conditions 3 4 re stisfied, the Hmiltonin of the theory cn be put into its reduced form, which cn be written s H x,p, ij,πijtt = d 3 x φ x,p, ij,πijtt. 7 The reduced Hmiltonin contins the full informtion for the dynmicl evolution of the cnonicl field nd mtter vribles 8, 8, 9. We expnd the constrint equtions nd in powers of /c where we tke into ccount tht see, e.g., Ref. 0 m O c, φ O c ij O, p O c 3 c 4, π ijtt O, π ij O c 3 c 5., 8

4 4 To clculte the reduced Hmiltonin 7 up to 3PN order we hve to expnd the Hmiltonin constrint eqution up to /c 0. Mking use of Eqs. 3, 4, nd 8, fter long clcultions, we obtin φ = 8 φ 64 φ 5 φ φ4 5 6 φ 5 8 φ 35 p 04 φ3 m φ 45 p 5 φ m p 8 φ 3 m 6 5 p 4 8 m φ 4 p pi p j m m ij 6 ij m δ π 8 φ ij 4 φ π ij π ijtt π ijtt 4 φ φ φ,ij φ φ,i φ,j π ik π jk ij φ ij,k φhtt ij htt ik,j 4 φ,k 6 φ ij,k htt ik,j 7 6 φ ij,k ij O c. 9 We lso need to expnd the momentum constrint equtions up to /c 7. Using Eqs. 3, 4, nd 8, we get π ij,j = 4 φ 5 64 φ p i δ 6 φ φ,j π ij φ,jπ ijtt π jk,k htt ij π jk htt jk,i htt ij,k O c 8. 0 All functions entering the right-hnd sides of Eqs. 9 nd 0 cn be written s sums of terms of different orders in /c. To the orders needed in our clcultions, they red φ = φ φ 4 φ 6 φ 8 O c 9, π ij = π ij 3 πij 5 πij 7 O c 8, ij = 4ij htt 5ij htt 6ij O c 7, 3 π ijtt = π ijtt 5 O c 6. 4 Using Eqs. 4 we extrct from Eq. 9 the Hmiltonin constrint equtions vlid t individul even orders in /c, up to the order /c 8. They re φ = m δ, 5

5 φ 4 φ 6 = = 8 φ p m 64 φ 8 φ φ m δ, 6 p m 8 p m 4 m δ 5 φ 8 = π ij 3 φ,ij 4ij, 7 5 φ3 3 φ φ 4 8 φ φ p m 4 p 3 6 m φ 5 p 6 φ 4 m p i p j m 4ij m δ 8 φ π ij 3 π 3 πij ij 5 πij 3 πijtt φ,iφ,j 4 φ φ,ij φ 4,ij 4ij 4ij 4 4ij,k htt 4ij,k htt 4ik,j φ,ij 6ij. 8 Explicit solutions of the equtions 5, 6, nd 7 for the functions φ, φ 4, nd φ 6, respectively, re shown in Appendix A. The full solution of the eqution 8 for the function φ 8 is not known. We split φ 8 into two prts φ 8 = φ 8 φ 8. 9 The function φ 8 is explicitly clculble nd cn be found in Appendix A. The unknown prt of φ 8 we cll φ 8. It is given by φ 8 = S 8, 0 where S 8 := 8 φ π ij 3 π 3 πij ij 5 πij 3 πijtt φ,iφ,j 4 φ φ,ij 4ij 4ij,k 4 htt 4ij,k htt 4ik,j φ,ij 6ij. Using the Eqs. 4 we extrct from Eq. 0 the momentum constrint equtions vlid up to the order /c 7 : π ij 3,j = p i δ, π ij 5,j = φ π ij 3, 3,j π ij 7,j = Γ ij 7,j, 4

6 6 where Γ ij 7 := φ π ij φ φ 4 π ij 3 φ π ijtt 5 π jk 3 htt 4ik π3 k π3,kl l 4jk,i. 5 Explicit solutions of the equtions nd 3 for the functions π ij 3 nd πij 5, respectively, re given in Appendix A, where one cn lso find the formul for the function π3 i connected with π ij 3 by mens of Eq. 6. The TT-prt of the metric in the leding order fulfils the eqution see, e.g., Eq. 4 of 3 4ij = δttkl ij A 4kl, 6 where A 4ij := p i p j m δ 4 φ,iφ,j, 7 nd where δ TTkl ij is the TT-projection opertor defined by see, e.g., Eq. 4 of 3 δij TTkl := δ il i l δ jk j k δ ik i k δ jl j l δ kl k l δ ij i j. 8 The explicit formul for the function 4ij cn be found in Appendix A. In our clcultions we lso need the explicit formul for tht prt of the metric function which diverges linerly t infinity. So we write 6ij where div 6ij r nd conv 6ij see, e.g., Eq. in 3 6ij = httdiv 6ij conv 6ij, 9 /r s r. The function div 6ij is given by the integrl div 6ij x,t = 8π δttkl ij d 3 x A 4kl x t,t x x = δij TTkl t A 4kl x,t, 30 where the function A 4ij is defined in Eq. 7. We hve clulted the integrl from the righthnd side of Eq. 30. The result is given in Appendix A. 3 The conservtive Hmilton function up to 3PN order In the following we re interested in conservtive 3PN Hmiltonin which depends on body vribles only. To chieve this gol we firstly trnsform our field-plus-mtter Hmilton functionl 7 into Routh functionl which is Hmilton function for the bodies but Lgrnge functionl for the field notice for the following the crucil difference between Hmiltonin nd Lgrngin: the functionl derivtive of the ltter is zero, the one of the former not. The Routh functionl up to 3PN order turns out to be R 3PN x,p, 4ij,ḣTT 4ij = H 3PN d 3 xπ ijtt 5 ḣ TT 4ij, 3

7 7 where H 3PN is the Hmiltonin up to 3PN order. The equtions of motion of the point msses determined by the Routh functionl 3 red We eliminte now in the Routh functionl 3 4ij ṗ = R 3PN x, ẋ = R 3PN p. 3 nd ḣtt 4ij by x, ẋ, p, nd ṗ through solving the field equtions which result from H 3PN see Eq. 6. After tht the Routh functionl 3 becomes higher order mtter Hmilton function denoted by tilde of the vribles x, ẋ, p, nd ṗ : H 3PN x,p,ẋ,ṗ = R 3PN x,p, 4ij x,p,ḣtt 4ij x,p,ẋ,ṗ. 33 The equtions of motion of the bodies determined by the Hmilton function 33 red ṗ = H 3PN d H 3PN, ẋ = H 3PN d H 3PN. 34 x dt ẋ p dt ṗ In these equtions the higher time derivtives my be eliminted by pplying lower order equtions of motion. The elimintion of the higher derivtives in the Hmilton function 33 would result in redefinition of the body vribles see, e.g., Ref.. The reduced Hmiltonin H 3PN of Eq. 3 cn be obtined from Eq. 9. After dropping full divergences including the term π ij π ijtt, which cn be written s π j π,jk k πijtt,i, cf. Eq. 6, it reds H 3PN = { d 3 x 8 φ 64 φ 5 φ φ4 5 6 φ 5 8 φ 35 p 04 φ3 m φ 45 p 5 φ m p 8 φ 3 m 6 5 p 4 8 m φ 4 p pi p j m m ij 6 ij π 8 φ ij 4 φ πij π ijtt π ijtt m δ 4 φ φ φ,ij φ φ,i φ,j π ik π jk ij 4 7 } 3 φ ij,k 6 φ ij,k ik,j. 35 Le us note tht the integrls of the full divergences we hve dropped in Eq. 35 do not contribute to the Hmiltonin becuse they fll off t infinity t lest s /r 4. This cn be inferred from the following symptotic behviour of the functions entering on the right-hnd side of Eq. 9 cf. Ref. 9: φ r, htt ij r, πij r, πijtt r, for r ; 36

8 8 see lso our discussion below Eq. 69. The higher order Hmiltonin H 3PN of Eq. 33 cn be split s follows H 3PN x,p,ẋ,ṗ = H 0 H N x,p H PN x,p H PN x,p H 3PN x,p,ẋ,ṗ, 37 where the Hmiltonins H 0 through H 3PN cn be extrcted from Eqs. 3 nd 35 by mens of the expnsions 4. The Hmiltonins H 0 through H PN re known. We re-clculte them below for completeness. Our im, however, is to clculte the 3PN Hmilton function H 3PN. The clcultion of the Hmiltonins H 0 through H PN cn be performed directly s we explicitly know ll functions needed to perform the integrtions. The Hmiltonins cn be written in the form H 0 = d 3 x m δ, 38 H N = d 3 x 8 φ p m m δ, 39 H PN = H H, 40 H = H = d 3 x d 3 x 64 φ 8 φ φ p m 8 p m 4 m δ, 4 { } π ij 3, 4 H PN = H H H 3, 43 H = H = H 3 = d 3 x 5 φ3 3 φ φ 4 8 φ φ p m 4 p 3 6 m 6 { d 3 x 8 φ d 3 x 5 8 φ 5 p 6 φ 4 m p i p j m 4ij m δ, 44 π ij 3 π 3 πij ij 5 6 φ,iφ,j 4ij } 4ij,k, 45 4 { } 4 φ φ,j 4ij. 46,i The Hmiltonin H 3PN, s extrcted from Eqs. 3 nd 35, reds { H 3PN = d 3 x 4096 φ4 3 5 φ φ 4 64 φ 4 3 φ φ 6 8 φ φ φ φ φ 6 p m

9 φ 9 p 64 φ 4 m φ p 3 m 6 5 p 4 8 m φ p pi p j 4ij 4 4ij m δ 6 8 φ 4 m m π ij 3 4 φ πij 3 πij 5 π ij 5 π 3 πij ij 7 π ijtt 5 4 φ πij 3 πijtt 5 π 3 πjk ik 3 htt 4ij 5 8 φ φ,i φ,j 5 64 φ φ,ij 4 φ φ 4,ij 3 8 φ,iφ 4,j 4 φ,ijφ 4 4ij 7 3 φ 4ij,k 6 φ 4ij,k htt 4ik,j htt 4ij,k htt 6ij,k πijtt 5 ḣ TT 4ij p i p j 3 6ij m δ 6 φ,iφ,j 4 φ φ,ij 6ij }. 47 For the clcultion of the Hmiltonin H 3PN we perform some mnipultions which llow us to do the integrtions in Eq. 47 without explicit knowledge of ll functions entering on the right-hnd side of Eq. 47. There re three such functions: the prt φ 8 of the function φ 8 cf. Eq. 9, the function π ij 7, nd the prt httconv 6ij of the function 6ij. The function φ 8 we eliminte by mens of the identity φ φ 8 = φ S 8 φ,i φ 8 φ φ 8,i. 48,i To eliminte the unknown function π ij 7 we use the following reltion, which cn be proved by mens of the Eqs. 4 nd 5, nd by the id of the Eq. 6: π ij 3 πij 7 = π3 i π3,ik k π ij 7 Γ ij 7,j 3 3 φ π ij 3 φ 4 π ij 3 φ πij 3 πij 5 φ πij 3 πijtt 5 π3 i π3,il l,k πjk 3 htt 4ij π3 i π3,il l π3,j k π3,jkm m 4ij,k. 49 All the terms which depend on the unknown function 6ij we were be ble to write s full divergence see Eq. 60 below. To do this we hve used the following reltion, which cn be derived using the explicit formul 8 for the δij TTkl opertor nd the trceless property of the function 6ij : δ TTkl ij A 4kl A 4ij 6ij = { A 4kk,j A 4jk,k A 4kl,jkl 6ij },i. 50

10 0 In the lst stge we eliminte the field momentum π ijtt 5 by x, ẋ, p, nd ṗ. We use the field eqution for the field momentum in the leding order. It reds see, e.g., Eq. 3 of 3: One cn check tht the following reltion holds Substituting Eq. 5 into Eq. 5 we obtin π ijtt 5 = ḣtt 4ij φ π 3 ij TT. 5 φ π ij TT 3 = φ π ij 3 πij 5. 5 π ijtt 5 = ḣtt 4ij φ πij 3 πij Using Eqs. 9,, 48, 49, 50, nd 53 we rewrite the Hmilton function H 3PN given by Eq. 47 s sum of terms 6 H 3PN = H 3I, 54 where the H 3I re defined s follows H 3 := d 3 x { 4096 φ4 3 5 φ φ 4 64 φ 4 3 φ φ 6 8 φ φ φ φ 4 5 p 6 φ 6 m 45 5 φ 9 p 64 φ 4 m 4 I= 3 8 φ p 3 m 6 5 p m 8 6 φ 4 H 3 := H 33 := H 34 := H 35 := { d 3 x 7 8 φ 4 d 3 x d 3 x p pi p j m m 4ij } 4ij m δ,55 6 π ij } 3 8 φ,iφ 4,j 4ij φ π ij TT 3, 56 { ḣtt } 4ij ḣtt 4ij φ 4 π ij TT 3, 57 { φ π ij 3 π ik 3 π jk 3 htt 4ij π 3 i π3,il l,k πjk 3 htt 4ij π3 i π3,il l π3 k π3,km m,j htt 4ij,k 5 64 φ φ,i φ,j 4ij } 4 φ 4ij,k, 58 d 3 x { 7 64 φ φ,j 4ij φ 4 φ 4,j htt 4ij π j 5 π5,jk k ḣ TT 4ij 3 φ π ij TT 3

11 π j 3 π3,jk k π ji 7 Γ ji 7 } φ 8,i φ 8 φ φ 8,i, 59,i H 36 := d 3 x { 5 6 φ φ,j 6ij htt 4jk,i htt 6jk A 4kk,j A 4jk,k } A 4kl,jkl 6ij. 60,i 4 Results of regulriztion procedures To diminish the number of terms we perform clcultions in the center-of-mss reference frme, so we cn use the reltions p p = 0, ṗ ṗ = 0. 6 We lso use, t the 3PN level, the Newtonin reltions between the coordinte velocities ẋ, ẋ of the bodies nd their reltive velocity v := ẋ ẋ : ẋ = m M v, ẋ = m v, 6 M where M := m m is the totl mss of the system. To shorten the formule we introduce the following reduced vribles r := 6π M x x, r := r, n := r r, p := p µ = p ṗ, q := µ 6πν = ṗ 6πν, 63 where µ := m m M, ν := µ M. 64 We lso introduce the reduced Hmilton function Ĥ := H µ. 65 Ĥ depends on msses of the binry system only through the prmeter ν. From now on the ht will indicte division by the reduced mss µ. To clculte the Hmiltonins H 0 through H PN we proceed s follows. The terms H 0, H N, H, nd H given by Eqs. 38, 39, 4, nd 44, respectively we regulrize by mens of the Hdmrd s procedure described in Appendix B.. The terms H nd H given by Eqs. 4 nd 45, respectively we regulrize using the procedure from Appendix B.. The term H 3 from Eq. 46 is full divergence. The integrnd in H 3 flls off t infinity s /r 4 so H 3 does not contribute to the Hmiltonin. We hve checked tht the regulriztion procedure of Appendix B. pplied to H 3 gives zero. The finl result for the Hmiltonin H 0 H N H PN H PN coincides with tht known in the literture see, e.g., Eq. 3. in 7. The clcultion of the Hmiltonin H 3PN is much more complicted. The term H 3 of Eq. 55 we regulrize using the Hdmrd s procedure described in Appendix B.. After long clcultions we obtin Ĥ 3 = 5 35ν 70ν 35ν 3 p 4 8 { 7 4ν 53ν 6ν 3 p 3 νν n p p 3n p 4 p } 6 r

12 { 08 55ν 84ν } p 46 35ννn p p 45 55ννn p 4 64 r 60 { ν 37ν } p ννn p r 3 6 ν r4. 66 The clcultion of H 3, H 33, nd H 34 defined by Eqs we do by mens of the regulriztion procedure described in Appendix B.. The performnce of the procedure differs considerbly in the cse of H 3 nd H 33 compred to H 34. In H 3 nd H 33 ech term cn be regulrized seprtely. This is not the cse for H 34. For ech term of H 34 tken seprtely the limit 30 of Appendix B. does not exist. By series expnsion with respect to ε one cn see tht for ech term entering H 34 the sum from Eq. 30 contins prt proportionl to /ε similr divergences rise in the dimensionl regulriztion procedure in quntum field theory, see, e.g.,. After collecting ll /ε terms they re of the type 75 only the prts cncel ech other. Furthermore, fter regulriztion, there re terms of the type 75 proportionl to ln r. After collecting ll such terms the logrithms cncel out, too. Unfortuntely, the result of ppliction of the regulriztion procedure to H 34 is not unique. The mbiguity which rises cn be expressed in terms of exctly one unknown number nd will be discussed in detil in the next section. The result of pplying the regulriztion procedure of Appendix B. to Eqs , fter tedious clcultions, is Ĥ 3 = 3 { 4 64 ν 3νp 3νn p p } r { 5 ν 603π π ν p } π 5π 44 ν n p r 3 4 ν r4, 67 Ĥ 33 = { ν 4n pn qp q 5n q p 5n p q n p n q 3p q p q }r {n 8 ν vp q 5p n p n p p q v p qp v 3 n p3 q v n q p n p } n pn v p v { 48 ν 47 0νn pn qn v 5 νn p n q 5 8νn qp 46 5νn pp q 4 νn vp q 47 ν n pq v n qp v 3ν 8n pn vp p v 8n p 3 n vp v n p p v 5p v n p 4 v 5n v p 6n p n v p 5n p 4 n v 4n p q v 4p q v } r { 48 ν 55 7νn p 3 n v 03 5νn p n v 37 35νn pn vp 83 8νn pn vp v

13 3 5 3νn p p v 4 77νn p v 9 9νn v p } 34 9νp p v 37νp v 46 7νp v r { 536 ν 3 97π π ν n pn v 6 36 π 45π 464ν n v 76 3π 676 5π ν π 09π 80ν p v } v r3, 68 Ĥ 34 = { } 9 ν 99 3νp 83 43νn p p 5 303νn p 4 r { 90 ν π π ν p } 7 60π 3 6 5π ν n p r π ν r 4 ω p 3n p ν r The number ω in Eq. 69 is unknown due to the non-uniqueness of the results of the regulriztion procedure, s discussed in the next section. The integrls given by H 35 nd H 36 defined in Eqs. 59 nd 60 re integrls of full divergences. To discuss them let us first observe tht the 3PN Hmiltonin given by Eqs not only describes the system of two point msses but cn be pplied lso to the system of two dusty bodies provided the following substitutions re mde: m δ ρ, p i δ P i, p /m δ P /ρ, etc., where ρ nd P re the mss nd the liner momentum density of the th body, respectively. It is obvious tht for dusty bodies the integrnds in H 35 nd H 36 re loclly integrble, so these terms cn contribute to the Hmiltonin only if the integrnds fll off t infinity slower thn /r 4. Now we proceed s follows. We check whether the integrnds in H 35 nd H 36 fll off t infinity t lest s /r 4. If yes, the terms do not contribute to the Hmiltonin, if no, we must study them in more detil. Becuse the symptotic behviour of ll functions entering H 35 nd H 36 is the sme for both point msses nd dusty bodies, we use the bove prescription lso for the system of point msses, treting them s limiting cse of extended dusty bodies. One cn check tht the whole integrnd in 59 flls off t infinity s /r 4, so formlly H 35 does not contribute to the Hmiltonin. For checking the consistency of our regulriztion procedures we hve integrted out ll terms in H 35 except the lst two ones which contin the unknown functions π ij 7 nd φ 8. To do this we hve used the regulriztion of Appendix B.. We hve obtined zeros for ll terms but the first one. This nonzero result is connected with the mbiguity of the term H 34, s explined in the next section. The integrnd in H 36 does not fll off t infinity fst enough not to possibly contribute to the Hmiltonin. The reson is tht the function 6ij is sum httdiv 6ij conv 6ij, where div 6ij r nd conv 6ij /r s r cf. Eqs. 9 nd 30. The integrnd in H 36 for 6ij = httconv 6ij flls off t infinity s /r 4 but it flls off only s /r for 6ij = httdiv 6ij. Therefore we hve clculted H 36 for 6ij = httdiv 6ij. We hve used the regulriztion procedures from Appendixes B. nd B. the need to use the Hdmrd s procedure of Appendix B. rises becuse 4ij contins some terms with Dirc delts. To clculte 4ij properly we lso hd to employ the rule of differentition of homogeneous functions from the Appendix B.4. The result is zero.

14 4 On the bsis of the bove discussion we put the integrls given by H 35 nd H 36 equl to zero, H 35 = H 36 = 0, 70 nd djust the regulrized expression different from zero to the Hmiltonin H 34. Collecting the Eqs we finlly obtin the utonomous higher order 3PN Hmilton function Ĥ3PN. It reds Ĥ 3PN r,p,v,q = 5 35ν 70ν 35ν 3 p 4 8 { 7 4ν 53ν 6ν 3 p 3 νν n p p 3n p 4 p } 6 r { ν 96ν } p 68 5ννn p p 35 67ννn p 4 48 r { π ν 736ν } p π 4ν νn p r { π ν } r 4 Ĥ33 r,p,v,q ω p 3n p ν r 3, 7 where the number ω is unknown nd the function Ĥ33 r,p,v,q is given by the right-hnd side of Eq Ambiguity The source of mbiguity given by the regulriztion procedure described in Appendix B. cn be explined s follows. Vi integrtion by prts, some terms in H 34 cn be represented in different wys. The regulriztion method pplied to both representtions give different results. As n exmple let us consider the integrl of φ φ,i φ,j 4ij which, by mens of integrtion by prts, cn be written s integrl of φ φ,ij 4ij in the Appendix B.4 nother such term is treted. Appliction of the regulriztion procedure to the difference between these integrls the integrnd is full divergence gives in the reduced vribles d 3 x µ φ φ,j 4ij = 3,i 5 p 3n p ν r 3. 7 Only if the result 7 would be zero, ppliction of the regulriztion procedure to the integrls of φ φ,i φ,j 4ij nd φ φ,ij 4ij would give the sme results. We cn lso obtin the result 7 in different wy. Let us denote the integrnd from the left-hnd side of 7 by F i,i nd let us consider the volume integrl d 3 xf i,i, 73 B0,R\Bx,ε Bx,ε where B x,ε =, is bll of rdius ε round the position x of the th body nd B 0, R is bll of rdius R centered t the origin of the coordinte system. We pply Guss s theorem to the integrl 73 nd then we clculte the limit lim dσ i F i lim dσ i F i lim dσ i F i 74 ε 0 Bx,ε ε 0 Bx,ε R B0,R

15 5 with the norml vectors pointing inwrds the spheres B x,ε nd outwrds the sphere B 0,R. It is esy to check tht the integrl over the sphere B 0,R vnishes in the limit R wheres the integrls over the spheres B x,ε diverge s ε 0, so to clculte them we use the Hdmrd s procedure from Appendix B.. The result coincides with tht given on the right-hnd side of Eq. 7. We hve checked tht for ll these terms entering H 34 for which one cn use different, vi integrtion by prts, representtions the mbiguity is lwys multiple of the quntity p 3n p ν r 3 = νp i i r, 75 or, written in non-reduced form nd with the moment substituted through the velocities prt from fctor 4, m v i vj c rs i j Gm r m v i vj c rs i j Gm r, 76 where r s nd r s denote the Schwrschild rdii r s = Gm /c, in isotropic coordintes of the bodies nd, respectively. This expression indictes tht the rdius of the bodies might come into ply lredy t 3PN s indicted e.g. in Ref. 0. The interction described by the terms 76 is the interction of the Newtonin kinetic energy tensor of ech body with the Newtonin tidl potentil of the other body scled to the respective Schwrzschild rdius. We hve tried to resolve the mbiguity 75 in severl wys. Firstly we hve extrcted the locl non-integrbilities in H 34. Let us denote by F the totl integrnd of the H 34 nd by F s =, the non-integrble prt i.e., of the order /r 3 or higher of the Lurent expnsion of the F round the position x of the th point mss. Then we replce F by the loclly integrble expression F F s F s. We hve checked tht the results of pplying the regulriztion procedure of Appendix B. re the sme for F nd F F s F s so the regulrized vlue of the singulr prt F s F s of the integrnd is zero. This is so becuse by subtrcting the locl non-integrbilities we hve trnsferred the integrnd, initilly loclly non-integrble but integrble t infinity, to n integrnd loclly integrble but non-integrble t infinity. We hve studied two further possibilities to overcome the mbiguity. The first one, bsed on the Riesz s kernel representtion of the Dirc delt distribution, is described in detil in Appendix B.3. The Riesz s kernel regulriztion of the divergence 7 gives zero, i.e. it tkes into ccount lso the contributions coming from the points x. We were yet not ble to compute ll terms in H 34 using the Riesz s kernel regulriztion becuse of serious clcultionl problems, s described in Appendix B.3 only the first nd the moment dependent prt of the fifth term out of six terms in H 34 we could clculte, but with ill-defined /ε nd lnr multipliction fctors for the terms of type 75. In the second cse we hve tried to employ the rule of differentition of homogeneous functions coming from the distribution theory. This is described in Appendix B.4. Here we lso hven t fully succeeded in removing the mbiguity. We conjecture tht the mbiguity hs its origin in the zero extension of the bodies. We strted with point-like bodies but the formlism rected in such mnner tht the Schwrzschild rdii of the bodies got introduced. They, however, re fr beyond the pplicbility of the post- Newtonin pproximtion scheme nd thus, the result turned out to be mbiguous. In the future we shll investigte our conjecture in more detil. We close this section with nother observtion. According to the pplied cnonicl formlism one is llowed to put ij of the initil stte is given by H 3PN x,p, ij nd π ijtt equl to zero t ny freely chosen initil time. The energy = 0,π ijtt = 0. It exctly contins the first

16 6 term of H 34. But this term is ill-defined even with respect to the Riesz s kernel regulriztion procedure Appendix B.3. This is further indiction tht the binry point-mss model is ill-defined t the 3PN level. 6 Comprison with the known results The 3PN Hmiltonin derived in Section 4 we now compre with the results known in the literture. We hve found only two such results: the test body limit nd the sttic prt of the full Hmiltonin. Let us stress tht these comprisons do not fix the mbiguity present in our pper. The Hmiltonin expressed in the reduced vribles defined in the beginning of Section 4 describing test body orbiting round Schwrzschild blck hole reds 0 we restored the explicit dependence on the speed of light c Ĥ test = c c r c 4 p r c r c. 77 Expnding 77 with respect to /c nd tking the /c 6 contribution we obtin the 3PN Hmiltonin for test body: Ĥ3PN test = 5p 4 8 7p 3 7p 6r 6r 5p 8r 3 8r4. 78 To obtin the test-body limit of the 3PN Hmiltonin of our two-body point-mss system we hve to substitute ν = 0 into the right-hnd side of Eq. 7 note tht the mbiguous term nd the term Ĥ33 vnish for ν = 0. The result coincides with 78. The sttic prt of the Hmiltonin is defined by mens of the conditions p = 0, π ijtt = 0, ij = In the pper 3 one cn find the following expression for the sttic prt of the mny-body point-mss 3PN Hmiltonin: 3PN = 8 6π m m b m c m d m e H sttic... b c d e d 3 x...d 3 x 5 δ x x δ x x b δ x 3 x c δ x 4 x d δ x 5 x e r r 3 r 34 r 45 4 r r 3 r 34 r 35 r r 3 r 4 r 5, 80 where r mn := x m x n. We hve clculted H3PN sttic from 80 for two-body system. For this we hve used the Hdmrd s regulriztion procedure of Appendix B.. The result is in reduced vribles Ĥ3PN sttic = 8 ν r4. 8 To clculte the sttic prt of the 3PN Hmiltonin in our pproch, we hve to implement the conditions 79. After tht the quntities H 3 H 36 from Eqs vnish. Only some prt

17 7 of the expression H 3 of Eq. 55 survives. Applying for this prt the regulriztion procedure from Appendix B. one obtins the formul 8. Acknowledgments The uthors grtefully cknowledge discussions with Sergei Kopeikin, Jn Pwlowski, nd Prof. Hns Triebel. P. J. is thnkful to Prof. Gernot Neugebuer for invittions to the Friedrich- Schiller-University nd kind hospitlity. This work ws supported by the Mx-Plnck-Gesellschft Reserch Unit Theory of Grvittion, by the Mx-Plnck-Gesellschft Grnt No TG74, nd by the KBN Grnt No. P303D 0 P. J.. A Some explicit solutions In this ppendix we cite the known pproximte solutions of the constrint equtions nd. The functions φ, φ 4, φ 6 implicitly, π ij 3, nd htt 4ij cn be found e.g. in Refs. 7 nd 4. The functions φ 8, π ij 5, nd httdiv 6ij hve been clculted by us. In most of the Poisson equtions given below there re source terms of the form fxδ, where the function f is singulr t x = x. The Poisson integrls for these terms we hve clculted s follows { } { } fxδ = f reg x δ = f reg x δ = f reg x, 8 4π r where the regulrized vlue f reg of the function f is defined in Eq. of Appendix B.. The Poisson integrls which do not contin Dirc delts we hve computed using the explicit formule for inverse Lplcins listed t the end of this ppendix. The Hmiltonin constrint eqution is explicitly solved up to the /c 6 order. The solutions of Eqs. 5 nd 6 for the functions φ nd φ 4 re known for n-body point-mss system. They red φ = m, 83 4π r φ 4 = 6π b m m b r b r 8π p m r. 84 The solution of Eq. 7 for the function φ 6 is fully known only for two-body point-mss systems. It cn be written s φ 6 = φ 6 φ 6 φ 63, 85 where φ 6 φ 6 := { 64 φ 8 φ φ p m 8 p } m 4 m δ, 86 { } := π ij 3, 87 φ 63 := { } φ,ij 4ij = 8π { m r,ij 4ij }. 88

18 8 Computtion of the Poisson integrls from Eqs. 86 nd 87 gives the results φ 6 = p 3π m 3 r 6π 6π 3 b m m b r b φ 6 = 9 8 6π 3 p r 6π b b r r b p i p j ln r,ij m b p 5 m r b r r b, 89 { 4 p p b b b r r b 8p p b b 3p p b b 8p b p b r r b 4p p b b r }. 90 r b The Poisson integrl needed to clculte φ 63 for two-body point-mss systems reds here b { } m m b 4ij = r,ij 36π r 3r brb 3 3r 4 r3 r b 8r r b r rb 3 3r4 b 8r 3 r b r r b r b r r br b 8r 3 br b 30r r b 60r r b r b 30r b r b 36r r 3 b 36r br 3 b 35r4 b 3 64π m r p 3n p 3π m b p b b 4p b p b b b 4p b p b b b 6 p b b b r3 b r r b r Only prt of the solution of Eq. 8 for the function φ 8 is known, nmely } { p b b r b r 4 p b r r b 8p b r r b. 9 φ 8 = S 8 { } φ 4,ij 4ij 4ij, 9 where S 8 := { 5 φ3 3 φ φ 4 8 φ φ φ p m

19 9 64 φ p m 4 p 3 6 m 6 9 } p i p j m 4ij m δ. 93 We hve computed the Poisson integrls from Eq. 9. They re given by p 3 { S 8 = m p 64π m 5 r 46π m 3 b r 9 } b r rb 6π b b 8 p nb p p 3n b p 4 m m b r b r { m 47 86π 3 r p 6 5n b p 6 } b b r r b r rb 5 m p 6π 3 m b rb 3 r r b b { } φ 4,ij 4ij = { 5m m b 6π 4 rb 3 r b m3 m b r 3 b p { } 4ij 6π m r,ij 6π b r r b }, { } m m b 4ij. 94 r b r,ij The momentum constrint equtions 4 cn be written in the form π ij,j = Si. 95 Mking use of the decomposition 6 it is not difficult to write the solution of Eq. 95 in the form π ij = Θij k Sk, 96 where the opertor Θ ij k is defined s follows Θ ij k := δ ij k δ ik j δ jk i i j k. 97 Using Eq. 96 we hve obtined the solutions of Eqs. nd 3 vlid for n-body pointmss systems. The solution of Eq. reds π ij 3 = { p k δ ij δ ik δ jk 6π r r } r r,ijk. 98,k The function π3 i, connected with the function πij 3 by mens of Eq. 6, is equl to π i 3 = 8π p i r 3π,j,i p j r,ij. 99

20 0 The explicit formul for the function π ij 5 reds π ij 5 = 3 6π m r 3 6π p i n j p jn i n p δ ij n p n i nj { m b δ ij p k r b,k r b 4δ ij p 3 i j p 4p i j j i 4p j i p i j r r b δ ij p i j j i p r r b 4 i j p r r b i j p r r b }. 00 The TT-prt of the metric in the leding order /c 4 is the solution of Eq. 6. It is given by here s b := r r b r b ; be wre tht in Eq. 0 in 3 there is misprint 4ij = { p 5n p δ ij p i p j 4 6π m r } 3n p 5p n i nj 6n p n i p j n j p i 8 6π r r b r 3 b b { m m b 3 n i b s b r b s nj b b s n i n j b 6 b s b rb n i n j b nj n i b 5 r b r b r rb 3 3r 8 n i r s b r s nj b 5 r r 3 b r δ ij}. 0 r b r b r r r b s b r r b The prt of the TT-metric t /c 6 which diverges linerly t infinity is given by Eq. 30. After computing the double Poisson integrl in 30 one finds 6ij = 4 6π div t r { 7n p p δ ij 6p i p j m 7p n p } n i nj 0n p n i p j n j p i 6π t b m m b { 8 r 3 b r 3 r r b 5r δ ij r b

21 i bj j bi i j j i r3 r b 44rb 3 i j r 5 r3 r b 7r b δ ij } r r b 8 i bj i j. 0 r r b Below we list explicit formule for the inverse Lplcins which hve been used throughout this ppendix: r r b = ln s b, 03 = r r r b 36 3r r b rb 3r r b 3r b r b rb r r b r b ln s b, 04 r = r 3r r b rb 3r r b 3r b r b rb r rb rb ln s b, 05 r b 8 6 r r b = r 4 30r r b 6rr b 90r rb 3 63rb 4 30rr 3 b 30rr b r b 90r r br b 90r 3 br b 6r r b 90r r b r b 6r br b 90r r 3 b r r b = r3 r b = 90r b r 3 b 63r 4 b r 4 r 40 rb 3rb 4 rr b 6rbr b 3rb 4 ln s b, 06 63r 4 90r r b 6rr b 30r rb 3 rb 4 90rr 3 b 90rr b r b 30r rbr b 30rbr 3 b 6rr b 90r r b rb 6rbr b 90r rb 3 90r b rb 3 63rb 4 3r 4 0 r r b r4 b 6r r b r b r b 3r4 b ln s b, 07 63r 4 90r 3 00 r b rr b 90r rb 3 63rb 4 50rr 3 b 90r r b r b 90r r br b 50r 3 br b 6r r b 90r r b r b 6r br b 90r r 3 b 90r b r 3 b 63r 4 b 40 3r 4 rr b 3rb 4 6rr b 6rbr b 3rb 4 ln s b. 08 B Regulriztion procedures In this ppendix we describe techniques which cn be used to regulrize integrls which pper in our pper. More detils re given in Ref. 5.

22 B. Hdmrd s prtie finie regulriztion Let f be rel vlued function defined in neighbourhood of point x o R 3, excluding this point. At x o the function f is ssumed to be singulr. We define the fmily of uxiliry complex vlued functions f n lbelled by unit vectors n in the following wy: We expnd f n into Lurent series round ε = 0: f n : C ε f n ε := f x o εn C. 09 f n ε = m= N m nε m. 0 The coefficients m in this expnsion depend on the unit vector n. We define the regulrized vlue of the function f t x o s the coefficient of ε 0 in the expnsion 0 verged over ll directions: f reg x o := dω 0 n. 4π We use formul to compute ll integrls which contin Dirc delt distribution. It mens tht we define d 3 xf x δ x x := f reg x. The procedure described here ws used by one of the uthors cf. Appendix B in 7 in the clcultion of the PN nd.5pn ADM Hmiltonins for mny-body point-mss systems. More detils on the pplictions of the Hdmrd s regulriztion cn be found in 5. Also the reltion of the Hdmrd s procedure to the regulriztion described in Appendix B. nd the rule of differentition of homogeneous functions of Appendix B.4 is discussed in 5. B. Riesz s formul bsed regulriztion The following formul, firstly derived by Riesz see Eqs. 7 nd 0 in Chpter of 6, cn serve s tool to regulrize clss of divergent integrls of r α rβ using the nlytic continution rguments: α3 d 3 xr α r β := π 3/Γ Γ β3 Γ αβ3 reg Γ α Γ β Γ αβ6 r αβ3. 3 In deriving the 3PN Hmiltonin for two-body point-mss system more generl integrls thn those which cn be regulrized by mens of the formul 3 pper. We succeeded in deriving the generliztion of the formul 3 for those integrls. The generlized formul reds d 3 xr α rβ r r r γ := Rα,β,γr αβγ3, 4 where Γ α Γ β Γ α β γ 4 Rα,β,γ := π Γ γ reg I / α, α γ I / β, β γ I / α β 4, α β γ 4. 5

23 3 The function I / in Eq. 5 is defined s follows I / x,y := B / x,y B x,y, 6 where B stnds for the bet function Euler s integrl of the first kind nd B / is the incomplete bet function; it cn be expressed in terms of the Guss hypergeometric function F : B / x,y = x x F y,x;x ;. 7 The regulriztion procedure bsed on the formule 4 nd 5 consists of severl steps. We enumerte them now.. By mens of Eqs. 6 nd 6 we replce p by p, ṗ by ṗ, nd ẋ nd ẋ by v.. We eliminte the unit vector n by the reltion Then we expnd the sclr product n n : n = r r n r r n. 8 n n = r r r r r All integrnds we reduce to integrnds which depend on r nd r only. To do this we use set of substitutions. For the integrnds of the type n p k f r,r, k =,...,6, we use the rules obtined by considering the integrtion in prolte spheroidl coordintes n p f r,r f r,r, 0 n p f r,r b f r,r, n p 3 f r,r n p 4 f r,r n p 5 f r,r n p 6 f r,r 3 3 b f r,r, 4 3 b 3 8 b4 f r,r, b 5 8 b4 f r,r, b 45 8 b b6 f r,r, 5 where := n p r r r r r, 6 b := p n p r r r r r r r r. 7

24 4 After this step the integrnd becomes the liner combintion of the type N I= c I r α I rβ I r r r γ I, 8 where α I, β I, nd γ I re integers; the constnt coefficients c I my depend on r, p, ṗ, v, n p, n ṗ, n v, v p, v ṗ, nd p ṗ. 4. We define the uxiliry function J ε µ,ν α,β,γ := Rα µε,β νε,γr αβγµνε3 b 9 nd for the combintion 8 we clculte the limit lim N ε 0 I= c I J µ,ν ε α I,β I,γ I. 30 It turns out tht for ll terms considered in the pper this limit is of the form where A, B, C do not depend on µ nd ν. A B µ ν C ν µ, 3 5. As the regulrized vlue of the integrl of the sum 8 we tke the number A, i.e. d 3 x N I= B.3 Riesz s kernel bsed regulriztion c I r α I rβ I r r r γ I Three-dimensionl Dirc delt distribution cn be represented s the limit reg := A. 3 δx x = lim α 0 I α x,x, 33 where I α x,x is the Riesz s kernel defined s follows: 3 α Γ I α x,x := x x α α π 3/ α Γ The regulriztion of clss of integrls bsed on the kernel 34 nd its reltion to Hdmrd s procedure of Appendix B. is discussed by Riesz 6. For us, the Riesz s kernel regulriztion consists in using the kernel 34 insted of the Dirc delt in the source terms of the constrint equtions, solving these equtions, nd performing the regulriztion of integrls needed to obtin the 3PN dynmics in the wy described below. For point-mss two-body system we need two fmilies of Riesz s kernels, ech ssocited with the different mss. Let the first fmily be lbelled by the prmeter α, nd the second one by β. After replcing the Dirc delts by Riesz s kernels the right-hnd sides of the constrint equtions depend on α nd β. The first step of the procedure is to solve constrint equtions itertively keeping α nd β s unspecified complex prmeters.

25 5 This is very difficult tsk. Out of ll equtions for the functions φ n we hve fully solved only the simplest one for the function φ. Also only the leding order momentum constrint eqution for π ij 3 we hve solved. The TT-prt of the metric is not known even in the lowest order only the prt depending qudrticlly on the moment hs been clculted. The solutions obtined by us re enough to clculte tht prt of the full divergence 7 discussed in Section 5 which depends qudrticlly on the moment. The regulriztion of the full divergence by mens of the procedure described in Appendix B. gives non-zero result. If one uses the Riesz s kernel generliztion of the function φ nd of the moment dependent prt of the function 4ij the result of regulriztion of the full divergence using the procedure described below is zero. The procedure is similr to tht described in Appendix B.. The first three steps of the regulriztion procedure from Appendix B. pply here. After tht ll integrnds to be regulrized re functions of r, r nd the prmeters α nd β. So ny integrnd cn be written s f r,r,α,β. 35 We need now to use formule nlogous to the generlized Riesz s formul 3 from Appendix B. to perform integrtion of 35 for ny α nd β. Let s denote the result of integrtion by R α,β. Then the following limit is clculted: lim Rα µε,β νε. 36 ε 0 For the divergence 7, to perform the integrtion it is enough to use the formul 3 from Appendix B. for γ = 0. The limit 36 comes out to be zero. For ll terms for which we were ble to pply the bove procedure we obtined the limits 36 of the form given in Eq. 3. However, we were not ble to check the cncelltion of ll divergent terms in H 34 fter regulriztion becuse we were not ble to clculte ll terms using the Riesz s kernel procedure. But gin, the ill-defined terms we hve obtined were of the type 75. B.4 The rule to differentite homogeneous functions. In distribution theory there exists rule to differentite some homogeneous nd loclly nonintegrble functions under the integrl sign 7. We hve studied how useful this rule is in our clcultions. Let f be rel vlued function defined in neighbourhood of the origin in R 3. f is sid to be positively homogeneous function of degree λ, if for ny number > 0 f x = λ f x. 37 Let k := λ. If λ is n integer nd if λ i.e. k is non-negtive integer, then the prtil derivtive of f with respect to the coordinte x i hs to be clculted by mens of the formul cf. Eq. 5.5 in 5 i f = i f ord k k! k δx x i... x i k Σ dσ i f x i...x i k, 38 where i f ord mens the derivtive computed using the stndrd rules of differentitions, Σ is ny smooth close surfce surrounding the origin nd dσ i is the surfce element on Σ.

26 6 As n exmple of pplying the rule 38 let us consider the full divergence connected with the fourth term in the H 34 prt of the 3PN Hmiltonin, see Eq. 58 d 3 x { } π3 i π3,il l π3 k π3,km m,j htt 4ij. 39,k Applying the regulriztion procedure from Appendix B. to the divergence 39 gives result much more complicted thn tht given by Eq. 75. After performing in 39 differentition with respect to x k one finds tht only in one term the rule 38 cn be pplied. The term reds The rule 38 pplied to π k 3,jk yields d 3 x 3 π3 i π3,il l 4ij πk 3,jk. 40 { } π3 i 4 π3,il l 4ij p j δ. 4 The integrl 4 is clculted by mens of Hdmrd s procedure from Appendix B.. After dding the result of regulriztion of the integrl 4 to the result of regulriztion of the integrl 39 we obtin multiple of the quntity 75. Using similr considertions we were lwys ble to restrict the mbiguity to multiple of the quntity 75. Let us lso stress tht the H 34 prt of the 3PN Hmiltonin is written in such form tht there is no need to use the rule 38 for individul derivtives ppering in H 34. A method of pplying the rule 4 in regulriztions of integrls cn be found in Section 6 of 5. We hve yet found tht this method without modifictions is not ble to give zero for ll full divergences one meets in our clcultions. References B. R. Iyer nd C. M. Will, Phys. Rev. D 5, L. Blnchet, Phys. Rev. D 55, P. Jrnowski nd G. Schäfer, Phys. Rev. D 55, T. Dmour, in Grvittion in Astrophysics, edited by B. Crter nd J. B. Hrtle Plenum Press, New York, 987, p A. Gopkumr, B. R. Iyer, nd S. Iyer, preprint 997, gr-qc/ T. Dmour, in Three Hundred Yers of Grvittion, edited by S. Hwking nd W. Isrel Cmbridge University Press, Cmbridge, 987, p T. Dmour nd G. Schäfer, Nuovo Cimento B 0, R. Arnowitt, S. Deser, nd C. W. Misner, in Grvittion: An Introduction to Current Reserch, edited by L. Witten John Wiley, New York, 96, p L. E. Kidder, C. M. Will, nd A. G. Wisemn, Clss. Quntum Grv. 9, L5 99; Phys. Rev. D 47,

27 7 0 N. Wex nd G. Schäfer, Clss. Quntum Grv. 0, G. Schäfer nd N. Wex, in Perspectives in Neutrinos, Atomic Physics, nd Grvittion, edited by J. Trân Thnh Vân, T. Dmour, E. Hinds, nd J. Wilkerson Editions Frontières, 993, p. 53. T. Dmour, B. R. Iyer, nd B. S. Sthyprksh, preprint 997, gr-qc/ T. Dmour, in Grvittionl Rdition, edited by N. Deruelle nd T. Pirn North-Hollnd, Amsterdm, 983, p T. Dmour nd G. Schäfer, Gen. Rel. Grv. 7, S. M. Kopeikin, Sov. Astron. 9, S. Wolfrm, Mthemtic. A System for Doing Mthemtics by Computer, nd ed. Addison-Wesley, Redwood City, G. Schäfer, Annls of Physics NY 6, T. Kimur, Prog. Theor. Phys. 6, T. Regge nd C. Teitelboim, Annls of Physics NY 88, G. Schäfer, in Symposi Gussin, Proceedings of the nd Guss Symposium, Conference A: Mthemticl nd Theoreticl Physics, Munich, edited by M. Behr, R. Fritsch, nd R. G. Lintz Wlter de Gruyter, Berlin, 995, p T. Dmour nd G. Schäfer, J. Mth. Phys. 3, C. Itzykson nd J. B. Zuber, Quntum Field Theory McGrw-Hill, New York, 980, p T. Kimur nd T. Toiy, Prog. Theor. Phys. 48, T. Oht, H. Okmur, T. Kimur, nd K. Hiid, Prog. Theor. Phys. 5, 0 974; Nuovo Cimento B 7, P. Jrnowski, in Mthemtics of Grvittion. Prt II. Grvittionl Wve Detection, Bnch Center Publictions, Vol. 4, Prt II Institute of Mthemtics, Polish Acdemy of Sciences, Wrszw, 997, p M. Riesz, Act Mth. 8, I. M. Gel fnd nd G. E. Shilov, Generlized Functions, Vol. Acdemic Press, New York, 964.