arxiv: v1 [astro-ph] 19 Nov 2007

Transcription

1 Long wvelength limit of evolution of cosmologicl perturbtions in the universe where sclr fields nd fluids coexist Tkshi Hmzki 1 Kmiyugi Hchioji-city Tokyo Jpn rxiv: v1 [stro-ph] 19 Nov 2007 Abstrct We present the LWL formul which represents the long wvelengh limit of the solutions of evolution equtions of cosmologicl perturbtions in terms of the exctly homogeneous solutions in the most generl cse where multiple sclr fields nd multiple perfect fluids coexist. We find the conserved quntity which hs origin in the dibtic decying mode, nd by regrding this quntity s the source term we determine the correction term which corrects the discrepncy between the exctly homogeneous perturbtions nd the k 0 limit of the evolutions of cosmologicl perturbtions. This LWL formul is useful for investigting the evolutions of cosmologicl perturbtions in the erly stge of our universe such s reheting fter infltion nd the curvton decy in the curvton scenrio. When we extrct the long wvelength limits of evolutions of cosmologicl perturbtions from the exctly homogeneos perturbtions by the LWL formul, it is more convenient to describe the corresponding exctly homogeneous system with not the cosmologicl time but the scle fctor s the evolution prmeter. By pplying the LWL formul to the reheting model nd the curvton model with multiple sclr fields nd multiple rdition fluids, we obtin the S formul representing the finl mplitude of the Brdeen prmeter in terms of the initil dibtic nd isocurvture perturbtions PACS number(s):98.80.cq 1 emil ddress: yj4t-hmzk@shi-net.or.jp

2 1 Introduction nd summry Recently we come to be required to investigte the evolution of cosmologicl perturbtions intheveryerlyuniverse[26],[6]. Accordingtotheinfltionryscenriosndthecurvton scenrio, in the erly universe the wvelength of cosmologicl perturbtions responsible for the present cosmic structures such s glxies nd clusters of glxies is much lrger thn the horizon scles. Therefore the methods for reserching the cosmologicl perturbtions on superhorizon scles hve been sought. In this context, Nmbu nd Truy pointed out tht there exists n LWL formul representing the k 0 limit of the cosmologicl perturbtions in terms of the exctly homogeneous perturbtions [32 ]. Soon lter in the multiple sclr fields system [15 ] [30 ], the complete LWL formule were constructed. Since the evolution equtions of the corresponding exctly homogeneous universe look simpler thn the evolution equtions of cosmologicl perturbtions, the LWL formul brought bout gret simplifiction. In ddition, the viewpoint tht the evolutions of the cosmologicl perturbtions on superhorizon scles re governed by the stbility nd instbility of the corresponding exctly homogeneous universe [9 ], [10 ] is useful for physicl interprettion. In this context, the phse spce of the corresponding exctly homogeneous system ws investigted in detil nd the role of the fixed points in the phse spce in the stbility nd instbility of cosmologicl perturbtions ws discussed [10 ]. For these resons, the LWL formul ws used for investigting the evolution of cosmologicl perturbtion on superhorizon scles by severl uthors [14 ], [30 ], [32 ], [9 ], [10 ]. In this pper, in section 2, in order to investigte the evolutionry behviors of cosmologicl perturbtions during reheting nd the curvton decys, we construct the complete LWL formule for the most generl system where the multiple sclr fields nd the multiple perfect fluids coexist. In the erly universe, the cosmologicl perturbtions on superhorizon scles responsible for the cosmic structures nd the CMB temperture nisotropies experience the reheting nd/or the curvton decys. In these processes, the multiple sclr fields such s infltons nd curvtons oscillte coherently, grdully decying into rdition fluid. By replcing the oscillting sclr fields with the dust fluids, the evolution of the cosmologicl perturbtions during reheting fter the infltion [8 ] nd in the curvton decy [20 ] were investigted. These uthors treted the system dominted by dust-like sclr field fluid nd rdition nd investigted the influence of the entropy perturbtion originting from the multicomponent property to the evolution of the totl curvture perturbtion vribles such s the Brdeen prmeter. The purpose of these nlyses ws to determine the initil perturbtion of the present Friedmnn universe in terms of the erly stge seed perturbtion. Although it ws shown prtilly tht this replcement is physiclly resonble [14 ], [8 ], we re required to tret the decying oscilltory sclr fields directly. In fct, the instbilities chrcteristic to the rpidly oscillting sclr fields were pointed out [3 ], [5 ], [37 ], [10 ], [36 ] nd investigted [33 ], [16 ], [31 ], [17 ]. In order to tret the oscilltory sclr fields directly, the ction ngle vribles were introduced nd the verging method representing the verging over the fst chnging ngle vribles ws pplied [9 ],[10 ]. In this pper, in order to investigte the evolutionry behviors of cosmologicl perturbtions during reheting nd/or the curvton decys, in section 3 the ction ngle vribles nd the ction ngle perturbtion vribles re introduced, nd in section 5, the verging method ws pplied to the decying oscilltory sclr fields. In the ppers [9 ], [10 ], the verging method ws used to investigte the corresponding 1

3 exctly homogeneous system, nd by the LWL formul the evolution of the cosmologicl perturbtion in the long wvelength limit ws constructed from the corresponding exctly homogeneous perturbtion. By using the LWL formul nd by pplying the verging method, in single oscilltory sclr field system [14] nd in nonresonnt multiple oscilltory sclr fields system [9 ] it ws shown tht the Brdeen prmeter is conserved, nd in resonnt multiple oscilltory sclr fields system [10 ] it ws shown tht the cosmologicl perturbtion including the Brdeen prmeter cn grow. In this pper, in section 6, 7, by using the LWL formul nd by pplying the verging method to the decying oscilltory sclr fields, we construct S formule representing the finl mplitude of the Brdeen prmeter in terms of the initil seed dibtic nd entropic perturbtions, in the reheting nd in the curvton decys, respectively. The orgniztion nd the summry of the pper is explined s follows. In section 2, we construct the LWL formul s for such sclr-fluid composite system bsed on the the philosophy of the pper [15 ]. The discrepncy exists between the evolution equtions of the cosmologicl perturbtions in the k 0 limit nd the evolution equtions of the exctly homogeneous perturbtions becuse the former contins k 2 Φ = O(1) terms nd the ltter does not, therefore this discrepncy should be corrected by the correction term which contributes the dibtic decying mode, but ny generl methods for determining such correction term hve not been presented yet, nd only in the multiple sclr fields system such correction term ws determined. We show in the k 0 limit the existence of the conserved quntity which hs origin in the dibtic decying mode nd which is relted with k 2 Φ. By regrding this conserved quntity s the source term, nd obtining the specil solution A, we correct the exctly homogeneous perturbtion A nd we obtin the complete LWL formul A = A +A in the most generl sclr-fluid composite system. In section 3, we point out tht it is more pproprite to use the scle fctor rther thn the cosmologicl time t s the evolution prmeter when we use the LWL formul. As for the sclr quntity T, we use the perturbtion vrible DT, D is the opertor which mps the exctly homogeneous sclr quntity T to the guge invrint perturbtion vrible representing thet fluctutioninthe fltslice. D defined inthis wy cnbeinterpreted s kind of derivtive opertor. In fct, the exctly homogeneous prt (DT) cn be expressed s the derivtive of T with respect to the solution constnt with the scle fctor fixed. In order to investigte the exctly homogeneous system contining oscilltory sclr fields, we use the ction ngle vribles I, θ. By using D defined in this wy, we cn define the ction ngle perturbtion vribles DI, Dθ whose exctly homogeneous prts re given s the derivtives of I, θ with solution constnt C with the scle fctor fixed. When we use the derivtive opertor D nd the LWL formule, it is essentil to use the scle fctor s the evolution prmeter. In section 4, we pply the LWL formule to the nonintercting multicomponents system nd discuss the long wvelength limit of the evolution of the Brdeen prmeter. In section 5, we pply the verging method by which the system is verged over the fst chnging ngle vribles to the decying sclr fields which hve been discussed in the reheting model nd the curvton model. By evluting the corrections produced by the verging process nd the errors produced by the trunction of the sufficient reduced ngle vribles dependent prt, the vlidity of the verging method is estblished. In section 6, 7, we pply the LWL formul nd the verging method to the intercting multicomponents model such s the reheting model, the curvton model, 2

