arxiv:gr-qc/ v2 10 Jun 2003

Transcription

1 The binry blck-hole dynmics t the third-nd--hlf post-newtonin order in the ADM-formlism Christin Königsdörffer, Guillume Fye, nd Gerhrd Schäfer Theoretisch-Physiklisches Institut, Friedrich-Schiller-Universität Mx-Wien-Pltz, 043 Jen, Germny Dted: November, 08 rxiv:gr-qc/ v 0 Jun 003 We specilize the rdition-rection prt of the Arnowitt-Deser-Misner ADM Hmiltonin for mny non-spinning point-like bodies, clculted by Jrnowski nd Schäfer, to third-nd--hlf post-newtonin pproximtion to generl reltivity, to binry systems. This Hmiltonin is used for the computtion of the instntneous grvittionl energy loss of binry to PN rective order. We lso derive the equtions of motion, which include PN rective terms vi Hmiltonin nd Euler- Lgrngin pproches. The results re consistent with the expressions for rective ccelertion provided by Iyer-Will formlism in Ref. in generl clss of guges. PACS numbers: 04.5.Nx, Db, 9.60.Jd, 9.60.Lf I. INTRODUCTION Compct binries re the most promising sources of detectble grvittionl wves for ground nd spce bsed interferometric experiments, TAMA300, GEO600, VIRGO, LIGOLser Interferometer Grvittionl Wve Observtory nd LISA Lser Interferometer Spce Antenn. Becuse of grvittionl wve emission, the binry seprtion decreses dibticlly in time while the frequency increses ccording to the reltivistic version of the Kepler lw. This process ends s soon s the first instble orbit is reched before the colescence. During the inspirling phse the binry system produces most likely detectble signl. The construction of wve ptterns for dt nlysis requires the knowledge of the theoreticl wve form to very high post-newtonin orders see, e.g., Refs. 3, 4. This ssumes good understnding of the compct-binry dynmics. The seculr effects governing the long-time evolution re of prticulr importnce. They result from the rection force due to the qudrupole dmping rising t the.5 order in power of v /c, v representing typicl orbitl velocity nd c stnding for the speed of light in vcuum. For this reson, they re often referred to s the.5 post- Newtonin PN dmping. The mesurements resulting from the observtion of the Hulse-Tylor binry pulsr PSR 936 show excellent greement with the vlues predicted by Einstein theory 5. A while go, Iyer nd Will 6 derived the equtions of motion for two non-spinning point-like objects including the first dissiptive terms ppering t the.5 nd 3.5PN level, i.e. the qudrupole dmping nd its PN correction. They proceeded by postulting blnce eqution between the instntneous flux of energy nd ngulr momentum in the fr zone on the one hnd, nd the instntneous loss of energy nd ngulr momentum in the Electronic ddress: C.Koenigsdoerffer@tpi.uni-jen.de system snerzoneontheotherhnd. Inlterpper, they detiled their derivtion nd proposed n lterntive method, bsed on the speciliztion of the Blnchet s rdition-rection potentils to binry systems, which led to the sme results. More recently, the dynmics of two-body system t the 3.5PN order ws derived by Pti nd Will in hrmonic coordintes 8 using the integrtion method for the Einstein field equtions developed in Ref. 9. This rticle dels with the 3.5PN grvittionl dmping of compct binries regrded s point-msses. We bse our investigtion on the N-body point-mss Hmiltonin clculted by Jrnowski nd Schäfer within the Arnowitt-Deser-Misner ADM Hmiltonin formlism of generl reltivity 0, which hs proved to be very efficient in determining the pproximte 3.5PN dynmics,,, 3. A remrkble feture of the ADM formlism t this level is the bsence of symptotic mtching between the ner field nd the fr field given in ner-zone nd fr-zone-defined coordinte system respectively. We cn sty in one single globl coordinte systeminllourclcultions,likein Ref. 8. Letusmention tht the ADM formlism does not provide directly the blnce equtions between instntneous losses nd fluxes, even when ssuming qusi-sttionrity in the rdition emission. Although the ADM formlism gives us the conserved totl energy, the blnce between the lost energy of the mtter system nd the energy flux in the wve zone is proved up to 3.5PN order only see, e.g., Ref. 4. Themin purposesofthis workrei to specilizethe N-body 3.5PN ADM Hmiltonin of Ref. to N =, ii to deduce from it the energy loss of two-body system in generl orbits, nd iii to clculte the equtions of motions in the Hmiltonin form nd in the Euler- Lgrngin form in order to determine the rditionrection force. The expression obtined for the energy loss is different from tht derived in Jrnowski nd SchäferndinIyerndWill, butthe time-verged expressions coincide in the cse of qusi-elliptic orbitl

