If indirect detection searches are to be used to discriminate between dark matter particle models, it is crucial to understand the expected energy spectra of secondary particles such as neutrinos, charged antiparticles and gamma-rays emerging from dark matter annihilations in the local Universe. In this work we study the effect that both the choice of event generator and the polarisation of the final state particles can have on these predictions. For a variety of annihilation channels and dark matter masses, we compare yields obtained with Pythia8 and Herwig7 of all of the aforementioned secondary particle species. We investigate how polarised final states can change these results and do an extensive study of how the polarisation can impact the expected flux of neutrinos from dark matter annihilations in the centre of the Sun. We find that differences between the event generators are larger for yields of hadronic end products such as antiprotons, than for leptonic end products. Concerning polarisation, we conversely find the largest differences in the leptonic spectra. The large differences in the leptonic spectra point to the importance of including polarisation effects in searches for neutrinos from dark matter annihilations in the Sun. However, we find that these differences are ultimately somewhat washed out by propagation effects of the neutrinos in the Sun.

It was recently suggested that dark matter consists of ~GeV particles that carry baryon number and mix with the neutron. We demonstrate that this could allow for resonant dark matter-neutron oscillations in the early universe, at finite temperature, leading to low-scale baryogenesis starting from a primordial dark matter asymmetry. In this scenario, the asymmetry transfer happens around 30 MeV, just before big bang nucleosynthesis. We illustrate the idea using a model with a dark U(1)' gauge interaction, which has recently been suggested as a way of addressing the neutron lifetime anomaly. The asymmetric dark matter component of this model is both strongly self-interacting and leads to a suppression of matter density perturbations at small scales, allowing to mitigate the small-scale problems of cold dark matter cosmology. Future CMB experiments will be able to consistently probe, or firmly exclude, this scenario.

The existence of supermassive black holes at high redshifts ($z\sim7$) is difficult to accommodate in standard astrophysical scenarios. It has been shown that dark matter models with a subdominant self-interacting component are able to produce early seeds for supermassive black holes through the gravothermal catastrophe. Previous studies used a fluid equation approach, requiring some limiting assumptions. Here we reconsider the problem using $N$-body gravitational simulations starting from the formation of the initial dark matter halo. We consider both elastic and dissipative scattering, and elucidate the interplay between the dark matter microphysics and subsequent accretion of the black hole needed to match the properties of observed high redshift supermassive black holes. We find a region of parameter space in which a small component of self-interacting dark matter can produce the observed high redshift supermassive black holes.

Searches for supersymmetric electroweakinos have entered a crucial phase, as the integrated luminosity of the Large Hadron Collider is now high enough to compensate for their weak production cross-sections. Working in a framework where the neutralinos and charginos are the only light sparticles in the Minimal Supersymmetric Standard Model, we use gambit to perform a detailed likelihood analysis of the electroweakino sector. We focus on the impacts of recent ATLAS and CMS searches with 36 fb$^{-1}$ of 13 TeV proton-proton collision data. We also include constraints from LEP and invisible decays of the $Z$ and Higgs bosons. Under the background-only hypothesis, we show that current LHC searches do not robustly exclude any range of neutralino or chargino masses. However, a pattern of excesses in several LHC analyses points towards a possible signal, with neutralino masses of $(m_{\tilde{\chi}_1^0}, m_{\tilde{\chi}_2^0}, m_{\tilde{\chi}_3^0}, m_{\tilde{\chi}_4^0})$ = (8-155, 103-260, 130-473, 219-502) GeV and chargino masses of $(m_{\tilde{\chi}_1^{\pm}}, m_{\tilde{\chi}_2^{\pm}})$ = (104-259, 224-507) GeV at the 95% confidence level. The lightest neutralino is mostly bino, with a possible modest Higgsino or wino component. We find that this excess has a combined local significance of $3.3\sigma$, subject to a number of cautions. If one includes LHC searches for charginos and neutralinos conducted with 8 TeV proton-proton collision data, the local significance is lowered to 2.9$\sigma$. We briefly consider the implications for dark matter, finding that the correct relic density can be obtained through the Higgs-funnel and $Z$-funnel mechanisms, even assuming that all other sparticles are decoupled. All samples, gambit input files and best-fit models from this study are available on Zenodo.

