Testing an Inflation Model with Nonminimal Derivative Coupling in the Light of PLANCK 2015 Data

We study the dynamics of a generalized inflationary model in which both the scalar field and its derivatives are coupled to the gravity. We consider a general form of the nonminimal derivative coupling in order to have a complete treatment of the model. By expanding the action up to the second order in perturbation, we study the spectrum of the primordial modes of the perturbations. Also, by expanding the action up to the third order and considering the three point correlation functions, the amplitude of the non-Gaussianity of the primordial perturbations is studied both in equilateral and orthogonal configurations. Finally, by adopting some sort of potentials, we compare the model at hand with the Planck 2015 released observational data and obtain some constraints on the model's parameters space. As an important result, we show that the nonminimal couplings help to make models of chaotic inflation, that would otherwise be in tension with Planck data, in better agreement with the data. This model is consistent with observation at weak coupling limit.


Introduction
The existence of an inflationary stage in the history of the early universe was proposed to solve some problems of the standard big bang cosmology, such as the flatness, horizon, and relics problems. A successful inflationary paradigm can also provide a causal mechanism for generating the initial density perturbations needed to seed the formation of structures in the universe [1][2][3][4][5][6][7][8]. A simple single field inflationary scenario (in which the early universe was dominated by the potential energy of the scalar field dubbed inflaton) predicts the dominant mode of the primordial density fluctuations to be highly adiabatic, scale invariant, and Gaussian [9]. But recent observational data have detected a level of scale dependence in the primordial perturbations [10,11]. On the other hand, although there is no direct signal for primordial non-Gaussianity in observation, Planck team has obtained some tight limits on primordial non-Gaussianity [12]. Some inflationary models also predict a level of non-Gaussianity in the perturbation's mode [9,[13][14][15][16][17][18][19][20]. This is really an important and interesting point, because a large amount of information on the cosmological dynamics deriving the initial inflationary expansion of the universe can be found in the primordial scalar non-Gaussianities [21,22]. In this regard, the extended inflationary models, which can show the non-Gaussianity and scale dependence of the primordial perturbation, are more favorable.
One successful class of the extended inflationary models is the one with a nonminimal coupling between the scalar field and the curvature scalar . This coupling usually is shown by the Lagrangian term as ( ) . There are some motivations behind including an explicit nonminimal coupling in the action. It is necessary for the renormalization term arising in quantum field theory in curved space [23]. Considering the quantum corrections to the scalar field theory, this term would arise at the quantum level [24,25]. With nonminimal coupling term, it is possible to have gravity as a spontaneous symmetry-breaking effect [26,27]. NMC term allows having an oscillating universe [28]. Also, with this term, it is possible to solve the graceful exit problem of the old inflation model by slowing the false-vacuum expansion [29,30]. So, it seems that the presence of nonminimal coupling 2 Advances in High Energy Physics between the scalar field and Ricci scalar is a useful extension of the theory and in this regard many authors have studied such models [17,[31][32][33][34][35][36][37][38][39][40].
It has been shown that the presence of NMC makes the inflation model more viable. For instance, the author of [41] has considered an inflation model with nonminimal coupling and noncanonical kinetic terms leading to a chaotic inflation model. It is shown that this chaotic inflation model is favored by the BICEP2 observation. But this model is in tension with recent observational data released by Planck 2015. The authors of [42] have considered an inflation model with a nonminimal coupling between the scalar field and gravity as ( 2 + 4 ) . They compared their model with Planck 2015 data in Einstein frame to explore the viability of the model. It is shown that this model just for = 65 and with some maximum values of is consistent with Planck 2015 data. In [43], a nonminimally coupled inflationary model with quadratic potential in Einstein frame has been considered. The authors of this paper have shown that their model for some values of the NMC parameter is consistent with observational data. For other literature dealing with these issues, we refer the reader to [44][45][46][47][48][49][50][51][52][53].
Another extension of the inflation theory arises from considering a nonminimal coupling between the derivatives of the scalar field and the curvature, which leads to the interesting cosmological behaviors [54]. This coupling is usually given by the Lagrangian term as Some authors have shown that, during inflation, with a nonminimal coupling between the scalar field and gravity, the unitary bound of the theory is violated [55][56][57]. However, it has been shown that the model with nonminimal couplings between the kinetic term of the inflaton (derivatives of the scalar field) and the Einstein tensor preserves the unitary bound during inflation [58]. Also, the presence of nonminimal derivative coupling is a powerful tool to increase the friction of an inflaton rolling down its own potential [58]. Some authors have considered the model with this coupling term and have studied the early time accelerating expansion of the universe as well as the late time dynamics [59][60][61].
As an important extension of the nonminimal models, in addition to the standard coupling between the scalar field and curvature, it is possible also to have coupling between the derivatives of the scalar field and gravitational sector which makes the inflationary model more powerful to be a successful scenario.
In this work, we consider a generalized inflationary model in which both the scalar field and its derivatives are coupled to the curvature. We also consider an extension of the nonminimal derivative term in the Lagrangian as F( ) ] ∇ ∇ ] ( ). In this term, F is a general function of which, for F ∼ (1/2) , leads to the simple cases introduced earlier in the literature. In Section 2, we present some main equations of this generalized inflationary model. In Section 3, we use the ADM metric and study the linear perturbations of the model. By expanding the action up to the second order in perturbation and considering the 2-point correlation functions, we obtain the amplitude of the scalar perturbation and its spectral index. Also, by considering the tensor part of the perturbed metric, we obtain the tensor perturbation and its spectral index as well. In Section 4, by expanding the action up to the cubic order in perturbation, we explore the nonlinear perturbation in this model. By considering the 3-point correlation functions, we study the non-Gaussian modes of the primordial perturbations. In this section, we obtain the amplitude of the non-Gaussianity in the equilateral and orthogonal configuration and in 1 = 2 = 3 limit. Then, in Section 5, we test our generalized inflationary model in confrontation with the recently released observational data. To this end, we consider two types of nonminimal coupling parameter ( 2 and 4 ), several types of nonminimal derivative coupling function, and also several types of potential. By these functions, we study the behavior of the perturbative parameters in the background of the Planck 2015 observational data. Our main goal is to see whether this nonminimal model has the potential to make models of chaotic inflation in better agreement with Planck 2015 data or not.
In this regard, the model with nonminimal coupling between the scalar field and its derivative to gravity is a subset of Horndeski theory which satisfies the necessary conditions for being a good theoretical inflationary model. Also, this model is observationally viable, since this model with several inflationary potentials is consistent with Planck 2015 observational data. As we can see from Planck team's released paper, a simple inflationary model, with potential corresponding to 2 , 3 , and 2/3 , is not consistent with observation. In [53], it has been shown that a nonminimal model (a model with nonminimal coupling between the scalar field and gravity) is viable with potentials they have considered. However, their model has been studied in Einstein frame. We will show that adding a nonminimal coupling between the gravity and derivatives of the scalar field leads to a viable inflationary model even in Jordan frame.
The main properties of the cosmological perturbations are described by the tensor-to-scalar ratio and the scalar spectral index. In this regard, several observational teams are trying to find some constraints on these perturbative parameters. WMAP team has found the constraints < 0.13 and = 0.9636 ± 0.0084 from the combined WMAP9 + eCMB + BAO + H 0 data [62]. The constraints obtained by the joint Planck 2013 + WMAP9 + BAO data are as < 0.12 and = 0.9643 ± 0.0059 [63]. Then, in 2014, BICEP team found a surprising result. The joint Planck 2013 + WMAP + BICEP2 + BAO data has set the constraints = 0.2096 +0.0443 −0.0608 and = 0.9653 ± 0.0129 on the perturbative parameters [64,65]. This large value of the tensor-to-scalar ratio, which shows the large amplitude of the gravitational wave modes generated during inflation, was a really challenging result. But Planck collaboration has performed a new study on the polarized dust emission in our galaxy and has shown that the part of the sky observed by BICEP2 is contaminated significantly by galactic-dust emission. So, the polarized dust emission might contribute a significant part of the B-mode detected by BICEP2. The Planck collaboration released the constraints < 0.099 and 0.9652 ± 0.0047 from Planck TT, TE, and EE + lowP + WP data (note that Planck TT, TE, and EE + lowP refer to the combination of the likelihood at > 30 using TT, TE, and EE spectra and the low multipole temperature polarization likelihood) [10][11][12]. In this regard, to compare our model with observational data, we study the behavior of the tensor-to-scalar ratio versus the scalar spectral index in the background of the Planck TT, TE, and EE + lowP data and deduce some constraints on the parameters space of the model. We also analyze the non-Gaussianity feature of the model by studying the behavior of the orthogonal configuration versus the equilateral configuration in the background of the observational data.

