1. Introduction
The final release of the
Planck Cosmic Microwave Background (CMB) data [
1] has severely constrained many primordial observables, such as the scalar spectral index,
, the tensor-to-scalar ratio,
r, the primordial non-Gaussianities and the primordial isocurvature spectrum, leaving us with the simplest vanilla inflationary models as the most preferred ones. However, the presence of the running,
, and the running of running,
, of the scalar spectral index that were inferred by the
Planck data may hint otherwise. For the ΛCDM model, the
Planck 2018 TT(TT,TE,EE)+lowE+lensing data constrains the scalar spectral index, its running and the running of its running as [
1]
all at
CL. According to
Planck, such values for the running and running of running of the scalar spectral index yield a better fit to the low-
ℓ deficit of the TT spectrum. It is interesting to note that in these constraints there is a slight preference for a positive and rather large
, which is also larger than the running of the spectral index (
). This is unexpected, when compared, for instance, to the values derived from standard inflationary models [
2], where a hierarchy such as
is expected. It has been shown in Refs. [
3,
4] that if we consider the empirical relation
, where
N is the number of
foldings, then in these vanilla models one expects
The
Planck-preferred vanilla inflationary models in general also produce negative
and
. Therefore, such inflation models have been considered to be in tension with the current observations by
Planck [
4,
5]. Going beyond these vanilla models does not seem to improve the scenario either. It was shown in Ref. [
5] that single field inflation models and non-interacting two-field models are incapable of producing a
larger than
. Slow-roll violating models and those with non-trivial evolution of sound speed could be able to provide an exception to this trend, but with a considerable amount of fine-tuning [
5].
The
Planck-preferred vanilla inflaton potential models in general consider the traditional cold inflation (CI) scenario, where the dynamics of the inflaton field
is assumed to be independent of its couplings with other fields during inflation. However, such interaction terms play a major role at the end of inflation, when the inflaton field needs to release its energy density in the form of radiation through decay processes to reheat the Universe and, hence, leading to a radiation dominated Universe as required by the big bang cosmology. On the other hand, in the warm inflation (WI) scenario [
6], the couplings between the inflaton and other fields are strong enough such that their effects on the inflaton dynamics cannot be ignored. In WI (for reviews on WI, see, e.g., Refs. [
7,
8]), the inflaton field is able to keep dissipating its energy such that a non-negligible radiation bath throughout inflation can be produced, while preserving the flatness of the inflaton potential required for slow-rolling of the inflaton field. The presence of such a quasi-equilibrium thermal radiation bath during WI helps transiting smoothly from the inflationary accelerated phase to the radiation dominated phase after inflation ends, avoiding a phase of reheating in between. In WI, the inflaton dynamics are modified with respect to the one in CI due to the presence of an extra friction term,
, that accounts for the energy transfer between the inflaton field and the radiation bath present during inflation. Due to such modified dynamics, the inflationary observables, such as
r,
and the non-Gaussianity parameter
, are also modified with respect to those obtained in CI. These changes have certain advantages. For example, some of the inflaton potentials excluded by data in the context of CI can be made in tune with the observations in the WI context. One such example is the quartic chaotic potential
[
9]. Moreover, it has recently been shown that, while CI fails to comply with the recently proposed Swampland Conjectures in String Theory [
10,
11,
12], WI can easily accommodate those criteria [
13,
14,
15,
16,
17,
18]. Hence, WI provides a way to construct inflationary models that can be consistent as effective models that could descend from an ultraviolet complete quantum gravity, despite the swampland conjectures barring them from constructing de Sitter vacua in String Landscapes.
The running and the running of the running of the scalar spectral index in the context of WI were first studied in Ref. [
19], where several inflationary potentials were analyzed with two different forms of dissipative terms considered in WI,
and
(
T being the temperature of the thermal bath). The WI models studied in Ref. [
19] were, however, treated in the weak dissipative regime (
, where
Q is the ratio of the thermal friction term
and the Hubble expansion friction term
present in the inflaton equation of motion in WI). This was because the models studied in that reference could only lead to consistent observables, e.g., values for
r and
, in that specific dissipation regime of WI. In particular, it was shown in Ref. [
19] that in all the models studied there, in the weak dissipative regime, there was still a large hierarchy between the values of
and
.
The aim of this paper is to study the running and the running of the running of the scalar spectral index in the context of WI, where WI is realized in a strong dissipative regime
. This is motivated by the recent results in WI in the context of the swampland conjectures [
13,
14,
15,
16] and also on the solution of the so-called
-problem [
20], which exactly favors WI being realized in the strong dissipative regime. For this purpose, we will focus mostly on a certain WI model, dubbed the Minimal Warm Inflation (MWI) [
21], where WI has been shown to be possible to be realized in the strong dissipative regime and in which the dissipative term is proportional to the cubic power of the temperature of the thermal bath
. However, our results will be kept as general as possible, so that they can be extended to other models. In the original paper of MWI [
21], a hybrid potential was used to produce a red-tilted scalar power spectrum. Later, MWI was studied with generalized exponential potentials in Ref. [
22], and it was shown that this model is not only in accordance with the current observations (yielding the appropriate values for
and
r), but also are in tune with the swampland conjectures. We will analyze both the original MWI model, with the hybrid potential and the one with the generalized exponential potential to derive
and
that are produced in these two scenarios. Note that although the running and the running of the running of the scalar spectral index in the context of WI were first studied in Ref. [
19], they were estimated numerically. Here, however, we aim to produce explicit analytical expressions for these quantities in WI. To the best of our knowledge, this is the first time that such analytical analysis and derivation of
and
in the context of WI are presented. Given the advent of new generations of cosmological observatories probing both the cosmic microwave background (CMB) measurements [
23,
24], the distribution of matter at low-redshift from optical, near-infrared and 21 cm intensity surveys [
25], it is expected that theories of cosmic inflation will be further constrained with the more precise cosmological data. Thus, it is important to have such analytical expressions to help in our ability in finding possible models of interest and also in model building in WI.
The remainder of the paper is organized as follows. In
Section 2, we briefly review WI and present the general expressions for
,
and
. In
Section 3, we analyze the case of MWI with a hybrid potential. In
Section 4, we turn our attention to the determination of
,
and
in the recently studied generalized exponential potentials for the inflation in WI. Then, in
Section 5, we discuss our findings and conclude.
2. Running and Running of the Running of Scalar Spectral Index in WI—General Expressions
The equation of motion of the inflaton field
in WI is
where
Q is the ratio of the two frictional terms, the friction due thermal bath and the friction due to Hubble expansion, present in the theory,
where the dissipation coefficient
is in general a function of the temperature of the thermal radiation bath generated and can also depend on the inflaton amplitude (for examples of different forms of dissipation coefficients derived in quantum field theory and used in WI, see, e.g., Ref. [
26]). In Equation (
3),
, where
V is the potential of the inflaton field. As the inflaton dissipates part of its energy density, it can sustain a radiation bath, with energy density
and whose evolution equation is given by
We define the slow-roll parameters
and
in the usual way,
where
GeV is the reduced Planck mass. In WI, another set of slow-roll parameters are also defined,
due to the fact that WI actually ends when
, or, similarly, when
, while in CI the usual condition for the end of the accelerated inflationary regime is simply
.
Besides
and
, it is also useful to define the higher order slow-roll coefficients
and
, as [
27]
and
which will appear in the equations for the running and for the running of the running to be derived later on.
The primordial scalar curvature power spectrum of a typical WI model can be written as [
19,
28]
where the subindex * means that the quantities are evaluated at the Hubble radius crossing of the pivot scale
. Here, the function
is given as
where
is the thermal distribution of the inflaton field due to the presence of the radiation bath and
accounts for the effect of the coupling of the inflaton and radiation fluctuations [
29,
30,
31].
The primordial scalar power spectrum
can be expanded in terms of the small scale-dependence as
where
is the scalar amplitude,
is the scalar tilt,
is the running, and
is the running of the running (or second running).
From Equation (
14), the scalar spectral index
can be calculated at the horizon crossing (
) as
Likewise, the expressions for
and
are determined, respectively, by
and
Using Hubble radius crossing
, we can write
, with
being the number of
e-folds and
, then
The definitions given by Equations (
15)–(
17) can now be applied directly to the scalar curvature in WI, Equation (
12). Once an appropriate functional form for the dissipation coefficient is given, explicit expressions for
,
and
can be derived. Here, we work with the well-motivated functional form for the dissipation coefficient
1
where
is a constant,
T is the temperature, and
is some appropriate mass scale. For specific examples of microscopic quantum field theory derivations of such dissipation coefficients in WI, see, e.g., Refs. [
7,
34,
35,
36,
37].
When taking derivatives of Equation (
12) such as to obtain
,
and
, we are naturally faced with accounting for derivatives on
Q and
with respect to
N. These expressions depend on the specific form of the dissipation coefficient used. For completeness, let us quote them here for the generic dissipation term Equation (
19). We use the slow-roll approximation for the background dynamical equations for the inflaton field and for the radiation energy density in WI, e.g.,
From this set of equations, one can easily deduce for instance that
Thus, from Equations (
20)–(
22), together with Equation (
19), we obtain that
and
where we have also introduced the quantity
, defined as
. From Equations (
24) and (
25), we find in particular that
where we have defined
as
The expressions for
and
are obtained by taking further derivatives of Equation (
26) using Equation (
18). Though this is a lengthy but straightforward exercise, the expressions are, however, too long to quote here in full. However, since we are mostly interested in models of WI in the strong dissipative regime, we can expand Equation (
26) for
and obtain, at least for
, a relatively shorter expression. Hence, expanding Equation (
26) for
, we obtain that
where we have also used Equations (
24) and (
25). Approximate expressions for
and
can be obtained analogously. For instance, from Equation (
28) we obtain for
that
while for
we find
In the above expressions,
and
. The above approximate expressions for
,
and
applies to any WI model with a generic dissipation coefficient of the form of Equation (
19) and once the primordial potential for the inflaton is specified.
Below, as an example, we will work with two explicit models in WI, which have been shown to be of interest recently. Then, specific expressions for , and will be derived for those models, along with their respective analysis.
3. Running and Running of the Running of the Scalar Spectral Index in MWI
The minimal WI (MWI) model proposed in Ref. [
21] (see also Ref. [
38]) successfully realizes WI in the strong dissipative regime. In such a model, the inflaton is being treated as an axionic field with coupling to the non-Abelian gauge fields. Due to such gauge couplings, sphaleron transitions between gauge vacua at high temperatures lead to a friction term of the form
where
is the sphaleron rate,
f is the axion decay rate,
T is the temperature,
, with
g being the Yang–Mills gauge coupling, and
is a dimensionless quantity depending on the dimension of the gauge group
, the representation of the fermions
and on the gauge coupling
g through
. For instance, for a quantum chromodynamics (QCD) type of axion and using typical values of parameters for QCD, we have that
.
This model is particularly attractive in the context of WI. Because of the axionic shift symmetry, the inflaton is protected from any perturbative backreaction and, hence, from acquiring a large thermal mass. Similar symmetry properties allowing the inflaton to be coupled directly to the radiation fields have also been studied in Refs. [
20,
37]. Models based on pseudo-Goldstone bosons for the inflaton [
39] are quite reminiscent of these ideas and have gained increased attention recently in the context of WI [
40].
The scalar power spectrum in the MWI model that was considered in Ref. [
21] was given by
where
. At the end of this section and in the next one, we will discuss a more accurate form for the scalar spectrum in this model for a dissipation coefficient of the form
as it has been considered here. However, for now, let us use the above expression Equation (
32) as considered by the authors in Ref. [
21].
Since we are also deriving expressions in the high dissipative regime,
, let us note that in this case, using
and
in Equation (
19), that
Hence, the scalar spectral index, in the strong dissipative regime can be calculated as
For a detailed derivation of the scalar spectral index see, e.g., Ref. [
22].
It is evident from Equation (
37) that to obtain the observed red tilt of the scalar spectral index
in MWI, one requires a potential which yields
. In Ref. [
21], for this objective a potential of a hybrid inflation model was considered. The effective potential of hybrid inflation contains the inflaton field
along with a waterfall field
, which can be written as
2 The
field has a squared mass term as
in this potential. When
, which is the region where inflation takes place,
has only one minimum at
, and the effective potential becomes
3 During this stage, the constant term
drives the expansion. Below the threshold
, the waterfall field
quickly rolls down to its minimum
and puts an end to inflation. The effective potential, given in Equation (
39), yields the slow-roll parameters as
4
which shows that
Therefore, such a potential is suitable to produce a red-tilted scalar spectral index in MWI. Hence, with the hybrid inflation potential given in Equation (
39), the scalar spectral index in Equation (
37) can be approximated as
Above all, as we can see from Equation (
40) that
is effectively constant in this model, we denote
, where
is a constant. Hence, to yield the observed value of the red-tilt,
, we require
We now estimate the running
and the running of the running
of the scalar spectral index in this model. The running of the scalar spectral index,
, can be calculated as
where derivative has been taken with respect to the number of
e-foldings
N, and in MWI we have [
22]
However, as the effective hybrid potential given in Equation (
39) yields
, we can approximate the above equation as
This yields the running of the scalar spectral index in MWI as
Now, the value of
, which yields the observed red-tilt of the scalar spectral index given in Equation (
43), produces a running of the scalar spectral index as
The running of the running of the scalar spectral index then turns out to be
Again, the value of
given in Equation (
43), which yields the observed red-tilt of the scalar spectral index, gives rise to a running of the running of the scalar spectral index as
Therefore, we note that the MWI model with hybrid-type potential generates negative
and
, and the running of the running is two orders of magnitude smaller than the running of the scalar spectral index. These results come from the simplified analysis made in particular with the approximated scalar power spectrum given by Equation (
32).
Let us now check whether an improved form for the scalar power spectrum might somehow improve on the estimates made above for
,
and
. Using the previously derived Equations (
28)–(
30), we obtain for instance that
and
where
is given by Equation (
39), and the function
is derived using the growth function
as defined in the next section and given by Equation (
57).
We solve the complete system of dynamical equations for the hybrid inflation model and compare the numerical results with the analytical ones coming from Equations (
51)–(
53). As in Ref. [
21], we consider the case of the model consisting of a pure
, with relativistic degrees of freedom
, consisting of two polarizations per eight gauge bosons plus one for the axion. The gauge coupling is assumed to be
and
for the coefficient in the dissipation coefficient in Equation (
31). As also shown in Ref. [
21], the strong dissipative regime (
in the model with the hybrid type of potential Equation (
38) is favored when the parameters
M,
g and axion term
f satisfies
. We also recall that when assuming a QCD axion, astrophysical constraints [
43] put a lower bound
GeV on the QCD axion decay constant. In our numerical experiments, we have checked that both the strong dissipative regime and acceptable values for
f satisfying the astrophysical bounds can be achieved by taking, for example,
GeV and choosing appropriate values for
g and for the other constant,
, appearing in the potential Equation (
38).
We have studied two representative cases, and in the
Table 1 we summarize the relevant parameters and quantities obtained from the numerical analysis. Note that for the parameters considered, we have the temperature at the end of the WI slightly smaller than the axion decay constant,
, which indicates an axion symmetry breaking still in the inflationary regime and, hence, prior to the end of WI. As an aside regarding the hybrid inflation model studied here, we note that the inflaton field excursions for examples shown in
Table 1 are all sub-Planckian,
. Hence, the hybrid WI model studied here satisfies the swampland distance conjecture, while it marginally satisfies the (refined) de Sitter conjecture [
10,
11] by having
, but
.
The number of
e-folds
at Hubble radius crossing is found from the relation obtained in WI (see, e.g., Ref. [
22] for details)
where
is the temperature at the end of WI,
is the Hubble parameter today and for which we assume the
Planck result,
s
[from the
Planck Collaboration [
44], TT,TE,EE-lowE+lensing+BAO 68% limits,
s
],
is the CMB temperature today,
, while the for the pivot scale
we take the
Planck value
, and we also use the convention
. For
, we assume the MWI value,
.
The numerical and analytical results obtained for
,
and
are given in
Table 2. We note from the results from
Table 2 that the analytical and numerical results for
and
agree quite well. There are discrepancies between the values for
, especially for the lower dissipative case studied,
, but that can be attributed to the lack of numerical precision in the determination of
in that case, recalling that to obtain
numerically, it requires three derivatives of the scalar power spectrum with respect to the number of e-folds. In any case, these results also corroborate the previous ones using the approximated form of the power spectrum, with
and
having the same order of magnitudes as obtained previously, Equations (
48) and (
50), respectively. We have not quoted the results for the tensor-to-scalar ratio
r, but as typical for WI in the strong dissipative regime,
r is too small to be accessible by any future observation. For the dissipation values quoted in
Table 1,
.
We now study another type of WI models that can allow for large dissipation and are still consistent with the Planck observations. The objective is to check whether with these models we can break the large hierarchy between and .
4. Running and Running of the Running of the Scalar Spectral Index in MWI with Exponential Potentials
In CI, steep exponential potentials of the form
lead to power-law type inflation [
45], where inflation does not exit gracefully in standard general relativity. However, it was shown in Ref. [
46] that WI with dissipative coefficient of the form
, with
, can gracefully exit inflation with such an exponential potential. The MWI model [
21] is, thus, such a model where the dissipative coefficient scales with the cubic power of the temperature of the thermal bath,
.
A MWI model with such an exponential potential was studied in Ref. [
47]
5. Though inflation ends gracefully in such a model, it was shown there that such a combination yields way too large red-tilt in the scalar spectral index to be compliant with the current observation. Therefore, such models are not viable models of inflation.
A generalized form of the exponential potential, given as
with
, was then studied in Ref. [
22], which shows that this model not only exits inflation gracefully, but also fully satisfies all the observational constraints (comply with the observed values of scalar spectral index
and the tensor-to-scalar ratio
r) and is in tune with the swampland conjectures (see, e.g., Table 1 of Ref. [
22]). This same model was also considered in the context of WI for quintessential inflation [
49], displaying interesting features as far as the late Universe physics is concerned. Therefore, here we will further analyze such a model to determine the running and running of the running of the scalar spectral index that are yielded by this model.
In the model that we are presently considering, the explicit form for the function
appearing in the scalar of curvature power spectrum in WI, Equation (
12), has been derived in Ref. [
22], and it is given as
6We will first derive the scalar spectral index
, its running
and running of its running
analytically in this model. Note that, although the values of
in this model have been reported previously in Table 1 of Ref. [
22], they have been obtained numerically. To the best of our knowledge, this is the first time anyone has determined
,
and
in any WI model analytically.
If we ignore the thermal distribution factor
of the inflaton field in Equation (
13), which can be justified when working in the strong dissipative regime
(see, e.g., Ref. [
28]), then the scalar power spectrum simplifies to the form
This is the form of the power spectrum that we will use to calculate the spectral index and its running and running of its running analytically. Such an approximation of the power spectrum will later be justified by obtaining
,
and
numerically and matching them with the analytically obtained results.
First, we calculate the scalar spectral index
. Using the scalar spectrum given in Equation (
58), we obtain
From the equations in Equation (
36), after some straightforward algebra, we obtain the scalar spectral index as
where
has been defined in Equation (
27).
We now analytically determine
. Using the form of
from Equation (
60), we obtain
Using now that
the expression for the running of the scalar index can be computed to be
To determine
, one requires
Using the form of
from Equation (
64),
can be written as
The equations for
,
and
can also be expressed explicitly in terms of the inflaton field
using the potential Equation (
56), but that only complicates the form of the expressions more, and we do not need to perform that explicitly here. In
Table 3, we will explicitly furnish the numerical values for the required slow-roll parameters.
We now estimate the values of
,
and
using the analytical forms of these quantities that we derived above. To obtain these values, we need to fix
n,
and
, the values of which have been taken from Table 1 of Ref. [
22] which produce
and
r that are in tune with the observations
7. In the table, we have also quoted the temperature at the end of WI,
and the value obtained for the axion decay constant. To obtain
f, we made the same considerations regarding the gauge sector coupled to the scalar field as in the previous section for the MWI model with the hybrid potential. We have considered QCD like values for the parameters. With relativistic degrees of freedom
, the gauge coupling is assumed to be
and
for the coefficient in the dissipation coefficient in Equation (
31). Note that in all the cases considered here, we have that
, pointing to cases where the axion is already in a symmetry broken phase at the end of WI.
In
Table 3, we give the values of the required slow-roll parameters required in the derivation of
,
and
. We also give the normalization
of the potential for each of the models and the scale of the inflaton potential at Hubble radius crossing,
. The analytically estimated values for
,
and
that are obtained for each of those models are then shown in
Table 4. In the same table, the numerically obtained values for
,
and
are given, where the full scalar power spectrum, given in Equation (
12), has been used to determine them.
From the results shown in
Table 4, we notice that the analytical results match well with the numerical results. Therefore, the approximation we made to the scalar power spectrum in Equation (
58) is justified. We also note that the MWI model with generalized exponential potentials also produces negative
and
. However, in comparison to the MWI model with hybrid potential, it produces one order higher
and two-order higher
. These results are promising in the sense that they indicate that the large hierarchy observed between
and
, which appears in the vanilla CI models, can be broken in the context of WI.