Source author record

Peter Athron

Peter Athron appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

29works
4topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

29 published item(s)

preprint2022arXiv

Global fits of SUSY at future Higgs factories

In this work, we study the impact of electroweak and Higgs precision measurements at future electron-positron colliders on several typical supersymmetric models, including the Constrained Minimal Supersymmetric Standard Model (CMSSM), Non-Universal Higgs Mass generalisations (NUHM1, NUHM2), and the 7-dimensional Minimal Supersymmetric Standard Model (MSSM7). Using publicly-available data from the \textsf{GAMBIT} community, we post-process previous SUSY global fits with additional likelihoods to explore the discovery potential of Higgs factories, such as the Circular Electron Positron Collider (CEPC), the Future Circular Collider (FCC) and the International Linear Collider (ILC). We show that the currently allowed parameter space of these models will be further tested by future precision measurements. In particular, dark matter annihilation mechanisms may be distinguished by precise measurements of Higgs observables.

preprint2022arXiv

GM2Calc-2 for the 2HDM

GM2Calc is a leading tool for calculating precise contributions to $a_μ$ in the Minimal Supersymmetric Standard Model. In this proceeding we detail GM2Calc version 2 where it is extended so it can calculate two-loop contributions to $a_μ$ in the Two-Higgs Doublet Model (2HDM), based on the work in Ref. [1]. The 2HDM is a simple model, yet it is one of the few single field extensions of the Standard Model which is able to explain the muon $g-2$ anomaly. We demonstrate the powerful and flexible 2HDM capabilities of GM2Calc2, which include the most precise contributions in the literature and allow the user to work in their favourite type of the 2HDM as well as use complex and lepton flavour violating couplings. With its multiple interfaces and input flexibility, GM2Calc2 is a powerful tool both as a standalone code and as part of a larger code toolchain.

preprint2022arXiv

Likelihood analysis of the flavour anomalies and $g-2$ in the general two Higgs doublet model

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\to c l \overlineν$ and $b\to 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)$ ratio at $1σ$, but it can not accommodate the $R(D^{*})$ and $F_{L}(D^{*})$ observables. We find that the fitted values for the angular observables in $b\to 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 $\mathrm{BR}(t\to ch)$, $\mathrm{BR}(h\to bs)$, $\mathrm{BR}(h\to τμ)$, $\mathrm{BR}(B_{s}\rightarrowτ^{+}τ^{-})$, $\mathrm{BR}(B^{+}\rightarrow K^{+}τ^{+}τ^{-})$ and the flavour violating decays of the $τ$ lepton, $\mathrm{BR}(τ\rightarrow3μ)$ and $\mathrm{BR}(τ\toμγ)$. The model predicts values of $\mathrm{BR}(t\to ch)$, $\mathrm{BR}(B_{s}\rightarrowτ^{+}τ^{-})$ and $\mathrm{BR}(B^{+}\rightarrow K^{+}τ^{+}τ^{-})$ that are out of reach of future experiments, but its predictions for $\mathrm{BR}(h\to bs)$ and $\mathrm{BR}(h\to τμ)$ are within the future sensitivity of the HL-LHC or the ILC. We also find that the predictions for the $τ\rightarrow3μ$ and $τ\toμγ$ 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.

preprint2022arXiv

Simple and statistically sound recommendations for analysing physical theories

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.

preprint2022arXiv

Two-loop Prediction of the Anomalous Magnetic Moment of the Muon in the Two-Higgs Doublet Model with GM2Calc 2

We present an extension of the GM2Calc software to calculate the muon anomalous magnetic moment ($a_μ^{\text{BSM}}$) in the Two-Higgs Doublet Model. The Two-Higgs Doublet Model is one of the simplest and most popular extensions of the Standard Model. It is one of the few single field extensions that can give large contributions to $a_μ^{\text{BSM}}$. It is essential to include two-loop corrections to explain the long standing discrepancy between the Standard Model prediction and the experimental measurement in the Two-Higgs Doublet Model. The new version GM2Calc 2 implements the state of the art two-loop calculation for the general, flavour violating Two-Higgs Doublet Model as well as for the flavour aligned Two-Higgs Doublet Model and the type I, II, X and Y flavour conserving variants. Input parameters can be provided in either the gauge basis or the mass basis, and we provide an easy to use SLHA-like command-line interface to specify these. Using this interface users may also select between Two-Higgs Doublet Model types and choose which contributions to apply. In addition, GM2Calc 2 also provides interfaces in C++, C, Python and Mathematica, to make it easy to interface with other codes.

