|
|
SPheno
This is the development page for the SPheno project.
The current version
is 3.1.10 (4 May 2012) (Download); for older version see the
old SPheno webpage.
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.
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 . Moreover, starting with version 3.1.10 also parts of
the proposal to include seesaw models has been implemented, see the corresponding
chapter of the
Les Houches Proceedings 2011.
Detailed comparisons between SPheno and other public codes have been performed
in:
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 11.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 SPheno-3.1.10.tar.gz
> tar -xvf SPheno-3.1.10.tar
then type (if necessary change the definitions for the f90 compiler)
> cd SPheno-3.1.10
> 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).
The compilation of the RGEs for the seesaw tpye II and type III models
is rather time consuming and can take up to one hour in the optimized mode.
For this reason they are not included in the compilation using the
default options. If you want to include, include the option -DSEESAWIII
in the Makefile that is located in the src directory. In this way
both set of RGEs are included.
- 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
- 8) Three body decays of sleptons via virtual neutralinos have been added
- 9) Including the possibility to give within mSUGRA scenarios
m_A and mu at the electroweak scale. An example for a
corresponding input file can be found in the directory
input and is called LesHouches.in.mNUHM
file is called. An example where the Higgs mass parameters
are non-universal is given in LesHouches.in.mSUGRA_NUHM
- 10) Implementing the SLHA proposal for seesaw scenarios as given
in arXiv:1203.1488 (Les Houches 2011 proceedings)
- 11) Improving speeds of the RGE calculation in case of complex parameters,
this is in
particular relevant for the seesaw models.
- Some features of the SLHA conventions are not yet implemented, see
appendix B of the new manual.
- 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.
- Correcting input for low energy parameters in case of bilinear
R-parity violation.
- Correcting a bug in case of Z decays into two charged leptons with different
flavour.
- Correcting a bug in case of the Higgs-mass calculation
related to the use of the loop-corrected mu-parameter. This can lead to shifts of
up to several GeV for m_h in case of large A_t (at the electroweak scale).
The heavier Higgs bosons are hardly affected. The standard examples are hardly
affected.
- Correcting a bug in case of neutrino mass calculation for bilinear R-parity
violation: in some cases the seesaw formula had been implemented incorrectly
taking only the real part of the inverses neutralino mass matrix.
- In case of complex phases there has been a missing relative
complex conjugation between the A-parameter and the mu-parameter
in the LR-sector in case of generation mixing. This affects mainly
the EDM calculation but leaves the spectrum essentially unchanged.
- Correcting a bug in case of input parameters are the
electroweak scale and where on-shell mass for the
pseudoscalar Higgs boson is used as input.
- 1) the values of EXTPAR are also given in SPheno.spc
- 2) the leptonic EDMs had been too small by a factor 100 due to a typo
in an overall factor. Moreover, there has been a convention problem
in the sfermion mass problems which mainly affects the electron EDM
and only slightly the other EDMs and at most at the per-mille the
Delta(M_Bd) and Delta(M_Bs).
- 1) input for left slepton mass squared parameters fixed
- 2) fixing numerical problem in expansion for loop function DB_0
- 1) a problem related to the input of slepton mass parameters in the
super-PMNS basis at the electroweak scalae has been fixed.
- 2) a bug in the charged Higgs boson contribution to the fermion
self-energies has been fixed, little impact on the numerical results
- 3) the coefficient of the U(1) coupling in the running of the dim-5
neutrino mass operator has been corrected, small numerical impact
- 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
|