Micromegas v_2.2 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 (June 2008) is micromegas 2.2

Changes from micromegas 2.1
 
1. Bug removed in  nucleusRecoil routine, increase of factor 3 of  the  spin-dependent contribution. 
All WIMP-nucleon cross sections were evaluated correctly.

2. For MSSM and NMSSM we have implemented  an option to use as input  an SLHA file. For  this, comment the following lines  
#define SUGRA  
#define EWSB
#define AMSB  // for MSSM only

The main program main.c/main.F will then try to read the   SLHA file specified as a parameter 
    ./main  <SLHA file name>

3.  SuSpect(2.41) and  NMSSMTools(2.0.1) are included


4. The parameters of the routine
   cutRecoilResult(dNdE, Emin,Emax) 
are now double Emin, Emax which specify the region of recoil energy in KeV.
In the previous version these parameters were presented by integer numbers.

5. We have renormalized outgoing parameters of  
   nucleonAmplitude(LoopFactor,pA0,pA5,nA0,nA5)
routine responsible for spin-dependent interaction.  Now the  Wimp-proton cross sections are

      
SI= SCoeff*pA0[0]^2 for the spin-independent cross section and 

SD= 3*SCoeff*pA5[0]^2  for the spin-dependent one

such  that  pA5[0]  directly corresponds to the coupling $\xi_n$ defined in the
manual.


where 
     Scoeff= 4/M_PI(Nnucleon*Mwimp/(Mnucleion+Mwimp))2 



  
6. We added 2 new  routines 

extern double nucleusRecoilAux(
      double rho, double(*vfv)(double),
      int A, int Z, double J,
      double(*S00)(double),double(*S01)(double),double(*S11)(double),
      double Mwimp,
      double cs_SI_P,double cs_SI_N,  double cs_SD_P, double cs_SD_N,
      double * dNdE);
                                                                                
extern double nucleusRecoil0Aux(
      double rho, double(*vfv)(double),
      int A, int Z, double J,
      double Sp, double Sn,
      double Mwimp,
      double cs_SI_P,double cs_SI_N,  double cs_SD_P, double cs_SD_N,
      double * dNdE);

which are similar to nucleusRecoil and nucleusRecoil0 but use as input
Mwimp   - LSP mass in [GeV]
cs_SI_P - proton  spin-independent cross section in [pb]
cs_SI_N - neutron spin-independent cross section in [pb]
cs_SD_P - proton  spin-dependent   cross section in [pb]
cs_SD_N - neutron spin-dependent   cross section in [pb]

These routines depend only on the parameters specified as their arguments.  
In particular the couplings of WIMP-nucleon interactions are extracted from  the cross section. 
In order  to simulate a sign difference between  proton and 
neutron couplings one has to enter a negative cross sections. 

7. A CalcHEP online calculator is available.  From within any model directory it can  be called as 
../calc "2*2"
../calc "PI"
../calc "sqrt(0.7)*(1E32 + 100)"
Names of mathematical functions are the same as in C.

For more information see the README files and also

If you are using micromegas 2.x, 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)