4 respectively. We ssume tht the multiple sclr fields nd the multiple rdition fluid components exist. In these models, we construct the S formule representing the finl mplitude of the Brdeen prmeter in terms of the initil dibtic nd isocurvture perturbtions. In our previous pper [10 ], the evolutionry behviors of cosmologicl perturbtions in the erly universe where multiple oscilltory sclr fields interct with ech other hve been investigted. This S formul give the informtion bout how the cosmologicl perturbtions which grew in such erly er re trnsmitted into the rdition energy density perturbtions through the energy trnsfer from the sclr fields into the rdition fluids. We present the necessry condition for the initil entropic perturbtions produced in the erly er to survive until the lte rdition dominnt universe. Section 8 is devoted to discussions contining non-liner generliztion of our LWL formilism nd comment of the cse where the decy rte depends on other physicl quntities. In ppendices, the proofs of the propositions presented in section 5 nd the evlutions of the useful mthemticl formule used in section 6 re contined. In this pper, we consider the cse where the homogeneous sclr fields obey the phenomenologicl evolution equtions s φ+3h φ+ U φ +S = 0. (1.1) The interctions between sclr fields re described by the interction potentil U, while the interction between sclr fields nd fluids re described by S. This nlysis includes the well known cse S = Γ φ [11 ], [21 ] nd the generl cse where S is n rbitrry nlytic function of φ, φ which ws discussed in the pper [35 ] but whose perturbtions hve not been investigted yet. Another supplementl purpose of our pper is to present the evolution equtions of cosmologicl perturbtions corresponding to the homogeneous system with vrious source term S especilly dependent on φ. The nottion used in this pper re bsed on the the review [12 ] nd the pper [15 ]. 2 Derivtion of the LWL formul We give the definitions nd the evolution equtions s for the bckground nd the perturbtion vribles. Bsed on these nottions, in the most generl model where the multiple sclr fields nd the multiple perfect fluid components interct, we give the LWL formul representing the evolutions of the perturbtions vribles in terms of the exctly homogeneous solutions. We consider perturbtions on sptilly flt Robertson-Wlker universe given by ds 2 = (1+2AY)dt 2 2BY j dtdx j + 2 [(1+2H L Y)δ jk +2H T Y jk ]dx j dx k, (2.1) where Y, Y j nd Y jk re hrmonic sclr, vector nd tensor for sclr perturbtion with wve vector k on flt three-spce: Y := e ik x, Y j := i k ( j 1 k Y, Y jk := 3 δ jk k ) jk k Y. (2.2) k 2 3