2 motion. The pln of the pper is to reduce the originl Hmiltonin for two point-msses. Next, we perform prtil time differentition to get the energy loss. On the other side, we compute the rdition-rection force from the Euler-Lgrngin equtions of motion. The corresponding energy loss turns out to be equivlent to the preceding one fter verging. Our results re compred with those of Iyer nd Will nd the ccessory guge coefficients re determined explicitly. Hdmrd regulriztion is systemticlly employed. We use units in which 6πG = = c, where G is the Newtonin grvittionl constnt, but G nd c re restored in the presenttion of our finl results. In our nottion, t is the coordinte time. Chrcters in bold, like x = x i, represents point in the 3-dimensionl Eucliden spce R 3 endowed with stndrd Eucliden metric. The sclr product is denoted by dot. Letters,b,... re prticle lbels, so tht x R 3 denotes the position of the th prticle. We lso define: r := x x, r := r, n := r /r ; nd for b: r b := x x b = r b r, r b := r b, n b := r b /r b ;... stnds here for the Eucliden length of vector. The momentum vector of the prticle with mss prmeter m is denoted by p = p i. A dot over symbol, like in ẋ, represents its totl time derivtive. The prtil differentition with respect to x i is denoted by i, or equivlently by comm, i.e. i φ φ,i. II. 3.5PN FIELD EQUATIONS The ADM formultion is bsed on 3 splitting of spce-time. The cnonicl vribles re, on one hnd, the projection g ij of the 3-metric on the spce-like hypersurfces x 0 = t = const. in suitble coordinte grid nd, on the other hnd, its conjugte momentum π ij. In the so-clled ADM guge,, 5, the sptil metric g ij decomposes into digonl prt 8 φ 4 δij, plus prt h TT ij trnsverse nd trceless with respect to the Eucliden metric δ ij, so tht g ij = 8 φ 4 δ ij h TT ij. with j h TT ij = 0. In ddition, the field momentum π ij is trceless: π ii = 0. Let us consider now h TT ij s the new field vrible, nd denote its conjugte momentum by π ijtt. The degrees of freedom ssocited with π ij = π ij π ijtt nd the potentil φ re removed by solving the constrint equtions. After inserting the vlues of π ij nd φ into the generl reltivistic Hmiltonin of the isolted system under considertion in n symptoticlly flt spce-time, the ltter simplifies to surfce integrl often referred to s the reduced Hmiltonin H. It is functionl of the independent degrees of freedom h TT ij nd π ijtt 0. For system of N point-like bodies with position vectorsx ndmomentp =,...,N, the Hmiltonin tkes the form H = H x,p,h TT ij,πijtt.. The Hmiltonin equtions of motion for body re ṗ = H x nd ẋ = H p,.3 nd the field equtions for the independent degrees of freedom red t πijtt = δ TTij δh kl δh TT, kl t htt ij = δij TTkl δh δπ kltt,.4 where δ.../δ... denotes the Fréchet derivtive see, e.g., Ref. 6 nd where δ TTij kl is the TT-projection opertor see, e.g., Eq.. of. The Hmiltonin H cn be uniquely decomposed into mtter prt, H mt, depending on the mtter vribles only, field prt, H field, depending on the field vribles only, nd n interction prt, H int, depending on both sets of vribles, mtter nd field, in wy tht it vnishes if one set is put to zero. In short, the full content of the field-plus-mtter dynmics t the 3.5PN order is included in the Hmiltonin H 3.5PN = H mt 3.5PN Hfield 3.5PN Hint 3.5PN..5 The nottion 3.5PN indictes tht ll orders lower or equl thn 3.5PN re included. In wht follows, only the interction prt of the Hmiltonin.5 will be needed. It hs been explicitly computed t the dissiptive 3.5PN order in Ref. by Jrnowski nd Schäfer. Its trunction t the.5pn level, H.5PN int, tkes the wellknown form presented in Refs., with H int.5pnx,p,t = 5π χ 4ij tχ 4ij x,p.6 χ 4ij x,p := 60π 6π p m δ ij 3p i p j b m m b r b 3n i bn j b δ ij.. Thefunction χ 4ij tcomesfromthepost-newtoninexpnsion of the retrded integrl giving h TT ij s function of x nd p. As H/ x nd H/ p re differentited for constnt h TT ij nd π ijtt, the positions nd moment ppering in Eq.. re not ffected by the opertors / x nd / p. In order not to confuse the ltter prticle vribles with the ones tht re lredy present in the originl Hmiltonin. nd re ffected by / x nd / p, we shll mrk them with prime symbol: x, p. Thus, in our nottion, χ 4ij t denotes the time derivtive of χ 4ij x,p.

3 The 3.5PN Hmiltonin H3.5PN int reds H int 3.5PNx,p,t =5πχ 4ij x,p H int 3.5PN Π ij t Π ij t Π 3ij t 5π χ 4ij tπ 3ij x,p χ 4ij t 5π χ 4ij t Q ij x,p,tq ij x,t Π ij x,p Π ij x,t 3 t 3 R x,p,tr x,t..8 is prmetrized by the following functions Π ij x,p := 4 p 5 6π m 3 5 6π p δ ij 3p i p j 8 5 r b b 6π b 9p p b 3n b p n b p b δ ij 4p i p bj } 3 5p p b n b p n b p b n i bn j b 6n b p b n i bp j n j b p i π b m m b r b π 80 r b s bc 0 s bc 4 b c,b s bc r b r bc r c δ ij 3n i bn j b 8 m m b m c δ ij 3n i b r b r nj b c n i b nj b b n i b s nj bc bc r b r bc r c r b r bc r c m b p m r δ ij 5p i p j p ni b nj b b } δ ij,.9 with s bc := r b r bc r c, Π ij x,p := m b 5n b p p δ ij p i p j 5p 5 6π m r 3n b p n i b nj b b Π 3ij x,p := 5 } 6n b p n i bp j n j b p i π 3 0 6π 3 r b r 3 c b c,b r bc r b rc 3 5 r b r c rbc 3 6 6π b r b r c 40 r b b m m b m c 5 r c rb 3 r c r bc s bc r c r b 88 s bc s bc n i bn j c m m b r b n i bn j b 3n i bn j b δ ij 3 40 δ ij r b r c r b s bc },.9b m b 5n b p δ ij n b p n i b nj b n i b p j n j b p i,.9c 3 We observe tht there re few misprints in the expressions for H3.5PN int, ppering in. We mke pproprite chnges in those expressions nd their positions re mrked by....

4 Π ij x,t := m 5 6π m r 5n p p δ ij p ip j 5p } 3n p n i nj 6n p n i p j n j p i 3 m m m b 0 6π b 6 r b s b r r r 3 3 b r b 3 s b n i nj b nj ni b 5 r b r 8 s b r s b 5 r r 3 r 4 b r b r b r r r b n i r b s b nj b b r r b r 3 b n i nj s b n i nj b 8 } 4 δ ij,.9d s b r r b with s b := r r b r b. VriblesorigintingfromhTT ij nd π ijtt reprimed. It isindeed cruciltodistinguish the prticle positions nd moment inside nd outside the trnsverse-trceless TT quntities t this stge. The other four functions which enter H3.5PN int re Q ij x,p,t := m p i p j n p n i 6 6π m r p j 5p ni nj 3n p n i nj,.9e Q ij x,t := 3 m m b m 3 6π s b b 8 s b r s b with s b := r b r r b, R x,p,t := 05 6π n i nj r r b r 3 b n ib r b s njb 3 r b rb 3 5 r b r b r rb 3 r s n i nj b 3 b r } n ib b s, nj b.9f r 5p m m p p p 4n p p 4n p p n p n p p p 05 6π r r b r b nd R x,t := m m b 05 6π m b m m b m b n b p 4 r 6r n p 3 r b rb r r4 r 3 r r b b r 3 5 r p 4 r n p b r b r b 6 r3 r r n p n b p b 5 r r b r r b r 3 b 3 m m b m m b 0 6π b b r4 r rb 3 p r b n b p r b n b p n p r r b r 3 b r r b r r 3 b r b r r b 4.9g

