
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 be easily 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 3BITHIT, a
commandline utility for computing mass spectra, couplings, decays and
precision observables in the MSSM, which shows how the three modules can be
easily used independently of GAMBIT.

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 NonUniversal 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 coannihilation is now ruled out in the CMSSM at more than 95\%
confidence. Stop coannihilation 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 highlikelihood
regions of parameter space featuring light stops and charginos, making them
potentially detectable in the near future at the LHC. We also show that
tonnescale direct detection will play a largely complementary role, probing
large parts of the remaining viable parameter space, including essentially all
models with multiTeV neutralinos.

We study the sevendimensional Minimal Supersymmetric Standard Model (MSSM7)
with the new GAMBIT software framework, with all parameters defined at the weak
scale. Our analysis significantly extends previous weakscale, 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 highlyefficient
statistical sampling methods to scan the parameter space. We find regions of
the MSSM7 that exhibit coannihilation of neutralinos with charginos, stops and
sbottoms, as well as models that undergo resonant annihilation via both light
and heavy Higgs funnels. We find highlikelihood 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.

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 highenergy neutrinos from
dark matter annihilation in the Sun with IceCube, and searches for gamma rays
from annihilation in dwarf galaxies with the FermiLAT. 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 lowmass 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 finetuned.

An excess of gamma rays has been identified at the centre of the Milky Way,
and annihilation of dark matter has been posited as a potential source. This
hypothesis faces significant challenges: difficulty characterizing
astrophysical backgrounds, the need for a nontrivial adiabatic contraction of
the inner part of the Milky Way's dark matter halo, and recent observations of
photon clustering, which suggest that the majority of the excess is due to
unresolved point sources. Here we point out that the apparent pointlike nature
of the emission rules out the dark matter interpretation of the excess
entirely. Attempting to model the emission with dark matter point sources
either worsens the problem with the inner slope, requires an unrealistically
large minihalo fraction toward the Galactic Centre, or overproduces the
observed emission at higher latitudes.

The radiativelyinduced mass splitting between components of an electroweak
multiplet is typically of order 100 MeV. This is sufficient to endow the
charged components with macroscopicallyobservable lifetimes, and ensure an
electricallyneutral 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 noniterative calculation at the same
loop order. By estimating the uncertainties of the two oneloop results, we
show that the iterative procedure is significantly more sensitive to the choice
of renormalisation scale and gauge parameter than the noniterative 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 scaledependence regardless of the procedure, but that the
leading scaledependent terms cancel when computing the mass splitting if and
only if the noniterative procedure is employed. We show that this behaviour
persists at twoloop order: the precision of the mass splitting improves in the
noniterative approach, but our results suggest that higherorder corrections
do not reduce the uncertainty in the iterative calculation enough to resolve
the problem at twoloop order. We conclude that the iterative procedure should
not be used for computing pole masses in situations where electroweak mass
splittings are phenomenologically relevant.

The radiativelyinduced splitting of masses in electroweak multiplets is
relevant for both collider phenomenology and dark matter. Precision twoloop
corrections of $\mathcal{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 twoloop selfenergy
calculations for the wino model by obtaining consistent input parameters to the
calculation via twoloop renormalisationgroup running, and including the
effect of finite light quark masses. We also present the first twoloop
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 doublycharged
fermions in this model. We find that the twoloop mass splittings in the MDM
quintuplet are not constant in the largemass 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.

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 gammaray
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 spinindependent and
spindependent 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 opensource global fitting package GAMBIT: the Global And
Modular BeyondtheStandardModel 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 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 introduce ScannerBit, the statistics and sampling module of the public,
opensource 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 TWalk) 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
bestfit region, number of likelihood evaluations, and total runtime. For
Bayesian posterior estimation at high resolution, TWalk 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 and TWalk;
in ten or more dimensions, Diver substantially outperforms the other three
samplers on all metrics.

I give a brief review of a few recent developments and future directions in
the search for dark matter using highenergy neutrinos from the Sun. This
includes the ability to recast neutrino telescope limits on nuclear scattering
of dark matter to arbitrary new theories, and new calculations of the solar
atmospheric background relevant to such searches. I also touch on applications
to global searches for new physics, and prospects for improving searches for
asymmetric dark matter in the Sun.

