We delineate the allowed parameter and mass range for a wino-like dark matter particle containing some Higgsino admixture in the MSSM by analysing the constraints from diffuse gamma-rays from the dwarf spheroidal galaxies, galactic cosmic rays, direct detection and cosmic microwave background anisotropies. A complete calculation of the Sommerfeld effect for the mixed-neutralino case is performed. We find that the combination of direct and indirect searches poses significant restrictions on the thermally produced wino-Higgsino dark matter with correct relic density. For $\mu>0$ nearly the entire parameter space considered is excluded, while for $\mu<0$ a substantial region is still allowed, provided conservative assumptions on astrophysical uncertainties are adopted.

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.

This paper concludes the presentation of the non-relativistic effective field theory formalism designed to calculate the radiative corrections that enhance the pair-annihilation cross sections of slowly moving neutralinos and charginos within the general minimal supersymmetric standard model (MSSM). While papers I and II focused on the computation of the tree-level annihilation rates that feed into the short-distance part, here we describe in detail the method to obtain the Sommerfeld factors that contain the enhanced long-distance corrections. This includes the computation of the potential interactions in the MSSM, which are provided in compact analytic form, and a novel solution of the multi-state Schr\"odinger equation that is free from the numerical instabilities generated by large mass splittings between the scattering states. Our results allow for a precise computation of the MSSM neutralino dark matter relic abundance and pair-annihilation rates in the present Universe, when Sommerfeld enhancements are important.

We present a detailed discussion of Sommerfeld enhancements in neutralino dark matter relic abundance calculations for several popular benchmark scenarios in the general MSSM. Our analysis is focused on models with heavy wino- and higgsino-like neutralino LSP and models interpolating between these two scenarios. This work is the first phenomenological application of effective field theory methods that we have developed in earlier work and that allow for the consistent study of Sommerfeld enhancements in non-relativistic neutralino and chargino co-annihilation reactions within the general MSSM, away from the pure-wino and pure-higgsino limits.