arxiv:astro-ph/ v1 5 Jul 2005

Transcription

1 Cosmologicl nonliner hydrodynmics with post-newtonin corrections Ji-chn Hwng, Hyerim Noh b, nd Dirk Puetzfeld c Deprtment of Astronomy nd Atmospheric Sciences, Kyungpook Ntionl University, Tegu, Kore b Kore Astronomy nd Spce Science Institute, Dejon, Kore c Deprtment of Physics nd Astronomy, Iow Stte University, Ames, IA, 500, USA rxiv:stro-ph/ v 5 Jul 2005 The purpose of this pper is to present generl reltivistic cosmologicl hydrodynmic equtions in Newtonin-like forms using the post-newtonin PN method. The PN pproximtion, bsed on the ssumptions of wek grvittionl fields nd slow motions, provides wy to estimte generl reltivistic effects in the fully nonliner evolution stge of the lrge-scle cosmic structures. We extend Chndrsekhr s first order PNPN hydrodynmics bsed on the Minkowski bckground to the Robertson-Wlker bckground. We ssume the presence of Friedmnn s cosmologicl spcetime s bckground. In the bckground we include the three-spce curvture, the cosmologicl constnt nd generl pressure; we show tht our PN pproch is successful only for sptilly flt cosmologicl bckground. In the Newtonin order nd PN order we include generl pressure, stress, nd flux. The Newtonin hydrodynmic equtions pper nturlly in the 0PN order. The sptil guge degree of freedom is fixed in unique mnner nd the bsic equtions re rrnged without tking the temporl guge condition. In this wy we cn conveniently try lterntive temporl guge conditions. We investigte number of temporl guge conditions under which ll the remining vribles re equivlently guge-invrint. Our im is to present the fully nonliner PN equtions in form suitble for implementtion in conventionl Newtonin hydrodynmic simultions with miniml extensions. The PN terms cn be considered s reltivistic corrections dded to the well known Newtonin equtions. The proper rrngement of the vribles nd equtions in combintion with suitble guge conditions would llow the possible future PN cosmologicl simultions to become more trctble. Our equtions nd guges re rrnged for tht purpose. We suggest wys of controlling the numericl ccurcy. The typicl PN order terms re bout times smller thn the Newtonin terms. However, we cnnot rule out possible presence of cumultive effects due to the time-delyed propgtion of the reltivistic grvittionl field with finite speed, in contrst to the Newtonin cse where chnges in the grvittionl field re felt instntneously. The quntittive estimtion of such effects is left for future numericl simultions. Contents I Introduction 2 II Bsic quntities 3 A Curvture B Energy-momentum tensor C Fluid-frme kinemtic quntities D Norml-frme kinemtic quntities E ADM quntities III Derivtions 0 A Equtions of motion B Einstein s equtions C ADM pproch IV PN Equtions 4 A Complete equtions to PN order B Idel fluid cse V Guge issue 7 A Guge trnsformtion B Chndrsekhr s guge C Uniform-expnsion guge D Trnsverse-sher guge

2 E Hrmonic guge F Trnsformtion between two guges VI Discussion 23 I. INTRODUCTION The nonliner evolution of lrge-scle cosmic structures is usully investigted within the frmework of Newtonin theory either in nlytic studies or numericl simultions. Without investigting generl reltivistic effects, however, it is not cler whether the Newtonin theory is sufficient to hndle such lrge-scle cosmologicl structures. If we regrd Einstein s grvity theory to be the correct frmework on such cosmic scles, the reltivistic effects should exist lwys. The point is to which level we cn prcticlly ignore reltivistic correction terms considering currently vilble levels of observtions nd experiments. In cse reltivistic or nonliner effects re not lrge, we hve two widely known wys to estimte reltivistic effects: one is the perturbtion pproch 3, nd the other one is the post-newtonin PN pproximtion 4 9. In the perturbtion study, we recently hve investigted the wekly nonliner regimes bsed on Einstein s grvity. We showed tht, except for the grvittionl wve coupling, the reltivistic sclr-type perturbtions coincide exctly with the Newtonin ones up to the second order 0. The pure reltivistic correction terms pper in the third order. These terms turn out to be independent of the horizon scle nd smll 0 5 order compred with the second-order Newtoninn/reltivistic terms. Our studies hve shown, from the view point of Einstein s grvity, tht Newtonin grvity is prcticlly relible ner the horizon scle where structures re supposed to be wekly nonliner. Therefore, justifying its use in current lrge-scle numericl simultions which cover the Hubble volume. Considering the ctiont--distnce nture of Newton s grvity our result is rther surprising one becuse the reltivistic perturbtion theory is pplicble on ll scles including the super-horizon scle. Thus, our perturbtion study ensures tht to the wekly nonliner order Newtonin theory is relible even ner nd beyond the horizon scle. Now, how bout sitution in the nonliner regime on scles smller thn the current horizon? On such scles the structures could be in fully nonliner stge. However, if the reltivistic grvity effect is smll we cn pply nother pproximtion scheme which is well developed to hndle isolted bodies in Einstein s theory: the PN pproximtion. The PN pproximtion hs been importnt to test Einstein s grvity. It presents the reltivistic equtions in form similr to the well known Newtonin equtions. In the PN method, by ssuming tht the reltivistic effects re smll, we expnd reltivistic corrections in powers of v/c. In nerly virilized systems we hve GM/Rc 2 v/c 2. Thus, the PN pproximtion is pplicble in the slow motion nd wek grvittionl field regime. The first order PN PN pproximtion corresponds to dding reltivistic correction terms of the order v/c 2 to the Newtonin order terms. In this pproch we recover the well known Newtonin hydrodynmics equtions in the 0PN order 5. The PN pproximtion is suitble to study systems in which Newtonin grvity hs dominnt role, while the reltivistic effects re smll but non-negligible. Well known pplictions include the precession of Mercury s perihelion, nd vrious solr system tests of Einstein s grvity theory like the light deflection 2. Recent successful pplictions include the genertion of grvittionl wves from compct binry objects, nd the wekly reltivistic evolution stges of isolted systems of celestil bodies 3. Another importnt ppliction is reltivistic celestil mechnics which is required by recent technology driven precise mesurements of the solr system bodies 4,5. Notice tht ll these pplictions re bsed on the PN pproximtion of isolted systems ssuming the Minkowski bckground spcetime. In this work on the cosmologicl PN pproximtion we ssume the presence of Robertson-Wlker bckground. Thus, we hve the Friedmnn world model built in s the bckground spcetime. We will show tht the Newtonin hydrodynmic equtions come out nturlly in the 0PN order which is the Newtonin limit. In the context of lrgescle cosmic structures PN effects could be essentilly importnt in hndling grvittionl wves nd grvittionl lensing effects. It is well known tht Newtonin order tretment of the grvittionl lensing shows only hlf of the result from Einstein s theory which is due to ignoring the PN correction in the metric. We lso nticipte tht PN correction terms could ffect the dynmic evolution of lrge-scle structures, especilly considering the ctiont--distnce nture of the Newtonin theory. This cn be compred with the reltivistic sitution in which the grvittionl field propgtes with finite speed. In this work we present the complete set of equtions to study PN effects in cosmologicl context. The PN pproch cn be compred with our previous studies of the wekly nonliner regime bsed on the perturbtion pproch which is pplicble in the fully reltivistic regime nd on ll scles. In comprison, lthough the PN equtions re pplicble in wek grvity regions inside the horizon, these re pplicble in fully nonliner sitution. Thus, the two pproches re complementry in enhncing our understnding of the reltivistic evolutionry spects of the lrge-scle structure in the universe. In the PN pproximtion we ttempt to study the equtions of motion nd the field equtions in Einstein s grvity in Newtonin wy s closely s possible. The Newtonin nd post-newtonin equtions of motion follow from 2

