arxiv:astro-ph/ v2 16 May 2005

Transcription

1 Growth of perturbtions in drk mtter coupled with quintessence Tomi Koivisto Helsinki Institute of Physics,FIN-4 Helsinki, Finlnd nd Deprtment of Physics, University of Oslo, N-36 Oslo, Norwy (Dted: September, 28) rxiv:stro-ph/5457v2 6 My 25 We consider the evolution of liner perturbtions in models with nonminiml coupling between drk mtter nd sclr field drk energy. Growth of mtter inhomogeneities in two exmples of such models proposed in the literture re investigted in detil. Both of these models re bsed on low-energy limit of effective string theory ction, nd hve been previously shown to nturlly led to lte ccelertion of the universe. However, we find tht these models cn be ruled out by tking properly into ccount the impct of the sclr field coupling on the formtion of structure in the drk mtter density. In prticulr, when the trnsition to ccelertion in these models begins, the interction with drk energy enchnces the smll scle clustering in drk mtter much too strongly. We discuss lso the role of n effective smll scle sound speed squred nd the issue of dibtic initil conditions in models with coupled drk sector. I. INTRODUCTION Cosmologicl observtions suggest tht bout ninetysix percent of the present energy density in the universe resides in forms unknown to the stndrd model of prticle physics. Mesurements of the fluctutions in the cosmic microwve bckground(cmb) rdition re consistent with nerly criticl density universe[], s predicted by the erly infltion. Assuming the vlidity of generl reltivity, roughly two thirds of this density nowdys must be due to negtivepressure component clled drk energy[2](de) to explin the ccelertion of the universe indicted by distnce-redshift curves of the high redshift supernove[3]. The rest of the density is in grvittionlly clustered mtter with negligible pressure on cosmologicl scles. However, only frction of this mtter consists of the (prtly luminous) bryonic mtter known to the stndrd model. Drk mtter (DM) is needed for consistency with convincing number of strophysicl observtions, for exmple the flt rottion curves of glxies[4]. It hs been proposed tht the two drk components do not evolve independently but re non-minimlly coupled. Due to their different pressures, the densities of these components decy with different rtes. This rises the disturbing question why they pper to be of the sme order of mgnitude right now. An explicit coupling between DE nd DM[5] could provide phenomenologicl solution to this so clled cosmic coincidence problem [6, 7, 8, 9]. On the other hnd, there is theoreticl motivtion[] to couple quintessence fields[, 2] to DM (see, however[3]). In prticulr, if it is held tht sclr-tensor theories should be described in the Einstein frme, it follows tht sclr fields with nonminiml coupling to other mtter re generl predictions of ll dimensionlly reduced theories[4, 5]. In the context of string theory, the dilton could be such field, nd indeed hs been proposed s cndidte for DE [6, 7]. Since generlly the dilton coupling to different mtter species could be different[8] nd such couplings to bryonic mtter re tightly constrined by high precision solr system tests of generl reltivity[9], it is of interest to concentrte on simpler models where the dilton is only coupled to DM, while the possible coupling to bryons is ssumed negligible. This pper concerns the liner perturbtions in such models. The purpose is to find how the formtion of cosmologicl structure in pressureless, dust-like mtter is influenced by coupling to sclr field, when this coupling is significnt lredy t the bckground order s required if it is to ddress the coincidence problem. First the cosmologicl field equtions for sclr field intercting with DM re derived to liner order in section II. These equtions hve been presented in the literture before (for exmple, see [2]), nd the reder fmilir with them my skip this section. However, there we hope lso to somewht clrify the reltions of different possible formultions nd theoreticl frmeworks for the coupled drk sector. In the two following sections we investigte specific cses. The simplest models of intercting DM nd DE cn be excluded by comprison to the CMB observtions[2], but more sophisticted scenrios exist[, 6, 7]. In section III we study model feturing two DM components, only the subdominnt of them intercting with quintessence[22]. A model where the DE-DM coupling is given by polynomil function of the DE field[7] is considered in section IV. These two models re bsed on low-energy limit of effective string theory ction nd hve been previously shown to nturlly led to lte-time ccelertion of the universe. Electronic ddress: tomikoiv@pcu.helsinki.fi