5 4 r br r 3 5 r r 3 b b r b r b r 3 r b r b r b r b r b r b r b r b r b n b n b 4 r br r bb r 3 n n bb 6 r4 b r n b n b r b 34r b r n b r b n b..9h b 5 Now, we re in position to reduce the N-body 3.5PN Hmiltonin.8 to N =. The result will serve s bse for the derivtion of the grvittionl energy loss of compct binry moving in generl orbits, like in Ref., nd provide reference to check the equtions of motion tht will be clculted fterwrd from this two-body Hmiltonin. III. ENERGY LOSS OF A TWO-BODY SYSTEM Let us denote by H 3.5PN the Hmiltonin resulting from H 3.5PN, Eq..5 fter writing the field prt s function of x,p,ẋ,ṗ nd subtrcting its.5pn nd 3.5PN contribution H field By definition, we hve.5pn nd Hfield 3.5PN from Hfield 3.5PN. H 3.5PN :=H mt 3.5PN H field 3.5PN H field.5pn H field 3.5PN H int 3.5PN, 3. which cn be split into conservtive nd dissiptive prt H 3.5PN x,p,ẋ,ṗ,t =H cons 3PNx,p,ẋ,ṗ H diss 3.5PN x,p,t, 3. where H 3PN cons :=HN mt HPN mt HPN mt HPN field HPN int H3PN mt H3PN field H3PN int, 3.3 H diss 3.5PN :=Hint.5PN Hint 3.5PN. 3.3b The conservtive Hmilton function corresponding to the expression 3.3 is given in Refs. 8, 9. With the help of the coordinte trnsformtion described in Ref. 3, it is possible to eliminte the dependence on ẋ nd ṗ. The difference between H 3PN cons nd Hdiss 3.5PN in Eq. 3. All lgebric computtions re performed using the commercil softwre Mthemtic 0. reltes to their behviors with respect to time reversl: H 3PN cons is symmetric wheres Hdiss 3.5PN is ntisymmetric. Thetotltime derivtiveof H 3.5PN isequltoitsprtil time derivtive, the other terms cnceling by virtue of the equtions of motion. Only the dissiptive prt of H 3.5PN depends explicitly on time, through the substitution of the trnsverse-trceless vribles tht hve been replced by their ctul vlues. We get therefore d dt H 3.5PN x,p,ẋ,ṗ,t = t H 3.5PN x,p,ẋ,ṗ,t = t Hdiss 3.5PN x,p,t. 3.4 The instntneous energy loss of the mtter system due to the grvittionl wve emission is defined s L 3.5PN t := d dt H 3.5PN x,p,ẋ,ṗ,t. 3.5 Out tsk is now to compute L 3.5PN for system of two point-like bodies with msses m nd m. We strt by rewriting the right-hnd side of Eq. 3.5 by mens of Eqs. 3.3b nd 3.4 in the form L 3.5PN = t H int.5pn Hint 3.5PN. 3.6 Next, weinserttheformuls.6nd.8intoeq.3.6, nd specilize ll sums of Eqs.. nd.9 for N =. Notice tht the triple sum disppers in reltions.9 nd.9b. For clrity, we lso reintroduce b c,b the Newtonin grvittionl constnt G nd the speed of light c. In differentiting the Hmiltonins.6 nd.8 with respect to time, we replce the time derivtives of the primed coordintes nd moment ccording to the PN equtions of motion: ẋ = p m c m 3 p p G c 6 m r m ẋ = p m c m 3 p p G c r p p n p n, 3. 6 m p p n m p n, 3.b

6 ṗ = Gm m r n G m m m m c r 3 n G c r n p p n p p 3 m } p m 3m p m p p 3n p n p n, 3.c ṗ = Gm m r n G m m m m c r 3 n G c r 3 m m n p p n p p } p 3m p m p p 3n p n p n. 3.d The clcultion of the third time derivtive of R nd R is prticulrly tedious. The prticle moment p / is expressed in terms of the prticle coordinte velocities v / up to the PN order see, e.g., Eq. 4. of by: p = m v c m v v Gm m c 6v v n v n, 3.8 r p = m v c m v v Gm m c 6v v n v n. 3.8b r Anlogous reltions hold for p nd p. The time derivtives of the primed coordintes x nd x coming from the differentition of H 3.5PN in Eq. 3.5 equtethe primed coordintevelocitiesv nd v. We elimintell ccelertionsby mkinguseofthe PNequtions of motion see, e.g., Eq..5 of 5: v = Gm Gm r n c r 4n v 3n v v v v v 4v v } 3 n v G 5m 4m n, 3.9 r v = Gm Gm r n c r 4n v 3n v v v v v 4v v } 3 n v G 5m 4m n r. 3.9b 6 At this stge, we cn identify positions nd velocities of prticles inside nd outside the trnsverse-trceless vribles i.e. the primed nd unprimed quntities. As the limit x x, x x is singulr, we shll resort to the Hdmrd regulriztion procedure described in Appendix Afor more detiled investigtion, see pper nd the references therein. The computtions re extremely long, but strightforwrd. To write the finl formul for the energy loss in more compct wy, we introduce the totl mss of the system M := m m, the reduced mss µ := m m /M, s well s the rtio prmeter ν := µ/m. The individul msses m nd m re given in terms of µ nd ν by: m = µ ν 4ν, 3.0 m = µ ν 4ν, 3.0b ssuming m m. We lso express the velocities of the bodies in the center-of-mss frme s function of their reltive velocity v = v := v v nd of their seprtion r = r := x x = r r with v = ṙ. This is chieved by mens of the reltions linking the positions ofthe individul prticles to r t the PNordersee, e.g., Eq..4 of 3: 3 x = µ µm m m c M v GM r r, 3. 3 At the PN order the Eqs. 3. in ADM guge nd in hrmonic guge re identicl.

