Hamilton-Jacobi approach for quasi-exponential inflation: predictions and constraints after Planck 2015 results

In the present work we study the consequences of considering an inflationary universe model in which the Hubble rate has a quasi-exponential dependence in the inflaton field, given by $H(\phi)=H_{inf}\exp \left[\frac{\frac{\phi}{m_p}}{p\left(1+\frac{\phi}{m_p}\right)}\right]$. We analyze the inflation dynamics under the Hamilton-Jacobi approach, which allows us to consider $H(\phi)$, rather than $V(\phi)$, as the fundamental quantity to be specified. By comparing the theoretical predictions of the model together with the allowed contour plots in the $n_s-r$ plane and the amplitude of primordial scalar perturbations from the latest Planck data, the parameters charactering this model are constrained. The model predicts values for the tensor-to-scalar ratio $r$ and for the running of the scalar spectral index $dn_s/ d\ln k$ consistent with the current bounds imposed by Planck, and we conclude that the model is viable.


I. INTRODUCTION
Inflation has become the most acceptable paradigm that describes the physics of the very early universe. Besides of solving most of the shortcomings of the hot big-bang scenario, like the horizon, the flatness, and the monopole problems [1][2][3][4][5][6], inflation also generates a causal mechanism to explain the large-scale structure (LSS) of the universe [7][8][9][10][11] and the origin of the anisotropies observed in the cosmic microwave background (CMB) radiation [12][13][14][15][16][17][18][19], since primordial density perturbations may be sourced from quantum fluctuations of the inflaton scalar field during the inflationary expansion.
Several representative inflationary models have been studied within the framework of the so-called slow-roll approximation [20], where the kinetic term of the inflaton field is much smaller than the potential energy, i.e.φ 2 V (φ), together with the approximation φ H φ . Moreover, in this approach the full shape of the inflaton potential is considered in order to identify the value of the scalar field at the end of inflation and hence the value of the scalar field when the largest scales observable today cross the Hubble radius.
Upon comparison to the current cosmological and astronomical observations, specially those related with the CMB temperature anisotropies, it is possible to constrain several inflation models. Particularly, the constraints in the n s − r plane give us the predictions of a number of representative inflationary potentials. Recently, the Planck collaboration has published new data of enhanced precision of the CMB anisotropies [19]. Here, the Planck full mission data has improved the upper bound on the tensor-to-scalar ratio r 0.002 < 0.11(95% CL) which is similar to obtained from [17], in which r < 0.12 (95% CL). From the particle physics point of view, it is natural to begin by specifying the functional form of the potential. However, even for simple choices, such as exponential [21], constant [1] or power-law potentials [6], it is not possible to go further analytically. An alternative way is to specify the time-dependence of the scale factor a(t). Following Refs. [22][23][24]. exact solutions can also be found in the scenario of intermediate inflation. In this inflationary model the scale factor evolves as a(t) ∼ exp At f , where A and f are two constant parameters such that A > 0 and 0 < f < 1. The expansion rate of this scale factor is slower than de Sitter inflation [1], for which a(t) ∼ exp(Ht), where H is the Hubble rate, which is a contant, but faster than power-law inflation, a(t) ∼ t n [21] , where n > 1.
Alternative to the slow-roll approximation, there is another method for studying inflation known as the Hamilton-Jacobi approach [25,26]. This formulation is a powerful way of rewriting the equations of motion for single-field inflation. It can be derived by considering the scalar field itself to be the time variable, which is possible during any epoch in which the scalar field evolves monotonically with time. It allows us to consider the Hubble rate or Hubble function H(φ) (from now on, not confuse with with the Hamiltonian function H), rather than the inflaton scalar potential V (φ), as the fundamental quantity to be specified.
Because H(φ), unlike V (φ), is a geometric quantity, inflation is described more naturally in that language. The advantage of such an approach is that the form of the potential is readily deduced. As it was suggested in Refs. [27][28][29], H(φ) should be viewed as the solution generating function when analysing inflationary cosmologies. For instance, H(φ) ∼ exp(φ) gives the power-law inflation model [21]. Furthermore, this formalism has been considered by Planck collaboration in order to reconstruct the inflaton potential beyond slow-roll approximation [17,19]. For a representative list of recent inflation models studied under Hamilton-Jacobi formalism where several expressions for H(φ) have been considered, see Refs. [30][31][32][33][34][35].
Following Ref. [31], a phenomenological quasi-exponential Hubble rate H(φ) yielding an inflationary solution was proposed to be where p is a dimensionless parameter, m p denotes the Planck mass, and H inf is a parameter with dimensions of Planck mass. It is interesting to mention that this model presents an improvement in comparison to power-law inflation model [21], because the first one addresses the graceful-exit problem of inflation and the value predicted for the tensor-to-scalar ratio was compatible with Seven-Year WMAP [14], being supported by the current data available at that time.
The main goal of the present work is to study the realization of inflation by reconsidering the expression for the Hubble rate given by Eq.(1), in the light of the recent Planck results.
We stress that our work is different to previous work [31] in three ways. Firstly, in this work we restrict ourselves only to the inflationary predictions of this models. Secondly, in the previous paper the authors did not used the contour plots in the n s − r and n s − dn s /d ln k planes to constrain the parameters of the model they studied. Finally, in our work here we make use the latest data from Planck, not available at that time, to put bounds on the parameters of the model. We will show that our results are modified compared to [31] using the Planck results. By comparing the theoretical predictions of the model together with the allowed contour plots in the n s − r plane, the model predicts a value for the tensor-to-scalar ratio r detectable by Planck, and we conclude that the model is viable.
We organize our work as follows: After this introduction, in the next section we summarize the dynamics of inflation in the Hamilton-Jacobi formalism. In the third section we analyze the inflation dynamics of the Hubble rate given by Eq.(1) in the Hamilton-Jacobi framework, obtaining expressions for the scalar power spectrum, scalar spectral index, and tensor-toscalar ratio in terms of the free parameters characterizing the model which are constrained by considering the Planck 2015 results, through the allowed contour plots in the r − n s plane and the amplitude of the scalar power spectrum. In section IV we discuss a little further how in the warm inflation scenario the radiation-dominated phase is achieved without introducing the reheating phase for this quasi-exponential Hubble function. In the last section we finish with our conclusions. We choose units so that c = = 1.