2 Though the detils of the mechnisms re different, the DE domintion in ech cse is triggered by the time-evolving coupling. Lrge scle structure in both of the two scenrios hve previously been climed to fit observtions s well s or better thn the stndrd ΛCDM model. However, here closer inspection shows tht these models re ruled out becuse of the excessive DM structure they predict in the lte-time universe. In subsection IVB we consider the dibtic initil conditions in for sclr field coupled to DM. In section V we discuss the reltion of our results to other models nd conclude. 2 II. EQUATIONS We consider n ction of the form A. Field equtions nd Binchi identities S = d 4 x g[ κ2 2 R 2 µ φ µ φ V(φ)+γ χ (φ)l (χ) +L (m) ], () where χ stnds for n intercting mtter species, nd m for non-intercting ones. The constnt κ (8πG) /2. We use the definition T µν (i) 2 δ( gl i ) g δ(g µν, (2) ) nd the Lgrngin of the field φ is tken to consist of the two first terms in the ction integrl Eq.(). The Einstein equtions resulting from this ction cn be written s nd the eqution of motion for the field φ is The Binchi identities then imply G µν = κ 2(T(φ) µν +γ χ (φ)t µν (χ) +T µν (m) ), (3) φ = V (φ) γ χ (φ)l(χ). (4) µ T ν µ (χ) = (δ ν µ L(χ) T ν µ µ T ν µ (m) =, (5) (χ) ) γ χ (φ) γ χ (φ) µφ. (6) Thus the energy momentum tensor for the intercting components re not seprtely conserved. From now on we will ssume tht L (χ) = gµν ẋ µ ẋ ν m δ(x x (t)), (7) g where ẋ i, i.e. the intercting component behves s dust. This is in fct the cse most commonly studied in cosmology(formodelswithdecoupledtoneutrinos,seerefs.[23, 24]). Onenotesthtthetermγ χ (φ)l (χ) inthection Eq.() cnbe nturllyinterpretedslgrnginforprticleswith field-dependent msses, m (φ) = m γ χ (φ). Such Lgrngin cn follow from conforml trnsformtion. Consider sclr-tensor theory defined by grvittionl Lgrngin proportionl to R/ γ(φ). Then the metric in the Einstein frme is ĝ µν g µν / γ(φ). In terms of this metric, the ction for the Lgrngin Eq.(7) reds S χ = d 4 x ĝγ(φ)ˆl (χ). (8) (At this point one my hve to redefine the field in order to hve the cnonicl kinetic term for the sclr field s in the ction Eq.().) A more stndrd wy, however, to describe sclr-tensor theories in the Einstein frme is to define the component T (X) µν = γ χ (φ)t (χ) µν. (9)

3 This energy momentum tensor then depends on the field φ, which turns out to be coupled to the trce T (X). For some purposes the χ-formultion is more convenient, since there the divergence in Eq.(6) vnishes except for the sptil prt in the liner order. However, we will lso use the X-formultion for some clcultions, especilly in subsection IIIB. With the vrible mss prticle interprettion, the difference is simply tht ρ χ is proportionl to the number density of the DM mtter prticles, while ρ X is their grvitting energy density. The fct tht the ction Eq.() is equivlent to the cse of trce-coupled sclr field nd the possibility to mke direct contct with conformlly trnsformed sclr tensor nd vrible mss prticle theories re independent of the form of the coupling. However, it is required tht the Lgrngin is effectively of the form Eq.(7). From tht follows (when ẋ i ) These simple identities would not, of course, hold for more generl Lgrngin. L (χ) = T (χ) = ρ χ. () 3 B. Cosmologicl bckground nd perturbtions We consider flt FRW universe. Including smll perturbtions bout the bckground, we cn write the line element s ds 2 = 2 (τ)[ dτ 2 +(δ ij +h ij )dx i dx j ]. () This prmeteristion of the metric perturbtions is clled the synchronous guge. There we cn chrcterize the sclr modes of h ij in the Fourier spce s h ij (x,τ) = d 3 ke [ˆk ik x iˆkj h(k,τ)+(ˆk iˆkj ] 3 δ ij)6η(k,τ), (2) when k = kˆk. In the following, prime is derivtive with respect to φ, n overdot represents derivtive with respect to conforml time τ, nd H is the conforml Hubble prmeter H ȧ/. These nd mostly ll other nottions nd conventions re those of M nd Bertschinger[25]. With these, the Friedmnn eqution reds H 2 = 2 3κ 2 ( φ V(φ)+γ χρ χ +ρ m while evolution of the sclr field is governed by the eqution ), (3) φ+2h φ+ 2 V(φ) = 2 γ χ (φ)ρ χ, (4) following from Eq.(4). In ddition to the intercting species nd the sclr field, we include stndrd CDM, bryons nd rdition in the universe. Then the energy densities scle s usully, nd therefore, with the scle fctor normlized to unity tody, we my write ρ χ = ρ χ 3, ρ c = ρ c 3, ρ b = ρ b 3, ρ r = ρ r 4. (5) Here superscript denotes the density tody. However, insted of ρ χ the grvitting energy density is now the density ρ X = γ χ ρ χ from definition Eq.(9), s implied by Eq.(3). The equtions for the perturbed metric re k 2 η 2 Hḣ = 2 2κ 2δT, (6) k 2 η = 2 2κ 2ik iδt i, (7) ḧ+2hḣ 2k2 η = 2 κ 2δTi i, (8)