7 x = µ µm m m c M v GM r r, 3.b νv v Gµ r vn vn 3.3 where r = r. By differentiting Eqs. 3. with respect to time nd eliminting the temporl derivtives with the id of the Newtonin eqution of motion where n := r/r, we re led to: v = GM n, 3. r v = νv 4ν 4ν c νv v Gµ vn vn. 3.3b r v = νv 4ν 4ν c Finlly, we rrive t the instntneous energy loss 4: L 3.5PN = 4 5 G M 3 ν c 5 r 3 G M r 05 G M 3 ν c r 3 v 9n v GM v 4 60n v v 45n v 4} r 5 04ν G3 M 3 r 3 656νv 339νn v G M 3083νv νn v v νn v 4 GM r 64064νv νn v v νn v 4 v 058 νn v 6}. 3.4 r The comprison of the expression bove for the instntneous grvittionl energy loss with the known instntneous fr-zone flux see Eq. 5 of 5 corrected by n errtum nd eqution below Eq.3.40 of6 shows tht they re not either identicl. This is nothing to worry bout since the form of the instntneous grvittionl energy loss in n isolted system is s well coordinte s representtion dependent. We lso observe tht the expression 3.4 is different from formul 3 in the work of Jrnowski nd Schäfer t the 3.5PN order. However, the difference my be expressed s totl time derivtive. To verify tht, we computed L 3.5PN nd LJS 3.5PN, the orbitl verge of L 3.5PN nd LJS 3.5PN respectively, using the prescription given in Ref. 6, for qusi-ellipticl orbitl motionnd found L JS 3.5PN = L 3.5PN. So, werecover Eqs. 4.0 nd 4. of Ref. 6. The importnt point is tht ll verged losses coincide. In the following, we shll mke the ssumption tht the binry hs entered the inspirling phse, so tht its reltive motion is qusi-circulr. The investigtion ppers to be simpler in this cse, since ll terms involving the fctor n v vnish. The verging is performed by expressing L 3.5PN s function of E. According to eqution 3.35 in Ref. 6, we hve for the energy reduced by the fctor µ E = v GM r GM rc 3 8c 3νv4 3νv νn v GM r. 3.5 The right-hnd side depends on the reltive velocity v nd the reltive distnce r, wheres the symbol E in the left-hnd side my be viewed s n dibtic prmeter. By solving Eq. 3.5 for v, we find v =E GM r c E 39ν GM r E 6ν G M r 0ν. 3.6 Next, we substitute Eq. 3.6 to v in the instntneous energy loss 3.4, nd eliminte the reltive distnce r in fvor of the reduced energy with the id of Eq. 43 in Ref. : r = GM E νec. 3.

8 8 As consequence, the time-verged energy loss specilized to the cse of circulr orbits c.o. reds c.o. L 3.5PN = 04 ν E 5 t 5 c 5 E 3 G c 68 5ν 3.8 We recover formul 4.0 in Ref. 6 with e R = 0 due to the fct tht the orbits re qusi-circulr. This is in greement with Eq.8 of Ref.5. Processing Eq.4.6 or Eq. 3 in Ref. in similr wy ends in the sme result. The instntneous grvittionl energy loss L 3.5PN given in Eq. 3.4 will be tken s reference to check the equtions of motion tht will be determined in the next section.. IV. HAMILTONIAN EQUATIONS OF MOTION In this section we compute the rection prt of the equtions of motion ṗ = Hint 3.5PN nd ẋ = Hint 3.5PN, 4. x p in fully explicit form in terms of x,x,p nd p. For this gol, we hve to differentite the Hmiltonins.6 nd.8 with respect to the prticle coordintes x nd the prticle moment p with =, nd pply the Hdmrd regulriztion procedure described in Appendix A. We find 4: ẋ = c 5 G r 4 5 n p n p m 4 m 5 n p 6 5 p 6 5 p p n 8m n p 3m 8 5 n p p 6 } } 5 n p p G m 36 c r m 3 n p 4 60 n p p 4 35 p4 m 48 n p 3 n p n p n p p n p p p 35 p p p 04 m m n p n p n p p n p n p p p p p 3 35 n p p 05 p p 8 m n p n p n p p p 6 35 n p n p p 5 05 p p p m n m n p n p p 64 m n p n p n p p n p p p 6 m m n p n p n p p p n p p m n p n p p p m 5 n p n p p 99 m m 35 n p n p n p p 400 n p p p 9 m 35 n p n p n p p p 5 } 05 n p p p G3 m r 3 86 m 05 n p p m 48 5 n p n p n p p p p m 05 n p n p 58 5 n p 484 p p 9 6m 588 n p m 5m 05 n p 43 n p m n p p 5 p n

9 469 } 35 m n p m 484 n p 03 5 n p p } G4 r 4 m m m m 3 n, 4. 9 ẋ =. 4.b For the time derivtive of the prticle moment, we get ṗ = c 5 c G 3 56 r 4 5 m m n p 56 5 m m n p n 6 5 m m p 6 } } 5 m m p G m r 3 m 3 0n p 5 6 n p 3 p n p p 4 30n p 4 n p 80 n p n p p 8 m n p p p p m m n p p 4 84 n p 3 p p 30n p 3 n p 6n p n p p 44n p n p p p 8 35 n p p p p n p p p 30 n p 3 p n p p p m 4 n p n p p p n p p p 30 m m 0n p n p 3 n p 3 p n p n p p n p p p 4 35 n p p p p m n m 3 6 n p 4 4n p p 6 40 m n p 3 n p 6 5 n p n p p 8 5 n p p p 3 05 p p p n p n p n p p 8 5 n p n p p p 6 p p n p p 05 p p 4 5 n p n p p 4 5 p p p m m m m p m 05 p4 n p n p n p p p 6 n p 4 5 n p p 35 p4 n p 3 n p 35 n p n p p n p p p 8 5 p p p 6 n p n p n p p n p p 3 p p p } G3 m r 4 m 35 n p n p p p 40 p p 5 5 n p n p p 0 m n p n p n p 4 05 n p p n p p n p p p m 0 n p n p n p n p n p p 4 35 n p p p n p p p 4 05 n p p m 48 m n p 3 34 m 35 n p p n 6 m n p 8 05 p

