Dynamics of warm power-law plateau inflation with a generalized inflaton decay rate: predictions and constraints after Planck 2015

In the present work, we study the consequences of considering a new family of single-field inflation models, called power-law plateau inflation, in the warm inflation framework. We consider the inflationary expansion is driven by a standard scalar field with a decay ratio Γ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Gamma $$\end{document} having a generic power-law dependence with the scalar field ϕ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\phi $$\end{document} and the temperature of the thermal bath T given by Γ(ϕ,T)=CϕTaϕa-1\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Gamma (\phi ,T)=C_{\phi }\frac{T^a}{\phi ^{a-1}}$$\end{document}. Assuming that our model evolves according to the strong dissipative regime, we study the background and perturbative dynamics, obtaining the most relevant inflationary observable as the scalar power spectrum, the scalar spectral index and its running and the tensor-to-scalar ratio. The free parameters characterizing our model are constrained by considering the essential condition for warm inflation, the conditions for the model evolves according to the strong dissipative regime and the 2015 Planck results through the ns\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$n_\mathrm{s}$$\end{document}–r plane. For completeness, we study the predictions in the ns\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$n_\mathrm{s}$$\end{document}–dns/dlnk\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathrm{d}n_\mathrm{s}/\mathrm{d}\ln k$$\end{document} plane. The model is consistent with a strong dissipative dynamics and predicts values for the tensor-to-scalar ratio and for the running of the scalar spectral index consistent with current bounds imposed by Planck and we conclude that the model is viable.


