Warm Inflation with Nonminimal Derivative Coupling

We study the effects of the nonminimal derivative coupling on the dissipative dynamics of the warm inflation where the scalar field is nonminimally coupled to gravity via its kinetic term. We present a detailed calculation of the cosmological perturbations in this setup. We use the recent observational data from the joint data set of WMAP9 + BAO + H 0 and also the Planck satellite data to constrain our model parameters for natural and chaotic inflation potentials. We study also the levels of non-Gaussianity in this warm inflation model and we confront the result with recent observational data from the Planck satellite.


Introduction
Although the hot big bang scenario has considerable successes in explaining the general properties of our universe, this model cannot address properly some issues such as flatness, horizon, and relics problems.The latest observations show that a period of accelerating expansion (inflation) during the early history of the universe with ä > 0 ( < −/3) can solve these problems.Inflation also provides the primordial perturbations which seed the formation of structure in the universe [1].Despite the great successes of the inflation paradigm, there is a problem for realization of this scenario that we do not know how to integrate it with ideas of particle physics.Also problems such as unexpected low power spectrum at large scales and egregious running of the spectral index [2] are some yet unsolved problems in this paradigm.From a thermodynamic viewpoint, there are two different dynamical realizations of inflation which are called cold and warm inflation.In the standard inflationary model, called supercooled or isentropic inflation, the universe rapidly supercools during the inflation and its temperature decreases, which is a result of the exponential expansion.When the inflation ends, a reheating phase occurs, which explains the subsequent processes and evolution.The fluctuations created during inflation are zero-point ground state fluctuations and evolution of the inflaton field is governed by the ground state evolution equations.The other picture of inflation dynamics is the warm or nonisentropic inflation [3][4][5].In this picture, the dissipative effects are significant.So, that radiation production and inflationary expansion occur concurrently.Several mechanisms for implementing such a dissipation during inflation have been proposed [6,7].In warm inflationary models, the sources of density fluctuations are the thermal fluctuations.Matter fields interact with particles that are in a thermal bath with mean temperature smaller than the grand unified theories (GUT) critical temperature [8].The most important candidate for the inflaton is the Standard Model's Higgs field.Nevertheless, if Higgs field is minimally coupled to gravity, it cannot realize a good inflation stage [9].To preserve the Higgs field as an inflaton candidate, a model has been postulated in which the Higgs field is nonminimally coupled to gravity [10].It has been argued in [11,12] that this nonminimal coupling breaks the unitary bound by giving rise to nonrenormalizabe operators.However, the authors in [13] have shown that this is not actually the case and the unitary bound is preserved in this setup.On the other hand, if gravity is nonminimally coupled to derivatives of the scalar field, the unitary bound is preserved during inflation.So, this model can explain the inflationary phase properly [14,15].In recent years, the issue of inflation and also the late time accelerated expansion with nonminimal derivative coupling are studied in detail [16][17][18][19][20]. Nevertheless, the warm inflation dynamics in the context of nonminimal derivative coupling and perturbations in this setup have not been studied yet.Therefore, the present study is devoted to filling this gap.
With these preliminaries, the goal of the present study is to investigate the effects of the nonminimal derivative coupling on the dissipative dynamics of the warm inflation.We investigate the effects of the nonminimal derivative coupling and dissipation on the inflationary dynamics of the model.Then, we present a detailed calculation of the perturbations in this framework.By using the recent observational data from WMAP9 + BAO +  0 and also the Planck satellite data, we constrain our model parameters for natural and chaotic inflation potentials.Also we study the levels of non-Gaussianity in this scenario and we confront the obtained results with the recent observational data from the Planck satellite.

