arxiv: v1 [astro-ph.co] 16 Feb 2010

Transcription

1 Collpse of Smll-Scle Density Perturbtions during Preheting in Single Field Infltion Krsten Jedmzi Lbortoire de Physique Théorique et Astroprticules, UMR 507-CNRS, Université de Montpellier II, F Montpellier, Frnce rxiv: v1 [stro-ph.co] 16 Feb 010 Mrtin Lemoine nd Jérôme Mrtin Institut d Astrophysique de Pris, UMR 7095-CNRS, Université Pierre et Mrie Curie, 98bis boulevrd Argo, Pris, Frnce (Dted: August, 013) After cosmic infltion nd before the trnsition to rdition domintion, the cosmic energy density my hve been dominted during n extended period by n oscillting mssive sclr condenste. We show tht during this period, sub-hubble scle perturbtions re subject to metric preheting instbility in the nrrow resonnce regime. We nlyze in detil both, qudrtic nd qurtic potentils. The instbility leds to the growth of density perturbtions which in mny cses become non-liner lredy before the beginning of rdition dominted Universe. This is prticulrly the cse when requiring phenomenologiclly preferred low rehet temperture. These erly structures my led to the emission of grvittionl wves nd the production of primordil blc holes. Furthermore, it is not cler if they could modify the prediction of liner curvture perturbtions on very lrge scles. I. INTRODUCTION The Universe hving pssed through n erly infltionry epoch is by now widely ccepted prdigm for n explntion of pprent cosmic lrge-scle homogeneity s well s the perturbtions observed in the cosmic microwve bcground rdition (CMBR). One of the erliest incrntions of infltion ws given in Ref. [1] where mssivesclrfield with potentil V(φ) = m φ /, strting from super-plnc vlues, slowly rolls down towrds its minimum thereby driving period of ccelerted expnsion. Though mny lternte scenrios of infltion hve been discussed over the yers, the bove scenrio of chotic infltion remins one of the most successful in terms of reproducing the observed perturbtions []. In prticulr, the predicted spectrl index n S 0.96 in this model is consistent with tht inferred from the WMAP- 7 nlysis of the CMBR inhomogeneities, which gives n S 0.963±0.01 (68% CL) [3]. After infltion hs ended, cosmic energy density is dominted by n lmost homogeneous oscillting inflton field. A trnsition towrds rdition dominted Universe, s required to exist t the ltest shortly before the epoch of Big Bng nucleosynthesis (BBN) which strts t cosmic temperture T 1MeV, is usully ssumed to occur by the perturbtive decy of mssive inflton prticles into rdition. This epoch is referred to s the epoch of reheting. Prmetric resonnce phenomen, inducing energy trnsfer to other sclrs cou- Electronic ddress: jedmzi@lpt.univ-montp.fr Electronic ddress: lemoine@ip.fr Electronic ddress: jmrtin@ip.fr pled to the inflton, referred to s preheting, though rich in their physics, hve been shown to fil to provide complete trnsition towrds rdition dominted Universe (see, e.g. [4, 5] for review of preheting models). In cse the inflton would decy immeditely fter the end of chotic infltion, nd ssuming m φ / infltion, the cosmic rehet temperture would be exceedingly lrge, i.e. T rh GeV. Such lrge T rh is phenomenologiclly problemtic, s liely leding to the re-genertion of unwnted relics, such s mgnetic monopoles t GUT trnsition, or mssive grvitinos in supergrvity [6, 7]. Modern constrints from BBN on unstble grvitinos constrin this reheting temperture down to T rh 10 7 GeV [8 10]. It is thus possible (nd desirble) tht between the end of chotic infltion nd reheting cosmic energy density is dominted for very long expnsion period by n lmost homogeneous inflton condenste. For exmple, for T rh s low s 1MeV the Universe would hve expnded n immense fctor 10 5 between the end of infltion nd reheting. It is such cosmic er which is the topic of the current pper. An inflton oscilltory phse fter infltion hs received some ttention with respect to the evolution of super-hubble size perturbtions [11 13]. It hs been shown tht, though one finds prmetric resonnce in the m φ /cse, super-hubble size perturbtionsthoughtto be responsible for the nisotropies in the CMBR sty constntinmplitude, sisthecseinrditiondominted Universe. It ws therefore concluded tht in single-field models metric preheting does not led to non-trivil chnges in the evolution of perturbtions. However, s we will see in this pper (see lso Ref. [11]), prmetric resonnce occurs even for sub-hubble size regions, i.e. for ll modes with physicl wve vectors phys < 3Hm (for qudrtic potentil), where H is the Hubble constnt.