3 the energy-momentum conservtion. The Newtonin order potentil nd PN order metric vribles re determined by Einstein s equtions in terms of the Newtonin mtter nd potentil vribles. Thus, the metric or reltivistic contributions re reinterpreted s smll correction terms to the well known Newtonin hydrodynmic equtions. The Newtonin equtions nturlly follow in the 0PN order. The form of PN equtions is ffected by our choice of the guge conditions. In this work we tke unique sptil guge conditions which fix the sptil guge modes completely; under these sptil guge conditions the remining vribles re ll equivlently sptilly guge-invrint, see Sec. V. The temporl guge condition slicing condition will be deployed to hndle the mthemticl tretment of the equtions conveniently. We show how to choose the temporl guge condition which lso removes the temporl guge mode completely. Under ech of such temporl guge conditions the remining vribles re equivlently sptilly nd temporlly guge-invrint. We rrnge the equtions nd guges so tht the finl form of equtions is suitble for numericl implementtion in conventionl cosmologicl hydrodynmic simultions. We present bsic quntities nd the derivtion of our PN equtions in Sec. II nd III. If the reder is more interested in the ppliction of the PN equtions, the bsic set of cosmologicl PN equtions is summrized in Sec. IV. The guge issue is expounded in Sec. V. We discuss our results nd severl unresolved issues in Sec. VI. II. BASIC QUANTITIES A. Curvture As the metric we tke g 00 c 22U + 2U 2 c 4 4Φ +O 6, g 0i c 3P i +O 5, g ij 2 + c 22V γ ij +O 4, where x 0 ct, nd t is the cosmic scle fctor of the bckground Friedmnn world model. Indices,b,c,... indicte spcetime, nd i, j, k,... indicte spce. Tildes indicte spcetime covrint quntities, i.e., spcetime indices of quntities with tilde re rised nd lowered with the spcetime metric g b. The sptil index of P i is bsed on γ ij in rising nd lowering indices with γ ij, n inverse of γ ij. γ ij is the comoving time-independent sptil prt of the Robertson-Wlker metric, see Eq. 2 in 6. In flt Robertson-Wlker bckground γ ij becomes δ ij if we tke Crtesin coordintes. Compred with our previous nottions used in perturbtion studies in 7, we set the comoving prt of the three-spce bckground metric g 3 ij γ ij nd use prt of the Ltin indices to indicte the spce. We re following Chndrsekhr nd Nutku s nottion 5,7 extended to the cosmologicl sitution, see lso Fock 4. The 2U/c 2 term in g 00 gives the Newtonin limit, nd if we ignore ll the Newtonin nd post-newtonin correction terms we hve the Robertson-Wlker spcetime. Thus, our PN formultion is built on the cosmologicl bckground spcetime. O n indictes v/c n nd higher order terms tht we ignore. The expnsion in Eq. is vlid to PN order 5. Dimensions re s follows g b = g b =, γ ij = γ ij =, =, c = LT, U = V = c 2, P i = P i = c 3, Φ = c 4, 2 where L nd T indicte the length nd the time dimensions, respectively. In our metric convention in Eq. we hve ignored the possible presence of c 2 2C,i j +C i j +C j i like terms in g ij with C i i 0 by choosing the sptil C-guge conditions C 0 C i which remove the sptil guge mode completely; this will be explined in Sec. V; verticl br indictes the covrint derivtive bsed on γ ij. Under such sptil guge conditions, we cn regrd ll our remining PN vribles s equivlently sptilly guge-invrint ones, see Sec. V. We still hve freedom to tke the temporl guge condition which will be chosen lter depending on the mthemticl simplifiction or fesibility of physicl interprettion of the problem under considertion, see Sec. V: in the perturbtion pproch we termed this guge-redy pproch 8. We lso hve ignored the sptilly trcefree-trnsverse prt of the metric C ij with Ci i 0 Cj i j becuse grvittionl wves re known to show up only from the 2.5PN order 8. The inverse metric becomes 3

4 The determinnt of the metric tensor g is g 00 = + c 22U + 2U 2 c 4 +4Φ +O 6, g 0i = c 3 Pi +O 5, g ij = 2 c 22V γ ij +O 4. 3 g = 3 γ where γ is the determinnt of γ ij. The connection is Γ 0 00 = c 3 U + c 5 2 Φ+ Pi U,i +L O 7, + c 2 3V U+O 4, 4 Γ 0 0i = c 2U,i c 4 2Φ,i +ȧp i +L O 6, { Γ 0 ij = 2 ȧ c γ ij + V c 3 +2ȧ U +V γ ij + } P i j +L O 5, Γ i 00 = { 2 c 2U,i + 2U c 4 +VU,i 2Φ,i i } P +L O 6, Γ i 0j = ȧ c δi j + Vδ i c 3 j P i i j P j +L O 5, 2 Γ i jk = Γ γi jk + c 2 V,k δ i j +V,j δ i k V,i γ jk +L O 4, 5 where Γ γi jk is the connection bsed on γ ij; we introduce P i j 2 U,0 = U c t U. c The Riemnn curvture is R 0 00i = c 5 äp i + U,i jp +L j 2 O 7, R 0 0ij = L 2 O 6, R 0 i0j = äγij c 2 +U,i j + { 2 c 4 2U,i V,j U,i U,j +U,k V,k γ ij +2Φ,i j + P i j Pi j +P j i nd Pi j 2 V +2ȧ U +V γ ij + ȧ U +2ȧ 2 U +V γ ij } +L 2 O 6, R 0 ijk = 2 V c 32 + ȧ U γ ki,j P j ki + 2K γ ijp k +L 2 O 5, R i 00j = ä c 2 δi j + 2U,i j + c 4 { V + ȧ +U,i U,j +U,i V,j +U,j V,i U,k V,k δ i j 2Φ,i j R i 0jk = 2 V c 3 + ȧ U,j R i j0k = c 3 V,j δ i k V,i γ jk + ȧ δ i k P i j k U +2 V δ i j 2 2U +VU,i j +L 2 O 5, U,j δk i U,i γ jk + P i j 2 R i jkl = 2Kγ jkδ i l + c 22 ȧ 2 γ jk δ i l +V,j kδ i l V,i k γ lj P i j i P j k Pi j P j i. We hve } i +P j +L 2 O 6, +L 2 O 5, +L 2 O

5 It is convenient to hve P i jk Pi kj Rγi ljk Pl, P i jk = P i kj +R γl ijk P l, R γi jkl = K δk i γ jl δl i γ γ jk, R ij = 2Kγ ij, R γ = 6K, 7 where K indictes the comoving time-independent prt of bckground sptil curvture with dimension K = L 2. The Ricci curvture nd the sclr curvture become R 0 0 = c 2 3ä + + c { V 2U 4 3 V +3ȧ U +2 +6ä U } 2 U,i U V,i +2VU 2Φ P i i +L 2 O 6, R i 0 = 2 V c 3 + ȧ U +,i 2 R 0 i = c V + ȧ U R j i = 2K 2 δi j + ä +2ȧ2 c K 2 R = 6K 2 + ä c ȧ P j ji P i 2KP i +L 2 O 5,,i + P j i 2 j P i +2KP i +L 2 O 5, V δj i + 2 U V i j +L 2 O 4, 2U 4+3K 2 V +L 2 O 4, 8 where is Lplcin opertor bsed on γ ij. The Weyl curvture is introduced s C bcd R bcd δc 2 R bd δd R bc + g bd R c g bc R d + 6 Rδ c g bd δd g bc. 9 To PN order we hve C 0 00i = c 3 K P i +L 2 O 5, C 0 0ij = L 2 O 6, C 0 i0j = c 2 2 U +V,i j 6 U +Vγ ij +L 2 O 4, C 0 ijk = P l c 3 lk 2 P k 2KP k γ ji 2P j ki C i 00j = c 2 2 C i 0jk = c 3 2 C i j0k = c U +V i j 6 U +Vδi j P l lk P k +2KP k P l i l P i +2KP i γ jk +L 2 O 4, δ i j 2P i j k +L 2 O 5, +L 2 O 5, P l lj P j +2KP j δ i k +2 C i jkl = c U +Vδi k γ lj +U +V,j k δ i l U +V,i k γ lj P i j i P j +L 2 O 5, k +L 2 O 4. 0 We cn check the trcefree nture of the Weyl tensor: Cb bcd = 0 = C c bcd. Here we encounter non-vnishing trces involving K terms, like wheres C b b0i = c 3 K P i +L 2 O 5, C c 0ci = c 3 2K P i +L 2 O 5, C c 0c0 = L 2 O 4, Cb bij = L 2 O 4, Cc ic0 = L 2 O 5, Cc icj = L 2 O 4. 2 This indictes tht the presence of K terms in the PN order ppers to be not relible. Lter this point will be resolved s we show tht the K terms indeed cn be relted to the higher order terms in the PN expnsion, see below Eq. 88. Until we rech such conclusion bout the effect of K term we will keep it in our equtions. 5

