Reheating mechanism of Curvaton with Nonminimal Derivative Coupling to Gravity

In this paper, we continue our study on the curvaton model with nonminimal derivative coupling (NDC) to Einstein gravity proposed in our previous work, focusing on the reheating mechanism. We found that according to whether the curvaton has dominated the background after the end of inflation, it will have two different behaviors of evolution, which should be the general property of curvaton with nonminimal couplings. This will cause diffferent rates of particle creation, which goes on via the parametric resonance process. The reheating temperature is estimated for both cases in which reheating completes before and after curvaton domination, and the constraints are quite loose compared to that of overproduction of gravitino. Finally we investigated the evolution of curvature perturbation during reheating. We have shown both analytically and numerically that the curvature perturbation will not blow up during the resonance process.

As an interesting application, in previous works [19,20], we have proposed a curvaton model with such a NDC term. This term possesses a "self-modulation" mechanism that can phenomenologically make the power spectrum of perturbations scale-invariant independent of background evolutions. This is because that in this case the kinetic term of the field couples directly to the geometric variables such as Hubble parameter, so that it can compensate the deviation of spacetime from de-Sitter and make the field "feel" itself in the inflationary background, even if it actually is not. In [19], we calculated the background as well as both scalar and tensor perturbations generated by this model, and provided several constraints on parameters considering its transfer to the curvature perturbations. While in [20] we investigated the full types of non-Gaussianities up to the 3-rd order.
As a consequent work, in this paper we consider another region, namely the reheating process of this curvaton. Reheating is a very important region in the evolu-tion of the universe, for it explains how particles and light elements can be generated, after the dilution of the inflation. The oldest discussions of reheating can be pursued to the last decades of last century, see [21][22][23][24] for original ideas and see [25,26] where the great breakthrough of parametric resonance mechanism was proposed for efficient reheating of single field inflation. As more and more models for early universe come out, more and more reheating mechanisms are also proposed, such as geometric reheating [27], curvaton reheating [28], modulated reheating [29], bounce reheating [30], and so on. Moreover, Recently there has been some discussions on the reheating mechanism for the inflaton with such non-minimal derivative coupling (NDC) terms [31][32][33][34][35][36]. These works raised a couple of questions and discussions about reheating in such models, which is interesting and related to the topic of this paper. For reviews on reheating mechanisms, see [37,38].
In the following, we will investigate in detail the reheating process caused by the NDC curvaton model. For simplicity, we assume that the background is still given by inflation. We will mainly focus on the following questions: • How will the curvaton evolve after the end of inflation? • Can it give efficient mechanism for particle creation (namely, does parametric resonance exist)?
• How is the reheating temperature constrained? • How will the curvature perturbation be affected by reheating process? and by our study, we wish to discover interesting properties of nonminimal coupling curvaton reheating, and find the difference from other kinds of reheating processes, such as that caused by inflaton itself, or other minimal coupling curvatons.
Our paper is organized as follows: In Sec. II we briefly review the NDC curvaton model. In Sec. III we study the reheating process of our model. In Sec. IIIA we show the background evolution both before and after curvaton domination, in Sec. IIIB we study the process of parametric resonance, and in Sec. IIIC we estimate the reheating temperature. In Sec. IV we discuss about the evolution of curvature perturbation in reheating process. Sec. V comes the conclusion.

II. THE NDC CURVATON MODEL
Following the preceding work [19], one has the action including the NDC curvaton as where ξ is the dimensionless coupling constant of the NDC term. L bg is the Lagrangian of the background, which drives inflation, and we don't need to specify its detailed form. It is straightforward to write down the equation of motion for the curvaton field ϕ, such as 6ξ and its energy density and pressure can be expressed as The background evolution of the curvaton field with various types of potential as well as the linear perturbations during inflation has been classified and briefly analyzed in [19]. As has been shown there, if the curvaton field is massless, namely V (ϕ) = 0, exactly scaleinvariant power spectrum will be obtained due to the nonminimally kinetic coupling. However, since exact scale-invariance is not favored by today's data, a nonzero potential is needed which could give a mass term to ϕ. 1 For the simplest choice, we choose as its potential. Then the equation of motion (2) reduces toφ where the slow-roll parameter ǫ is defined as ǫ ≡ −Ḣ/H 2 . Moreover, one can define an effective mass of the curvaton field,m 1 For varying Hubble parameter in non-inflationary case, a nearlyconstant correction could also be obtained, with a varying mass of ϕ.