2 In fct, these modes lie in the nrrow resonnce regime nd the red-shifting is such tht ny mode which hd once been resonnt will lwys remin resonnt. Due to this resonnce the curvture perturbtions of sub-hubble modes sty constnt rther thn decresing. One cn redily show tht this implies tht the density contrst grows with scle fctor δρ /ρ implying non-liner cosmic structure formtion fter n expnsion fctor of This fct hs so fr received little ttention in the literture. Whether the Universe indeed entered n er with nonliner structures, lredy before the epoch of reheting, depends on the rehet temperture nd the scle of infltion. Attempting to eep the rehet temperture low enough to void overproduction of grvitinos, i.e. T rh < GeV [9, 10] implies tht essentilly ll lrge-field infltionry models re followed by phse of non-liner structure formtion. Figure 1 provides n illustrtion of this fct ssuming qudrtic potentil during the oscillting phse. It is seen tht even for H end s low s 10 7 GeV non-linerities typiclly develop unless T rh 10 8 GeV. Therefore, nonliner cosmic phse before reheting ppers to be generic prediction. Note, however, tht these considertions ssume tht fter infltion energy density is dominted by single oscillting sclr field. These findings my not be surprising since it is nown tht the effective eqution of stte during such phse is tht of Universe filled with pressure-less dust. One would thus expect the inflton condenste to be grvittionlly unstble, which indeed, hs been shown in the Newtonin limit [14, 15] nd in n verged reltivistic limit [16, 17]. However, our nlysis nd the effect tht we report go beyond this Newtonin nlogue, becuse this instbility is found to be intrinsiclly some form of metric preheting. And, indeed, we find tht similr effect lso tes plce in the cse of qurtic inflton potentil, for which the effective eqution of stte during oscilltions is tht of rdition. In this ltter cse, the resonnce is (nd remins) shrply peed round the physicl scle ssocited to the Hubble scle t the end of infltion, but the growth of the density contrst is exponentil. The outline of the pper is s follows. In Section II the metric preheting of sclr perturbtions is nlyzed in detil. The primordil power spectrum is clculted numericlly over wide rnge of scles. In Section III, we exmine in detil the growth of density perturbtions in qudrtic inflton potentil, while in Section IV, we discuss the cse of qurtic inflton potentil. Conclusions re drwn in Section V. II. GENERAL CONSIDERATIONS The evolution of sclr (density) perturbtions is controlled by the so-clled Muhnov vrible. If mtter is described by sclr field (s it is the cse during in- FIG. 1: Contour plot of the r.m.s. density fluctution δρ /ρ on the scle corresponding to the horizon t the end of infltion s mesured t the onset of the rdition er in the T rh (reheting temperture) H end (Hubble constnt t the end of infltion). The solid line tht is lbelled 1 delinetes prmeter spce where non-liner structures form between the end of infltion nd reheting (the region below this line) from prmeter spce where the Universe remins in the liner regime of density perturbtions, i.e. δρ /ρ 1, before reheting (the blue region bove the line). The normliztion for the curvture perturbtions on the Hubble scle t the end of infltion hs been ssumed to be P ζ in ccordnce with Fig.. The two dotted lines indicte where the power spectrum of the fluctutions is lrger nd smller thn one or, equivlently, how the solid line would move when the normliztion of the curvture perturbtions on the Hubble scle fter infltion is higher, P ζ = 10 9 (line lbeled 0.01), or lower, P ζ = (line lbeled 100). The solid red line shows pproximte upper limits on T rh in order to void grvitino overproduction. Note tht the figure is bsed on Eq. (5), nd mes the firly generic ssumptions tht energy density fter infltion is dominted by single oscillting sclr field with potentil m φ /, nd tht φ M Pl t the end of infltion. fltion nd preheting), then the eqution of motion of the Fourier trnsform of Muhnov s vrible is given by [18 1] v + [ ( ) ] ǫ1 v = 0. (1) ǫ 1 Inthisrticleprimedenotesderivtivewithrespectto conforml time while dot mens derivtive with respect to cosmic time. The quntity is the comoving wve numbernd ǫ 1 Ḣ/H is the first slow-rollprmeter. The quntity H = ȧ/ is the Hubble prmeter nd denotes the Friedmn-Lemître-Robertson-Wler scle fctor. Muhnov s vrible is relted to the curvture perturbtion, ζ = κ v, () ǫ 1

3 3 where κ 8π/m Pl, m Pl being the (non-reduced) Plnc mss. The importnce of the curvture perturbtion lies in the fct tht, on lrge scles, it is conserved quntity [] nd, therefore, it cn be used to propgte the infltionry spectrum from the end of infltion to the post-infltionry er. As consequence of the bove reltion, the spectrum of ζ cn be expressed s P ζ () 3 π ζ = 3 πm Pl v ǫ 1. (3) In order to clculte P ζ (), one needs to integrte Eq. (1), which requires the nowledge of the initil conditions for the mode function v. Since, t the beginning of infltion, ll the modes of strophysicl interest tody were much smller thn the Hubble rdius, the initil condition re chosen to be the Bunch-Dvis vcuum which mounts to lim v = 1 e iη, (4) /(H) + where η denotes conforml time. The bove formultion is well-dpted to the clcultion of the power spectrum during infltion. However, in order to study the metric perturbtions during the phse where the inflton field oscilltes it is more convenient to wor in terms of cosmic time nd to use the rescled vrible defined by ṽ = 1/ v. Then, Eq. (1) tes the form [13] [ ṽ + + d V dφ +3κ φ κ H φ 4 + 3κ 4 + κ φ ] dv H dφ ( ) φ V ṽ = 0. (5) This formultion is prticulrly convenient s ll the terms in the bove eqution remin non-singulr even when the field oscilltes t the bottom of its potentil, tht is to sy when φ, φ nd ǫ1 periodiclly vnish [1, 3]. To proceed further, one needs to choose the potentil V(φ). The form of the inflton potentil is presently unnown (lthough it is constrined, see Refs. [4 6]). A convenient wy to tret the most generl cse is to consider the following potentil V(φ) = V 0 ± 1 m φ + λ 4 φ4 +, (6) where the dots refer to higher order terms. As it is wellnown, this expnsion cn be difficult to control, most notbly so for lrge field models where, during infltion, the vcuum expecttion vlue of the field is lrge compred to the Plnc mss. However, one cn lwys ssume tht some terms in the bove expnsion re bsent, for instnce due to prticulr symmetries of the underlying theory. Different choices of coefficients in the bove eqution (6) correspond to different models. With the three terms displyed in Eq. (6), one cn describe lrge field models, smll field models nd hybrid infltion (t lest in the infltionry vlley). In the next section, we consider the mssive cse. III. THE QUADRATIC CASE A. Density Perturbtions In this section, we investigte model with simple qudrtic potentil, V(φ) = 1 m φ. (7) We hve in mind lrge field infltionry model, but s discussed bove, this qudrtic potentil hs more generl pplicbility when it comes to discussing the physics of inflton oscilltions fter the ctul er of infltion. It my lso describe the physics of curvton oscilltions. We will therefore eep wherever needed, the dependence of our results on the mss m (nd Hubble scle H end t the end of infltion), lthough we will use s fiducil vlues those corresponding to lrge field infltionry model. In lrge field models, the mss m is fixed by the normliztion of the mplitude of density fluctutions to tht mesured by COBE nd WMAP. Stndrd clcultions led to m 180π N +1 Q rms PS m Pl (8) T where N is the number of e-folds before the end of infltion t which the scle of strophysicl interest left the Hubble rdius during infltion (we te N = 50) nd where Q rms PS /T This gives m m Pl. The power spectrum P ζ () of density fluctutions t the end of infltion is clculted numericlly in the following, vi the integrtion of Eq. (1) throughout the infltionry er. Although one my find excellent pproximtions to this power spectrum on CMBR scles through the use of stndrd slow-roll clcultions, we indeed find tht up to second order in the slow-roll prmeters, these methods cnnot predict relibly the power spectrum on the shorter scles tht we re interested in. More specificlly, one finds to first order in the slow-roll prmeter [7]

