micrOMEGAs: Version 1.3 Original Research Article

We present the latest version of micrOMEGAs, a code that calculates the relic density of the lightest supersymmetric particle (LSP) in the minimal supersymmetric standard model (MSSM). All tree-level processes for the annihilation of the LSP are included as well as all possible coannihilation processes with neutralinos, charginos, sleptons, squarks and gluinos. The cross-sections extracted from CalcHEP are calculated exactly using loop-corrected masses and mixings as specified in the SUSY Les Houches Accord. Relativistic formulae for the thermal average are used and care is taken to handle poles and thresholds by adopting specific integration routines. The input parameters can be either the soft SUSY parameters in a general MSSM or the parameters of a SUGRA model specified at the GUT scale. In the latter case, a link with Suspect, SOFTSUSY, Spheno and Isajet allows one to calculate the supersymmetric spectrum, Higgs masses, as well as mixing matrices. Higher-order corrections to Higgs couplings to quark pairs including QCD as well as some SUSY corrections (image) are implemented. Routines calculating image, image and image are also included. In particular the image routine includes an improved NLO for the SM and the charged Higgs while the SUSY large image effects beyond leading-order are included. This new version also provides cross-sections for any image process as well as partial decay widths for two-body final states in the MSSM allowing for easy simulation at colliders.