6 B. Energy-momentum tensor The normlized fluid four-vector ũ with ũ ũ gives, to PN order, ũ 0 = + c 2 2 v2 +U + c v4 +v + 2 U +V 2 U2 +2Φ v i P i +O 6, ũ i c vi ũ 0, { ũ 0 = + c 2 2 v2 U { ũ i = c v i + c 3 + c v4 +v 2 U +V + } 2 U2 2Φ +O 6, } v i 2 v2 +U +2V P i +O 5, 3 where the index of v i is bsed on γ ij ; v i corresponds to the peculir velocity field. The energy-momentum tensor is decomposed into fluid quntities s follows T b = c + 2 c Π 2 ũ ũ b + pũ ũ b + g b + q ũ b + q b ũ + π b. 4 The energy density µ c 2 + Π is decomposed into the mteril energy density c 2 nd the internl energy density Π; p, q, nd π b re the isotropic pressure, the flux, nd the nisotropic stress, respectively. We hve q ũ 0, π b ũ b 0, π c c 0, nd π b π b, thus We introduce q 0 = c q iv i, π 0i = c π ijv j, π 00 = c 2 2 π ijv i v j. 5 q i c Q i, π ij 2 Π ij, 6 where indices of Q i nd Π ij re bsed on γ ij. We tke q nd π b to hve post-newtonin orders s introduced in Eq. 6, see for exmple 9. This will be justified by the energy conservtion eqution in the Newtonin context which will be derived lter, see Eq. 02. In strictly single component sitution we cn lwys remove Q i by following the fluid element. However, there exist situtions where we hve the dditionl flux terms present even when we follow the fluid elements; the fundmentl origin of such flux terms cn be trced to the presence of dditionl components in the energy-momentum. Thus, in our cse it is more generl nd convenient to keep the flux terms seprtely. Up to the PN order the condition π c c 0 gives We lso set Dimensions re s follows Π i i = c 2Π ijv i v j. 7, Π Π, p p. 8 ũ = ũ =, T b = T b = T b = p = q = q = π b = c 2 = ML T 2, v i = v i = c, Π = c 2, Q i = Q i = c q = c 3, Π ij = Π i j = Πij = π b = c 2. 9 We hve T 00 = c 2 { + + c 2 v 2 2U +Π + c 4 2Qi v i +Π ij v i v j +O 6 }, v 2 v 2 +2V +Π+ p 2UΠ+2U 2 4Φ 6

7 { T 0i = c 2 c v i + c 3 v i v 2 +2V +Π+ p P i + Qi +Π ij v j } +O 5, { T ij = 2 c 2 v i v j + p γij + Πij c 2 + } c 4 v i v j v 2 +2U +4V +Π+ p +2V p γij 2v i P j + 2 Qiv j +O 6, T = c 2 + c 2 Π 3 p +O From Eq. 4 covrintly we hve T = c + 2 Π c +3 p. 2 C. Fluid-frme kinemtic quntities The projection tensor h b bsed on the fluid-frme four-vector ũ is defined s h b g b +ũ ũ b. 2 The kinemtic quntities re defined s 20 θ b h c h d bũc;d, θ ũ ;, σ b θ b b, ω b 3 θ h θ b, ã ũ ũ ;b ũ b, 22 where we hve ũ θ b 0, ũ σ b 0, ũ ω b 0, ũ ã 0, nd θ θ. The kinemtic quntities θ, ã, σ b, nd ω b re the expnsion sclr, the ccelertion vector, the sher tensor, nd the vorticity tensor, respectively. To PN order, the projection tensor becomes h 0 0 = c 2v2 c 4 v 2 v 2 +2U +2V v i P i +O 6, { h 0 i = c v i + vi v 2 c 3 +2U +2V } P i +O 5, h i 0 = { c vi + c 2v2 + v 2 c 4 v 2 +2U +2V v j } P j +O 7, h i j = δi j + c 2vi v j + c 4 v i v j v 2 +2U +2V v i P j +O The kinemtic quntities become { θ ij = c v i j +ȧγ ij + c 3 v j v i +v i j 2 v2 +U +2V +v k v i v k j +v j v i k } + V + ȧ U +2V+ vk V,k + ȧ γ 2 v2 ij +2v i U +V,j P i j +L O 5, 24 θ = 3ȧ c + vi i + c 3 3 V +3ȧ U + 3 V,iv i + Uvi i v 2 + v 2 v i +L O 5, 25 2 i { σ ij = v c i j 3 vk k γ ij + c 3 v i vj +vi j 2 v2 +U +2V +v k v k i v j +v i v j k v i v j ȧ+ 3 vk k 3 γ ij 2 v 2 +U +2Vv k k + v 2 v k } +L O 5, 26 2 k ω ij = c v i j + {v c 3 j vi +vi j 2 v2 +U +2V +v k } v k j v i +v j v i k +2vi U +V,j P i j +L O 5, 27 7

8 ã i = c 2U,i + c 2v2 c 42Φ,i + + c 2 2 v2 +U t + v + v 2 cũi c 4 V,i +v j P j i +L O 6. From ũ c ã c 0 nd ũ b θ b 0 we hve ã 0 = c vi ã i, θ0i = c vj θij, θ00 = c 2 2vi v j θij, 29 nd similrly for σ 0 nd ω 0. The electric nd the mgnetic prts of the Weyl conforml tensor re introduced s 28 Ẽ b C cbd ũ c ũ d, Hb 2 η ef c C efbd ũ c ũ d, 30 where η bcd g ǫ bcd is the totlly ntisymmetric tensor density with ǫ bcd, the totlly ntisymmetric Levi-Civit symbol with ǫ 023. Ẽ b nd H b re symmetric, trcefree nd ũ Ẽ b = 0 = ũ Hb. To PN order we hve Ẽj i = C i 00j = c 2 2 H i j = 2ũ0 ũ 0 η i k 0l = { c 3 2 3ηikl 2 2 U +V i j 6 U +Vδi j +L 2 O 4, 3 C l k0j c vl C0 k0j + η i lk 0 c vl Ck 00j P m ml P l +2KP l + } 3 v lu +V γ kj +P l kj v l U +V kj +L 2 O 5, 32 where we introduced η ijk γ ǫ ijk with ǫ ijk, the totlly ntisymmetric Levi-Civit symbol with ǫ 23. Thus η 0ijk = γ g ηijk. Indices of η ijk re bsed on γ ij. Ẽ 0 follows from ũ b Ẽ b = 0 which gives Ẽ0 = ũ i /ũ 0 Ẽi = v i /cẽi, thus Ẽ0i L 2 O 3 nd Ẽ00 L 2 O 4. For nonvnishing K, Hb is not trcefree to L 2 O 3, but s we mentioned erlier the term involving K is lredy O 2 order higher, thus consistent. D. Norml-frme kinemtic quntities The normlized norml-frme four-vector ñ, which is norml to the spce-like hypersurfces, is defined s ñ i 0 with ñ ñ. To PN order we hve ñ 0 = + c 2U + c 4 2 U2 +2Φ +O 6, ñ i = c 3 Pi +O 5, ñ 0 = + c 2U c 4 2 U2 2Φ +O 6, ñ i By setting ũ i 0 in Eq. 3 we lso recover the norml vector, thus v i = c 2P i +LT O 4, 34 where v i is, sy, the velocity of the norml frme vector. The projection tensor bsed on ñ becomes h ij = g ij, h0i = g 0i, h00 = O 6, hi j = δ i j, hi 0 = c 3 Pi +O 5, h0 i = 0 = h Kinemtic quntities bsed on ñ become θ = c 3ȧ + c 3 3 V +3ȧ σ ij = c 3 P i j 3 Pk k γ ij U + Pi i +L O 5, 36 +L O 5, 37 ω ij = 0, 38 ã i = c 2U,i c 42Φ,i +L O

