Quintessential inflation at low reheating temperatures

We have tested some simple quintessential inflation models, imposing the requirement that they match with the recent observational data provided by the BICEP and Planck team and leading to a reheating temperature, which is obtained via gravitational particle production after inflation, supporting the nucleosynthesis success. Moreover, for the models coming from supergravity one needs to demand low temperatures in order to avoid problems such as the gravitino overproduction or the gravitational production of moduli fields, which are obtained only when the reheating temperature is due to the production of massless particles with a coupling constant very close to its conformal value.


Introduction
Quintessential inflation [1] is a good candidate to unify the early and late time acceleration of our universe. These models are a combination of an inflationary potential, used to explain the early acceleration of the universe, and a quintessential one -which could only be a cosmological constant -, which takes into account the current cosmic acceleration. At early times the inflationary acceleration is the one that dominates and it ceases to be dominant in a phase transition where the universe enters in a kinetic regime [2]. At this moment, particles coupled with gravity are produced to reheat the universe and match with the current hot one. Finally, at very late times, the quintessential potential dominates and the universe starts to accelerate again.
In order to ensure the viability of these models, they are required to fit well with the recent observational data provided by the BICEP and Planck's teams [3,4], but also the reheating temperatures have to be compatible with the nucleosynthesis and baryogenesis bounds, preventing the increase a e-mails: llibert.areste@estudiant.upc.edu; Llibert.Salo@campus.lmu.de b e-mails: jaime.haro@upc.edu of entropy due to the decays of gravitational relics such as moduli fields and gravitinos. It is well known that the success of nucleosynthesis constrains the reheating temperature to be between 1 MeV and 10 9 GeV [5][6][7][8][9][10][11], whereas in order to overcome the overproduction of very light gravitinos (particles which appear in supersymmetric gravitational theories), i.e., with masses in the range 10 −3 MeV m 3/2 10 −1 MeV, it is stated that the reheating temperature has to be less than 10 2 GeV [12]. It is more complicated to avoid the increase of entropy due to the production of moduli fields, which also appears in supergravity models, during the inflationary period; in that case, as it was showed in [13], the reheating temperature has to be less than 1 GeV. However, this low temperature seems to be problematic with baryogenesis, according to which it is commonly assumed that the reheating temperature is at least of the order of the electroweak scale (10 2 GeV). But, fortunately, there are mechanisms that explain the baryon asymmetry at very low temperatures [14]. For all these reasons, we will consider that those supergravity models with a reheating temperature in the MeV regime, i.e., between 1 MeV and 1 GeV, succeed in solving all these problems.
Having this in mind, the main goal of the present work is to study the viability of some well-known inflationary potentials, some of them coming from supersymmetric theories, adapted to quintessence. To do it in a simple way, first of all we consider a universe with a cosmological constant, we choose positive inflationary potentials that vanish at some value of the scalar field and then we extend them to zero for the other values of the field. Hence, we obtain a potential with a phase transition that models an inflationary universe at early times, which reheats the universe via gravitational particle production after inflation and it is finally dominated by the cosmological constant. Once we have these potentials, we calculate their spectral parameters, the number of e-folds, which has to be in quintessential inflation between 63 and 73 as we show, and its reheating temperature in three cases that will be analytically calculated: via the production of heavy massive particles conformally coupled with gravity, massless nearly conformally coupled particles and particles far from the conformal coupling. This last case has also been studied in [15] in the context of braneworld inflation for exponential and power-law potentials. Note that there is another way to reheat the universe, via the so-called "instant preheating" [13], which has been applied to quintessential inflation in [16][17][18][19] for exponential and power-law potentials. This kind of reheating deserves future investigation for different potentials such as the ones studied in this paper, which we will deal with in a future work.
The paper is organized as follows: in Sect. 2 we build a non-singular model based on a universe filled with a barotropic fluid with a non-linear equation of state where the universe starts and finishes in a de Sitter phase. Via the reconstruction method we find the potential that mimics these dynamics when the universe is filled with a scalar field, we find the analytic solution that mimics the fluid model, and we see that it provides spectral quantities such as the spectral index, its running and the tensor to scalar ratio that match with the current observational data. However, this solution appears to be unstable at early times, in the sense that any small perturbation leads to a solution which is singular at early times, in the same way as happens with the Starobinsky model [20,21]. Moreover, with our criterion, the model only supports a reheating via the creation of massless particles nearly conformally coupled with gravity. In Sect. 3, we consider a universe with a small cosmological constant and we adapt three simple inflationary potentials [exponential SUSY inflation potential (ESI), the Higgs inflation potential in Einstein frame (HI), and the power-law inflation potential (PLI)] to quintessence, showing that the HI potential leads to an unacceptable number of e-folds, the ESI potential, which comes from a supersymmetric potential and thus suffers the gravitino and moduli problems, only supports a reheating via the production of massless particles nearly conformally coupled with gravity, and the PLI, which is only viable in the quadratic case, supports the three kind of reheating studied in the work because, since it does not come from any supergravity potential, gravitinos or relic moduli fields do not appear in this theory. In the last section we adapt other well-known potentials [Witten-O'Raifeartaigh inflation (WRI), Kähler moduli inflation I (KMII), open string tachyonic inflation (OSTI), brane inflation (BI) and loop inflation (LI)] coming from inflation, to quintessence and we study its viability. Finally, in the appendix we review from a critical viewpoint the overproduction of gravitational waves in quintessential inflation.
The units used throughout the paper areh = c = 1 and, with these units, is the reduced Planck's mass.

A non-singular model
In this section we are going to study inflation coming from fluids [22,23] and to develop the idea of a non-singular universe proposed in [24], where the universe was filled with a barotropic fluid whose non-linear equation of state (EoS), namely P = P(ρ) P being the pressure and ρ the energy density, satisfies P + ρ = 0 at two different scales, leading, at early and late times, to two de Sitter eras.
The simplest realization of this model, which could be generalized introducing viscosity [25], is to consider three parameters H i H E H f , where H i and H f are two fixed points of the dynamical system and H E is the value of the Hubble parameter when the transition from inflation to kinetic regime is produced. To simplify, we assume that H i = M pl , i.e., the universe starts, in a de Sitter phase, at Planck scales, and in order to reproduce the current cosmic acceleration we have to take H f ∼ H 0 . The simplest way to obtain the model is to consider the following differential equation: where we have chosen k = to ensure the continuity ofḢ . Regarding the value of k, it will be determined in the next section so that our model matches with the current observational data. The dynamical system can be analytically solved leading to the following Hubble parameter: and the corresponding scale factor is Moreover, the effective EoS parameter, namely w eff , which is defined as w eff ≡ P ρ = −1 − 2Ḣ 3H 2 , for our model is given by which shows that for H ∼ M pl one has w eff (H ) ∼ = −1, meaning that we have an early inflationary quasi-de Sitter period. At the phase transition, i.e. when H ∼ = H E , the EoS parameter satisfies w eff (H ) ∼ = 1 and the universe enters in a kinetic or deflationary period [26,27], and finally, for H ∼ = H f one also has w eff (H ) ∼ = −1 depicting the current cosmic acceleration.

Cosmological perturbations
In order to study the cosmological perturbations, one needs to introduce the slow-roll parameters [28] which allow us to calculate the associated inflationary parameters, such as the spectral index (n s ), its running (α s ) and the ratio of tensor to scalar perturbations (r ) defined by where the star ( * ) means that the quantities are evaluated when the pivot scale crosses the Hubble radius. In our case, these quantities become From the theoretical [28] and the observational [29] value of the power spectrum and the observational value of the spectral index obtained by Planck2015 data, n s = 0.968 ± 0.006, we can find the value of H * and k that fits in our model with these data by inserting the value of k = into the first equation in (7). This is equivalent to finding a root to Hence, for the corresponding range of values of n s , we obtain that 3 × 10 −5 < H * M pl < 4 × 10 −5 . Thus, it is straightforward that 7 × 10 −12 < k < 10 −11 and, by using the equations in (7), one gets 0.1040 < r < 0.1520 and −7×10 −4 < α s < −3 × 10 −4 . Therefore, some values of n s in the 1dimensional marginalized 2σ C.L. (0.970 < n s < 0.974) correspond to a ratio of tensor to scalar perturbations which agrees with the constraint r < 0.12, provided by the BICEP and Planck Collaborations [3]. It fulfills, as well, the limits obtained for α s in Planck2015 [4] (α s = −0.003 ± 0.007).
Finally, the number of e-folds is given by where (end) stands for the end of inflation, i.e., end = 1, . Thus, by using the calculated values of k and H * , we obtain a number of e-folds satisfying 52 < N < 77. Now, we are going to compare this value with the one that we will obtain from the following equation [30] in an analogous way to [31]: where R and M symbolize the beginning of radiation era and the beginning of the matter domination era and we have used We use that H 0 ≈ 2 × 10 −4 Mpc −1 and, as usual, we choose a 0 = 1 taking as a physical value of the pivot scale k phys = 0.02 Mpc −1 (value used by Planck 2015 [4]). Then one sees that, in co-moving coordinates, the pivot scale will be k * ≡ a 0 k phys = 0.02 Mpc −1 . Moreover, we know that the process after reheating is adiabatic, i.e. T 0 = a M a 0 T M , and that we have the relations ρ M ≈ π 2 15 g M T 4 M and ρ R ≈ π 2 30 g R T 4 R (where {g i } i=R,M are the relativistic degrees of freedom [32]). Hence, We use H 0 ∼ 6 × 10 −61 M pl and, from Eq. (8), we infer that H * ∼ 4 × 10 4 √ * M pl . We know as well that T 0 ∼ 2 × 10 −13 GeV and g M = 3.36 [32]. Also, g R = 107, 90 and 11 for T R ≥ 135 GeV, 175 GeV ≥ T R ≥ 200 MeV and 200 MeV ≥ T R ≥ 1 MeV, respectively [32]. On the other hand, Thus, we obtain Therefore, with the values in our model and with the range 1 MeV ≤ T R ≤ 10 9 GeV required in order to have a successful nucleosynthesis [33], we find that 63 N 73. With Eq. (9), this range is verified in our values for 0.969 < n s < 0.973.
In conclusion, with the intersection of the bounds obtained from the constraints r < 0.12 and 63 N 73, in our model we have 0.970 < n s < 0.973 with the corresponding values of k ∼ 8 × 10 −12 and H * ∼ 3 × 10 −5 M pl (corresponding to H E ∼ 2 × 10 −6 M pl ), thereby having a number of e-folds of 66 ≤ N ≤ 73. Hence, our model satisfies all the values obtained from recent observations. Moreover, one can see in Fig. 1 that our model also provides theoretical values that enter in the marginalized 2σ C.L. contour in the plane (n s , r ).

The scalar field
In this section we are going to mimic the perfect fluid that fills the FLRW universe with a scalar field. If we represent the energy density and the pressure by the notations ρ ϕ , p ϕ , respectively, then they assume the following simplest forms: Now, using Eq. (14) and the Raychaudhuri equationḢ = −φ 2 2M 2 pl , we find which can be analytically solved, where Now, we aim to analyze the stability of the solution to our model. Using the potential found in Eq. (17), we will study the dynamics of the equation: where . The numerical results are presented in Fig. 2. For clarity in the understanding of the dynamical system, we have used a value of k which differs in various orders of magnitude from the one in our model. Yet the behavior is exactly the same as with the real values. Therefore we observe that the analytical orbit is the only one with two de Sitter points. All the other orbits start with an infinite energy at t → −∞, they experience a phase transition at ϕ = ϕ E and then asymptotically approach the analytical orbit for t → ∞, reaching at t → ∞ the de Sitter point H = H f .
An important remark is in order: From the phase space portrait we can see that there is only one non-stable nonsingular solution, i.e., if one takes initial conditions near the analytical solution one obtains a past singular orbit. This is exactly the same that happens in the Starobinsky model [20] (see Figure 3 of [21] for details).

Reheating constraints
Before studying the reheating constraints in our model, we are going to carry out some useful approximations that will help us in our further calculations. Since k ∼ =  hand, it is also fulfilled that 0 ∼ = α M 2 pl − 2H 2 * and, thus, we obtain all the following approximative expressions: and, hence, r = 16 * ∼ = 16k Moreover, the number of e-folds is but using Eq. (13) we also obtain Equating both quantities one obtains the equation where we have introduced the notation Y = 6 1−n s . Hence, from this equation we obtain the following relation between the cosmological parameters (n s , r ) and T R , showing that the suitable constraints will be verified for 1 MeV ≤ T R ≤ 10 5 GeV, as we can see in Fig. 3. Now, let us consider the gravitational production of χparticles. Firstly, we treat the case when the produced particles are very massive and conformally coupled with gravity. Since a classical picture of the universe is only possible at energy densities less than the Planck one, it seems natural to choose that at this scale the quantum field χ is in the vacuum. In fact, we will choose as a vacuum state the adiabatic one, defined by the modes [34] where k satisfies the equation Then, to obtain an approximate expression of k , one can use the WKB approximation, which holds when m ≥ H . For this reason, since at the Planck epoch H ∼ M pl , in order to obtain the WKB solution one has to assume that the mass of the field, namely m, satisfies m ≥ M pl . However, since for m > √ 4π M pl the produced particles are micro black holes [35], whose thermodynamical description is unknown [36], one has to choose m ∼ M pl so as to avoid this black hole production.
The reheating temperature caused by the decay of these particles into lighter ones, with a thermalization rate -the same used in [5,37] [26], is of the order (see for details [31,38]) (25) and, inserting this expression in (22), one obtains which leads to a spectral index n s ∼ = 0.9716, a tensor/scalar ratio r ∼ = 0.1136 and a reheating temperature T R ∼ 11 GeV. Since this model does not come from supergravity, we do not need to take into account the gravitino and moduli fields problems. Hence, these results lead to a viable model. On the other hand, when considering massless particles nearly conformally coupled with gravity, the reheating temperature becomes [37,39] with ξ the coupling constant, N = 1 dτ H E , and where R is the scalar curvature. It has been numerically computed that N ∼ 1. So as to verify the bounds 1 MeV ≤ T R ≤ 10 5 GeV (coming from the restriction r ≤ 0.12 in our model), the coupling constant ξ must satisfy If we consider now massless particles far from the conformal coupling with gravity, using the results of [40,41] which we will discuss in the appendix, the reheating temperature was calculated in [38], leading to In this case combining Eqs. (22) and (29), one obtains whose solution is n s = 0.970, corresponding to r = 0.119, coinciding thus with the limit where our model starts being valid. Hence, our model supports the presence of these particles.

Simple quintessential inflation models
In this section, we are going to proceed in an inverse way from the one we have just done. We are going to take some simple and well-known potentials, some of them with a similar form from the one already studied, and we are going to verify whether we can adjust the parameters so as to fulfill all the corresponding constraints for the cosmological parameters. To simplify, we will consider some well-known positive inflationary potentials that vanish at some value of the scalar field, namely ϕ E , and we will extend it to zero for the other values of the field. Moreover, we introduce a cosmological constant ∼ 3H 2 0 (H 0 being the current Hubble cosmological constant) to ensure the current cosmic acceleration.

Exponential SUSY inflation (ESI)
The first potential we are going to study is an exponential SUSY inflation (ESI) style potential [42,43], λ being a dimensionless positive parameter. By using the following approximate expressions of the slow-roll parameters as a function of the potential: we obtain * = M pl . Hence, we can compute the associated inflationary parameters, analogously as we did in Sect. 3: where we have usedḢ = −φ 2 2M 2 pl for the calculation of α s . It is also straightforward to calculate the power spectrum: where we have used V (ϕ * ) ≈ 3M 2 pl H 2 * . We will guarantee that the value P ∼ 2×10 −9 is verified in our model by taking λ = . Finally, regarding the number of e-folds, So as to fit our model with the recent observational values, we will initially take, as before, the spectral index obtained by Planck2015 data, n s = 0.968±0.006 and we will find the corresponding value ϕ * from the first equation in (33), which is −4.37 ≤ ϕ * M pl ≤ −3.99. Therefore, the other cosmological parameters are 0.0013 ≤ r ≤ 0.0028 and −7×10 −4 ≤ α s ≤ −3 × 10 −4 . Finally, the corresponding number of e-folds is 49 N 73.
Hence, all the values of r and α s fit to our restrictions and, regarding the number of e-folds, in order to have 63 N 73 we need to take the range of 0.970 ≤ n s ≤ 0.974 for the spectral index, obtaining thus 0.0013 < r < 0.0017 and −5 × 10 −4 < α s < −3 × 10 −4 . Moreover, the theoretical values provided by our model enters in the marginalized 1σ C.L. contour in the plane (n s , r ) as one can see in Fig. 4.
Since the approximative expression of the cosmological parameters is only valid during the inflation process, we cannot directly calculate the Hubble constant at the transition, namely H E . As we can verify in the former model where we had the analytical form of the dynamical system, it is of the same order as at the end of the inflation (H end ). Since end = 1 and ϕ end = ln , meaning that for our model one has H end ∼ 3 × 10 −6 M pl . So as to study the reheating constraints for this model, we are going to proceed as in the former case, i.e., with the approximative expressions r ≈ 2(1 − n s ) 2 ,  . Now we cannot reproduce the same calculation as we did for the former case so as to obtain Eq. (13), because we do not have an analytic expression to perform the calculations given in (12), but since this value is negligible we can continue using it.
Therefore, by combining Eq. (13) and the just obtained approximate value of N , we have again using Y = 6 1−n s . Hence, from solving this equation, one can see in Fig. 5 that the constraints on r and n s are fulfilled, as we had already pointed out. Now, with regards to the case when the produced particles are very massive and conformally coupled with gravity (again choosing m ∼ M pl for the same argument as in previous section), we have to proceed as follows: The square modulus of the β-Bogoliubov coefficient is given by the expression [31] where a E is the value of the scale factor at the transition phase, andḦ ± E are the value of the second derivative of the Hubble parameter after and before the phase transition.
pl ∼ = 32H 6 end . A simple integration leads to an energy density for the produced particles of ρ χ ∼ = 5H 6 end 16 3 π m 2 , and following the calculations made in [31] one gets the following reheating temperature: which leads, by combination with (36), to obtaining that n s ∼ = 0.9720, r ∼ = 1.56 × 10 −3 and T R ∼ 5 GeV. This result again satisfies the gravitinooverproduction problem but not the moduli fields one. In order to study massless particles nearly conformally coupled with gravity, we will use the same formula as used in (27) and, since we do not have the analytical expression of the Hubble constant throughout all the time, we will assume that N ∼ 1 as in the other case. Hence, the reheating temperature is obtaining, thus, the result that for the bounds of T R coming from the nucleosynthesis it should be verified that 10 −7 ξ − 1 6 16. This only gives us a restriction for the lower bound, since given that we have considered the particles to be nearly conformally coupled with gravity ξ − 1 6 10 −1 , which adds the constraint that T R ≤ 7×10 5 GeV. By considering a reheating temperature in the MeV scale, one would obtain 10 −7 ξ − 1 6 10 −5 . Finally, regarding massless particles far from the conformal coupling with gravity, the reheating temperature becomes which leads, using also Eq. (36), to whose solution is n s = 0.970, corresponding to r = 0.0018, which now clearly falls within the bounds imposed in our model so as to satisfy the observational data, but again fails to fulfill the bound T R ≤ 1 GeV suggested in the introduction. Finally we are going to briefly study the dynamics of the equationφ + 3H (ϕ,φ)φ + V ϕ = 0. In Fig. 6 we have represented in blue some orbits which are solution of these dynamical system, observing that, as in the former case, they start with infinite energy at t → −∞, they experience a phase transition at ϕ = 0 and then evolve towards a de Sitter phase for t → ∞. Moreover, we have drawn in black the slow-roll We observe that there is a long period in inflation where orbits coincide with the slow-roll one. Furthermore, as happened with the analytical solution in the former case we studied, the slow-roll solution is a limit orbit that starts with a de Sitter phase at t → −∞, since all the other orbits start with infinite energy.

Higgs inflation (HI)
Another potential that could work would be the following Higgs inflation (HI) style potential in the Einstein frame [44]: λ being a dimensionless positive parameter. In this case, the slow-roll parameters are * = 2s 2 * and η * = −s * (1 − s * ), s * being the same as above. In this case, the running, power spectrum and number of e-folds are However, in this model now we have N ∼ 1 s * ∼ = 1 1−n s , which is bounded by 50 for the allowed values of the spectral index. Thus, the potential is not a viable quintessential inflation model because it leads to an insufficient number of e-folds.

Power-law inflation (PLI)
Now, we are going to study a power-law inflation (PLI) [45] potential, adapted to quintessence With the same procedure as in the former cases, we obtain * = 2n 2 M pl ϕ * 2 and η = 2n(2n − 1) 1−n s and N = n+1 1−n s − n 2 . Hence, we can easily verify that, so as to obtain a number of e-folds 63 ≤ N ≤ 73 with a spectral index n s = 0.968 ± 0.006, we need 0.65 < n < 1.35.
Therefore, taking n = 1 is a good choice in order to match our model with the observational results. Thus, if we consider the range of spectral index 0.9685 ≤ n s ≤ 0.973, we obtain the desired number of e-folds. Regarding the ratio of tensor to scalar perturbations r = 16 * = 4(1 − n s ), the constraint r ≤ 0.12 is verified for n s ≥ 0.97. So does the running, which becomes α s = − 8 resulting that −3.7 × 10 −3 ≤ α s ≤ −3.1 × 10 −3 for 0.97 ≤ n s ≤ 0.973. Finally, the power spectrum obeys the following expression: (46) and, thus, since P ∼ 2 × 10 −9 , we can determine λ, namely λ ∼ 10 −11 .
Remark We can see that the power-law potential for n = 2 is equivalent during the inflation, i.e. for ϕ −|ϕ E | to the double well inflation (DWI) [46] potential, namely Therefore, we have seen that the DWI is not a valid model that matches with the corresponding observational results. Now, we are going to proceed analogously as for the potential previously studied. Therefore, we approximate the Hubble constant at the transition point by H E ∼ H end = M pl √ λ ≈ 10 −4 (1 − n s )M pl . Now, combining Eq. (13) and the number of e-folds for this particular potential, we obtain where, as before, Y = 6 1−n s . Thus, as we can see in Fig. 7, the bounds coming from the observational values used throughout all this paper are satisfied for 1 MeV ≤ T R ≤ 10 4 GeV (corresponding to 0.970 ≤ n s ≤ 0.973).  Fig. 7 Evolution of the tensor/scalar ratio r (left) and the spectral index n s (right) versus the reheating temperature T R If we start considering the production of massless particles nearly conformally coupled to gravity, then the reheating temperature becomes Hence, for 1 MeV ≤ T R ≤ 10 4 GeV, we obtain 10 −7 ξ − 1 6 5 × 10 −3 . On the other hand, in the case of heavy massive particles (m ∼ M pl ), since the first derivative of the potential is continuous at the transition phase, one has to use, for the β-Bogoliubov coefficient the expression given in [39], To obtain the third derivative of the Hubble parameter first of all we use the formula ...
. Using the con-  [31] and, using the thermalization rate introduced in Sect. 2.3, one gets the following reheating temperature: Remark 3.1 Note that the formula of the reheating temperature is a little bit different from the one obtained in [39], since there another thermalization rate has been used.
By combining it with Eq. (48), namely one obtains that n s ∼ = 0.9721, r ∼ = 0.1117 and T R ∼ 11 MeV, which means that this potential supports the production of heavy massive particles. If we consider massless particles far from the conformal coupling, then the reheating temperature becomes Therefore, it falls out of the bounds that we have previously found (1 MeV ≤ T R ≤ 10 4 GeV) and, thus, our model does not support the production of massless particles far from the conformal coupling.
Regarding the dynamics of the equationφ +3H (ϕ,φ)φ + V ϕ = 0, we have also proceeded analogously as with the other potential. In Fig. 8 we have represented some blue orbits solution of the dynamical system which start with infinite energy at t → −∞ and then evolve to the de Sitter phase for t → ∞ after having gone through a phase transition at ϕ = 0. In this case, the slow-roll solution (in black) iṡ 3 . We note that in this case the slow-roll solution does not start with a de Sitter phase.
We can compare this power-law potential with the one studied in [38], namely Firstly we note that, for ϕ ϕ E , i.e. during all the inflation, the behavior will be analogous. Regarding the reheating temperature, the upper bound T R ≤ 10 4 GeV coincides. It only differs the fact that for heavy massive particles nearly conformally coupled with gravity the reheating temperature is in the GeV regime, instead of the MeV one, for the model coming from the potential in Eq. (54).

Other quintessential inflation potentials
In this last section we are going to study some other potentials that appear in [44], studying whether they can be implemented in the quintessential inflation.

Open string tachyonic inflation (OSTI)
We consider the following adapted form of the OSTI potential [47]: where |ϕ 0 | M pl . We note that this potential is only nonnegative for ϕ > −|ϕ 0 |. Therefore, we will always move through this domain. Now, let us compute the cosmological parameters.
Since cancels at ϕ = −e −1/2 |ϕ 0 |, the inflation period will happen between this value and ϕ end (corresponding to = 1). Using that |ϕ 0 | M pl , we obtain ϕ end ≈ √ 2M pl , ϕ * ≈ showing, as in Fig. 9, that the valid bounds of the reheating temperature for this potential are 1 MeV ≤ T R ≤ 10 5 GeV. As in the case of the power-law potential with n = 1, this potential has a continuous derivative at the transition phase. Thus, when considering the production of massless particles nearly conformally coupled to gravity, the reheating temperature is Thus, for 1 MeV ≤ T R ≤ 10 5 GeV, we obtain 10 −8 ξ − 1 6 3 × 10 −3 . However, dealing with heavy massive particles (m ∼ M pl ), one can see that the second derivative of the potential, and thus the third derivative of the Hubble parameter, diverges at the transition phase, which means that we cannot use the WKB solution to approximate the modes. Therefore, we are not able to compute, in this case, the reheating temperature. Regarding the case of massless particles far from the conformal coupling, the reheating temperature becomes which shows that this model neither supports the production of these particles. Finally, as done with the other potentials, we study the dynamics of the equationφ+3H (ϕ,φ)φ+V ϕ = 0. In Fig. 10 we see again that the blue orbits (solution of the dynamics equation) start with infinite energy at t → −∞ and then evolve to the de Sitter phase after a phase transition at ϕ = 0. As happened with the PLI potential, the slow-roll inflation,

Witten-O'Raifeartaigh inflation (WRI)
In this case, the version of the WRI [48,49] potential is where |ϕ E | M pl . The cosmological parameters are * = 2 With regards to the number of e-folds, it can be exactly integrated as N = |ϕ E | 2 8M 2 pl x 2 (2 ln(x) − 1) x end x * (where x = −ϕ |ϕ E | ), whose expression up to order 2 in which results in a too high number of e-folds, namely N 100.

Kähler moduli inflation I (KMII)
The expression of the Kähler moduli inflation I (KMII) potential [44] is where α is a positive dimensionless constant such that α ≥ e and ϕ E is the value of ϕ where ϕ M pl e −ϕ/M pl = 1/α such that ϕ E ≥ M pl . We note that, in contrast to the former potentials considered, in this one we are going to assume thatφ < 0. The cosmological parameters are where x = ϕ/M pl . By considering that x * 1, we obtain that n s − 1 ≈ 2αx * e −x * . Regarding the number of e-folds, they can be exactly integrated, namely where x end is the value of x where = 1 and E i is the exponential integral function, which verifies that, for x * 1, x * , this being the dominant term in the previous equation. Hence, by using x * e −x * ≈ n s −1 2α , one obtains which does not fulfill our bounds for the number of e-folds and the spectral index. Therefore, we have proved that this potential is not a viable quintessential inflation model.

Brane inflation (BI)
The brane inflation (BI) potential behaves as [50] V where μ and p are positive dimensionless parameters and ϕ E ≡ −μM pl . We will consider for simplicity only p ∈ N. Thus, the cosmological parameters are

Loop inflation (LI)
In this case the potential behaves as [51] V where λ and α are positive dimensionless constants. The slow-roll parameters are where we have introduced the parameter x ≡ ϕ M pl . We consider to different asymptotic cases: In this case one has * ∼ = and thus, n s −1 ∼ = −2η * ∼ = 2α . For the number of e-folds one has which leads, as in the case of HI, to a not high enough number of e-folds. 2. α 1: Now the slow-roll parameters become * = and the spectral index and the tensor/scalar ratio will be as a function of x * 1 − n s = 1 Then, at 2σ C.L., for the allowed values of the spectral index, we can see, after some numerics, that x * ranges in the domain 6.94 ≤ x * ≤ 7.98. On the other hand, the number of e-folds is Using the range of values for x * one finds that 34 N 50, which results from the nucleosynthesis bounds.

Conclusions
We have adapted some inflationary potentials to quintessence inflation, extending them to zero after they vanish and adding a small cosmological constant. Once we have done it, we have tested the models imposing the requirements that: (1) they fit well with the current observational data provided by BICEP and Planck teams. (2) The number of e-folds must range between 63 and 73, this number is larger than the usual one used for potentials with a deep well, due to the kination phase after inflation. (3) The reheating temperature due to the gravitation particle production during the phase transition from inflation to kination has to be compatible with the nucleosynthesis success, i.e., it has to range between 1 MeV and 10 9 GeV, although if one wants to remove the gravitino and moduli problems, which appear in supergravity, this temperature has to be less than 1 GeV. Our study shows that the potentials WRI and KMII lead to a too high number of e-folds, while for HI and LI potentials this number is too small. Other potentials such as ESI, PLI (only when the potential is quadratic), OSTI and BI satisfy the prescriptions 1 and 2. Dealing with the reheating temperature, all the four models, namely ESI, PLI, OSTI and BI, lead to a temperature compatible with the nucleosynthesis bounds. However, since the ESI potential comes from a supersymmetric gravitational one, the problems related with the gravitino and moduli overproduction will appear, which can only be removed when the reheating temperature is in the MeV regime. Then, for this potential, when reheating is due to the production of massless nearly conformally coupled particles (the coupling constant is very close to 1/6), these problems are removed because the reheating temperature is very low (in the MeV regime). On the contrary, reheating due to the production of very heavy massive particles conformally coupled with gravity or massless particles far from the conformal coupling leads to a too high reheating temperature, which does not avoid the gravitino and moduli problems for the ESI potential.
Fortunately, the other three potentials, namely PLI, OSTI and BI, do not come from supersymmetric potentials and, thus, do not suffer these problems because gravitino and moduli fields only appear in supergravity, meaning that the reheating via very massive particles conformally coupled with gravity or via massless particles far from the conformal coupling is also viable for these three potentials. and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. Funded by SCOAP 3 .

Appendix: A critical review of the gravitational wave particle production
Gravitons satisfy the same equation as scalar fields minimally coupled with gravity, so its energy density at the end of inflation is twice that of a single scalar field due to the two graviton polarization states. Experimental observations show that, at the end of inflation, the ratio of the energy density of the relic gravitational waves to the one of the produced particles is less than 10 −2 [16]. In [37], using the results of [40] and [41], the authors take as the energy density of the where H E is the value of the Hubble parameter at the transition time, so assuming that there are N s different light scalar fields with a coupling constant near zero, one deduces that ρ g ρ χ = 2 N s , where ρ χ denotes the energy density of the produced particles. Then, according to the observational bound, this means that the number of different scalar fields must be greater than 100, which only happens in some supersymmetric theories [37].
Note that in this situation the reheating temperature is of the order 10 3 N

3/4 s
GeV ∼ 3 × 10 4 GeV, allowing the overproduction of gravitinos and moduli fields. The problem is even worse when one considers reheating via the production of heavy massive particles because this mechanism is less efficient than the one due to the creation of massless particles. The situation could be addressed when one considers models whose potential does not come from supergravity, which avoids the gravitino and moduli problems, and reheating is produced via oscillations of an auxiliary field in hibrid quintessential models [52] or via instant preheating [15][16][17][18][19] because these kinds of reheating are more efficient than the one due to gravitational particle production.
However, a more detailed investigation shows that the results used in [37] about the production of gravitons may not be correct and the conclusions could be wrong. Here, we summarize the problems that appear in this calculation 1. Gravitational particle production of massless particles far from the conformal coupling is only considered in toy models, where there is a transition from the exact de Sitter phase to a kination or radiation one. This happens because the vacuum modes are represented by Hankel functions and one can perform analytic calculation, but, as is well known, the inflationary period is not exactly de Sitter. 2. In the calculations, only the energy density of the particles is considered and vacuum polarization effects (see [53] for a detailed calculation) are disregarded, but as has been showed in [13] they could be very important in some inflationary scenarios leading to undesirable effects. 3. Performing the calculations, an UV divergence appears but it is removed by disregarding the modes inside the Hubble horizon, because they "do not feel" gravity and, thus, do not produce particles, but in the case of gravitons we obtain an IR divergence which is ignored and produces important consequences. Let us explain it in detail: When one considers a transition from the de Sitter phase to radiation, in [40] it was showed that the density energy of the produced massless particles coupled with gravity, in the long wave-length limit, i.e., for modes that leave the Hubble radius after the phase transition (k ≤ a E H E , H E being the Hubble parameter during the de Sitter phase and a E the scale factor at the transition time), satisfies where ν 2 = 9 4 − 12ξ , ξ being the coupling constant. This quantity is IR divergent for ν ≥ 3/2, that is, for ξ ≤ 0, and in this case, to obtain the energy density of the produced particles, one cannot simply integrate in the domain [0, a E H E ]. One can only do that when ξ > 0, obtaining the well-known result ρ ∼ 10 −2 H 4 E a E a(t) 4 . Thus, the case of gravitons has to be considered in another way, because for them ξ = 0. This was done in [54], where the IR singularity is removed assuming an early radiation phase before inflation (see [55]). Then, let a i be the scale factor when the universe passes from the early radiation phase to inflation. In this situation, if one is only interested in the gravitons produced by the modes that leave the Hubble radius during the de Sitter phase, i.e., in modes satisfying a i H E ≤ k ≤ a E H E , one obtains [54] where N ≡ ln a E a i is the number of e-folds that the de Sitter phase lasts. On the other hand, dealing with massless particles with a small negative coupling constant, one has and thus where we have used the fact that for negative small values of the coupling constant ν = 3 2 +4|ξ |. Since N has to be greater than N (the number of e-folds from the leaving of the pivot scale to the end of inflation), one can safely take N ≥ 100, then choosing for example N = 100, and ξ = −0.012 the bound is reached. So, a unique massless quantum field with a small negative coupling constant satisfies the bound ρ g ρ χ ≤ 10 −2 , meaning that gravitational supersymmetric potentials are not needed and the gravitino and moduli problems do not appear, allowing a greater reheating temperature. In this case, the reheating temperature, taking as usual to obey H E ∼ 10 −6 M pl , is given by T R ∼ 2 × 10 4 e 8|ξ |N −1 |ξ |

3/4
GeV, leading for the values N = 100 and ξ = −0.012 to the reheating temperature T R ∼ 7 × 10 8 GeV, which enters in the viable range.
Summing up, the problem of the production of massless particles far from the conformal coupling and in particular the overproduction of gravitational waves in quintessential inflation is a problem that deserves future investigation for many reasons, but basically because nowadays we do not have any analytic method which allows us to calculate, in realistic models, the backreaction of a massless quantum field not conformally coupled with gravity, and only results that are valid for toy models which may not apply to viable models are used to perform analytic calculations.