A. Dynamics of inflation
In the simplest model of inflation in Einstein's General Relativity is a classical homogeneous scalar field φ = φ(t), named the inflaton field, which is introduced into the action.
The properties of the scalar potential determine how inflation evolves. For a flat Friedman-Lemaître-Robertson-Walker (FLRW) metric, the Friedmann and acceleration equations be- respectively, where m p = 1/G corresponds to the Planck mass.
Besides the Einstein equations, the field satisfies the Klein-Gordon equation in this FLRW universeφ where prime indicates derivative with respect to φ, and dot a derivative with respect to cosmic time.
The Friedmann and Klein-Gordon equations are the basis to construct the Hamilton-Jacobi formulation. By combining Eqs. (2) and (4), we obtain the following expressioṅ which gives the relation between φ and cosmic time t. This allows us to write the Friedmann equation in a first-order form, from which the inflaton potential V (φ) becomes This last equation is the Hamilton-Jacobi equation [20]. It allows us to consider H(φ), rather than V (φ), as the fundamental quantity to be specified. On the other hand, a relatiom da dφ = a Ḣ φ with Eq.(5) yields a differential equation for a(φ), whose integration becomes This equation implies that, once the functional form of a geometrical quantity H(φ) has been specified, the cosmological dynamics is determined. The advantage of such an approach is that the form of the potential is readily deduced from Eq.(6).
We can use the Hamilton-Jacobi formalism to write down a slightly different version of the slow-roll approximation, defining the Hubble hierarchy parameters H and η H as [20] In the slow-roll limit, H → and η H → η − , where and η are the usual slow-roll parameters. By using Eq.(8), the acceleration equation (3) is rewritten as During inflation H satisfies the condition H < 1, and the inflationary expansion ends when H becomes one.
On the other hand, the number of e-folds between the Hubble-radius crossing and the end of inflation yields where φ * and φ end are the values of the scalar field when the cosmological scales cross the Hubble-radius and at the end of inflation, respectively. The last value is found by

