Warm Little Inflaton becomes Dark Energy

TL;DR

This paper proposes a model where the warm inflation inflaton field transitions into dark energy, explaining late-time cosmic acceleration through a thawing quintessence scenario with specific observational signatures.

Contribution

It introduces a unified framework where warm inflation naturally evolves into dark energy without requiring separate mechanisms.

Findings

The model predicts a modified tensor-to-scalar ratio consistency relation.

It shows a variation in dark energy equation of state at low redshifts.

Dark energy isocurvature perturbations are correlated with primordial fluctuations.

Abstract

We present a model where the inflaton field behaves like quintessence at late times, generating the present phase of accelerated expansion. This is achieved within the framework of warm inflation, in particular the Warm Little Inflaton scenario, where the underlying symmetries guarantee a successful inflationary period in a warm regime sustained by dissipative effects without significant backreaction on the scalar potential. This yields a smooth transition into a radiation-dominated epoch, at which point dissipative effects naturally shut down as the temperature drops below the mass of the fermions directly coupled to the inflaton. The post-inflationary dynamics is then analogous to a thawing quintessence scenario, with no kination phase at the end of inflation. Observational signatures of this scenario include the modified consistency relation between the tensor-to-scalar ratio and tensor spectral index typical of warm inflation models, the variation of the dark energy equation of state at low redshifts characteristic of thawing quintessence scenarios, and correlated dark energy isocurvature perturbations.

Figures

Equations

$${\varphi }_1=\frac{M}{\sqrt{2}}{e}^{i\varphi /M},{\varphi }_2=\frac{M}{\sqrt{2}}{e}^{-i\varphi /M}.$$

$$-{\mathcal{L}}_{\varphi \psi }$$

$${\rm Y}\approx \frac{2g^2{h}^2}{(2\pi {)}^3}T\left(\frac{1}{5}-\ln \left(\frac{m_{\psi }}{T}\right)\right),m_{\psi }/T\ll 1,$$

$${\rm Y}\approx \frac{2g^2{h}^2}{(2\pi {)}^3}T\sqrt{\frac{\pi }{2}}\sqrt{\frac{m_{\psi }}{T}}{e}^{-\frac{m_{\psi }}{T}},m_{\psi }/T\gg 1.$$

$$\ddot{\varphi }+(3H+{\rm Y})\dot{\varphi }+V^{\prime }(\varphi )=0,$$

$${\dot{\rho }}_R+4H{\rho }_R={\rm Y}{\dot{\varphi }}^2.$$

