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.

For particle physics observables at colliders such as the LHC at CERN, it has been common practice for many decades to estimate the theoretical uncertainty by studying the variations of the predicted cross sections with a priori unpredictable scales. In astroparticle physics, this has so far not been possible, since most of the observables were calculated at Born level only, so that the renormalization scheme and scale dependence could not be studied in a meaningful way. In this paper, we present the first quantitative study of the theoretical uncertainty of the neutralino dark matter relic density from scheme and scale variations. We first explain in detail how the renormalization scale enters the tree-level calculations through coupling constants, masses and mixing angles. We then demonstrate a reduction of the renormalization scale dependence through one-loop SUSY-QCD corrections in many different dark matter annihilation channels and enhanced perturbative stability of a mixed on-shell/$\bar{\rm DR}$ renormalization scheme over a pure $\bar{\rm DR}$ scheme in the top-quark sector. In the stop-stop annihilation channel, the Sommerfeld enhancement and its scale dependence are shown to be of particular importance. Finally, the impact of our higher-order SUSY-QCD corrections and their scale uncertainties are studied in three typical scenarios of the phenomenological Minimal Supersymmetric Standard Model with eleven parameters (pMSSM-11). We find that the theoretical uncertainty is reduced in many cases and can become comparable to the size of the experimental one in some scenarios.

We analyse the prospects for constructing hybrid models of inflation that provide a dynamical realisation of the apparent closeness between the supersymmetric GUT scale and the possible scale of cosmological inflation. In the first place, we consider models based on the flipped SU(5)$\times$U(1) gauge group, which has no magnetic monopoles. In one model, the inflaton is identified with a sneutrino field, and in the other model it is a gauge singlet. In both cases we find regions of the model parameter spaces that are compatible with the experimental magnitudes of the scalar perturbations, $A_s$, and the tilt in the scalar perturbation spectrum, $n_s$, as well as with an indicative upper limit on the tensor-to-scalar perturbation ratio, $r$. We also discuss embeddings of these models into SO(10), which is broken at a higher scale so that its monopoles are inflated away.