The Setup
The action of a 4-dimensional model in the presence of the nonminimal derivative coupling between the scalar field and gravity is given by where  ] =  ] − (1/2)  ] is the Einstein tensor and  is a constant having the dimension of mass and   = 2.4 × 10 18 GeV is the reduced Planck mass. is the Ricci scalar and () is the potential of the scalar field .Variation of action (1) with respect to the scalar field gives the following equation of motion: If we consider phenomenologically a dissipation coefficient Γ that is responsible for the decay of the scalar field into the radiation during the inflationary regime, the scalar field equation of motion will be as follows: Γ can be assumed to be a constant or a function of the scalar field  or the temperature  or both of them.In the slow-roll approximation where φ ≪ |(3 + Γ) φ | and | Ḣ| ≪  2 , the scalar field equation of motion takes the following form: Note that the sufficient condition for canonical normalization of the scalar field is given by  2 ≫  2 [14,15].The energy density and the pressure of the scalar field are given by (5) Also, the energy conservation equation of the model is given by During the inflaton evolution, the dissipation effect leads to the production of entropy.The entropy density of the radiation (, ) in this nonminimal setup is defined by  = − ∫(()/).Now the total energy density of the system is given by and using the relation (6), the rate of entropy production can be deduced.In this setup, entropy production is affected by the nonminimal derivative coupling between the scalar field and gravity.This can be seen more explicitly in the following relation: where by (4), φ is directly related to the nonminimal derivative coupling.The basic idea of the warm inflation is that radiation production is occurring concurrently with inflationary expansion due to dissipation from the inflaton field system.The equation of state for radiation is given by   =   /3.Therefore, the continuity equation for the radiation yields the following result: We assume a quasi-stable radiation production during the warm inflation phase [21,22]; that is, ρ  ≪ 4  and ρ  ≪ Γ φ 2 .So we can write Substituting (4) into this relation, we get Now we define the dissipation factor as follows: which is a dimensionless parameter.With this definition, we can write During the slow-roll stage, the scalar field evolution is damped.For the high (or weak) dissipation regime we have  ≫ 1 (or  ≪ 1), respectively.Now, in the spatially flat FRW background, the Friedmann equation can be written as follows: The slow-roll regime can be parametrized by a set of slowroll parameters , , and , defined by Inflation can be realized only if {, , } < 1; once one of these parameters reaches the unity, the inflation phase terminates.We need to calculate these parameters in our nonminimal derivative coupling model.For  we get and the second slow-roll parameter  becomes And finally, the third slow-roll parameter takes the following form: where, as usual, a prime marks the differentiation with respect to .The number of e-foldings during inflation is defined as which in the slow-roll approximation can be written as where   (  ) is the value of the scalar field at the end (beginning) of the inflationary phase.

