Search

Abstract We present PEANUTS (Propagation and Evolution of Active NeUTrinoS), an open-source Python package for the automatic computation of solar neutrino spectra and active neutrino propagation through Earth. PEANUTS is designed to be fast, by employing analytic formulae for the neutrino propagation through varying matter density, and flexible, by allowing the user to input arbitrary solar models, custom Earth density profiles and general detector locations. It provides functionalities for a fully automated simulation of solar neutrino fluxes at a detector, as well as access to individual routines to perform more specialised computations. The software has been extensively tested against the results of the SNO experiment, providing excellent agreement with their results. In addition, the present text contains a pedagogical derivation of the relations needed to compute the oscillated solar neutrino spectra, neutrino propagation through Earth and nadir exposure of an experiment.

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

Abstract 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}}}$$ G ~ -EWMSSM, exhibits quite varied collider phenomenology due to its many permitted electroweakino production processes and decay modes. Characteristic $${{\tilde{G}}}$$ 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}}}$$ 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}}}$$ G ~ -EWMSSM parameter samples and GAMBIT input files generated for this work is available via Zenodo.

Abstract We present a likelihood analysis of the general two Higgs doublet model, using the most important currently measured flavour observables, in view of the anomalies in charged current tree-level and neutral current one-loop rare decays of B mesons in b → cl ν ¯ $$ \overline{\nu} $$ and b → sμ + μ − transitions, respectively. We corroborate that the model explains the latter and it is able to simultaneously fit the experimental values of the R(D) charged current ratio at 1σ, but it can not accommodate the D * charmed meson observables R(D *) and F L (D *). We find that the fitted values for the angular observables in b → sμ + μ − transitions exhibit better agreement with the general two Higgs double model in comparison to the SM. We also make predictions for future collider observables BR(t → ch), BR(h → bs), BR(h → τμ), BR(B s → τ + τ − ), BR(B + → K + τ + τ − ) and the flavour violating decays of the τ lepton, BR(τ → 3μ) and BR(τ → μγ). The model predicts values of BR(t → ch), BR(B s → τ + τ − ) and BR(B + → K + τ + τ − ) that are out of reach of future experiments, but its predictions for BR(h → bs) and BR(h → τμ) are within the future sensitivity of the HL-LHC or the ILC. We also find that the predictions for the τ → 3μ and τ → μγ decays are well within the projected limits of the Belle II experiment. Finally, using the latest measurement of the Fermilab Muon g − 2 Collaboration, we performed a simultaneous fit to ∆a μ constrained by the charged anomalies, finding solutions at the 1σ level. Once the neutral anomalies are included, however, a simultaneous explanation is unfeasible.

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

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

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

Abstract The extension of the Standard Model by right-handed neutrinos can not only explain the active neutrino masses via the seesaw mechanism, it is also able solve a number of long standing problems in cosmology. Especially, masses below the TeV scale are of particular interest as they can lead to a plethora of signatures in experimental searches. We present the first full frequentist analysis of the extension of the Standard Model by three right-handed neutrinos, with masses between 60 MeV and 500 GeV, using the Global and Modular BSM (beyond the Standard Model) Inference Tool GAMBIT. Our analysis is based on the Casas-Ibarra parametrisation and includes a large range of experimental constraints: active neutrino mixing, indirect constraints from, e.g., electroweak precision observables and lepton universality, and numerous direct searches for right-handed neutrinos. To study their overall effect, we derive combined profile likelihood results for the phenomenologically most relevant parameter projections. Furthermore, we discuss the role of (marginally) statistically preferred regions in the parameter space. Finally, we explore the flavour mixing pattern of the three right-handed neutrinos for different values of the lightest neutrino mass. Our results comprise the most comprehensive assessment of the model with three right-handed neutrinos model below the TeV scale so far, and provide a robust ground for exploring the impact of future constraints or detections.

Abstract 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 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 = (8–155, 103–260, 130–473, 219–502) GeV and chargino masses of = (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 $$ 3.3σ , 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.

Abstract We analyse the gravitational wave and low energy signatures of a Pati-Salam phase transition. For a Pati-Salam scale of M PS ∼ 105 GeV, we find a stochastic power spectrum within reach of the next generation of ground-based interferometer experiments such as the Einstein Telescope, in parts of the parameter space. We study the lifetime of the proton in this model, as well as complementarity with low energy constraints including electroweak precision data, neutrino mass measurements, lepton flavour violation, and collider constraints.

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

Abstract In Ref. (GAMBIT Collaboration: Athron et. al., Eur. Phys. J. C. arXiv:1705.07908, 2017) we introduced the global-fitting framework GAMBIT. In this addendum, we describe a new minor version increment of this package. GAMBIT 1.1 includes full support for Mathematica backends, which we describe in some detail here. As an example, we backend SUSYHD (Vega and Villadoro, JHEP 07:159, 2015), which calculates the mass of the Higgs boson in the MSSM from effective field theory. We also describe updated likelihoods in PrecisionBit and DarkBit, and updated decay data included in DecayBit.

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

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

Abstract In this paper, we embed the model of flipped GUT sneutrino inflation — in a flipped SU(5) or SO(10) set up — developed by Ellis et al. in a supergravity framework. The GUT symmetry is broken by a waterfall which could happen at early or late stage of the inflationary period. The full field dynamics is thus studied in detail and these two main inflationary configurations are exposed, whose cosmological predictions are both in agreement with recent astrophysical measurements. The model has an interesting feature where the inflaton has natural decay channels to the MSSM particles allowed by the GUT gauge symmetry. Hence it can account for the reheating after the inflationary epoch.