B. Attractor behavior
The Hamilton-Jacobi approach is usefull to show that all possible inflationary trajectories will rapidly converge to a common attractor solution, if they are sufficiently close to each other initially. This is exactly the behaviour that one expect within the slow-roll approximation, but the proof do not use of that approximation. Suppose that H 0 (φ) is any solution to Eq.(6), inflationary or not. If we add to this solution a linear homogeneous perturbation δH(φ), the attractor behaviour will be satisfied if δH(φ) H 0 (φ) tends quickly to zero as φ evolves [36]. Replacing H(φ) = H 0 (φ) + δH(φ) in Eq.(6) and linearizing, we have that Integrating last expression we get where δH(φ i ) is the initial value of the perturbation at φ = φ i . Knowing H(φ), it is possible to study the behaviour of perturbation δH(φ).
In the next section we will give a review of cosmological perturbations and use Hubble hierarchy parameters for describing scalar and tensor perturbations.

C. Cosmological perturbations
We consider the gauge invariant quantity ζ = −ψ − H δρ ρ . Here, ζ is defined on slices of uniform density and reduces to the curvature perturbation R at super-horizon scales. A fundamental feature of ζ is that it is nearly constant on super-horizon scales [37], and in fact this property does not depend on the gravitational field equations [38]. Therefore, at super-horizon scales we have that R = H δφ φ , where |δφ| = H/2π. In this way, the power spectrum of scalar perturbations is given by [20,39] This perturbation is evaluated at Hubble radius crossing k = aH during inflation.
Important observational quantities are not only the amplitude of the primordial curvature perturbations but also the scalar spectral index which represents the scale dependence of the power spectrum, defined by Thus, the scalar spectral index of the power spectrum (14) is given by where H and η H are the Hubble hierarchy parameters, given by Eqs. (8) and (9), respectively.
We also introduce the running of the scalar spectral index, which represents the scale dependence of the spectral index, by n run = dns d ln k , yielding where ξ 2 H is a third Hubble hierarchy parameter, defined by [40] On the other hand, the power spectrum of tensor perturbations generated from inflation is given by [20,39] As the cosmological parameter related to the primordial tensor perturbation, the ratio between the amplitude of the primordial tensor perturbation and that of the primordial curvature perturbation, the so-called tensor-to-scalar ratio, defined by r ≡ P R P T , becomes Additionally, by combining Eqs. (5) and (20), we obtain the Lyth bound [41], which relates the tensor-to-scalar ratio and the evolution of the scalar inflaton field This means that the tensor-to-scalar ratio measures (up to order-one constants) the distance that the inflaton field φ traveled in field space during inflation. For detectable r, this implies ∆φ ∼ m p .
Up to now, the basis of the Hamilton-Jacobi formalism has been presented. In next section, in order to get a specific result, we are going to introduce the quasi-exponential form for Hubble rate given by Eq. (1).

A. Dynamics of inflation
In this section we describe an inflationary model by using the quasi-exponential generating function H(φ) given by Eq.(1). By combining Eqs. (1) and (5) we get that where Ei [x(φ)] denotes the Exponential Integral function [42], given by the integral with . This latter expression yields the inflaton field as function of cosmic time. Additionally, the scale factor a(t) turns out to be where a i denotes the value of the scale factor when the inflaton field has the value φ i , i.e., a i = a(φ i ). In order to have an inflationary solution, the condition φ < φ i must be satisfied, which means that the inflaton starts to rolling down the potential at large values of φ i .
The form of the potential is readily deduced from Eqs. (1) and (6), which results to be given by For sake of comparison, in the slow-roll approximation,φ 2 V (φ) and φ H φ , the inflaton potential becomes and respectively.
From the condition H (φ end ) = 1, we obtain the value of the inflaton field at the end of the inflationary expansion, yielding Restricting ourselves only to positive incursion of the inflaton field trough the potential, the allowed range for p becomes 0 < p < 1 2 √ π ≈ 0.282. The number of inflationary e-folds between the values of the scalar field when a given perturbation scale leaves the Hubble-radius and at the end of inflation, can be computed from Eqs.(1), (11), and (29), resulting in By solving Eq.(30) for φ * , we may obtain the value of the scalar field at the time of Hubbleradius crossing, giving As we shall see later on, the several inflationary observables will be evaluated at the value of the inflaton field given by Eq. (31).

