Efficient numerical integration of thermal interaction rates

In many problems in particle cosmology, interaction rates are dominated by ${2}\leftrightarrow{2}$ scatterings, or get a substantial contribution from them, given that ${1}\leftrightarrow{2}$ and ${1}\leftrightarrow{3}$ reactions are phase-space suppressed. We describe an algorithm to represent, regularize, and evaluate a class of thermal ${2}\leftrightarrow{2}$ and ${1}\leftrightarrow{3}$ interaction rates for general momenta, masses, chemical potentials, and helicity projections. A key ingredient is an automated inclusion of virtual corrections to ${1}\leftrightarrow{2}$ scatterings, which eliminate logarithmic and double-logarithmic IR divergences from the real ${2}\leftrightarrow{2}$ and ${1}\leftrightarrow{3}$ processes. We also review thermal and chemical potential induced contributions that require resummation if plasma particles are ultrarelativistic.