5 By using the guge dependent vribles R nd σ g representing the sptil curvture perturbtion nd the sher, respectively: R := H L H T, σ g := kḣt B, (2.3) we cn define two independent guge invrint vribles: ( ) R A := A, Φ := R H H k σ g. (2.4) In order to define the mtter perturbtion vribles, we will consider the sclr quntity perturbtion vribles generlly. As for covrint sclr quntity T = T +δty, we define the guge invrint perturbtion vrible representing the T fluctution in the flt slice: T DT := δt R, (2.5) H Next we consider the covrint sclr quntity T 2 whose bckground quntity is the time derivtive of T: T. The extension of T into the covrint sclr quntity T 2 is not unique. For exmple, ( )[ 1/2, T 2 = sgn 0 T g µν µ T ν T] (2.6) nd where ñ µ is n rbitrry vector field stisfying T 2 = ñ µ µ T, (2.7) ñ µ ñ µ = 1, (2.8) hve the sme T s the bckground prt. But these defferent T 2 s give the unique perturbtion prt: DT 2 = (DT) TA. Therefore we cn define DT by D T := (DT) TA. (2.9) We consider the the universe where the sclr fields φ, (1 N S ) nd the fluids ρ α,p α (1 α N f ) coexist, whose energy momentum tensor is divided into A = (S,f) prts where S represents the multiple sclr fields, f represents the multiple fluids. The energy momentum tensor of f prt re further divided into individul fluids prts α. On the other hnd, the energy momentum tensor of S prt cnnot be divided into individul sclr fields prts, since the interction potentil U contins the terms consisting of plurl sclr fields φ : T µ ν = 0 = ) ( Tµ ν + S ) ( Qµ S + ) ( Tµ ν f ) ( Qµ f = = ) ( Tµ ν + S α ) ( Qµ S + α where the energy-momentum trnsfer vector Q Aµ is defined by T µ αν, (2.10) Q αµ, (2.11) ν Tν Aµ = Q Aµ = Q A ũ µ + f Aµ, (2.12) 4

6 where ũ µ is the four velocity of the whole mtter system nd the momentum trnsfer f Aµ stisfies ũ µ faµ = 0. For the sclr perturbtion, the energy-momentum tensor nd the energy-momentum trnsfer vector of ech individul component re expressed s nd T 0 A0 = (ρ A +δρ A Y), (2.13) T 0 Aj = (ρ A +P A )(v A B)Y j, (2.14) T j Ak = (P A δ j k +δp AYδ j k +Π TAY j k ), (2.15) Q A0 = [Q A +(Q A A+δQ A )Y], (2.16) Q Aj = [Q A (v B)+F ca ]Y j, (2.17) where ρ A, P A nd Q A re the bckground quntities of the energy density, the pressure nd the energy trnsfer of the individul component A, respectively. The nisotropic pressure perturbtion Π TA nd the momentum trnsfer perturbtion F ca re lredy guge invrint. As for the sclr quntities T = (ρ A,P A,Q A ), we use DT s the guge invrint perturbtion vribles. As for the guge invrint velocity perturbtion vrible, we use Z A := R H k (v A B). (2.18) The energy-momentum tensor of sclr fields prt is given by ) ( Tµ ν S = µ φ ν φ 1 2 ( λ φ λ φ+2 Ũ) δ µ ν. (2.19) Since divergence of the energy momentum tensor is given by ( ) ) ( µ Tµ ν = φ Ũ ν φ, (2.20) S φ in order tht the phenomenologicl equtions of motion of the sclr fields become we ssume tht ( Qν ) φ Ũ φ = S, (2.21) S = S ν φ, (2.22) By using the sclr fields bckground vribles φ, φ, S nd the corresponding perturbtion vribles Dφ, D φ, DS, the bckground prt of the fluid vribles re given by ρ S = 1 2 ( φ) 2 +U, (2.23) P S = 1 2 ( φ) 2 U, (2.24) h S = ) 2, ( φ (2.25) Q S = S φ, (2.26) 5

7 nd the perturbtion prt of fluid vribles re given by (Dρ) S = ρ S φ Dφ+ ρ S φ D φ, (2.27) (DP) S = P S φ Dφ+ P S φ D φ, (2.28) (hz) S = H φ Dφ, (2.29) (Π T ) S = 0, (2.30) (DQ) S = S D φ φ DS, (2.31) (F c ) S = S ( kdφ k ) H φ Z. (2.32) When the source of the sclr field φ, S is given s functions of the covrint sclr quntities T nd T 2 whose bckground prt is T, tht is S = S ( T, T 2 ), DS is given by In such cse, (DQ) S cn be written s DS = S T DT + S T D T. (2.33) (DQ) S = Q S T DT + Q S T D T, (2.34) which is ssumed from now on. In the sme wy s the individul components T µ Aν, s for the totl energy-momentum tensor T µ ν = T µ A Aν, we cn define the guge invrint perturbtion vribles such s Dρ, DP, hz nd Π T. From (2.10),(2.11), we obtin the bckground equtions s ρ = ρ S + α P = P S + α h = h S + α 0 = Q S + α ρ α, (2.35) P α, (2.36) h α, (2.37) Q α, (2.38) 6

8 nd perturbtion equtions s Dρ = Dρ S + α DP = DP S + α Dρ α, (2.39) DP α, (2.40) hz = (hz) S + α Π T = (Π T ) S + α 0 = (DQ) S + α 0 = (F c ) S + α h α Z α, (2.41) Π Tα, (2.42) DQ α, (2.43) F cα. (2.44) This Z is known s the Brdeen prmeter [2 ] [12 ], [23 ]. In the long wvelength limit, the Brdeen prmeter is conserved in the cse where the entropy perturbtions re negligible. But in vrious systems it ws reported tht the entropy perturbtions cnnot be neglected [3], [5], [6],[10], sointhepresentpperwewillinvestigtetheevolutionrybehviorofthe Brdeen prmeter more crefully. Until now, s for the guge invrint sclr quntity perturbtion vribles, we use D. But trditionly most sclr quntity perturbtion vribles hve been written without using D: Y := Dφ, ρ α gα := Dρ α, P α Π Lα := DP α, Q α E gα := DQ α. (2.45) This Y hs been clled the Sski-Mukhnov vrible [28 ] [22 ]. In terms of the guge independent vribles defined bove, we give the evolution equtions of cosmologicl perturbtions. From (2.21), the bckground nd the perturbtion prts cn be written s where φ+3h φ+ U φ +S = 0, (2.46) L 1 (DT,A) = k2 k2 φ 2Dφ Φ, (2.47) 2 H L 1 (DT,A) = (Dφ) +3H(Dφ) + 2 U Dφ+DS φ A+2( U φ φ φ As for the fluid components, µ Tµ αν = Q αν gives the bckground equtions s nd the perturbtion equtions s +S)A. (2.48) ρ α = 3Hh α +Q α, (2.49) L 2α (DT,A) = k2 2 H h α(φ Z α ), (2.50) ( ) hα Z α +3h α Z α +h α A+DP α 2 H 3 Π Tα = k F cα + Q α Z, (2.51) H 7

9 where L 2α (DT,A) = (Dρ α ) +3HDρ α +3HDP α Q α A DQ α. (2.52) G µ ν = κ 2 Tµ ν gives the bckground equtions s H 2 = κ2 ρ, (2.53) 3 ρ = 3Hh, (2.54) Ḣ = 3 2 (1+w)H2, (2.55) nd the perturbtion equtions s where L 3 (DT,A) = 2ρ k2 3 2 H2Φ, (2.56) L 4 (DT,A) = κ2 3 Π T k2 2Φ, (2.57) A+ 3 (1+w)Z = 0, (2.58) 2 A+ 1 ( ) H Φ κ 2 = Π k 22 T, (2.59) L 3 (DT,A) = 2ρA+Dρ, (2.60) L 4 (DT,A) = HA+2ḢA κ2 (Dρ+DP). (2.61) 2 The dynmicl perturbtion vribles re clssified into two groups, tht is, wht hs nlogy with the exctly homogeneous perturbtions nd wht is not relted with the exctlty homogeneous perturbtions t ll. The dynmicl perturbtion vribles of the former type re DT representing the sclr quntity T = (ρ,p,φ,q,s) perturbtion in the flt slice, DT nd the metric perturbtion vrible A. The dynmicl perturbtion vribles of the ltter type re the Newtonin grvittionl potentil Φ nd Z A, F ca, Π TA which hve vector or tensor origin. In the bove L i (i = 1,,4) equtions, the former type dynmicl perturbtion vribles re contined in the left hnd side while the ltter type perturbtion vribles re collected in the right hnd side. The exctly homogeneous perturbtions DT nd A corresponding to DT nd A, respectively re constructed s ( ) T (DT) T := C t H R, (2.62) ( R A ) :=, (2.63) H R := 1 ( ), (2.64) C t where C is the solution constnt of the bckground solution nd the subscript t implies tht the derivtive with respect to C is performed with the cosmologicl time t fixed. On 8

10 the other hnd, the dynmicl perturbtion vribles of the ltter type such s Φ, Z A, F ca nd Π TA do not hve exctly homogeneous counterprts. The evolution equtions of cosmologicl perturbtions contining L i (i = 1,,4) hve nlogy in the exctly homogeneous perturbtion equtions. In fct, the vritions of the exctly homogeneous equtions (2.46), (2.49), (2.53) nd (2.55) give L i (DT,A ) = 0 (i = 1,,4), (2.65) respectively. The only difference between the exctly homogeneous perturbtion L i (i = 1,,4)equtions ndthectul k 0 cosmologicl perturbtion L i (i = 1,,4)equtions is tht k 2 Φ terms exist in the ltter but k 2 Φ terms do not exist in the former. Then the effect ofthesourcetermk 2 Φiscorrected inthefollowing wy. Inperformimgthecorrection process, it is importnt to notice tht the source terms k 2 Φ cn be represented in terms of conserved quntity which hs origin in the universl dibtic decying mode. In fct, s for f defined by ( f = 3 H A+ 1 ) 2 g = k2 Φ, (2.66) 3H using (2.50), (2.56), (2.57) yields When we ssume tht for k 0 limit df dt = 3 H 2 wπ T k2 (1+w)Z. (2.67) Π T 0, kz 0, (2.68) re stisfied, the quntity f is conserved, whose vlue is written s c. Therefore for k 0 limit, k 2 Φ 3H c = O(1). (2.69) This expression of Φ is well known s tht of the universl dibtic decying mode [15 ]. In the L i (i = 1,,4) equtions contining DT, A, the Newtonin potentil Φ ppers only in the form k 2 Φ, tht is, ccompnied by k 2. When we ssume tht DT = O(1), A = O(1), k 2 Φ behves s O(1). Since in the liner perturbtion, the scle of the perturbtion vribles is rbitrry, the fct tht Φ = O(1/k 2 ) does not imply the brekdown of the liner perturbtion. If one wnt to get Φ = O(1), one simply ssumes tht DT = O(k 2 ), A = O(k 2 ). But s explined lter, we cnnot ssume tht Φ is vnishing, since c defined by (2.69) must stisfy the constrint (2.80). Therefore in the k 0 limit where (2.68), (2.69) re stisfied, (2.47), (2.50), (2.56), (2.57) cn be written s L 1 (DT,A) = 3 φ c, (2.70) 3 L 2α (DT,A) = 3h α c, 3 (2.71) L 3 (DT,A) = 2ρ c, 3 H (2.72) L 4 (DT,A) = 3H c. (2.73) 3 9

11 It cn be verified tht bove four sets of equtions (2.70), (2.71), (2.72) nd (2.73) re stisfied by where A = 3 2 (1+w)g + ġ H, (2.74) T DT = g, H (2.75) g = c dt 1 t 0 3. (2.76) This specil solution for A = (DT,A) is written s A. Since the vrition of the exctly homogeneous solution A stisfies (2.65), the generl solutions of (2.70), (2.71), (2.72), (2.73) A = (DT,A) cn be expressed s A = A +A. (2.77) The perturbtion equtions except L i (i = 1,,4) equtions hve vector origin, tht is, they re derived from the spce component of the Einstein equtions. Therefore these perturbtion equtions do not hve ny nlogy with the exctly homogeneous perturbtion equtions. As explined in the pper [15 ], these perturbtion equtions determine the evolutions of the dynmicl perturbtion vribles which hve vector or tensor origin, tht is, which hve no correspondence with the exctly homogeneous pertubtions, or give the constrint which should be stisfied in order tht the exctly homogeneous perturbtions become the k 0 limit of evolutions of cosmologicl perturbtions. Therefore (2.51), (2.58) cn be interpreted s the decision of the evolution of the vribles Z α which is not relted to the exctly homogeneous solution t ll in terms of A = (DT,A), the constrint to the exctly homogeneous perturbtions, respectively. Integrting (2.51) yields h α Z α H [C 3 α + dt 3 t 0 ( h α A DP α k F cα + Q α H Z, )]. (2.78) By summing (2.78) with respect to ll the fluid components, we obtin [ (hz) f = H ( ) ] C 3 α + dt 3 h f A (DP) f S Dφ α t 0 [ = H ( C 3 α + 2 ) ( 3 3H ρa+3 φ Dφ 2 ) ] 3 3H ρa+3 φ Dφ (2.79) α 0 Therefore(2.58)givestheconstrintbetween C α, cdefinedby(2.69), nd2n S +N f solution constnts of the exctly homogeneous perturbtion s C α + 2 ( κ 2c 2 ) 3 3H ρa+3 φ Dφ = 0. (2.80) Integrting (2.59) gives α Φ = H 0 ( ) C t dta, (2.81) t 0 10

12 where the first term contining C t is well known universl dibtic decying mode [15 ] nd by compring with (2.69) we obtin C t = 3 k2c. (2.82) If we ssume c = 0, since (2.80) gives one constrint reltion, we obtin 2N S + 2N f 1 solutions nd the Newtonin potentil is obtined by (2.81) with C t = 0: Φ = H dta, (2.83) t 0 If we ssume tht c is nonvnishing, since c = O(1), A O(1), therefore A C t O(k 2 ), (2.84) we obtin the k 0 limit of the universl dibtic decying mode [15 ]: 1 3 k2 Φ H c, (2.85) which is consistent with (2.69). Then we hve obtined the long wvelength limit of ll the solutions to the evolution equtions of cosmologicl perturbtions. 3 Use of the scle fctor s the evolution prmeter As the guge invrint vrible representing the fluctution of the sclr quntity T, we dopt DT defined by (2.5) which represents the T fluctution in the flt slice, since it is the esiest to see the correspondence with the exctly homogeneous perturbtion of T. While until nowwe described theexctly homogeneous vriblessfunctions oft, C where tisthe cosmologicl time nd C s re solution constnts, we cn describe the exctly homogeneous vribles s functions of, C where is the scle fctor. For n rbitrry sclr quntity such s ρ, P, S, Q, φ, from (2.62), (2.64), (DT) cn be written s the prtil derivtive of the corresponding exctly homogeneous sclr quntity T with respect to solution constnt C with the scle fctor fixed: ( ) T (DT) =. (3.1) C Since ( ) 1 = 1 ( ) ρ = 1 ȧ Cȧ 2ρ C 2ρ (Dρ) = A, (3.2) this property of D lso holds s for the time derivtive of the sclr quntity T: ( DT ) ( = T C ). (3.3) Therefore the opertor D defined by (2.5) cn be interpreted s kind of derivtive opertor, tht is, the derivtive with respect to the solution constnt C with fixed. Becuse 11

13 of this derivtive property of D, s for the sclr quntities T = (ρ,p,s,q) which re functions of φ, φ, we cn understnd DT = T φ Dφ+ T φ D φ, (3.4) esily. The contribution to DT from the dibtic decying mode is given by the sme form s tht of DT: ( DT ) T = g, (3.5) H where g is defined by (2.76). As seen from the bove duscussion, in order to derive the long wvelength limit of cosmologicl perturbtions from the corresponding exctly homogeneous system by using the LWL formul, it is more pproprite to use s the evolution prmeter the scle fctor thn the cosmologicl time t. For exmple, s for the sclr-fluid composite system, the corresponding exctly homogeneous expressions re obtined by solving the first order differentil equtions setting φ, p := 3 φ, ρ α s independent vribles nd the scle fctor s the evolution prmeter: d d φ = 1 p H 3, d d p = 3 H U φ 3 H S, replcing H with the right hnd side of the Hubble lw: [ H 2 = κ2 1 p U(φ)+ α d d ρ α = 3h α + Q α H, (3.6) ρ α ]. (3.7) While under use of t s the evolution prmeter our system is the constrined system with the Hmiltonin constrint (3.7), under use of s the evolution prmeter our system becomes the unconstrined system with respect to independent vribles φ, p := 3 φ, ρ α. The corresponding first order perturbtion vribles re Dφ, P := 3 D φ, Dρ α. For some time, we consider the system consisting of multiple sclr fields φ only. Since the evolutions of φ, p := 3 φ cn be described in terms of the Hmilton equtions of motion, the evolutions of the corresponding perturbtion vribles Y = Dφ P := 3 D φ cn lso be written in terms of the Hmilton equtions of motion: whose Hmiltonin is given by dy dt = H P, dp dt H = 1 2 3P P V b Y Y b + κ2 V b = 2 U φ φ b + 3κ2 2 φ φb + κ2 2H 12 = H Y, (3.8) 2H φ φb P Y b, (3.9) ( U φb + φ φ U )+ k2 φ b 2δ b. (3.10)

14 When we discuss the quntiztion of the fluctutions, the other set of the cnonicl perturbtion vribles Ỹ := Y, P := 3 Ẏ hs been used [23 ]. But in the viewpoint of the LWL formlism, the set of the cnonicl vribles Y := Dφ, P := 3 D φ is more nturl thn the set of the cnonicl vribles Ỹ := Y, P := 3 Ẏ, becuse the long wvelength limit of the former set is generted from the derivtive of the homogeneous vribles φ, p := 3 φ with respect to the solution constnts with the scle fctor fixed. The connection between the old cnonicl vribles Ỹ := Y, P := 3 Ẏ nd the new cnonicl vribles Y := Dφ, P := 3 D φ re given by the cnonicl trnsformtion defined by the generting function: W = ỸP H 4 ρ φ φb Ỹ Ỹ b. (3.11) When we tret the oscilltory sclr fields, the ction ngle vribles I, θ re useful [9 ], [10 ]: φ = 1 2I cosθ 3/2, m p = 3/2 2m I sinθ, (3.12) where m is the mss of the sclr field φ. The ction ngle vribles obey the evolution eqution s d d I = 3 U int + 3/2 2I sinθ S +3I cos2θ, (3.13) H θ H m d d θ = m H + 3 H U int + 3/2 I H 1 2m I cosθ S 3 2 sin2θ. (3.14) In order to investigte the cosmologicl perturbtions in the universe contining oscilltory sclr fields, by using D defined in the bove we define the ction ngle perturbtion vribles DI, Dθ strting from Y := Dφ, P := 3 D φ. In the LWL formlism, the perturbtion vribles corresponding with the ction ngle vribles I, θ re DI, Dθ defined by the following expressions: [ ] 1 2I Y = D cosθ 3/2, m P = D [ 3/2 ] 2m I sinθ. (3.15) where D in the right hnd side is interpreted s D = DI I + Dθ θ. (3.16) The expressions obtined from vritions of I, θ with fixed in the previous ppers [9 ], [10 ] re the long wvelength limits of DI, Dθ defined by (3.15). In fct, the LWL formule ( ) ( ) I I DI = + C H 3I g, (3.17) Dθ = ( θ C ) 13 + θ g, (3.18) H

15 where g is defined in (2.76) hold. While the third term in the right hnd side of (3.17) ppers becuse of the scle fctor dependence of the trnsformtion lw from φ, φ to I, θ, the prts of (3.17) nd (3.18) reflect the fct tht D is the derivtive opertor with respect to the solution constnt with the scle fctor fixed. In order to solve the dynmics of the system contining the oscilltory sclr fields, we re required to perform the verging over the fst chnging ngle vribles θ [9 ], [10 ]. If we use the cosmologicl time t s the evolution prmeter, our system is constrined system. Therefore we must check tht our verging procedure is consistent with the constrint nd this process is rther cumbersome. But if we use the scle fctor s the evolution prmeter, our system becomes unconstrined system, so the definition of the verging procedure becomes rther simple. We cn conclude tht use of the scle fctor s the evolution prmeter brings bout the two merits. One is tht it becomes esier to see the correspondence between the exctly homogeneous solution nd the long wvelength limit of the first order perturbtion nd tht the LWL formule become more simple. The other is the more simple definition of the verging process. We consider the evolution of the Brdeen prmeter Z. Following the pper [34 ], we define ζ s the guge invrint vrible representing the curvture perturbtion in the uniform density slice: ζ := Ḣ ρ Dρ = R Ḣ δρ. (3.19) ρ From (2.56) (2.58) we cn see tht the Bredeen prmeter Z nd ζ re closely connected s Z = ζ 2 1 k 2 91+w 2 H2Φ, (3.20) whose k 0 limit is Z = ζ 2 1 c 31+w 3 H, (3.21) where c is constnt relted with the dibtic decying mode defined by (2.69). While the Brdeen prmeter Z is expressed s the weighted sum of Z S, Z α which do not relted with the exctly homogeneous quntity t ll nd whose evolution is written in the rther cumbersome integrl form (2.78), ζ is connected with the corresponding exctly homogeneous quntity ndζ evolution cnbe written in terms of thederivtive of the totl energy density ρ with respect to solution constnt. Then we consider much esier ζ evolution. In the pper [19 ], s the nonliner generliztion of the Brdeen prmeter ζ, ζ(t,x) ws introduced. When P = P(ρ), ζ(t,x) is reduced to ζ(t,x) = ln(t,x)+ 1 3 ρ(t,x) In fct, the first order quntity of ζ(t,x) is given by ζ 1 (t,x) = δ(t,x) dρ ρ+p(ρ). (3.22) + 1 δρ(t, x) 3ρ+P(ρ), (3.23) which grees with the Brdeen prmeter ζ. In the viewpoint of the LWL formlism, we doptthezerocurvtureslice (t,x)/ x i = 0. WessumetheequtionofstteP r = ρ r /3, 14

16 since the finl stte of reheting nd of the curvton decy is rdition dominnt. In this cse, the bove ζ(t,x) is reduced to ζ(,x) = 1 4 lnρ r(,c(x)). (3.24) ρ r (,C(x)) is the expression of ρ r obtined by solving the loclly homogeneous system (see the seprte universe pproch [34 ]) with use of s the evolution prmeter. The solution constnts C(x) hve sptil dependence. By considering C(x) = C + δc(x) nd expnding with respect to δc(x), we cn obtin the perturbtion of n rbitrry order up to the decying modes of perturbtions. For exmple, the n-th order perturbtion is given by 1 n! i n C 1 C 2 C n [ ] 1 4 lnρ r(,c) δc 1 (x)δc 2 (x) δc n (x). (3.25) When we obtin the exctly homogeneous exression ρ r = ρ r (,C), we cn know the long wvelength limits of perturbtions of rbitrry orders. Lter we determine the exressions ρ r = ρ r (,C) in the reheting nd in the curvton decy. In the reheting nd in the curvton decy, δc(x) s re given by the ction ngle vribles in the initil time: δi (x) := δi ( = 0,x), δθ (x) := δθ ( = 0,x). (3.26) We discuss how to determine the sttisticl properties of δi (x), δθ (x). δi (x), δθ (x) re given t the time when the slow rolling phse ends nd the coherent oscilltion begins. At this time, φ = 0, nd s for the perturbtions δφ (x) = d 3 ke ik x e (k), δ φ (x) = 0, (3.27) where e (k) is the Gussin rndom vrible stisfying < e (k)e b (k ) >= P (k)δ b δ 3 (k+k ), k := k. (3.28) By solving I θ in terms of φ φ nd by Tylor expnding, we obtin δi (x) = 3 0 m φ δφ (x) m [δφ (x)] 2, (3.29) δθ (x) = 0, (3.30) where we use φ = 0, δ φ (x) = 0. The bove fct tht δθ ( = 0,x) = 0 does not imply δθ (,x) = 0, since θ () is function depending on not only θ ( 0 ) but lso I ( 0 ). Therefore it is often importnt to consider the role of the perturbtions of ngle vribles. Determining the mny point correltion function of fluctutions is reduced to the evlution of < e 1 (k 1 )e 2 (k 2 ) e n (k n ) >. (3.31) This quntity cn be determined by pplying the differentil opertion defined by [ ] exp d 3 δ δ kp (k). (3.32) δe (k) δe ( k) e=0 15

17 4 Appliction of the LWL formul to the non-intercting multicomponent system Bsed on the results obtined in sections 2, 3, we consider the long wvelength limit of the evolutions of cosmologicl perturbtions in the universe consisting of multiple cosmic components. For simplicity, we consider the cse where ech component do not interct with ech other. We consider the w α fluid where w α is constnt nd which does not interct with other components Q α = 0. ρ α is solved s ρ α = A α 3(1+wα), (4.1) where A α is solution constnt nd by differentiting with respect to A α with the scle fctor fixed, we obtin (Dρ α ) = δa α 3(1+wα), (4.2) where δa α is perturbtion solution constnt corresponding with A α. Therefore we obtin gα = δa α A α = const. (4.3) Next we consider the cse tht the oscilltory sclr field φ does not interct with other cosmic components, tht is U int = 0, S = 0. Since in such cse the right hnd side of (3.13) is oscilltory function depending on the ngle vrible θ with vnishing men vlue, by tking the verging over θ, we obtin I = A, (4.4) where A is constnt. The estimte bout the effects of the oscilltions due to the fst chnging ngle vribles θ ws discussed in the previous ppers [9 ] [10 ]. Therefore by tking the derivtive with respect to solution constnt A with the scle fctor fixed, we obtin DI = δa, (4.5) where δa is perturbtion constnt corresponding with A. Since ρ = m I 3, Dρ = m DI 3, (4.6) where the bove second expression is given by the D opertion to the bove first expression, we obtin g = DI I = δa A = const. (4.7) As seen from the bove expressions, s for the oscilltory sclr field φ the energy density ρ nd the energy density perturbtion ρ g behve like those of the dust fluid. We cn summrize tht for non-intercting cosmic components, gα, g re conserved. 16

18 As for the multicomponent non-intercting fluids system, ζ evolution is given by ( ) ( ζ δa α / = 3 ) A α (1+w 3(1+wα) α ). (4.8) 3(1+wα) α α We cn see tht ζ is exctly conserved for the dibtic growing mode defined by δa α A α 1 1+w α = αindependent. (4.9) For some time, we consider the two components system consisting of dust nd rdition. In this cse, (4.8) is reduced to ζ = δa d +δa r 3A d +4A r, (4.10) where the suffix d, r imply dust, rdition, respectively. We obtin ζ init = 1 δa r, ζ fin 4 A = 1 δa d, (4.11) r 3 A d in the limit 0,, respectively. We dopt different, more physicl prmeriztion: δa d = 3ξA d ηa d, δa r = 4ξA r. (4.12) ξ represents the dibtic growing mode nd η represents the isocurvture mode defined by ζ init = 0, S rd = 3 4 gr gd = const =: η. (4.13) Then we obtin ζ fin = ζ init 1 η, (4.14) 3 which is the fmous formul.[13 ] Although in the pper [34 ], (4.10) hs lredy been derived essentilly without using the LWL formul, our result is more rigorous in the point tht we tret the contribution from the dibtic decying mode chrcterized by c defined by (2.69) more ppropritely, while the pper [34 ] simply ssumes tht k 2 Φ is vnishing. 5 Appliction of the verging method to the decying sclr fields We derive the evolution equtions of the multiple sclr fields decying into the multiple rdition fluids. By solving these evolution equtions nd tking the exctly homogeneous perturbtions, we cn obtin the informtion of the evolutionry behviors of cosmologicl perturbtions during reheting nd in the curvton model. We ssume tht the source S is given by S = Γ φ. (5.1) 17

19 We nondimensionlize the dynmicl quntities s 0, I I 0 I, ρ α ρ 0 ρ α, U int ρ 0 U int = O(ν), (5.2) nd the prmeters s m m 0 m, Γ Γ 0 Γ, (5.3) where ν is the smll prmeter implying the rtio of the interction energy to the free prt energy defined by Then we obtin the dimensionless prmeters ρ 0 := m 0I 0. (5.4) 3 0 ǫ := κ ρ 1/2 0, γ := 3 m 0 3Γ0, (5.5) κρ 1/2 0 which imply the rtio of the Hubble prmeter to the mss of the sclr fields H/m, the rtioofthedecy rtetothehubbleprmeter Γ/H ttheinitiltime = 0, respectively. By using the bove dynmicl vribles nd prmeters, the evolution equtions cn be expressed s d d I = 1 3 U int I γγ ǫρ 1/2 θ ρ (1 cos2θ )+3I 1/2 cos2θ, (5.6) d d θ = 1 ǫ d d σ α = γ m ρ + 1 1/2 ǫ ρ 1/2 3 ρ 1/2 U int I 1 2 γ Γ ρ 1/2 sin2θ 3 2 sin2θ, (5.7) Γ α m I (1 cos2θ ), (5.8) where σ α is defined by ρ α = σ α 4, (5.9) nd Γ α is the decy rte from the sclr field I to the rdition component σ α nd therefore Γ is given by Γ = Γ α. (5.10) α In this pper, we investigte the evolutionry behvior of cosmologicl perturbtions during the period when the decy rtes from the sclr fields to the rdition fluids re lrge compred to the interction between the sclr fields, tht is γ ν/ǫ, while in the pper [10 ] the evolutions of cosmologicl perturbtions during the period when the interction between the sclr fields is dominnt ν/ǫ γ were discussed, which is thought to give the initil conditions for the present studies. Next we show tht there exists trnsformtion such tht in system obtined by tht trnsformtion, the dynmics of the ction vribles I nd the rdition energy densities σ α cn be determined independently of the ngle vribles θ. In order to show this sttement, we put severl ssumptions. 18

20 (i) The interction energy of the sclr fields U int is nlytic with respect to the dynmicl vribles I, θ nd is 2π periodic with respect to θ. U int is bounded s U int ν 9/2 I, (5.11) which implies tht there exists positive constnt M such tht 9/2 U int νm I, (5.12) where I := I. (5.13) (ii) As for the totl energy density ρ = m I 3 + α σ α 4 +U int, (5.14) the ction vribles I nd the rdition energy densities σ α stisfy c 1 m I + α σ α c 2, (5.15) for some positive constnts c 1, c 2. Note tht the evolution equtions of the system cn be written s d d I = F (I,σ,θ,), (5.16) d d σ α = F α (I,σ,θ,), (5.17) d d θ = 1 ǫ ω (I,σ,)+G (I,σ,θ,), (5.18) where F F α G re nlytic with respect to the dynmicl vribles I σ α θ nd 2π periodic with respect to the ngle vribles θ α. We sy tht the evolution equtions is of the type C k : (i) The verged prts of F F α G re bounded s < F > 2 I, (5.19) < F α > 3 I, (5.20) < G > 2, (5.21) where < A > implies the verging over the ngle vribles θ : 1 < A >:= (2π) N S 2π 0 d N S θa. (5.22) In cse of the resonnt cse, it is prescribed tht the verging is performed with respect to the fst ngle vribles only [10 ]. 19

21 (ii) The oscilltory prts of F F α G re bounded s F ǫ k 2 I, (5.23) F α ǫ k 3 I, (5.24) G ǫ k 2, (5.25) where à implies the residul prt fter the verging over the ngle vribles θ : Under this nottion, the following proposition holds. à := A < A >. (5.26) Proposition 1 Let k be some non-negtive integer nd consider the evolution equtions s d d I(k) d d σ(k) α d d θ(k) = F (k) (I (k),σ (k),θ (k),), (5.27) = F (k) α (I (k),σ (k),θ (k),), (5.28) = 1 ǫ ω(k) (I(k),σ (k),)+g (k) (I(k),σ (k),θ (k),). (5.29) Suppose tht this set of evolution equtions is of the type C k, there exists trnsformtion I (k) σ (k) α θ (k) = I (k+1) = σ (k+1) α = θ (k+1) stisfying the following conditions: +u (k) +u (k) (I (k+1),σ (k+1),θ (k+1),), (5.30) α (I(k+1),σ (k+1),θ (k+1),), (5.31) +v (k),σ (k+1),θ (k+1),), (5.32) (i) u (k) u (k) α v (k) re nlytic with respect to the dynmicl vribles I (k+1) re 2π periodic with respect to the ngle vribles θ (k+1) σ (k+1) α, nd re bounded s θ (k+1), u (k) ǫ k+1 I (k+1), (5.33) u (k) α ǫ k+1 I (k+1), (5.34) v (k) ǫ k+1. (5.35) (ii) The evolution equtions of the trnsformed vribles I (k+1) σ α (k+1) type C k+1 nd the chnges of < F > < F α > < G > re bounded s θ (k+1) re of the where < A > is defined by < F > ǫ k+1 2 I (k+1), (5.36) < F α > ǫ k+1 3 I (k+1), (5.37) < G > ǫ k+1 2, (5.38) < A >:=< A (k+1) (I (k+1),σ (k+1),θ (k+1),) > < A (k) (I (k+1),σ (k+1),θ (k+1),) >. (5.39) 20