4 4 [ ( )] P ζ () = H πǫ 1 m 1 (C +1)ǫ 1 Cǫ +( ǫ 1 ǫ )ln, (9) Pl while, t second order, it is given by [7] H { ) P ζ () = πǫ 1 m 1 (C +1)ǫ 1 Cǫ + (C +C + )ǫ (C π Pl C + 7π 1 7 ǫ 1 ǫ ( C + )ǫ ( + π C + π 4 ( ) ln + 1 ( 4ǫ 1 +ǫ 1 ǫ +ǫ ) ǫ ǫ 3 ln ( ) ǫ ǫ 3 + [ ǫ 1 ǫ +(C +1)ǫ 1 +(C 1)ǫ 1ǫ +Cǫ Cǫ ǫ 3 ] )}. (10) In these expressions, ǫ nd ǫ 3 re slow-roll prmeters belonging to the Hubble flow hierrchy [7, 8] ǫ n+1 dlnǫ n /(dn) (N being the number of e-folds) nd C is numericl constnt C The quntity is the so-clled pivot scle the loction of which is rbitrry but is usully chosen in the middle of the rnge of scles probed by the CMBR. These pproximtions fil to reproduce the power spectrum on smll scles becuse the lever rm is then too importnt for the ccurcy of the slow-roll clcultion. This is rther unusul, s the slow-roll method generlly provides n excellent pproximtion. This bre-down of the slow-roll clcultion on smll scles is clerly seen in Fig., which compres the full numericl clcultion of the power spectrum with the bove pproximtions for lrge field model qudrtic infltionry model. As one might expect, second-order slow-roll clcultions re more ccurte on smll scles thn the first order (which, mong others justifies the techniques tht hve been developed for this purpose [9, 30]). 1. Evolution of Density Perturbtions fter Infltion We re now interested in the behvior of perturbtions fter the end of infltion. On time scles t H 1 end (the subscript end denotes herefter quntities evluted t the end of infltion), the inflton vcuum expecttion vlue evolves pproximtely s ( end ) 3/sin(mt) φ(t) φ end, (11) for qudrtic inflton potentil. Let us recll tht φ end m Pl /( π) nd H end m/ 3 in the prticulr cse of lrge field qudrtic infltionry model. Plugging Eq. (11) into Eq. (5), it is esy to see tht the dominnt term of the lst four grvittionl terms is the lst one s it scles s 1/ 3/ while the others decy s 1/ 3. Neglecting the other three terms one finds n eqution FIG. : Primordil power spectrum of curvture perturbtions (solid line) for the lrge field model V(φ) = m φ /. Here N exit() N end is ln[ cross()/ end ], with cross() the scle fctor t Hubble rdius crossing of mode during infltion (i.e. / cross πh), nd end the scle fctor t the end of infltion. The dotted nd dshed lines respectively represents the first nd second order slow-roll results s given by Eqs. (9) nd (10) with ǫ 1 1/[(N + 1/)], ǫ = ǫ 3 = 1/(N + 1/) nd N 50. The slow-roll result is ccurte in the CMBR region but bres down on smller scles where the predicted power spectr become negtive. The upturn is relted to modes tht hve not crossed out the Hubble rdius nd hve remined in their vcuum stte, see text for detils. of motion of the form d ṽ ( [1+ dz + m end ) ] 3/cos(z) 6κφ end ṽ = 0, (1) where we hve defined z mt + π/4. This eqution is similr to Mthieu eqution d ṽ dz +[ A qcos(z) ] ṽ = 0 (13)