Introduction
The inflationary universe has become in the most acceptable framework in describing 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][7], inflation also genera e-mail: abduljawad@ciitlahore.edu.pk; jawadab181@yahoo.com b e-mail: nelson.videla@ing.uchile.cl c e-mail: fazi.gull@yahoo.com ates a mechanism to explain the large-scale structure (LSS) of the universe [8][9][10][11][12][13] and the origin of the anisotropies observed in the cosmic microwave background (CMB) radiation [14][15][16][17][18][19][20][21], since primordial density perturbations may be sourced from quantum fluctuations of the inflaton scalar field during the inflationary expansion. The standard cold inflation scenario is divided into two regimes: the slowroll and reheating phases. In the slow-roll period the universe undergoes an accelerated expansion and all interactions of the inflaton scalar field with other field degrees of freedom are typically neglected. Subsequently, a reheating period [22][23][24][25] is invoked to end the brief acceleration. After reheating, the universe is filled with relativistic particles and thus the universe enters in the radiation big-bang epoch.
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 [21]. 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 [19], in which r < 0.12 (95% CL). In particular, some representative models, as chaotic inflation, which predict a large value of the tensor-to-scalar ratio r , are ruled out by the data. As reported in Ref. [26], the Planck data tends to support plateau-like inflaton scalar potentials, which are asymptotically constant. Starobinsky R 2 [1] and Higgs inflation [27] are the most representative models with this class of potentials and, more recently, the α-attractors models [28,29]. For this last class of models, the approach to the inflationary plateau is exponential, being indistinguishable from Starobinsky and Higgs inflation models. In this direc-tion, starting from global supersymmetry and considering a superpotential, there was proposed in [30,31] a new class of models called shaft inflation. As opposed to Starobinsky and Higgs inflation, the approach to the plateau is of powerlaw type. In a subsequent work [32], another new family of inflationary models is studied, proposing a phenomenological potential where n and q are real parameters, V 0 is a constant density scale and α ≡ M M p , with M being a mass scale and M p = 2.43 × 10 18 GeV denotes the reduced Planck mass. In addition, it is required that φ ≥ M or equivalently φ ≥ α M p ; otherwise this model is indistinguishable from monomial inflation where V ∝ φ nq . It was demonstrated in the r -n s plane that the predictions of power-law plateau inflation are distinct and testable compared to several inflation models. In particular, the case n = 2 and q = 1 corresponds to the best model choice in Ref. [32]. Despite this attractiveness, in order to be a realistic model, it needs to be embedded in a convenient theoretical framework.
On the other hand, some classes of inflaton models excluded by current data in the standard cold inflation scenario which can be saved in the warm inflation scenario, which is an alternative mechanism for having successful inflation. The warm inflation scenario, as opposed to standard cold inflation, has the essential feature that a reheating phase is avoided at the end of the accelerated expansion due to the decay of the inflaton into radiation and particles during the slow-roll phase [33][34][35]. During warm inflation, the temperature of the universe does not drop dramatically and the universe can smoothly enter into the decelerated, radiation-dominated period, which is essential for a successful big-bang nucleosynthesis. In the warm inflation scenario, dissipative effects are important during the accelerated expansion, so that radiation production occurs concurrently with the accelerated expansion. For a representative list of recent references see Refs. [36][37][38][39][40][41][42][43][44][45]. The dissipative effect arises from a friction term or dissipative coefficient , which describes the processes of the scalar field dissipating into a thermal bath via its interaction with other field degrees of freedom. The effectiveness of warm inflation may be parametrized by the ratio R ≡ /3H . The weak dissipative regime for warm inflation is for R 1 while for R 1, it is the strong dissipative regime for warm inflation. It is important to emphasize that the dissipative coefficient may be computed from first principles in quantum field theory considering that encodes the microscopic physics resulting from the interactions between the inflaton and the other fields that can be present. For instance, by considering different decay mechanisms, it is possible to obtain several expressions for the dissipative coefficient . In particular, in Refs. [46][47][48], a supersymmetric model containing three superfields , X , and Y has been studied, with a superpotential W = f ( ) + g 2 X 2 + h 2 XY 2 , where the scalar components of the superfields are φ = √ 2 , χ , and y respectively. The inflaton scalar potential is given by 2 |, which spontaneously breaks supersymmetry (SUSY). By coupling the inflaton to the bosonic and fermionic X fields and their subsequent decay into Y scalars and fermions, which form the thermal bath, and for the case of low-temperature regime, when the mass of the catalyst field m χ is larger than the temperature T , the resulting dissipation coefficient can be well described by the expression where C φ is a dimensionless parameter related to the dissipative microscopic dynamics. For this particular case 1 and N X,Y denote the multiplicity of chiral superfields. In this direction, SUSY ensures that quantum and thermal corrections to the effective potential are under control [49,50]. As another example, in Ref. [51], it was demonstrated that a form = C φ φ 2 /T , in principle, with problems of large thermal corrections for the models studied in Refs. [52,53], may produce a consistent warm inflation scenario for a specific model. On the other hand, it was shown for the first time in Ref. [39] that warm inflation can be realized by directly coupling the inflaton to a few light fields, instead of considering indirect couplings to light fields through heavy mediator fields, as in Refs. [46][47][48]. Then the expression obtained for turns out to be = C φ T . Following Refs. [47,54,55], a general parametrization of the dissipative coefficient (T, φ) can be written as This expression for includes the several cases mentioned above. Specifically, for the value a = 3, this case corresponds to a low temperature regime, when the mass of the catalyst field m χ is larger than the temperature T [46][47][48].
Additionally, thermal fluctuations during the inflationary scenario may play a fundamental role in producing the primordial fluctuations [56][57][58]. In 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 [33][34][35][56][57][58][59][60][61][62]. When the universe heats up and becomes radiation dominated, inflation ends and the universe smoothly enters in the radiation bigbang phase [33,34]. For a comprehensive review of warm inflation; see Refs. [49,50]. In this direction, there are many phenomenological models of warm inflation, but more interesting are the first principles models of warm inflation in which the dissipative coefficient and effective potential are computed from quantum field theory. For instance, in Refs. [49,50], the following inflaton superpotential W ( ) was considered: which reduces to chaotic inflation models for p > 1. In particular, in Ref. [48] the authors studied the quartic potential V (φ) = λφ 4 , which corresponds to a superpotential W (φ) = λ 3 /3 together with the dissipative coefficient corresponding to a = 3, i.e. ∝ T 3 /φ 2 , was confronted with current data available at that time. On the other hand, when p = 0 and λ < 0, we have the standard supersymmetric hybrid inflation. Given the attractiveness of the power-law plateau inflation models as a new class of candidates in describing inflation, the main goal of this work is study the consequences of considering this new family of single-field inflation models in the warm inflation scenario in order to avoid the reheating phase. We would like to emphasize that our analysis is phenomenological in the sense that, in order to describe the dissipative effects during the inflationary expansion, we consider the generalized parametrization for the inflaton decay rate given by Eq. (2) and without considering a first principles construction for our model. However, in Ref. [32], the authors presented a toy-model in supergravity (SUGRA) which can produce the scalar potential of Plateau inflation (1) for n = 2 and q = 1. Specifically, they considered only global supersymmetry (SUSY) and sub-Planckian fields with the superpotential where S, 1 , 2 are chiral superfields and m is a large, but sub-Planckian scale. An interesting approach could be specify a decaying mechanism for the scalar components of the superfields 1 and 2 in light fields which form the thermal bath and compute the corresponding dissipative coefficient starting form first principles. However, this further considerations go beyond the scope of our present work, but these may be regarded as basis of a future work. On the other hand, we will restrict ourselves only to the study of the strong dis-sipative regime, R 1. For this dissipative regime, in the slow-roll approximation, we study the background as well as the perturbative dynamics. The free parameters characterizing our model are constrained by considering the essential condition for warm inflation, T > H , the condition for the model evolves according to strong dissipative regime and the 2015 Planck results through the n s -r plane. For completeness, we study the predictions of our model regarding the running of the scalar spectral index through the n s -dn s /d ln k plane.
This paper is organized as follows: In the next section, we present the basic setup of warm inflation. In Sect. 3, we study the background and perturbative dynamics when our model evolves according to strong regime. Specifically, we find explicit expressions for the most relevant inflationary observable as the scalar power spectrum, scalar spectral index, the running of the scalar spectral index and the tensor-to-scalar ratio. In order to establish a direct comparison between the power-law plateau inflation in the cold and warm scenarios, we will restrict ourselves to the case n = 2 and q = 1. For this particular case, we obtain the predictions in the n s -r and n s -dn s /d ln k plane.
Finally, Sect. 4 summarizes our finding and exhibits our conclusions. We have chosen units such that c =h = 1.