9 σ 0c follows from σ c ñ c 0, thus σ 0i L O 6 nd σ 00 L O 9 ; we hve ω b = 0. Similrly, ã 0 follows from ã c ñ c 0, thus ã 0 L O 5. The electric nd the mgnetic prts of Weyl tensor bsed on ñ give Ẽj i = C i 00j = c U +V i j 6 U +Vδi j +L 2 O 4, 40 H i j = 2ñ0ñ 0 η 0ikl C0 jkl = c 3 2 3ηikl P m ml 2 P l +2KP l γ jk +P l kj +L 2 O 5. 4 Thus, Ẽi j is the sme in both frmes to the PN order; compred with Eq. 32 H j i differs. The nonvnishing K term which cuses H b to be not trcefree to L 2 O 3 cn be ignored becuse the K term is of the O 2 order. E. ADM quntities In the Arnowitt-Deser-Misner ADM nottion the metric nd fluid quntities re 2 g 00 N 2 +N i N i, g 0i N i, g ij h ij, 42 ñ 0 N, ñ i 0, ñ 0 = N, ñ i = N N i, 43 E ñ ñ b Tb, J i ñ b Tb i, S ij T ij, S h ij S ij, Sij S ij 3 h ijs, 44 where N i, J i nd S ij re bsed on h ij s the metric, nd h ij is the inverse metric of h ij. The extrinsic curvture is K ij 2N N i:j +N j:i h ij,0, K h ij K ij, Kij K ij 3 h ij K, 45 where indices of K ij re bsed on h ij ; colon : denotes the covrint derivtive bsed on h ij with Γ hi jk 2 hil h jl,k +h lk,j h jk,l. The intrinsic curvtures re bsed on the metric h ij N i = c 3 Pi R hi jkl Γhi jl,k Γhi jk,l +Γhm jl Γhi km Γhm jk Γhi lm, R h ij R hk ikj, Rh h ij R h ij, Rh ij R h ij 3 h ijr h. 46 A complete set of ADM equtions will be presented in Sec. IIIC. To the PN order we hve N = c 2U + c 4 2 U2 2Φ +O 6, c 2U +O 7, N i = c 3P i + c 2 2V U +O 7, h ij = 2 + c 22V γ ij +O 4, h ij = 2 c 22V γ ij +O 4, 47 Γ hi jk = Γγi jk + V,k c 2 δj i +V,j δk i V,i γ jk +L O 4, 48 R h ij = R γ ij V,i j c 2 +γ ij V +L 2 O 4, R h = 2 6K +L c 24+3KV 2 O 4, 49 K ij = c ȧγ ij { V c 32 + ȧ U +2V γ ij + } P i j +L O 5, K = c 3ȧ c 3 3 V + ȧ U + Pi i +L O 5, 50 9

10 { E = c 2 + v 2 c 2 +Π + c 4 v 2 v 2 +2U +2V +Π+ p } +2UΠ 2P i v i + 2Qi v i +Π ij v i v j +O 6, { J i = c 2 c v i + c 3 v i v 2 +U +2V +Π+ p P i + Qi +Π ij v j } +O 5, { S ij = 2 c 2 v i v j + p γij + Πij c 2 + } c 4 v i v j v 2 +2U +4V +Π+ p +2 p Vγ ij 2v i P j + 2 Qiv j +O 6, { S = c 2 c 2 v 2 +3 p + c v 2 4 v 2 +2U +2V +Π+ p 2P i v i + 2Qi v i +Π ij v i v j } +O 6. 5 III. DERIVATIONS A. Equtions of motion Using Eq. 20 the energy nd momentum conservtion equtions give 0 = c T b 0;b = 3 3 σ + σv i + Q i c 2 +Π i j vj + i c 2 V + vi V U,i + ȧ v2 ṗ + T O 4, 52 0 = T i;b b = { 4 4 σv i + j} Qi c 2 +Π ij v + {σv i v j +Π ji + c } 2 Q j v i +Q i v j 2U +VΠ j i c σu 22U,i + c 2 v2 V,i j + p,i + {v c 2 i t + v U +3V+ 2 U +VU,i Φ,i P i 2 vj P i j + } U +3V,j Πj i + LT 2 O 4, 53 where σ + c 2 v 2 +2V +Π+ p. 54 For Q i = 0 = Π ij, V = U, = nd γ ij = δ ij these equtions reduce to Eqs. 64, 67 in Chndrsekhr 5. Eqution 52 cn be written in nother form s 3 + v i 3 i + c 2 Q i i +Πi jv j i + t + v Π+ 3ȧ + v p + T O 4 = 0, 55 where g 3 γũ0 + c 2 2 v2 +3V +O 4, 56 with 3 γ the bckground prt of g. This corresponds to Eq. 7 in Chndrsekhr 5; see lso Eq. 44 in 6. The mss conservtion continuity eqution 0 = ũ c ;c = + θ gives 0

11 0 = c ũ c ;c = c gũ c g,c = 3 + v i 3 i + T O 4 = 3 + i v 3 i + 2 c t + v 2 v2 +3V + T O Thus, if we ssume the mss conservtion, Eq. 55 gives t + v Π+ 3ȧ + p v = Q i i +Πi j vj i +T O The specific entropy is introduced s Td S 2 = d Π/c + p/c 2 d/, thus long the flow we hve T S = Π+ p/ c 2 = c t 3 + v Π+ 3ȧ + p v +T O This shows tht for n dibtic fluid flow the LHS of Eq. 58 vnishes; this lso mkes the RHS to vnish, which is nturlly required for n idel fluid. According to Chndrsekhr6: the conservtion of mss nd the conservtion of entropy re not independent requirements in the frmework of generl reltivity. And the reson for their independence in the Newtonin limit is tht in this limit c 2 our Eq. 57 reduces simply to the eqution of continuity. In 22 Blnchet, et l. suggested to use the following new combintion of vrible insted of v i vi g c 4 γ T i 0, 60 which becomes vi = v i + c 2 v i 2 v2 +U +2V +Π+ p P i + Qi +Π ij v j +LT O 4. 6 Notice tht v i is directly relted to the ADM flux vector in Eq. 5 Eqution T i;b b = 0 cn be written s follows g g T0 c γ i + Tj γ i = g γ J i = cv i. 62 g γ Tb g b i = 2 j {U,i + c U 2,i v 2 +Π +V,i v 2 +3 p v j P j i +2Φ,i +L 2 T 2 O 4 }, 63 nd we hve g Tj γ i = 3 v j v g i + pδ j γ i +Πj i + c 2 Q j v i 2VΠ j i Π ikv j v k +LT O Thus, Eq. 63 cn be rrnged in the following form v i + v i j vj = {U,i + c 2 U,i 2 v2 U +Π +V,i v 2 +3 p } v j P j i +2Φ,i {+ c pδ 2 3V U ji +Πji + c 2 Q j v i 2VΠ j i Π ikv j v k} j + c 2 v i Q + j +Π jk vk j t + v Π+ 3ȧ + v p +LT 2 O 4, 65 where the lst line vnishes in n dibtic idel fluid.

