Warm $\frac{\lambda}{4}\phi^{4}$ inflationary universe model in light of Planck 2015 results

In the present work we show that warm chaotic inflation characterized by a simple $\frac{\lambda}{4}\phi^{4}$ self-interaction potential for the inflaton, excluded by current data in standard cold inflation, and by an inflaton decay rate proportional to the temperature, is in agreement with the latest Planck data. The parameters of the model are constrained, and our results show that the model predicts a negligible tensor-to-scalar ratio in the strong dissipative regime, while in the weak dissipative regime the tensor-to-scalar ratio can be large enough to be observed.


I. INTRODUCTION
The inflationary universe has become one of the central paradigms in modern cosmology. This is due to the fact that many long-standing problems of the Big Bang model, such as the horizon, flatness, homogeneity and monopole problems, find a natural explanation in the framework of the inflationary universe [1][2][3][4][5][6]. However, the essential feature of inflation is that it generates a mechanism to explain the Large-Scale Structure (LSS) of the universe [7][8][9][10][11] and provides a causal interpretation of the origin of the anisotropies observed in the Cosmic Microwave Background (CMB) radiation [12][13][14][15][16][17], since primordial density perturbations may be produced from quantum fluctuations during the inflationary era.
The original "old inflation" scenario assumed the inflaton was trapped in a metastable false vacuum and had to exit to the true vacuum via a first-order transition [1,2]. However, the exit could occur neither gracefully nor completely. The revised version of inflation was proposed by A. Linde [3,4], and A. Albrecht and J. Steinhardt [5] in 1982 referred as "new inflation". However, these scenarios suffer from theoretical problems about the duration of inflation and initial conditions. In 1983, A. Linde considered the case that the initial conditions for scalar field driving inflation may be chaotic, which is called "chaotic inflation" [6]. This inflation model can solve the remaining problems, where the potential was chosen to be cuadratic or quartic form, i.e. m 2 2 φ 2 or λ 4 φ 4 , terms that are always present in the scalar potential of the Higgs sector in all renormalizable gauge field theories [19] in which the gauge symmetry is spontaneously broken via the Englert-Brout-Higgs mechanism [20].
Such models are interesting for their simplicity, and has become one of the most favored, because they predict a significant amount of tensor perturbations due to the inflaton field gets across the trans-Planckian distance during inflation [21]. After that, many kinds of inflationary scenarios have been proposed, related to supersymmetry (SUSY) theory, brane world, string theory, etc. (for review, see [22][23][24][25] ).
On the other hand, with respect to the dynamical mechanisms of inflation, 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 [26]. 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 further relation to these dissipative effects, the dissipative coefficient Γ is a fundamental quantity, which has been computed from first principles in the context of supersymmetry. In particular, in Ref. [27], a supersymmetric model containing three superfields Φ, X, and Y has been studied, with a superpotential W = g where the scalar components of the superfields are φ, χ, and y respectively. For a scalar field with multiplets of heavy and light fields, and different decay mechanisms, it is possible to obtain several expressions for the dissipative coefficient Γ, see e.g., [27][28][29][30][31][32].
Following Refs. [30,31], a general parametrization of the dissipative coefficient Γ(T, φ) can be written as where the parameter a is related with the dissipative microscopic dynamics and the exponent m is an integer. This expression for the dissipative coefficient includes different cases studied in the literature, depending of the values of m (see Refs. [30,31]). Specifically, for the value m = 3, i.e., Γ ∝ T 3 /φ 2 , the parameter a corresponds to 0.02h 2 N Y , where a generic supersymmetric model with chiral superfields Φ, X, and Y i , i = 1, ...N Y has been considered.
This case corresponds to a low temperature regime, when the mass of the catalyst field m χ is larger than the temperature T [32]. On the other hand, m = 1, i.e., Γ ∝ T corresponds to a high temperature regime, where the thermal corrections to the catalyst field mass start to be important, where a = 0.97g 2 /h 2 [27]. For m = 0, the dissipative coefficient represents an exponentially decaying propagator in the high temperature regime. Finally, For m = −1, i.e., Γ ∝ φ 2 /T , agrees with the non-SUSY case [28,40]. Additionally, thermal fluctuations during the inflationary scenario may play a fundamental role in producing the primordial fluctuations [34,35]. During the 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 [26,34,35]. When the universe heats up and becomes radiation dominated, inflation ends and the universe smoothly enters in the radiation Big-Bang phase [26].
For a comprehensive review of warm inflation, see Ref. [36].
Upon comparison to the current cosmological and astronomical observations, specially those related with the CMB temperature anisotropies, it is possible to constrain the inflationary models. In particular, 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 [18] . 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). In particular, the λ 4 φ 4 model, which predicts a large value of the tensor-to-scalar ratio r, lies well outside of the joint 99.7% CL region in the n s − r, so it is ruled out by the data. This result confirms previous findings from e.g., Hinshaw et al. [15] in which this model is well outside the 95% CL for the WMAP 9-year data and is further excluded by CMB data at smaller scales.
In this way, the goal of the present work is to study the possibility that the λ 4 φ 4 model can be rescued in the warm inflation scenario and be able to agree with the latest observational data. In order to achieve this, we consider an inflaton decay rate Γ proportional to the temperature, which has been computed in the context of a high temperature supersymmetric model [27]. We stress that, in previous works (see Ref. [32]), the authors have also studied the quartic potential in the framework of warm inflation. However, our work is different in two ways. First, contrary to the standard cold inflation where the dynamics is determined only by the inflaton potential, in warm inflation also the dissipative coefficient plays an important role, and here we have considered an expression for it not studied in the previous works. Furthermore, in none of these papers the authors used the contour plots in the r and n s plane to constrain the parameters of the model they studied. On the contrary, in our work here we have used the latest data from Planck, not available at that time, to put bounds on the parameters of the model we have considered.
The outline of the paper is as follows: The next section presents a short review of the basics of warm inflation scenario. In Sect. III we study the dynamics of warm inflation for our quartic potential, in the weak and strong dissipative regimes; specifically, we obtain analytical expressions for the slow-roll parameters and the dissipative coefficient. Immediately, we compute the cosmological perturbations in both dissipative regimes, obtaining expressions for the inflationary observables such as the scalar power spectrum, the scalar spectral index, and the tensor-to-scalar ratio. Finally, Sect.IV summarizes our results and exhibits our conclusions. We choose units so that c = = 1.