22 (For the proof, see Appendix A.) This proposition implies tht we cn mke the prt depending on the ngle vribles rbitrrily smll by tking the originl set of the evolution equtions s the strting point k = 0 nd pplying trnsformtions given in the proposition itertively. Therefore it cn be expected tht the evolution of the dynmicl vribles cn be described by the truncted system obtined by discrding the ngle vrible dependent prt with sufficiently good ccurcy, if we tke sufficiently lrge k. By estimting the errors produced by the trunction, we show tht the bove expecttion is correct. We use the symbol to represent the difference of quntity for the exct system nd corresponding quntity for the truncted system. For A = (I,σ α,θ ), the errors of the bckground vribles re written s A := A A tr, (5.40) nd the errors of the perturbtion vribles re written s δa := δa δa tr, (5.41) where A, δa represent quntities of the exct system nd A tr, δa tr represent quntities of the corresponding truncted system. For function f(), let us write when f() is bounded s f() = E( 2 ), (5.42) f() p()exp( λ 2 ) (5.43) for polynomil of : p(), nd for positive number λ. For function f(), let us define f() by f() := sup f( ). (5.44) 1 When we write ll the inequlities, it is prescribed tht ll the coefficients of order unity re omitted. The trunction error for the m-th order system cn be estimted s follows. Proposition 2A Let m be n integer lrger thn or equl to 2. For the m-th order system, the trunction errors of the bckground vribles re given by I E( 2 )ǫ m, (5.45) σ () ǫ m, (5.46) θ 2 ǫ m 1, (5.47) nd the trunction errors of the perturbtion vribles re given by δi E( 2 )ǫ m 1 δa 1 (1), (5.48) δσ ǫ m 1 δa 1 (1), (5.49) δθ 2 ǫ m 2 δa 1 (1)+ǫ m 1 exp( 2 ǫ m ) [ 2 ǫ δθ(1) + 4 δa m (1) ], (5.50) where δa m (1) := δi(1) + δσ(1) +ǫ m δθ(1), (5.51) 21