5 the resonnce bnd during infltion (blue lines in Fig. 3) nd those tht hve entered it during the inflton oscilltion er (green lines in Fig. 3). Of course, the preheting instbility is only opertive fter infltion nd before reheting. Hving identified the scles ffected by the resonnce, we now discuss their mplifiction. For the first instbility bnd, the Floquet index is given by µ = q/. Then the mode function evolves s ṽ e µz. However, in the present cse, we hve time-dependent Floquet index, for which the corresponding solution cn be written s [13] ( ṽ exp 5 ) µdz 3/. (17) FIG. 3: Hubble rdius l H (red solid line), lower bound of resonnce bnd l C (red dshed line) nd physicl wvelengths (blue nd green dotted lines), in Plnc units, versus the e- fold number during the finl stge of infltion nd the inflton oscilltion er for lrge field m φ / model. The reheting temperture is T rh = 10 7 GeV hence the number of e-folds during the inflton oscilltion er is N rh 5.6. The inset shows the detiled behviors of l H nd l C t the trnsition between infltion nd inflton oscilltion. with A = 1+ m, (14) 6κ ( q = φ end ) 3/ end. (15) Since q 1, we re in the nrrow resonnce regime. The first instbility bnd is given by 1 q < A < 1+q which mounts to 0 < < 3Hm. (16) From the bove inequlity, we note the emergence of newchrcteristicsptillengthgivenbyl C 1/ 3Hm, which sets the lower bound to the rnge of resonnt sptil scles. During infltion, in lrge field model, l C is lmost constnt nd l C l H, with l H H 1. After the end of infltion, during the inflton oscilltion er, with H 3/ one hs l C 3/4 nd, therefore, l C grows more slowly thn the Hubble rdius l H 3/. The sitution is summrized in Fig. 3 where both scles re displyed (solid nd dshed red lines). Their evolution hs been computed by n exct numericl integrtion of the equtions of motion in the cse of lrge field model. Since the physicl wvelength of comoving Fouriermode, λ, growsfsterthnl C, the sclesmy cross l C during infltion or during the mtter-dominted er nd thereby enter the instbility bnd. Furthermore, mode tht hs entered this bnd will lwys remin in it, t lest until reheting is completed. One cn thus distinguish two types of modes: those tht hve entered Note tht, ting the integrl of the time-dependent µ does not correspond to the WKB pproximtion, but rther is n educted guess confirmed by numericl simultion, see below. This implies tht v nd therefore ζ is constnt ccording to its definition Eq. (). Note lso tht the curvture perturbtion usully remins constnt only on super-hubble scles while decying on sub- Hubble scles, unless there is grvittionl instbility. We re now interested in the behvior of the frctionl mss-energydensityperturbtion, δ (η) δρ /ρ, where ρisthebcgroundenergydensityofthesclrfield. The perturbed Einstein eqution implies tht δ = ( ) 3 H +3 Φ Φ H, (18) where Φ is the Brdeen potentil [18, 31]. The Brdeen potentil is relted to the quntity ζ by ζ = H 1 Φ +Φ +Φ, (19) 3 1+p/ρ where p denotes the pressure. For constnt p/ρ, if Φ is constnt so is ζ. This is generic result on lrge scles in one fluid models. If Φ evolves s power-lw in the scle fctor (lso generic result), Φ s, then H 1 Φ = sφ nd consequently, Eq. (19) indictes tht ζ follows the sme power-lw behvior. Therefore, the constncy of ζ in the resonnce bnd implies tht Φ lso remins constnt. This is confirmed by numericl integrtion of the equtions of motion, see Fig. 4. Then ζ = 5Φ /3, s the pressure vnishes on verge during the preheting stge. Consequently, one my express δ s function of ζ s follows δ = ( ) 5 H +3 ζ. (0) This reltion remins vlid in both the lrge scle nd the smll scle limit (with respect to the Hubble rdius) provided the modes re in the instbility bnd (i.e. ζ is constnt). Given tht H 1/ in inflton oscilltion dominted Universe, one then finds tht δ grows s

6 which enters the resonnce bnd t the time of completion of reheting. It corresponds to the green dotted line lbelled mx in Fig. 3. Interestingly enough, this scle hs never been super-hubble. From the condition mx / rh = 3H rh m, where H rh is the Hubble prmeter t the end of the reheting epoch or t the beginning of the rdition dominted phse, one obtins ( ) 1/ ˆ mx = m m Pl ( ) 1/3 ( ) 1/3 Trh Hend 10 7 GeV () GeV 6 FIG. 4: Evolution of the density perturbtions δρ /ρ, ζ nd Φ towrds the end of infltion nd during the ensuing period of metric preheting fter infltion. The nottion is s in the text with phys /; N denotes the e-fold number of the scle fctor. The reduced Plnc mss M Pl m Pl / 8π is used in this plot. (t) on sub-hubble scles but remins constnt on super- Hubble scles. This behvior is lso confirmed by the numericl integrtion whose results re shown in Fig. 4. The preheting instbility for ll sub-hubble modes with / < 3Hm my thus be understood s the grvittionl instbility of the oscillting sclr condenste. Let us summrize our min findings. We hve found tht ll modes in the instbility bnd nd outside the Hubble rdius re such tht ζ nd δ re constnt. This is of course the cse for modes on scles probed by the CMBR nd, therefore, one recovers the fct tht the infltionry power spectrum propgtes through the reheting stge without being ffected. However, for those modes which re in the instbility bnd nd inside the Hubble rdius, one finds tht ζ is constnt but tht δ grows s (t). Therefore, nd contrry to the stndrd lore, the preheting stge cn ffect the behvior of the perturbtions even if there is only one non-intercting fluid.. Spectrum of Density Perturbtions t Reheting From the bove considertions, we see tht the mplifiction described bove occurs only for rnge of scles, [ min, mx ] in terms of comoving wvenumber. We nowdiscussthispointinmoredetilndspecify min nd mx. Henceforth, we introduce the short-hnd nottion: ˆ end H end, (1) where the lbel end refers s before to the end of infltion. The smllest unstble sptil scle is tht The lst eqution ssumes g = 30 degrees of freedom t reheting. All modes with ˆ < ˆ mx hve entered the resonnce bnd during infltion or during the preheting er. We hve lso seen tht necessry condition for the growth of δ on given scle is tht the corresponding mode lies within the Hubble rdius t some time before the reheting epoch. As consequence, the lrgest sptil scle tht is unstble with respect to the growth of δ is tht which re-enters the Hubble rdius t the end of the pre-reheting epoch, i.e. min rh H rh. This corresponds to the blue dotted line lbeled min in Fig. 3 nd ( ) /3 ( ) 1/3 ˆ min = Trh Hend 10 7 GeV (3). GeV As the density contrst increses in direct proportion to the scle fctor for modes in the rnge [ min, mx ], non-liner structures my eventully form. The mplitude of the density contrst t the end of the preheting epoch is given by the density contrst t the end of infltion nd by the e-fold number of growth while the mode is sub-horizon sized nd in the resonnce bnd. At this point, one should distinguish the modes tht enter the resonnce bnd during infltion (the dotted blue lines in Fig. 3) nd those tht hve entered the bnd during the preheting er but tht hve never exited the Hubble rdius (the dotted green lines in Fig. 3). We first ddress the former ones, for which [ min, end ], with ˆ end 1. Given tht, in the resonnce bnd, the Fourier density contrst remins constnt on super-hubble scles, the behvior of δ on sub-hubble scle cn be pproximted s δ = δ hc = δ end, (4) hc hc where the lbel hc refers to the epoch of Hubble rdius crossingof mode nd δ end is the density contrst clculted t the end of infltion, see Fig.. Let us precise the vlue of hc : in order to mtch both symptotic behviorsofδ in thesuper-horizon,nmelyδ 6ζ /5ndin the sub-horizon regimes, nmely δ [/(H)] ζ /5, one defines hc s tht t which = 3 hc H hc. This