$$V(\phi)=\left\{\begin{array}[]{c}\lambda(\phi^{4}+m^{4})\\ \frac{\lambda m^{4}}{1+\left(\frac{\phi}{m}\right)^{\alpha}}\\ \end{array}\right.~{}\begin{array}[]{l}\phi\leq 0\\ \phi>0\end{array}\,~{},$$

$$\frac{{\rho }_R}{V(\varphi )}\simeq \frac{1}{2}\frac{{ϵ}_{\varphi }}{1+Q}\frac{Q}{1+Q},$$

$${\varphi }^{\prime \prime }+\left(3+\frac{H^{\prime }}{H}\right){\varphi }^{\prime }=0,$$

$$\varphi ={\varphi }_0^{\prime }\left(1-e^{-(N_e-N_0)}\right),$$

$$w_{\varphi }=\frac{K(\varphi )-V(\varphi )}{K(\varphi )+V(\varphi )}\approx -1+\frac{K(\varphi )}{V(\varphi )},$$

Full text

Since the discovery of cosmic acceleration by the observations of Type Ia supernovae (SNe) [1, 2], dark energy (DE) has become an essential part of the standard cosmological model. Whether DE is a complement to the well-known energy budget of the Universe described by General Relativity or generated by modifications of the latter remains to be seen.

Microscopic/first principles descriptions of either scenario have remained elusive, but these appear crucial if one is to have a more complete understanding of the problem of cosmic acceleration. We think, however, that this is not emphasized enough: in the particular case of quintessence, a scalar field description of dynamical dark energy, one does not usually motivate neither the origin nor the characteristics of such a homogeneous scalar field.

Interestingly, it is also hypothesized that a homogeneous scalar field, the inflaton $\phi$, generates a phase of accelerated expansion in the early Universe, thereby solving the homogeneity and horizon problems and providing the seeds for structure formation and Cosmic Microwave Background (CMB) anisotropies [3]. One could then ask: can the same field generate these two similar stages of accelerated expansion, happening at such different energy scales?

Inflaton-quintessence unification was originally suggested in [4, 5] and investigated in more general settings [6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 38, 39, 37]. It may have noteworthy cosmological consequences as the generation of a baryon asymmetry [40, 41] and the modification of the constants of nature [42].

These models rely on reheating mechanisms, as gravitational particle production [43], which occurs due to the changing background geometry (sudden transition between inflation and kination) or instant preheating [44], which occurs through a two-stage mechanism $\phi\rightarrow\chi\rightarrow\psi$, where $\chi$ is either a boson or fermion of mass $m_{\chi}=g\phi$ (with $10^{-7}<g<1$) and $\psi$ are light fermions. It was pointed out in [5, 44] that gravitational particle production is very inefficient. Furthermore, it requires a large number of beyond the Standard Model scalar fields to dilute the contribution of gravitons to the number of relativistic degrees of freedom during Big Bang Nucleosynthesis (BBN) [5] and may lead to other cosmological problems, as large isocurvature fluctuations or overproduction of dangerous relics [44]. Instant preheating appears a better suited solution to the "graceful exit" problem but its near-future detection appears unlikely.

In this work, we present a model where the inflaton accounts for the current DE density within the warm inflation paradigm [45, 46]. This shares various similarities with the scenario that we have proposed in [47] where the inflaton accounted for the all the dark matter in the Universe, within the Warm Little Inflaton model [48], which constitutes the simplest successful realization of warm inflation in quantum field theory. This scenario exhibits three main features: (i) the inflaton is generated by the collective spontaneous breaking of a U(1) symmetry [49], which, in combination with a discrete interchange symmetry, removes all corrections to the inflaton effective potential (both at zero and finite temperature), (ii) the inflaton interacts with other fields only during inflation, becoming stable and sterile afterwards and (iii) these interactions sustain a non-negligible (albeit sub-dominant) radiation bath during inflation through dissipative effects, leading to a smooth transition into a radiation-dominated era with no need for a reheating process [50]. Here we show that these same features can naturally lead to quintessential inflation.

Given the similarities of this model with that of [47], here we will only briefly discuss its main technical features and refer the interested reader to [47, 48] for more details. The inflaton is generated when two complex scalar fields, $\phi_{1}$ and $\phi_{2}$, with identical charges under a U(1) gauge symmetry, develop an identical non-zero vacuum expectation value, leading to collective spontaneous symmetry breaking parametrically close to the grand unification scale. While the overall phase of the fields is the Goldstone boson that becomes the longitudinal polarization of the massive U(1) vector field, the relative phase between the two complex fields, the only remaining light field is the inflaton $\phi$. The vacuum manifold can be written, in the unitary gauge, as:

[TABLE]

The complex scalars also interact with chiral fermion fields $\psi_{1L,R}$ and $\psi_{2L,R}$, and we impose a discrete interchange symmetry under which $\phi_{1}\leftrightarrow\phi_{2}$ and $\psi_{1L,R}\leftrightarrow\psi_{2L,R}$, as well as a sequestering of the "1" an "2" sectors (e.g. via additional global symmetries or localization in extra-dimensions). This yields a Lagrangian density:

[TABLE]

where in the last line we used the vacuum parametrization of Eq. (1) and $\psi_{1,2}$ denote the resulting Dirac fermions, of mass $gM$ (at zero temperature). These fermions then remain light provided that the temperature during inflation is $gM\lesssim T\lesssim M$. The fact that the physical mass of the fermions is independent of the inflaton field ensures that no quantum and thermal corrections can spoil the flatness of the inflaton potential, which holds both during and after inflation. This corresponds to a modification of the original Warm Little Inflaton scenario [48], where only the leading thermal corrections to the scalar potential were absent as a result of the discrete interchange symmetry.

Note that the interchange symmetry corresponds to a $\mathbb{Z}_{2}$ reflection of the inflaton, $\phi\leftrightarrow-\phi$, composed with an interchange $\psi_{1}\leftrightarrow\psi_{2}$, which prevents the inflaton from decaying into any field besides the $\psi_{1,2}$ fermions. When $T<gM$, these mediators become non-relativistic and decouple from the dynamics, such that the inflaton becomes stable at late times.

During inflation, however, $T>gM$, and the fermions can be produced through non-equilibrium dissipative effects sourced by the slow motion of the inflaton field. These arise from non-local corrections to the inflaton’s effective action, which are non-zero at finite temperature despite the absence of local corrections to the scalar potential. Physically, such dissipative effects are a consequence of the inflaton’s interaction with the fermions in the thermal bath, as characteristic of any dynamical system interacting with the surrounding environment.

We assume that the fermions $\psi_{1,2}$ are kept close to a local thermal equilibrium distribution, $\Gamma_{\psi}>H$, through Yukawa interactions with an additional light scalar field $\sigma$ and light fermion field $\psi_{\sigma}$ in the thermal bath, of the form $h\sigma\overline{\psi}_{1,2}\psi_{\sigma}+\text{H.c.}$111We do not consider direct couplings between the inflaton and the light fields since their renormalized values can be chosen to be arbitrarily small at the high energies relevant during inflation. After inflation, the fermions $\psi_{1,2}$ decouple from the effective field theory and hence do not regenerate such couplings through quantum corrections.. Dissipative effects then lead to an additional friction term $\Upsilon\dot{\phi}$ in the inflaton’s equation of motion, where the dissipation coefficient, $\Upsilon$, can be computed from first principles using quantum field theory at finite temperature222More details are available at https://journals.aps.org/prl/supplemental/

10.1103/PhysRevLett.122.161301/Sup$\_$Mat$\_$Letter.pdf [47],

[TABLE]

where $m_{\psi}=\sqrt{g^{2}M^{2}+h^{2}T^{2}/8}$ is the thermally corrected $\psi_{1,2}$ mass. We thus have $\Upsilon\propto T$ in the high-temperature regime relevant during inflation. The warm inflation dynamics is then described by the coupled differential equations for the inflaton and radiation energy density:

[TABLE]

The Hubble parameter is determined by the Friedmann equation $H^{2}=(\rho_{\phi}+\rho_{R})/(3M_{\text{P}}^{2})$, where $M_{\text{P}}^{2}=1/(8\pi G)$ is the reduced Planck mass, and the inflaton and radiation energy densities are given by $\rho_{\phi}\equiv\dot{\phi}^{2}/2+V(\phi)$ and $\rho_{R}=(\pi^{2}/30)g_{*}T^{4}$, with $g_{*}$ denoting the number of relativistic degrees of freedom (DOF). Note that, if $\Upsilon\rightarrow 0$, i.e. if interactions between $\phi$ and $\psi_{1,2}$ become negligible, whatever radiation density $\rho_{R}$ exists in the early Universe is diluted away by the inflationary expansion as in cold inflation models. This would then require an extra reheating mechanism at the end of inflation to recover the Hot Big Bang cosmology.

To completely determine the dynamics, one has to choose the inflaton potential $V(\phi)$. In the case where this is given by the simplest renormalizable potential [47], one obtains inflaton cold dark matter. To obtain inflaton dark energy, one has to choose a potential with no minima at finite field values that is flat enough at late times. As we will see, the required flatness of the potential depends on the initial conditions for the quintessence behavior, which are set at the end of the inflationary period. As a proof of concept, we adopt the same potential as [5]:

[TABLE]

where $\lambda\sim 10^{-15}$ is fixed by the CMB anisotropies power spectrum [51], $\alpha\geq 0$ is a tunable parameter that determines the steepness of the quintessence potential and $m$ is a mass scale, fixed by present dark energy density. For a dimensionless Hubble parameter $h\approx 0.68$ and $\Omega_{\Lambda}\approx 0.69$, consistent with [52], the example in Fig. 1, where $\alpha=1$, has $m\approx 3\text{ MeV}$. Although the $\mathbb{Z}_{2}$ symmetry is broken in the second branch of the potential (and in fact in any quintessence model with a finite field value), the effective mass of the inflaton field is $\lesssim H_{0}\sim 10^{-33}\text{ eV}$, and its interactions are mediated by heavy fields, $\psi_{1,2}$, so that $\phi$ becomes effectively stable and sterile at late times.

For large negative values of $\phi$, the quartic term dominates and slow-roll inflation occurs. For large positive values of $\phi$, the potential becomes an inverse power law, $V(\phi)\approx\lambda m^{4+\alpha}/\phi^{\alpha}$. Provided that $\phi$ is large enough after inflation, the inflaton will eventually reenter a slow-roll regime, causing a second stage of accelerated expansion.

A representative example of the cosmological evolution is illustrated in Fig. 1 from right to left. Since initially $K(\phi)\ll V(\phi)$, the slowly rolling inflaton field acts as an effective cosmological constant for at least $\sim 60$ e-folds, at the same time sustaining a non-negligible radiation density through the aforementioned dissipative effects, as displayed in Fig. 2: while $H$, $\Upsilon$ and $\phi$ are slowly varying, there is a non-zero steady-state solution for $\rho_{R}\approx\Upsilon\dot{\phi}^{2}/4H$. Dissipative effects become strong ($Q\equiv\Upsilon/(3H)\gtrsim 1$) at the end of inflation, allowing radiation to take over as the dominant component of the energy density of the Universe at a temperature $T_{R}\sim 10^{13}-10^{14}$ GeV, since

[TABLE]

where we used the steady-state solution above and the slow-roll equations $3H(1+Q)\dot{\phi}=-V^{\prime}(\phi)$, $H^{2}\approx V(\phi)/(3M_{\text{P}}^{2})$, with $\epsilon_{\phi}\equiv M_{\text{P}}^{2}(V^{\prime}/V)^{2}/2$. During inflation, $\epsilon_{\phi}\ll 1+Q$ and radiation is sub-dominant, while at the end of the slow-roll regime $\epsilon_{\phi}\sim 1+Q$ and radiation can smoothly take over if $Q\gtrsim 1$ [50]. This is in contrast with [5, 40, 41] where inflation was followed by a period of kination. Concomitantly, the $\dot{\rho}_{R}$ term in (4) becomes significant and makes the temperature drop below the mass of the $\psi_{1,2}$ fermions. Therefore, these decouple from the dynamics, making the inflaton effectively stable. Around the same time, $T/H$ becomes large enough to excite all the Standard Model DOF [53], resulting in a temperature drop by a factor $\sim$2 and further suppressing $Q$. This is represented by the first dashed line (from right to left) in the top panel of Fig. 1.

Since the dissipation coefficient $\Upsilon$, given in Eq. (3), becomes exponentially suppressed for $m_{\psi}>T$, the inflaton field will then quickly roll from $|\phi_{\text{eq}}|\lesssim M_{\text{P}}$ at inflaton-radiation equality towards $\phi=0$ in an approximately quartic potential (check first dashed line from right to left in the middle panel of Fig. 1). A precise description of this rolling phase requires a numerical solution of (4) (see, however, [54]333The parametric shape of $\phi(t)$ can be obtained by considering the inflaton field rolling down the quartic potential with negligible dissipation in a radiation-dominated universe. This is similar to the solution obtained in [54], where $\rho_{\phi}$ is dominant but its equation of state is close to that of relativistic matter.), but the important point is to determine the velocity at which the field crosses the origin, $\phi^{\prime}_{0}\equiv\left.\dot{\phi}/H\right|_{\phi=0}$, since this determines how far it can move in the second branch of the potential, as we will see below.

An instantaneous conversion of the inflaton’s potential energy at the end of inflation into kinetic energy at the origin would then yield $\phi^{\prime}_{0}\simeq\sqrt{3}M_{\text{P}}$. However, we find that the finite duration of this rolling period, where both Hubble and dissipative friction play a non-negligible role, yields somewhat smaller values. For instance, in the example shown in Fig. 1, we find $\phi^{\prime}_{0}\simeq 1.2M_{\text{P}}$. More generally, we find numerically that this quantity exhibits a mild dependence on the value of the dissipative ratio $Q_{*}$ when the relevant CMB scales become super-horizon approximately 60 e-folds before the end of inflation, yielding $\phi^{\prime}_{0}\simeq(0.5-1.5)M_{\text{P}}$ for $Q_{*}\sim 10^{-3}-10^{-1}$.

Although here we will not discuss in detail the spectrum of primordial density fluctuations in the Warm Little Inflaton scenario, referring instead the reader to the detailed analyses in [48, 55, 56] and references therein, we note that for a quartic potential consistency with CMB observations requires weak dissipation at horizon-crossing, in the above mentioned range. This ensures that the growth of thermal inflaton fluctuations during inflation is mild enough to yield a sufficiently red-tilted power spectrum, at the same time giving the necessary suppression of the tensor-to-scalar ratio through the thermal enhancement of scalar inflaton fluctuations when $T\gtrsim H$ at horizon-crossing. We emphasize that $Q$ is dynamical, with $Q\gtrsim 1$ at the end of inflation even if $Q_{*}\ll 1$, as illustrated in Fig. 2.

Once the field rolls beyond $\phi>0$, its potential energy becomes negligible and the equation of motion becomes:

[TABLE]

with primes denoting derivatives with respect to the number of e-folds of expansion, $dN_{e}=Hdt$. Since radiation becomes dominant at the end of the slow-roll regime, $H^{\prime}/H=-2$, so that Eq. (7) is $\phi^{\prime\prime}+\phi^{\prime}=0$ and admits the simple solution:

[TABLE]

where $N_{0}$ is the number of e-folds at which $\phi=0$. This implies that the field value quickly approaches $\phi_{\infty}=\phi^{\prime}_{0}$, as anticipated above. We note that, more accurately $H^{\prime}/H<-2$, since the inflaton’s kinetic energy is not completely negligible, yielding slightly larger asymptotic field values, for instance $\phi_{\infty}\simeq 1.3M_{\text{P}}$ in the example shown in Fig. 1. This nevertheless generically yields $\phi\sim M_{P}$ in this kinetic energy-dominated phase.

Since $H^{2}\propto a^{-4}$ and $\phi^{\prime}\propto\exp(-N_{e})\propto a^{-1}$, the inflaton behaves as stiff matter ($w_{\phi}=1$), with $\rho_{\phi}\simeq H^{2}\phi^{\prime 2}/2\propto a^{-6}$. The universe then undergoes a long radiation-dominated phase with a sub-leading inflaton stiff matter component (check $w_{\phi}$ in the bottom panel of Fig. 1). Eventually, the kinetic energy becomes diluted enough and the potential dominates the inflaton energy density after a number of e-folds since horizon-crossing $N_{\text{t}}\approx N_{e}+\log\left(K(0)/V(\phi_{\infty})\right)/6$, where $K(0)=K(\phi=0)$ and $V(\phi_{\infty})=\lambda m^{4+\alpha}/\phi_{\infty}^{\alpha}$. $N_{\text{t}}$, the total number of e-folds from horizon crossing, is $\sim 100$ ($\log_{10}(1+\text{z}_{\text{t}})\approx 10$), i.e. shortly after BBN. From this point onwards, $V(\phi)\gg K(\phi)$, $w_{\phi}\approx-1$ and the inflaton fluid is indistinguishable from a cosmological constant because the Hubble friction term, $H$, is much larger than the effective mass of the field $m_{\phi}^{2}=V^{\prime\prime}(\phi_{\infty})=\alpha(\alpha+1)V(\phi_{\infty})/\phi_{\infty}^{2}\sim\alpha(\alpha+1)H_{0}^{2}$. As a result, the field is effectively frozen, and only when $H\sim H_{0}$ can $\phi$ start to evolve significantly, as usual in thawing quintessence models444This occurs even though the inflaton potential is an inverse power law, usually associated with tracking freezing models of quintessence [57]. This shows the importance of initial conditions in quintessence models of dark energy. This behavior had been pointed out earlier in [5].. This means that the cosmological evolution of this model is somewhat independent of the particular quintessence potential, as long as it decreases rapidly once the field crosses the origin. It could, therefore, be replicated in more well-motivated potentials [58] and our analysis is not restricted to the particular Peebles-Vilenkin form.

Due to the dilution of radiation and matter, the inflaton field eventually becomes the dominant component, causing the present accelerated expansion. Since the slow-roll parameters $\epsilon_{\phi}=\alpha^{2}M_{\text{P}}^{2}/(2\phi^{2})$ and $\eta_{\phi}=\alpha(\alpha+1)M_{\text{P}}^{2}/\phi^{2}$ are non-negligible when the field becomes dynamical, $m_{\phi}\sim H$, the inflaton equation of state deviates from $w_{\phi}\approx-1$ due to an increase of $K(\phi)$:

[TABLE]

where we have expanded the denominator for $K(\phi)\ll V(\phi)$. It relaxes back to $w_{\phi}\approx-1$ as the field dominates the energy budget of the Universe, according to the de Sitter no-hair theorem [59], as the field moves towards larger values where $\epsilon_{\phi}\ll 1$.

Choosing different values of $\alpha$ or $Q_{*}$ (and therefore $\phi^{\prime}_{0}$) will result in different late time behaviors, as represented in Fig. 3. The $(w_{0},w_{a})$ parameterization [60], $w_{\phi}=w_{0}+w_{a}(1-a)$ with $a=1/(1-z)$, conveniently encapsulates the variation of $w_{\phi}$, yielding in this example $(-0.73,-0.33)$, $(-0.85,-0.20)$, $(-0.93,-0.09)$ for $\alpha=2,~{}3/2,~{}1$, respectively555It was mentioned in [60] that a large range of thawing models is well fitted by $w_{a}\approx-1.58(1+w_{0})$. We observe a slightly smaller slope, $w_{a}\approx-1.26(1+w_{0})$, for the examples of Fig. 3. . As expected, smaller values of $\alpha$ correspond to shallower quintessence potentials (smaller $\epsilon_{\phi}$, $\eta_{\phi}$), which lead to smaller deviations from $w_{\phi}\approx-1$ at late times. This has some degeneracy with the dissipation at horizon crossing, $Q_{*}$, as larger values of this quantity lead to a larger $\phi^{\prime}_{0}$ and consequently a value of $w_{\phi}$ closer to a cosmological constant. For instance, choosing $Q_{*}=0.07$, we find $\phi^{\prime}_{0}\approx 1.3M_{\text{P}}$ and so $(w_{0},w_{a})$ becomes $(-0.77,-0.29)$, $(-0.87,-0.17)$ , $(-0.94,-0.08)$ for $\alpha=2,~{}3/2,~{}1$. Comparing the above results with bounds on $(w_{0},w_{a})$ from the Planck collaboration [61, 52] suggests that $\alpha=2$ ($\alpha=3/2$) is (marginally) outside the $68\%$ confidence interval for the Planck+SNe+BAO data, with $\alpha=1$ inside it, but these conclusions should be taken at face value, given the particular details of the model and the data-set dependence of the $(w_{0},w_{a})$ bounds (see e.g. [62]).

In summary, in this Letter, we have shown that the inflaton can account for the present accelerated expansion of the Universe in the context of a concrete quantum field theory model, the Warm Little Inflaton scenario [48], as a result of the smooth transition between a period of warm inflation and a radiation-dominated era and the fact that the symmetries of the model make the inflaton both stable and sterile at late times. Alongside our earlier work on inflaton-dark matter [47] and previous works on baryogenesis/mattergenesis during warm inflation [63, 64, 65], this shows that this alternative paradigm may have a very significant impact not only on the inflationary dynamics itself but also on the other major open problems in cosmology.

The dynamics provides several interesting observational features: warm inflation generically predicts a suppressed tensor-to-scalar ratio and a modified consistency relation between the latter and the tensor spectral index [66, 67, 68], which will be further probed with next generation CMB experiments [69, 70]. Furthermore, the current data already constrains the consistent $(w_{0},w_{a})$ space of models of quintessence severely, causing them to be very close to a cosmological constant (this depends somewhat on the analysis, see [62, 71]). This causes smaller values of $\alpha$ to be favored in the particular scenario displayed here. Note, however, that this conclusion is quite model-dependent, as scenarios which predict larger values of $Q_{*}$ will consequently have larger asymptotic field values, allowing for a larger range for the consistent values of $\alpha$.

We also emphasize that our main results are not restricted to the particular Peebles-Vilenkin form of the quintessence part of the potential nor on how it is "glued" to the quartic inflationary potential, since the dynamics is insensitive to the potential function for $0<\phi\lesssim M_{\text{P}}$, and the exponent $\alpha$ simply characterizes the slope of the potential at $\phi\sim\phi_{\infty}\sim M_{\text{P}}$.

Finally, a potentially distinctive feature of our scenario are dark energy isocurvature fluctuations, which would be (anti-)correlated with the main adiabatic curvature perturbations given their common origin in inflaton fluctuations. The thermal nature of such fluctuations, with a distinctive spectrum dependent on the form of the dissipation coefficient during inflation ($\Upsilon\propto T$ in this case), could also in principle be used to probe the warm vs. cold nature of the quintessential inflation scenario. Detecting dark energy isocurvature modes may, however, be quite hard in general, given that they only affect the CMB spectrum on large scales through the integrated Sachs-Wolfe effect at late times (see e.g. [72, 73, 74]), but we hope that the present work motivates further exploration of this possibility.

Acknowledgements

L.B.V. is supported by FCT grant PD/BD/140917/2019. J.G.R. is supported by the FCT Investigator Grant No. IF/01597/ 2015. This work was partially supported by the H2020-MSCA-RISE-2015 Grant No. StronGrHEP-690904 and by the CIDMA Project No. UID/MAT/04106/2019. L.B.V. thanks Tiago Barreiro and Juan P. Beltrán Almeida for private correspondence.

During the completion of this work, we became aware of the recent analysis of warm quintessential inflation in [75], which shares some similarities with this Letter, albeit using a different but nevertheless complementary approach. Here we consider a concrete quantum field theory model of inflation, which allows us to follow in detail the dynamical evolution of dissipative effects. This affects not only the duration of inflation and associated observational predictions but also the quintessence phase at late times. In particular, in our scenario the non-negligible dissipation at the end of inflation prevents an intermediate kination era, so that the field cannot evolve towards large super-Planckian values, thus requiring shallower quintessence potentials.