The Setup
The four-dimensional action for a cosmological model where a scalar field is nonminimally coupled to the Ricci scalar, in the presence of a nonminimal derivative coupling between the scalar field and gravity, is given by the following expression: where is the Ricci scalar, is a scalar field (inflaton) with the potential ( ), and ( ) and F( ) are general functions of the scalar field. Also, in the nonminimal derivative term of action (1), ] = ] − (1/2) ] is the Einstein tensor. Note that, in some papers, there is a coefficient * /2 in front of the nonminimal derivative term, where the constant parameter * has length-squared dimension. In this paper, we absorb such a coefficient in function F( ). By assuming the FRW metric, variation of action (1) with respect to the metric leads to the Friedmann equation of the model as where a dot refers to a time derivative of the parameter and a prime denotes a derivative with respect to the scalar field. By varying action (1) with respect to the scalar field, we get the following equation of motion: Now, by defining the slow-roll parameters as ≡ −̇/ 2 and = −(1/ )(̈/̇), we find these parameters in our setup as follows: where parameter is defined as During the inflationary era, the evolution of the Hubble parameter is so slow, so that in this era the conditions ≪ 1 and ≪ 1 are satisfied. As one of these two slow-varying parameters reaches unity, the inflation phase terminates. The number of e-folds during inflation which is defined as = ∫ hc (6) in our setup and within the slow-roll limit (̈≪ |3| anḋ 2 ≪ ( )) takes the following form: where hc shows the value of the inflaton at the horizon crossing of the universe scale and denotes the value of the field when the universe exits the inflationary phase. In the next section, we study the linear perturbation of the model which is the key test of any inflationary model. In this regard, we calculate the spectrum of perturbations produced due to quantum fluctuations of the fields about their homogeneous background values.