We present global analyses of effective Higgs portal dark matter models in the frequentist and Bayesian statistical frameworks. Complementing earlier studies of the scalar Higgs portal, we use GAMBIT to determine the preferred mass and coupling ranges for models with vector, Majorana and Dirac fermion dark matter. We also assess the relative plausibility of all four models using Bayesian model comparison. Our analysis includes up-to-date likelihood functions for the dark matter relic density, invisible Higgs decays, and direct and indirect searches for weakly-interacting dark matter including the latest XENON1T data. We also account for important uncertainties arising from the local density and velocity distribution of dark matter, nuclear matrix elements relevant to direct detection, and Standard Model masses and couplings. In all Higgs portal models, we find parameter regions that can explain all of dark matter and give a good fit to all data. The case of vector dark matter requires the most tuning and is therefore slightly disfavoured from a Bayesian point of view. In the case of fermionic dark matter, we find a strong preference for including a CP-violating phase that allows suppression of constraints from direct detection experiments, with odds in favour of CP violation of the order of 100:1. Finally, we present DDCalc 2.0.0, a tool for calculating direct detection observables and likelihoods for arbitrary non-relativistic effective operators.

We present the most comprehensive global fits to date of three supersymmetric models motivated by grand unification: the Constrained Minimal Supersymmetric Standard Model (CMSSM), and its Non-Universal Higgs Mass generalisations NUHM1 and NUHM2. We include likelihoods from a number of direct and indirect dark matter searches, a large collection of electroweak precision and flavour observables, direct searches for supersymmetry at LEP and Runs I and II of the LHC, and constraints from Higgs observables. Our analysis improves on existing results not only in terms of the number of included observables, but also in the level of detail with which we treat them, our sampling techniques for scanning the parameter space, and our treatment of nuisance parameters. We show that stau co-annihilation is now ruled out in the CMSSM at more than 95\% confidence. Stop co-annihilation turns out to be one of the most promising mechanisms for achieving an appropriate relic density of dark matter in all three models, whilst avoiding all other constraints. We find high-likelihood regions of parameter space featuring light stops and charginos, making them potentially detectable in the near future at the LHC. We also show that tonne-scale direct detection will play a largely complementary role, probing large parts of the remaining viable parameter space, including essentially all models with multi-TeV neutralinos.

We study the seven-dimensional Minimal Supersymmetric Standard Model (MSSM7) with the new GAMBIT software framework, with all parameters defined at the weak scale. Our analysis significantly extends previous weak-scale, phenomenological MSSM fits, by adding more and newer experimental analyses, improving the accuracy and detail of theoretical predictions, including dominant uncertainties from the Standard Model, the Galactic dark matter halo and the quark content of the nucleon, and employing novel and highly-efficient statistical sampling methods to scan the parameter space. We find regions of the MSSM7 that exhibit co-annihilation of neutralinos with charginos, stops and sbottoms, as well as models that undergo resonant annihilation via both light and heavy Higgs funnels. We find high-likelihood models with light charginos, stops and sbottoms that have the potential to be within the future reach of the LHC. Large parts of our preferred parameter regions will also be accessible to the next generation of direct and indirect dark matter searches, making prospects for discovery in the near future rather good.

Scalar singlet dark matter is one of the simplest and most predictive realisations of the WIMP (weakly-interacting massive particle) idea. Although the model is constrained from all directions by the latest experimental data, it still has viable regions of parameter space. Another compelling aspect of scalar singlets is their ability to stabilise the electroweak vacuum. Indeed, models of scalar dark matter are not low-energy effective theories, but can be valid all the way to the Planck scale. Using the GAMBIT framework, we present the first global fit to include both the low-energy experimental constraints and the theoretical constraints from UV physics, considering models with a scalar singlet charged under either a $\mathbb{Z}_2$ or a $\mathbb{Z}_3$ symmetry. We show that if the model is to satisfy all experimental constraints, completely stabilise the electroweak vacuum up to high scales, and also remain perturbative to those scales, one is driven to a relatively small region of parameter space. This region has a Higgs-portal coupling slightly less than 1, a dark matter mass of 1 to 2 TeV and a spin-independent nuclear scattering cross-section around 10$^{-45}$ cm$^2$.

One of the simplest viable models for dark matter is an additional neutral scalar, stabilised by a $\mathbb{Z}_2$ symmetry. Using the GAMBIT package and combining results from four independent samplers, we present Bayesian and frequentist global fits of this model. We vary the singlet mass and coupling along with 13 nuisance parameters, including nuclear uncertainties relevant for direct detection, the local dark matter density, and selected quark masses and couplings. We include the dark matter relic density measured by Planck, direct searches with LUX, PandaX, SuperCDMS and XENON100, limits on invisible Higgs decays from the Large Hadron Collider, searches for high-energy neutrinos from dark matter annihilation in the Sun with IceCube, and searches for gamma rays from annihilation in dwarf galaxies with the Fermi-LAT. Viable solutions remain at couplings of order unity, for singlet masses between the Higgs mass and about 300 GeV, and at masses above $\sim$1 TeV. Only in the latter case can the scalar singlet constitute all of dark matter. Frequentist analysis shows that the low-mass resonance region, where the singlet is about half the mass of the Higgs, can also account for all of dark matter, and remains viable. However, Bayesian considerations show this region to be rather fine-tuned.

