Elsevier Science Home
Computer Physics Communications Program Library
Full text online from Science Direct
Programs in Physics & Physical Chemistry
CPC Home

[Licence| Download | New Version Template] aezl_v1_0.tar.gz(1188 Kbytes)
Manuscript Title: GravitinoPack and decays of supersymmetric metastable particles
Authors: Helmut Eberl, Vassilis C. Spanos
Program title: GravitinoPack
Catalogue identifier: AEZL_v1_0
Distribution format: tar.gz
Journal reference: Comput. Phys. Commun. 202(2016)310
Programming language: Fortran 77 and C for the MathLink.
Computer: Workstation, PC, Mac.
Operating system: Linux, Mac OSX.
Keywords: Supersymmetry, Gravitino Physics, BBN.
PACS: 11.30.Pb, 12.60.Jv, 04.65.+e, 14.80.Ly.
Classification: 11.1.

Nature of problem:
If Supersymmetry (SUSY) is realized in nature, each Standard Model (SM) particle is accompanied by its superpartner. The gravitino, the spin 3/2 superpartner of the graviton, belongs to the spectrum of the most important supersymmetric models. The gravitino can be the Lightest Supersymmetric Particle (LSP), in the so-called gravitino Dark Matter (DM) models. Thus, the Next to the Lightest Supersymmetric Particle (NLSP) is long-lived, since it decays gravitationally to the LSP (gravitino) and SM particles. Such decays produce electromagnetic energy and hadrons which can affect the primordial Big-Bang Nucleosynthesis (BBN) predictions for the abundance of the light nuclei. Similarly, when the gravitino is not the LSP, it decays into the LSP, which can be another SUSY particle, like the neutralino, and SM particles. GravitinoPack provides important results on all these decays, making the application of the BBN data more accurate.

Solution method:
GravitinoPack is a package for the evaluation of processes with gravitino interaction. The version GravitinoPack1.0 includes all two-body decays of the gravitino and all three-body decays of the gravitino to a neutralino and a pair of two particles. In the case that the gravitino is the LSP, all two- and three-body decays of the lightest neutralino, the light stop or stau NLSP are included. GravitinoPack works at Fortran77 level and has a Mathematica interface.

Unusual features:
MathLink works properly only up to version Mathematica 9. There is still a conflict between Mathematica 10 and the used integration routine CUBA.

Running time:
All two-body decays widths are calculated in a time far below one second. The three-body decay calculations can be fast, especially the stop and stau decays, or slow with times up to a few minutes, especially when the phase space becomes large and no approximation is chosen.