arxiv:gr-qc/ v2 10 Jun 2003
|
|
- Louise Spencer
- 5 years ago
- Views:
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.
arxiv:gr-qc/ v1 14 Mar 2000
The binry blck-hole dynmics t the third post-newtonin order in the orbitl motion Piotr Jrnowski Institute of Theoreticl Physics, University of Bi lystok, Lipow 1, 15-2 Bi lystok, Polnd Gerhrd Schäfer Theoretisch-Physiklisches
More information4 The dynamical FRW universe
4 The dynmicl FRW universe 4.1 The Einstein equtions Einstein s equtions G µν = T µν (7) relte the expnsion rte (t) to energy distribution in the universe. On the left hnd side is the Einstein tensor which
More informationClassical Mechanics. From Molecular to Con/nuum Physics I WS 11/12 Emiliano Ippoli/ October, 2011
Clssicl Mechnics From Moleculr to Con/nuum Physics I WS 11/12 Emilino Ippoli/ October, 2011 Wednesdy, October 12, 2011 Review Mthemtics... Physics Bsic thermodynmics Temperture, idel gs, kinetic gs theory,
More informationSUMMER KNOWHOW STUDY AND LEARNING CENTRE
SUMMER KNOWHOW STUDY AND LEARNING CENTRE Indices & Logrithms 2 Contents Indices.2 Frctionl Indices.4 Logrithms 6 Exponentil equtions. Simplifying Surds 13 Opertions on Surds..16 Scientific Nottion..18
More informationWeek 10: Line Integrals
Week 10: Line Integrls Introduction In this finl week we return to prmetrised curves nd consider integrtion long such curves. We lredy sw this in Week 2 when we integrted long curve to find its length.
More informationMath 1B, lecture 4: Error bounds for numerical methods
Mth B, lecture 4: Error bounds for numericl methods Nthn Pflueger 4 September 0 Introduction The five numericl methods descried in the previous lecture ll operte by the sme principle: they pproximte the
More information1.9 C 2 inner variations
46 CHAPTER 1. INDIRECT METHODS 1.9 C 2 inner vritions So fr, we hve restricted ttention to liner vritions. These re vritions of the form vx; ǫ = ux + ǫφx where φ is in some liner perturbtion clss P, for
More informationReview of Calculus, cont d
Jim Lmbers MAT 460 Fll Semester 2009-10 Lecture 3 Notes These notes correspond to Section 1.1 in the text. Review of Clculus, cont d Riemnn Sums nd the Definite Integrl There re mny cses in which some
More informationReview of basic calculus
Review of bsic clculus This brief review reclls some of the most importnt concepts, definitions, nd theorems from bsic clculus. It is not intended to tech bsic clculus from scrtch. If ny of the items below
More informationarxiv:gr-qc/ v1 17 Dec 1997
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
More informationChapter 4 Contravariance, Covariance, and Spacetime Diagrams
Chpter 4 Contrvrince, Covrince, nd Spcetime Digrms 4. The Components of Vector in Skewed Coordintes We hve seen in Chpter 3; figure 3.9, tht in order to show inertil motion tht is consistent with the Lorentz
More informationRecitation 3: More Applications of the Derivative
Mth 1c TA: Pdric Brtlett Recittion 3: More Applictions of the Derivtive Week 3 Cltech 2012 1 Rndom Question Question 1 A grph consists of the following: A set V of vertices. A set E of edges where ech
More informationA5682: Introduction to Cosmology Course Notes. 4. Cosmic Dynamics: The Friedmann Equation. = GM s
4. Cosmic Dynmics: The Friedmnn Eqution Reding: Chpter 4 Newtonin Derivtion of the Friedmnn Eqution Consider n isolted sphere of rdius R s nd mss M s, in uniform, isotropic expnsion (Hubble flow). The
More informationConservation Law. Chapter Goal. 5.2 Theory
Chpter 5 Conservtion Lw 5.1 Gol Our long term gol is to understnd how mny mthemticl models re derived. We study how certin quntity chnges with time in given region (sptil domin). We first derive the very
More information8 Laplace s Method and Local Limit Theorems
8 Lplce s Method nd Locl Limit Theorems 8. Fourier Anlysis in Higher DImensions Most of the theorems of Fourier nlysis tht we hve proved hve nturl generliztions to higher dimensions, nd these cn be proved
More informationName Solutions to Test 3 November 8, 2017
Nme Solutions to Test 3 November 8, 07 This test consists of three prts. Plese note tht in prts II nd III, you cn skip one question of those offered. Some possibly useful formuls cn be found below. Brrier
More informationg i fφdx dx = x i i=1 is a Hilbert space. We shall, henceforth, abuse notation and write g i f(x) = f
1. Appliction of functionl nlysis to PEs 1.1. Introduction. In this section we give little introduction to prtil differentil equtions. In prticulr we consider the problem u(x) = f(x) x, u(x) = x (1) where
More informationSection 6.1 INTRO to LAPLACE TRANSFORMS
Section 6. INTRO to LAPLACE TRANSFORMS Key terms: Improper Integrl; diverge, converge A A f(t)dt lim f(t)dt Piecewise Continuous Function; jump discontinuity Function of Exponentil Order Lplce Trnsform
More informationMath 8 Winter 2015 Applications of Integration
Mth 8 Winter 205 Applictions of Integrtion Here re few importnt pplictions of integrtion. The pplictions you my see on n exm in this course include only the Net Chnge Theorem (which is relly just the Fundmentl
More informationODE: Existence and Uniqueness of a Solution
Mth 22 Fll 213 Jerry Kzdn ODE: Existence nd Uniqueness of Solution The Fundmentl Theorem of Clculus tells us how to solve the ordinry differentil eqution (ODE) du = f(t) dt with initil condition u() =
More informationCHM Physical Chemistry I Chapter 1 - Supplementary Material
CHM 3410 - Physicl Chemistry I Chpter 1 - Supplementry Mteril For review of some bsic concepts in mth, see Atkins "Mthemticl Bckground 1 (pp 59-6), nd "Mthemticl Bckground " (pp 109-111). 1. Derivtion
More informationThe Regulated and Riemann Integrals
Chpter 1 The Regulted nd Riemnn Integrls 1.1 Introduction We will consider severl different pproches to defining the definite integrl f(x) dx of function f(x). These definitions will ll ssign the sme vlue
More informationEuler, Ioachimescu and the trapezium rule. G.J.O. Jameson (Math. Gazette 96 (2012), )
Euler, Iochimescu nd the trpezium rule G.J.O. Jmeson (Mth. Gzette 96 (0), 36 4) The following results were estblished in recent Gzette rticle [, Theorems, 3, 4]. Given > 0 nd 0 < s
More informationLecture 13 - Linking E, ϕ, and ρ
Lecture 13 - Linking E, ϕ, nd ρ A Puzzle... Inner-Surfce Chrge Density A positive point chrge q is locted off-center inside neutrl conducting sphericl shell. We know from Guss s lw tht the totl chrge on
More informationCMDA 4604: Intermediate Topics in Mathematical Modeling Lecture 19: Interpolation and Quadrature
CMDA 4604: Intermedite Topics in Mthemticl Modeling Lecture 19: Interpoltion nd Qudrture In this lecture we mke brief diversion into the res of interpoltion nd qudrture. Given function f C[, b], we sy
More informationNumerical Integration
Chpter 5 Numericl Integrtion Numericl integrtion is the study of how the numericl vlue of n integrl cn be found. Methods of function pproximtion discussed in Chpter??, i.e., function pproximtion vi the
More informationPhysics 116C Solution of inhomogeneous ordinary differential equations using Green s functions
Physics 6C Solution of inhomogeneous ordinry differentil equtions using Green s functions Peter Young November 5, 29 Homogeneous Equtions We hve studied, especilly in long HW problem, second order liner
More informationProperties of Integrals, Indefinite Integrals. Goals: Definition of the Definite Integral Integral Calculations using Antiderivatives
Block #6: Properties of Integrls, Indefinite Integrls Gols: Definition of the Definite Integrl Integrl Clcultions using Antiderivtives Properties of Integrls The Indefinite Integrl 1 Riemnn Sums - 1 Riemnn
More informationImproper Integrals, and Differential Equations
Improper Integrls, nd Differentil Equtions October 22, 204 5.3 Improper Integrls Previously, we discussed how integrls correspond to res. More specificlly, we sid tht for function f(x), the region creted
More informationand that at t = 0 the object is at position 5. Find the position of the object at t = 2.
7.2 The Fundmentl Theorem of Clculus 49 re mny, mny problems tht pper much different on the surfce but tht turn out to be the sme s these problems, in the sense tht when we try to pproimte solutions we
More informationChapter 5. , r = r 1 r 2 (1) µ = m 1 m 2. r, r 2 = R µ m 2. R(m 1 + m 2 ) + m 2 r = r 1. m 2. r = r 1. R + µ m 1
Tor Kjellsson Stockholm University Chpter 5 5. Strting with the following informtion: R = m r + m r m + m, r = r r we wnt to derive: µ = m m m + m r = R + µ m r, r = R µ m r 3 = µ m R + r, = µ m R r. 4
More informationLecture 6: Singular Integrals, Open Quadrature rules, and Gauss Quadrature
Lecture notes on Vritionl nd Approximte Methods in Applied Mthemtics - A Peirce UBC Lecture 6: Singulr Integrls, Open Qudrture rules, nd Guss Qudrture (Compiled 6 August 7) In this lecture we discuss the
More informationNUMERICAL INTEGRATION. The inverse process to differentiation in calculus is integration. Mathematically, integration is represented by.
NUMERICAL INTEGRATION 1 Introduction The inverse process to differentition in clculus is integrtion. Mthemticlly, integrtion is represented by f(x) dx which stnds for the integrl of the function f(x) with
More informationTheoretical foundations of Gaussian quadrature
Theoreticl foundtions of Gussin qudrture 1 Inner product vector spce Definition 1. A vector spce (or liner spce) is set V = {u, v, w,...} in which the following two opertions re defined: (A) Addition of
More informationTHE EXISTENCE-UNIQUENESS THEOREM FOR FIRST-ORDER DIFFERENTIAL EQUATIONS.
THE EXISTENCE-UNIQUENESS THEOREM FOR FIRST-ORDER DIFFERENTIAL EQUATIONS RADON ROSBOROUGH https://intuitiveexplntionscom/picrd-lindelof-theorem/ This document is proof of the existence-uniqueness theorem
More informationQuantum Mechanics Qualifying Exam - August 2016 Notes and Instructions
Quntum Mechnics Qulifying Exm - August 016 Notes nd Instructions There re 6 problems. Attempt them ll s prtil credit will be given. Write on only one side of the pper for your solutions. Write your lis
More informationarxiv: v2 [gr-qc] 7 Jan 2016
Derivtion of locl-in-time fourth post-newtonin ADM Hmiltonin for spinless compct binries Piotr Jrnowski Fculty of Physics, University of Bi lystok, Cio lkowskiego L, 5 45 Bi lystok, Polnd Gerhrd Schäfer
More informationNew Expansion and Infinite Series
Interntionl Mthemticl Forum, Vol. 9, 204, no. 22, 06-073 HIKARI Ltd, www.m-hikri.com http://dx.doi.org/0.2988/imf.204.4502 New Expnsion nd Infinite Series Diyun Zhng College of Computer Nnjing University
More informationThe Velocity Factor of an Insulated Two-Wire Transmission Line
The Velocity Fctor of n Insulted Two-Wire Trnsmission Line Problem Kirk T. McDonld Joseph Henry Lbortories, Princeton University, Princeton, NJ 08544 Mrch 7, 008 Estimte the velocity fctor F = v/c nd the
More informationJack Simons, Henry Eyring Scientist and Professor Chemistry Department University of Utah
1. Born-Oppenheimer pprox.- energy surfces 2. Men-field (Hrtree-Fock) theory- orbitls 3. Pros nd cons of HF- RHF, UHF 4. Beyond HF- why? 5. First, one usully does HF-how? 6. Bsis sets nd nottions 7. MPn,
More informationThe Properties of Stars
10/11/010 The Properties of Strs sses Using Newton s Lw of Grvity to Determine the ss of Celestil ody ny two prticles in the universe ttrct ech other with force tht is directly proportionl to the product
More informationThe Riemann-Lebesgue Lemma
Physics 215 Winter 218 The Riemnn-Lebesgue Lemm The Riemnn Lebesgue Lemm is one of the most importnt results of Fourier nlysis nd symptotic nlysis. It hs mny physics pplictions, especilly in studies of
More informationP 3 (x) = f(0) + f (0)x + f (0) 2. x 2 + f (0) . In the problem set, you are asked to show, in general, the n th order term is a n = f (n) (0)
1 Tylor polynomils In Section 3.5, we discussed how to pproximte function f(x) round point in terms of its first derivtive f (x) evluted t, tht is using the liner pproximtion f() + f ()(x ). We clled this
More informationA REVIEW OF CALCULUS CONCEPTS FOR JDEP 384H. Thomas Shores Department of Mathematics University of Nebraska Spring 2007
A REVIEW OF CALCULUS CONCEPTS FOR JDEP 384H Thoms Shores Deprtment of Mthemtics University of Nebrsk Spring 2007 Contents Rtes of Chnge nd Derivtives 1 Dierentils 4 Are nd Integrls 5 Multivrite Clculus
More informationPartial Derivatives. Limits. For a single variable function f (x), the limit lim
Limits Prtil Derivtives For single vrible function f (x), the limit lim x f (x) exists only if the right-hnd side limit equls to the left-hnd side limit, i.e., lim f (x) = lim f (x). x x + For two vribles
More informationVyacheslav Telnin. Search for New Numbers.
Vycheslv Telnin Serch for New Numbers. 1 CHAPTER I 2 I.1 Introduction. In 1984, in the first issue for tht yer of the Science nd Life mgzine, I red the rticle "Non-Stndrd Anlysis" by V. Uspensky, in which
More informationMatFys. Week 2, Nov , 2005, revised Nov. 23
MtFys Week 2, Nov. 21-27, 2005, revised Nov. 23 Lectures This week s lectures will be bsed on Ch.3 of the text book, VIA. Mondy Nov. 21 The fundmentls of the clculus of vritions in Eucliden spce nd its
More informationWe partition C into n small arcs by forming a partition of [a, b] by picking s i as follows: a = s 0 < s 1 < < s n = b.
Mth 255 - Vector lculus II Notes 4.2 Pth nd Line Integrls We begin with discussion of pth integrls (the book clls them sclr line integrls). We will do this for function of two vribles, but these ides cn
More informationDuality # Second iteration for HW problem. Recall our LP example problem we have been working on, in equality form, is given below.
Dulity #. Second itertion for HW problem Recll our LP emple problem we hve been working on, in equlity form, is given below.,,,, 8 m F which, when written in slightly different form, is 8 F Recll tht we
More informationPHYS 601 HW3 Solution
3.1 Norl force using Lgrnge ultiplier Using the center of the hoop s origin, we will describe the position of the prticle with conventionl polr coordintes. The Lgrngin is therefore L = 1 2 ṙ2 + 1 2 r2
More informationThe solutions of the single electron Hamiltonian were shown to be Bloch wave of the form: ( ) ( ) ikr
Lecture #1 Progrm 1. Bloch solutions. Reciprocl spce 3. Alternte derivtion of Bloch s theorem 4. Trnsforming the serch for egenfunctions nd eigenvlues from solving PDE to finding the e-vectors nd e-vlues
More information7.2 The Definite Integral
7.2 The Definite Integrl the definite integrl In the previous section, it ws found tht if function f is continuous nd nonnegtive, then the re under the grph of f on [, b] is given by F (b) F (), where
More information221B Lecture Notes WKB Method
Clssicl Limit B Lecture Notes WKB Method Hmilton Jcobi Eqution We strt from the Schrödinger eqution for single prticle in potentil i h t ψ x, t = [ ] h m + V x ψ x, t. We cn rewrite this eqution by using
More informationMath 32B Discussion Session Session 7 Notes August 28, 2018
Mth 32B iscussion ession ession 7 Notes August 28, 28 In tody s discussion we ll tlk bout surfce integrls both of sclr functions nd of vector fields nd we ll try to relte these to the mny other integrls
More informationINTRODUCTION. The three general approaches to the solution of kinetics problems are:
INTRODUCTION According to Newton s lw, prticle will ccelerte when it is subjected to unblnced forces. Kinetics is the study of the reltions between unblnced forces nd the resulting chnges in motion. The
More informationThe Active Universe. 1 Active Motion
The Active Universe Alexnder Glück, Helmuth Hüffel, Sš Ilijić, Gerld Kelnhofer Fculty of Physics, University of Vienn helmuth.hueffel@univie.c.t Deprtment of Physics, FER, University of Zgreb ss.ilijic@fer.hr
More informationThe Algebra (al-jabr) of Matrices
Section : Mtri lgebr nd Clculus Wshkewicz College of Engineering he lgebr (l-jbr) of Mtrices lgebr s brnch of mthemtics is much broder thn elementry lgebr ll of us studied in our high school dys. In sense
More informationMATH34032: Green s Functions, Integral Equations and the Calculus of Variations 1
MATH34032: Green s Functions, Integrl Equtions nd the Clculus of Vritions 1 Section 1 Function spces nd opertors Here we gives some brief detils nd definitions, prticulrly relting to opertors. For further
More informationImproper Integrals. Type I Improper Integrals How do we evaluate an integral such as
Improper Integrls Two different types of integrls cn qulify s improper. The first type of improper integrl (which we will refer to s Type I) involves evluting n integrl over n infinite region. In the grph
More informationA HELLY THEOREM FOR FUNCTIONS WITH VALUES IN METRIC SPACES. 1. Introduction
Ttr Mt. Mth. Publ. 44 (29), 159 168 DOI: 1.2478/v1127-9-56-z t m Mthemticl Publictions A HELLY THEOREM FOR FUNCTIONS WITH VALUES IN METRIC SPACES Miloslv Duchoň Peter Mličký ABSTRACT. We present Helly
More informationLecture XVII. Vector functions, vector and scalar fields Definition 1 A vector-valued function is a map associating vectors to real numbers, that is
Lecture XVII Abstrct We introduce the concepts of vector functions, sclr nd vector fields nd stress their relevnce in pplied sciences. We study curves in three-dimensionl Eucliden spce nd introduce the
More informationLecture 1. Functional series. Pointwise and uniform convergence.
1 Introduction. Lecture 1. Functionl series. Pointwise nd uniform convergence. In this course we study mongst other things Fourier series. The Fourier series for periodic function f(x) with period 2π is
More informationf(x) dx, If one of these two conditions is not met, we call the integral improper. Our usual definition for the value for the definite integral
Improper Integrls Every time tht we hve evluted definite integrl such s f(x) dx, we hve mde two implicit ssumptions bout the integrl:. The intervl [, b] is finite, nd. f(x) is continuous on [, b]. If one
More informationINTERNATIONAL CENTRE FOR THEORETICAL PHYSICS THE ALGEBRAIC APPROACH TO THE SCATTERING PROBLEM ABSTRACT
IC/69/7 INTERNAL REPORT (Limited distribution) INTERNATIONAL ATOMIC ENERGY AGENCY INTERNATIONAL CENTRE FOR THEORETICAL PHYSICS THE ALGEBRAIC APPROACH TO THE SCATTERING PROBLEM Lot. IXARQ * Institute of
More informationpotentials A z, F z TE z Modes We use the e j z z =0 we can simply say that the x dependence of E y (1)
3e. Introduction Lecture 3e Rectngulr wveguide So fr in rectngulr coordintes we hve delt with plne wves propgting in simple nd inhomogeneous medi. The power density of plne wve extends over ll spce. Therefore
More informationHigher Checklist (Unit 3) Higher Checklist (Unit 3) Vectors
Vectors Skill Achieved? Know tht sclr is quntity tht hs only size (no direction) Identify rel-life exmples of sclrs such s, temperture, mss, distnce, time, speed, energy nd electric chrge Know tht vector
More informationNote 16. Stokes theorem Differential Geometry, 2005
Note 16. Stokes theorem ifferentil Geometry, 2005 Stokes theorem is the centrl result in the theory of integrtion on mnifolds. It gives the reltion between exterior differentition (see Note 14) nd integrtion
More informationCalculus of Variations
Clculus of Vritions Com S 477/577 Notes) Yn-Bin Ji Dec 4, 2017 1 Introduction A functionl ssigns rel number to ech function or curve) in some clss. One might sy tht functionl is function of nother function
More informationSufficient condition on noise correlations for scalable quantum computing
Sufficient condition on noise correltions for sclble quntum computing John Presill, 2 Februry 202 Is quntum computing sclble? The ccurcy threshold theorem for quntum computtion estblishes tht sclbility
More informationPressure Wave Analysis of a Cylindrical Drum
Pressure Wve Anlysis of Cylindricl Drum Chris Clrk, Brin Anderson, Brin Thoms, nd Josh Symonds Deprtment of Mthemtics The University of Rochester, Rochester, NY 4627 (Dted: December, 24 In this pper, hypotheticl
More informationExam 2, Mathematics 4701, Section ETY6 6:05 pm 7:40 pm, March 31, 2016, IH-1105 Instructor: Attila Máté 1
Exm, Mthemtics 471, Section ETY6 6:5 pm 7:4 pm, Mrch 1, 16, IH-115 Instructor: Attil Máté 1 17 copies 1. ) Stte the usul sufficient condition for the fixed-point itertion to converge when solving the eqution
More informationSummary: Method of Separation of Variables
Physics 246 Electricity nd Mgnetism I, Fll 26, Lecture 22 1 Summry: Method of Seprtion of Vribles 1. Seprtion of Vribles in Crtesin Coordintes 2. Fourier Series Suggested Reding: Griffiths: Chpter 3, Section
More information221A Lecture Notes WKB Method
A Lecture Notes WKB Method Hmilton Jcobi Eqution We strt from the Schrödinger eqution for single prticle in potentil i h t ψ x, t = [ ] h m + V x ψ x, t. We cn rewrite this eqution by using ψ x, t = e
More informationAdvanced Calculus: MATH 410 Notes on Integrals and Integrability Professor David Levermore 17 October 2004
Advnced Clculus: MATH 410 Notes on Integrls nd Integrbility Professor Dvid Levermore 17 October 2004 1. Definite Integrls In this section we revisit the definite integrl tht you were introduced to when
More information1.2. Linear Variable Coefficient Equations. y + b "! = a y + b " Remark: The case b = 0 and a non-constant can be solved with the same idea as above.
1 12 Liner Vrible Coefficient Equtions Section Objective(s): Review: Constnt Coefficient Equtions Solving Vrible Coefficient Equtions The Integrting Fctor Method The Bernoulli Eqution 121 Review: Constnt
More informationAn approximation to the arithmetic-geometric mean. G.J.O. Jameson, Math. Gazette 98 (2014), 85 95
An pproximtion to the rithmetic-geometric men G.J.O. Jmeson, Mth. Gzette 98 (4), 85 95 Given positive numbers > b, consider the itertion given by =, b = b nd n+ = ( n + b n ), b n+ = ( n b n ) /. At ech
More informationTHE QUADRATIC RECIPROCITY LAW OF DUKE-HOPKINS. Circa 1870, G. Zolotarev observed that the Legendre symbol ( a p
THE QUADRATIC RECIPROCITY LAW OF DUKE-HOPKINS PETE L CLARK Circ 1870, Zolotrev observed tht the Legendre symbol ( p ) cn be interpreted s the sign of multipliction by viewed s permuttion of the set Z/pZ
More informationGravitational scattering of a quantum particle and the privileged coordinate system
rxiv:1712.02232v1 [physics.gen-ph] 5 Dec 2017 Grvittionl scttering of quntum prticle nd the privileged coordinte system A.I.Nishov July 20, 2018 Abstrct In grvittionl scttering the quntum prticle probes
More informationChapter 3. Vector Spaces
3.4 Liner Trnsformtions 1 Chpter 3. Vector Spces 3.4 Liner Trnsformtions Note. We hve lredy studied liner trnsformtions from R n into R m. Now we look t liner trnsformtions from one generl vector spce
More informationWHEN IS A FUNCTION NOT FLAT? 1. Introduction. {e 1 0, x = 0. f(x) =
WHEN IS A FUNCTION NOT FLAT? YIFEI PAN AND MEI WANG Abstrct. In this pper we prove unique continution property for vector vlued functions of one vrible stisfying certin differentil inequlity. Key words:
More informationPhysics 201 Lab 3: Measurement of Earth s local gravitational field I Data Acquisition and Preliminary Analysis Dr. Timothy C. Black Summer I, 2018
Physics 201 Lb 3: Mesurement of Erth s locl grvittionl field I Dt Acquisition nd Preliminry Anlysis Dr. Timothy C. Blck Summer I, 2018 Theoreticl Discussion Grvity is one of the four known fundmentl forces.
More informationLecture 14: Quadrature
Lecture 14: Qudrture This lecture is concerned with the evlution of integrls fx)dx 1) over finite intervl [, b] The integrnd fx) is ssumed to be rel-vlues nd smooth The pproximtion of n integrl by numericl
More informationNumerical integration
2 Numericl integrtion This is pge i Printer: Opque this 2. Introduction Numericl integrtion is problem tht is prt of mny problems in the economics nd econometrics literture. The orgniztion of this chpter
More informationMath 520 Final Exam Topic Outline Sections 1 3 (Xiao/Dumas/Liaw) Spring 2008
Mth 520 Finl Exm Topic Outline Sections 1 3 (Xio/Dums/Liw) Spring 2008 The finl exm will be held on Tuesdy, My 13, 2-5pm in 117 McMilln Wht will be covered The finl exm will cover the mteril from ll of
More informationSemigroup of generalized inverses of matrices
Semigroup of generlized inverses of mtrices Hnif Zekroui nd Sid Guedjib Abstrct. The pper is divided into two principl prts. In the first one, we give the set of generlized inverses of mtrix A structure
More informationIn-Class Problems 2 and 3: Projectile Motion Solutions. In-Class Problem 2: Throwing a Stone Down a Hill
MASSACHUSETTS INSTITUTE OF TECHNOLOGY Deprtment of Physics Physics 8T Fll Term 4 In-Clss Problems nd 3: Projectile Motion Solutions We would like ech group to pply the problem solving strtegy with the
More informationMAA 4212 Improper Integrals
Notes by Dvid Groisser, Copyright c 1995; revised 2002, 2009, 2014 MAA 4212 Improper Integrls The Riemnn integrl, while perfectly well-defined, is too restrictive for mny purposes; there re functions which
More informationReview of Gaussian Quadrature method
Review of Gussin Qudrture method Nsser M. Asi Spring 006 compiled on Sundy Decemer 1, 017 t 09:1 PM 1 The prolem To find numericl vlue for the integrl of rel vlued function of rel vrile over specific rnge
More informationBest Approximation in the 2-norm
Jim Lmbers MAT 77 Fll Semester 1-11 Lecture 1 Notes These notes correspond to Sections 9. nd 9.3 in the text. Best Approximtion in the -norm Suppose tht we wish to obtin function f n (x) tht is liner combintion
More informationRiemann Sums and Riemann Integrals
Riemnn Sums nd Riemnn Integrls Jmes K. Peterson Deprtment of Biologicl Sciences nd Deprtment of Mthemticl Sciences Clemson University August 26, 203 Outline Riemnn Sums Riemnn Integrls Properties Abstrct
More informationJackson 2.26 Homework Problem Solution Dr. Christopher S. Baird University of Massachusetts Lowell
Jckson 2.26 Homework Problem Solution Dr. Christopher S. Bird University of Msschusetts Lowell PROBLEM: The two-dimensionl region, ρ, φ β, is bounded by conducting surfces t φ =, ρ =, nd φ = β held t zero
More informationSOLUTIONS FOR ADMISSIONS TEST IN MATHEMATICS, COMPUTER SCIENCE AND JOINT SCHOOLS WEDNESDAY 5 NOVEMBER 2014
SOLUTIONS FOR ADMISSIONS TEST IN MATHEMATICS, COMPUTER SCIENCE AND JOINT SCHOOLS WEDNESDAY 5 NOVEMBER 014 Mrk Scheme: Ech prt of Question 1 is worth four mrks which re wrded solely for the correct nswer.
More information5.7 Improper Integrals
458 pplictions of definite integrls 5.7 Improper Integrls In Section 5.4, we computed the work required to lift pylod of mss m from the surfce of moon of mss nd rdius R to height H bove the surfce of the
More informationSpecial Relativity solved examples using an Electrical Analog Circuit
1-1-15 Specil Reltivity solved exmples using n Electricl Anlog Circuit Mourici Shchter mourici@gmil.com mourici@wll.co.il ISRAE, HOON 54-54855 Introduction In this pper, I develop simple nlog electricl
More informationDIRECT CURRENT CIRCUITS
DRECT CURRENT CUTS ELECTRC POWER Consider the circuit shown in the Figure where bttery is connected to resistor R. A positive chrge dq will gin potentil energy s it moves from point to point b through
More informationHomework # 4 Solution Key
PHYSICS 631: Generl Reltivity 1. 6.30 Homework # 4 Solution Key The metric for the surfce of cylindr of rdius, R (fixed), for coordintes z, φ ( ) 1 0 g µν = 0 R 2 In these coordintes ll derivtives with
More information1 1D heat and wave equations on a finite interval
1 1D het nd wve equtions on finite intervl In this section we consider generl method of seprtion of vribles nd its pplictions to solving het eqution nd wve eqution on finite intervl ( 1, 2. Since by trnsltion
More informationLecture 19: Continuous Least Squares Approximation
Lecture 19: Continuous Lest Squres Approximtion 33 Continuous lest squres pproximtion We begn 31 with the problem of pproximting some f C[, b] with polynomil p P n t the discrete points x, x 1,, x m for
More informationRecitation 3: Applications of the Derivative. 1 Higher-Order Derivatives and their Applications
Mth 1c TA: Pdric Brtlett Recittion 3: Applictions of the Derivtive Week 3 Cltech 013 1 Higher-Order Derivtives nd their Applictions Another thing we could wnt to do with the derivtive, motivted by wht
More informationIntroduction to the Calculus of Variations
Introduction to the Clculus of Vritions Jim Fischer Mrch 20, 1999 Abstrct This is self-contined pper which introduces fundmentl problem in the clculus of vritions, the problem of finding extreme vlues
More information