New decay channels for the neutron into dark matter plus other particles have been suggested for explaining a long-standing discrepancy between the neutron lifetime measured from trapped neutrons versus those decaying in flight. Many such scenarios are already ruled out by their effects on neutron stars, and the decay into dark matter plus photon has been experimentally excluded. Here we explore the decay into a dark Dirac fermion $\chi$ and a dark photon $A'$, which can be consistent with all constraints if $\chi$ is a subdominant component of the dark matter. Neutron star constraints are evaded if the dark photon mass to coupling ratio is $m_{A'}/g' \lesssim (45-60)\,$MeV, depending upon the nuclear equation of state. $g'$ and the kinetic mixing between U(1)$'$ and electromagnetism are tightly constrained by direct and indirect dark matter detection, supernova constraints, and cosmological limits.

We introduce DarkBit, an advanced software code for computing dark matter constraints on various extensions to the Standard Model of particle physics, comprising both new native code and interfaces to external packages. This release includes a dedicated signal yield calculator for gamma-ray observations, which significantly extends current tools by implementing a cascade decay Monte Carlo, as well as a dedicated likelihood calculator for current and future experiments (gamlike). This provides a general solution for studying complex particle physics models that predict dark matter annihilation to a multitude of final states. We also supply a direct detection package that models a large range of direct detection experiments (DDcalc), and provides the corresponding likelihoods for arbitrary combinations of spin-independent and spin-dependent scattering processes. Finally, we provide custom relic density routines along with interfaces to DarkSUSY, micrOMEGAs, and the neutrino telescope likelihood package nuLike. DarkBit is written in the framework of the Global And Modular Beyond the Standard Model Inference Tool (GAMBIT), providing seamless integration into a comprehensive statistical fitting framework that allows users to explore new models with both particle and astrophysics constraints, and a consistent treatment of systematic uncertainties. In this paper we describe its main functionality, provide a guide to getting started quickly, and show illustrative examples for results obtained with DarkBit (both as a standalone tool and as a GAMBIT module). This includes a quantitative comparison between two of the main dark matter codes (DarkSUSY and micrOMEGAs), and application of DarkBit's advanced direct and indirect detection routines to a simple effective dark matter model.

We describe the open-source global fitting package GAMBIT: the Global And Modular Beyond-the-Standard-Model Inference Tool. GAMBIT combines extensive calculations of observables and likelihoods in particle and astroparticle physics with a hierarchical model database, advanced tools for automatically building analyses of essentially any model, a flexible and powerful system for interfacing to external codes, a suite of different statistical methods and parameter scanning algorithms, and a host of other utilities designed to make scans faster, safer and more easily-extendible than in the past. Here we give a detailed description of the framework, its design and motivation, and the current models and other specific components presently implemented in GAMBIT. Accompanying papers deal with individual modules and present first GAMBIT results. GAMBIT can be downloaded from gambit.hepforge.org.

In this conference paper, I give an overview of the capabilities of DarkBit, a module of the GAMBIT global fitting code that calculates a range of dark matter observables and corresponding experimental likelihood functions. Included in the code are limits from the dark matter relic density, multiple direct detection experiments, and indirect searches in gamma-rays and neutrinos. I discuss the capabilities of the code, and then present recent results of GAMBIT scans of the parameter space of the minimal supersymmetric standard model, with a focus on sensitivities of future dark matter searches to the current best fit regions.