preprint2021arXiv

The GAMBIT Universal Model Machine: from Lagrangians to Likelihoods

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.

preprint2021arXiv

Thermal WIMPs and the Scale of New Physics: Global Fits of Dirac Dark Matter Effective Field Theories

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.

preprint2020arXiv

Global analyses of Higgs portal singlet dark matter models using GAMBIT

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.

preprint2020arXiv

PhaseTracer: tracing cosmological phases and calculating transition properties

We present a C++ software package called PhaseTracer for mapping out cosmological phases, and potential transitions between them, for Standard Model extensions with any number of scalar fields. PhaseTracer traces the minima of effective potential as the temperature changes, and then calculates the critical temperatures, at which the minima are degenerate. PhaseTracer is constructed with modularity, flexibility and practicality in mind. It is fast and stable, and can receive potentials provided by other packages such as FlexibleSUSY. PhaseTracer can be useful analysing cosmological phase transitions which played an important role in the very early evolution of the Universe. If they were first order they could generate detectable gravitational waves and/or trigger electroweak baryogenesis to generate the observed matter anti-matter asymmetry of the Universe. The code can be obtained from https://github.com/PhaseTracer/PhaseTracer.

preprint2020arXiv

Reinterpretation of LHC Results for New Physics: Status and Recommendations after Run 2

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.

preprint2019arXiv

