There is an increasing interest in accurate dark matter relic density predictions, which requires next-to-leading order (NLO) calculations. The method applied up to now uses zero-temperature NLO calculations of annihilation cross sections in the standard Boltzmann equation for freeze-out, and is conceptually problematic, since it ignores the finite-temperature infrared (IR) divergences from soft and collinear radiation and virtual effects. We address this problem systematically by starting from non-equilibrium quantum field theory, and demonstrate on a realistic model that soft and collinear temperature-dependent divergences cancel in the collision term. Our analysis provides justification for the use of the freeze-out equation in its conventional form and determines the leading finite-temperature correction to the annihilation cross section. This turns out to have a remarkably simple structure.

The relic density of TeV-scale wino-like neutralino dark matter in the MSSM is subject to potentially large corrections as a result of the Sommerfeld effect. A recently developed framework enables us to calculate the Sommerfeld-enhanced relic density in general MSSM scenarios, properly treating mixed states and multiple co-annihilating channels as well as including off-diagonal contributions. Using this framework, including on-shell one-loop mass splittings and running couplings and taking into account the latest experimental constraints, we perform a thorough study of the regions of parameter space surrounding the well known pure-wino scenario: namely the effect of sfermion masses being non-decoupled and of allowing non-negligible Higgsino or bino components in the lightest neutralino. We further perform an investigation into the effect of thermal corrections and show that these can safely be neglected. The results reveal a number of phenomenologically interesting but so far unexplored regions where the Sommerfeld effect is sizeable. We find, in particular, that the relic density can agree with experiment for dominantly wino neutralino dark matter with masses ranging from 1.7 to beyond 4 TeV. In light of these results the bounds from Indirect Detection on wino-like dark matter should be revisited.