Thermal history after inflation is studied in a chaotic inflation model with supersymmetric couplings of the inflaton to matter fields. Time evolution equation is solved in a formalism that incorporates both the back reaction of particle production and the cosmological expansion. The effect of the parametric resonance gives rise to a rapid initial phase of the inflaton decay followed by a slow stage of the Born term decay. Thermalization takes place immediately after the first explosive stage for a medium strength of the coupling among created particles. As an application we calculate time evolution of the gravitino abundance that is produced by ordinary particles directly created from the inflaton decay, which typically results in much more enhanced yield than what a naive estimate based on the Born term would suggest.

Large amplitude oscillation of cosmic field that may occur right after inflation and in the decay process of weakly interacting fields gives rise to violent particle production via the parametric resonance. In the large amplitude limit the problem of back reaction against the field oscillation is solved and the energy spectrum of created particles is determined in a semi-classical approximation. For large enough coupling or large enough amplitude the resulting energy spectrum is broadly distributed, implying larger production of high energy particles than what a simple estimate of the reheating temperature due to the Born formula would suggest.