III. REHEATING MECHANISM
A. the behavior of curvaton during oscillation In order to investigate the reheating mechanism of the curvaton field ϕ, let's assume that after the inflation ends, the background field (inflaton) decays rapidly without oscillation. A typical example is that the inflaton field is a canonical scalar field dominated by its kinetic energy, such as the so-called "quintessential inflation" [40]. Therefore one has ρ bg ∝ a −6 (t) , w bg = 1 .
However, as the curvaton field has a potential (5), it will fall down and oscillate around its minimum, and if the energy density of curvaton decrease slower than that of inflaton, the curvaton can dominate the universe, unless it reheats before domination. First of all, one can determine the averaged equation of state during its oscillation, which will be used for later analysis. To do this, one can parametrize the scale factor a(t) and Hubble parameter H(t) during the rapid deflation of the universe. One can have: and w should be equal to w bg before curvaton domination, and w ϕ after that. For convenience, we define t eq as the time when the energy densities of inflaton and curvaton are equal. For the former case, p = 1/3 as w bg = 1, while the latter case depends on the behavior of ϕ which is the solution of the equation of motion (2). Considering (11), Eq. (2) can be solved as: where Φ(t) = Φ 0 t (1−3p)/2 is the oscillation amplitude of ϕ with the initial condition of Φ 0 , and J is the Bessel function. Note that in reheating era where ϕ begins to oscillate, one hasm > H, which roughly gives mMt 2 2 √ 6ξp > 1. Thus the Bessel function can be mimicked by the trigonometric functions, namely During rapid oscillation, the kinetic and potential terms in the energy density (3) are of the same order, therefore one roughly has For t < t eq , p = 1/3, so we have ρ ϕ ∝ t 0 ∼ a 0 . This means that before curvaton domination, the energy density of curvaton field will behave as a constant, giving its equation of motion to be w ϕ = −1.
Since the background energy density is decreasing, the curvaton will exceed the background (at the time t eq ) and dominate the universe. For t > t eq , the same solution of (13) applied but with p be related to w ϕ , so one cannot determine p, or the behavior of ρ ϕ , solely by the relation (14). However, since the curvaton dominates the universe, we have another equation, namely the Friedmann equation, This equation, together with (11), tells us that ρ ϕ ∝ t −2 .
Comparing it with Eq. (14), one gets p = 1, which gives In order to confirm our analysis, we performed numerical calculations for both the two cases. One can see that in Fig. 1, both the field value ϕ and its energy density ρ ϕ oscillate with a stable amplitude, giving a constant averaged value. However, in Fig. 2, the field value oscillates with a damping amplitude. The Hubble parameter is therefore also oscillating, so the value of p is also an oscillating function with its averaged value p = H t. From the figure it is clearly seen that p is at the position of about unity. We also plot velocities of ϕ and H for later use.
As a side remark, we remind the reader that our result is different from that for minimal coupling curvatons [28]. In the latter case, if the curvaton oscillates (let's still take a mass-squared potential for curvaton as an example), the solution of curvaton will be as ϕ = t −3p/2 cos(mt), and ρ ϕ ∼ t −3p ∼ a −3 , accordingly. One can see that, no matter what value p will be given, or no matter whether curvaton is dominant or not, the curvaton will behave like an ordinary matter, with w ϕ = 0 averagely, as long as the mass-squared potential is applied. This means the energy density will evolve synchronically with the scale factor, in terms of t. This result is consistent with [28]. However when curvaton is nonminimally coupled, this synchronicity will be violated, and the behavior before and after domination will be different. As will be seen later, this different will directly cause the different rate of particle creation. To our knowledge, this has not been given much notice in previous studies.

B. parametric resonance
In this section, we begin to discuss about how the oscillating behavior of the curvaton ϕ, as shown above, can be responsible for creating particles (we denote as χ), namely, the preheating process. One assume the freefield Lagrangian of χ has a canonical form, and that ϕ can interact with χ through some interacting term, such as gϕ 2 χ 2 , with a dimensionless coefficient of g. It gives the total form of Lagrangian of χ field as: It is true that one can also have the interaction term of χ and the background field (the inflaton), however here we drop this out for simplicity. Actually, as the inflaton decays quickly after inflation ends, it will have much less effects on χ than ϕ. Since χ is a quantum field, one can write down the Heisenberg presentation of χ as [25]: k is the annihilation and creation operators respectively, and k is the comoving wavenumber. Using this, one can also get the equation of motion for χ by simply varying the Lagrangian (16): For the rest of this paper, we assume m χ = 0 for simplicity. Since in our case, ϕ(t) behaves as in Eq. (13), in the static limit where a = const., H = 0, Eq. (18) will reduce to the Mathieu-like equation, the solution of which is the well-known Floquet solution [25]. In expanding phase such as inflation, however, one can define a new variable: X k (t) ≡ a 3/2 (t)χ k , and Eq. (18) can be rewritten as: One can apply the WKB approximation to get the solution of Eq. (19): with α k and β k being the time-dependent coefficient of positive and negative frequency parts, respectively. Moreover, the α k and β k satisfy the equations: with additional normalization relation Furthermore, the comoving occupation number of χ particles in the mode k is defined by Using Eq.s (20), (21) and (22), one get a simple expression of n k = |β k | 2 , and the vacuum expectation value of particle number density of χ per comoving volume is: It is straightforward to know how the χ particles increase with time just by integrating |β k | 2 , which is the solution of Eq. (21). However, it is rather difficult, if not impossible, to have analytical solution of (21) due to the complication of the dispersion relation in (20). In order to solve the problem, Ref. [25] introduced a method, that is, to consider the solution of each period around the time point at which the source field of particle creation is zero, namely t j , j = 1, 2, 3, ... where ϕ(t j ) = 0 in our case. This is because that only around t j the particle number can be dramatically changed. If we Taylor expand the function gϕ 2 (t) with ϕ given in (13) around t j , we can get: So the leading order of gϕ 2 (t) gives a parabolic potential for χ, and the particle creation at t j could be viewed as the scattering of χ particles through this potential for the j-th time. We assume the solution of X k at the period of t j−1 < t < t j is with α j k and β j k approximately constants. Following [25] and noticing the relation of n k = |β| 2 , one can get a recursive relation of n j k and n j+1 k as: with µ j k ≡ (2π) −1 ln(1+2e −πκ 2 +...), where the ellipsis denotes the random (stochastic) terms. Here κ ≡ k/(ak * ) and k * = g/(6ξ)mM t 3(1−p)/2 j /p. In [25], k * is independent on t j , but as is mentioned before, in our case there are two differences. One is that p does not have the same value before and after curvaton domination, and the other is that the scaling of ϕ on t j is different from the minimal coupling case, as was shown in (13). These causes k * ∝ t 1/2 j for t j < t eq , while k * ∝ t 0 j for t > t eq . Considering also the dependence of a on t, namely a ∝ t 1/3 before t eq and a ∝ t after, one gets κ ∝ t −5/6 before t eq and κ ∝ t −1 after. Therefore one gets: In the limit of κ 2 ≪ π −1 , µ k reduces to the trivial case in [25]: µ j k ≃ (2π) −1 ln 3 ≃ 0.18 where we neglected the random terms.
The total number during the whole process of resonance yields in mode k. In principle one can get n k by integrating µ k according to the expressions (28), however unfortunately it might not be analytical, so we will refer to numerical calculations later on. The particle number density in real space is: while the increasing rate of the particle number in the comoving volume is roughly estimated as: In Fig.s 3 and 4 we numerically calculated the creation of particles before and after curvaton domination, and plot the values of X k and ln n k in terms of t. The numerical results are in good agreement with the theoretical analysis. In our plot we can see that in both cases X k oscillates rapidly due to the parametric resonance, and with the amplitude get higher and higher. The number density n k (in Logarithm) of both cases increase rapidly, though in some instances it decreases, due to the stochastic random process. Nevertheless, one can see that there are quantitative difference in the shapes as well as the increasing rates of X k and n k in the two cases, which is due to the different behavior of the source field ϕ before and after curvaton domination. As a side remark, we roughly estimate the time when back-reaction of the produced particle χ becomes important, which can prevent further reheating process. To see this, one could define the energy density of the χ as: where n χ is calculated in (30) while the effective mass of the χ particle is from Eq. (19). The condition of back-reaction becoming important is ρ χ ≈ ρ ϕ , where ρ ϕ ≈ m 2 ϕ 2 . Since it is difficult to have an analytical solution of the time t, we perform the numerical calculation, and plot both ρ χ and ρ ϕ (after neglecting the oscillation/resonance effect) in

C. constraint on reheating temperature
When reheating process completed, the curvaton will decay into the relativistic products, which scales as T 4 where T is the temperature. Therefore the reheating temperature will be related to the final state of the curvaton field. Knowing the initial state of the curvaton field at where inflation has just ended, and the scaling of the curvaton field during reheating, one can estimate the reheating temperature and compare it to various data constraints. This procedure is initially done for normal curvaton mechanism given in [28]. Now we also use this method to estimate the reheating temperature in our model.
The curvaton may reheat the universe in two ways. Since we know from above that the energy density of curvaton scales as a 0 before domination, we have: where the subscripts osc and eq denotes the value when ϕ begins to oscillate and when it has the same amount of energy density with background, respectively. On the other hand, at the time when reheating completes we have: Therefore, if the curvaton reheats after its domination, then we have ρ eq ϕ > ρ rh ϕ , and the amount of difference depends on how long reheating will last after curvaton domination. From Eq.s (34) and (35) we have: At the time when ϕ ends slow-rolling and begins oscillating, the kinetic and potential term of ϕ will be of the same order, namely 9ξH 2φ2 /M 2 ∼ m 2 ϕ 2 /2, therefore where y osc is the value of y(≡ ξφ 2 /M 2 ) at the oscillation time. So the reheating temperature is: The constraints on the Hubble parameter and y can be obtained from the Planck constraints on scalar spectrum: |P ζ | ≃ 2.2065 +0.0763 −0.0738 × 10 −9 (TT, TE, EE+lowP, 68% CL) and the tensor/scalar ratio: r 0.10 (TT, TE, EE+lowP, 95% CL) [39]. In our model, from [19] we have: for this case, where " * " denotes the value evaluated at horizon crossing. These first give the constraints: from where one can see, if we set ǫ * ∼ O(10 −2 ), we can get y * to be as low as 10 −7 M 2 p . As y and H varies slowly during inflation, one can roughly have H * ≃ H osc , y * ≃ y osc , therefore the reheating temperature can be set to be: which actually gives very loose constraint on the reheating temperature, comparing to that from the gravitino producing [41].
On the other hand, if curvaton decays and reheats before dominating the universe, one has ρ rh ϕ ≃ ρ osc ϕ ≃ ρ eq ϕ . From Eq.s (35), (37), (38) we have: For perturbation in this case, we have: where r ≡ ρ ϕ /ρ r is the ratio of energy densities of ϕ field and the radiation. From Planck constraints, we have so we have which gives a high reheating temperature. So if the constraints on gravitino are trustable, then this case will be ruled out 3 . Note that recently, in Ref. [42], the authors obtained a more stringent upper bound for reheating temperature of T rh < 5 × 10 4 GeV for minimal coupling curvaton models by Bayesian inference method 4 . Moreover, one can also constrain some of the parameters for the model, using that on Hubble parameter. From Eq. (9), one has H * ∼ 10mM/ √ 6ξ during inflation while H osc ≃ (m 2 M 2 /6ξ) 1/4 when ϕ begins to oscillate, which is given by H osc ≃m. Since H * ≃ H osc , these two conditions actually do not differ much. From either (40) or (44), one can get the relation of mM/ √ 6ξ 10 −10 M p . That means, if we choose m to be of 10 −6 M p , we will get M/ √ 6ξ 10 −4 M p , which is consistent with our numerical calculations.

IV. EFFECTS ON CURVATURE PERTURBATIONS
In this section, we discuss about how the resonance process could affect on the final curvature perturbations in our model, especially whether the points ofφ = 0 will make the curvature perturbations diverge. This issue has been in debate in early works of canonical inflation reheating [43,44] (see also [45][46][47] for further discussions), and has been revisited recently for NDC inflation reheating [36]. As we know, since there are more than one component in the universe, the perturbations generated by curvaton is basically isocurvature, which will transfer to curvature perturbation after curvaton domination. According to the original analysis on curvaton in [48], the curvature perturbation is defined as: in spatial flat gauge. Assuming that the perturbations of the background field (inflaton) are negligible, we have δρ bg ≈ 0. From continuity equations for both components,ρ = −3H(ρ + P ). For background whose equation of state is assumed to be unity, we haveρ bg = −6Hρ bg . Therefore for the time before curvaton dominating, one roughly have: Since ρ bg are non-zero, one can easily see that ζ will not diverge.
On the other hand, for the time after curvaton dominating, one have: where we have made use of Eqs. (3), (4) as well as the relation: So it seems that whenφ passes through 0 during oscillation, ζ might diverge. However, this is not true. From equation of motion for δϕ [19]: where z = a √ Q/c s ∼ aH, and Q, c s as well as m ef f have all been defined in [19]. Making use of (48), we have: (51) On the other hand, from Eq. (2), one has: where V ,ϕ = m 2 ϕ, and one more time derivative gives: Expanding Eq. (51) and making use of the above two equations, one finally have: whenφ = 0. It has the unique solution:ζ = (2Ḣ/H)ζ. However, since now the curvaton is dominant, we havė which is equal to zero whenφ = 0. This can also be seen in Fig 2. So we still haveζ = 0, and ζ will not diverge. Note that similar conclusion has been obtained in the early work of canonical single field inflation [44] 5 .
We performed numerical calculations for both the two cases and plot the evolutions of ζ in Figs. 6 and 7. One can see that, in Fig. 6, ζ oscillates with an increasing amplitude. This can be explained as that, the perturbation generated before curvaton domination is of isocurvature type, which can be a source of ζ. However in Fig. 7, ζ behaves as a nearly constant. This is because when curvaton dominates and the isocurvature perturbation has been transformed into the curvature one, ζ will be a conserved quantity. Although there does have some features on ζ, which may be due to the secondary effects of the resonance, numerical fluctuation or other unknown reasons, it will not be so bad to make it diverge. We also plotζ in Fig. 8, so one can see more clearly that, at the range of typical amplitude of ζ (∼ 10 −5 ), the variation of ζ can hardly be seen. 5 However, since the parameter sound speed squared c 2 s contains ǫ ∼Ḣ [19,35], in the reheating era when H is oscillating,Ḣ will also oscillate between positive and negative values, which causes negative c 2 s and gradient instability for large-k modes during reheating [35]. Since now we're interested the fluctuations outside horizon, the instability will not affect too much. We thank the anonymous referee for pointing this to us.

V. CONCLUSIONS
In this paper, we investigated the reheating mechanism of the curvaton model nonminimally derivative coupled to gravity. We assume that the curvaton has a mass squared potential term, which can not only explain the observationally required tilt of the perturbation spectrum from scale invariance, but also provide a minimum around which the curvaton can oscillate. This gives an efficient way of particle creation, namely parametric resonance process. What is different from the reheating process of inflaton is that, the background evolutions before and after curvaton domination are not the same, and if the curvaton is nonminimally coupled, the averaged decaying rate of curvaton may also not be the same, which can affect the rate of particle creation. Though we didn't provide a rigid proof (which will be left for future work), this may be a general phenomenon.
For our model, as a specific case, we found that the averaged value of curvaton scales as a constant before curvaton domination, while proportional to time inverse after that. The particles, whose field interacts with the curvaton, will indeed be effectively created due to the resonance process. To confirm our analysis, we also numerically calculated the background evolution for the curvaton field and the increase of created particle numbers by the resonance.
Moreover, we also estimated the constraints on reheating temperature of our model, for both cases which the reheating completes before and after curvaton domination. We only got a upper limit on the temperature, which is looser than that given by the overproduction of gravitino. Finally, we investigated the curvature perturbation generated in our model. We showed that the curvature perturbation will be regular, although it is apparently divergent whenφ goes to zero in each oscillation. Unlike that the amplitude of curvature perturbation is increasing before curvaton domination, which should be due to the source of isocurvature perturbation, after the curvaton domination it will remain an approximately constant.
One can also investigate the back reaction of the created particle χ on the reheating process, which we leave for a future study. As an application, our work can be used either to compare cases of single and multiple field theories nonminimally coupled to gravity, or to compare cases of minimal and nonminimal couple theories in fixed degrees of freedom.