23 under the initil conditions I(1) = σ(1) = θ(1) = 0, (5.52) δi(1) = δσ(1) = δθ(1) = 0. (5.53) (For the proof, see Appendix A.) By the trnsformtion lws, the errors of the m-th order vribles ffect the originl vribles s shown by the next proposition. Proposition 2B The difference between A (0) obtined from A (m) by the trnsformtion lws nd A (0) tr obtined from A (m) tr by the sme trnsformtion lws hs upper bound I (0) E( 2 )ǫ m, (5.54) σ (0) ǫ m, (5.55) θ (0) 2 ǫ m 1, (5.56) nd the corresponding difference s for the perturbtion vribles hs upper bound δi (0) E( 2 )ǫ m 1 δa (m) 1 (1), (5.57) δσ (0) ǫ m 1 δa (m) 1 (1), (5.58) δθ (0) 2 ǫ m 2 δa (m) 1 (1)+ǫ m 1 exp( 2 ǫ m ) [ 2 ǫ δθ (m) (1) + 4 δa m (m) (1)], (5.59) under the sme initil condition s in the previous proposition. (For the proof, see Appendix A.) According to the bove proposition, we cn conclude tht we cn mke the trunction errors s for the originl vribles rbitrrily smll if we truncte the system t the rbitrily lrge m-th order system. From Proposition 1, we cn see tht the prt independent of the ngle vribles of the evolution equtions re shifted fter the trnsformtions reducing the prt dependent on the ngle vribles. By the trunction, our system become much simpler thn the originl system. But it is still difficult to solve the truncted evolution equtions with such correction terms becuse the evolution equtions re complictedly entngled with ech other. In order to solve the evolution equtions nlyticlly, we wnt to discrd such correction terms. The errors produced by discrding such corrections re evluted in the following proposition. Proposition 3 The difference between the system with the correction terms produced by the trnsformtion nd the system obtined by discrding such correction terms is evluted in the following wy. As for the bckground vribles, the discrd errors re evluted s I E( 2 )ǫ, (5.60) σ ǫ, (5.61) θ 2, (5.62) 22

