We study a dark matter production mechanism based on decays of a messenger WIMP-like state into a pair of dark matter particles that are self-interacting via exchange of a light mediator. Its distinctive thermal history allows the mediator to be stable and therefore avoid strong limits from the cosmic microwave background and indirect detection. A natural by-product of this mechanism is a possibility of a late time, i.e., after recombination, transition to subdominant dark radiation component through three-body and one-loop decays to states containing the light mediator. We examine to what extent such a process can help to alleviate the $H_0$ tension. Additionally, the mechanism can provide a natural way of constructing dark matter models with ultra-strong self-interactions that may positively affect the supermassive black hole formation rate. We provide a simple realization of the mechanism in a Higgs portal dark matter model and find a significant region of the parameter space that leads to a mild relaxation of the Hubble tension while simultaneously having the potential of addressing small-scale structure problems of $\Lambda$CDM.

Calculating the abundance of thermally produced dark matter particles has become a standard procedure, with sophisticated methods guaranteeing a precision that matches the percent-level accuracy in the observational determination of the dark matter density. Here, we point out that one of the main assumptions in the commonly adopted formalism, namely local thermal equilibrium during the freeze-out of annihilating dark matter particles, does not have to be satisfied in general. We present two methods for how to deal with such situations, in which the kinetic decoupling of dark matter happens so early that it interferes with the chemical decoupling process: i) an approximate treatment in terms of a coupled system of differential equations for the leading momentum moments of the dark matter distribution, and ii) a full numerical solution of the Boltzmann equation in phase-space. For illustration, we apply these methods to the case of Scalar Singlet dark matter. We explicitly show that even in this simple model the prediction for the dark matter abundance can be affected by up to one order of magnitude compared to the traditional treatment.

We provide an updated and improved study of the prospects of the H.E.S.S. and Cherenkov Telescope Array (CTA) experiments in testing neutralino dark matter in the Minimal Supersymmetric Standard Model with nine free parameters (p9MSSM). We include all relevant experimental constraints and theoretical developments, in particular the calculation of the Sommerfeld enhancement for both present-day annihilation and the relic abundance. We perform a state-of-the-art analysis of the CTA sensitivity with a log-likelihood test ratio statistics and apply it to a numerical scan of the p9MSSM parameter space focusing on a TeV scale dark matter. We find that, assuming Einasto profile of dark matter halo in the Milky Way, H.E.S.S. has already been able to nearly reach the so-called thermal WIMP value, while CTA will go below it by providing a further improvement of at least an order of magnitude. Both H.E.S.S. and CTA are sensitive to several cases for which direct detection cross section will be below the so-called neutrino floor, with H.E.S.S. being sensitive to most of the wino region, while CTA also covering a large fraction of the ~1 TeV higgsino region. We show that CTA sensitivity will be further improved in the monochromatic photon search mode for both single-component and underabundant dark matter.

e-ASTROGAM (enhanced ASTROGAM) is a breakthrough Observatory space mission, with a detector composed by a Silicon tracker, a calorimeter, and an anticoincidence system, dedicated to the study of the non-thermal Universe in the photon energy range from 0.3 MeV to 3 GeV - the lower energy limit can be pushed to energies as low as 150 keV for the tracker, and to 30 keV for calorimetric detection. The mission is based on an advanced space-proven detector technology, with unprecedented sensitivity, angular and energy resolution, combined with polarimetric capability. Thanks to its performance in the MeV-GeV domain, substantially improving its predecessors, e-ASTROGAM will open a new window on the non-thermal Universe, making pioneering observations of the most powerful Galactic and extragalactic sources, elucidating the nature of their relativistic outflows and their effects on the surroundings. With a line sensitivity in the MeV energy range one to two orders of magnitude better than previous generation instruments, e-ASTROGAM will determine the origin of key isotopes fundamental for the understanding of supernova explosion and the chemical evolution of our Galaxy. The mission will provide unique data of significant interest to a broad astronomical community, complementary to powerful observatories such as LIGO-Virgo-GEO600-KAGRA, SKA, ALMA, E-ELT, TMT, LSST, JWST, Athena, CTA, IceCube, KM3NeT, and LISA.

Kinetic decoupling of dark matter typically happens much later than chemical freeze-out. In fact, local thermal equilibrium is an important assumption for the usual relic density calculations based on solving the Boltzmann equation (for its 0-th moment) describing the dark matter number density. But is this assumption always justified? Here we address this question and discuss the consequences of more accurate treatments. The first treatment is relying on the inclusion of higher moments of the Boltzmann equation and the second on solving the evolution of the phase-space distribution function fully numerically. For illustration, these methods are applied to the Scalar Singlet model, often referred to as the simplest WIMP DM possibility from a model-building perspective. It is explicitly shown that even in this simple model the prediction for the dark matter abundance can be affected by as much as one order of magnitude.

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 perform a multichannel analysis of the indirect signals for the Wino Dark Matter, including one-loop electroweak and Sommerfeld enhancement corrections. We derive limits from cosmic ray antiprotons and positrons, from continuum galactic and extragalactic diffuse gamma-ray spectra, from the absence of gamma-ray line features at the galactic center above 500 GeV in energy, from gamma-rays toward nearby dwarf spheroidal galaxies and galaxy clusters, and from CMB power-spectra. Additionally, we show the future prospects for neutrino observations toward the inner Galaxy and from antideuteron searches. For each of these indirect detection probes we include and discuss the relevance of the most important astrophysical uncertainties that can impact the strength of the derived limits. We find that the Wino as a dark matter candidate is excluded in the mass range bellow ~800 GeV from antiprotons and between 1.8 and 3.5 TeV from the absence of a gamma-ray line feature toward the galactic center. Limits from other indirect detection probes confirm the main bulk of the excluded mass ranges.

We compute the present-day Wino dark matter annihilation cross-section including the one-loop radiative corrections together with the fully treated electroweak Sommerfeld effect. We discuss what is the consistent way of incorporating these two corrections simultaneously and why simply using the running coupling constants values at the Wino mass scale is not correct. The results show that up to a few TeV scale the full one-loop computation makes the cross-section smaller up to about 30% with respect to the Sommerfeld enhanced tree level result and are considerably larger than the tree or one-loop level without the Sommerfeld effect.