# Extended Warm Inflation

## I Introduction

Since the advent of inflationary models, the roles played by scalar fields at different epochs of the cosmic history have been extensively investigated. Such fields have been invoked in a variety of disparate scenarios with rather different goals. Some known examples are: (i) the inflaton, the field that drives inflation , (ii) the axion, a cold dark matter candidate  and (iii) the dilaton, the field appearing in the low energy string action  which addresses the same issues of inflation and may provide a solution to the singularity problem . More recently, inspired by the existing observational data and theoretical speculations, some authors have also suggested scalar fields (sometimes called “quintessence”) as the sough non-baryonic dark matter. These “remnant” fields, might have important consequences on the formation of the large scale structure, as well as be responsible by the present day accelerated phase , as indicated by the latest type Ia Supernovae observations .

Despite their generalized use in the cosmological framework, the physical situations in which the scalar fields are commonly considered are rather particular. For example, in the new inflation case there is not a fundamental justification for “turning off” all the possible couplings during the slow rollover phase and doing the opposite just at the onset of the thermalisation phase of reheating (as well as afterwards, if some potential energy is still available). This can be achieved only by considering very strict initial conditions, which weakened the viability of such a scenario. But this is not an isolated case. Neglecting possible thermal couplings between scalar fields and the other constituents of the universe is a general feature assumed in nearly all scalar field models presented in the literature.

One important exception is the explosive reheating period required to succeed any version of isentropic inflation. Along this process, as the field coherently oscillates about the minimum of the potential, its energy is drained to the matter and radiation components. Either in its earlier  or in its modern version based on parametric resonance (sometimes called preheating), the reheating mechanism is a relatively fast process, and virtually all the entropy of the present universe may be generated in this way. This is certainly not the most general case. In principle, a permanent or temporary coupling of the scalar field $`\varphi `$ with other fields might also lead to dissipative processes producing entropy at different eras of the cosmic evolution. It is expected that progresses in non-equilibrium statistics of quantum fields will provide the necessary theoretical framework for discussing dissipation in more general cases (see for example  and references therein). Another possibility is the so-called “instant preheating” . In this process, the inflaton decays continuously into another scalar field as it rolls down the potential. This second field is very short-lived and rapidly decays into fermions thus furnishing a sustained entropy generation, including for quintessence-like models.

Although a justification from first principles for dissipative effects has not been firmly achieved, such effects should not be ruled out only on readiness basis. Much work can be done in phenomenological grounds as, for instance, by applying nonequilibrium thermodynamic techniques to the problem or even studying particular models with dissipation. An interesting example of the latter case is the warm inflationary picture recently proposed . Like in new inflation, a phase transition driving the universe to an inflationary period dominated by the scalar field potential is assumed. However, a standard phenomenological friction-like term $`\mathrm{\Gamma }\dot{\varphi }^2`$ is inserted into the scalar field equation of motion to represent a continuous energy transfer from $`\varphi `$ to the radiation field. This persistent thermal contact during inflation is so finely adjusted that the scalar field evolves all the time in a damped regime generating an isothermal expansion. As a consequence, the subsequent reheating mechanism is not needed and thermal fluctuations produce the primordial spectrum of density perturbations (see also reference ).

Warm inflation was originally formulated in a phenomenological setting, but some attempts of a fundamental justification has also been presented . Furthermore, a dynamical systems analysis  showed that a smooth transition from inflationary to a radiation phase is attained for many values of the friction parameter, thereby showing that the warm scenario may be a workable variant to inflation. As it appears, its unique negative aspect is closely related to a possible thermodynamic fine-tunning, because an isothermal evolution of the radiative component is assumed from the very beginning in some versions of warm inflation (for comments on this issue, see ). In other words, the thermal coupling acting during inflation is so powerful and finely adjusted that the scalar field decays ensuring a constant temperature even considering the exponential expansion of the universe.

