SPheno

[ Latest version | manual for version 2.3 (pdf-file) | manual for additional features of version 3.1 (pdf-file) ]

The current version is 3.2.1 (10 October 2012) which can be found at the new SPheno webpage; for older version see below. This webpage will also be transfered to the new site in the next weeks.

SPheno stands for S(upersymmetric) Pheno(menology). The code calculates the SUSY spectrum using low energy data and a user supplied high scale model as input. The spectrum is used to calculate two- and three body decay modes of supersymmetric particle as well as of Higgs bosons. In addition the production cross sections for supersymmetric particle and Higgs bosons in e^+ e^- annihilation is calculated. Moreover, the branching of the decay $b \to s \gamma$, the SUSY contribution to anomalous magnetic moment of the muon as well as the SUSY contributions to the rho parameter due to sfermions are calculated. The code is written in F90 with an emphasis on easy generalisability. The structure is set such that complex phases as well as the extension to include the flavour structure can be done in a straight forward way. The 2-loop renormalization group equations as well as the one-loop finite corrections a la Bagger, Matchev, Pierce and Zhang are included. In addition the two-loop corrections to the neutral Higgs boson masses (a la Brignole, Degrassi, Slavich and Zwirner) and to the mu-parameter (a la Dedes and Slavich) are included. Starting with version 2.2.2 the SUSY Les Houches Accord is supported as well as the SPA conventions (for details see hep-ph/0511344). Starting with version 3.0.beta the complete flavour and CP structure of the MSSM as well as bilinear R-parity violation has been implemented. The input and output is done via he SUSY Les Houches Accord 2 .

If you use SPheno to write a paper, please cite W. Porod, Comput. Phys. Commun. 153 (2003) 275 [arXiv:hep-ph/0301101] and W. Porod and F. Staub arXiv:1104.1573 , which are the SPheno manuals. The version on the electronic hep-ph archive will be updated regularly with more recent versions.

Detailed comparisons between SPheno and other public codes have been performed in:

Precise determination of the neutral Higgs boson masses in the MSSM by B. C. Allanach, A. Djouadi, J. L. Kneur, W. Porod and P. Slavich, JHEP 0409 (2004) 04 [hep-ph/0406166]. The input files for the mh_max scenario are HighScale.in (SUSY parameters) and StandardModel.in (SM parameters). The resulting output file is SPheno.out.
Theoretical uncertainties in sparticle mass predictions from computational tools by B. Allanach, S. Kraml and W. Porod JHEP 0303 (2003) 016 [hep-ph//0302102]
Comparison of SUSY mass spectrum calculations by B. Allanach, S. Kraml and W. Porod, hep-ph/0207314
Online comparison of various codes maintained by S. Kraml. (Note that the various programs have in general slightly different Standard Model input except for the top-mass.)

Please send an email to porod@physik.uni-wuerzburg.de with the subject 'subscribe SPheno' if you want to be automatically informed on new releases.

Installation

To install SPheno, you will need a F95 (F90) compiler. No additional library is necessary. The code has been successfully compiled using NAG F95, Lahey/Fujitsu lf95, GNU g95, and Intel ifort (version 10.1) on linux PC system [from this list, the g95-compiler is completely free]. In case of the gfortran compiler depedning on the compiler version and platform deviations from the results of the other compilers, which agree among themselves perfectly, have been found up to a few per-cent. The code is supposed to be standard ANSI compatible F95.

The present release contains the main program SPheno.f90 and the Code for the SPheno library (libSPheno.a, link with -lSPheno). A Makefile is also included: the default compiler is ifort, but the options are now also included for other compilers. In linux, just unpack the files with
> gunzip SPheno3.1.4.tar.gz
> tar -xvf SPheno3.1.4.tar
then type (if necessary change the definitions for the f90 compiler)
> cd SPheno3.1.4
> make (or gmake)
The program is then run by
> bin/SPheno
The output file SPheno.spc should be identical (except for small numerical deviations depending on the compiler) to SPheno.spc (SLHA format).

Release 3.1.4: 14.07.11

gzipped tarred source, manual (ps-file) , manual addendum (ps-file)
input/output according to the SUSY Les Houches Accord and SUSY Les Houches Accord 2
supports the conventions of Supersymmetry Parameter Analysis (SPA) project
implemented high scale models: mSUGRA, GMSB, AMSB, various string effective field theories, general interface for abitrary boundary conditions at the GUT scale. In all models but the one with general boundary conditions at the GUT scale one can switch on flavour violation.
implemented models at the electroweak scale: general MSSM and bilinear R-parity violation.

New features

1) Implementation of Seesaw type I, II and III models. In case you want to use type II and/or type III you have to change the compiler options in the Makefile of the src directory
2) Implementation of rare lepton decays and additional B-physics observables
3) Enhanced speed in case of difficult scenarios
4) New blocks in the output which serve as input for the program HiggsBounds
5) taking into account all Yukawas in for the running of tan(beta)
6) improving the routine which determines the R-parity violating parameters when fitting the neutrino data
7) in case of SPA conventions one has to use entry 25 of EXTPAR to specify tan(beta) at m_Z and not as up to now entry 3 of MINPAR where tan(beta) is given at m_Z

Solved bugs (14.07.2011)

1) wrong value for mu in the output
2) in case of seesaw I models the manual contained the wrong ordering of the indices of the Yukawa couplings
3) correcting an error in the routine for calculating the non-universal contributions to the gauge couplings and m_W at m_Z in case of generation mixing and phases in the lepton/slepton sector; can lead to numerical differences which are significant
4) in case of flavour violating 3-body decays of charginos a bug has been fixed for the final states t q neutralino_j, where q=d,s

Known problems (24.05.2011)

1) Resummation of tan(beta) enhanced corrections to the bottom und tau Yukawa couplings cannot be switched of
2) Some features of the SLHA conventions are not yet implemented, see appendix B of the new manual.
3) Bugs have been reported for sparticle and Higgs production in e+ e- annihilation in case of R-parity violation; the corresponding calculations are switched of..

Model files plus the corresponding output files (default values of the program are used). The original SPheno input/output is no longer supported and everything is done using SLHA conventions, e.g.:

[ LesHouches.in.mSugra | SPheno.spc.mSugra ] minimal SUGRA model
[ LesHouches.in.GMSB | SPheno.spc.GMSB ] gauge mediated SUSY breaking
[ LesHouches.in.AMSB | SPheno.spc.AMSB ] anomaly mediated SUSY breaking

Pre-compiled versions

Due to increasing problems in the past, there are currently no pre-compiled versions.