7 7 implies tht hc = 3 endˆ. Now, using Eq. (0), one cn write δ end = (6/5) π P ζ ()/ 3. Then, the power spectrum of δ t the time of reheting is evluted s ( Pδ rh () = 3 δ end rh π hc ( = ˆ4 Pζ ) ( Trh ) 8/ GeV ( ) 4/3 Hend (5) GeV The vrince of the density fluctution enclosing the scle defined by is directly expressed t reheting s δ = 0 ) dln P rh δ ( ), (6) nd the condition δ 1 thus roughly mounts to Pδ rh 1. This singles out specific sptil scle (corresponding to wvenumber nl ), such tht scles with < nl hve not hd sufficient time to become nonliner, while those with > nl hve become non-liner by reheting. Strightforwrd mnipultions led to ) /3 ( ) 1/4 ( ˆ nl Pζ Trh GeV ( ) 1/3 Hend (7) GeV Let us now ddress the unstble modes tht hve lwys been sub-hubble, nmely [ end, mx ]. For these modes, Muhnov s vrible is still given by the vcuum solution, Eq. (4), until it enters the instbility bnd. This leds to ζ = e i(η η end) ǫ1 ( end ) ζ end, (8) where ǫ 1 3/ in the preheting er. In order to derive this expression, we hve set ǫ 1 = 1 t the end of infltion nd we hve ssumed tht Muhnov s vrible hs remined constnt through the (instntneous) trnsition from infltion to the preheting er; this is expected insofr s the modes re well sub-hubble t tht time, H, hence Muhnov s vrible oscilltes freely. The quntity ζ end represents the vlue of ζ t the end of infltion or t the onset of the oscilltory phse. Using Eq. (19) with p = 0, one finds ζ = Φ since ζ nd Φ must shre the sme power-lw behvior 1 (s discussed bove), hence H 1 Φ = Φ. Then, Eq. (18) leds to δ = 3 H Φ = ˆ e i(η η end) ζ end 3 ǫ1 = C()e iη. (9) The density contrst then oscilltes with constnt mplitude until the mode enters the resonnce bnd, in which δ = C()/ rc where rc denotes the scle fctor t the entry into the resonnce bnd. This ltter is given by rc = 4 /(9m 3 end H end ) which implies tht rh / rc = (ˆ mx /ˆ) 4. The power spectrum of density fluctutions t the time of preheting cn now be computed s Pδ rh 4 8ˆ () = 7 (ˆmx ˆ ) 8 P end ζ (). (30) The quntity Pζ end is the power spectrum t the end of infltion for the modes tht hve never exited the Hubble rdius, Pζ end () = ˆ ( ) Hend. (31) π m Pl We cn chec tht it givesthe correctorder ofmgnitude bycompringwiththeexctresultinfig.forthecseof lrge field model. Inserting this expression in Eq. (30), we see tht the spectrum now decreses s ˆ while it wsproportionlto ˆ 4 in theotherrngeofwvenumbers [ min, end ] tht we discussed previously. Finlly, for those modes tht hve never entered the resonnce bnd nd never crossed out the Hubble rdius, tht is to sy ˆ > ˆ mx, the density contrst lwys behves ccording to Eq. (9) nd, therefore, the corresponding power spectrum is given by Pδ rh 4 8ˆ () = 7 P ζ(). (3) This spectrum is proportionl to 6. Of course, these modes re still in vcuum stte nd therefore, despite the fct tht the vcuum fluctutions re lrge (which is lielyto be lined tothe pthologiclbehviorofthe theory in the ultr-violet regime), there is no energy density ssocited to this phenomenon. One wy to see this is to compute the evolution of the modes, nd in prticulr the Bogoliubov coefficient β of the negtive frequency WKB brnch t lte times, on sub-horizon scles (ssuming the vcuum corresponds to the positive frequency brnch). For these vcuum modes, β = 0 becuse the WKB pproximtion hs lwys remined vlid. The energy density stored in these modes, which is given by dδρ dlnˆ 4 π 4 β, (33) consequently vnish. The sitution is summrized in Fig. where the spectrum P δ () is displyed in the different regimes studied before. The striing feture of this plot is of course tht, for > nl, the power spectrum is non-liner t the time of reheting; let us stress tht of course, the plotted vlues only represent the extrpoltion of the liner clcultion.

8 The miniml mss, ssocited to mx is given by ( ) M ( mx ) Trh 10 7 GeV ( ) 3/ m g. (36) m Pl The mximum mss is ssocited to the scle nl, i.e. 8 FIG. 5: Power Spectrum P rh δ (), computed t reheting s givenbyeqs.(5),(9)nd(3). Inthernge [ min, end ], the spectrum is evluted numericlly, see Fig.. The smll discontinuity observed t end is just numericl rtefct, the ctul power spectrum being of course continuous. The non-linerity of modes with > end is rther intriguing, s these modes hve never crossed out the Hubble rdius, hence one might question their clssicl nture. However, it is rther strightforwrd to show tht forthese modes, β 1, hence energydensityisindeed stored on those scles. From this point of view, the non-linerities observed in the power spectrum should priori be ten seriously nd viewed s sign tht structure formtion hs strted on those scles. In order to ssess the nture of these modes, one should pply the discussion of the so-clled quntum to clssicl trnsition of modes [1, 3] to the present cse. However, this question is beyond the scope of this pper. 3. Associted Mss Scles Finlly, it is interesting to quntify the rnge of msses ssocited with the collpsing structures. In terms of the wvenumber of the perturbtion, the enclosed mss cn be written s M = 4π ρ 3 (π)3 end ˆ 3 Hend 3. (34) According to the previous considertions, the mss corresponding to mximum growth of the density fluctutions corresponds to the scle end. Using the bove expression, strightforwrd mnipultions led to the following expression ( ) 1 M end Hend g. (35) GeV ( ) 3/4 ( ) M nl Pζ Trh g. (37) GeV Quite interestingly, the bundnce of primordil blc holes with mss 10 9 g is constrined by erly Universe cosmology. Such blc holes indeed evporte before BBN hs strted nd the blc hole temperture is high enough to llow the production of copious mounts of lightest supersymmetric prticles, in excess of the closure density [33], s well s grvitinosnd other lte time decying mssive prticles, whose decy products would ruin the success of BBN [34, 35]. If such blc holes leve Plnc-mss relics t their finl point of evportion, dditionl constrints cn be obtined from closure density rguments [36]. For blc hole msses 10 9 g, stringent constrints cn be derived from the impct of the evported rdition on BBN, see [37] for detiled discussion of this nd other observtionl constrints. B. Tensor modes In the previous section, we hve studied the mplifiction of density perturbtions during the reheting epoch. Nturlly, one my wonder whether tensor modes cn be mplified through similr instbility. Their Fourier mplitude h = µ / obeys the eqution [38] µ + ( ) µ = 0. (38) The effective potentil is thus much simpler thn in the cse of density perturbtions s it only involves the scle fctor nd its time derivtives up to second order. Defining µ 1/ µ nd woringin termsofthe cosmictime, one obtins ( ) µ + 3Ḣ 9 4 H µ = 0. (39) Using the explicit form of the inflton field given by Eq. (11), the eqution of motion tes the form d µ [ dz + m + 3κ ( end ) ] 3cos(z) 8 φ end µ = 0, (40) where z mt. Therefore, one lso obtins Mthieu eqution with A = m, (41) q = 3κ ( end ) 3 16 φ end. (4)

9 9 These equtions should be compred to Eqs. (14) nd (15). We re still in the regime q 1 but the crucil point is tht there is no term 1 in the expression of A in Eq. (41). This mens tht, for given scle, the condition 1 q < A < 1+q cnnot be mintined since A does not scle s q does. Therefore the mplifiction does not occur for grvittionl wves. IV. THE QUARTIC CASE A. A qurtic potentil V(φ) = λφ 4 /4 In the previous section, we hve studied lrge field model with qudrtic potentil. This corresponds to pretty generic sitution since, close to its minimum, the inflton potentil cn lwys be pproximted by prbol. However, it might lso be tht, for some unspecified reson (genericlly, the presence of prticulr symmetry), the qudrtic term in the Tylor expnsion (6) is bsent. Then, one hs to repet the previous nlysis in the cse where the potentil is given by V(φ) = λ 4 φ4. (43) In the following, we find tht growth lso occurs for this cse. This result is prticulrly interesting, s it shows tht the growth of density perturbtions goes beyond the simple nlogue of Newtonin collpse in dust Universe. Indeed, in the cse of λφ 4 /4, the effective eqution of stte is tht of rdition. As before, we ssume tht the potentil is pproximted by the bove qurtic form in the preheting er without ming ny specific ssumption on the ctul potentil during the infltionry er. We note tht qurtic infltion cn be normlized to COBE nisotropies with λ 10 1, but tht t the sme time such n infltionry potentil is disfvored by detiled nlysis of CMB dt [4]. For qurtic potentil, the evolution of the inflton field round its minimum cn be expressed in the lrge time limit t H 1 end s [4, 5, 39, 40] ( end ) ( ) 1 φ(η) = φ end cn x x end ;, (44) where x end φ end λη, η being the conforml time nd cn denotes the Jcobi cosine function [41, 4]. This solutionshouldbeinsertedineq.(5). Then, eepingonlythe dominnt terms, the eqution tht controls the evolution of Muhnov s vrible reduces to [ ṽ ( end ) + +3λφ end cn (x x end ; 1 ) ] ṽ = 0. (45) It is more convenient to chnge bc to the vrible v nd to express this eqution in terms of the vrible x defined bove: d v dx + [ φ end end λ +3cn ( x x end ; 1 )]v = 0. (46) As emphsized in Refs. [5, 39, 40], this eqution is the so-clled Lmé eqution for which the Floquet chrt is exctly nown. There is mplifiction if 3/ < /(φ end end λ) < 3. Using Hend κλφ4 end /1, the bove inequlity cn be cst into the following form ( 3 m Pl < π φ ˆ 3 ) 1/ 3 mpl <. (47) end π φ end For the prticulr cse of qurtic lrge field infltion, φ end /m Pl 1/ π, hence the lower nd upper bound of the bove bnd reduce respectively to 3/ nd (3 3/) 1/. For the bove instbility, the mximum growth occurs with Floquet index µ t /(φ end end λ) (i.e. ˆ 1.56 for lrge field qurtic infltionry model) [5]. The Floquet index is constnt, therefore growth becomes exponentil, unlie wht hs been found for the qudrtic cse. Severl remrs re in order t this point. Firstly, the resonnce bnd corresponds to fixed nrrow rnge of comoving scles, contrry to the qudrtic cse for which the bnd widens with time. Secondly, the modes tht re subject to prmetric resonnce re of sptil scle comprble with the Hubble rdius t the end of infltion. Thirdly, the effect found here is similr to tht investigted in Ref. [4, 5] except tht it concerns Muhnov s vrible rther thn the perturbed inflton field or the Fourier mplitude of sclr field χ coupled to the inflton s g φ χ /, where g is the coupling constnt. The cse of density perturbtions cn be viewed s specil cse with specific vlue of g. But, of course, the underlying physics is different in the sense tht metric fluctutions re different from n externl sclr field coupled to φ. In order to chec the effect discussed bove, we hve numericlly integrted the eqution of motion for lrge field infltionry qurtic potentil, see Fig. 6. The numericl nlysis confirms the previous nlyticl nlysis, lthough we find tht the instbility bnd is shifted by 1% to higher vlues in comprison with Eq. (47); this effect is liely due to our ssumption of instntneous trnsition from infltion to (rdition dominted) preheting er in the nlyticl clcultions. In detil, one finds the bnd t < ˆ < When the mode is outside the resonncebnd (left pnel in Fig. 6), the perturbtions remin constnt while when it is inside (right pnel in Fig. 6), we observe violent mplifiction of Φ, ζ nd δ. To summrize, the qurtic cse presents substntil differences in comprison with the qudrtic cse. There

10 10 FIG. 6: Evolution of the density perturbtions δρ /ρ, ζ nd Φ towrds the end of qurtic infltion nd during the ensuing period of metric preheting fter infltion. The dimensionless prmeter is λ The red solid line represents the Hubble rdius while the blue dotted line represents the physicl wvelength of the Fourier mode. Left pnel: the Fourier mode is chrcterized by ˆ 1.47 nd lies outside the resonnce bnd. Right pnel: ˆ nd the mode is inside the instbility bnd s the exponentil growth indictes. exists nrrow instbility bnd t scles of order of the Hubble rdius t the end of infltion, of fixed width (in terms of comoving wvenumber) nd with ssocited exponentil growth. B. The mixed qudrtic-qurtic cse The purely qurtic cse rises only if the qudrtic term is bsent in the potentil. Here we briefly ddress the more generl cse, in which both qudrtic nd qurtic terms re present: V(φ) = 1 m φ + λ 4 φ4. (48) At lrge vcuum expecttion vlues of the inflton, the qurtic term lwys domintes while t smll vlues of φ, the qudrtic term is the most importnt one. For which vlue of φ these terms re of equl importnce depends on the two free prmeters m nd λ. If the potentil were purely qurtic, then infltion would stops t φ end = m Pl / π. For this vlue, the two terms re equl if m = λ/(π)mpl. It is therefore convenient to prmeterize the potentil (48) by λ nd the dimensionless prmeter Υ defined by Υ m m Pl π λ. (49) The regime where Υ 1 reduces to the qudrtic cse while the cse Υ 1 is equivlent to the qurtic cse. The most interesting sitution is of course when Υ is of order one. In order to study this cse, we hve numericlly integrted the equtions of motion, ssuming tht the potentil form remins the sme during infltion nd the prehetinger. Asfirstexmple, wehvetenλ = 10 1, tht is to sy the sme vlue s in the previous section for numericl pplictions nd Υ = We hve lso chosen ˆ = 1.738, i.e. vlue which is inside the resonnce bnd in the pure qurtic cse, see right pnel in Fig. 6. The result is displyed in Fig. 7 (left pnel). Severl comments re in order. Firstly, one notices tht the Hubble rdius nd the chrcteristic scle l C no longer coincide. This is of course due to the fct tht the mss m is no longer COBE normlized. In fct, in this model, the end of infltion is delyed becuse t φ = m Pl / π, the qudrtic term still plys role. Regrding the perturbtions themselves, despite the fct tht the Fourier mode is inside the qurtic resonnce bnd of Eq. (47), there is no mplifiction. It is of course due to the fct tht, for smll vcuum expecttion vlues of the field, the dominnt term is the qudrtic one. We lso notice tht δ does not grow s either. This is becuse, since the end of infltion hs been postponed nd the vlue of l C incresed(in comprisonwith the correspondingvlue in the pure qudrtic cse), the time t which the mode crosses l C is lso delyed (nd therefore lies outside of the plot). It is cler tht if we decrese the vlue of Υ, nd eep the sme vlue for, we should observe the growth described in the lst section in the pure qurtic cse. This is indeed the cse s demonstrted in Fig. 7 (right pnel) where Υ = Just fter the end of infltion, the Fourier mode is inside the pure qurtic resonnce bnd nd the perturbtions grow. But s the mplitude of the inflton field decreses, the influence of the qudrtic term becomes more nd more importnt. As consequence, the exponentil growth stops t some point. Then, the mplitude of the fluctutions remin constnt until the mode exceeds the chrcteristic scle l C ssocited to the qudrtic prt of the potentil. As before, this time is delyed due to the smll vlue of the mss m

11 11 FIG. 7: Left pnel: evolution of the density perturbtions δρ /ρ, ζ nd Φ towrds the end of infltion with potentil given by Eq. (48) nd during the ensuing period of of metric preheting. The dimensionless prmeter re λ 10 1 nd Υ = The red solid line represents the Hubble rdius, the red dshed line the scle l C while the blue dotted line represents the physicl wvelength of the Fourier mode. The Fourier mode is chrcterized by ˆ Right pnel: sme s left pnel but with Υ = nd ˆ = The vlue of is the sme for the two pnels but the vlue of ˆ sightly differs becuse the end of infltion does not occur exctly t the sme time in the two cses. (smll in comprison with the COBE normlized vlue) nd lies outside of the plot. Then, the density contrst grows proportionl to the scle fctor s discussed in Section III. V. CONCLUSION In this pper we hve investigted the evolution of subhorizon perturbtions between the end of cosmic infltion nd the beginning of rdition dominted erly Universe. This study ssumes tht energy density of the Universe during this phse is dominted by the oscilltions of single sclr field round its minimum in qudrtic or qurtic potentil. We hve found tht perturbtions re subject to preheting instbility in the nrrow resonnce regime, leding to the growth of the density contrst δρ /ρ on sub-horizon scles. In the cse of qudrtic potentil, the domin of instbility (in terms of comoving wvenumbers) widens with time. This implies tht wvenumber tht hs entered the bnd t some time remins inside for the whole durtion of the preheting er. The density contrst of sub-horizon unstble modes then increses s the scle fctor, reching nonlinerity for wide rnge of scles. Mximum growth occurs for sptil scles of the order of the Hubble rdius t the end ofinfltion, to which is ssocitedmss scle M end g ( H end /10 13 GeV ) 1. A qurtic potentil leds to different phenomenology, lthough the min effect of preheting instbility remins. Notbly, the rnge of unstble wvenumbers is now fixed in terms of comoving wvenumbers nd peed round the Hubble scle t the end of infltion. On the other hnd, the density contrst of sub-horizon unstble modes increses exponentilly fst. The growth of structure on sub-horizon scles is nown to occur in mtter dominted er, hence the growth of the density contrst for qudrtic inflton potentil cn thus be well understood in this respect. However, the result obtined for qurtic potentil confirms tht the effect goes beyond simple nlogy with Newtonin collpse of self-grvitting dust cloud, s the effective eqution of stte in this ltter cse is tht of rdition. These results open interesting venues of reserch in the serch for observtionl signtures of the dynmics of the very erly Universe. In prticulr, s we hve noted, the bundnce of blc holes of mss comprble to the bove is constrined by cosmologicl rguments. The formtion of inflton hlos my lso led to copious production of grvittionl wves, s discussed in compnion pper. Finlly, the very existence of nonlinerities on smll scles csts some shdow on the relibility of the predictions for curvture perturbtions on those lrge scles responsible for the present lrge-scle structure. These ltter could be modified if effects of mode-mode coupling become importnt. [1] A. D. Linde, Phys. Lett. B19, 177 (1983). [] A. D. Linde, Lect. Notes Phys. 738, 1 (008), [3] E. Komtsu et l., (010), [4] L. Kofmn, A. D. Linde, nd A. A. Strobinsy, Phys.

12 1 Rev. D56, 358 (1997), hep-ph/ [5] B. A. Bssett, S. Tsujiw, nd D. Wnds, Rev. Mod. Phys. 78, 537 (006), stro-ph/ [6] J. R. Ellis, A. D. Linde, nd D. V. Nnopoulos, Phys. Lett. B118, 59 (198). [7] D. V. Nnopoulos, K. A. Olive, nd M. Srednici, Phys. Lett. B17, 30 (1983). [8] K. Jedmzi, Phys. Rev. D74, (006), hepph/ [9] M. Kwsi, K. Kohri, T. Moroi, nd A. Yotsuyngi, Phys. Rev. D78, (008), [10] S. Billy, K.-Y. Choi, K. Jedmzi, nd L. Roszowsi, JHEP 05, 103 (009), [11] Y. Nmbu nd A. Truy, Prog. Theor. Phys. 97, 83 (1997), gr-qc/ [1] T. Hmzi nd H. Kodm, Prog. Theor. Phys. 96, 113 (1996), gr-qc/ [13] F. Finelli nd R. H. Brndenberger, Phys. Rev. Lett. 8, 136 (1999), hep-ph/ [14] M. Khlopov, B. A. Mlomed, nd I. B. Zeldovich, Mon. Not. Roy. Astron. Soc. 15, 575 (1985). [15] Y. Nmbu nd M. Ssi, Phys. Rev. D4, 3918 (1990). [16] B. Rtr, Phys. Rev. D44, 35 (1991). [17] J.-C. Hwng, Phys. Lett. B401, 41 (1997), stroph/ [18] V. F. Muhnov, H. A. Feldmn, nd R. H. Brndenberger, Phys. Rept. 15, 03 (199). [19] J. Mrtin, Brz. J. Phys. 34, 1307 (004), stroph/ [0] J. Mrtin, Lect. Notes Phys. 669, 199 (005), hepth/ [1] J. Mrtin, Lect. Notes Phys. 738, 193 (008), [] J. Mrtin nd D. J. Schwrz, Phys. Rev. D57, 330 (1998), gr-qc/ [3] H. Kodm nd T. Hmzi, Prog. Theor. Phys. 96, 949 (1996), gr-qc/ [4] J. Mrtin nd C. Ringevl, JCAP 0608, 009 (006), stro-ph/ [5] L. Lorenz, J. Mrtin, nd C. Ringevl, JCAP 0804, 001 (008), [6] L. Lorenz, J. Mrtin, nd C. Ringevl, Phys. Rev. D78, (008), [7] S. M. Lech, A. R. Liddle, J. Mrtin, nd D. J. Schwrz, Phys. Rev. D66, (00), stro-ph/ [8] D. J. Schwrz, C. A. Terrero-Esclnte, nd A. A. Grci, Phys. Lett. B517, 43 (001), stro-ph/ [9] J.-O. Gong nd E. D. Stewrt, Phys. Lett. B510, 1 (001), stro-ph/ [30] J. Mrtin nd D. J. Schwrz, Phys. Rev. D67, (003), stro-ph/ [31] J. M. Brdeen, Phys. Rev. D, 188 (1980). [3] J. Lesgourgues, D. Polrsi, nd A. A. Strobinsy, Nucl. Phys. B497, 479 (1997), gr-qc/ [33] A. M. Green, Phys. Rev. D60, (1999), stroph/ [34] M. Lemoine, Phys. Lett. B481, 333 (000), hepph/ [35] M. Y. Khlopov, A. Brru, nd J. Grin, Clss. Qunt. Grv. 3, 1875 (006), stro-ph/ [36] B. J. Crr, J. H. Gilbert, nd J. E. Lidsey, Phys. Rev. D50, 4853 (1994), stro-ph/ [37] B. Crr, K. Kohri, Y. Sendoud, nd J. Yooym, (009), [38] L. P. Grishchu, Sov. Phys. JETP 40, 409 (1975). [39] D. I. Kiser, Phys. Rev. D56, 706 (1997), hepph/ [40] D. I. Kiser, Phys. Rev. D57, 70 (1998), hepph/ [41] M. Abrmowitz nd I. A. Stegun, Hndboo of mthemticl functions with formuls, grphs, nd mthemticl tbles, ninth ed. (Ntionl Bureu of Stndrds, Wshington, US, 1970). [4] I. S. Grdshteyn nd I. M. Ryzhi, Tble of Integrls, Series, nd Products (Acdemic Press, New Yor nd London, 1965).