In brief, the aim of this paper is to relax this hyphothesis. However, instead of proposing another particular inflationary model, we discuss how the differences between the isentropic and the isothermal inflationary scenarios can be depicted in a convenient parameter space. As we shall see, these models are only two extreme cases of an infinite two-parametric family. Hopefully, this unified view may indicate ways to a consistent phenomenological treatments of these models based on the methods of nonequilibrium thermodynamics. We also discuss how the standard slow roll conditions are modified due to the scalar field decay.

## II Scalar Field with Dissipation

We will limit our analysis to homogeneous and isotropic universes, described by the flat Friedmann-Robertson-Walker (FRW) line element

$$ds^2=dt^2a^2(t)\left(dr^2+r^2d\theta ^2+r^2\mathrm{sin}^2\theta d\varphi ^2\right),$$
(1)   
where $`a(t)`$ is the scale factor (in our units $`\mathrm{}=c=1`$). The source of this spacetime is a mixture of a real and minimally coupled scalar field interchanging energy with a perfect fluid representing all the other fields. The Lagragian density for the scalar field is

$$=\frac{1}{2}^\mu \varphi _\mu \varphi V(\varphi )+_{\mathrm{int}},$$
(2)   
where the interaction is implied by the term $`_{\mathrm{int}}`$ and $`V(\varphi )`$ is the scalar field potential. This field has the stress-energy tensor given by

$$T_\varphi ^{\mu \nu }=^\mu \varphi ^\nu \varphi g^{\mu \nu }.$$
(3)   
The other component of the mixture is a simple fluid with energy-stress tensor

$$T_m^{\mu \nu }=(\rho +p)u^\mu u^\nu pg^{\mu \nu },$$
(4)   
where energy density and pressure are given respectively by $`\rho `$ and $`p`$. The total energy stress tensor of the system $`T_t^{\mu \nu }=T_\varphi ^{\mu \nu }+T_m^{\mu \nu }`$ obeys Einstein’s field equations, from which we obtain the equations of motion

$$3H^2=\frac{8\pi }{m_{\mathrm{pl}}^2}\left(\frac{\dot{\varphi }^2}{2}+V(\varphi )+\rho \right),$$
(5)   
$$3H^2+2\dot{H}=\frac{8\pi }{mpl^2}\left(\frac{\dot{\varphi }^2}{2}V(\varphi )+p\right),$$
(6)   
where a dot means time derivative, $`H=\dot{a}/a`$ is the Hubble parameter, $`m_{\mathrm{pl}}^2=1/G`$ is the Planck mass, and we have used that the scalar field energy density and pressure are, respectively

$`\rho _\varphi ={\displaystyle \frac{1}{2}}\dot{\varphi }^2+V(\varphi )`$ (7)   
$`p_\varphi ={\displaystyle \frac{1}{2}}\dot{\varphi }^2V(\varphi ).`$ (8)   
Now, assuming that the perfect fluid complies with the “gamma law” equation of state,

$$p=(\gamma 1)\rho ,$$
(9)   
the energy conservation law for this interacting selfgravitating mixture can be cast in the form

$$\dot{\varphi }(\ddot{\varphi }+3H\dot{\varphi }+V^{}(\varphi ))=\dot{\rho }3\gamma H\rho .$$
(10)   
where the prime denotes derivative with respect to the field $`\varphi `$. In order to decouple the two sides of the above equation and allow for the decay of the scalar field into ordinary matter and radiation, we will assume the usual phenomenological friction-like term (3$`\mathrm{\Gamma }\dot{\varphi }^2`$), thereby spliting (10) in two equations

$$\ddot{\varphi }+3H\dot{\varphi }+V^{}(\varphi )=3\mathrm{\Gamma }\dot{\varphi }$$
(11)   
and

