Search

We perform a global fit of dark matter interactions with nucleons using a non-relativistic effective operator description, considering both direct detection and neutrino data. We examine the impact of combining the direct detection experiments CDMSlite, CRESST-II, CRESST-III, DarkSide-50, LUX, LZ, PandaX-II, PandaX-4T, PICO-60, SIMPLE, SuperCDMS, XENON100, and XENON1T along with neutrino data from IceCube and ANTARES. While current neutrino telescope data lead to increased sensitivity compared to underground nuclear scattering experiments for dark matter masses above 100 GeV, our future projections show that the next generation of underground experiments will significantly outpace solar searches for most dark matter-nucleon elastic scattering interactions., Comment: 12+9 pages, 26 figures, Likelihoods available at https://zenodo.org/records/10032218

Using the GAMBIT global fitting framework, we constrain the MSSM with an eV-scale gravitino as the lightest supersymmetric particle, and the six electroweakinos (neutralinos and charginos) as the only other light new states. We combine 15 ATLAS and 12 CMS searches at 13\,TeV, along with a large collection of ATLAS and CMS measurements of Standard Model signatures. This model, which we refer to as the $\tilde G$-EWMSSM, exhibits quite varied collider phenomenology due to its many permitted electroweakino production processes and decay modes. Characteristic $\tilde G$-EWMSSM signal events have two or more Standard Model bosons and missing energy due to the escaping gravitinos. While much of the $\tilde G$-EWMSSM parameter space is excluded, we find several viable parameter regions that predict phenomenologically rich scenarios with multiple neutralinos and charginos within the kinematic reach of the LHC during Run 3, or the High Luminosity LHC. In particular, we identify scenarios with Higgsino-dominated electroweakinos as light as 140 GeV that are consistent with our combined set of collider searches and measurements. The full set of $\tilde G$-EWMSSM parameter samples and GAMBIT input files generated for this work is available via Zenodo., Comment: 38 pages, 12 figures, 5 tables

Global fits explore different parameter regions of a given model and apply constraints obtained at many energy scales. This makes it challenging to perform global fits of simplified models, which may not be valid at high energies. In this study, we derive a unitarity bound for a simplified vector dark matter model with an $s$-channel vector mediator, and apply it to global fits of this model with \GB in order to correctly interpret missing energy searches at the LHC. Two parameter space regions emerge as consistent with all experimental constraints, corresponding to different annihilation modes of the dark matter. We show that although these models are subject to strong validity constraints, they are currently most strongly constrained by measurements less sensitive to the high-energy behaviour of the theory. Understanding when these models cannot be consistently studied will become increasingly relevant as they are applied to LHC Run 3 data., Comment: 15 pages, 7 figures, 3 tables, v2: matches EPJC version

Quantum computing has the potential to offer substantial computational advantages over conventional computing. Recent advances in quantum computing hardware and algorithms have enabled a class of classically parallel quantum workloads, whereby individual quantum circuits can execute independently on many quantum processing units. Here, we present the full-stack software framework developed at Quantum Brilliance to enable multi-modal parallelism for hybrid quantum workloads. Our software provides the capability to distribute quantum workloads across multiple quantum accelerators hosted by nodes of a locally-networked cluster, via the industry-standard MPI (Message Passing Interface) protocol, or to distribute workloads across a large number of cloud-hosted quantum accelerators.

Simplified models provide a useful way to study the impacts of a small number of new particles on experimental observables and the interplay of those observables, without the need to construct an underlying theory. In this study, we perform global fits of simplified dark matter models with GAMBIT using an up-to-date set of likelihoods for indirect detection, direct detection and collider searches. We investigate models in which a scalar or fermionic dark matter candidate couples to quarks via an s-channel vector mediator. Large parts of parameter space survive for each model. In the case of Dirac or Majorana fermion dark matter, excesses in LHC monojet searches and relic density limits tend to prefer the resonance region, where the dark matter has approximately half the mass of the mediator. A combination of vector and axial-vector couplings to the Dirac candidate also leads to competing constraints from direct detection and unitarity violation., Comment: 21 pages, 11 figures, 4 tables, v2: matches EPJC version

