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] adqr_v4_0.tar.gz(9577 Kbytes)
Manuscript Title: MicrOMEGAs4.1: Two dark matter candidates
Authors: G. Bélanger, F. Boudjema, A. Pukhov, A. Semenov
Program title: MicrOMEGAs4.1
Catalogue identifier: ADQR_v4_0
Distribution format: tar.gz
Journal reference: Comput. Phys. Commun. 192(2015)322
Programming language: C and Fortran.
Computer: PC, Mac.
Operating system: UNIX (Linux, Darwin).
RAM: 50MB depending on the number of processes required.
Keywords: Dark matter, Relic density, Indirect detection, MSSM, Beyond Standard Model.
PACS: 95.35.+d, 14.80.Ly, 12.60.Jv.
Classification: 1.9, 11.6.

External routines: CalcHEP, SuSpect, NMSSMTools, CPSuperH, LoopTools, HiggsBounds

Does the new version supersede the previous version?: Yes

Nature of problem:
Calculation of the relic density and direct and indirect detection rates of the lightest stable particle in particle physics models with at most two stable dark matter candidates.

Solution method:
In the case where the two dark matter particles have very different masses, we find that the equations for the evolution of the density of dark matter behave as stiff equations. To solve these we use the backward scheme and the Rosenbrock algorithm. The standard solution based on the Runge-Kutta method is still used for models with only one dark matter candidate.

Reasons for new version:
There are many experiments that are currently searching for the remnants of dark matter annihilation and the relic density is determined precisely from cosmological measurements. In this version we generalize the Boltzmann equations to take into account the possibility of two dark matter candidates. Thus, in solving for the relic density we include interactions between the two dark matter sectors as well as semi-annihilation. The dark matter signals in direct and indirect detection are computed as well.

Summary of revisions:
  • Generalization of the Boltzmann equations to include two dark matter candidates, their interactions and semi-annihilations, the relative density of the two dark matter components is taken into account when computing direct/indirect detection rates.
  • Upgrade of the numerical method for solving the Boltzmann equations
  • Include sample extensions of the standard model with extra doublet and singlets which contain two stable neutral particles

Unusual features:
Depending on the parameters of the model, the program generates additional new code, compiles it and loads it dynamically.

Running time:
4 sec