$$\dot{\rho }+3\gamma H\rho =3\mathrm{\Gamma }\dot{\varphi }^2,$$
(12)   
where the dissipative coefficient $`\mathrm{\Gamma }`$ is the decay width of the scalar field and the factor 3 has been introduced for mathematical convenience. For models endowed with an isentropic inflationary stage, we know that $`\rho 0`$, at the onset of the coherent oscillations, that is, when the field $`\varphi `$ oscillates about the minimum of its potential $`V(\varphi )`$. In these cases, the term $`\mathrm{\Gamma }\dot{\varphi }^2`$ is generally inefficient for describing the first stages of the reheating process  and is not acting during exponential inflation. However, in the presence of a nonnegligible thermal component, a friction term can be justified under special conditions during a de Sitter regime . Aplications of this theory lead to complicated but promising warm inflationary models, with an unreasonably large quantity of light fields coupled to the inflaton , although some of these models might be physically interpreted in terms of string theory.

In what follows, it will be assumed the validity of this friction term as representing the thermal contact between $`\varphi `$ and the other fields for any epoch. We also consider the “thermal decay width” $`\mathrm{\Gamma }`$ as a generic function of the temperature, or equivalently, of the cosmic time.

Let us now define the following dimensionless parameters:

$$x\frac{\dot{\varphi }^2}{\dot{\varphi }^2+\gamma \rho },$$
(13)   
and

$$\alpha \frac{3H\dot{\varphi }^2}{3H\dot{\varphi }^2+|\dot{\rho }|}.$$
(14)   
The adoption of these parameters is somewhat natural and may be understood as follows. The warm picture departs quantitatively from new inflation because during the inflationary stage the energy density of the material component is not negligible. More precisely, it is not negligible only in comparison to the scalar field kinetic term since the potential $`V(\varphi )`$ must be strongly dominant in order to generate exponential expansion. This means that one needs to compare the kinetic term $`\dot{\varphi }^2=\rho _\varphi +p_\varphi `$ with a quantity involving the energy of the material component written in a convenient form. For the mixture, the (total) equivalent quantity is $`\rho _t+p_t=\dot{\varphi }^2+\gamma \rho `$. In this way, recalling the energy conservation law, $`x`$ is defined by the ratio between two “inertial mass” terms which redshift away due to the universe expansion. Additionally, it is not enough to compare the values of $`\dot{\varphi }^2`$ and $`\gamma \rho `$ since we also need to quantify how they evolve in the course of the expansion. This explains the introduction of the parameter $`\alpha `$ involving $`3H\dot{\varphi }^2`$ and $`\dot{\rho }`$. The presence of $`|\dot{\rho }|`$ in this parameter is also reasonable. It comes into play because we are adopting the original isothermal scenario proposed by Berera as a limiting case. Indeed, it seems to be an extreme theoretical situation where the cooling rate of the radiation due to expansion is fully compensated by the transference of particles from the scalar field to the ordinary constituents. It is implicit in the definition of $`\alpha `$ that $`\dot{\rho }0`$, since a reheating phase is unecessary in warm inflation. However we should point out that a positive $`\dot{\rho }`$ is possible and has been investigated .

The convenience of these parameters is apparent since by construction they are dimensionless and constrained on the intervals $`0x1`$ and $`0\alpha 1`$. In particular, for the isothermal inflation  we have $`\alpha =1`$ and $`x0`$, because $`\dot{\rho }=0`$ and $`\gamma \rho \dot{\varphi }^2`$. For isentropic inflation one has $`x=1`$ and $`\alpha =1`$ since $`\gamma \rho \dot{\varphi }^2`$ and $`\dot{\rho }3H\dot{\varphi }^2`$ (the radiation becomes exponentially negligible). Similarly, noninteracting quintessence-like models lie at some intermediary value of $`\alpha =x`$ between 0 and 1. When both parameters are equal to zero, we have the standard model (with a possible cosmological constant). Therefore, the most common solutions, with or without thermal couplings, can be portrayed in this bidimensional parameter space (see Fig. 1).