Perturbations
In this section we study scalar and tensor perturbations in our model.The scalar perturbations can be decomposed to entropy or isocurvature perturbations and adiabatic or curvature perturbations.In the framework of warm inflation, the scalar and radiation fields are interacting.Therefore, isocurvature perturbations are generated besides the adiabatic ones.This occurs because warm inflation can be considered as an inflationary model with two basic fields.
Dissipative effects can produce a variety of spectral indexes, ranging between red and blue and thus producing the running blue to red spectral index suggested by WMAP data [23][24][25].In the longitudinal gauge, the scalar metric perturbations of the FRW background are given by [26,27] where () is the scale factor and Φ = Φ(, ) and Ψ = Ψ(, ) are the metric perturbations.The spatial dependence of perturbed quantities is in the form of plane waves  ⋅ , where  is the wave number.Therefore, the perturbation of the metric leads to the perturbation in the energy-momentum tensor through Einstein's field equations.The perturbed Einstein's field equations can be obtained as follows: where V appears from the decomposition of the velocity field as   = −(  /)V ⋅ ( = 1, 2, 3) [16] and we have omitted the subscript .We perturb the scalar field equation of motion and the continuity equation for the radiation field as We consider these equations in the slow-roll approximation and on large scales where  ≪  [28,29].In this situation, (26) become, respectively, and ( 24) takes the following form: where we have used the relation of the velocity field as V ≃ −(/4)(Φ +   /4  + (3Γ φ /4  )).
We can solve the above set of equations by substituting ( 29) into (27), where we get Now, by defining a function  as (30) can be rewritten as A solution of this equation is given as where  is an integration constant.So, from (31), we find For simplicity we define the following quantity: With this definition, we can rewrite (34) as and, therefore, the density perturbation is given by The fluctuations of the scalar field are generated by thermal interaction with the radiation field, instead of quantum fluctuations where the wave number   is defined by   = √Γ + 3 2 (1 + 3( 2 / 2 )) and corresponds to the freezeout scale at which dissipation damps out the thermally excited fluctuations.So, the density perturbation amplitude becomes The scalar spectral index is given by The interval in wave number is related to the number of efoldings by the relation So we obtain ) . ( The running of the spectral index which is defined as in our model takes the following form: Now we pay attention to the tensorial perturbations.As has been mentioned in [30], the generation of tensor perturbations during inflation period produces stimulated emission in the thermal background of gravitational waves.This process changes the power spectrum of the tensor modes by an extra, temperature-dependent factor given by coth(/2).Therefore, the spectrum of tensor perturbations is given by The ratio between the amplitudes of tensor and scalar perturbations (tensor-to-scalar ratio) is given by where   = (25/4) 2  .This is the complete treatment of perturbations in this setup.In what follows, we consider two specific inflaton potentials to study the inflationary aspects of this model and confrontation with recent observational data from Planck + WMAP9 + BAO +  0 [23][24][25].

Specific Examples
where , , and  are defined as follows: Note that if we set  = 1 for the end of the inflation, the dimensionless field at the end of the inflation attains the following form: Substituting (49) into (47), we find  2 as follows: Advances in High Energy Physics 7 Using ( 42) and ( 46), the scalar spectral index and the tensor-to-scalar ratio in the warm chaotic inflation are given by respectively, where  is given by the relation (50).The left panel of Figure 1 shows the behavior of the tensor-to-scalar ratio versus the scalar spectral index with  = 1000.For  = 50, the model with 4 <  < 53.5, for  = 60, the model with 3.94 <  < 52.73, and for  = 70, the model with 3.90 <  < 51.86 lie inside the 95% CL of the Planck + WMAP9 + BAO +  0 joint data.We have also plotted the evolution of the running of the scalar spectral index versus the scalar spectral index in the right panel of Figure 1.For this case, the running is negative and in some range of  compatible with recent observational data.From another perspective, we can fix  at, for instance,  = 1000 and let  vary in order to see the effect of nonminimal derivative coupling.In Figure 2, we survey the effects of the nonminimal derivative decoupling.The left panel of this figure shows the behavior of the tensorto-scalar ratio versus the scalar spectral index with  = 1000.
For  = 50, the model with 23 <  < 141.3, for  = 60, the model with 13.67 <  < 140.05, and for  = 70, the model with 16 <  < 127.83 are compatible with observational data.Also, evolution of the running of the scalar spectral index versus the scalar spectral index in the background of Planck + WMAP9 + BAO +  0 joint data is shown in the right panel of Figure 2.For this case, the running is negative in some range of  and lies inside the 95% CL of the Planck + WMAP9 + BAO +  0 .

Natural Inflation.
Natural inflation is characterized by the potential where Λ and  are constants having the dimension of mass.
The number of e-foldings at the end of inflation with the above potential is given by where So, from ( 16) the value of  at the end of inflation takes the following form: Advances in High Energy Physics Substituting (55) into (53),  is determined fully for given values of , , and .Then   and  can be obtained, respectively, as follows: ) . (56) Now we consider the observational constraints on natural inflation parameters in our setup.As in the previous example, we consider three values of the number of e-folds.In Figure 3(a), we have shown evolution of the tensor-to-scalar ratio versus the scalar spectral index with  = 100 and variable .For  = 50, the model with  > 176.55, for  = 60, the model with  > 177.122, and for  = 70, the model with  > 177.05 are compatible with Planck + WMAP9 + BAO +  0 data.The right panel of Figure 3 shows the behavior of the running of the spectral index versus the scalar spectral index in the background of the Planck + WMAP9 + BAO +  0 joint data.The Figure has been plotted with  = 100.In our setup with natural inflation potential, the running of the scalar spectral index is negative and lies inside of the 95% CL of the Planck + WMAP9 + BAO +  0 joint data.Now we consider the opposite case; we fix  and let  to be variable.We fix  at  = 100.The left panel of Figure 4 shows the behavior of the tensor-to-scalar ratio versus the scalar spectral index with  = 100 and variable .For  = 50, the model with 156 <  < 273.6, for  = 60, the model with 156.03 <  < 273.04, and for  = 70, the model with 156.08 <  < 273.31 are compatible with data.The right panel of Figure 4 shows evolution of the running of the spectral index versus the scalar spectral index.For  = 50, the model with 112.37 <  < 211.71, for  = 60, the model with 113.5 <  < 213.38, and for  = 70, the model with 114.24 <  < 216.88 lie inside the 95% CL of the Planck + WMAP9 + BAO +  0 joint data.
In Tables 1 and 2, we have summarized the results of numerical analysis for different values of  for the above mentioned potentials in our warm inflation model with nonminimal derivative coupling.We have specified the range of , , and  parameters in which the model is compatible with Planck + WMAP9 + BAO +  0 joint data.