Basics of warm inflation scenario
In this section, we introduce the basic setup of warm inflation

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 equation reads is the reduced Planck mass. The dynamics of ρ φ and ρ γ is described by the equations [33][34][35] anḋ 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. As we have mentioned in the introduction, the parametrization given by Eq.
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. [33][34][35][56][57][58][59][60][61]. 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 [60]. Combining Eqs. (8) and (9) 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 con-ditions must be satisfied for a prolonged period of inflation to take place. For warm inflation, the slow-roll parameters are [46,60] The slow-roll conditions for warm inflation can be expressed as [46,60,61] When one these conditions is not longer satisfied, either the motion of the inflaton is no longer overdamped and slowroll 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.
The number of e-folds in the slow-roll approximation, using (5) and (8), yields where φ * and φ end are the values of the scalar field when the cosmological scales crosses the Hubble-radius and at the end of inflation, respectively. As it can be seen, the number of e-folds is increased due to an extra term of (1 + R). This implies a more amount of inflation, between these two values of the field, compared to cold inflation.

Cosmological 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 [49,50,57,60,61], 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 to leading order in the slow-roll approximation is given by [60,61] We also introduce the running of the scalar spectral index, which represents the scale dependence of the spectral index, by n run = dn s d ln k . In particular, for the strong dissipative regime, this expressions is given by [60,61] dn s d ln k where ζ 2 and γ are second-order slow-roll parameters defined by and respectively.
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 [62,63]. However, the tensor-to-scalar ratio r is modified with respect to standard cold inflation, yielding [49,50] 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 with a generic power-law dependence with the scalar field φ and the temperature of the thermal bath T influences the inflationary dynamics for the power-law plateau potential. We will restrict ourselves to the strong dissipation regime.

