The impact of electroweak Sudakov logarithms on the endpoint of the photon spectrum for wino dark matter annihilation was studied intensively over the last several years. In this work, we extend these results to Higgsino dark matter $\chi_1^0$. We achieve NLL' resummation accuracy for narrow and intermediate spectral energy resolutions, of order $m_W^2 / m_\chi$ and $m_W$, respectively. This is the most accurate prediction to date for the yield of high-energy $\gamma$-rays from $\chi_1^0 \chi_1^0 \to \gamma + X$ annihilation for the energy resolutions realized by current and next-generation telescopes. We also discuss for the first time the effect of power corrections in $m_W/m_\chi$ in this context and argue why they are not sizeable.

We calculate the SU(2)$\times$U(1) electroweak static potential between a fermionic triplet in the broken phase of the Standard Model in the one-loop order (NLO). The one-loop correction provides the leading non-relativistic correction to the large Sommerfeld effect in the annihilation of wino or wino-like dark matter particles $\chi^0$. We find sizeable modifications of the $\chi^0\chi^0$ annihilation cross section and determine the shifts of the resonance locations due to the loop correction to the wino potential.

The annihilation cross section of TeV scale dark matter particles $\chi^0$ with electroweak charges into photons is affected by large quantum corrections due to Sudakov logarithms and the Sommerfeld effect. We calculate the semi-inclusive photon energy spectrum in $\chi^0\chi^0\to \gamma+X$ in the vicinity of the maximal photon energy $E_\gamma = m_\chi$ with NLL' accuracy in an all-order summation of the electroweak perturbative expansion adopting the pure wino model. This results in the most precise theoretical prediction of the annihilation rate for $\gamma$-ray telescopes with photon energy resolution of parametric order $m_W^2/m_\chi$ for photons with TeV energies.

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.

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.

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.

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 estimate the B-polarisation induced in the Cosmic Microwave Background by the non-linear evolution of density perturbations. Using the second-order Boltzmann code SONG, our analysis incorporates, for the first time, all physical effects at recombination. We also include novel contributions from the redshift part of the Boltzmann equation and from the bolometric definition of the temperature in the presence of polarisation. The remaining line-of-sight terms (lensing and time-delay) have previously been studied and must be calculated non-perturbatively. The intrinsic B-mode polarisation is present independent of the initial conditions and might contaminate the signal from primordial gravitational waves. We find this contamination to be comparable to a primordial tensor-to-scalar ratio of $r\simeq10^{-7}$ at the angular scale $\ell\simeq100\,$, where the primordial signal peaks, and $r\simeq 5 \cdot 10^{-5}$ at $\ell\simeq700\,$, where the intrinsic signal peaks. Therefore, we conclude that the intrinsic B-polarisation from second-order effects is not likely to contaminate future searches of primordial gravitational waves.

B-mode polarization of the cosmic background radiation is induced from purely scalar primordial sources at second order in perturbations of the homogeneous, isotropic universe. We calculate the B-mode angular power spectrum C_l^{BB} sourced by the second-order scattering term in the full second-order Boltzmann equations for the polarized radiation phase-space density, which have recently become available. We find that at l\approx 200 the second-order effect is comparable to the first-order effect for a tensor-to-scalar ratio of r=10^{-6}, and to about 2\cdot 10^{-4} at l\approx 1000. It is always negligible relative to the weak-lensing induced contribution.

Within the Closed Time Path (CTP) framework, we derive kinetic equations for particle distribution functions that describe leptogenesis in the presence of several lepton flavours. These flavours have different Standard-Model Yukawa couplings, which induce flavour-sensitive scattering processes and thermal dispersion relations. Kinetic equilibrium, which is rapidly established and maintained via gauge interactions, allows to simplify these equations to kinetic equations for the matrix of lepton charge densities. In performing this simplification, we notice that the rapid flavour-blind gauge interactions damp the flavour oscillations of the leptons. Leptogenesis turns out to be in the parametric regime where the flavour oscillations are overdamped and flavour decoherence is mainly induced by flavour sensitive scatterings. We solve the kinetic equations for the lepton number densities numerically and show that they interpolate between the unflavoured and the fully flavoured regimes within the intermediate parametric region, where neither of these limits is applicable.

Non-gaussianity and B-mode polarization are particularly interesting features of the cosmic microwave background, as -- at least in the standard model of cosmology -- their only sources to first order in cosmological perturbation theory are primordial, possibly generated during inflation. If the primordial sources are small, the question arises how large is the non-gaussianity and B-mode background induced in second-order from the initially gaussian and scalar perturbations. In this paper we derive the Boltzmann hierarchy for the microwave background photon phase-space distributions at second order in cosmological perturbation theory including the complete polarization information, providing the basis for further numerical studies. As an aside we note that the second-order collision term contains new sources of B-mode polarization and that no polarization persists in the tight-coupling limit.