4 ḧ+6 η +2H(h +6 η) 2k 2 η = 32 κ 2 (ˆk iˆkj 3 δ ij)σ i j, (9) where the components of the perturbed stress energy tensor re δt = ρ r δ r ρ χ (γ χϕ+γ χ δ χ ) ρ b δ b ρ c δ c 2 φ ϕ V ϕ, (2) 4 ik i δt i = 4 3 ρ rδ r θ r +γ χ ρ χ θ χ +ρ b θ b + φk 2 ϕ, (2) δt i i = ρ rδ r +3c 2 b ρ bδ b φ ϕ 3V ϕ. (22) We hve grouped together photons nd neutrinos with subscript r, since their contribution is s usul, nd used ϕ to denote δφ/φ s in[26]. In ddition to the field equtions, we need the equtions of motion for the individul components. For the perturbtions of the coupled sclr field, we derive the Klein-Gordon eqution from Eq.(4) to find ϕ+2h ϕ+( 2 V +k 2 )ϕ+ 2 φḣ = 2 ρ χ (γ χ δ χ +γ χ ϕ). (23) For the χ-component, the Binchi identities Eq.(6) in the liner order re δ χ = θ χ (24) 2ḣ, θ χ = Hθ χ + γ χ γ χ (k 2 ϕ φθ χ ). (25) Conservtion equtions for uncoupled mtter sty unmodified from the ones given in[25]. III. MODEL WITH TWO COLD DARK MATTER SPECIES A. Bckground solution In this section we investigte the model proposed by Huey nd Wndelt[22]. In this model quintessence with n exponentil potentil, is exponentilly coupled to χ-mtter V(φ) = V e καφ, (26) γ χ (φ) = e κβ(φ φc). (27) There is lso present n uncoupled species of cold DM. At erly times, the field behves s stndrd exponentil quintessence, but when it reches the minimum of its effective potentil, V eff = V +γ χ ρ χ, it evolves to potentildominted stge driving the ccelertion of the universe. There is motivtion for the dditionl complexity of introducing two kinds of DM. It is well known tht with minimlly coupled exponentil quintessence the trcking solutions do not yield n ccelerting phse[27]. If the field is exponentilly coupled to mtter, there exists n ccelerting ttrctor[28], but obviously the field could not hve reched it until tody lest there would not hve been mtter dominted er t ll. Moreover, it is known tht in the such models the coupling between DE nd DM is tightly constrined by the CMB spectrum[29] nd other cosmologicl dt[3]. The mount of quintessence during BBN must be smll, constrining α > 5 [27]. The onset of ccelertion should occur fter mtter-rdition equlity, which is chieved when κφ c > 8/α. Finlly, consistent lte-time bckground evolution requires α/β < 3/5 [22]. As representtive prmeter combintion we will use α = 5, β = 5, φ c = 25 in numericl clcultions. Fig.() shows then the evolution of the frctionl energy densities Ω i κ 2 ρ i /(3H 2 ). We hve rrnged the intercting nd stndrd CDM components to be of equl importnce tody. We see some trnsient

5 oscilltions t the lte epoch when the coupling term begins to hold the field t the minimum of its effective potentil. However, symptoticlly the coupling will keep the rtio of quintessence nd intercting DM densities constnt nd the field evolution is gin chrcterized with constnt scling. Consider the component ρ X γ χ ρ χ. Since during the trcking phse, ssuming constnt bckground eqution of stte w B, κφ = 3(+w B) α log() + constnt, (28) then ρ X 3[ + β α (+wb)]. Thus this energy density is growing in phntom-like mnner. Its effective eqution-of stte prmeter, defined by ρ X +3Hρ X (+w (eff) X ) =, (29) is shown together with the totl eqution of stte of the universe in the left pnel of Fig.(2). Similrly we define the effective eqution of stte for the field φ: Then we hve from Eq.(4) tht We plot both w φ nd w (eff) φ ρ φ +3Hρ φ (+w (eff) φ ) =. (3) ρ φ +3H(ρ φ +p φ ) = γ χ(φ)ρ χ φ w (eff) φ in the right pnel of Fig.(2). = p φ +γ χρ χ φ/(3h) ρ φ. (3) Ω FIG. : Evolution of the frctionl energy densities in model of section III. The solid line is Ω φ, the dshed line Ω X, the dotted line is Ω r nd the dsh-dotted line (Ω c +Ω b ).

