We present MadDM v.3.0, a numerical tool to compute particle dark matter observables in generic new physics models. The new version features a comprehensive and automated framework for dark matter searches at the interface of collider physics, astrophysics and cosmology and is deployed as a plugin of the MadGraph5_aMC@NLO platform, inheriting most of its features. With respect to the previous version, MadDM v.3.0 can now provide predictions for indirect dark matter signatures in astrophysical environments, such as the annihilation cross section at present time and the energy spectra of prompt photons, cosmic rays and neutrinos resulting from dark matter annihilation. MadDM indirect detection features support both $2\to2$ and $2 \to n$ dark matter annihilation processes. In addition, the ability to compare theoretical predictions with experimental constraints is extended by including the Fermi-LAT likelihood for gamma-ray constraints from dwarf spheroidal galaxies and by providing an interface with the nested sampling algorithm PyMultinNest to perform high dimensional parameter scans efficiently. We validate the code for a wide set of dark matter models by comparing the results from MadDM v.3.0 to existing tools and results in the literature.

Weakly-coupled TeV-scale particles may mediate the interactions between normal matter and dark matter. If so, the LHC would produce dark matter through these mediators, leading to the familiar "mono-X" search signatures, but the mediators would also produce signals without missing momentum via the same vertices involved in their production. This document from the LHC Dark Matter Working Group suggests how to compare searches for these two types of signals in case of vector and axial-vector mediators, based on a workshop that took place on September 19/20, 2016 and subsequent discussions. These suggestions include how to extend the spin-1 mediated simplified models already in widespread use to include lepton couplings. This document also provides analytic calculations of the relic density in the simplified models and reports an issue that arose when ATLAS and CMS first began to use preliminary numerical calculations of the dark matter relic density in these models.

We find new and universal relations for the properties of dark matter particles consistent with standard relic abundances. Analysis is based on first characterizing the $s$-channel resonant annihilation process in great detail, keeping track of all velocity-dependence, the presence of multiple scales and treating each physical regime above, below, and close to thresholds separately. The resonant regime as well as extension to include non-resonant processes are then reduced to analytic formulas and inequalities that describe the full range of multi-dimensional numerical work. These results eliminate the need to recompute relic abundance model by model, and reduce calculations to verifying certain scale and parameter combinations are consistent. Remarkably simple formulas describe the relation between the total width of an $s$-channel intermediate particle, the masses and the couplings involved. Eliminating the width in terms of the mass produces new consistency relations between dark matter masses and the intermediate masses. The formulas are general enough to test directly whether new particles can be identified as dark matter. Resonance mass and total width are quantities directly observable at accelerators such as the LHC, and will be sufficient to establish whether new discoveries are consistent with the cosmological bounds on dark matter.

We find model-independent upper limits on rates of dark matter annihilation in galactic halos. The Born approximation generally fails, while exotic threshold enhancements akin to "Sommerfeld factors" also turn out to be baseless The most efficient annihilation mechanism involves perturbatively small decay widths that have largely been ignored. Widths that are very small compared to TeV mass scales suffice to effectively saturate unitarity bounds on annihilation rates. Bound state formation in weakly coupled theories produces small effects due to wave function normalizations. Unitarity shows the Sommerfeld factor cannot produce large changes in cross sections, and serves to identify where those approximations break down.