24 nd s for the perturbtion vribles, the discrd errors re evluted by δi E( 2 )ǫδb(1), (5.63) δσ ǫδb(1), (5.64) δθ 2 δb(1), (5.65) where under the initil conditions δb(1) := δi(1) + δσ(1), (5.66) I(1) = σ(1) = θ(1) = 0, (5.67) δi(1) = δσ(1) = δθ(1) = 0. (5.68) (For the proof, see Appendix A.) By the bove propositions, it cn be understood tht we cn obtin the informtion of the originl system with sufficiently good ccurcy by investigting the evolution equtions by simply dropping the prt dependent on the ngle vribles, becuse the errors produced by dropping re sufficiently smll nd in prticulr σ, δσ re bounded. The reson why these errors re mild is tht the finl stte in which ll the energy of the sclr fields is completely trnsferred into tht of rdition fluids is the ttrcting equilibrium round which the perturbtions do not grow. 6 Appliction of the LWL formul to the multicomponent reheting model In thissection, wepply thelwl formul tothereheting where theenergy ofthemultiple sclr fields is trnsferred into tht of the multiple rdition fluids. The decy rte from the sclr field φ to the rdition fluid ρ α is given by Γ α. When the interctions between the sclr fields φ : U int re neglected, the bckground quntities re solved s m I = A exp { γγ d 1 }, (6.1) 1 ρ 1/2 σ α = B α + d γ Γ ρ 1/2 α m I, (6.2) 1 where A, B α re integrtion constnts. As long s we do not give the expression of the totl energy density ρ in the integrls, the bove solutions do not give ny physicl informtion of reheting. But it is difficult to solve the evolution equtions in the form where the exct expression of ρ is explicitly described, becuse the evolution equtions of I, σ α re complicted nd highly nonliner. Then we expect tht the contribution to the integrtions owes minly to the period when the energy of the sclr fields is dominnt, tht is ρ cn be pproximted s ρ = A + B 3 4, (6.3) 23