10 m 8 05 n p n p n p 8 35 p p p 8 m 35 n p n p 58 5 n p 8 05 p p p p 34 m n p p m 05 n p n p n p 5 05 p 9 05 p p m 30 n p 4 5 p p } G4 r m3 m n p m m 69 m m m m 3 ṗ =. m n p 38 p 05 n p n 6 35 m3 m 64 5 m m } } p, m m 3 n p b To check the correctness of the expressions for ẋ, ṗ, we employ them to compute the instntneous energy loss of binry, t.5pn nd 3.5PN orders by differentiting PN ccurte conservtive Hmiltonin expression see, Hmiltonin.5 of Ref. 5 truncted t the PN level. 4 The contribution to the Hmiltonin bove comes purely from the mtter prt nd it reds H PN mt x,p = p p Gm m m m r 8c Differentition with respect to time leds to d dt Hmt p p m m n p n p m m PN = Hmt PN x PN ẋ Hmt x p 4 m 3 p4 m 3 Gm m p rc 3 m p m G m m m m r c. 4.4 PN ẋ Hmt p ṗ Hmt PN p ṗ. 4.5 We then replce the time derivtives of x, x, p nd p occurring in Eq. 4.5 by their expressions 4. nd 4.3. At lst, we write the prticle moment in terms of the prticle coordinte velocities with the help of reltions 3.8. We introduce the reltive velocity, the reltive position, the totl mss nd the prmeter ν, nd multiply by ccording to definition 3.5. In this wy, we obtin for the luminosity 4: } L 3.5PN = d dt Hmt PN = 4 G M 3 ν 5 c 5 r 3 4v n v GM r G M 3 ν 05 c r 3 05 G3 M 3 r νv 4880νn v G M r 94 84νv νn v v νn v 4 GM r 6 344νv νn v v νn v 4 v 050 νn v 6}. 4.6 Though this result differs from the instntneous energy loss 3.4, one will see by using the method presented 4 The PN terms do not ply ny role in the computtion of the.5pn nd 3.5PN energy loss, since they first couple with the rection force t the.5pn = 4.5PN order. in Ref. 6 for verging tht it leds to the sme verged vlue, in greement with Eqs. 4.0 nd 4. of

11 Ref. 6. V. EULER-LAGRANGIAN EQUATIONS OF MOTIONS The im of the present section is to compute the PN rective corrections to ẍ employing two slightly different methods. We will see tht the computtion of these Euler-Lgrngin equtions of motion does not involve Lgrngin in either of these methods. However, we will require PN ccurte reltions between prticle moment nd its coordinte velocities. We turn to derive these reltions in the next subsection. A. Legendre trnsformtion For the present computtion, we require PN,.5PN nd 3.5PN corrections to the fmilir p = m v reltions. So, we hve to invert the equtions 4. giving ẋ = fx,x,p,p nd ẋ respectively, in order to determine p = fx,x,ẋ,ẋ = fx,x,v,v nd p. This will provide the link between the prticle moment nd the prticle coordinte velocities up to the order /c. As the first step, we dd Newtonin nd PN contributions to ẋ to the right-hnd side of Eqs. 4.. These corrections re obtined from Eqs. 3. nd 3.b fter replcing primed quntities by their unprimed counterprts. Through the expression presented below refer to the first body, it should be noted tht similr expressions hold true for the second body. It is the full coupled system tht must be inverted. The time derivtive of the position x hs the form ẋ = v = p m c A c 5A 5 c A, 5. wheres the structure of p fter inversion is p = m v c B c 5B 5 c B. 5. The coefficients B,B 5 nd B in Eq. 5. re computed itertively in terms of A,A 5 nd A of by mens of Eq. 5.. The first-order coefficient B is red directly from Eq We get 4: p = m v c m v v Gm } m n v } n6v v c G m m r 5 r 4 5 n v 4 5 n v n v 6 5 v v v n 3 n v 8 5 n v v 6 } } 5 n v v c v4 G m m r 36 n v 4 48 n v 3 n v 04 n v n v 8 n v n v 3 n v n v n v 9 35 n v 3 35 n v 5 n v n v 44 v v n v v v 05 v v 6 n v n v n v 3 8 n v n v v v G3 r 3 n v v v v v n v 6 n v n v n v n v n v 35 n v n v v n v 36 5 n v v v 5 n v n v n v 9 35 n v n v n v n v v v 44 } n v v v m 3 m 68 5 n v 04 n v n v 9 v 4 05 v v 9 5 n v n v v

12 336 m m 05 n v 8 05 n v n v 6 5 n v v 0 35 v v 3 m m 5 5 n v 3 n v m 3 m n v 3 05 n v v m3 m n v m m 05 n v 3 n v v } G4 r 4 5 v n m 3 m m m3 n} }, 5.3 p =. 5.3b B. First method to clculte the Euler-Lgrngin equtions of motions The clcultion of the ccelertions mkes the connection with the Euler-Lgrngin pproch. It is chieved in first stge by evluting the totl time derivtives of the functions ẋ = fx,x,p,p nd ẋ = fx,x,p,p s given by the Eqs. 4.. It is worth stressing tht this method does not need the explicit knowledge of the Lgrngin. We cn write indeed: ẍ = d dtẋ = ẋ x ẋ ẋ x ẋ ẋ p ṗ ẋ p ṗ, ẍ = d dtẋ = ẋ x ẋ ẋ x ẋ ẋ p ṗ ẋ p ṗ b Inserting pproprite PN contributions to ẋ nd ṗ given in Eqs. 4. nd 4.3, we obtin PN,.5PN nd 3.5PN corrections to the fmilir Newtonin expression for ẍ nd it reds prt from the PN nd 3PN terms 4: ẍ = Gm r n c G r m m p 4 p p 3 m m n p p 3 n p p 4 n p 3 n p m m m G m c 5 r 3 m 4n p n p p p m m 64 m 8 5 p G3 6 r 4 c 56 5 n p p p } G r 3 n 4 m 5m m 4m n m } n p m 48n p n p 5 n p p 5 n p p p n p 4n p 4 5 p n 5 n p 88 5 p m 56 5 n p n p 6 3 p p m 8 5 n p p m 5 n p 3 5 p m 5 n p n p 3 } 5 p p p 8 5 m n p m m 5 n p n p n 5 m m 8 } 5 m p G m r 3 m 4 46n p 5 9n p 3 p n p p 4 m 3 n p n p p n p p n p 3 p p 45 m m n p n p p p 44n p p p p 3 n p p 0 p p 0p p 4n p 4 n p 35 n p p p p 8n p 46 n p m m 4n p 3 p p n p 3

13 n p 06n p 3 86 n p p 35 p p p p n p 03 n p p p 3 35 p p p n p 84 m 3 n p 3 p p n p p p p 8n p 4 6 n p p 9 35 p4 n p n m 334 m 4 n p n p p 386 p4 m 3 5 n p 3 n p 88 n p n p p 64 5 n p p p p p p m 3 08 n p n p p 6 5 p4 40 p p 0 p p p p p n p p p4 4 n p p 96 m m 35 n p n p p p 6n p p n p m m 60 n p p p 584 n p n p p n p p m 3 8n p m m n p 3 n p n p n p p n p p p 396 p p p 54 m m 44 n p n p p p 48 p p 6 35 p p 5 m 3 35 n p p p 44 p p p G3 r n p n p m 8084 m n p p 394 n p n p n p n p p n p p 35 p p n p } n p p p p n p m 346 n p n p n p 384 m n p p p p p n p n p 95 3 p p p n p m 6464 m 35 n p n p p n p 84p p p n p n n p n p 68 5 n p p p p m m n p p m 064 m 05 n p 64 5 n p n p 8 5 p 3 05 p p m 8 m 5 n p p p 05 n p n p n p 348 p p p m 3005 m 05 n p p m 84n p n p 548 m 05 n p 04 5 p p 384 } 05 p p G4 r m3 n p m m 35 n p n p m m 8 35 n p 35 5 n p n 3

