Micromegas v_2.1 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 drak matter candidate at zero velocity, relevant for indirect detection of dark matter, are also computed automatically.
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).


1. The MSSM model is as described in micromegas_1.3
  The package also includes routines for g-2, b to s gamma, b_s to mu mu, experimental direct LEP constraints,..

2. The NMSSM model now uses the Susy Les Houches Accord 2 , for a description and physics studies see this paper.
      The NMSSM file uses NMHDECAY and includes LEP constraints.

3. The CPV   model  uses CPsuperH. A   routine for the electron and thallium edm (at two-loop) is included. For a description and physics studies see this paper.

4. The RHN model is a generic model with a Dirac right-handed neutrino inspired by extra-dimensional models that contains also extra gauge bosons and fermions. For a description and physics studies see this paper.

5. For a physics study using the  little Higgs model see this paper.

6. Implementation of new models first requires to write a new CalcHEP model file.

7. If you are using this program please refer to the standard micromegas manual   and the 2.1 update document.

Present version (March 2008) is micromegas 2.1

For more information see the README files and also

If you are using micromegas 2.1, please cite
G. Bélanger, F. Boudjema, A. Pukhov, A. Semenov, arXiv:0803.2360 [hep-ph]
and
G. Bélanger, F. Boudjema, A. Pukhov, A. Semenov, Comput.Phys.Commun.176:367-382,2007 hep-ph/0607059

Moreover if you are using the MSSM module, please also cite
G. Bélanger, F. Boudjema, A. Pukhov and A. Semenov, Comput. Phys. Commun. 174 (2006) 577; hep-ph/0405253.
and
G. Bélanger, F. Boudjema, A. Pukhov and A. Semenov, Comput. Phys. Commun. 149 (2002) 103; hep-ph/0112278.


If you are using the CPVMSSM please also cite
G. Bélanger, F. Boudjema,  S. Kraml, A. Pukhov, A. Semenov, Phys. Rev. D73:115007 (2006); hep-ph/0604150

If you are using the NMSSM please also cite
G. Bélanger, F. Boudjema,  C. hugonie, A. Pukhov, A. Semenov, JCAP0509:001(2005); hep-ph/0505142


If you are using the RHNM please also cite
G. Belanger, A.Pukhov, G.Servant, JCAP 0801:009 (2008); arXiv:0706.0526 [hep-ph].

If you are using the LHM please also cite
 A. Belyaev, C.R.Chen, K. Tobe, C.P. Yuan,  Phys.Rev. D74:115020 (2006)