25 where A := A, B := α B α. (6.4) We ssume tht the initil rdition energy density B is negligibly smll B A. We ssume tht ll Γ := α Γ α ( = 1,2,,N S ) re of the sme order of mgnitude. By substituting the bove ρ expression to the solutions (6.1)(6.2), by expnding the solutions with respect to B round B = 0, we obtin { m I = A exp 2 γγ 3 A 1/23/2 } +γγ A B A 3/21/2 exp { 2 γγ 3 A 1/23/2 } +O(B 2 ), 4 ρ r = 4 α ρ α where = B {1+ 32 G(2) 12 } G(1) + Γ := α ( ) 2/3 3 G(5/3) A1/3 2 γ 2/3 A +O(B 2 ), (6.5) Γ 2/3 Γ α, (6.6) nd G(t) is the Gmm function defined by G(t) := 0 dxx t 1 e x, (6.7) which is convergent for t > 0. By tking the exctly homogeneous perturbtion of (6.5) defined by ( ) D := δa A +δb, (6.8) B B=0 we cn obtin the Brdeen prmeter in the finl stte : ζ fin 1Dρ = r 4 ρ r = 1 ( ) 2/3 2 γ 2/ { 1 3 δa A A 1/3δB A Γ 2/3 { G(2) 1 2 G(1) } / + δa Γ 2/3 } / A Γ 2/3 A Γ 2/3 G(5/3). (6.9) From now on, we nme the expressions representing the finl mplitude of the Brdeen prmeter ζ fin in terms of the initil perturbtion mplitudes such s δa δb α, S formul fter S mtrix in the quntum mechnics. As for the initil energy density perturbtions of the sclr fields δa, by dopting more physicl prmetriztion introduced by δa cn be written s 1δA 3 A =: ξ, S b = δa δa b =: η b, (6.10) A A b δa = 3A ξ + b 24 A A b A η b. (6.11)

