Perturbative reheating and thermalization of pure Yang-Mills plasma

We investigate the thermalization of high-energy particles injected from the perturbative decay of inflaton during the pre-thermal phase of reheating in detail. In general, thermalization takes a relatively long time in a low-temperature plasma; therefore, the instantaneous thermalization approximation is not justified, even for the reheating of the Standard Model (SM) sector. We consider a pure Yang-Mills (YM) theory as an approximation of the SM sector or a possible dark sector, considering the Landau-Pomeranchuk-Migdal effect, a quantum interference effect in a finite temperature plasma. We perform the first numerical calculation to solve the time evolution of the system, including the redshift due to the expansion of the Universe, and show the details of the temperature evolution near the maximum and the behavior of the quasi-attractors at later times. The maximal temperature T max and time scale t max are determined quantitatively, such as T max similar or equal to 0.05 x Gamma I M PI 2 / m I 3 2 / 5 m I \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ {\left({\Gamma}_I{M}_{\textrm{PI}}<^>2/{m}_I<^>3\right)}<^>{2/5}{m}_I $$\end{document} and t max similar or equal to 2 x 103 x Gamma I M PI 2 / m I 3 - 3 / 5 m I - 1 \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ {\left({\Gamma}_I{M}_{\textrm{PI}}<^>2/{m}_I<^>3\right)}<^>{-3/5}{m}_I<^>{-1} $$\end{document} in the SM-like system, where m I and Gamma I are the mass and decay rate of inflaton. We also provide a similar formula for pure SU(N) and SO(N) YM theories for general values of N and coupling constant alpha, including T max proportional to alpha 4/5 and t max proportional to N -2 alpha -16/5 behaviors and their numerical coefficients. The thermalization occurs in a finite time scale, resulting in a lower maximal temperature of the Universe after inflation than that under the instantaneous thermalization approximation.