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 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.

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.

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.

Dedicated underground experiments searching for dark matter have little sensitivity to GeV and sub-GeV masses of dark matter particles. We show that the decay of B mesons to K(K^*) and missing energy in the final state can be an efficient probe of dark matter models in this mass range. We analyze the minimal scalar dark matter model to show that the width of the decay mode with two dark matter scalars B\to KSS may exceed the decay width in the Standard Model channel, B\to K\nu\bar\nu, by up to two orders of magnitude. Existing data from B physics experiments almost entirely exclude dark matter scalars with masses less than 1 GeV. Expected data from B factories will probe the range of dark matter masses up to 2 GeV.

The cosmic microwave background radiation provides unique constraints on cosmological models. In this Letter we present a summary of the spatial properties of the cosmic microwave background radiation based on the full 4 years of COBE DMR observations, as detailed in a set of companion Letters. The anisotropy is consistent with a scale-invariant power law model and Gaussian statistics. With full use of the multi-frequency 4-year DMR data, including our estimate of the effects of Galactic emission, we find a power-law spectral index of $n=1.2\pm 0.3$ and a quadrupole normalization $Q_{rms-PS}=15.3^{+3.8}_{-2.8}$ $\mu$K. For $n=1$ the best-fit normalization is $Q_{rms-PS}\vert_{n=1}=18\pm 1.6$ $\mu$K. These values are consistent with both our previous 1-year and 2-year results. The results include use of the $\ell=2$ quadrupole term; exclusion of this term gives consistent results, but with larger uncertainties. The 4-year sky maps, presented in this Letter, portray an accurate overall visual impression of the anisotropy since the signal-to-noise ratio is ~2 per 10 degree sky map patch. The improved signal-to-noise ratio of the 4-year maps also allows for improvements in Galactic modeling and limits on non-Gaussian statistics.

The Differential Microwave Radiometers (DMR) instrument aboard the Cosmic Background Explorer (COBE) has mapped the full microwave sky to mean sensitivity 26 uK per 7 deg field of view. The absolute calibration is determined to 0.7% with drifts smaller than 0.2% per year. We have analyzed both the raw differential data and the pixelized sky maps for evidence of contaminating sources such as solar system foregrounds, instrumental susceptibilities, and artifacts from data recovery and processing. Most systematic effects couple only weakly to the sky maps. The largest uncertainties in the maps result from the instrument susceptibility to the Earth's magnetic field, microwave emission from the Earth, and upper limits to potential effects at the spacecraft spin period. Systematic effects in the maps are small compared to either the noise or the celestial signal: the 95% confidence upper limit for the pixel-pixel rms from all identified systematics is less than 6 uK in the worst channel. A power spectrum analysis of the (A-B)/2 difference maps shows no evidence for additional undetected systematic effects.

The first two years of COBE DMR observations of the CMB anisotropy are analyzed and compared with our previously published first year results. The results are consistent, but the addition of the second year of data increases the precision and accuracy of the detected CMB temperature fluctuations. The two-year 53 GHz data are characterized by RMS temperature fluctuations of DT=44+/-7 uK at 7 degrees and DT=30.5+/-2.7 uK at 10 degrees angular resolution. The 53X90 GHz cross-correlation amplitude at zero lag is C(0)^{1/2}=36+/-5 uK (68%CL) for the unsmoothed 7 degree DMR data. A likelihood analysis of the cross correlation function, including the quadrupole anisotropy, gives a most likely quadrupole-normalized amplitude Q_{rms-PS}=12.4^{+5.2}_{-3.3} uK (68% CL) and a spectral index n=1.59^{+0.49}_{-0.55} for a power law model of initial density fluctuations, P(k)~k^n. With n fixed to 1.0 the most likely amplitude is 17.4 +/-1.5 uK (68% CL). Excluding the quadrupole anisotropy we find Q_{rms-PS}= 16.0^{+7.5}_{-5.2} uK (68% CL), n=1.21^{+0.60}_{-0.55}, and, with n fixed to 1.0 the most likely amplitude is 18.2+/-1.6 uK (68% CL). Monte Carlo simulations indicate that these derived estimates of n may be biased by ~+0.3 (with the observed low value of the quadrupole included in the analysis) and {}~+0.1 (with the quadrupole excluded). Thus the most likely bias-corrected estimate of n is between 1.1 and 1.3. Our best estimate of the dipole from the two-year DMR data is 3.363+/-0.024 mK towards Galactic coordinates (l,b)= (264.4+/-0.2 degrees, +48.1+/-0.4 degrees), and our best estimate of the RMS quadrupole amplitude in our sky is 6+/-3 uK.