|
|
Athron, P., Bach, M., Fargnoli, H. G., Gnendiger, C., Greifenhagen, R., Park, J. H., et al. (2016). GM2Calc: precise MSSM prediction for (g – 2) of the muon. Eur. Phys. J. C, 76(2), 62–16pp.
Abstract: We present GM2Calc, a public C++ program for the calculation of MSSM contributions to the anomalous magnetic moment of the muon, (g-2)(mu). The code computes (g -2)(mu) 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 beta resummation so that it is valid for arbitrarily high values of tan beta, 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.
|
|
|
|
Athron, P., Park, J. H., Stockinger, D., & Voigt, A. (2015). FlexibleSUSY-A spectrum generator generator for supersymmetric models. Comput. Phys. Commun., 190, 139–172.
Abstract: 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, E6SSM, USSM, R-symmetric models and models with right-handed neutrinos. Program Summary Program title: FlexibleSUSY Catalogue identifier: AEVIv10 Program summary URL: http://cpc.cs.qub.ac.uk/summaries/AEVIv10.html obtainable from: CPC Program Library, Queen's University, Belfast, N. Ireland Licensing provisions: GNU General Public License, version 3 No. of lines in distributed program, including test data, etc.: 129406 No. of bytes in distributed program, including test data, etc.: 854831 Distribution format: tar.gz Programming language: C++, Wolfram/Mathematica, FORTRAN, Bourne shell. Computer: Personal computer. Operating system: Tested on Linux 3.x, Mac OS X. Classification: 11.1, 11.6, 6.5. External routines: SARAH 4.0.4, Boost library, Eigen, LAPACK Nature of problem: Determining the mass spectrum and mixings for any supersymmetric model. The generated code must find simultaneous solutions to constraints which are specified at two or more different renormalization scales, which are connected by renormalization group equations forming a large set of coupled first-order differential equations. Solution method: Nested iterative algorithm and numerical minimization of the Higgs potential. Restrictions: The couplings must remain perturbative at all scales between the highest and the lowest boundary condition. FlexibleSUSY assumes that all couplings of the model are real (i.e. CP-conserving). Due to the modular nature of the generated code, adaption and extension to overcome restrictions in scope is quite straightforward. Running time: 0.06-0.2 seconds per parameter point.
|
|
|
|
GAMBIT Collaboration(Athron, P. et al), & Ruiz de Austri, R. (2017). Global fits of GUT-scale SUSY models with GAMBIT. Eur. Phys. J. C, 77(12), 824–44pp.
Abstract: 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 darkmatter in all threemodels, 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.
|
|
|
|
Staub, F., Athron, P., Basso, L., Goodsell, M. D., Harries, D., Krauss, M. E., et al. (2016). Precision tools and models to narrow in on the 750 GeV diphoton resonance. Eur. Phys. J. C, 76(9), 516–57pp.
Abstract: 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.
|
|