B. Attractor behavior
As final step of the analysis of background dynamic for this model, the attractor behavior of the solution is considered. From Eqs. (1) and (13), the solution for the perturbation δH(φ) In order to determine the attractor behaviour quantitatively, we consider the initial value of the inflaton field to be φ i = αφ * , with α > 1 and φ * given by Eq. (31). Then, replacing φ i into Eq.(32), the perturbation δH(φ) δH(φ i ) becomes For this quasi-exponential form of the Hubble function, the inflaton field decreases as time increases, therefore the exponential term on the right-hand side of Eq.(33) decreases by passing time and tends to zero rapidly, then the perturbation of the Hubble function vanishes, and the model has an attractive behavior.
Regarding the cosmological perturbations, the amplitude of the primordial curvature perturbation, using Eqs. (1) and (14), is found to be The scalar spectral index, using Eqs. (16), (27), and (28), becomes Additionally, the running of the scalar spectral index n run is found to be Finally, the tensor-to-scalar ratio can be obtained from Eqs. (20) and (27), yielding r = 1 After evaluating these inflationary observables at the value of the scalar field when a given perturbation scale leaves the Hubble-radius, given by (31), we may compare the theoretical predictions of our model with the observational data in order to obtain constraints on the parameters that characterize it.
The amplitude of the primordial curvature perturbation, the scalar spectral index, the running of the scalar spectral index, and the tensor-to-scalar ratio, evaluated at the Hubbleradius crossing k = aH, become The trajectories in the n s − r plane for the model studied here may be generated by plotting Eqs. (39) and (41) parametrically. In particular, we have obtained three different curves by fixing the number of e-folds to N = 55, 60, and 65, and plotting with respect to the parameter p in the range 0 < p < 1 2 √ π , obtained by considering a positive incursion of the inflaton field trough the potential, which gives an upper bound for p. The Fig.(2) shows the plot of the tensor-to-scalar ratio r versus the scalar spectral index n s for the quasi-exponential Hubble rate. Here, we have considered the two-dimensional marginalized joint confidence contours for (n s , r), at the 68% and 95% CL, from the latest Planck data [19]. We can determinate numerically from Eq.(41) that, by fixing N , the tensor-to-scalar  Table (I) summarizes the constraints obtained on p and H inf using the last data of Planck.
As we can see, using the latest Planck results, through the allowed contour plots in the n s −r plane and the amplitude of the scalar power spectrum, we were able to find the allowed range for p and H inf . Particularly, the allowed contour plots in the n s − r plane impose a strong constraint on the lower bound for p, which was not considered by the authors in the previous work [31].
After replacing Eq.(42) into Eq.(25), we can plot the scalar potential V as function of the inflaton field, as is shown in Fig.(3). We have plotted the inflaton potential for 3 different values for the number of e-folds N : the solid, dashed, and dotted lines correspond to N = 55, 60, and 65, respectively, and fixing the value p = 0.15, which lies in the allowed range for each value of N already obtained by using the Planck data. It is interesting to mention that this quasi-exponential form of the Hubble rate presents a graceful-exit of inflation, however, the inflaton potential does not present a minimum, which raises the issue of how to address the problem of reheating in this model. A way to address this problem may be to study this model in the warm inflation scenario [43][44][45], which has the attractiveness that it avoids the reheating period at the end of the accelerated expansion.
In such as scenario, the dissipative effects are important, and radiation production takes place at the same time as the expansion of the universe. When the universe heats up and becomes radiation dominated, inflation ends and the universe smoothly enters the radiation Big Bang phase. In section IV we discuss a little further how in the warm inflation scenario the radiation-dominated phase is achieved without introducing the reheating phase for this quasi-exponential Hubble function.
In order to determine the prediction of this model regarding the running of the spectral index, the trajectories in the n s − dn s /d ln k plane may be generated by plotting Eqs. (40) and (41) parametrically. In particular, we have obtained three different curves by fixing the number of e-folds to N = 55, 60, and 65, and plotting with respect to the parameter p in the allowed range obtained for each value of N , which are shown in left panel of Fig.(4). In order to compare the previous predictions with the observational data, the right panel of Fig.(4), shows the two-dimensional marginalized joint confidence contours for (n s , dn s /d ln k), at the confidence contours for (n s , dn s /d ln k), at the 68% and 95% CL, in the presence of a non-zero tensor contribution, from the latest Planck data [19] and the dns d ln k (n s ) curve for N = 60 (red-dashed line).
68% and 95% CL, in the presence of a non-zero tensor contribution, from the latest Planck data [19]. Given the indistinguishability of the curves in left panel, for the right panel we have only considered the curve corresponding to N = 60 (red-dashed line) to compare with Planck data. The thin black line in right panel shows the prediction of single-field monomial inflation models with 50 < N < 60. From both panels, we observe that all three dns d ln k (n s ) curves lie inside the 68% as well as 95% CL regions from Planck.
Finally, by replacing the expression that we have found for the tensor-to-scalar ratio r(N ), expressed by Eq.(41), into Eq. (21), the incursion of the inflaton field is found to be In particular, by considering the constraint on p for each value of N already obtained, we get that 2.623 < ∆φ mp < 3.428, 2.727 < ∆φ mp < 3.852, and 2.826 < ∆φ mp < 4.267, for N = 55, N = 60, and N = 70, respectively. We note that the incursion of the inflaton field decreases as p increases.