Background evolution
Assuming that the inflationary dynamics takes place in the strong dissipative regime, i.e., R 1 (or 3H ), by using Eqs. (2) and (11), the temperature of the thermal bath as a function of the inflaton field is found to be Replacing the last expression into Eq. (2), both the inflaton decay rate and the ratio R = /3H expressed in terms of the inflaton field becomes and respectively. In this way, by combining Eqs. (1), (8), and (22), the inflaton field as a function of cosmic time may be obtained from the following expression: where g 0 ≡ 24−8nq−8a+anq 2n(4+a) and 2 F 1 denotes the hypergeometric function [64].
For this model the set of slow-roll parameters become For the strong dissipative regime, the slow-roll conditions (13) become As we mentioned in the previous section, inflation ends when one of these parameters become the order of R.
On the other hand, 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) and (23) into (14), yielding where .
Since Eq. (30) has a complicated dependence in the inflaton field, it is not possible to express φ * as a function of N analytically. Instead, for numerical purposes, from Eq. (1) we may express the inflaton field as a function of the amplitude of the potential and evaluate this expression at the Hubbleradius crossing, obtaining The last expression will be useful to evaluate the several inflationary observables and put the observational bounds on our model.

Cosmological perturbations
Now, we shall study the cosmological perturbations for our model in the strong dissipative regime R = /3H > 1. For this regime, the amplitude of the scalar power spectrum (32) becomes then, by replacing Eqs. (1), (21), and (23), the power spectrum as a function of the inflaton field is found to be By considering the strong dissipative regime, the expression for the scalar spectral index becomes in this way, by replacing Eqs. (23) and (25)- (27), the inflaton field dependence of the scalar spectral index is given by Regarding the running of the scalar spectral index, the second-order slow-roll parameters ζ 2 and γ for this model become and γ = a × α n φ M p n + 4 28 + n(4 − 3q) + n 2 (4 + 3q) respectively. Then, by replacing last expressions together with Eqs. (25)- (27) into (17), the running of the spectral index is completely determined (not shown).
Regarding the tensor perturbations, the tensor-to-scalar ratio for the strong regime becomes The inflaton field dependence of the tensor-to-scalar ratio is determined by replacing Eqs. (1), (11), and (25) into (38), yielding In order to find observational constraints on our model, we will study the particular case n = 2 and q = 1, which corresponds to the best choice of model in the power-law plateau inflation in Ref. [32]. In addition, we consider α and V 0 , from the potential (1), and C φ , from the generalized inflaton decay ratio (2), as free parameters characterizing the model.

Special case n = 2 and q = 1
To compare the predictions of power-law plateau inflation in the cold and warm scenarios, we will restrict ourselves to the case n = 2 and q = 1, corresponding to the best choice of model in Ref. [32]. Moreover, for the generic parametrization of the inflaton decay rate (2), we consider the cases a = 3, 1, 0, and −1, which correspond to several dissipative ratios studied in the literature, but we consider α, V 0 and C φ to be free parameters. To put observational constraints on the parameters characterizing our model, we consider the essential condition for warm inflation, T > H , the condition for which the model evolves according to the strong regime, R 1, and finally the two-dimensional marginalized joint confidence contours for n s and r , at the 68 and 95% CL and the amplitude of the scalar power spectrum by Planck 2015 data [21]. In addition, we will try to ascertain whether the predictions for the running of the scalar spectral index dn s /d ln k are consistent with the current bounds imposed by Planck.

