A bstract In many problems in particle cosmology, interaction rates are dominated by 2 ? scatterings, or get a substantial contribution from them, given that 1 and 3 reactions phase-space suppressed. We describe an algorithm to represent, regularize, evaluate class of thermal for general momenta, masses, chemical potentials, helicity projections. key ingredient is automated inclusion virtual co...