Precision calculations of dark matter relic abundance

The dark matter annihilation channels sometimes involve sharp resonances. In such cases the usual momentum averaged approximations for computing the DM abundance may not be accurate. We develop an easily accessible momentum dependent framework for computing the DM abundance accurately and efficiently near such features. We apply the method to the case of a singlet scalar dark matter $s$ interacting with SM through higgs portal $\lambda_{\rm hs}s^2 h^2$ and compare the results with different momentum averaged methods. The accuracy of the latter depend strongly on the strength of the elastic interactions and corrections are large if WIMP has negligible interactions beyond the main annihilation channel. In the singlet scalar model however, the standard model scatterings induce an efficient kinetic equilibrium that validates the momentum averaged computation to 20 per cent accuracy. We update the current extent of the allowed region in the light singlet scalar dark matter to $m_{\rm S} \in [56,62.5]$ GeV.