In the context of the ESA M5 (medium mission) call we proposed a new
satellite mission, Theia, based on relative astrometry and extreme precision to
study the motion of very faint objects in the Universe. Theia is primarily
designed to study the local dark matter properties, the existence of Earthlike
exoplanets in our nearest star systems and the physics of compact objects.
Furthermore, about 15 $\%$ of the mission time was dedicated to an open
observatory for the wider community to propose complementary science cases.
With its unique metrology system and "point and stare" strategy, Theia's
precision would have reached the sub microarcsecond level. This is about 1000
times better than ESA/Gaia's accuracy for the brightest objects and represents
a factor 1030 improvement for the faintest stars (depending on the exact
observational program). In the version submitted to ESA, we proposed an optical
(3501000nm) onaxis TMA telescope. Due to ESA Technology readiness level, the
camera's focal plane would have been made of CCD detectors but we anticipated
an upgrade with CMOS detectors. Photometric measurements would have been
performed during slew time and stabilisation phases needed for reaching the
required astrometric precision.

The existence of substructure in halos of annihilating dark matter would be
expected to substantially boost the rate at which annihilation occurs.
Ultracompact minihalos of dark matter (UCMHs) are one of the more extreme
examples of this. The boosted annihilation can inject significant amounts of
energy into the gas of a galaxy over its lifetime. Here we determine the impact
of the boost factor from UCMH substructure on the heating of galactic gas in a
Milky Waytype galaxy, by means of Nbody simulation. If $1\%$ of the dark
matter exists as UCMHs, the corresponding boost factor can be of order $10^5$.
For reasonable values of the relevant parameters (annihilation cross section
$3\times10^{26} ~\textrm{cm}^3~ \textrm{s}^{1}$, dark matter mass 100 GeV,
10% heating efficiency), we show that the presence of UCMHs at the 0.1% level
would inject enough energy to eject significant amounts of gas from the halo,
potentially preventing star formation within $\sim$1 kpc of the halo centre.

We describe ColliderBit, a new code for the calculation of high energy
collider observables in theories of physics beyond the Standard Model (BSM).
ColliderBit features a generic interface to BSM models, a unique parallelised
Monte Carlo event generation scheme suitable for largescale supercomputer
applications, and a number of LHC analyses, covering a reasonable range of the
BSM signatures currently sought by ATLAS and CMS. ColliderBit also calculates
likelihoods for Higgs sector observables, and LEP searches for BSM particles.
These features are provided by a combination of new code unique to ColliderBit,
and interfaces to existing stateoftheart public codes. ColliderBit is both
an important part of the GAMBIT framework for BSM inference, and a standalone
tool for efficiently applying collider constraints to theories of new physics.

Flavour physics observables are excellent probes of new physics up to very
high energy scales. Here we present FlavBit, the dedicated flavour physics
module of the globalfitting package GAMBIT. FlavBit includes custom
implementations of various likelihood routines for a wide range of flavour
observables, including detailed uncertainties and correlations associated with
LHCb measurements of rare, leptonic and semileptonic decays of B and D mesons,
kaons and pions. It provides a generalised interface to external theory codes
such as SuperIso, allowing users to calculate flavour observables in and beyond
the Standard Model, and then test them in detail against all relevant
experimental data. We describe FlavBit and its constituent physics in some
detail, then give examples from supersymmetry and effective field theory
illustrating how it can be used both as a standalone library for flavour
physics, and within GAMBIT.

Dark matter with momentum or velocitydependent interactions with nuclei has
shown significant promise for explaining the socalled Solar Abundance Problem,
a longstanding discrepancy between solar spectroscopy and helioseismology. The
bestfit models are all rather light, typically with masses in the range of 35
GeV. This is exactly the mass range where dark matter evaporation from the Sun
can be important, but to date no detailed calculation of the evaporation of
such models has been performed. Here we carry out this calculation, for the
first time including arbitrary velocity and momentumdependent interactions,
thermal effects, and a completely general treatment valid from the optically
thin limit all the way through to the optically thick regime. We find that
depending on the dark matter mass, interaction strength and type, the mass
below which evaporation is relevant can vary from 1 to 4 GeV. This has the
effect of weakening some of the betterfitting solutions to the Solar Abundance
Problem, but also improving a number of others. As a byproduct, we also
provide an improved derivation of the capture rate that takes into account
thermal and optical depth effects, allowing the standard result to be smoothly
matched to the wellknown saturation limit.