12 To the O 0 order, Eqs. 52, 53 give the Newtonin mss nd momentum conservtion equtions, respectively i i v = 0, v i + vj j v i + i p+ j Π j i iu = Thus, we nturlly hve the Newtonin equtions to 0PN order, e.g., see 23. This cn be compred with the sitution in perturbtion pproch where the Newtonin correspondence cn be chieved only fter suitble choice of different guges for different vribles, thus being non-trivil result even to the liner order in perturbtion pproch 3,24. Equtions 66, 67, s well s ll the PN equtions in this section, re fully nonliner. In the Friedmnn bckground we set the PN vribles U, Φ, P i, V, v i, Q i, nd Π ij equl to zero, nd set = b, Π = Π b, nd p = p b, with µ c 2 +Π nd σ b b + c Π 2 b +p b / b. Equtions 52, 57 give By subtrcting the bckground prt, Eq. 52 becomes 3 3 σ σ b + σv i + Q i c 2 +Π i jv j + i c 2 µ b +3ȧ µ b +p b = 0, 68 b +3ȧ b = V + vi V U,i + ȧ v2 ṗ ṗ b + T O 4 = We tke Einstein s equtions in the form Dimensions re s follows R b = 8πG c 4 B. Einstein s equtions T b 2 Tδ b +Λδb. 7 g b =, R b = R b = R b = L 2, T b = T = c 2 = ML T 2, Λ = L 2, G = T 2, 72 where M indictes the dimension of mss. To PN order, Eqs. 8, 20 give Λ+ { c 2 3ä + 2U +4πG + c V 4 3 V +3ȧ U +2 +6ä U U,i 2 U V,i +2VU 2Φ P i i +8πG v Π+ 3 } p = 0, 73 2 c 3 2 V + ȧ U + P j ji,i 2 P i 2KP i 8πG v i = 0, 74 2K 2 Λ δj i + ä +2ȧ2 c K 2 V 4πG δj i + 2 U V i j = 48πG v i v j + Π p δj i c 2 + Πi j, 75 where Eqs. 73, 74, nd 75 re R 0 0, R i 0, nd R j i prts of Eq. 7, respectively. We kept the O 4 term on the RHS of Eq. 75 in order to get the correct Friedmnn bckground eqution. To be consistent, we lso need to keep the LHS to O 4 which will explicitly involve 2PN order vribles. But we do not need such efforts for our bckground subtrction process becuse ll the O 4 order terms involve PN order vribles, 2

13 which do not ffect the Friedmnn bckground. To get the Friedmnn bckground we set U, Φ, P i, V, v i, Q i, nd Π ij equl to zero, nd set = b, Π = Π b, nd p = p b. Then, Eqs. 73, 75 give {ä 3 c 2 + 4πG b + Π b 3 c 2 + 3p } b Λc2 c 2 = 0, 76 3 {ä +2ȧ2 c 2 2 4πG b + Π b c 2 p } b + 2Kc2 c 2 2 Λc 2 δj i = Thus, we hve ä = 4πG b + Π b 3 c 2 ȧ 2 2 = 8πG 3 b + Π b c 2 + 3p b c 2 + Λc2 3, 78 Kc2 2 + Λc2 3, 79 which re the Friedmnn equtions. By subtrcting the bckground equtions, Eqs give { c 2 2U +4πG b + c V 4 3 V +3ȧ U +2 +6ä U 2 U,i U V,i +2VU 2Φ P i i +8πG v Π bπ b + 3 b } 2 p p = 0, 80 c 3 2 V + ȧ c 2 { U + P j ji,i 2 P i 2KP i 8πG v i = 0, 8 +4K V +4πG b δj i + } 2 U V i j = To O 2 order, Eq. 80 gives 2U = 4πG b, 83 which is Poisson s eqution in Newton s grvity. Notice tht in the Newtonin limit of our PN pproximtion the homogeneous prt of density distribution is subtrcted. This revised form of Poisson s eqution is consistent with Newton s grvity, nd in fct, is n improved form voiding Jens swindle 25. The decomposition of Eq. 82 into trce nd trcefree prts gives Thus, compred with Eq. 83, for K = 0 we hve +3K 2 V = 4πG b, 84 U V i j = KVδi j. 85 V = U, 86 where we ignored ny surfce term S with S 0 nd S,i j 0. Notice tht we hve V = U even in the presence of nisotropic stress; this differs from the sitution in the perturbtion pproch where U nd V in the zero-sher guge setting P i i = 0, see lter re different in the presence of the nisotropic stress 3. For lter use we tke divergence of Eq. 8 which gives V 2 + ȧ U = 4πG i v i + K 3Pi i. 87 Here, we reched point to ddress the issue relted to the bckground curvture in our PN pproch. By tking divergence of Eq. 8 nd using Eqs. 83, 84, nd 66 we hve 3

14 c 3 6K V 2K Pi i = 0, 88 which is pprently inconsistent. This is becuse the K term in Eq. 79 is relted to the c 2 higher order terms in the bckground Friedmnn eqution. We cn check this point by expnding the eqution to 2PN order; in tht order we cn find tht the bove term cn be removed by subtrcting the bckground eqution. Thus, it looks like our PN pproximtion is properly pplicble only for K = 0. In this sense, our K terms do not hve significnce becuse these re relted to O 2 order higher terms through the bckground eqution. At the moment, it is uncler whether this limittion of our PN pproch to nerly flt bckground is due to our subtrction process of the bckground equtions which spred over different PN orders, or intrinsic to the whole PN pproximtion. Menwhile, recent observtions fvour the flt bckground world model with non-vnishing cosmologicl constnt. We include the cosmologicl constnt in our PN pproximtion which ppers only in the bckground Friedmnn equtions. Ignoring K term we simply hve V = U. The ADM equtions re 2,3,7 C. ADM pproch R h = K ij Kij 2 3 K 2 + 6πG c 4 E +2Λ, 89 K j i:j 2 3 K,i = 8πG c 4 J i, 90 K,0 N K,i N i N +N :i i N K ij Kij 3 K2 4πG c 4 E +S+Λ = 0, 9 K j,0n i K j:kn i k N + K jk N i:k N K kn i k :jn = K K j N i :i j 3 δi j N:k k N + R hi j 8πG c 4 Si j, 92 E,0 N E,i N i N E K + 3 S S ij Kij +N 2 N 2 J i = 0, 93 :i J i,0 N J i:j N j N J j N j :i N KJ i +EN,i N +S j i:j +Sj i N,jN = Using the ADM quntities presented in Sec. IIE we cn show tht Eqs give the sme equtions we lredy hve derived from Einstein s equtions nd the energy nd momentum conservtions. Eqution 89 gives Eqs. 79, 84. Eqution 90 gives Eq. 8. Eqution 9 gives Eq. 73. Eqution 92 gives the trcefree prt of Eq. 82. Eqution 93, 94 give Eqs. 52, 53, respectively. The ADM equtions re often used in the PN pproch 22,28,3,26. We present the ADM pproch becuse the ADM equtions show the fully nonliner structure of Einstein s grvity in form suitble for numericl tretment. In fct, Eq. 65 cn be derived from Eq. 94 using the identifiction mde in Eq. 62. IV. PN EQUATIONS A. Complete equtions to PN order Here we summrize the complete set of equtions vlid to PN order in the cosmologicl sitution. The bckground vribles, b, Π b nd p b re provided by solving Eqs. 68, 69, 78, nd 79 ä = 4πG b + Π b 3 c 2 + 3p b + Λc2 c 2 3, 95 ȧ 2 2 = 8πG 3 b + Π b c 2 Kc2 2 + Λc2 3, 96 µ b +3ȧ µ b +p b = 0, 97 b +3ȧ b = 0, 98 4