Combined collider constraints on neutralinos and charginos

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χ_1^0}, m_{\tildeχ_2^0}, m_{\tildeχ_3^0}, m_{\tildeχ_4^0})$ = (8-155, 103-260, 130-473, 219-502) GeV and chargino masses of $(m_{\tildeχ_1^{\pm}}, m_{\tildeχ_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σ$, 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$σ$. 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.

preprint2019arXiv

New Physics Interpretations with GAMBIT

I present recent results from the Global and Modular Beyond-the-Standard-Model Inference Tool (GAMBIT) collaboration. Global fits with GAMBIT have been carried out on a variety of models including supersymmetric models, scalar singlet dark matter, fermionic and vector Higgs portal dark matter and axions. Here I focus on a recent GAMBIT study interpreting collider constraints on electroweakinos (arXiv:1809.02097). We show that when the neutralinos and charginos are the only light states of the MSSM, there are scenarios which evade LHC constraints for any mass of the lightest neutralino and the lightest chargino, i.e. the profile likelihood shows no constraint in this plane when one only considers the possibility of excluding new physics. Intriguingly, in addition we also find that excesses in the data can lead to closed confidence level contours, indicating a preference for light neutralinos and charginos over the Standard Model. We find the excess has a local significance of 3.3 sigma when combining ATLAS and CMS 13 TeV searches, but this drops to 2.9 sigma when including 8 TeV searches as well.

preprint2019arXiv

SpecBit, DecayBit and PrecisionBit: GAMBIT modules for computing mass spectra, particle decay rates and precision observables

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.

preprint2019arXiv

Strong first-order phase transitions in the NMSSM --- a comprehensive survey

Motivated by the fact that the Next-to-Minimal Supersymmetric Standard Model is one of the most plausible models that can accommodate electroweak baryogenesis, we analyze its phase structure by tracing the temperature dependence of the minima of the effective potential. Our results reveal rich patterns of phase structure that end in the observed electroweak symmetry breaking vacuum. We classify these patterns according to the first transition in their history and show the strong first-order phase transitions that may be possible in each type of pattern. These could allow for the generation of the matter-antimatter asymmetry or potentially observable gravitational waves. For a selection of benchmark points, we checked that the phase transitions completed and calculated the nucleation temperatures. We furthermore present samples that feature strong first-order phase transitions from an extensive scan of the whole parameter space. We highlight common features of our samples, including the fact that the Standard Model like Higgs is often not the lightest Higgs in the model.

preprint2018arXiv

A global fit of the MSSM with GAMBIT

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.

preprint2018arXiv

Global fits of GUT-scale SUSY models with 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 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.

preprint2018arXiv

Impact of vacuum stability, perturbativity and XENON1T on global fits of $\mathbb{Z}_2$ and $\mathbb{Z}_3$ scalar singlet dark matter

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

preprint2018arXiv

Status of the scalar singlet dark matter model

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.

preprint2017arXiv

GAMBIT: The Global and Modular Beyond-the-Standard-Model Inference Tool

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.

preprint2017arXiv

SUSY Global Fits

We present comprehensive global fits of supersymmetric (SUSY) models from the Global and Modular Beyond-the-Standard-Model Inference Tool (GAMBIT) collaboration, based on arXiv:1705.07935 and arXiv:1705.07917. We investigate several variants of the minimal supersymmetric standard model, a fully constrained version (CMSSM) with universal scalar ($m_0$), gaugino ($m_{1/2}$) and trilinear masses ($A_0$) at the gauge coupling unification scale, a similar model that is relaxed by adding an extra parameter for the soft Higgs masses (NUHM1), another where the soft Higgs masses are also split (NUHM2) and finally a weak scale MSSM7 model. We use the public GAMBIT global fitting framework and take into account all relevant data to reveal the regions of parameter space with the highest likelihood. Our results reveal that all models have very heavy scenarios that are well out of reach of the LHC, but will be probed by forthcoming dark matter experiments, as well as a stop-co-annihilation region which has better prospects for detection in collider experiments. The stau co-annihilation region is excluded from the CMSSM at $2 σ$ but is present in the NUHM1 and NUHM2 variants. Finally by relaxing constraints in the NUHM1, NUHM2 and MSSM7 we see additional regions appear: lighter chargino co-annihilation region, sbottom co-annihilation and $h/Z$ funnels.

preprint2016arXiv

Higgs mass predictions of public NMSSM spectrum generators

The publicly available spectrum generators for the NMSSM often lead to different predictions for the mass of the standard model-like Higgs boson even if using the same renormalization scheme and two-loop accuracy. Depending on the parameter point, the differences can exceed 5 GeV, and even reach 8 GeV for moderate superparticle masses of up to 2 TeV. It is shown here that these differences can be traced back to the calculation of the running standard model parameters entering all calculations, to the approximations used in the two-loop corrections included in the different codes, and to different choices for the renormalization conditions and scales. In particular, the importance of the calculation of the top Yukawa coupling is pointed out.

preprint2016arXiv

Precision tools and models to narrow in on the 750 GeV diphoton resonance

The hints for a new resonance at 750 GeV from ATLAS and CMS have triggered a significant amount of attention. Since the simplest extensions of the standard model cannot accommodate the observation, many alternatives have been considered to explain the excess. Here we focus on several proposed renormalisable weakly-coupled models and revisit results given in the literature. We point out that physically important subtleties are often missed or neglected. To facilitate the study of the excess we have created a collection of 40 model files, selected from recent literature, for the Mathematica package SARAH. With SARAH one can generate files to perform numerical studies using the tailor-made spectrum generators FlexibleSUSY and SPheno. These have been extended to automatically include crucial higher order corrections to the diphoton and digluon decay rates for both CP-even and CP-odd scalars. Additionally, we have extended the UFO and CalcHep interfaces of SARAH, to pass the precise information about the effective vertices from the spectrum generator to a Monte-Carlo tool. Finally, as an example to demonstrate the power of the entire setup, we present a new supersymmetric model that accommodates the diphoton excess, explicitly demonstrating how a large width can be obtained. We explicitly show several steps in detail to elucidate the use of these public tools in the precision study of this model.

preprint2015arXiv

FlexibleSUSY -- A spectrum generator generator for supersymmetric models

We introduce FlexibleSUSY, a Mathematica and C++ package, which generates a fast, precise C++ spectrum generator for any SUSY model specified by the user. The generated code is designed with both speed and modularity in mind, making it easy to adapt and extend with new features. The model is specified by supplying the superpotential, gauge structure and particle content in a SARAH model file; specific boundary conditions e.g. at the GUT, weak or intermediate scales are defined in a separate FlexibleSUSY model file. From these model files, FlexibleSUSY generates C++ code for self-energies, tadpole corrections, renormalization group equations (RGEs) and electroweak symmetry breaking (EWSB) conditions and combines them with numerical routines for solving the RGEs and EWSB conditions simultaneously. The resulting spectrum generator is then able to solve for the spectrum of the model, including loop-corrected pole masses, consistent with user specified boundary conditions. The modular structure of the generated code allows for individual components to be replaced with an alternative if available. FlexibleSUSY has been carefully designed to grow as alternative solvers and calculators are added. Predefined models include the MSSM, NMSSM, E$_6$SSM, USSM, R-symmetric models and models with right-handed neutrinos.

preprint2015arXiv

GM2Calc: Precise MSSM prediction for $(g - 2)$ of the muon

We present GM2Calc, a public C++ program for the calculation of MSSM contributions to the anomalous magnetic moment of the muon, $(g-2)_μ$. The code computes $(g-2)_μ$ precisely, by taking into account the latest two-loop corrections and by performing the calculation in a physical on-shell renormalization scheme. In particular the program includes a $\tanβ$ resummation so that it is valid for arbitrarily high values of $\tanβ$, as well as fermion/sfermion-loop corrections which lead to non-decoupling effects from heavy squarks. GM2Calc can be run with a standard SLHA input file, internally converting the input into on-shell parameters. Alternatively, input parameters may be specified directly in this on-shell scheme. In both cases the input file allows one to switch on/off individual contributions to study their relative impact. This paper also provides typical usage examples not only in conjunction with spectrum generators and plotting programs but also as C++ subroutines linked to other programs.

preprint2014arXiv

Bayesian naturalness of the CMSSM and CNMSSM

The recent discovery of the 125.5 GeV Higgs boson at the LHC has fueled interest in the next-to-minimal supersymmetric standard model (NMSSM) as it may require less fine-tuning than the minimal model to accommodate such a heavy Higgs. To this end we present Bayesian naturalness priors to quantify fine-tuning in the (N)MSSM. These priors arise automatically as Occam razors in Bayesian model comparison and generalize the conventional Barbieri-Giudice measure. In this paper we show that the naturalness priors capture features of both the Barbieri-Giudice fine-tuning measure and a simple ratio measure that has been used in the literature. We also show that according to the naturalness prior the constrained version of the NMSSM is less tuned than the CMSSM.

preprint2014arXiv

FlexibleSUSY - a meta spectrum generator for supersymmetric models

FlexibleSUSY is a software package that takes as input descriptions of (non-)minimal supersymmetric models written in Wolfram/Mathematica and generates a set of spectrum generator libraries and executables, with the aid of SARAH. The design goals are precision, reliability, modularity, speed, and readability of the code. The boundary conditions are independent C++ objects that are plugged into the boundary value problem solver together with the model objects. This clean separation makes it easy to adapt the generated code for individual projects. The current status of the interface and implementation is sketched.

preprint2013arXiv

Fine Tuning in the Constrained Exceptional Supersymmetric Standard Model

Supersymmetric unified models in which the Z' couples to the Higgs doublets, as in the E6 class of models, have large fine tuning dominated by the experimental mass limit on the Z'. To illustrate this we investigate the degree of fine tuning throughout the parameter space of the Constrained Exceptional Supersymmetric Standard Model (cE6SSM) that is consistent with a Higgs mass m_h ~ 125 GeV. Fixing tan beta = 10, and taking specific values of the mass of the Z' boson, with M_Z' ~ 2-4 TeV. We find that the minimum fine tuning is set predominantly from the mass of Z' and varies from ~ 200-400 as we vary M_Z' from ~ 2-4 TeV. However, this is significantly lower than the fine tuning in the Constrained Minimal Supersymmetric Standard Model (cMSSM), of O(1000), arising from the large stop masses required to achieve the Higgs mass.

preprint2012arXiv

Threshold Corrections in the Exceptional Supersymmetric Standard Model

We calculate threshold corrections to the running gauge and Yukawa couplings in the Exceptional Supersymmetric Standard Model (E6SSM) and analyse the more precise and reliable mass spectra in a constrained model (CE6SSM). Full expressions for the corrections are provided and the implementation into a spectrum generator is described. We find a dramatic reduction in the matching scale dependency of the masses of many states and observe a significant adjustment of the correlation of low-scale physical masses and high-scale parameters. Still, in substantial regions of parameter space the mass of the lightest Higgs is compatible with the new boson discovered at the LHC and the model satisfies limits from collider searches for squark, gluinos and Z' bosons. We study the implications for gauge coupling unification from a new dependency of the spectrum on so-called survival Higgs fields which cannot be addressed without the inclusion of the threshold corrections.