14 m m 6 05 m m m3 p m m 6 } } 35 m m p, ẍ = b We now go on to present ẍ in terms of x,v. As Eqs. 5.5 involve p t PN,.5PN nd 3.5PN orders, we use only Newtonin, PN,.5PN contributions to p, given in Eqs. 5.3 rther thn Eqs. 3.8 truncted t the PN pproximtion. After some hevy lgebr, we get ẍ = fx,x,v,v nd ẍ = fx,x,v,v 4: ẍ = Gm r n c G m r 3 Gm 3 r 5m 4m n n v v 4v v v }c5 4n v 96 5 v 4v v n v 88 5 v 6 3 v v v v c v } G m m r 3 4n v 3 48n v 5 v 5 v v G3 m m 8 r 4 5 m n v m 64 n v n 5 n v } n 4n v 3n v v v 5 v 5 v v n v 5 n v 56 5 n v n v 6 5 n v 8 5 n v 5 n v n v 3 G m m r 3 46n v 5 4n v 4 n v 06n v n v 3 4v 84 v v n v 3 38 v 03 n v 3 5 v 8 n 5 m 8 } 5 m v } v v 84 v n v n v 8n v 40 v 45 v v v 56v 536 v v 46 v 86 v v v v v v v v v v v n v n v n v 8n v v v v v 0 v v 38 5 v v v v v 5 v4 n v n 334 n v 4 5 n v 3 n v 584 n v n v 3 08 n v v v v v 40 v v 6n v v v v 548 v n v n v n v v 64 5 v v v 8n v n v 3 n v 84 n v n v v v v v v G3 m m n v n v 05 v v v v v v 60 v v v 5 35 v v 65 v v 0 v v v v v 394 n v v v v v v v v 48 r 4 m 35 n v n v n v 35 8 n v v 4 3 v v v 088 n v m n v 3 35 v 5 v4 n v v n v 35 v 44 v v } n v v n v

15 36 35 n v n v 3 35 n v v v v v n v 5 5 n v 84 5 v 6 05 v v v n v v v v v n m n v 58 n v n v m n v 04 n v n v n v v v v v v m 05 n v n v n v v v v m 05 n v 54 n v n v n v v v v 6 } 05 v v G4 m m 8 r 5 35 m n v m 5 5 n v 64 5 n v m m 96 n v n v n 8 05 m m m m v 35 m m 38 } } 05 m v, 5 n v 5.6 ẍ =. 5.6b In the next subsection we show the second wy to compute the equtions bove. C. Second method to clculte the Euler-Lgrngin equtions of motions The second route is less strightforwrd, but leds to the confirmtion of the preceding result. We strt from the Hmiltonin equtions.3 for ṗ : ṗ = H x, 5. where we replce ṗ by the derivtive of the right-hnd side of Eq The forml time derivtives of ẋ nd ẋ, denoted by nd respectively, re kept unevluted. The right-hnd side of Eq.5. involve the rection contributions Eqs. 4.3, s well s the Newtonin nd post-newtonin prts 3.. We systemticlly eliminte the prticle moment in fvor of the prticle coordinte velocities by mens of Eqs. 5.3 without /c. To clculte ẍ, we consider the coefficient of in ṗ, nd invert it up to the order /c. The eqution of motion 5. in its present form is multiplied by the ltter result. Next, we gther ll terms but in the right-hnd side. Simultneously, the ccelertion ẍ is isolted in the sme wy. We rrive t the coupled system: =fx,x,ẋ = v,ẋ = v,ẍ =,ẍ =, 5.8 =fx,x,ẋ = v,ẋ = v,ẍ =,ẍ =. 5.8b The vrious contributions re sorted in powers of /c, nd lbeled with two different indices. The second index refers to the /c exponent, wheres the first index identifies the body; N stnds for the Newtonin prt. =N 0 c f, c 5f 5, c f,, 5.9 =N 0 c f, c 5f 5, c f,. 5.9b The decoupling is chieved by itertive substitution of Eq. 5.9 into Eq. 5.9b nd vice vers. We get =N 0 c g c 5g 5 c g, 5.0 =N 0 c g c 5g 5 c g. 5.0b Next, we compute nd itertively, nd insert the result into the right-hnd sides of Eqs The ccelertion hs finlly the following structure =fx,x,ẋ = v,ẋ = v, =fx,x,ẋ = v,ẋ = v b We recover Eqs In order to confirm the correctness of expressions bove for ẍ nd ẍ, we shll check it yields the correct energy loss 3.8. For this gol, we begin by deriving the rection forces with the help of Eqs. 5.6 F rec = m ẍ, 5.