26 ξ represents the dibtic growing mode nd η b represent the isocurvture modes. Since η b stisfy η b = η b, η b +η bc = η c, (6.12) theindependent quntitiesregivenbyξ η 12 η 23 η NS 1N S. Byusingthisprmetriztion, δa dependent prt of ζ fin is written s ( ζ fin ξ b 1 Γ 2/3 ) 1 A A b / Γ 2/3 A η b b A Γ 2/3, (6.13) where A B implies tht B is contined by A, tht is A = B +. From this S formul, we cn conclude tht for the dibtic growing mode ξ, the Brdeen prmeter is conserved, nd tht the initil entropy perturbtions survive in the cse tht the decy rtes re dependent on the sclr field φ from which the rdition energy comes, tht is Γ Γ b ( b). In the cse where multiple sclr fields exist, there is no reson why the perturbtion hs only the dibtic component, nd it is nturl to think tht in the perturbtion the dibtic components nd the entropic components coexist. In such mixed cses, so clled conservtion of the Brdeen prmeter does not hold nd the bove S formul gives useful tool for clculting the finl Brdeen prmeter. We cn consider the cse where the energy trnsfer rtes Γ α fluctutes [29 ], which is clled s the modulted reheting scenrio [4 ]. For simplicity, we consider the one sclr field cse. By tking the derivtive of (6.5) with respect to Γ with B vnishing, we obtin ζ fin 1 δγ 6 Γ, (6.14) which is well known formul derived in the pper [4 ], nd where the coefficient is successfully determined in this pper. We consider the influence of the resonnt interction between sclr fields on the finl mplitudeofthebrdeenprmeter. WetketheinterctionU int intoccount byitertion. As the first order correction from the interction term s we obtin 1 ǫ σ α where (m I ) := exp { γγ d 1 } d m {γγ 1 ρ 1/2 1 exp 1 3 ρ 1/2 U int θ d d I, (6.15) d γ Γ ρ 1/2 α (m I ), (6.16) d 1 1 ρ 1/2 Since s the zeroth order pproximtion I obeys (6.1), we cn write U int = U { int =1 1 exp γγ d 1 θ θ 3n/2 1 ρ 1/2 }( 1 ǫ ) 3 Uint. ρ 1/2 θ (6.17) }, (6.18) 25

27 where we ssumed tht U int contins n n-th order interction term nd tht Γ is the pproprite sum of Γ /2. By substituting the bove expression nd by expnding the correction term with respect to B round B = 0, we obtin σ [ Γ α m U int θ =1 ( 1 ǫ G(5/3, n+3,γ,γ)+ )( 2 3 γ A 1/2 ( γ ) 2/3 B A 1/2 A ) n 8/3 1 γ { (1 ) 1/3 (Γ Γ )G(5/3, n+10/3,γ,γ) 2 1 ( ) 2/3 2 G(5/3, n+7/3,γ,γ) ( ) 2/3 2 G(1, n+3,γ,γ)}] 2 3 ( 3 2 ) 1/3 Γ G(2, n+3,γ,γ) (6.19) where G(n 1,n 2,Γ 1,Γ 2 ) := dxx n1 1 exp( Γ 1 x) x 0 x x 0 dyy n 2 1 exp{(γ 1 Γ 2 )y}, (6.20) where x := γ 2 A 1/2 3 3/2, x 0 := γ 2 (6.21) A 1/2 3 The evlution of the double Gmm function defined by (6.20) is treted in Appendix B. We consider the concrete exmple defined by In this cse, we use the independent vribles defined by U int = λφ 2 1 φ2 2, m 1 = m 2. (6.22) θ 1 = q 0 θ 2 = q 0 +q 1 I 1 = p 0 p 1 I 2 = p 1, (6.23) where (q 0,p 0 ) nd (q 1,p 1 ) re clled fst nd slow ction-ngle vribles, respectively [10 ]. Becuse of the resonnt reltion stisfied by the msses of the sclr fields, the slow ngle vrible q 1 moves much more slowly thn the fst ngle vrible q 0. The verging over the slow ngle vrible q 1 cnnot be justified, nd therefore the slow ction-ngle vribles (q 1,p 1 ) cn hve evolutions. In the previous pper [10 ], we investigted the influences of the resonnt interction on the evolution of the cosmologicl perturbtions before the energy trnsfer from the sclr fields to the rdition fluids begins. According to this study, the slow ction-ngle vribles cn suffer from the instbility ner the hyperbolic fixed point in the phse spce of the slow ction-ngle vribles. Since the initil dibtic perturbtion ξ nd the initil isocurvture perturbtion η 12 re given by nd ξ = 1 δp 0 (6.24) 3 p 0 η 12 = p 1δp 0 p 0 δp 1 (p 0 p 1 )p 1, (6.25) 26