Linear Perturbation
In this section, we study the linear perturbation arising from the quantum behavior of both the scalar field and the spacetime metric, ] , around the homogeneous background solution.
To this end, we should expand action (1) up to the second order in small fluctuations. In this regard, it is convenient to work in ADM metric formalism [66]: where is the lapse function and is the shift vector. To obtain a general perturbed form of metric (8), we should expand the shift and laps functions as ≡ = + V and = 1 + Φ, respectively. Φ and are 3-scalar and V is a vector satisfying the condition V , = 0 [67,68]. Note that it is sufficient to compute or up to the first order. This is because the second-order perturbation is multiplied by a factor which is vanishing using the first-order solution. Since the third-order term is multiplied by a constraint equation at the zeroth-order obeying the equations of motion, the contribution of this term also vanishes [9,15,69]. ℎ should be written as Advances in High Energy Physics the spatial curvature perturbation and T is a spatial shear 3tensor which is symmetric and traceless. So, perturbed metric (8) takes the following form: In order to study the scalar perturbation of the theory, we work within the uniform-field gauge with = 0 and also neglect the gravitational wave (T = 0). So, by considering the scalar part of the perturbations at the linear level, we rewrite perturbed metric (9) as [67,68,70] Replacing the metric (10) in action (1) and expanding the action up to the second order in perturbations give In this gauge, the equations of motion for Φ and , obtained by varying action (11), lead to the following constraint equations: By substituting (12) into action (11) and integrating by parts, we reach where the parameters W and 2 (which is known as the sound velocity) are defined as Advances in High Energy Physics 5 respectively. To see more details of driving this type of equations, one can refer to [15,16,22,71]. Now, we calculate the quantum perturbations of Ψ. To this end, by variation of action (14), we find the equation of motion of Ψ asΨ To solve this equation, it is convenient to rewrite it in terms of the conformal time. By introducing the new variable as = Ψ , where = √ 2W, (17) becomes where a prime refers to derivative with respect to the conformal time. By considering the slow variation of W, we have / ≃ 2/ 2 in the de Sitter background. In this regard, the solution of (18) is given by the following expression: By choosing the Bunch-Davis vacuum characterized by the mode function = exp(− )/√2 in the asymptotic past, when the -mode is far inside the horizon ( → −∞), we obtain In this way, the solution of (17), up to the lowest order of the slow-roll approximation, is given by the following expression: There are two other treatments to obtain the constants 1 and 2 which have been mentioned in [72]. To study the power spectrum of the curvature perturbation, we should compute the two-point correlation function in our setup. We find the two-point correlation function, sometimes after the horizon crossing, by obtaining the vacuum expectation value of Ψ at = 0 (corresponding to the end of inflation): where A is called the power spectrum and is given by The scalar spectral index of the perturbations at the Hubble crossing is defined as which in our setup takes the following form: Any deviation of from the unity shows that the evolution of perturbations is scale dependent.
Other important parameters in an inflationary model are amplitude of the tensor perturbation and its spectral index. To study the tensor perturbation, we use the tensor part of perturbed metric (9). We can write T , in terms of the two polarization tensors, as follows: where (+,×) are symmetric, traceless, and transverse. The normalization condition imposes the constraints and the reality condition gives Now, we can write the second-order action for the tensor mode (gravitational waves) as follows: where W and 2 are defined with the following expressions, respectively: By following the strategy as applied for the scalar mode case, we can find the amplitude of the tensor perturbations as follows: which, by using the definition of the tensor spectral index Advances in High Energy Physics leads to the following expression for : The ratio between the amplitudes of the tensor and scalar perturbations (tensor-to-scalar ratio) is another important inflationary parameter which in our setup is given by