16 6 F rec = m ẍ. 5.b After multiplying these reltions by ẋ = v or ẋ = v respectively, we sum to obtin m ẍ v m ẍ v = hx,x,v,v 5.3 This is not yet the instntneous energy loss due prticulrly to the presence of the Newtonin nd post- Newtonin terms coming from Eqs We then must serch for the function g = gx,x,v,v to be dded to h s defined in Eq. 5.3 in order to isolte the rditive.5pnnd 3.5PNprts, which does not only mount to removing the Newtonin nd post-newtonin contributions. We shll then be ble to find the verge power by substituting to v its vlue 3.6. The correct energy loss is given by the time derivtive of E = Hx,x,p x,v,p x,v truncted t the PN level. The prticle moment of the PN Hmiltonin 4.4 re expressed by mens of the prticle coordinte velocities 5.3 up to the order /c 5 only. When performing the differentition of the resulting energy with respect to time, ẍ = nd ẍ = coming from the derivtives of ẋ = v nd ẋ = v re replced by expressions 5.6 except in the Newtonin term m v m v where they re left unevluted. The function g = gx,x,v,v is just the difference de/dt m v m v. By construction, hg contins no Newtonin or post- Newtonin terms. It ctully represents the instntneous luminosity L 3 3.5PN = de/dt, we intended to check, given by the negtive instntneous energy loss. 5 Theresultofthevergingprocedureinthecseofqusielliptic orbitl motion is in greement with our preceding results nd with Eqs. 4.0 nd 4. of Ref. 6. VI. COMPARISM WITH IYER-WILL FORMALISM In this section, we shll show tht the dissiptive prts of Eqs. 5.6 re comptible with the generic rection 5 The full result for L 3 3.5PN is presented in Ref. 4. force of Iyer nd Will depending on 6 guge prmeters, s the ones given in the Eqs d in Ref.. The first two prmeters t the.5pn order nd the six coefficients t the 3.5PN order correspond to the guge freedom. As mtter of fct, it is sufficient to prove tht for certin vlues of the unknowns, the ltter rection force coincides with ours. This requires to compute the reltive ccelertion = := = ẍ ẍ 6. from Eq. 5.5 for ẍ = fx,x,p,p nd ẍ = fx,x,p,p respectively nd to go to the center-ofmss frme p p = As the ccelertions we strt with re function of p nd p rther thn v nd v, we must simply set p to p nd p to p in = x,x,p,p see Eq. 6.. The next step consists in estblishing the link between the prticle moment p nd the reltive velocity v = v := v v in orderto eliminte the prticle moment from = x,x,p. We thus restore the Newtonin nd post-newtonin prt in Eqs. 4. nd clculte the reltive velocity v = v := v v = ẋ ẋ. 6.3 We mke the substitutions p p nd p p in v = vx,x,p,p. The reltion v = vx,x,p cn be inverted with the id of the itertive method described bove to get the momentum s function of the reltive velocity p = px,x,v. Inserting p = px,x,v in = x,x,p leds to = x,x,v. After introducing the totl mss M s well s the dimensionless mss prmeter ν through Eqs. 3.0 nd µ = νm, the reltive ccelertion reds 4: = GM r n c G M c 5 r 3 G M c r 3 GM r 3νv 3 ṙ ν n4ṙ ṙνv } G M 4ṙ 3 ν 96 5 ṙv ν 46ṙ 5 ν 4ṙ 5 ν v ṙ4 ν 68 ṙ4 ν v 4 64 n ν ν 5 ṙ ν 88 5 v ν ṙν ṙν v v r 3 4νn } G3 M 3 6 r 4 5 ṙνn 8 } 5 νv v 56ṙ 3 ν 4ṙ3 ν n ṙ ν 48 5 ṙ ν } v

17 with the nottion G3 M r 4 ṙ 3 ν v ν 6 35 ν v ṙ3 ν v } G4 M 4 r ṙν ṙν 5 63 ṙν 5 05ṙν n n ν 6 5 ν 05 ṙ ν 44 ṙ ν v } 6.4 ṙ = n v = n v. 6.5 This result corresponds to the generic eqution of motion.4 nd Eqs. in Ref. by Iyer nd Will in the following equtions we tke G = = c: d r dt = M r n M r na PN A PN A 3PN ṙvb PN B PN B 3PN 8 5 ν M M r ṙna.5pn A 3.5PN vb.5pn B 3.5PN 6.6 r with see Eqs..,.b,.5 nd. in Ref. A PN = 3νv 3 νṙ ν M r, 6. B PN = ν, A.5PN = v M r 3ṙ, B.5PN =b v b M r b 3ṙ, A 3.5PN =c v 4 c v M r c 3v ṙ c 4 ṙ M r c 5 ṙ 4 c 6 M r, B 3.5PN =d v 4 d v M r d 3v ṙ d 4 ṙ M r d 5 ṙ 4 d 6 M r. 6.b 6.c 6.d 6.e 6.f Note tht we re not interested in the PN nd 3PN terms in the present work. Iyer nd Will estblish in Refs., 6 tht the energy nd ngulr momentum blnces hold if nd only if the coefficients i,b i,c i nd d i stisfies the following reltions: = 33β, = 3 3 α 3β, 3 = 5β, 6.8 b = α, b = α, b 3 = 3 3α 6.8b for the.5pn coefficients see Eqs.. of Ref., nd c = 8 3ν 3 β 3ν3δ 3δ 6, 6.9 c = ν 3α 4ν 3 β3ν δ 3δ 3δ 5 3δ 6, 6.9b c 3 = ν 5 β 3ν 5δ 5δ 4 5δ 6, 6.9c c 4 = ν 6αν β54ν δ 5δ 4 6δ 5, 6.9d c 5 = δ 4, 6.9e c 6 = ν α49ν3β4ν δ 3 3δ 5, 6.9f d = 3 3ν 3 α 3ν δ, 6.9g d = ν α5νδ δ 3, 6.9h d 3 = ν 3 3αβ 3ν 3δ 3δ 6, 6.9i d 4 = ν α38 ν 3β 3ν δ 4δ 3 3δ 6, 6.9j d 5 = ν 5β 3ν5δ 6, 6.9k d 6 = να3νδ 3 6.9l for the 3.5PN coefficients see Eqs..8 of Ref.. 6 The two prmeters α, β t the.5pn order nd the six coefficients δ i t the 3.5PN order correspond to the guge freedom nd hve no physicl mening. As mentioned in Ref. 8, t the.5pn order, the vlues 6 In Ref. the nottion ε 5 ws used in plce of δ 6.