15 where µ b bc 2 +Π b /c 2. In Eqs only two re independent, nd from Eqs. 97, 98 we hve b Πb + 3ȧ/p b = 0. To the Newtonin order we hve Eqs. 66, 67, nd i i v = 0, 99 v i + vj j v i + i p+ j Π j i iu = 0, 00 2U +4πG b = 0. 0 In the Newtonin cse energy conservtion nd mss conservtion provide n dditionl eqution which is Eq. 58 t + v Π+ 3ȧ + p v + Q i i +Πi jv j i = In the non-expnding bckground this eqution gives the well known energy conservtion eqution in the Newtonin theory, see for exmple Eq in 27. The Newtonin order hydrodynmic equtions re vlid in the presence of generl bckground curvture K. In the following we set K = 0 to the PN order, becuse our PN expnsion fils to pply in the presence of K term. Thus we hve V = U. The energy nd the momentum conservtion equtions vlid to PN order re derived in Eqs. 70, 53, nd these re 3 3 σ + σv i + Q i c 2 +Π i j vj + i c 2 U + ȧ v2 ṗ = 0, 03 { σv 4 i + c j} Qi 2 +Π ij v 4 + c 22U p,i + c 2 + σ + c 2 v2 σv i v j +Π ji + c 2 Q j v i +Q i v j 4UΠ j i U,i 4v i t + v U + 2 U 2 Φ,i P i 2 vj P i j + 4 U,jΠ j i j = 0, 04 where σ + c 2 v 2 +2U +Π+ p. 05 To the Newtonin order Eqs. 03, 04 for nd v i reduce to Eqs. 99, 00, respectively. From Eq. 55, 65 we cn derive lterntive forms 3 + v i 3 i = c 2 Q i i +Πi jv j i + t + v Π+ 3ȧ + v p, 06 v i + v i j vj = { U,i + 3 c 2 2 v2 U +Π+3 p } U,i +2Φ,i v j P j i {+ c pδ 22U ji +Πji + c 2 Q j v i 2VΠ j i Π ikv j v k} j + c 2 v i Q + j +Π jk vk j t + v Π+ 3ȧ + v p, 07 where + c 2 2 v2 +3U, v i v i + c 2 2 v2 +3U +Π+ p v i P i + Qi +Π ij v j. 08 5

16 The hydrodynmic nd thermodynmic vribles Π, p, Q i nd Π ij should be provided by specifying the equtions of stte nd the thermodynmic stte of the system under considertion. The metric vribles U, Φ nd P i cn be expressed in terms of the Newtonin vribles, v i, Π, p, Q i nd Π ij by using Einstein s equtions. Equtions 83, 84 give 2U = 4πG b, 09 which lredy determines U. Equtions 80, 8, nd 87 give { 2U +4πG b+ c 2 2 2Φ 2UU + P i i +3Ü +9ȧ U +6ä U +8πG v Π bπ b + 3 b } 2 p p = 0, 0 2P i = 6πG v i + Pj j +4 U +4ȧ U,,i U i 2 = 4πG v i + ȧ b. 2 Notice tht to PN order Einstein s equtions do not involve the nisotropic stress or flux term. To PN order Eq. 0 determines Φ. Our conservtion equtions 03, 04 contin U terms. In order to hndle these terms it ws suggested in 28 tht the Poisson-type eqution in Eq. 2 provides better numericl ccurcy. Notice tht, wheres only the sptil grdient of the potentil U ppers in the Newtonin limit, we hve bre U terms present in the PN order 29. In order to hndle these equtions we hve the freedom to tke one temporl guge condition. This corresponds to imposing condition on P i i or Φ. Guge relted issues will be ddressed in detil in Sec. V. There, we will show tht ll our vribles in this section re sptilly guge-invrint, nd re temporlly guge-redy. B. Idel fluid cse We consider n idel fluid, i.e., Q i 0 Π ij, nd ssume the dibtic condition Eq. 02 pplies. Thus, the internl energy is determined by the energy conservtion eqution t + v Π+ 3ȧ + p v = 0. 3 To the bckground order we hve Eqs To the Newtonin order we hve Eqs Equtions 03, 04 become 3 3 σ + σv i i + c 2 U + ȧ v2 ṗ = 0, σv i + σvi v j j σ + c 2 v2 U,i + c 22U p,i + 4v c 2 i t + v U + 2 U 2 Φ,i P i 2 vj P i j = 0, 5 where σ is given in Eq. 05. These equtions cn be written s 3 + i v 3 i + 2 c t + v 2 v2 +3U = 0, 6 v i + v i jv j U,i + p,i + c 2 v 2 +4U +Π+ p p,i +v i t + v 2 v2 +3U +Π+ p v2 U,i + 2 U 2 Φ,i P i 2 vj P i j =

17 Alterntive forms follow from Eqs. 06, v i 3 i = 0, v i + v i j vj = + c 22U p,i c 2 2 v2 U +Π+ p U,i + 2Φ,i c 2 v j P j i, 9 where + c 2 2 v2 +3U, vi v i + c 2 2 v2 +3U +Π+ p v i P i. 20 All these equivlent three sets of equtions re written in Eulerin forms. Einstein s equtions in 0-2 will provide the metric perturbtion vribles U, Φ nd P i in terms of the Newtonin fluid vribles, v i, Π, nd p, fter tking one temporl guge condition. The pressure term should be provided by n eqution of stte. In the cse of zero-pressure fluid with vnishing internl energy we cn set p = 0 = Π. We cn tke ny set of equtions of motion depending on the mthemticl convenience in numericl tretments. V. GAUGE ISSUE A. Guge trnsformtion We consider the following trnsformtion between two coordintes x nd ˆx ˆx x + ξ x e. 2 For tensor quntity we use the tensor trnsformtion property between x nd ˆx coordintes t b x e = ˆxc x ˆx d x bˆ t cd ˆx e. 22 Compring tensor quntities t the sme spcetime point, x, we cn derive the guge trnsformtion property of tensor quntity, see Eq. 226 in 7, ˆ t b x e = t b x e 2 t c ξc,b t b,c ξc +2 t c ξd,b ξc,d + t cd ξc, ξd,b + ξ d 2 ξ c, t bc,d +2 t c ξc,bd + 2 t b,cd ξc + t b,c ξc,d. 23 We considered ξ to the second perturbtionl order which will turn out to be sufficient for the PN order. In the following we will consider the guge trnsformtion properties of the metric nd the energy-momentum vribles. As the metric we consider the following more generlized form g 00 c 22U + 2U 2 c 4 4Φ +O 6, g 0i c 3P i +O 5, g ij 2 + c 22V γ ij + c 2 2C,i j +2C i j +2C ij +O 4, 24 where C i is trnsverse C i i 0, nd C ij is trnsverse nd trcefree C j i j 0 Ci i ; indices of C i nd C ij re bsed on γ ij. In Eq. 24 we introduced 0 independent metric components: U nd Φ together -component, P i 3-components, V -component, C -component, C i 2-components, nd C ij 2-components. It is known tht grvittionl wves show up in the 2.5PN order 8. Thus, we ignore the trnsverse nd trcefree prt, i.e., set C ij 0 to PN order. We wish to keep the metric in the form of Eq. 24 in ny coordinte system. Thus, we tke the trnsformtion vrible ξ to be PN-order quntity. We consider coordinte trnsformtions which stisfy 7