This is the Snowmass2021 Energy Frontier (EF) Beyond the Standard Model (BSM) report. It combines the EF topical group reports of EF08 (Model-specific explorations), EF09 (More general explorations), and EF10 (Dark Matter at Colliders). The report includes a general introduction to BSM motivations and the comparative prospects for proposed future experiments for a broad range of potential BSM models and signatures, including compositeness, SUSY, leptoquarks, more general new bosons and fermions, long-lived particles, dark matter, charged-lepton flavor violation, and anomaly detection., Comment: 108 pages + 38 pages references and appendix, 37 figures, Report of the Topical Group on Beyond the Standard Model Physics at Energy Frontier for Snowmass 2021. The first nine authors are the Conveners, with Contributions from the other authors

Axion-like particles (ALPs) decaying into photons are known to affect a wide range of astrophysical and cosmological observables. In this study we focus on ALPs with masses in the keV-MeV range and lifetimes between $10^4$ and $10^{13}$ seconds, corresponding to decays between the end of Big Bang Nucleosynthesis and the formation of the Cosmic Microwave Background (CMB). Using the CosmoBit module of the global fitting framework GAMBIT, we combine state-of-the-art calculations of the irreducible ALP freeze-in abundance, primordial element abundances (including photodisintegration through ALP decays), CMB spectral distortions and anisotropies, and constraints from supernovae and stellar cooling. This approach makes it possible for the first time to perform a global analysis of the ALP parameter space while varying the parameters of $\Lambda$CDM as well as several nuisance parameters. We find a lower bound on the ALP mass of around $m_a > 300\,\text{keV}$, which can only be evaded if ALPs are stable on cosmological timescales. Future observations of CMB spectral distortions with a PIXIE-like mission are expected to improve this bound by two orders of magnitude., Comment: 29+16 pages, 9 figures. V2 corresponds to the published version. Auxiliary material available on Zenodo at https://zenodo.org/record/6573347

In this work, we summarize a global fit study of Type-II two Higgs doublet models (2HDM), and explore the impact of future SM-like Higgs and Z-pole precision measurements on the allowed parameter space. The work is based on the study results of a global fit of 2HDMs with the tool GAMBIT, utilising various current constraints including theoretical constraints (unitarity, perturbativity and vacuum stability), Higgs searches at colliders, electroweak physics and flavour constraints. We further investigate the ability of future facilities, such as the HL-LHC, CEPC, ILC and FCC-ee to explore the 2HDM parameter space., Comment: 9 pages, 3 plots, contribution to Snowmass 2021

Sky survey telescopes and powerful targeted telescopes play complementary roles in astronomy. In order to investigate the nature and characteristics of the motions of very faint objects, a flexibly-pointed instrument capable of high astrometric accuracy is an ideal complement to current astrometric surveys and a unique tool for precision astrophysics. Such a space-based mission will push the frontier of precision astrometry from evidence of Earth-mass habitable worlds around the nearest stars, to distant Milky Way objects, and out to the Local Group of galaxies. As we enter the era of the James Webb Space Telescope and the new ground-based, adaptive-optics-enabled giant telescopes, by obtaining these high precision measurements on key objects that Gaia could not reach, a mission that focuses on high precision astrometry science can consolidate our theoretical understanding of the local Universe, enable extrapolation of physical processes to remote redshifts, and derive a much more consistent picture of cosmological evolution and the likely fate of our cosmos. Already several missions have been proposed to address the science case of faint objects in motion using high precision astrometry missions: NEAT proposed for the ESA M3 opportunity, micro-NEAT for the S1 opportunity, and Theia for the M4 and M5 opportunities. Additional new mission configurations adapted with technological innovations could be envisioned to pursue accurate measurements of these extremely small motions. The goal of this White Paper is to address the fundamental science questions that are at stake when we focus on the motions of faint sky objects and to briefly review instrumentation and mission profiles., Comment: arXiv admin note: substantial text overlap with arXiv:1910.08028, arXiv:1707.01348

Asymmetric dark matter (ADM) that is captured in stars can act as an efficient conductor of heat. Small ADM-induced changes in a star's temperature gradient are known to alter neutrino fluxes and asteroseismological signatures, erase convective cores and modify a star's main sequence lifetime. The Sun's proximity to us makes it an ideal laboratory for studying these effects. However, the two formalisms commonly used to parametrize such heat transport were developed over 30 years ago, and calibrated with a single set of simulations. What's more, both are based on assumptions that break down at the Knudsen transition, where heat transport is maximized. We construct a Monte Carlo simulation to exactly solve the Boltzmann collision equation, determining the steady-state distribution and luminosity carried in stars by ADM with cross sections that depend on velocity and momentum. We find that, although the established (Gould and Raffelt) formalism based on local thermal equilibrium does well for constant cross sections, the isothermal (Spergel and Press) method actually performs better across all models with a simple, universal rescaling function. Based on simulation results, we provide recommendations on the parametrization of DM heat transport in stellar evolution models., Comment: 38 pages, 14 figures. Updated to match published version