18 8 α =, β = 0 chrcterize the guge of Dmour nd Deruelle 8, while α = 4, β = 5 correspond to the soclled Burke-Thorne gugesee for exmple 36. of9 lso used by Blnchet. It is then non-trivil check to verify tht the 8 coefficients i, b i, c i nd d i prmeterizing our.5pn nd 3.5PN terms Eqs. 6.c 6.f specilized for Eq. 6.4 yield unique, self-consistent solution for the 8 guge coefficients. They re 4: for the.5pn order: for the 3.5PN order: α = 5, β = 3; δ = ν, δ = ν, 6.0b 4 δ 3 = ν, δ 4 = ν, 6.0c δ 5 = ν, δ 6 = ν. 6.0d 4 VII. DISCUSSION In this pper, we hve computed PN rective contributions to the ADM Hmiltonin for pir of compct objects moving in generl orbits, using existing results vilble for the post-newtonin N-body Hmiltonin, computed by Jrnowski nd Schäfer. The two-body Hmiltonin is employed to compute the first PN corrections to the instntneous grvittionl energy loss ssocited with inspirling compct binries. We lso derive PN corrections to rective equtions of motion vi Hmiltonin nd Euler-Lgrngin methods. The expressions for PN rective ccelertions, computed in the ADM guge, re consistent with those of Iyer nd Will. Employing the vrious rective equtions of motion we lso computed instntneous grvittionl wve luminosity for the compct binry in generl orbits. These luminosities re in totl greement with their orbitl verged counterprts in the cse of qusi-ellipticl orbitl motion. Acknowledgments We thnk P. Jrnowski for discussions nd useful checks. We re grteful to A. Gopkumr for criticl reding of the mnuscript. This work ws supported by the EU Progrmme Improving the Humn Reserch Potentil nd the Socio-Economic Knowledge Bse Reserch Trining Network Contrct HPRN-CT nd by the Deutsche Forschungsgemeinschft DFG through SFB/TR Grvittionswellenstronomie. APPENDIX A: HADAMARD REGULARIZATION Let f be smooth rel-vlued function defined on R 3 deprived of point x 0 R 3, where it my be singulr. We consider the fmily of uxiliry functions f n ε := f x 0 εn, lbeled by unit vectors n. We expnd f n into Lurent series round ε = 0: f n ε = m= N m nε m. A The coefficients m of this expnsion depend on the unit vector n. We define the regulrized vlue of the function f t x 0 s the coefficient of ε 0 in the expnsion A verged over ll directions: f reg x 0 := 4π 0 ndω, A where dω is the elementry solid ngle viewed from the point x 0. This procedure is clled Hdmrd regulriztion. To clculte the integrls in Eq. A we used the following equtions n i n i...n i l dω = 0, for odd l, A3 4π 4π n i n i...n ip dω =... δ ip i p} p!! δii }}, p!! terms A3b with p!! = p. The curly brces in Eq. A3b men the symmetric combintion of the indices: δ ij δ kl} = δ ij δ kl δ ik δ jl δ il δ } jk. A4 } 3!! terms During the identifiction process of the primed nd unprimed vribles, x x nd x x, in the clcultion of the derivtion of the Hmiltonins given by the Eqs..6 nd.8, some singulrities pper. The purpose of the regulriztion is the elimintion of these singulrities. It is mtter of choice whether we pply first the regulriztion for index nd then for index, or the other wy round. As n exmple, we shll pply Hdmrd regulriztion procedure on the term n n n p n p r r3 = n i ni nj pj nk pk r r 3. A5 We intend to identify index with s first step. We strtbyclcultingthe limit for i = ri 0 fterthe replcement: r i = ri ri = ri i. A6

19 9 in Eq. A5. To this purpose, we expnd r m = f r i i = l=0 L l! Lr m = l=0 l r l l! nl Lr m. A We now perform the first ngulr integrtion over the terms of Eq. A5 depending on the vrible lbeled by. Here is m = 3 nd l =. 4π = n i nj r r 3 dω l=0 l l! = br 3 4π = br 3 = 30 By using we get n i 4π nj nl dω n i nj n nb dω r l Lr 3 δ ij δ b δ i δ jb δ ib δ j 5 ij r 3 δij r 3. A8 ij r m = mm n i n j mδ ij r m, r m = mmr m, A9 A9b n i 4π nj r r 3 dω = ni nj r 5. A0 In further computtions we replce the terms on the right-hnd side of Eq. A0 with r i = ri r i = ri i, n i = r r A n i r n i r = r n i n i r r. Ab The second step will consist in identifying index with. We clculte the limit for i = ri 0. The second ngulr integrtion over the finite prt terms i.e. the terms tht re of order zero in power of i depending on the vrible lbeled by cn be performed immeditely becuse there re no singulrities in Eqs. A for i 0: n i nj pj pk r 5 The result is n i 4π nk dω = ni n j pj pk r 5 n n n p n p r r3 3 δik. A x x x x = n p n p 3r 5. A3 P. Jrnowski nd G. Schäfer, Phys. Rev. D 55, B. R. Iyer nd C. M. Will, Phys. Rev. D 5, T. Dmour, B. R. Iyer, nd B. S. Sthyprksh, Phys. Rev. D 63, T. Dmour, B. R. Iyer, nd B. S. Sthyprksh, Phys. Rev. D 66, J. M. Weisberg nd J. H. Tylor 00, to be published in Rdio Pulsrs, ASP Conference Series, stroph/0. 6 B. R. Iyer nd C. M. Will, Phys. Rev. Lett. 0, L. Blnchet, Phys. Rev. D 4, M. E. Pti nd C. M. Will, Phys. Rev. D 65, M. E. Pti nd C. M. Will, Phys. Rev. D 6, R. Arnowitt, S. Deser, nd C. W. Misner, in Grvittion: An Introduction to Current Reserch, edited by L. Witten Wiley, New York, 96, pp. 65. T. Oht, H. Okmur, T. Kimur, nd K. Hiid, Prog. Theor. Phys. 5, G. Schäfer, Ann. Phys. N.Y. 6, T. Dmour, P. Jrnowski, nd G. Schäfer, Phys. Rev. D 6, G. Schäfer, in Very High Energy Phenomen in the Universe, editedbyy.girud-hérudndj. T.T. VânEditions Frontières, Pris, 99, pp T. Dmour nd G. Schäfer, Nuovo Cimento B 0, L. Schwrtz, Anlyse II, Clcul différentiel et équtions différentielles Hermnn, Pris, 99. G. Schäfer, in Symposi Gussin, edited by M. Behr, R. Fritsch, nd R. G. Lintz Wlter de Gruyter, Berlin, 995, pp P. Jrnowski nd G. Schäfer, Phys. Rev. D 5, P. Jrnowski nd G. Schäfer, Ann. Phys. Leipzig 9, S. Wolfrm, Mthemtic: A System for Doing Mthemtics by Computer Addison-Wesley Publishing Compny, Inc., 993, nd ed. T. Oht, H. Okmur, T. Kimur, nd K. Hiid, Prog. Theor. Phys. 5, L. Blnchet nd G. Fye, J. Mth. Phys. 4, T. Dmour nd N. Deruelle, Ann. Inst. Henri Poincré Phys. Théor. 43, C. Königsdörffer, Mster s thesis, Friedrich-Schiller- Universität, Jen R. V. Wgoner nd C. M. Will, Astrophys. J. 0, 64 96, Errtum: 5, L. Blnchet nd G. Schäfer, Mon. Not. R. Astron. Soc.

20 0 39, N. Wex nd G. Schäfer, Phys. Lett. A 4, , Errtum:, T. Dmour nd N. Deruelle, Phys. Lett. 8A, C. W. Misner, K. Thorne, nd J. A. Wheeler, Grvittion W.H. Freemn nd Compny, 999.