## III Inflation, Dissipation and Slow Roll Conditions

A fundamental ingredient of many inflationary variants is the period of “slow roll” evolution of the scalar field during which the inflaton field evolves so slowly that its kinetic term remains always much smaller than $`V(\varphi )`$. A full description of the two-component mixture requires the introduction of a dynamical parameter which is function of the total energy density of the mixture. In order to discuss slow roll dynamical conditions, a suitable choice of such a paramater is $`\gamma _{\mathrm{eff}}`$, derived from manipulations on Eqs. (5) and (6)

$$\gamma _{\mathrm{eff}}=\frac{2\dot{H}}{3H^2}=\frac{\dot{\varphi }^2+\gamma \rho }{\rho _t}.$$
(15)   
Slow roll conditions are imposed to assure (nearly) de Sitter solutions for an amount of time, say, $`\mathrm{\Delta }t_I`$, which must be long enough to solve the problems of the hot big bang. By inspection of Eqs. (5), (6) and (15) one obtains that the condition for having a de Sitter universe is

$$3H^22|\dot{H}|\gamma _{\mathrm{eff}}0V(\varphi )\rho +\frac{\dot{\varphi }^2}{2}.$$
(16)   
Additionally, a sufficient amount of inflation requires the above relations to be valid along the interval $`\mathrm{\Delta }t_I`$. This usually implies constraints on the slope and curvature of $`V(\varphi )`$. However, as we shall see, the coupling between $`\varphi `$ and the $`\gamma `$-fluid during inflation may relax these constraints for a large range of intermediary situations between the isentropic and warm scenarios. In our case, these conditions are represented by (see also)

$$H^2\frac{8\pi }{3m_{\mathrm{pl}}^2}V(\varphi ),$$
(17)   
$$3(H+\mathrm{\Gamma })\dot{\varphi }V^{}(\varphi ).$$
(18)   
Notice that the discussion that follows is independent of constraints on the signal of $`\dot{\rho }`$, so that we are not limited to the definition of $`\alpha `$. The approximated expressions (17) and (18) give us the following constraints on the shape of $`V(\varphi )`$

$$\left(\frac{V^{}(\varphi )}{V(\varphi )}\right)^2\frac{m_{\mathrm{pl}}^2}{16\pi }\left(1+\frac{\mathrm{\Gamma }}{H}\right)^2\frac{\dot{\varphi }^2}{2V(\varphi )}$$
(19)   
and

$$\frac{V^{\prime \prime }(\varphi )}{V(\varphi )}\frac{m_{\mathrm{pl}}^2}{24\pi }\left(1+\frac{\mathrm{\Gamma }}{H}\right)\frac{\dot{\varphi }^2+\gamma \rho }{2V(\varphi )}\frac{1}{3H}\frac{d}{dt}\left(\frac{\mathrm{\Gamma }}{H}\right),$$
(20)   
We recall that the standard slow roll conditions appearing in isentropic inflation are recovered for $`\mathrm{\Gamma }=0`$ (decoupled mixture). They imply that inflation is possible only if the potential is extremely flat, by the last inequalty in (16). However, in the limit ($`\mathrm{\Gamma }H`$), one may see that the first and the second derivative of the potential are not necessarily small, in order to guarantee the continued accelerated expansion. Thus, the extremely flat potential of usual inflation assuring the slow down of the scalar field by an enough amount of time, may be replaced by a large friction-like term with no extreme behavior of the derivatives of $`V(\varphi )`$. The field does not accelerate, thereby $`\dot{\varphi }^2`$ becoming comparable to $`V(\varphi )`$, because the friction may be really large and not because the potential is unusually flat. This explains why this class of extended scenarios may provide a solution for the fine tunning problems plaguing the new inflationary picture.

## IV A Toy Model