A. Background evolution
We start by considering a spatially flat Friedmann-Robertson-Walker (FRW) 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 where M p = 1 √ 8πG is the reduced Planck mass. The dynamics of ρ φ and ρ γ is described by the equations [26] andρ 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 [26].
During warm inflation, the energy density related to the scalar field predominates over the energy density of the radiation field, i.e., ρ φ ρ γ [26,34,[37][38][39], but even if small when compared to the inflaton energy density it can be larger than the expansion rate with 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. [26,34,[37][38][39]. Then, the equations of motion reduce to where , φ denotes differentiation with respect to inflaton, and where R is the dissipative ratio defined as In warm inflation, we can distinguish between two possible scenarios, namely the weak and strong dissipative regimes, defined as R 1 and R 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.
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. In the Minimal Supersymmetric Standard Model (MSSM), g = 228.75 and C γ 70 [34]. Combining Eqs. (5) and (6) with ρ γ ∝ T 4 , the temperature of the thermal bath becomes On the other hand, the consistency conditions for the approximations to hold imply that a set of slow-roll conditions must be satisfied for a prolonged period of inflation to take place. For warm inflation, the slow-roll parameters are [27,34] = M 2 The slow-roll conditions for warm inflation can be expressed as [27,34] 1 + R, η 1 + R, β 1 + R, σ 1 + R When one these conditions is not longer satisfied, either the motion of the inflaton is no longer overdamped and slow-roll ends, or the radiation becomes comparable to the inflaton energy density. In this way, inflation ends when one of these parameters become the order of 1 + R.
From first principles in quantum field theory, the dissipative coefficient Γ has been computed. As we have seen in the introduction, the parametrization given by Eq. ..N Y has been considered. In particular, this inflation ratio decay has been studied extensively in the literature [32,33], including the quartic potential [36] . For the special case m = 1, the dissipative coefficient Γ ∝ T is related with the high temperature supersymmetry (SUSY) case [27]. Finally, for the cases m = 0 and m = −1, Γ represents an exponentially decaying propagator in the high temperature SUSY model and the non-SUSY case, respectively [28,40].

B. Perturbations
In the warm inflation scenario, a thermalized radiation component is present with T > H, then the inflaton fluctuations δφ are predominantly thermal instead quantum. In this way, following [34][35][36], the amplitude of the power spectrum of the curvature perturbation is given by where the normalization has been chosen in order to recover the standard cold inflation result when R → 0 and T H.
By the other hand, the scalar spectral index n s is given by [34] n s = 1 + dP R d ln k where A R = R 1+7R . Regarding to tensor perturbations, these do not couple to the thermal background, so gravitational waves are only generated by quantum fluctuations, as in standard inflation [41]. However, the tensor-to-scalar ratio r is modified with respect to standard cold inflation, yielding [36] r H T We can see that warm inflation predicts a tensor-to-scalar ratio suppressed by a factor (T /H)(1 + R) 5/2 > 1 compared with standard cold inflation.
When a specific form of the scalar potential and the dissipative coefficient Γ are considered, it is possible to study the background evolution under the slow-roll regime and the primordial perturbations in order to test the viability of warm inflation. In the following we will study how an inflaton decay rate proportional to the temperature, corresponding to the case m = 1, influences the inflationary dynamics for the quartic potential. We will restrict ourselves to the weak and strong dissipation regimes.

III. DYNAMICS OF WARM λ 4 φ 4 INFLATION
Although inflation is widely accepted as the standard paradigm for the early universe, it is not a theory yet as we don't know how to answer the question that naturally arises, "what is the inflaton and what is its potential?". After the recent discovery of the Higgs boson at CERN [42], which showed that elementary scalars exist in nature, the most natural and simplest thing to assume is that inflation is driven by the Higgs boson (in the standard model or in some extension of it). Unfortunately it is well known that the quartic potential, which is the simplest Higgs potential provided by particle physics in renormalizable theories, has been excluded by current data [43] since it predicts too many gravity waves. Although the presence of a non-minimal coupling can make the quartic potential viable [44], warm inflation provides another solution that is simpler and at the same time, as we have already mentioned, avoids the discussion about reheating. If we look at the expressions for the observables in the framework of warm inflation, we see that the key ingredient that can in principle reduce the tensor-to-scalar ratio, and bring the predictions of the model inside the region allowed by observational data, is the suppression factors (T /H) and R 5/2 . And this is exactly what happens indeed as we will show in the discussion to follow.
Warm inflation with a quartic potential for the inflaton has also been studied in [31,45]. However there are some differences, as in these works the authors have used another expression for the dissipative coefficient, they have not derived the allowed range for the parameters of the model they studied, and finally in our work we have used the most recent data available today.
and the Hubble parameter is given by In this way, for the weak regime, the slow-roll parameters become It is easy to see that the end of inflation is determined by the condition η = 1, where the scalar field takes the value φ end = 2 √ 3M p .
By the other hand, the number of e-folds is given by the standard formula where we have assumed that φ * φ end .
In the following, we will study the scalar and tensor perturbations. In the weak dissipative regime, the amplitude of the power spectrum (11) becomes By using Eqs. (14), (15), and (17), it may we written in terms in the number of e-folds as The power spectrum constraint P 1/2 R ∼ 10 −5 [17,18] determines the dimensionless coupling λ in terms of a and N , while the scalar spectral index (12) turns out to be which may be expressed in terms of the number of the e-folds, obtaining while the tensor-to-scalar ratio (13) becomes so eventually we can obtain r as a function of n s . Using Eqs. (19), (21), and (22), the relation r(n s ) is given by In figure 1, the relation r(n s ) is shown for several values of a. In the same plot we also show the curve for standard inflation (a = 0) as well as the contours allowed by the Planck latest data. When a decreases the curve is shifted upwards and finally lies outside the allowed contours. This induces a lower bound on a. On the other hand, when a increases the curve is shifted downwards, but R also increases and eventually the condition for being in the weak dissipative regime is violated. This induces an upper bound on a, which is found to be 6.5 × 10 −5 < a < 3.4 × 10 −2 . This implies that the Eq. (19), evaluated when the cosmological scales cross the Hubble horizon during inflation at 60 e-folds, gives us the constraint on λ determined by 10 −15 < λ < 10 −13 . It is interesting to note that this result is in agreement with the value obtained for λ in the standard cold inflation using the COBE normalization [46], given by λ ∼ 10 −14 .
FIG. 1: Plot of the tensor-to-scalar ratio r versus the scalar spectral index n s in the weak dissipative regime, for the quartic potential and an inflaton ratio decay Γ = aT . 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 [18]. In this plot we have used four different values of the parameter a, where the value a = 0 corresponds to standard cold inflation.

B. The strong dissipative regime
Considering our model evolves in agreement with the strong dissipative regime, where R 1, under the slow-roll approximation, the temperature of the radiation field becomes and the Hubble parameter is given by Eq. (15) In this way, for the strong regime, the slow-roll parameters become For the strong regime, inflation ends when one of these slow-roll parameters becomes the order of R. In this case, the end of inflation is determined by the condition η = R, where the inflaton takes the value φ end = ( By the other hand, the number of e-folds is given by where we have assumed that φ * φ end . Now, the amplitude of the power spectrum (11) becomes Similarly to weak regime, the amplitude of the power spectrum may we written in terms of the number of e-folds. Using Eqs. (15), (24), and (26), we have that In this case, the power spectrum does not depend on a and then the constraint P 1/2 R ∼ 10 −5 [17,18] determines the inflaton self-interaction coupling λ.
For this regime, the scalar spectral index (12) turns out to be which expressed in terms of the number of the e-folds yields Finally, for the tensor-to-scalar ratio (13) we have that which may be expressed as function of n s . Using Eqs. (28), (30), and (31), the relation r(n s ) is given by r 8.5 × 10 −9 π 10/3 a 4 (1 − n s ).
In figure 2, the relation r(n s ) is shown for two different values of a. In the same plot, as in the weak regime, we also show the curve for standard cold inflation (a = 0) as well as the contours allowed by the latest Planck data. When a decreases the curve is shifted upwards, but R also decreases and eventually the condition for being in the strong dissipative regime is violated. This induces a lower bound on a. By the other hand, when a increases the curve is shifted downwards, but R also increases and the condition for being in the strong dissipative regime is always satisfied. This implies that there is only a lower bound for a found by the requirement of staying in the strong dissipative regime, and given by a > 3.4×10 −2 . Finally, the Eq.(28), evaluated at 60 e-folds, gives us the constraint on λ, determined by λ ∼ 10 −15 .
This value is almost the same order that obtained for λ in the standard cold inflation.

FIG. 2:
Plot of the tensor-to-scalar ratio r versus the scalar spectral index n s in the strong dissipative regime, for the quartic potential and an inflaton ratio decay Γ = aT . 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 [18]. In this plot we have used two different values of the parameter a, where the value a = 0 corresponds to standard cold inflation.

IV. CONCLUSIONS
In the present work we have studied warm inflation with a quartic inflaton potential V (φ) = (1/4)λφ 4 and an inflaton decay rate proportional to the temperature, namely Γ = aT . Warm inflation consists an alternative to the standard cold inflation, during which radiation is neglected and which requires two steps, a slow-roll phase followed by a reheating phase, about which very little is known. On the contrary, in warm inflation, which has the attractive feature that avoids reheating, radiation is also taken into account and it is coupled to the inflaton leading to testable predictions different than the predictions of standard inflation even in the weak dissipative regime. The model we have considered is characterized by two parameters, namely the dimensionless couplings a and λ. We have used the latest Planck data to constrain the parameters of the model, and the results we have obtained are shown in the figures 1 and 2 for the case of weak and strong dissipative regime respectively. In the weak regime first, where Γ 3H, the background equations look the same as in standard inflation, however the tensor-to-scalar ratio is suppressed by the factor T /H, which must always be larger than one in warm inflation. The power spectrum constraint determines λ in terms of a, and then the tensor-to-scalar ratio as a function of the scalar index n s changes according to a as follows: As a increases the theoretical curve is shifted downwards, and on the other hand as a decreases the theoretical curve is shifted upwards. We have obtained both an upper and a lower bound on a, since when a becomes too low the theoretical curve lies outside the contours allowed by data, and when a becomes too large the condition for the weak dissipative regime is not satisfied. In figure 1 we show the contours allowed by the data together with four theoretical curves, namely one for the standard inflation and for three different values of the coupling a in warm inflation, the minimum value, the maximum value and one intermediate value.
In the strong dissipative regime, where Γ 3H, the power spectrum does not depend on a and so the constraint determines the inflaton self-interaction coupling λ. In the figure 2, the r − n s plot is shown and there is only a lower bound for a, obtained by the requirement of staying in the strong dissipative regime. In this regime the tensor-to-scalar ratio is suppressed by the factor T /H as in the weak regime, but also by the factor R 5/2 . That is why in the strong regime the model always predicts a very low r. By the other hand, we observe that the constraints found on the coupling λ, in both dissipative regimes, are in agreement with the value obtained in standard cold inflation using the COBE normalization. In this way, we conclude that warm inflation can rescue the quartic potential that in standard inflation is ruled out by the data.