In recent years, a revised set of solar abundances has led to a discrepancy
in the soundspeed profile between helioseismology and theoretical solar
models. Conventional solutions require additional mechanisms for energy
transport within the Sun. Vincent et al. have recently suggested that dark
matter with a momentum or velocity dependent cross section could provide a
solution. In this work, we consider three models of dark matter with such cross
sections and their effect on the stellar structure. In particular, the three
models incorporate dark matter particles interacting through an electromagnetic
dipole moment: an electric dipole, a magnetic dipole or an anapole. Each model
is implemented in the \texttt{DarkStec} stellar evolution program, which
incorporates the effects of dark matter capture and heat transport within the
solar interior. We show that dark matter with an anapole moment of
$\sim1\mathrm{GeV}^{2}$ or magnetic dipole moment of $\sim10^{3}\mu_p$ can
improve the soundspeed profile, small frequency separations and convective
zone radius with respect to the Standard Solar Model. However, the required
dipole moments are strongly excluded by direct detection experiments.

This White Paper is an input to the ongoing discussion about the extension
and refinement of simplified Dark Matter (DM) models. Based on two concrete
examples, we show how existing simplified DM models (SDMM) can be extended to
provide a more accurate and comprehensive framework to interpret and
characterise collider searches. In the first example we extend the canonical
SDMM with a scalar mediator to include mixing with the Higgs boson. We show
that this approach not only provides a better description of the underlying
kinematic properties that a complete model would possess, but also offers the
option of using this more realistic class of scalar mixing models to compare
and combine consistently searches based on different experimental signatures.
The second example outlines how a new physics signal observed in a visible
channel can be connected to DM by extending a simplified model including
effective couplings. This discovery scenario uses the recently observed excess
in the highmass diphoton searches of ATLAS and CMS for a case study to show
that such a pragmatic approach can aid the experimental search programme to
verify/falsify a potential signal and to study its underlying nature. In the
next part of the White Paper we outline other interesting options for SDMM that
could be studied in more detail in the future. Finally, we discuss important
aspects of supersymmetric models for DM and how these could help to develop of
more complete SDMM.

Interpreting samples from likelihood or posterior probability density
functions is rarely as straightforward as it seems it should be. Producing
publicationquality graphics of these distributions is often similarly painful.
In this short note I describe pippi, a simple, publiclyavailable package for
parsing and postprocessing such samples, as well as generating highquality
PDF graphics of the results. Pippi is easily and extensively configurable and
customisable, both in its options for parsing and postprocessing samples, and
in the visual aspects of the figures it produces. I illustrate some of these
using an existing supersymmetric global fit, performed in the context of a
gammaray search for dark matter. Pippi can be downloaded and followed at
http://github.com/patscott/pippi .

Cosmological inflation generates primordial density perturbations on all
scales, including those far too small to contribute to the cosmic microwave
background. At these scales, isolated ultracompact minihalos of dark matter can
form well before standard structure formation, if the perturbations have
sufficient amplitude. Minihalos affect pulsar timing data and are potentially
bright sources of gamma rays. The resulting constraints significantly extend
the observable window of inflation in the presence of cold dark matter,
coupling two of the key problems in modern cosmology.

We estimate the sensitivity of the upcoming CTA gammaray telescope to DM
annihilation at the Galactic centre, improving on previous analyses in a number
of significant ways. First, we perform a detailed analyses of all backgrounds,
including diffuse astrophysical emission for the first time in a study of this
type. Second, we present a statistical framework for including systematic
errors and estimate the consequent degradation in sensitivity. These errors may
come from e.g. event reconstruction, Monte Carlo determination of the effective
area or uncertainty in atmospheric conditions. Third, we show that performing
the analysis on a set of suitably optimised regions of interest makes it
possible to partially compensate for the degradation in sensitivity caused by
systematics and diffuse emission. To probe dark matter with the canonical
thermal annihilation crosssection, CTA systematics like nonuniform variations
in acceptance over a single field of view must be kept below the 0.3% level,
unless the dark matter density rises more steeply in the centre of the Galaxy
than predicted by a typical NavarroFrenkWhite or Einasto profile. For a
contracted $r^{1.3}$ profile, and systematics at the 1% level, CTA can probe
annihilation to $b\bar{b}$ at the canonical thermal level for dark matter
masses between 100 GeV and 10 TeV.

