MicrOMEGAs4.1: Two dark matter candidates

Abstract

micrOMEGAs is a code to compute dark matter observables in generic extensions of the standard model. This version of micrOMEGAs includes a generalization of the Boltzmann equations to take into account the possibility of two dark matter candidates. The modification of the relic density calculation to include interactions between the two dark matter sectors as well as semi-annihilation is presented. Both dark matter signals in direct and indirect detection are computed as well. An extension of the standard model with two scalar doublets and a singlet is used as an example. Program summaryProgram title: MicrOMEGAs4.1Catalogue identifier: ADQR_v4_0Program summary URL:http://cpc.cs.qub.ac.uk/summaries/ADQR_v4_0.htmlProgram obtainable from: CPC Program Library, Queen’s University, Belfast, N. IrelandLicensing provisions: Standard CPC licence, http://cpc.cs.qub.ac.uk/licence/licence.htmlNo. of lines in distributed program, including test data, etc.: 738425No. of bytes in distributed program, including test data, etc.: 9807620Distribution format: tar.gzProgramming language: C and Fortran.Computer: PC, Mac.Operating system: UNIX (Linux, Darwin).RAM: 50MB depending on the number of processes required.Classification: 1.9, 11.6.Catalogue identifier of previous version: ADQR_v3.0Journal reference of previous version: Comput. Phys. Comm. 185 (2014) 960External routines: CalcHEP, SuSpect, NMSSMTools, CPSuperH, LoopTools, HiggsBoundsDoes the new version supersede the previous version?: YesNature 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