IV. A FIRST APPROACH TO WARM INFLATION WITH A QUASI-EXPONENTIAL HUBBLE FUNCTION
As we mentioned at previous section, the warm inflation scenario, as opposed to the standard cold inflation, has the attractive feature that it avoids the reheating period at the end of the accelerated expansion. During the evolution of warm inflation dissipative effects are important, and radiation production takes place at the same time as the expansion of the universe. The dissipative effects arise from a friction term Γ which accounts for the processes of the scalar field dissipating into a thermal bath. In addition, in warm inflationary scenario the density perturbations arise from thermal fluctuations of the inflaton and dominate over the quantum ones. In this form, an essential condition for warm inflation to occur is the existence of a radiation component with temperature T > H, since the thermal and quantum fluctuations are proportional to T and H, respectively [43][44][45]. When the universe heats up and becomes radiation dominated, inflation ends and the universe smoothly enters in the radiation Big-Bang.
We start by considering a spatially flat FLRW universe containing a self-interacting inflaton scalar field φ with energy density and pressure given by ρ φ =φ 2 /2 + V (φ) and P φ =φ 2 /2 − V (φ), respectively, and a radiation field with energy density ρ γ . The corresponding Friedmann equations reads and the dynamics of ρ φ and ρ γ is described by the equations [43][44][45] where the dissipative coefficient Γ > 0 produces the decay of the scalar field into radiation.
Recall that this decay rate can be assumed to be a function of the temperature of the thermal bath Γ(T ), or a function of the scalar field Γ(φ), or a function of Γ(T, φ) or simply a constant [43][44][45].
During warm inflation, the energy density related to the scalar field predominates over the energy density of the radiation field, i.e., ρ φ ρ γ [43][44][45], but even if small when compared to the inflaton energy density it can be larger than the expansion rate with ρ 1/4 γ > H.
Assuming thermalization, this translates roughly into T > H, which is the condition for warm inflation to occur.
When H, φ, and Γ are slowly varying, which is a good approximation during inflation, the production of radiation becomes quasi-stable, i.e.,ρ γ 4Hρ γ andρ γ Γφ 2 , see Refs. [43][44][45]. Then, the energy density of the radiation field becomes If we consider thermalization, then the energy density of the radiation field could be written as ρ γ = C γ T 4 , where the constant C γ = π 2 g * /30. Here, g * represents the number of relativistic degrees of freedom. By combining Eqs. (44) and (45), the time derivative of the Hubble function is given byḢ where Q is the dissipative ratio, defined as Q ≡ Γ/3H. In warm inflation, we can distinguish between two possible scenarios, namely the weak and strong dissipative regimes, defined as Q 1 and Q 1, respectively. In the weak dissipative regime, the Hubble damping is still the dominant term, however, in the strong dissipative regime, the dissipative coefficient Γ controls the damped evolution of the inflaton field.
By expressing the time derivative of the Hubble function in terms of inflaton field derivative, the time derivative of inflaton field becomeṡ which is the same expression found in [46], where warm inflation under the Hamilton-Jacobi formalism has been studied recently.
Introducing the dimensionless Hubble hierarchy parameter H , we write It is possible to find a relation between the energy densities ρ γ and ρ φ by combining Eqs.(47), (49), and (50), so that Warm inflation takes place when the parameter H satisfies H < 1. This condition given above implies that during inflation the energy density of the inflaton field satisfies ρ φ > 2(1+Q) Q ρ γ . Then, at the end of inflation, when H = 1, we have that ρ γ = Q 2(1+Q) ρ φ . The universe stops inflating and heats up to become radiation dominated at the time when ρ γ = ρ φ . This is one of the most attractive features of warm inflation, since provides a smooth transition to the radiation-dominated epoch without introducing a reheating epoch.
Given that for the quasi-exponential Hubble function the inflaton potential does not present a minimum, the dynamics for this model in warm inflation scenario provides a solution for the problem of reheating. The perturbation dynamics of warm inflation with the quasiexponential Hubble function deserves a more further analysis which goes beyond the scope of this work.