The standard procedure in scalar field cosmological model building is to assume a specific potential $`V(\varphi )`$, motivated or not from particle physics. This potential is used to solve the $`\varphi `$ equation of motion (11). Subsequently, if there is a nonegligible energy density stored in the other fields, the solution of $`\varphi `$ is inserted into (12), which is solved for the energy density $`\rho `$. As widely known, even for a given coupling, many solutions are possible just changing the potential $`V(\varphi )`$. To fix terminology, this method for generating solutions will be called dynamic approach.

On the other hand, the parameter space ($`\alpha ,x`$) can also be used as a guide in the search for new models. To show that a different (thermodynamic) route is also possible we first rewrite (12) in terms of $`x`$ and $`\alpha `$. We have

$$\frac{\mathrm{\Gamma }}{H}=\frac{1}{x}\frac{1}{\alpha }.$$
(21)   
It is worth noticing that (21) does not include the potential since it is representative only of the possible couplings, but not of the dynamics. In what follows we obtain a simple model example based on equations (13), (14) and (21). In principle, a more quantitative analysis will clarify the interesting features contained in the parameter space ($`\alpha ,x`$). In order to be as generic as possible, we also include the possibility of a decaying $`\varphi `$ during any era.

Since $`x`$ and $`\alpha `$ do not depend on $`V(\varphi )`$, one may show that equation (12) for the material medium can be rewritten as

$$\dot{\rho }=3\gamma H\rho \left(1\delta \right),$$
(22)   
where for short we have introduced the quantity

$$\delta \frac{\alpha x}{\alpha (1x)}.$$
(23)   
It is apparent that the simplest toy model is the one where the parameters $`\alpha `$ and $`x`$ are constants. In this case, the solution of (22) is

$$\rho (a)=\rho _\mathrm{I}\left(\frac{a}{a_\mathrm{I}}\right)^{3\gamma (1\delta )},$$
(24)   
where the integration constants, $`\rho _\mathrm{I}`$ and $`a_\mathrm{I}`$, are “initial” conditions just at the onset of the inflationary stage. It should be noticed that if $`\delta =1`$, or equivalently, if $`\alpha =1`$, the energy density $`\rho `$ remains constant ($`\rho =\rho _\mathrm{I}`$). In particular, if $`x0`$, that is, if $`\rho \dot{\varphi }^2`$, we recover the isothermal inflationary scenario proposed by Berera (see also Fig.1). This show more clearly why models with values of $`0<x<1`$ may also evolve isothermally during inflation (exponential or power-law). The other extreme situation is obtained if $`\delta =0`$, that is, if $`\alpha =x`$, with the energy density scaling as $`\rho a^{3\gamma }`$. This behavior is typical for an adiabatic expansion (very weak coupling), and in Fig. 1, it corresponds to the adiabatic line. Note that even in these circumstances, where the parameters are constants, one may expect that an intermediary situation ($`0<\delta <1`$) be physically more probable.

Under the above conditions, the energy density of the the scalar field may also be readily obtained. First we rewrite (11) as

$$\dot{\rho }_\varphi =3H\dot{\varphi }^2\left(1+\frac{\mathrm{\Gamma }}{H}\right).$$
(25)   
Now, since $`x=\mathrm{const}`$ we see from (13) that $`\dot{\varphi }^2=\frac{x}{1x}\gamma \rho `$, and using (21) and (24), this equation can be integrated in terms of the scale factor. The result is

$$\rho _\varphi =\rho _{\varphi _I}+\rho _I\left(\delta +\frac{x}{1x}\right)\left[\frac{\left(\frac{a}{a_I}\right)^{3\gamma (1\delta )}1}{1\delta }\right],$$
(26)   
where the constant $`\rho _{\varphi _I}`$ is the scalar field energy density when $`a=a_I`$.

In the adiabatic limit ($`\delta =0`$) the above expression reduces to