6 6 w 2 3 w φ FIG. 2: Left pnel: Totl eqution of stte w (solid line) nd the effective w (eff) X (dshed line) for intercting CDM in the model of section III. Since w < /3 tody, the universe ccelertes. Right pnel: Eqution of stte for the quintessence in the sme model. Solid line is the forml w φ, nd dshed line is the effective w (eff) φ. B. Perturbtion evolution Defining C log(γ χ ), the coupled perturbtion equtions for the component X cn be written s[2] δ X = θ X 2ḣ+C ϕ+c φϕ, (32) θ X = Hθ X +C (k 2 ϕ φθ X ). (33) For the lrgest scles the source term in Eq.(33) cn be neglected, nd thus no velocity perturbtions re generted. Since then δ c δ X = C ϕ, the coupling hs smll effect on δ X, which then follows closely the stndrd CDM overdensities. To find the smll scle evolution of perturbtions, we follow derivtions found in Refs.[2, 28, 3, 32]. Using liner combintion of Eqs.(6) nd (8), ḧ+hḣ = 2 κ 2 [ δt δt i i], (34) nd verging over the oscillting time derivtives ϕ, ϕ nd δ r, one cn combine Eqs. (32) nd (33) into δ +H δ 2 2κ 2ρ Xδ X C φθ = [ C k 2 +HC φ+ C 2 κ 2V ]ϕ+ 2 2κ 2(ρ cδ c +ρ b δ b ), (35) when c 2 b. Averging out the smll nd oscilltory ϕ nd ϕ lso from Eq.(23) yields ϕ = 2 ρ X C δ X + φḣ 2 2 (V +ρ X C )+k2. (36) This pproximtion is in fct resonbly good for ll scles. We hve tested it by compring the LHS nd RHS of this expression evluted by numericlly integrting the full perturbtion equtions[47]. Inserting the pproximtion in Eq.(35) nd using Eq.(33), we get in the smll scle limit ( δ X + H +C φ ) δx = 3 2 H2 [Ω X (+2 C 2 κ 2 )δ X +Ω c δ c +Ω b δ b ]. (37) Here we hve used C =. During trcking er, Ω X 3[(+wB)( β α +) ] vnishes to good pproximtion. If we ssume tht then the nonintercting mtter perturbtion scles s δ c τ m 2, we cn use the trcker solution Eq.(28)

7 to find tht [ δ X + + 3β ] α (+w 2 B) (+3w B )τ δ X = Aτ m, (38) with the constnt A determined by the mgnitude of nonintercting mtter perturbtions t given time. The full solution to this eqution is δ X = D +Bτ b + A (b+m+)(m+2) τm+2, (39) where B nd D re integrtion constnts nd b = [2+6β(+w B )/α]/(+3w B ). Thus the lst term gives the dominnt mode, which hs the sme scling s the nonintercting DM. Therefore the rtio δ X /δ c is constnt during trcking er. During mtter domintion, when m = 2, this mens tht δ X, nd δ X /δ c = 3( 3/α 2 )/( + 2β/α). When the field reches the minimum of its effective potentil, it evolves with the rte[22] κ φ = 3 H. (4) α+β The criticl period we re interested in is the trnsition to the joint domintion of quintessence nd intercting mtter tht is ongoing tody. The field derivtive is then oscillting bout the vlue given by Eq.(4), to which it will symptoticlly settle (see Fig.(2)). Since Ω X is now pproching Ω c, nd (+C 2 ), we pproximte Eq.(37) by dropping the source terms due to nonintercting mtter. This pproximtion is well justified since, s it turns out, the RHS of Eq.(37) is soon dominted by the rpidly growing perturbtion in the density of intercting DM. We therefore write ( δ X + + 3β ) H δ X = 3 ( α+β 2 H2 Ω X +2β 2 ) δ X. (4) To solve this eqution nlyticlly, we would need expressions for the bckground quntities H (which follows from w B ) nd Ω X. To tht end we use the symptotic vlue lso for the ltter Inserting these in Eq.(4) yields the solutions w B = β/(α+β), Ω X = α/(α+β). (42) δ n, n = (α+β) 4 (α+β) [ ± + 24α(α+β)(+2β2 ) (α+β) 2 ] 7. (43) The exponent of the growing solution (corresponding to the positive sign) is lrge for typicl vlues in this model. For α = 5 nd β = 5, n.2. This is in strong contrst to stndrd DM perturbtions growth rte of which is slowing fter the mtter domintion when δ c. Therefore θ X δ c δ X r, r = n 2β 3α 2(α+β), (44) which is, consistently, the growing solution of Eq.(33) in this limit. Numericl results re shown in Fig.(3)[48]. There we see tht the intercting DM density contrst dpts to this symptotic scling in fct s the trnsition to ccelertion begins, lthough the solution Eqs.(4),(42) hs not been reched. We conclude tht when the energy in the χ-mtter is becoming significnt contributor to the bckground expnsion, the perturbtion growth in this component is lmost explosively sped up by the interction pumping energy into it from the φ-field. As seen from Eq.(36), the DE is then inhomogenenous on smll scles. It hs been found lso tht in the nonliner regime DE tends to clump with DM if coupled[33, 34, 35]. The lrge coupling for which C = β in this model results in lrge effective grvittionl constnt in Eq.(4), G eff = ( + 2β 2 )G. In erlier cosmology this is not s importnt, since then the stndrd DM overdensities constitute the dominnt grvittionl sources. There the modified Hubble friction term in Eq.(37) cts to reduce the reltive mgnitude of δ X to δ c, but these evolve with the sme rte. However, in the lte universe the cosmic structure is drsticlly incresed, except t the lrgest scles where we see only little oscilltions. Consequently, the shpe of the mtter power spectrum in this model is highly incomptible with observtions. Also, s the increse in structure will cuse rpid deepening of the grvittionl wells, the resulting CMB spectrum is modified by strong lte integrted Schs-Wolfe effect due to time vrition of the grvittionl potentils. We hve confirmed these considertions by numericl integrtion for severl prmeter vlues stisfying the reltions given in the previous subsection. As non-negligible mount of χ-mtter nd reltively strong coupling re required in this model, we feel it sfe to stte lthough hving not performed n extensive explortion in the spce of ll cosmologicl prmeters, tht this model cnnot be succesfully confronted with observtions.