We introduce the GAMBIT Universal Model Machine (GUM), a tool for automatically generating code for the global fitting software framework GAMBIT, based on Lagrangian-level inputs. GUM accepts models written symbolically in FeynRules and SARAH formats, and can use either tool along with MadGraph and CalcHEP to generate GAMBIT model, collider, dark matter, decay and spectrum code, as well as GAMBIT interfaces to corresponding versions of SPheno, micrOMEGAs, Pythia and Vevacious (C++). In this paper we describe the features, methods, usage, pathways, assumptions and current limitations of GUM. We also give a fully worked example, consisting of the addition of a Majorana fermion simplified dark matter model with a scalar mediator to GAMBIT via GUM, and carry out a corresponding fit., Comment: 32 pages, 6 figures, 3 tables

We assess the status of a wide class of WIMP dark matter (DM) models in light of the latest experimental results using the global fitting framework $\textsf{GAMBIT}$. We perform a global analysis of effective field theory (EFT) operators describing the interactions between a gauge-singlet Dirac fermion and the Standard Model quarks, the gluons and the photon. In this bottom-up approach, we simultaneously vary the coefficients of 14 such operators up to dimension 7, along with the DM mass, the scale of new physics and several nuisance parameters. Our likelihood functions include the latest data from $\mathit{Planck}$, direct and indirect detection experiments, and the LHC. For DM masses below 100 GeV, we find that it is impossible to satisfy all constraints simultaneously while maintaining EFT validity at LHC energies. For new physics scales around 1 TeV, our results are influenced by several small excesses in the LHC data and depend on the prescription that we adopt to ensure EFT validity. Furthermore, we find large regions of viable parameter space where the EFT is valid and the relic density can be reproduced, implying that WIMPs can still account for the DM of the universe while being consistent with the latest data., Comment: 37 pages, 11 figures, 5 tables; v2: matches EPJC version

Capt'n General is a FORTRAN90 standalone package that can be used to compute the capture and heat transport of dark matter in stars. It can compute capture rates for constant, velocity and momentum-dependent DM-nucleon elastic scattering cross sections, as well as non-relativistic effective operator interactions. Capt'n General can be interfaced with the GAMBIT global fitting codebase as well as stellar evolution simulation codes such as MESA., Comment: 9 pages, 2 figures. Loosely based on a talk presented at TOOLS2020. Code available at: https://github.com/aaronvincent/captngen

Optimisation problems are ubiquitous in particle and astrophysics, and involve locating the optimum of a complicated function of many parameters that may be computationally expensive to evaluate. We describe a number of global optimisation algorithms that are not yet widely used in particle astrophysics, benchmark them against random sampling and existing techniques, and perform a detailed comparison of their performance on a range of test functions. These include four analytic test functions of varying dimensionality, and a realistic example derived from a recent global fit of weak-scale supersymmetry. Although the best algorithm to use depends on the function being investigated, we are able to present general conclusions about the relative merits of random sampling, Differential Evolution, Particle Swarm Optimisation, the Covariance Matrix Adaptation Evolution Strategy, Bayesian Optimisation, Grey Wolf Optimisation, and the PyGMO Artificial Bee Colony, Gaussian Particle Filter and Adaptive Memory Programming for Global Optimisation algorithms., Comment: Experimental framework publicly available at http://www.github.com/darkmachines/high-dimensional-sampling

Physical theories that depend on many parameters or are tested against data from many different experiments pose unique challenges to statistical inference. Many models in particle physics, astrophysics and cosmology fall into one or both of these categories. These issues are often sidestepped with statistically unsound ad hoc methods, involving intersection of parameter intervals estimated by multiple experiments, and random or grid sampling of model parameters. Whilst these methods are easy to apply, they exhibit pathologies even in low-dimensional parameter spaces, and quickly become problematic to use and interpret in higher dimensions. In this article we give clear guidance for going beyond these procedures, suggesting where possible simple methods for performing statistically sound inference, and recommendations of readily-available software tools and standards that can assist in doing so. Our aim is to provide any physicists lacking comprehensive statistical training with recommendations for reaching correct scientific conclusions, with only a modest increase in analysis burden. Our examples can be reproduced with the code publicly available at https://doi.org/10.5281/zenodo.4322283., Comment: 15 pages, 4 figures. extended discussions. closely matches version accepted for publication

