Micromegas
v_2.4 for a generic model for the calculation of e Relic density Direct detection rates Indirect detection rates
Code to calculate the properties of a stable massive particle in
a generic model. First developed to compute the relic density of
a stable massive particle, the code
also computes the rates for direct and indirect detection rates
of dark matter. It is assumed that a discrete symmetry like R-parity
ensures the stability of the lightest odd particle. All
annihilation and coannihilation channels are included in the
computation of the relic density. Specific examples of this general
approach include the MSSM and various extensions. Extensions to other
models can be implemented by the user. The New Physics model first
requires to write a new CalcHEP model file, a package for the automatic generation of squared matrix elements. This can be done through LanHEP. Once this is done, all annihilation and coannihilation channels are included automatically in any model.
The cross-sections for both spin dependent and spin independent
interactions of WIMPS on protons are computed automatically as well as
the rates for WIMP scattering on nuclei in a large detector.
Annihilation cross-sections of the dark matter candidate at zero
velocity, relevant for indirect detection of dark matter, are also
computed automatically. The propagation of charged particles in the Galactic halo is handled with a new module.
The package includes the minimal supersymmetric standard model (MSSM),
the NMSSM, the MSSM with complex phases (CPVMSSM), the little Higgs
model (LHM) and a model with right-handed neutrino DM (RHNM).
I) New options
1) The SM contribution in the bsgamma routine in the MSSM is given as an output
2) The value of the charm quark mass (McMc) is added as a parameter of the MSSM
II) Bugs removed
1) The annihilation cross section was repeated twice in the outut
2) A bug in the Spheno interface was removed
3) A mistake in the anti-proton propagator was corrected
4) In SHOWPLOTS, the electron spectrum was displayed instead of the electron flux
5) The System() command lead to file outputs repeated twice. fflush(NULL) was added
micromegas_2.4.1:
I) New Options
1) New functions that give the thermally averaged cross section as a function of
the temperature
2) New functions that specifies the contribution of each channel to the relic density or to the annihilation cross section
3) Function to calculate the width and branching ratios of any particle and write the result in SLHA format
4) Updated routine for the SM contribution to bsgamma in the MSSM
5) Interface to superISO
6) Possibility to set a time limit for execution of spectrum calculators in the MSSM
7) Facilities of SLHA reader were extended to work with new file formats used in superIso and HiggsBound.
II) Bugs removed
1) Compilation bug for Fedora 14 caused by extra spaces in the list of objects needed for archive file.
2) On Mac computers, sometimes Fortran and C create codes of different format (64 and 32 bits),
now we force the Fortran compiler to create codes in the same format as the C compiler.
micromegas_2.4.CPC:
I) New Options
The basic spectra tables relevant for the indirect detection module were redone.
II) Bugs removed
1) In MSSM/cycle2.c and MSSM/cycle5.c, "I" was replaced by "I_" as "I" is reserved for sqrt(-1)
2) In the routine readVarMSSM the number of free parameters was increased.
3) The PDG codes for ~o4 and ~g have been corrected in both the MSSM and the NMSSM.
Note that micrOMEGAs does not work with Cygwin now because Cygwin C compiler
does not support Complex numbers.
micromegas_2.4.O: 1. A bug that was introduced after optimisation
of version micromegas_2.4.L was removed. The value of some parameters were not properly reassigned.
This occurred only in models where some parameters depend on the scale Q.
2. For the NMSSM, the new version of NMSSMTools_2.3.1 is provided.
micromegas_2.4.Q: 1. Includes the SLHAplus library
2. Small bugs were removed
micromegas_2.4.R:
I) New Options
1. The input file MSSM/EWSB now needs parameters for first generation of sfermions. In previous versions we had assumed that the parameters of the first generation were identical to those of the second generation.
2. The new function
int slhaDecayPrint(char * name,FILE*f)
prints into a file all information about the decays of the particle 'name' in the SLHA format.
II) Bugs improved
1. The direct detection was zero when a model file was constructed in such a way that the DM particle was listed before the quarks in the particle table (prtcls1.mdl). The rate is now computed correctly in any case.
2. There was a bug with the calculation of constraints after getting a NAN in one of the constraints. Restoring the parameters used before getting NAN did not force micrOMEGAs to recalculate constraints and so led to an error. This is because micrOMEGAs does not recalculate constraints if one uses exactly the same parameters, This bug was not dangerous when performing scans over the parameter space.
3. On Mac computers it sometimes happens that one has 64-bit C and 32-bit Fortran compilers. This is fatal for C/Fortran interfaces, in particular one can not create mutual archives from C and Fortran codes. To solve this problem one has to set either the -m32 or -m64 options for both compilers. These flags are now set automatically by micrOMEGAs.
4. On Mac computers there was a conflict for define settings between micrOMEGAs and the C compiler, this lead to a problem in the compilation of MSSM/lib/gmuon.c.
5. We had a bug in the diagonalizing routines of the SLHAplus library, the rotation matrix was defined as the transpose of the standard SLHA rotation matrix. This had an impact
only for users working with a SLHAplus generated model file not for the model distributed with micrOMEGAs.
A web interface for the computation of the direct detection rates developed by Rachid Lemrani can be found here.