$$\rho _\varphi =\rho _{\varphi _I}+\rho _I\frac{x}{1x}\left[\left(\frac{a}{a_I}\right)^{3\gamma }1\right],$$
(27)   
The isothermal case ($`\delta =1`$) is also readily obtained using that $`lim_{q1}\frac{f^{1q}1}{1q}=\mathrm{ln}f`$. One finds

$$\rho _\varphi =\rho _{\varphi _I}\frac{3\gamma }{1x}\rho _I\mathrm{ln}\left(\frac{a}{a_\mathrm{I}}\right)$$
(28)   
If $`\delta 1`$, using (24) and (26), we have that the total energy density may be written as

$$\rho _t=\rho _{\varphi _I}B\rho _I+B\rho $$
(29)   
where $`B=\frac{x(1+\frac{\mathrm{\Gamma }}{H})}{(1x)(1\delta )}`$ is a constant.

The time dependence of the scale factor may also be obtained from (13) and (15), which give $`\gamma _{\mathrm{eff}}=\frac{\gamma \rho }{(1x)\rho _t}`$. We see that the possible evolution laws depend critically on the initial conditions. In particular, if $`B\frac{\rho _{\varphi _I}}{\rho _I}`$ we have $`\gamma _{\mathrm{eff}}\frac{\gamma \rho _I}{(1x)\rho _{\varphi _I}}`$. Thus, if $`\rho _{\varphi _I}\rho _I`$ we have exponential inflation, and if $`\rho _\varphi `$ dominates only moderately, the scale factor will evolve as a power law inflation with a coupled thermal component. As one may check, in this case the potential $`V(\varphi )`$ scales as $`e^{\lambda \varphi }`$, where $`\lambda (\delta )`$ is a positive parameter. Note still that if the first two terms in (29) do not cancel each other, we may have a “remnant” cosmological constant for large values of the cosmological time.

As we have seen, this simple model is representative of two different ways of extending the original warm inflation. One is keeping the isothermal condition ($`\alpha =\delta =1`$) for generic values of $`x`$. Another approach is to consider nonvanishing parameters $`\alpha `$ and $`x`$ ($`\alpha x`$), which allows for warm power law inflationary models. Hopefully, in a more general framework, where the pair of parameters ($`\alpha ,x`$) are time dependent functions, a consistent unified picture containing the new and warm inflation, as well as all the intermediary situations may be obtained.

## V Conclusions

Our analysis might be useful as an heuristic tool for building models with scalar fields coupled to a material medium. In principle, it can be applied to warm inflation-like models (with or without phase transitions involved), reheating (even for the old scenario) and quintessence-like models. As a rule, it could be tried as a first step, before attempting different potentials for the scalar fields. As we have shown, this happens because depending on the two parameters, the conditions for distinct dynamics of the scale factor (including inflationary regimes) are relatively independent of the shape of the potencial, since a large dissipative term may provide a natural slow rolling for the field. This may unconstrain the conditions on the flatness of the potential in such a way that even exotic models like the oscillating inflation of Damour and Mukhanov may provide the necessary number of e-foldings and enough post-inflationary radiation temperature. The toy model presented in the last section somewhat suggest that any inflationary interpolating solution between the isentropic and isothermal limits can be represented in the bidimensional parameter space ($`\alpha ,x`$). Particular examples will be discussed elsewhere .

As it appears, a more comprehensive phenomenological treatment of this matter should necessarily include thermodynamical constraints thus requiring the methods and techniques from nonequilibrium thermodynamics. As shown recently, a phenomenological coupling term explicitly dependent on the created particles (and not only on the scalar field) should be a natural outcome of these methods when the decay products thermalize with the heat bath. Such an approach may have interesting consequences on the old reheating and warm inflation models.

Acknowledgements The authors are grateful for the support of CAPES and CNPq (Brazilian research agencies). One of us (JASL) was also supported by the project PRONEX/FINEP (No. 41.96.0908.00).