8 8 5 δ FIG. 3: Intercting DM overdensities in the model of section III. The slope s of the curves in the logrithmic plot is s 2,,.2 erly in the rdition domintion (when the scle is outside the horizon), during mtter domintion nd during the trnsition to ccelerting universe, respectively. The k-vlues from bottom to top re k = 2.4 6, 2.4 5, 2.4 4,.24,.24 nd.24 Mpc. IV. MODEL WITH POLYNOMIAL COUPLING A. Bckground solution The model proposed by Ben nd Mgueijo[7] is nother exmple of modified exponentil quintessence scenrio. Agin field is coupled to DM in such wy tht quintessence is dominnt tody nd driving the ccelertion of the universe while the model still hs its prmeters of the order of the Plnck scle nd its outcome widely independent of the initil conditions for the field. This is chieved with only one DM species by using non-exponentil form for the coupling. To be definite, we hve in the ction Eq.() V(φ) = V e καφ, γ χ (φ) = +λ[κ(φ φ )] b. (45) For numericl vlues we will use α = 8, b = 8, λ = 5 nd κφ = The bckground evolution is shown is Fig.(4). In this model the field leves the stndrd exponentil trck long before mtter domintion. After initil trnsients, the field trcks the bckground rdition, until the coupling term pushes the field into kintion. This lsts until the field reches the minimum of its effective potentil. Then it will t due time drive the ccelertion. However, s the mtter density eventully dilutes, the coupling force will fde, leding the field bck to stndrd scling in the symptotic future. Thus mtter domintion is restored nd the ccelertion in this model is trnsient. We plot the effective eqution of stte prmeter for the field φ, s defined in Eq.(3), in the left pnel of Fig.(5). In the right pnel we plot the totl eqution of stte for the universe together with the w (eff) X s defined in Eq.(29). Since in this model ρ X governs the bckground expnsion lso in the mtter dominted er, it is importnt tht the coupling term hs then negligible effect on the scling of this density, since significnt devition from τ 2 would result in modified conforml distnce to the lst scttering surfce, possibly dngerous to the comptibility with observtions of the peks in the fluctution spectrum of the CMB rdition[2]. One lso expects tht then the growth rte of DM perturbtions will not be much modified. However, one would guess tht the nonstndrd evolution during trnsition to DE domintion, seen in Figs.(4) nd (5), hs some effects on the perturbtion growth s ws the cse in the model of section III. Indeed these effects prove ftl to this model, s will be seen below.

9 9.8.6 Ω FIG. 4: Evolution of the frctionl energy densities in the model of section IV. The solid line is Ω φ, the dshed line Ω X. Dotted line is Ω r nd the dsh-dotted line Ω b w w φ FIG. 5: Left pnel: Totl eqution of stte w (solid line) nd the effective w (eff) X (dshed line) for intercting DM in the model of section IV. Since w < /3 tody, the universe ccelertes. Right pnel: Eqution of stte for the quintessence for the sme model. Solid line is the forml w φ, nd dshed line is the effective w (eff) φ. B. Perturbtion evolution The linerized Einstein equtions hve in fct been solved in this model[36]. However, this ws not done tking fully into ccount the impct of the coupling to the mtter perturbtions, but insted by ssuming tht the fluctutions in ρ χ decouple from those in φ. This ws relised by setting the DM velocity perturbtion to zero, s conventionl in the synchronous guge. However, now the interction with DE will drive θ χ to evolve even if its initil vlue is chosen