We combine and extend the analyses of effective scalar, vector, Majorana and
Dirac fermion Higgs portal models of dark matter (DM), in which DM couples to
the Standard Model (SM) Higgs boson via an operator of the form
$\mathcal{O}_{\textrm{DM}}\, H^\dagger H$. For the fermion models, we take an
admixture of scalar $\overline{\psi} \psi$ and pseudoscalar $\overline{\psi}
i\gamma_5 \psi$ interaction terms. For each model, we apply constraints on the
parameter space based on the Planck measured DM relic density and the LHC
limits on the Higgs invisible branching ratio. For the first time, we perform a
consistent study of the indirect detection prospects for these models based on
the WMAP7/Planck observations of the cosmic microwave background, a combined
analysis of 15 dwarf spheroidal galaxies by FermiLAT and the upcoming
Cherenkov Telescope Array (CTA). We also perform a correct treatment of the
momentumdependent direct search cross section that arises from the
pseudoscalar interaction term in the fermionic DM theories. We find, in line
with previous studies, that current and future direct search experiments such
as LUX and XENON1T can exclude much of the parameter space, and we demonstrate
that a joint observation in both indirect and direct searches is possible for
high mass weakly interacting massive particles. In the case of a pure
pseudoscalar interaction of a fermionic DM candidate, future gammaray searches
are the only class of experiment capable of probing the high mass range of the
theory.

We study the effects of energy transport in the Sun by asymmetric dark matter
with momentum and velocitydependent interactions, with an eye to solving the
decadeold Solar Abundance Problem. We study effective theories where the dark
matternucleon scattering crosssection goes as $v_{\rm rel}^{2n}$ and $q^{2n}$
with $n = 1, 0, 1 $ or $2$, where $v_{\rm rel}$ is the dark matternucleon
relative velocity and $q$ is the momentum exchanged in the collision. Such
crosssections can arise generically as leading terms from the most basic
nonstandard DMquark operators. We employ a highprecision solar simulation
code to study the impact on solar neutrino rates, the sound speed profile,
convective zone depth, surface helium abundance and small frequency
separations. We find that the majority of models that improve agreement with
the observed sound speed profile and depth of the convection zone also reduce
neutrino fluxes beyond the level that can be reasonably accommodated by
measurement and theory errors. However, a few specific points in parameter
space yield a significant overall improvement. A 35 GeV DM particle with
$\sigma_{SI} \propto q^2$ is particularly appealing, yielding more than a
$6\sigma$ improvement with respect to standard solar models, while being
allowed by direct detection and collider limits. We provide full analytical
capture expressions for $q$ and $v_{\rm rel}$dependent scattering, as well as
complete likelihood tables for all models.

Smallscale dark matter structure within the Milky Way is expected to affect
pulsar timing. The change in gravitational potential induced by a dark matter
halo passing near the line of sight to a pulsar would produce a varying delay
in the light travel time of photons from the pulsar. Individual transits
produce an effect that would either be too rare or too weak to be detected in
30year pulsar observations. However, a population of dark matter subhalos
would be expected to produce a detectable effect on the measured properties of
pulsars if the subhalos constitute a significant fraction of the total halo
mass. The effect is to increase the dispersion of measured period derivatives
across the pulsar population. By statistical analysis of the ATNF pulsar
catalogue, we place an upper limit on this dispersion of $\log \sigma_{\dot{P}}
\leq 17.05$. We use this to place strong upper limits on the number density of
ultracompact minihalos within the Milky Way. These limits are completely
independent of the particle nature of dark matter.

Ultracompact Minihalos (UCMHs) have been proposed as a type of dark matter
substructure seeded by largeamplitude primordial perturbations and
topological defects. UCMHs are expected to survive to the present era, allowing
constraints to be placed on their cosmic abundance using observations within
our own Galaxy. Constraints on their number density can be linked to conditions
in the early universe that impact structure formation, such as increased
primordial power on small scales, generic weak nonGaussianity, and the
presence of cosmic strings. We use new constraints on the abundance of UCMHs
from pulsar timing to place generalised limits on the parameters of each of
these cosmological scenarios. At some scales, the limits are the strongest to
date, exceeding those from dark matter annihilation. Our new limits have the
added advantage of being independent of the particle nature of dark matter, as
they are based only on gravitational effects.