We determine the upper limit on the mass of the lightest neutrino from the most robust recent cosmological and terrestrial data. Marginalizing over possible effective relativistic degrees of freedom at early times ($N_\mathrm{eff}$) and assuming normal mass ordering, the mass of the lightest neutrino is less than 0.037 eV at 95% confidence; with inverted ordering, the bound is 0.042 eV. These results improve upon the strength and robustness of other recent limits and constrain the mass of the lightest neutrino to be barely larger than the largest mass splitting. We show the impacts of realistic mass models, and different sources of $N_\mathrm{eff}$., Comment: 5 pages, 2 figures + Appendix. Full dataset available at https://doi.org/10.5281/zenodo.4005381 (v3: Matches version published in PRD)

We introduce $\sf{CosmoBit}$, a module within the open-source $\sf{GAMBIT}$ software framework for exploring connections between cosmology and particle physics with joint global fits. $\sf{CosmoBit}$ provides a flexible framework for studying various scenarios beyond $\Lambda$CDM, such as models of inflation, modifications of the effective number of relativistic degrees of freedom, exotic energy injection from annihilating or decaying dark matter, and variations of the properties of elementary particles such as neutrino masses and the lifetime of the neutron. Many observables and likelihoods in $\sf{CosmoBit}$ are computed via interfaces to $\sf{AlterBBN}$, $\sf{CLASS}$, $\sf{DarkAges}$, $\sf{MontePython}$, $\sf{MultiModeCode}$, and $\sf{plc}$. This makes it possible to apply a wide range of constraints from large-scale structure, Type Ia supernovae, Big Bang Nucleosynthesis and the cosmic microwave background. Parameter scans can be performed using the many different statistical sampling algorithms available within the $\sf{GAMBIT}$ framework, and results can be combined with calculations from other $\sf{GAMBIT}$ modules focused on particle physics and dark matter. We include extensive validation plots and a first application to scenarios with non-standard relativistic degrees of freedom and neutrino temperature, showing that the corresponding constraint on the sum of neutrino masses is much weaker than in the standard scenario., Comment: 44 pages (including executive summary), 11 figures + appendices. Code available at https://github.com/GambitBSM/gambit_1.5 and all datasets available at https://doi.org/10.5281/zenodo.3975642 . V3: minor revisions, matches published version

Air showers, produced by the interaction of energetic cosmic rays with the atmosphere, are an excellent alternative to study particle physics at energies beyond any human-made particle accelerator. For that, it is necessary to identify first the mass composition of the primary cosmic ray (and its energy). None of the existing high energy interaction models have been able to reproduce coherently all air shower observables over the entire energy and zenith angle phase space. This is despite having tried all possible combinations for the cosmic ray mass composition. This proposal outlines a self-consistent strategy to study high energy particle interactions and identify the energy spectra and mass composition of cosmic rays. This strategy involves the participation of different particle accelerators and astrophysics experiments. This is important to cover the entire cosmic ray energy range and a larger phase-space of shower observables to probe the high energy interaction models., Comment: Snowmass2021 - Letter of Interest

The excess of electron recoil events seen by the XENON1T experiment has been interpreted as a potential signal of axion-like particles (ALPs), either produced in the Sun, or constituting part of the dark matter halo of the Milky Way. It has also been explained as a consequence of trace amounts of tritium in the experiment. We consider the evidence for the solar and dark-matter ALP hypotheses from the combination of XENON1T data and multiple astrophysical probes, including horizontal branch stars, red giants, and white dwarfs. We briefly address the influence of ALP decays and supernova cooling. While the different datasets are in clear tension for the case of solar ALPs, all measurements can be simultaneously accommodated for the case of a sub-dominant fraction of dark-matter ALPs. Nevertheless, this solution requires the tuning of several a priori unknown parameters, such that for our choices of priors a Bayesian analysis shows no strong preference for the ALP interpretation of the XENON1T excess over the background hypothesis., Comment: 40 pages, 9 figures. V3 corresponds to the accepted manuscript; added appendices on the DIC and MC simulations; supplementary material now available on Zenodo at https://doi.org/10.5281/zenodo.4384061