Non-Gaussianity
In this section we study non-Gaussianity of the primordial perturbations in our model.We calculate the three-point correlation function (bispectrum) of the density perturbation distribution in the Fourier space.Inflaton field equation of motion is given by the following expression: In order to study the fluctuations of the inflaton field (, ), it is assumed that the fluctuations are small and the full inflaton field is emerged as (, ) =  0 () + (, ), where  0 is homogeneous background field and (, ) ≪  0 ().Therefore, the full inflaton equation of motion with fluctuation in our setup will be as follows: where (, ) is a Gaussian noise term.The fluctuationdissipation theorem determines the properties of the noise, as shown in [31,32].So, this term in the momentum space has the following properties: In order to calculate the bispectrum, we consider the inflaton field fluctuation up to the second order (, ) =  1 (, ) +  2 (, ), where  1 is the first order of  and  2 is the second order.Thus, the equations of motion of the first Advances in High Energy Physics Therefore, the solution of the above equations will be as follows: where  and  are defined as follows: The second term on the right in both parts of (61) is memory terms.These memory terms within the Hubble time give the concept of freeze-out [33].By definition of freeze-out, for  ≥   this term damps away and for  ≤   the memory term reminds within a Hubble time.By definition,   can be obtained via the following condition: For warm inflation   ( 0 ) < Γ, so the freeze-out condition in this model can be simplified to Note that the second term in square root comes from the nonminimal derivative coupling.Now, we obtain the bispectrum before the end of inflation.The three-point correlation function (when the number of e-folding is equal to 60) becomes where where   represents the time at Hubble crossing of the smallest of the three inflation perturbation modes and   represents the time when the last of the three wave vectors thermalizes.The bispectrum for slow-roll, single field, supercooled inflation models, and also for the set of warm inflation models is given by [31,32] ⟨Φ ( 1 ) Φ ( 2 ) Φ ( 3 )⟩ =  inf (2) )) . (70) The  inf parameter for the pure warm inflation (without the nonminimal derivative coupling) is obtained as [19]  inf = − We see that the non-Gaussianity in the high-dissipation regime of warm inflation with nonminimal derivative coupling (70) is comparable to the high-dissipation regime of pure warm inflation (71) and the weak-dissipation regime of pure warm inflation which is given by

Conclusion
In this paper we have studied the cosmological inflation on the nonminimal coupling where the scalar field is nonminimally coupled to gravity via its kinetic term.We have studied the warm inflationary dynamics and dissipative effects of the model.We have calculated the inflation parameters and perturbations.Next we have surveyed observational constraints on two inflationary models with a field derivative coupling to the Einstein tensor.We also have explored the evolution of the inflationary parameters in the background of the Planck + WMAP9 + BAO +  0 data.For chaotic potential, as we see in Figure 1, the values of   , , and

4. 1 .
Chaotic Inflation.With () = (1/2) 2  2 , the number of e-foldings in the slow-roll approximation at the end of the inflation is given by

Figure 1 :Figure 2 :
Figure 1: Evolution of the tensor-to-scalar ratio (a) and running of the scalar spectral index (b) versus the scalar spectral index for chaotic potential in the background of Planck + WMAP9 + BAO +  0 data.The figure has been plotted for  = 50, 60, and 70 and with  = 1000 with varying .

Figure 3 :Figure 4 :
Figure 3: Evolution of the tensor-to-scalar ratio (a) and running of the scalar spectral index (b) versus the scalar spectral index for natural inflation in the background of Planck + WMAP9 + BAO +  0 data.The figure has been plotted for  = 50, 60, and 70, with  = 100 and  as a free parameter.

Table 1 :
Constraints on  =  2 / 2 for chaotic potential with  = 1000 and on  = ( 2 Λ 4 )/( 2  4  ) for natural potential with  = 100 for which the model of warm inflation with nonminimal derivative coupling is compatible with the Planck + WMAP9 + BAO +  0 joint data.

Table 2 :
Constraints on  with  = 1000 and  = 100 for which the model of warm inflation with nonminimal derivative coupling is compatible with the Planck + WMAP9 + BAO +  0 joint data.