Computer Physics Communications Program Library

Programs in Physics & Physical Chemistry



Manuscript Title: SuSpect: a Fortran Code for the Supersymmetric and Higgs Particle Spectrum in the MSSM
Authors: Abdelhak Djouadi, Jean-Loïc Kneur, Gilbert Moultaka
Program title: SuSpect
Catalogue identifier: ADYR_v1_0 Distribution format: tar.gz
Journal reference: Comput. Phys. Commun. 176(2007)426
Programming language: FORTRAN 77.
Computer: Unix machines, PC.
Operating system: Unix (or Linux).
RAM: approximately 2500 Kbytes
Keywords: Supersymmetry, renormalization group equations, loop calculations.
PACS: 12.60.Jv, 11.30.Qc, 12.10.-g, 12.15.Lk, 14.80.Cp, 14.80.Ly.
Classification: 11.6, 4.3, 4.12.
Nature of problem: SuSpect calculates the supersymmetric and Higgs particle spectrum (masses and some other relevant parameters) in the unconstrained Minimal Supersymmetric Standard Model (MSSM), as well as in constrained models (cMSSMs) such as the minimal Supergravity (mSUGRA), the gauge mediated (GMSB) and anomaly mediated (AMSB) Supersymmetry breaking scenarii. The following features and ingredients are included: renormalization group evolution between low and high energy scales, consistent implementation of radiative electroweak symmetry breaking, calculation of the physical particle masses with radiative corrections at the one- and two-loop level.
Solution method: The main methods used in the code are: an (adaptative fourth-order) Runge-Kutta type algorithm, (following a standard algorithm described in ``Numerical Recipes"), used to solve numerically a set of coupled differential equations resulting from the renormalization group equations at the two-loop level of the perturbative expansions; diagonalizations of mass matrices; some mathematical (Spence etc) functions resulting from the evaluation of one and two-loop integrals using the Feynman graphs techniques for radiative corrections to the particle masses; finally some fixed-point iterative algorithms to solve non-linear equations for some of the relevant output parameters.
Restrictions: The code is limited at the moment to real input parameters. It also does not include flavor non-diagonal terms which are possible in the most general soft supersymmetry breaking Lagrangian. There are some (mild) limitations on the possible range of values of input parameter, i.e. not any arbitrary values of some input parameters are allowed: these limitations are essentially based on physical rather than algorithmic issues, and warning flags and other protections are installed to avoid as much as possible execution failure if unappropriate input values are used.
Running time: between 1 and 3 seconds depending on options, with a 1 GHz processor.