The evaluation of higher-order cross-sections is an important component in the search for new physics, both at hadron colliders and elsewhere. For most new physics processes of interest, total cross-sections are known at next-to-leading order (NLO) in the strong coupling $\alpha_s$, and often beyond, via either higher-order terms at fixed powers of $\alpha_s$, or multi-emission resummation. However, the computation time for such higher-order cross-sections is prohibitively expensive, and precludes efficient evaluation in parameter-space scans beyond two dimensions. Here we describe the software tool $\textsf{xsec}$, which allows for fast evaluation of cross-sections based on the use of machine-learning regression, using distributed Gaussian processes trained on a pre-generated sample of parameter points. This first version of the code provides all NLO Minimal Supersymmetric Standard Model strong-production cross-sections at the LHC, for individual flavour final states, evaluated in a fraction of a second. Moreover, it calculates regression errors, as well as estimates of errors from higher-order contributions, from uncertainties in the parton distribution functions, and from the value of $\alpha_s$. While we focus on a specific phenomenological model of supersymmetry, the method readily generalises to any process where it is possible to generate a sufficient training sample., Comment: Accepted version

We report on the status of efforts to improve the reinterpretation of searches and measurements at the LHC in terms of models for new physics, in the context of the LHC Reinterpretation Forum. We detail current experimental offerings in direct searches for new particles, measurements, technical implementations and Open Data, and provide a set of recommendations for further improving the presentation of LHC results in order to better enable reinterpretation in the future. We also provide a brief description of existing software reinterpretation frameworks and recent global analyses of new physics that make use of the current data., Comment: 58 pages, minor revision following comments from SciPost referees

The Global and Modular Beyond-Standard Model Inference Tool (GAMBIT) is an open source software framework for performing global statistical fits of particle physics models, using a wide range of particle and astroparticle data. In this review, we describe the design principles of the package, the statistical and sampling frameworks, the experimental data included, and the first two years of physics results generated with it. This includes supersymmetric models, axion theories, Higgs portal dark matter scenarios and an extension of the Standard Model to include right-handed neutrinos. Owing to the broad spectrum of physics scenarios tackled by the GAMBIT community, this also serves as a convenient, self-contained review of the current experimental and theoretical status of the most popular models of dark matter., Comment: 108 pages, 28 figures, 3 tables; matches version accepted by Progress in Particle and Nuclear Physics; v3 includes a small correction to the acknowledgements

We present a two-field model that realises inflation and the observed density of dark energy today, whilst solving the fine-tuning problems inherent in quintessence models. One field acts as the inflaton, generically driving the other to a saddle-point of the potential, from which it acts as a quintessence field following electroweak symmetry breaking. The model exhibits essentially no sensitivity to the initial value of the quintessence field, naturally suppresses its interactions with other fields, and automatically endows it with a small effective mass in the late Universe. The magnitude of dark energy today is fixed by the height of the saddle point in the potential, which is dictated entirely by the scale of electroweak symmetry breaking., Comment: v2: 13 pages, minor improvements, version accepted in JHEP

We present global fits of cosmologically stable axion-like particle and QCD axion models in the mass range 0.1 neV to 10 eV. We focus on the case where the Peccei-Quinn symmetry is broken before the end of inflation, such that the initial value of the axion field can be considered to be homogeneous throughout the visible Universe. We include detailed likelihood functions from light-shining-through-wall experiments, haloscopes, helioscopes, the axion relic density, horizontal branch stars, supernova 1987A, white dwarf cooling, and gamma-ray observations. We carry out both frequentist and Bayesian analyses, with and without the inclusion of white dwarf cooling. We explore the degree of fine-tuning present in different models and identify parameter regions where it is possible for QCD axion models to account for both the dark matter in the Universe and the cooling hints, comparing them to specific DFSZ- and KSVZ-type models. We find the most credible parameter regions, allowing us to set (prior-dependent) upper and lower bounds on the axion mass. Our analysis also suggests that QCD axions in this scenario most probably make up a non-negligible but sub-dominant component of the dark matter in the Universe., Comment: 77 pages, 29 figures. Version v3: revised version of the paper due to a bug in the code (see Footnote 2); please use GAMBIT 1.3.1 or newer. Updated supplementary material available on Zenodo (version v3) at http://dx.doi.org/10.5281/zenodo.2615235

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., Comment: 38 pages, 16 figures, v3 is the version accepted by EPJC

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., Comment: 30 pages, 12 figures. v2 matches the published version