a = 3
In first place, for the special case a = 3, i.e., for ∝ T 3 /φ 2 , the scalar spectral index (35) becomes From last expression we see that n s is always greater than 1. Based on current observational data, for CDM +r + dn s /d ln k, the spectral index is measured to be n s = 0.9667± (68% CL, Planck TT + LowP). Hence, the inflaton decay ratio corresponding to a = 3 is not suitable to describe a strong dissipative warm inflation scenario consistent with current observations. It is interesting that, for other inflaton potentials, the inflaton decay ratio a = 3 describes a consistent warm inflationary dynamics (see Refs. [43,47,[65][66][67][68][69]). Figure 1 shows the ratio /3H and the tensor-to-scalar ratio r as functions of the scalar spectral index n s for the case a = 1, i.e., (φ, T ) = C φ T . To obtain the values to perform the plots, we have used three different values for C φ parameter, fixed the values C γ = 70 and M p = 1. For each value of C φ , we solve numerically the Eqs. (32) and (35) (after evaluating both equations at φ * given by Eq. (31), which is a function of V * ) for α and V 0 , considering the observational values P R 2 × 10 −9 , n s 0.9667 [21] and fixing V * = 10 −12 . In this way, for C φ = 4.39 × 10 −2 , we obtain the values α = 0.25 and V 0 = 1.03 × 10 −12 , whereas for C φ = 9.45 × 10 −2 , the solution is given by α = 0.45 and V 0 = 1.09 × 10 −12 .

a = 3
Finally, for C φ = 1.53 × 10 −1 , we found that α = 0.65 and V 0 = 1.19 × 10 −12 . In this way, the R(n s ) and r (n s ) curves of Fig. 1 may be generated by plotting Eqs. (23), (35), and (39) parametrically (after being evaluated at φ * given by Eq. (31)) with respect to V * . From the left panel, we observe that, for C φ ≥ 4.39 × 10 −2 , the model evolves according to the strong regime, R > 1. On the other hand, we find numerically that, for C φ ≥ 4.39 × 10 −2 , the ratio T H becomes T H 79 when n s 0.9667 (plot not shown). Hence, for C φ ≥ 4.39 × 10 −2 the essential condition for warm inflation, T H > 1, is always satisfied. Then the condition for which the model evolves in agreement with the strong regime gives us a lower limit on C φ . However, the essential condition for warm inflation does not impose any constraint on C φ . On the other hand, right panel of Fig. 1 shows the trajectories in the n s -r plane along with the two-dimensional marginalized constraints at 68 % and 95% C.L. on the parameters r and n s , by Planck 2015 data [21]. Here, we observe that, for C φ ≥ 4.39 × 10 −2 , the tensor-to-scalar ratio predicted by this model is always consistent with the observational bound found by Planck, given by r < 0.168 (95% CL, Planck TT + LowP). In order to determine the prediction of this model regarding the running of the spectral index, Fig. 2 shows the trajectories in the n s -dn s /d ln k plane. Again we note that for C φ 4.39×10 −2 the running of the spectral index predicted by the model is consistent with the bound found by Planck, given by dn s /d ln k = −0.0126 +0.0098 −0.0087 (68% CL, Planck TT + LowP). After the previous analysis, we only were able to find a lower limit on C φ as well as for α and V 0 , given by C φ = 4.39 × 10 −2 , α = 0.25 and V 0 = 1.03 × 10 −12 .
From the left panel of Fig. 3, the condition for the model evolves according to strong regime is satisfied for C φ ≥ 7.49 × 10 −7 , which gives us a lower limit for C φ . Additionally, for C φ ≥ 7.49 × 10 −7 the condition for warm inflation, T H > 1, is always satisfied. In particular, for C φ = 7.49 × 10 −7 , the ratio T H becomes T H 55 when n s 0.9667. Then, just like the case a = 1, for a = 0 the essential condition for warm inflation does not impose any constraint on C φ . Moreover, from the right panel, for C φ ≥ 7.49 × 10 −7 , the tensor-to-scalar ratio becomes r ∼ 10 −5 , but this value is still supported by the last data of Planck. For completeness, the running of the scalar spectral index becomes dn s /d ln k −0.001 at n s 0.9667 (plot not shown). Then, for the case a = 1, the previous analysis gives us only a lower limit for C φ as well as for α and V 0 , given by C φ = 7.49×10 −7 , α = 0.4, and V 0 = 1.02×10 −12 , respectively. Despite this result, it is interesting to mention that, for this power-law plateau potential, the inflaton decay ratio ∝ φ describes a strong dissipative warm inflation scenario compatible with current observations. In previous works [65][66][67][68][69] it was found that the decay rate ∝ φ is not able to describe a consistent strong dissipate dynamics, since the predicted scalar spectral index is always greater than unity.
Finally, for C φ = 2.53 × 10 −12 , we found that α = 2.4 and V 0 = 1.13 × 10 −12 . Figure 4 shows the plots of the ratios R = /3H and T /H as functions of the scalar spectral index n s . From the left panel, we see that for C φ ≥ 2.01×10 −12 the model takes place in the strong dissipative regime of warm inflation. Moreover, from the right panel, we observe that the essential condition for warm inflation T H > 1 is always guaranteed. In particular, for C φ = 2.01 × 10 −12 , this ratio takes the value T H 60 when n s 0.9667. Again, the T H plot does not impose any constraint on C φ . Regarding the predictions of this case in the r -n s plane, for C φ = 2.01 × 10 −12 , the tensor-to-scalar ratio becomes r ∼ 10 −5 , but this value is still supported by the last data of Planck by considering the two-dimensional marginalized joint confidence contours for (n s , r ), at the 68 and 95% C.L. (plot not shown). Finally, the predictions for the running of the spectral index are similar to previous ones, yielding dn s /d ln k −0.001 at n s 0.9667 for all the values considered for C φ (plot not shown). The result of this analysis yields only a lower limit for C φ as well as for α and V 0 , given by C φ = 2.01 × 10 −12 , α = 1.5, and V 0 = 1.07 × 10 −12 , respectively. Just like the case a = 0, the case a = −1 has an interesting feature, because yields a strong dissipative dynamics compatible with observations, since that in previous works [52,53,[65][66][67][68][69], the inflaton decay rate ∝ φ 2 T is not able to describe a consistent strong dissipative dynamics.

Discussion
From the analysis carried out in Ref. [32], and provided that the model be distinguishable from monomial inflation, i.e., φ > αM p , the authors found that the best model choice in the power-law plateau inflation family correspond to the values n = 2 and q = 1. In a first approach and ensuring a sub-Planckian excursion for the inflaton through the potential, the maximum value allowed for α was found to be α = 0.04, and the values for the scalar spectral index and the tensor-to-scalar ratio at the Hubble-radius crossing correspond to n s = 0.97 and r = 0.000157. For our warm power-plateau model, in order to produce a strong dissipative dynamics, for each value of a, all the values obtained for α are greater than 0.04, implying a trans-Planckian excursion of the inflaton field, but ensuring that φ > αM p . On the other hand, the predictions for the scalar spectral index are very similar for the cold and warm power-law plateau inflation models. Regarding the tensor-to-scalar ratio, in the warm inflation scenario, this quantity is suppressed by a factor (T /H )R 5/2 > 1 compared with standard cold inflation. In particular, for a = 1, the tensor-to-scalar ratio is almost the same order compared to cold power-law plateau inflation. However, for a = 0 and a = −1, the tensor-to-scalar ratio becomes smaller than for cold power-law plateau inflation.
In a second approach adopted in Ref. [32], the authors considered a trans-Planckian excursion of the inflaton field, obtaining values for α going from α = 1 up to α = 5, and the tensor-to-scalar ratio and the running of the scalar spectral index taking values from r = 0.004106 up to r = 0.024412, and from dn s d ln k = −0.00057 up to dn s d ln k = −0.00051, respectively. This implies that, for any value of a, the tensorto-scalar ratio for our warm power-law plateau inflation is always lower than the predicted by the cold scenario. On the other hand, the values predicted for the scalar spectral index in the cold and warm scenarios are very similar. In addition, it is interesting to mention that the running of the scalar spectral index dn s d ln k predicted by our warm power-law plateau model is almost two orders of magnitude greater than the predicted by the cold scenario.
After the analysis performed previously for each value of a, we only found a lower limit on C φ as well as for α and V 0 , which means that we have a larger range of parameter values to enter in accordance with the Planck results and consistent with a strong dissipative dynamics. This degeneracy could be broken combining these results with the constraints on the inflationary observables related with non-Gaussianities, particularly the f N L parameter, since in warm inflation scenario these have different features when comparing with cold inflation [36]. Despite this issue, the predictions of warm powerlaw plateau inflation are comparable to those of power-law plateau cold inflation, however, the difference between both scenarios is a way to address the problem of reheating in cold power-law plateau inflation is provided by the warm inflation scenario.

Conclusions
In the present work, we have studied the consequences of considering a new family of single-field inflation models, called power-law plateau inflation, in the warm inflation scenario. As far we know, this is the first work in studying the dynamics of warm inflation by using the power-law plateau potential. In order to describe the dissipative effects during the inflationary expansion, we considered a generalized expression for the inflaton decay ratio given by (φ, T ) = C φ T a φ a−1 , where a = 3, 1, 0, −1, denotes several inflaton decay ratios studied in the literature. We restricted ourselves only to the study of the strong dissipative regime, R 1. For this dissipative regime, under the slow-roll approximation, we have studied the background as well as the perturbative dynamics. In particular, we have found the expressions for the scalar power spectrum, scalar spectral index and its running as well as the tensor-to-scalar ratio. Contrary to the standard cold inflation, in the warm inflation scenario it is not sufficient to consider only the constraints on the r -n s plane, but we also have to consider the essential condition for warm inflation T > H and the conditions for the model evolves under strong dissipative regime R 1. For completeness, we study the predictions of our model regarding the running of the scalar spectral index, through the n s -dn s /d ln k plane.
To compare the predictions of power-law plateau inflation in the cold and warm scenarios, we restricted ourselves to the case n = 2 and q = 1, corresponding to the best choice of model in Ref. [32]. For this particular case, the inflaton decay a = 3, i.e. ∝ T 3 φ 2 , fails in describe a strong dissipative dynamics consistent with current data, since the predicted value for the scalar spectral index is always greater than unity. We recall that, for the more representative potentials studied in the literature, the inflaton decay rate a = 3 describes a warm inflationary dynamics consistent with current data. Regarding the predictions in the n s − r and n s -dn s /d ln k planes, for a = 1, the tensor-to-scalar ratio and the running of the spectral index becomes r 10 −4 and dn s d ln k 0.002, respectively, whereas for both the cases a = 0 and a = −1, these inflationary observables become r 10 −5 and dn s d ln k 0.001, being consistent with current bounds imposed by Planck for CDM +r + dn s /d ln k. It is interesting to mention that, for other kinds of potentials already studied in the warm inflaton scenarios, the decay ratios a = 0 and a = −1 predicted a scalar spectral index always greater than unity. On the other hand, for any value of m, the condition for the model evolves according to the strong dissipative regime sets the lower limit for the disipative parameter C φ as well for α and V 0 . However, the essential condition for warm inflation to occur, T > H neither the Planck data, by considering the two-dimensional marginalized constraints at 68% and 95% C.L. on the parameters r and n s , do not impose any constraints on the model for this dissipative regime, obtaining a lower limit on C φ as well as for α and V 0 . However, if we consider the observational constraints on the inflationary observables related with non-Gaussianities, particularly the f N L parameter, this degenerancy in the parameters could be broken.
Comparing our warm power-law plateau inflation model with the standard one, we found that the strong dissipative warm inflation dynamics is only consistent with a trans-Planckian incursion of the inflaton potential, according with second approach adopted in [32], ensuring that this powerlaw plateau potential be distinguishable from monomial inflation. For this trans-Planckian evolution of the inflaton, and for any value of a, the tensor-to-scalar ratio for our warm power-law plateau inflation is always lower than predicted by the cold scenario. On the other hand, the values predicted for the scalar spectral index in the cold and warm scenarios become similar, however, the running of the scalar spectral index dn s d ln k is almost two orders of magnitude greater than predicted by the cold scenario. We have shown that warm power-law plateau inflaton, with decay ratios parametrized by a = 1, 0, and −1, is consistent with a strong dissipative dynamics and predicts values for the scalar spectral index, the running of the scalar spectral index, and the tensor-to-scalar ratio consistent with current bounds imposed by Planck, for CDM +r + dn s /d ln k.