10 to be zero. One sees tht the source term in the Euler eqution Eq.(25) vnishes only in the bsence of coupling or if the DE is perfectly smooth. Neither is the cse now. We show the effect of physicl ssumptions in Fig.(6). Since now C =, Eq.(37) is not vlid. However, using the pproximtions tht led to tht eqution, we cn derive the corresponding evolution eqution for the quntity δ χ. It should cpture the qulittive behviour of the DM perturbtion, since δ X = δ χ +C ϕ δ χ for smll scles. The eqution we find reds ) δ χ + (H λbyb +λy bκ φ δ χ = 3 [ 2 H2 Ω X + λyb (λy b 2 (2b 2 y 2 ] ) ) (+λy b ) 2 δ χ, (46) where y κ(φ φ ). The effective grvittionl constnt G eff χ cn be thought to be the expression is squre brckets times G. We plot G (eff) /G s function of φ in Fig.(7) for choice of prmeters. The vlues of G eff becomes relevnt t the mtter dominted er, when the source term in the RHS of Eq.(46) is non-negligible. However, s the coupling hs driven φ ner to φ, we hve then G eff χ G (for our exmple prmeters,.5 < y < from = 4 until =.). This chnges when Ω φ climbs up, since then φ rolls forwrd. Similrly s in the previous cse, δ X is driven to rpid increse. We show its evolution in Fig.(8), where we lso plot the bryon overdensities δ b for comprison. The boost experienced by δ X is not s vst s we sw in the previous section, but still results in n order of mgnitude or two too lrge δ X. Agin we believe tht notble improvement is not found by djustments of prmeters, but by revision of the scenrio..5.4 C l l FIG. 6: The CMB spectrum in the model of section IV (solid line). Dsh-dotted line corresponds to suppressing the coupling of drk energy nd DM perturbtions (i.e. θ χ = ), nd dshed line corresponds to suppressing DE perturbtions (i.e. ϕ = ). Finlly, remrk bout the dibtic initil conditions we hve used in the clcultions. Adibticity of perturbtions[37] mens tht given ny two superhorizon regions A nd B in perturbed universe t n instnt τ, B is identicl to A t some other instnt τ +δτ. Evolution of the regions is the sme, but slightly synchronous. Then the rtio of fluctution nd time derivtive of sy density field is common to ll density fields. Thus we cn define the entropy between species i nd j s S ij δρ i δρ j ρ i ρ j = δ i +w (eff) i δ j +w (eff) j. (47)

11 2 5 G X (eff) /G κφ FIG. 7: The devition of the effective grvittionl constnt from the stndrd one, G (eff) /G, s function of the field κφ for the model of section IV. The solid curve corresponds to the G χ (eff) in Eq.(46), nd the dshed curve is 2C 2, corresponding to G eff X s it ppers in Eq.(37). For both definitions the devition is smll when φ φ, but becomes very rpidly lrge when the field φ evolves to lrger vlues. 5 δ FIG. 8: DM overdensities (solid curves) in the model of section IV. Dshed curves re δ b. The k-vlues from bottom to top re k = 2.4 6, 2.4 5, 2.4 4,.24,.24 nd.24 Mpc. The generl dibtic condition then dicttes tht ll S ij s together with their derivtives vnish. Note tht only in the bsence of interctions does this condition coincide with the usully cited δ i +w i = δ j +w j, θ i = θ j. (48) There cn be difference between these conditions, s seen for the cses we re now considering in Figs.(2) nd (5). For scle fr outside the horizon nd deep in the rdition dominted er, the growing dibtic mode[25, 38] of

12 the sclr field perturbtion stisfies 3φ 2 ( 2 V 6 γ χ ρχ) φτ + 2 ρ χ ϕ = (6 φ+γ χ 2 ρ χ τ)v γ χρ χ φ2 τ [ γ χ φ 2 γ χ 2 V γ 2 χ 2 ρ χ ] τ 2 2 φak 2 τ 2, (49) [ ϕ = (4 3 γ χρ χ 2 V V 2 2 γ χρ φ2 χ )τ + 2 ρ χ 2 γ χ(v + 3 γ χρ χ )+ 2 γ 2 V γ 2 χ 2 ρ χ V ] τ 2 (6 φ+γ χ 2 ρ χ τ)v γ χ ρ φ χ 2 τ φak 2 τ 2. (5) Normliztion A = /2 is convenient nd sets comoving curvture perturbtion equl to unity. Fortuntely the bove expressions cn in mny occsions be simplified by dropping the squre brcket terms, but we hve kept them here for the ske of completeness. The deprture of the usul dibtic condition Eq.(48) from the generlized one hs been pointed out lso in Refs.[39, 4]. If the field is trcker, the initil conditions my indeed be irrelevnt lso for the perturbed field. However, for minimlly coupled fields it hs been shown tht in some instnces this is not the cse[4, 42], while isocurvture perturbtions in coupled quintessence remin to be systemticlly studied. V. DISCUSSION Interestingly, in the two models we hve investigted, the growth of DM perturbtions is enchnced, wheres in minimlly coupled models the lte time effect of DE is the opposite. This my be understood by compring to the cse where the intercting DE nd DM stisfy the dibtic reltion S φχ = Ṡφχ =, where S φχ is defined by Eq.(47). Then the system of these components cn be treted s unified fluid without intrinsic entropy. It is known tht the effective sound speed squred c 2 S ṗ/ ρ of such fluid governs the lte dynmics of the perturbtions[43, 44]. If c 2 S >, they re driven to oscillte, if c2 S < they tend to blow up. This is the work of the pressure grdient sourcing the evolution eqution one cn derive for the unified fluid density contrst, δ +[ 3(2w c 2 S)]H δ 3 2 ( 6c2 S +8w 3w 2 )H 2 δ = k 2δp δ. (5) δρ This eqution is exct (to liner order, in the comoving guge). As w nd c 2 S re given by the bckground expnsion, one needs only informtion of the pressure perturbtion to solve the structure in the universe. In the dibtic cse, δp/δρ = c 2 S, but in the sclr field models considered here this ssumption is not vlid. However, using equtions from subsection IIB, the smll scle limit of the pressure grdient in the lte universe is redily found. It is k 2δp δρ = φ k2 ϕ V (φ)ϕ 2 ρ χ δ χ 2V (φ)ϕ 2 C (φ)v (φ), (52) (lso in the comoving guge) nd turns out positive for the models of section III nd IV. Thus one my regrd the enchnced smll scle growth rte there s consequence of negtive effective (non-dibtic) sound speed squred. Now the pressure grdient pproches k-independent expression nd the growth rte sturtes t given vlue when k is incresed beyond some criticl scle, unlike in the dibtic cse. Note lso tht the totl density perturbtion δ in Eq.(5), δ = ρ χγ χ (φ)δ χ +[ρ χ γ (φ)+v (φ)]ϕ+ φ ϕ 2 γ χ (φ)ρ χ + φ, (53) 2 +V(φ) 2 is directly relevnt to the integrted Schs-Wolfe effect in the CMB fluctution spectrum. Before concluding, we use the opportunity to comment on yet nother model with DE-DM interction[7]. There DE with constnt eqution of stte prmeter w φ is coupled to DM in such wy tht ρ φ +3Hρ φ = 3Hλ 2 (ρ X +ρ φ ), (54) ρ X +3Hρ X = 3Hλ 2 (ρ X +ρ φ ), (55) where λ 2 is constnt mesuring the intensity of the interction. It might be difficult to find field theory underlying this phenomenologicl model, but tht would be necessry to be ble to derive the corresponding equtions to the

13 liner order. Otherwise they must be put in by hnd. The liner DE continuity equtions presented in the recent pper[45] imply tht there µ T µ (φ) = 3Hρ φ λ 2 (δ φ +rδ X ), (56) 3 µ T µ (φ) i = [3Hρ φ λ 2 (+r)θ φ ],i, (57) where r ρ χ /ρ φ. Except perhps for the sign of the -component of the divergence[49], these expressions re not unresonble, but do not follow unmbiguously from Eqs.(54) nd (55). However, once the divergence of the DE stress tensor is given, the intercting DM continuity equtions re uniquely fixed. When Eqs.(56) nd (57) hold for the DE nd the DM Lgrngin is effectively of the form (7), the Binchi identities determine tht δ X = θ X 2ḣ 3Hλ2 ( r δ φ +δ X ), (58) θ X = Hθ X +3Hλ 2 (+ r )θ φ. (59) From these we recover the equtions used in Ref.[45] only by setting θ = in the first eqution nd forgetting bout the second one. But s detiled in the previous section, such pproximtion cn leve crucil fetures out from the evolution of perturbtions. Therefore we must regrd tht the results of the dt nlysis in Ref.[45], where it is nnounced tht observtions prefer intercting DE models, cll for reinvestigtion. As mentioned in the beginning of this pper, the origin of most of the energy density in the universe is unknown. Suggestions cn be bsed on specultions of string theory, but this is tenttive t the present. However, some inventive mechnism, s the ones considered here, is necessry if some light is required to be shed on the coincidence problem. Moreover, s the results of this pper demonstrte, considertions of cosmologicl structure provide powerful mens to constrin such mechnisms. Acknowledgments The uthor cknowledges useful discussions with Hnnu Kurki-Suonio, Dvid Mot nd Finn Rvndl. This work ws supported in prt by NorFA, Wldemr von Frenckels Stiftelse nd Emil Altosen Säätiö. [] D. N. Spergel et l. (WMAP), Astrophys. J. Suppl. 48, 75 (23), stro-ph/3229. [2] P. J. E. Peebles nd B. Rtr, Rev. Mod. Phys. 75, 559 (23), stro-ph/ [3] A. G. Riess et l. (Supernov Serch Tem), Astrophys. J. 67, 665 (24), stro-ph/4252. [4] G. Bertone, D. Hooper, nd J. Silk, Phys. Rept. 45, 279 (25), hep-ph/4475. [5] C. Wetterich, Astron. Astrophys. 3, 32 (995), hep-th/ [6] W. Zimdhl nd D. Pvon, Phys. Lett. B52, 33 (2), stro-ph/5479. [7] L. P. Chimento, A. S. Jkubi, D. Pvon, nd W. Zimdhl, Phys. Rev. D67, 8353 (23), stro-ph/3345. [8] R.-G. Ci nd A. Wng, JCAP 53, 2 (25), hep-th/425. [9] M. C. Bento, O. Bertolmi, nd A. A. Sen, Phys. Rev. D7, 8359 (24), stro-ph/ [] G. R. Frrr nd P. J. E. Peebles, Astrophys. J. 64, (24), stro-ph/3736. [] I. Zltev, L.-M. Wng, nd P. J. Steinhrdt, Phys. Rev. Lett. 82, 896 (999), stro-ph/9872. [2] R. R. Cldwell, R. Dve, nd P. J. Steinhrdt, Phys. Rev. Lett. 8, 582 (998), stro-ph/ [3] M. Dorn nd J. Jeckel, Phys. Rev. D66, 4359 (22), stro-ph/238. [4] A. Albrecht, C. P. Burgess, F. Rvndl, nd C. Skordis, Phys. Rev. D65, 2357 (22), stro-ph/7573. [5] P. Brx, C. vn de Bruck, nd A. C. Dvis, JCAP 4, 4 (24), stro-ph/ [6] M. Gsperini, F. Pizz, nd G. Venezino, Phys. Rev. D65, 2358 (22), gr-qc/86. [7] R. Ben nd J. Mgueijo, Phys. Lett. B57, 77 (2), stro-ph/799. [8] J. Grci-Bellido, Int. J. Mod. Phys. D2, 85 (993), hep-ph/ [9] T. Dmour nd A. M. Polykov, Nucl. Phys. B423, 532 (994), hep-th/9469. [2] L. Amendol, Phys. Rev. D69, 3524 (24), stro-ph/375. [2] M. B. Hoffmn (23), stro-ph/3735. [22] G. Huey nd B. D. Wndelt (24), stro-ph/4796. [23] R. D. Peccei, Phys. Rev. D7, (25), hep-ph/437.

14 [24] A. W. Brookfield, C. vn de Bruck, D. F. Mot, nd D. Tocchini-Vlentini (25), stro-ph/ [25] C.-P. M nd E. Bertschinger, Astrophys. J. 455, 7 (995), stro-ph/ [26] P. G. Ferreir nd M. Joyce, Phys. Rev. D58, 2353 (998), stro-ph/972. [27] E. J. Copelnd, A. R. Liddle, nd D. Wnds, Phys. Rev. D57, 4686 (998), gr-qc/9768. [28] L. Amendol, Phys. Rev. D62, 435 (2), stro-ph/ [29] L. Amendol, C. Quercellini, D. Tocchini-Vlentini, nd A. Psqui, Astrophys. J. 583, L53 (23), stro-ph/2597. [3] U. Frnc nd R. Rosenfeld, Phys. Rev. D69, 6357 (24), stro-ph/3849. [3] L. Amendol nd D. Tocchini-Vlentini, Phys. Rev. D66, (22), stro-ph/535. [32] L. Amendol nd D. Tocchini-Vlentini (22), stro-ph/ [33] D. F. Mot nd C. vn de Bruck, Astron. Astrophys. 42, 7 (24), stro-ph/454. [34] N. J. Nunes nd D. F. Mot (24), stro-ph/4948. [35] C. vn de Bruck nd D. F. Mot (25), stro-ph/5276. [36] R. Ben, Phys. Rev. D64, 2356 (2), stro-ph/4464. [37] A. R. Liddle nd D. H. Lyth, Cosmologicl infltion nd lrge-scle structure (Univ. Pr., 22), cmbridge, UK; 4 p. [38] M. Bucher, K. Moodley, nd N. Turok, Phys. Rev. D62, 8358 (2), stro-ph/ [39] K. A. Mlik, D. Wnds, nd C. Ungrelli, Phys. Rev. D67, 6356 (23), stro-ph/262. [4] K. A. Mlik nd D. Wnds, JCAP 52, 7 (25), stro-ph/473. [4] M. Kwski, T. Moroi, nd T. Tkhshi, Phys. Lett. B533, 294 (22), stro-ph/88. [42] L. R. W. Abrmo nd F. Finelli, Phys. Rev. D64, 8353 (2), stro-ph/4. [43] R. R. Reis, M. Mkler, nd I. Wg, Clss. Qunt. Grv. 22, 353 (25), stro-ph/563. [44] T. Koivisto, H. Kurki-Suonio, nd F. Rvndl, Phys. Rev. D7, 6427 (25), stro-ph/4963. [45] G. Olivres, F. Atrio-Brndel, nd D. Pvon, Phys. Rev. D7, (25), stro-ph/ [46] A. Lewis, A. Chllinor, nd A. Lsenby, Astrophys. J. 538, 473 (2), stro-ph/9977. [47] In numericl clcultions we use the public CAMB code[46], see For simplicity, we ssume scle-invrint initil power spectrum, no reioniztion nd three mssless neutrino species. In our clcultions, h =.7. Normliztion is such tht the primordil curvture perturbtion equls unity. [48] Since we hve normlized the primordil curvture perturbtion to unity, these perturbtions in fct exit the liner regime t bout δ 5. So these clcultions re untrustworthy beyond those vlues, but tht is unrelevnt to our conclusions. [49] The sign is different from the zero-order divergence. Possibly this is misprint. 4