V. CONCLUSIONS
To summarize, in this article we have studied an inflationary model the Hubble rate has an quasi-exponential dependence in the inflaton field. We have studied the inflation dynamics in the Hamilton-Jacobi formalism, in which the scalar field itself to be the time variable, which is possible during any epoch in which the scalar field evolves monotonically with time. It allows us to consider the Hubble function H(φ), rather than the inflaton potential V (φ), as the fundamental quantity to be specified. Because H(φ), unlike V (φ), is a geometric quantity, inflation is described more naturally in that language. This model is characterized by the dimensionless parameter p and H inf . In order to constraints our model, we have considered the amplitude of the primordial scalar perturbations, as well as the allowed contour plots in n s − r and n s − dn s /d ln k planes from Planck 2015 data. First, in the n s − r plane we show the theoretical predictions of the model for three different values of e-folds N = 55, 60, and 65. By finding the points where each r(n s ) curve enters the joint 95% CL region, the allowed range for the p parameter may be determined. After that, using the constraint for the amplitude of scalar perturbations we determined the allowed range for H inf . In addition, in the n s − dn s /d ln k plane we show that all the three dns d ln k (n s ) curves for N = 55, 60, and 65 lie inside the 68% as well as 95% CL regions from Planck. This model predicts values for the tensor-to-scalar ratio r and for the running of the scalar spectral index consistent with the current bounds imposed by Planck 2015, and we conclude that the model is viable. As we mention before, this quasi-exponential form of the Hubble rate presents a graceful-exit of inflation, however, the inflaton potential does not present a minimum, which raises the issue of how to address the problem of reheating. In order to address last problem and as a first approach to further research, in section IV we discussed how in the warm inflation scenario the radiation-dominated phase is achieved without introducing the reheating phase for this quasi-exponential Hubble function. We hope to return to this point in the near future.