18 ξ 0 c ξ20 + c 3ξ , ξi c 2 ξ2i +..., 25 where the indices of the ξ 2i s re bsed on γ ij. The guge trnsformtion of Eq. 23 pplied to the metric gives Û = U + ξ 20, ˆΦ = Φ+ 2 ξ 40 2 Uξ 20 2 U,iξ 2i 2 ξ20 ξ20 4 ξ 20 ξ20 2 ξ 20,i = 0, ˆP i = P i ξ40,i + ξ2 i +3KĈ C = +3K ξ2i i +2KĈi = +2K C i ξ2 i +3K ˆV = +3K ξ2i ξ 20,i, 27 U + ξ 20 ξ 20,i + ξ 20,i ξ ξ 2j ξ 20,j + ξ 20,i 4 2 ξ 20,i 3 3 ξ2j ji 2 + { 4 i 3 +3Kξ2j j + V ȧ ξ20 ξ 20,i ξ20,j 4 2 ξ 20,i ξ20,j j 28,i, 29 ij, 30 } kj ξ 20,k ξ20,j, 3 +2K 4 2 ξ 20,i ξ 20,i + ij 4 2 ξ 20,i,j ξ Equtions 26, 27 follow from the trnsformtion of g 00 ; Eqs. 28, 29 follow from g 0i ; Eqs follow from g ij. Eqution 28 shows tht ξ 20 is sptilly constnt, After this, Eq. 30 gives ξ 20 = ξ 20 t. 33 Ĉ = C ξ2i i. 34 Thus, by choosing C 0 C 0 is enough s the guge condition this implies tht we set C = 0 to be vlid in ny coordinte, we hve ξ 2i i Imposing the conditions in Eqs. 33, 35, Eq. 3 gives = Ĉ i = C i ξ2 i. 36 Thus, by choosing C i 0 to be vlid in ny coordinte, i.e., choosing C i 0 s the guge condition, we hve ξ 2 i = Under the conditions in Eqs. 33, 37 the remining metric perturbtion vribles trnsform s Û = U + ξ 20, ˆV = V ȧ ξ20, ˆΦ = Φ+ 2 ξ 40 2 Uξ 20 2 ξ20 ξ20 4 ξ 20 ξ20, ˆP i = P i ξ40,i. 38 Notice tht for non-vnishing ξ 20, U nd V trnsform differently under the guge trnsformtion even in the flt cosmologicl bckground. Since the sptilly constnt ξ 20 cn be bsorbed by globl redefinition of the time coordinte, without losing generlity we cn set ξ 20 equl to zero 8

19 ξ 20 0, 39 thus llowing V = U in ny coordinte. If we set ξ 20 = 0 but do not tke the sptil guge which fixes ξ 2 i, from Eqs , the metric vribles trnsform s Û = U, 40 ˆΦ = Φ+ 2 ξ 40 2 U,iξ 2i, 4 ˆP i = P i ξ40,i + ξ2 i, 42 Ĉ = C ξ2i i, 43 +2KĈi = +2K C i ξ2 i 3 ξ2j ji i +3Kξ 2j j, 44 ˆV = V. 45 Tking the sptil guge conditions C = 0 = C i we hve ξ 2 i = 0 = ξ 20, nd the remining metric perturbtion vribles trnsform s Û = U, ˆV = V, ˆΦ = Φ+ 2 ξ 40, ˆPi = P i ξ40,i. 46 In the perturbtion nlysis we cll C 0 C i the sptil C-guge 7. Apprently, under these guge conditions the sptil guge trnsformtion function to PN order, ξ 2 i, is fixed completely, i.e., ξ 2 i = 0. By tking such guge conditions, the only remining guge trnsformtion function to PN order is ξ 40. This temporl guge trnsformtion function ffects only Φ nd P i. If we tke Φ 0 s the temporl guge condition, we hve ξ 40 = 0, thus ξ 40 does not vnish even fter imposing the temporl guge condition nd hs generl dependence on sptil coordinte ξ 40 x. Thus, in this guge, even fter imposing the temporl guge condition the temporl guge mode is not fixed completely. Wheres, if we tke P i i 0 s the temporl guge condition, we hve ξ40 = 0, thus the temporl guge condition is fixed completely. In fct, from the guge trnsformtion properties of Φ nd P i we cn mke the following combintions Φ,i + 2 P i, Φ+ 2 P i i, 47 which re invrint under the temporl guge trnsformtion, i.e., temporlly guge-invrint. Menwhile, U nd V re lredy temporlly guge-invrint. Since our sptil C-guge lso hs removed the sptil guge trnsformtion function completely, we cn correspond ech remining vrible to unique sptilly nd temporlly guge-invrint combintion. Thus, in this sense, we cn equivlently regrd ll our remining vribles s guge-invrint ones. For exmple, for K = 0, from Eqs we cn show ˆP i +Ĉi +Ĉ,i = Pi +C i +C,i ξ40,i. 48 Thus, P i + C i +C,i is sptilly guge-invrint nd becomes P i under the sptil C-guge. This implies tht P i under the sptil C-guge is equivlent to unique guge-invrint combintion P i + C i +C,i. Similrly, in the cse of the temporl guge, from Eq. 47 we cn show tht Φ under the P i i = 0 guge is the sme s unique guge-invrint combintion Φ + 2 Pi i. In this sense, under guge conditions which fix the guge mode completely, the remining vribles cn be regrded s equivlently guge-invrint ones with corresponding guge-invrint vribles. In Eq. we begn our PN nlysis by choosing the sptil C-guge C 0 C i. 49 By exmining Eqs we notice tht the sptil C-guge is most economic in fixing the sptil guge mode completely without due lterntive. In this sense the sptil C-guge cn be regrded s unique choice nd we do not lose ny mthemticl convenience by tking this sptil guge condition. In the literture these conditions re often expressed in the following forms. The sptil component of the hrmonic guge condition sets the PN prt of 9

20 equl to zero. We cn lso set the PN prt of g b Γi b = g g g ic,c = c 2 2 C,i +C i +L O 4, 50 g jk g ij,k 3 g jk,i = c 2 43 C,i +C i +L O 4, 5 equl to zero; notice tht in Eqs. 50, 5 we ssumed K = 0. In either cse we rrive t Eq. 49. We hve shown tht under the conditions in Eq. 49, the sptil guge trnsformtion is fixed completely without losing ny generlity or convenience. Wheres, we hve not chosen the temporl guge condition which cn be best chieved by imposing condition on P i i. We cll this the guge-redy strtegy 8,7. It is convenient to choose this remining temporl guge condition depending on the problem we encounter; or to try mny different ones in order to find out the best suitble one or ones. Now, let us consider the guge trnsformtion property of the energy-momentum tensor. We set ξ The guge trnsformtion property of Eq. 23 pplied to the energy-momentum tensor in Eq. 20 gives ˆ = + c 2 2, ˆΠ = Π,i ξ2i 2, ˆv i = v i + c 2v2 i, 54 ˆQ i = Q i v 2 i ξ2 i + v jξ 2j i + v i jξ 2j, 55 ˆp = p c 2 p,i ξ 2i pξ2i i Πi jξ 2j i, 56 ˆΠ ij = Π ij c 2 2p ξ 2 i j 3 γ ijξ 2k k +Π ij k ξ 2k +2Π ki ξ 2k j 2 3 γ ijπ k l ξ2l k. 57 Equtions 52, 53 follow from the trnsformtion property of T 00 ; Eqs. 54, 55 follow from T 0i ; Eqs. 56, 57 follow from T ij. The trnsformtion functions 2 nd v 2 i re not determined, see below. If we tke the sptil C-guge, i.e., set ξ 2 i 0, we hve ˆ = + c 2 2, ˆΠ = Π 2 Thus, we hve, ˆv i = v i + c 2v2 i, ˆQi = Q i v 2 i, ˆp = p+o 4, ˆΠij = Π ij +O ˆ + c ˆΠ 2 = + c 2Π, ˆv i + c 2 ˆQ i = v i + c 2 59 Qi. Although the guge trnsformtion properties of nd Π re not determined individully, the energy density µ +Π/c 2 is guge-invrint under our C-guge. For vnishing Π we hve 2 = 0. Similrly, the guge trnsformtion properties of v i nd Q i re not determined individully. For vnishing flux term in ny coordinte, Q i = 0, we hve v 2 i = 0, thus ˆv i = v i +O 4. Thus, the guge trnsformtion property of v i is determined only for vnishing flux term. From Eq. 08 we hve ˆ = + 2 2, ˆv i = vi + 2 c c 2 ξ40,i. 60 We cn check tht up to the PN order Einstein s equtions nd the energy nd momentum conservtion equtions in Sec. IV re invrint under the guge trnsformtion. In the following we introduce severl temporl guge conditions ech of which fixes the temporl guge mode completely. In reltivistic grvity, guge conditions consist of four non-tensoril reltions imposed on the metric tensor or the energy-momentum tensor. Our purpose is to employ the temporl guge slicing condition to mke the resulting equtions simple for mthemticl/numericl tretment. Depending on the sitution we cn lso tke lterntive sptil guge conditions for tht purpose. In the following we impose the sptil C-guge in Eq. 49, nd set K = 0. vi 20

