documentation

Documentation for version 2.2

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 pape r.

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

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.x 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))² 



  
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)