The search for the particle nature of dark matter has given rise to a number of experimental, theoretical and statistical challenges. Here, we report on a number of these statistical challenges and new techniques to address them, as discussed in the DMStat workshop held Feb 26 - Mar 3 2018 at the Banff International Research Station for Mathematical Innovation and Discovery (BIRS) in Banff, Alberta., Comment: 32 pages, 8 figures, 331 references

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$., Comment: 25 pages, 13 figures. v2 matches the published version. Supplementary data available at https://doi.org/10.5281/zenodo.1298565

I provide a brief summary of the status of supersymmetric models with parameters defined at either the unification or weak scale, based on global fits using the GAMBIT framework., Comment: 4 pages, 3 figures, proceedings of invited contribution to Rencontres de Moriond QCD 2018

The radiatively-induced splitting of masses in electroweak multiplets is relevant for both collider phenomenology and dark matter. Precision two-loop 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 two-loop self-energy calculations for the wino model by obtaining consistent input parameters to the calculation via two-loop renormalisation-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., Comment: 31 pages, 10 figures, 2 Tables

I give a brief review of some of the implications of Fermi data for theories of the identity of dark matter, and their combination with data from other complementary probes. I also preview some of the prospects for probing such models with future data., Comment: 8 pages, 5 figures. To appear in the proceedings of the 7th International Fermi Symposium, Garmish, Oct 15-20 2017

I give a brief review of a few recent developments and future directions in the search for dark matter using high-energy 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., Comment: 7 pages, 4 figures, to appear in the Proceedings of Rencontres du Vietnam: Exploring the Dark Universe, Quy Nhon, Vietnam, July 23-28

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 radiatively-corrected 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., Comment: 21 pages, 6 figures, submitted to EPJC

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 Earth-like 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 micro-arcsecond level. This is about 1000 times better than ESA/Gaia's accuracy for the brightest objects and represents a factor 10-30 improvement for the faintest stars (depending on the exact observational program). In the version submitted to ESA, we proposed an optical (350-1000nm) on-axis 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., Comment: http://theia.phyip3.dur.ac.uk

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 and T-Walk; in ten or more dimensions, Diver substantially outperforms the other three samplers on all metrics., Comment: 51 pages, 25 figures, 2 tables, matches version accepted for publication in EPJC

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 global-fitting 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., Comment: 25 pages, 2 figures, 9 tables. v2: References added

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 3-BIT-HIT, a command-line 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., Comment: 67 pages, 9 figures, updated to match version published in EPJC

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., Comment: 19 pages, 5 figures, matches published version. v3 updates the Zenodo URL for supplementary data (available at https://doi.org/10.5281/zenodo.801510)

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., Comment: 31 pages, 14 figures; v2 matches version accepted by EPJC; v3 has corrected link to supplementary data on Zenodo

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 large-scale 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 state-of-the-art 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., Comment: 31 pages, 7 figures, 17 tables. Replacement version accepted by EPJC

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., Comment: 50 pages, 21 figures, 7 tables, v2 accepted for publication in EPJC, v3 update Zenodo link

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., Comment: 84 pages, 8 Figures. Code available at http://gambit.hepforge.org. v2: Main paper now matches version accepted by EPJC. This version also adds an addendum submitted to EPJC on Mathematica support in GAMBIT v1.1

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., Comment: 53 pages, 9 figures, 25 tables. v2: Updated to reflect new functionality in version 1.1.0, including the addition of Xenon1T and PICO-60 likelihoods. Minor corrections and clarifications throughout. Matches published version

Dark matter with momentum- or velocity-dependent interactions with nuclei has shown significant promise for explaining the so-called Solar Abundance Problem, a longstanding discrepancy between solar spectroscopy and helioseismology. The best-fit models are all rather light, typically with masses in the range of 3-5 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 momentum-dependent 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 better-fitting solutions to the Solar Abundance Problem, but also improving a number of others. As a by-product, 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 well-known saturation limit., Comment: 32 pages including appendices; 10 figures. References added and minor corrections made to match published version

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 non-trivial 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 point-like 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., Comment: 14 pages, 4 figures