21 B. Chndrsekhr s guge Chndrsekhr s temporl guge condition in his Eq. 24 of 5 corresponds to setting the PN prt of g g ij 0i,j 2 g ij,0 = c 3ȧ c 3 Pi i +3 U +L O 5, 6 equl to zero; in the literture this is often clled the stndrd PN guge 22. We tke Pi i +3 U +mȧ U = 0, 62 s Chndrsekhr s guge. We used the freedom to dd n rbitrry mȧ/u term with m, rel number. In this cse Eqs. 0, give 2P i = 6πG v i + U m 4 ȧ U, 63,i { 2U +4πG b + c 2 2 ȧ 2Φ m 3 U + 6 m ä mȧ2 2 U +8πG v Π bπ b +U b+ 3 b } 2 p p = Therefore, U, P i, nd Φ re determined by Eqs. 09, 63, nd 64. The vrible U cn be determined from Eq. 2. This completes our PN scheme bsed on Chndrsekhr s guge. When we hndle this complete set of PN equtions numericlly, we should monitor whether the chosen guge condition is stisfied lwys; this could be used to control the numericl ccurcy. C. Uniform-expnsion guge The expnsion sclr of the norml-frme vector, θ ñ c ;c, is given in Eq. 36. It is the sme s the trce of extrinsic curvture K with minus sign, see Eq. 50. Tking the PN prt of K equl to zero Pi i +3 U +3ȧ U 0, 65 cn be nturlly clled the uniform-expnsion guge. In the literture it is often clled the ADM guge 22, or the mximl slicing condition in numericl reltivity30. This condition corresponds to the m = 3 cse of Chndrsekhr s guge in Eq. 62. D. Trnsverse-sher guge The sher of the norml-frme vector is given in Eq. 37. Thus, the guge condition P i i 0, 66 cn be clled the trnsverse-sher guge. In this cse Eqs. 0, give 2P i = 6πG v i + 4 U + ȧ U, 67,i { 2U +4πG b + c 2 2 2Φ+3Ü +9ȧ U +6ä U +8πG v Π bπ b +U b+ 3 b } 2 p p = Therefore, U, P i, nd Φ re determined by Eqs. 09, 67, nd 68. U cn be determined from Eq. 2. We still hve Ü which cn be determined similrly by tking the time derivtive of Eq. 2 nd by using the bckground nd the Newtonin equtions in Eqs

22 E. Hrmonic guge We hve g b Γ0 b = g g 0 = g, c 3ȧ + c 3 Pi i +4 U +6ȧ U +L O The well known hrmonic guge condition sets the PN prt of Eq. 69 equl to zero, thus we tke Pi i +4 U +mȧ U 0, 70 where we used freedom to dd n rbitrry mȧ/u term. Combined with Eq. 50 the full hrmonic guge condition cn be expressed by setting the PN prts of g b Γc b equl to zero. In the perturbtion pproch the hrmonic guge often clled de Donder guge is known to be bd choice becuse the condition involves derivtives of the metric vribles; this leds to n incomplete guge fixing nd consequently we would hve to hndle higher order differentil equtions which is unnecessry, see the Appendix in 3. In this guge Eqs. 0, give 2P i = 6πG v i m 4 ȧ U,i, 7 { 2U +4πG b + c 2 2 ȧ 2Φ Ü m U + 6 m ä mȧ2 2 U +8πG v Π bπ b +U b+ 3 b } 2 p p = Therefore, U, P i, nd Φ re determined by Eqs. 09, 7, nd 72. The vrible U cn be determined from Eq. 2. We lso hve Ü term in Eq. 72 which might demnd for more involved numericl implementtion s explined below Eq. 68. However, its presence mkes the propgting nture of the PN order metric fluctutions pprent in this guge condition. This time-delyed propgtion of the grvittionl field in the PN pproximtion cn be compred with the ction-t--distnce nture of Poisson s eqution in the Newtonin order in Eq. 72. We nticipte tht the reltivistic time-delyed propgtion could led to seculr time cumultive effect. This could be importnt even in the cse in which the PN correction terms re smll compred with the Newtonin terms; in the lrge-scle clustering regions we would expect v/c 2 GM/Rc , see Sec. VI. Using g 00 +2U/c 2 we hve U U + c 2 U 2 +2Φ +c 2 O 4, 73 thus U U ;c c = 2U c 2 Ü+3ȧ U+2U 2U +T 2 O 4 = 2U c 2 2 U 2 2Φ +Ü +3ȧ U +2U 2U +T 2 O 4, 74 nd Eq. 72 cn be written s U+4πG b + c 2 { 2U2 m 4 ȧ U+ 6 m ä mȧ2 2 U+8πG v Π bπ b + 3 } 2 p p b = Thus, in the hrmonic guge condition the propgtion speed of the grvittionl potentil U is the sme s the speed of light c. When we mention the propgtion speed, we re considering the D Alembertin prt of the wve eqution ignoring the nonliner terms nd bckground expnsion. In this cse we hve U = 2 U c 2 Ü, nd for K = 0, becomes the ordinry Lplcin in flt spce. Apprently, the wve speed of the metric potentil cn tke n rbitrry vlue depending on the temporl guge condition we choose. For exmple, if we tke Pi i +n U +mȧ U 0, 76 22

23 s the guge condition, with n nd m rel numbers, Eqs. 0, give 2P i = 6πG v i n 4 U +m 4 ȧ U { 2U +4πG b + c 2 2 ȧ 2Φ n 3Ü 2n+m 9 U +, 77,i 6 m ä mȧ2 2 U +8πG v Π bπ b +U b+ 3 b } 2 p p = In this cse, for n 3, the speed of propgtion corresponds to c n 3, 79 which cn tke n rbitrry vlue depending on our choice of the vlue of n. It becomes c for n = 4 e.g., the hrmonic guge, nd infinity for n = 3 e.g., Chndrsekhr s guge nd the uniform-expnsion guge. In the cse of the trnsverse-sher guge we hve n = 0, thus Eq. 68 is no longer wve eqution. F. Trnsformtion between two guges Now, let us show how we cn relte the equtions nd solutions known in one guge condition to the ones in ny other guge condition. Our sptil C-guge condition lredy fixed the guge trnsformtion function ξ 2 i 0, nd we hve ξ Under the remining temporl guge trnsformtion we hve ˆt = t+ c 4ξ40, ˆx i = x i, 80 ˆΦ = Φ+ 2 ξ 40, ˆPi = P i ξ40,i, 8 nd ll the other vribles re guge-invrint. As n exmple, let us consider the two guge conditions used in Sec. VD nd Sec. VE. Let us ssume tht the x coordinte is the trnsverse-sher guge in Sec. VD, nd ˆx is the hrmonic guge in Sec. VE. Thus, we hve P i i 0 in x, nd 4 U +mȧ U + P i i 0 in ˆx. From Eq. 8 we hve thus 4 U +mȧ U = 0 2ξ40, 82 2ξ40 = 4 U +mȧ U. 83 Thus, when we trnsform the result known in the trnsverse-sher guge to the hrmonic guge, we use Eq. 8 with ξ 40 given in Eq. 83. All the other vribles re invrint under the guge trnsformtion. We cn show tht under this trnsformtion Eqs.67, 68 become Eqs.7, 72, respectively. In this wy, if we know the solution in ny guge condition the rest of the solutions in other guge conditions cn be simply derived without solving the equtions gin. Thus, s in other guge theories like Mxwell s or Yng-Mills theories the guge condition should be deployed to our dvntge in hndling the problem mthemticlly. VI. DISCUSSION Here we compre the PN pproch with the perturbtive one. The perturbtion nlysis is bsed on the perturbtion expnsion of the metric nd energy-momentum vribles in given bckground. All perturbtion vribles re ssumed to be smll. In liner order perturbtions we keep only the first-order devitions from the bckground 3, wheres 23