Thermal field theory is applied to particle production rates in inflationary models, leading to new results for catalysed or two-stage decay, where massive fields act as decay channels for the production of light fields. A numerical investigation of the Boltzmann equation in an expanding universe shows that the particle distributions produced during small amplitude inflaton oscillations or even alongside slowly moving inflaton fields can thermalise.
I. INTRODUCTIONInflationary models give a picture of the early universe that has shown spectacular agreement with observation [1,2,3]. All inflationary models require a mechanism for reheating the universe to take it from the vacuum dominated inflationary phase to the hot radiation-dominated universe which we know must follow. The details of particular reheating mechanisms depend on the interactions between the inflaton and other fields, but the underlying process is particle production which fills the universe with radiation.The original work on reheating in the 1980's introduced particle production in an ad hoc fashion, assuming the rate of particle production by the inflaton φ in the limit of smallφ was proportional toφ n [4, 5, 6]. Most authors settled for n = 2, which has the advantage of being equivalent to a simple friction term in the inflaton equation of motion. At around the same time, a less ad hoc approach was based on particle production caused by an inflaton field oscillating about the minimum of the inflaton potential after the end of inflation [7,8]. This latter approach appeared to be the more consistent, and it is still widely used today.Following this early work, there where several attempts to apply new ideas in thermal field theory to the reheating process. These usually focused on finding effective field equations for the inflaton. The small-φ equation of motion was derived first, using linear response theory [9,10]. This has been used in the theory of warm inflation [11,12,13]. More general forms of the inflaton field equation, which where not limited to small time derivatives and could be applied to the oscillating inflaton, followed later [14,15,16].Renewed interest in particle production was ignited by the discovery of preheating, a nonperturbative process of inflaton decay though parametric resonance [17,18,19,20,21]. Like the earlier work, preheating involves particle production from an oscillating inflaton field. Preheating is a result of very large amplitude oscillations in the inflaton field. Large amplitude oscillations are a feature, though not necessarily a desirable one, of most single-field inflationary models.Preheating may or may not occur depending on the details of the particle model. In this paper we shall focus on models with a mechanism which we call catalysed or two-stage decay [22]. Unlike in the case of pre-heating, we examine what happens when the fields coupled to the inflaton are too massive to be produced directly. Instead, the fields can act as decay channels for the production of light fields. This requires the kind ...