One of the simplest viable models for dark matter is an additional neutral scalar, stabilised by a 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 a cornellj@physics.mcgill.ca 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 ∼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.

We introduce ScannerBit, the statistics and sampling module of the public, open-source global fitting framework GAMBIT. ScannerBit provides a standardised interface to different sampling algorithms, enabling the use and comparison of multiple computational methods for inferring profile likelihoods, Bayesian posteriors, and other statistical quantities. The current version offers random, grid, raster, nested sampling, differential evolution, Markov Chain Monte Carlo (MCMC) and ensemble Monte Carlo samplers. We also announce the release of a new standalone differential evolution sampler, Diver, and describe its design, usage and interface to ScannerBit. We subject Diver and three other samplers (the nested sampler MultiNest, the MCMC GreAT, and the native ScannerBit implementation of the ensemble Monte Carlo algorithm T-Walk) to a battery of statistical tests. For this we use a realistic physical likelihood function, based on the scalar singlet model of dark matter. We examine the performance of each sampler as a function of its adjustable settings, and the dimensionality of the sampling problem. We evaluate performance on four metrics: optimality of the best fit found, completeness in exploring the best-fit region, number of likelihood evaluations, and total runtime. For Bayesian posterior estimation at high resolution, T-Walk provides the most accurate and timely mapping of the full parameter space. For profile likelihood analysis in less than about ten dimensions, we find that Diver and MultiNest score similarly in terms of best fit and speed, outperforming GreAT a

We describe the open-source global fitting package GAMBIT: the Global And Modular Beyond-the-StandardModel 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 a e-mail: benjamin.farmer@fysik.su.se b e-mail: anders.kvellestad@nordita.org c e-mail: p.scott@imperial.ac.uk d e-mail: c.weniger@uva.nl e Also Institut Universitaire de France, 103 boulevard Saint-Michel, 75005, Paris, France parameter scanning algorithms, and a host of other utilities designed to make scans faster, safer and more easilyextendible 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.

We present the GAMBIT modules SpecBit, DecayBit and PrecisionBit. Together they provide a new framework for linking publicly available spectrum generators, decay codes and other precision observable calculations in a physically and statistically consistent manner. This allows users to automatically run various combinations of existing codes as if they are a single package. The modular design allows software packages fulfilling the same role to be exchanged freely at runtime, with the results presented in a common format that can easily be passed to downstream dark matter, collider and flavour codes. These modules constitute an essential part of the broader GAMBIT framework, a major new software package for performing global fits. In this paper we present the observable calculations, data, and likelihood functions implemented in the three modules, as well as the conventions and assumptions used in interfacing them with external codes. We also present 3-BIT-HIT, a command-line utility for computing mass spectra, couplings, decays and precision observables in the MSSM, which shows a

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 or a 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–2 TeV and a spin-independent nuclear scattering cross-section around .

In this work we present a simple extension of the Standard Model that contains, as the only new physics component, a massive spin-one matter field in the adjoint representation of SU (2) L . In order to be consistent with perturbative unitarity, the vector field must be odd under a Z 2 symmetry. Radiative corrections make the neutral component of the triplet (V 0 ) slightly lighter than the charged ones. We show that V 0 can be the dark matter particle while satisfying all current bounds if it has a mass between 2.8 and 3.8 TeV. We present the current limit on the model parameter space from highly complementary experimental constraints including dark matter relic density measurement, dark matter direct and indirect detection searches, LHC data on Higgs couplings to photons and LHC data on disappearing track searches. We also show that the two-dimensional parameter space can be fully covered by disappearing track searches at a future 100 TeV hadron collider, which will probe, in particular, the whole mass range relevant for dark matter, thus giving an opportunity to discover or exclude the model. * a.belyaev@soton.ac.uk

The radiatively-induced splitting of masses in electroweak multiplets is relevant for both collider phenomenology and dark matter. Precision two-loop corrections of OðMeVÞ to the triplet mass splitting in the wino limit of the minimal supersymmetric standard model can affect particle lifetimes by up to 40%. We improve on previous two-loop self-energy calculations for the wino model by obtaining consistent input parameters to the calculation via two-loop renormalization-group running, and including the effect of finite light quark masses. We also present the first two-loop calculation of the mass splitting in an electroweak fermionic quintuplet, corresponding to the viable form of minimal dark matter (MDM). We place significant constraints on the lifetimes of the charged and doubly-charged fermions in this model. We find that the two-loop mass splittings in the MDM quintuplet are not constant in the large-mass limit, as might naively be expected from the triplet calculation. This is due to the influence of the additional heavy fermions in loop corrections to the gauge boson propagators.

The radiatively-induced mass splitting between components of an electroweak multiplet is typically of order 100 MeV. This is sufficient to endow the charged components with macroscopically-observable lifetimes, and ensure an electrically-neutral dark matter particle. We show that a commonly used iterative procedure to compute radiativelycorrected pole masses can lead to very different mass splittings than a non-iterative calculation at the same loop order. By estimating the uncertainties of the two one-loop results, we show that the iterative procedure is significantly more sensitive to the choice of renormalisation scale and gauge parameter than the non-iterative method. This can cause the lifetime of the charged component to vary by up to 12 orders of magnitude if iteration is employed. We show that individual pole masses exhibit similar scale-dependence regardless of the procedure, but that the leading scale-dependent terms cancel when computing the mass splitting if and only if the non-iterative procedure is employed. We show that this behaviour persists at two-loop order: the precision of the mass splitting improves in the non-iterative approach, but our results suggest that higher-order corrections do not reduce the uncertainty in the iterative calculation enough to resolve the problem at two-loop order. We conclude that the iterative procedure should not be used for computing pole masses in situations where electroweak mass splittings are phenomenologically relevant.