Equation of state and beginning of thermalization after preheating

We study the out-of-equilibrium nonlinear dynamics of fields after post-inflationary preheating. During preheating, the energy in the homogeneous inflaton is exponentially rapidly transfered into highly occupied out-of-equilibrium inhomogeneous modes, which subsequently evolve towards equilibrium. The infrared modes excited during preheating evolve towards a saturated distribution long before thermalization completes. We compute the equation of state during and immediately after preheating. It rapidly evolves towards radiation domination long before the actual thermal equilibrium is established. The exact time of this transition is a nonmonotonic function of the coupling between the inflaton and the decay products, and it varies only very weakly (around 10(-35) s) as this coupling changes over several orders of magnitude. This result is applied to refine the relation between the number of e-foldings N and the physical wavelength of perturbations generated during inflation. We also discuss the implications for the theory of modulated perturbations from preheating. We finally argue that many questions of the thermal history of the universe should be addressed in terms of prethermalization, illustrating this point with a calculation of perturbative production of gravitinos immediately after chaotic inflation. We also highlight the effects of three-legs inflaton interactions on the dynamics of preheating and thermalization in an expanding universe.