Nonlinear Perturbations and Non-Gaussianity
Non-Gaussianity of the primordial density perturbations is another important aspect of an inflationary model. To explore the non-Gaussianity of the density perturbations, one has to study the nonlinear perturbation theory. The twopoint correlation function of the scalar perturbations gives no information about the non-Gaussianity of the model.
The non-Gaussian property shows itself up in the threepoint correlation function. This is because, for a Gaussian perturbation, all odd -point correlators vanish and the higher even -point correlation functions can be expressed in terms of sum of products of the two-point functions. To calculate the three-point correlation function, we should expand action (1) up to the cubic order in the small fluctuations around the homogeneous background solution. The obtained cubic terms lead to a change in both the ground state of the quantum field and nonlinearities in the evolution. By expanding action (1) up to the third order in perturbation, we obtain a complicated expression which can be seen in the Appendix. We use (12) to eliminate the perturbation parameter Φ in the expanded action. Then, by introducing the new parameter X as we can rewrite the third-order action, up to the leading order, as follows: Now, we are in the position to calculate the three-point correlation function by using the interacting picture. In this picture, we obtain the vacuum expectation value of the curvature perturbation Ψ for the three-point operator in the conformal time interval between the beginning and end of the inflation as [9,22,71] where the interaction Hamiltonian, int , is equal to the Lagrangian of the third-order action. To solve the integral of (38), we can approximate the coefficients in the brackets of the Lagrangian (37) to be constants since these coefficients would vary slower than the scale factor. By solving the integral, we find the three-point correlation function of the curvature perturbation in the Fourier space as where A in (40) is given by (23). Also, the parameter E Ψ is defined as where As it can be seen from (40), a three-point correlator depends on the three momenta 1 , 2 , and 3 . To satisfy the translation Advances in High Energy Physics 7 invariance, these momenta should form a closed triangle, meaning that the sum of these momenta should be zero ( 1 + 2 + 3 = 0). Also, satisfying the rotational invariance makes the shape of the triangle important [21,[73][74][75][76]. Depending on the values of momenta, there are several shapes and each shape has a maximal signal in a special configuration of triangle. One of these shapes is a local shape [77][78][79][80] which has a maximal signal in the squeezed limit ( 3 ≪ 1 ≃ 2 ). Another shape which is corresponding to the equilateral triangle [81] has a signal at 1 = 2 = 3 . A linear combination of the equilateral and orthogonal templates, which is orthogonal [82] to equilateral one, gives a shape corresponding to folded triangle [83] with a pick in 1 = 2 2 = 2 3 limit. We mention that the orthogonal configuration has a signal with a positive peak at the equilateral configuration and a negative peak at the folded configuration. In order to measure the amplitude of the non-Gaussianity, we define the dimensionless parameter NL , called "nonlinearity parameter," as follows: In this paper, we study the non-Gaussianity in the equilateral and orthogonal configurations. To this end, we should find E Ψ in these configurations. In this regard, we should estimate the correlation between two different shapes. So, following [84][85][86], we define the following quantity: The region of integration is 0 < 1 < ∞, 0 ≤ 2 / 1 < 1, and 1 − 2 / 1 ≤ 3 / 1 ≤ 1. We note that two shapes satisfying the condition |C(B ΨBΨ )| ≃ 0 are almost orthogonal. Now, we introduce a shapeS equi as [20,85] By using (40)-(42) and (44)-(46), we can show that the following shape is orthogonal to (46): Note that, to normalize the quantitiesS equi andS ortho , we have considered that in the equilateral configuration limit (corresponding to 1 = 2 = 3 ) two shapes are equal and we haveS equi =S ortho = 3 (for more details, see [85]). In this regard, and by solving the condition |C(B ΨBΨ )| ≃ 0, we have ≃ 1.1967996. Now, bispectrum (41) can be expressed in terms of the equilateral and orthogonal basis as where the parameters C 1 and C 2 are given by the following expressions: By using (43)-(49), we can obtain the amplitude of the non-Gaussianity in the equilateral and orthogonal configuration as follows: Since the equilateral shape has a maximal signal at 1 = 2 = 3 limit and also an orthogonal configuration has a positive peak at this limit, we rewrite (50) So far, we have obtained the main equations of our setup. In the following section, we examine our model by using the recently released observational data by Planck collaboration and find some constraints on the model's parameters space.

Observational Constraints
In this section, we are going to find some observational constraints on the parameters space of the model in hand. To this end, first we define the extension of the nonminimal derivative term as F( ) ∼ (1/2 ) , which for = 1 leads to the usual nonminimal derivative case (as mentioned in Section 1). We set the values 1, 2, 3, and 4 for . Another general function in this setup is ( ), which shows the nonminimal coupling between the scalar field and the curvature scalar. For this function, we consider two cases: ( ) ∼ 2 and ( ) ∼ 4 . The last general function is the potential term in the action. We adopt various types of potential to explore the linear potential as ( ) ∼ [87], the quadratic potential as ( ) ∼ 2 , the cubic potential as ( ) ∼ 3 , the quartic potential as ( ) ∼ 4 , the axion monodromy's motivated potential as ( ) ∼ 2/3 [88], and the exponential  potential as ( ) ∼ − . With this specification of the general functions, we can analyze the model numerically and find some constraints on the nonminimal coupling parameter . For this purpose, we solve the integral of (7) by adopting the mentioned functions of the scalar field. By solving this equation, we find the value of the inflaton at the horizon crossing in terms of the e-folds number. By substituting hc into (25), (35), and (51), we can find the scalar spectral index, tensor-to-scalar ratio, and the amplitude of the orthogonal and equilateral configuration of the non-Gaussianity in terms of the e-folds number. Then, we can analyze these parameters numerically to see the observational viability of this setup in confrontation with recently released data. Plotting the figures and analyzing the model, we have rescaled all constant parameters to unity with the choice = = 1. In this regard, we study the behavior of the tensor-to-scalar ratio versus the scalar spectral index and the orthogonal configuration versus the equilateral configuration. The observational parameters are defined at 0 = 0.002 −1 where subscript 0 refers to the value of when it left the Hubble radius during inflation. The results are shown in Figures 1-8. We have found that an inflationary model with a nonminimal coupling between the scalar filed and Ricci scalar and also a nonminimal derivative coupling, in some ranges of nonminimal coupling parameter , is consistent with Planck 2015 dataset. In fact, these nonminimal couplings help to make models of chaotic inflation that would otherwise be in tension with Planck data, in better agreement with the data. The ranges of nonminimal coupling are shown in Tables      term and studied the model numerically. We can compare the result of our setting with their results here. For example, it is shown in [41] that the presence of nonminimal coupling and noncanonical kinetic terms lead to a chaotic inflation model, which is favored by the BICEP2 observation but is in tension with recent observational data. Here, we have shown that the presence of nonminimal coupling between the scalar field and Ricci scalar and also a general nonminimal derivative coupling between the scalar field and gravity make the model observationally more viable. This model with these couplings is consistent with Planck 2015 TT, TE, and EE + lowP data. The authors of [42] have considered a nonminimal coupling between a scalar field and gravity as ( 2 + 4 ) in Jordan frame. They perform a frame transformation to Einstein frame and compare their model with Planck 2015 data. Their model just for = 65 and with some maximum values of is consistent with observational data and for = 50 and 60 is observationally disfavored. In our setup, we have considered nonminimal coupling functions 2 and 4 separately. Our numerical analysis shows that the presence of NMC and NMDC makes the model observationally viable for both = 50 and = 60 in the Jordan frame (note that although we do not present here the results with = 50, we have analyzed the model with this value of number of e-folds and have found the viability of the model). In [43], a nonminimally coupled inflationary model with quadratic potential in Einstein frame is considered in the background of Planck 2015 observational data. Their numerical analysis shows that this model for some value of NMC parameter gives ≃ 0.96. Interestingly, our setup also, in some values of , gives ≃ 0.96. Nevertheless, we should emphasize that, for some other choices for model parameters, this would be not the case necessarily.
Our numerical analysis shows that this nonminimal model with small values of the NMC parameter is consistent with observational data and so this model in Jordan frame in weak coupling limit has cosmological viability. In week coupling limit, our setup is applicable to several inflationary potentials. In this limit and in the ranges obtained by numerical analysis, the scalar spectral index is red tilted. Another point is that, in our setup, in some ranges of NMC parameter, it is possible to have large non-Gaussianity, as can be seen in figures. Note that, to obtain the constraints on , we have focused on those ranges of the NMC parameter that lead to the values of the observable quantities allowed by the Planck 2015 observational data. For instance, we have focused on the values of leading to 0.95 < < 0.97 allowed by Planck 2015 TT, TE, and EE + lowP data. In the same way, we focused on the values of the NMC parameter leading to −147 < equi NL < 143 allowed by Planck 2015 TTT, EEE, TTE, and EET data. Also, (25), (35), and (51) show that the scalar spectral index, tensor-to-scalar ratio, and the amplitudes of non-Gaussianity have nonlinear dependence on the NMC parameter . The function ( ) has an explicit nonminimal Advances in High Energy Physics 13   coupling dependence. In the mentioned equations, there are different powers of the function ( ) and therefore different powers of the nonminimal coupling parameter. Because of complexity of equations, we have avoided rewriting the equations in terms of the NMC parameter . However, we note that the effects of these nonlinearities in have been included in our analysis and graphs. Since quantities such as the scalar spectral index, tensor-to-scalar ratio, and the amplitudes of non-Gaussianity have nonlinear dependence on the NMC parameter , any small change in the values of has considerable effect on the values of these quantities. In fact, there is a nonlinear dependence on model parameters that really demands four figures of accuracy in quantities such as to get two figures of accuracy in observable quantities such as . So, we have tried to be more accurate and found the values of with higher precision.