Multi-species thermalization cascade of energetic particles in the early universe

Heavy long-lived particles are abundant in BSM physics and will, under generic circumstances, get to dominate the energy density of the universe. The resulting matter dominated era has to end through the decay of the heavy matter component particles of mass M, and before the universe cools down to a temperature of a few MeV, in order not to upset the successful predictions of Big Bang Nucleosynthesis. The process of thermalization of their decay products, of initial energy & SIM;M/2, into a thermal bath of temperature T primarily involves near-collinear splittings of energetic particles into two particles with lower energy. The correct treatment of these processes requires the inclusion of coherence effects which suppress the splitting rate. We write down and numerically solve the resulting coupled Boltzmann equations including all gauge bosons and fermions of the unbroken phase of the Standard Model (SM). We then comment on the dependence of the nonthermal spectra on the ratio M/T, as well as on the matter decay rate and branching ratios into various SM particles.