micrOMEGAs 2.0.7: a program to calculate the relic density of dark matter in a generic model

详细信息    查看全文

• 作者：G. Bé ; langer ; F. Boudjema ; A. Pukhov ; A. Semenov
• 关键词：Perturbative quantum field theory ; Mellin– ; Barnes integrals ; Feynman integrals
• 刊名：Computer Physics Communications
• 出版年：2007
• 出版时间：1 December 2007
• 年：2007
• 卷：177
• 期：11
• 页码：894-895
• 全文大小：80 K

文摘

micrOMEGAs2.0.7 is a code which calculates the relic density of a stable massive particle in an arbitrary model. The underlying assumption is that there is a conservation law like R-parity in supersymmetry which guarantees the stability of the lightest odd particle. The new physics model must be incorporated in the notation of CalcHEP, a package for the automatic generation of squared matrix elements. Once this is done, all annihilation and coannihilation channels are included automatically in any model. Cross-sections at v=0, relevant for indirect detection of dark matter, are also computed automatically. The package includes three sample models: the minimal supersymmetric standard model (MSSM), the MSSM with complex phases and the NMSSM. Extension to other models, including non supersymmetric models, is described.

Program summary

Title of program:micrOMEGAs2.0.7

Catalogue identifier:ADQR_v2_1

Program summary URL: http://cpc.cs.qub.ac.uk/summaries/ADQR_v2_1.html

Program obtainable from: CPC Program Library, Queen's University of Belfast, N. Ireland

Licensing provisions:Standard CPC licence, http://cpc.cs.qub.ac.uk/licence/licence.html

No. of lines in distributed program, including test data, etc.:216 529

No. of bytes in distributed program, including test data, etc.:1 848 816

Distribution format:tar.gz

Programming language used:C and Fortran

Computer:PC, Alpha, Mac, Sun

Operating system:UNIX (Linux, OSF1, SunOS, Darwin, Cygwin)

RAM:17 MB depending on the number of processes required

Classification:1.9, 11.6

Catalogue identifier of previous version:ADQR_v2_0

Journal version of previous version:Comput. Phys. Comm. 176 (2007) 367

Does the new version supersede the previous version?:Yes

Nature of problem:Calculation of the relic density of the lightest stable particle in a generic new model of particle physics.

Solution method:In numerically solving the evolution equation for the density of dark matter, relativistic formulae for the thermal average are used. All tree-level processes for annihilation and coannihilation of new particles in the model are included. The cross-sections for all processes are calculated exactly with CalcHEP after definition of a model file. Higher-order QCD corrections to Higgs couplings to quark pairs are included.

Reasons for new version:The main changes in this new version consist, on the one hand, in improvements of the user interface and treatment of error codes when using spectrum calculators in the MSSM and, on the other hand, on a completely revised code for the calculation of the relic density in the NMSSM based on the code NMSSMTools1.0.2 for the computation of the spectrum.

Summary of revisions:

• The version of CalcHEP was updated to CalcHEP 2.4.

• The procedure for shared library generation has been improved. Now the libraries are recalculated each time the model is modified.

• The default value for the top quark mass has been set to 171.4 GeV.

Changes specific to the MSSM model.

• The deltaMb correction is now included in the B,t,H-vertex and is always included for other Higgs vertices.

• In case of a fatal error in an RGE program, micrOMEGAs now continues operation while issuing a warning that the given point is not valid. This is important when running scans over parameter space. However this means that the standard ˆC command that could be used to cancel a job now only cancels the RGE program. To cancel a job, use “kill -9 -N” where N is the micrOMEGAs process id, all child processes launched by micrOMEGAs will be killed at once.

• Following the last SLHA2 release, we use key=26 item of EXTPAR block for the pole mass of the CP-odd Higgs so that micrOMEGAs can now use SoftSUSY for spectrum calculation with EWSB input. The Isajet interface was corrected too, so the user has to recompile the isajet_slha executable. For SuSpect we still support an old “wrong” interface where key=24 is used for the mass of the CP-odd Higgs.

• In the non-universal SUGRA model, we set the value of M₀ (M_1/2,A₀) to the value of the largest subset of equal parameters among scalar masses (gaugino masses, trilinear couplings). In the previous version these parameters were set arbitrarily to be equal to MH2, MG2 and At respectively. The spectrum calculators need an input value for 33736723e960922a1078df30"" title=""Click to view the MathML source"">M₀,M_1/2 and A₀ for initialisation purposes.

• We have removed bugs in micrOMEGAs-Isajet interface in case of non-universal SUGRA.

• $(FFLAGS) is added to compilation instruction of suspect.exe. It was omitted in version 2.0.

• The treatment of errors in reading of the LesHouches accord file is improved. Now, if the SPINFO block is absent in the SLHA output it is considered as a fatal error.

• Instructions for calculation of Δ_ρ, (g−2)_μ, Br(b→sγ) and Br(B_s→μ⁺μ⁻) constraints are included in EWSB sample main programs omg.c/omg.cpp/omg.F.

• We have corrected the name of the library for neutralino–neutralino annihilation in our sample files MSSM/cs br.*.

Changes specific to the NMSSM model.

• The NMSSM has been completely revised. Now it is based on NMSSMTools_1.0.2.

• The deltaMb corrections in the NMSSM are included in the Higgs potential.

CP violation model.

• We have included in our package the MSSM with CP violation. Our implementation was described in Phys. Rev. D 73 (2006) 115007. It is based on the CPSUPERH package published in Comput. Phys. Comm. 156 (2004) 283.

Unusual features:Depending on the parameters of the model, the program generates additional new code, compiles it and loads it dynamically.

Running time:0.2 seconds