There are presently several discrepancies in $b \to s \ell^+ \ell^-$ decays of $B$ mesons suggesting new physics coupling to $b$ quarks and leptons. We show that a $Z'$, with couplings to quarks and muons that can explain the $B$-decay anomalies, can also couple to dark matter in a way that is consistent with its relic abundance, direct detection limits, and hints of indirect detection. The latter include possible excess events in antiproton spectra recently observed by the AMS-02 experiment. We present two models, having a heavy (light) $Z'$ with $m_{Z'}\sim 600\,(12)\,$GeV and fermionic dark matter with mass $m_\chi \sim 50\,(2000)\,$GeV, producing excess antiprotons with energies of $\sim 10\, (300)\,$GeV. The first model is also compatible with fits for the galactic center GeV gamma-ray excess.

The wide range of probes of physics beyond the standard model leads to the need for tools that combine experimental results to make the most robust possible statements about the validity of theories and the preferred regions of their parameter space. Here we introduce a new code for such analyses: GAMBIT, the Global and Modular BSM Inference Tool. GAMBIT is a flexible and extensible framework for global fits of essentially any BSM theory. The code currently incorporates direct and indirect searches for dark matter, limits on production of new particles from the LHC and LEP, complete flavor constraints from LHCb, LHC Higgs production and decay measurements, and various electroweak precision observables. Here we present an overview of the code's capabilities, followed by preliminary results from scans of the scalar singlet dark matter model.

Dark matter (DM) annihilations have been widely studied as a possible explanation of excess gamma rays from the galactic center seen by Fermi/LAT. However most such models are in conflict with constraints from dwarf spheroidals. Motivated by this tension, we show that p-wave annihilating dark matter can easily accommodate both sets of observations due to the lower DM velocity dispersion in dwarf galaxies. Explaining the DM relic abundance is then challenging. We outline a scenario in which the usual thermal abundance is obtained through s-wave annihilations of a metastable particle, that eventually decays into the p-wave annihilating DM of the present epoch. The couplings and lifetime of the decaying particle are constrained by big bang nucleosynthesis, the cosmic microwave background and direct detection, but significant regions of parameter space are viable. A sufficiently large p-wave cross section can be found by annihilation into light mediators, that also give rise to Sommerfeld enhancement. A prediction of the scenario is enhanced annihilations in galaxy clusters.

In supersymmetric theories without R-parity, the gravitino can play the role of a decaying Dark Matter candidate without the problem of late NLSP decays affecting Big Bang Nucleosynthesis. In this work, we elaborate on recently discussed limits on R-parity violating couplings from decays to antideuterons and discuss the implications for two classes of flavor symmetries: horizontal symmetries, and Minimal Flavor Violation. In most of the parameter space the antideuteron constraints on R-parity violating couplings are stronger than low-energy baryon-number-violating processes. Even in the absence of flavor symmetries, we find strong new limits on couplings involving third-generation fields, and discuss the implications for LHC phenomenology. For TeV scale superpartners, we find that the allowed MFV parameter space is a corner with gravitino masses smaller than O(10) GeV and small $\tan\beta$.

The recent discovery of a Higgs boson with mass of about 125 GeV, along with its striking similarity to the prediction from the Standard Model, informs and constrains many models of new physics. The Higgs mass exhausts one out of three input parameters of the minimal, five-dimensional version of universal extra dimension models, the other two parameters being the Kaluza-Klein (KK) scale and the cut-off scale of the theory. The presence of KK fermions with large coupling to the Higgs implies a short-lived electro-weak vacuum, unless the cut-off scale is at most a few times higher than the KK mass scale, providing an additional tight constraint to the theory parameter space. Here, we focus on the lightest KK particle as a dark matter candidate, and investigate the regions of parameter space where such particle has a thermal relic density in accord with the cosmological dark matter density. We find the paradoxical result that, for low enough cutoff scales consistent with vacuum stability, larger than previously thought KK mass scales become preferred to explain the dark matter abundance in the universe. We explain this phenomenon by pinpointing the additional particles which, at such low cutoffs, become close enough in mass to the dark matter candidate to coannihilate with it. We make predictions for both collider and direct dark matter searches that might soon close in on all viable theory parameter space.

The size of the smallest dark matter collapsed structures, or protohalos, is set by the temperature at which dark matter particles fall out of kinetic equilibrium. The process of kinetic decoupling involves elastic scattering of dark matter off of Standard Model particles in the early universe, and the relevant cross section is thus closely related to the cross section for dark matter scattering off of nuclei (direct detection) but also, via crossing symmetries, for dark matter pair production at colliders and for pair annihilation. In this study, we employ an effective field theoretic approach to calculate constraints on the kinetic decoupling temperature, and thus on the size of the smallest protohalos, from a variety of direct, indirect and collider probes of particle dark matter.

Dark matter kinetic decoupling involves elastic scattering of dark matter off of leptons and quarks in the early universe, the same process relevant for direct detection and for the capture rate of dark matter in celestial bodies; the resulting size of the smallest dark matter collapsed structures should thus correlate with quantities connected with direct detection rates and with the flux of high-energy neutrinos from dark matter annihilation in the Sun or in the Earth. In this paper we address this general question in the context of two widely studied and paradigmatic weakly-interacting particle dark matter models: the lightest neutralino of the minimal supersymmetric extension of the Standard Model, and the lightest Kaluza-Klein particle of Universal Extra Dimensions (UED). We argue and show that while the scalar neutralino-nucleon cross section correlates poorly with the kinetic decoupling temperature, the spin-dependent cross section exhibits a strong correlation in a wide range of models. In UED models the correlation is present for both cross sections, and is extraordinarily tight for the spin-dependent case. A strong correlation is also found, for both models, for the flux of neutrinos from the Sun, especially for fluxes large enough to be at potentially detectable levels. We provide analytic guidance and formulae that illustrate our findings.