micrOMEGAs is a code to compute dark matter observables in generic extensions of the standard model. This version of micrOMEGAs includes a generalization of the Boltzmann equations to take into ...account the possibility of two dark matter candidates. The modification of the relic density calculation to include interactions between the two dark matter sectors as well as semi-annihilation is presented. Both dark matter signals in direct and indirect detection are computed as well. An extension of the standard model with two scalar doublets and a singlet is used as an example.
Program title: MicrOMEGAs4.1
Catalogue identifier: ADQR_v4_0
Program summary URL:http://cpc.cs.qub.ac.uk/summaries/ADQR_v4_0.html
Program obtainable from: CPC Program Library, Queen’s University, 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.: 738425
No. of bytes in distributed program, including test data, etc.: 9807620
Distribution format: tar.gz
Programming language: C and Fortran.
Computer: PC, Mac.
Operating system: UNIX (Linux, Darwin).
RAM: 50MB depending on the number of processes required.
Classification: 1.9, 11.6.
Catalogue identifier of previous version: ADQR_v3.0
Journal reference of previous version: Comput. Phys. Comm. 185 (2014) 960
External routines: CalcHEP, SuSpect, NMSSMTools, CPSuperH, LoopTools, HiggsBounds
Does the new version supersede the previous version?: Yes
Nature of problem: Calculation of the relic density and direct and indirect detection rates of the lightest stable particle in particle physics models with at most two stable dark matter candidates.
Solution method: In the case where the two dark matter particles have very different masses, we find that the equations for the evolution of the density of dark matter behave as stiff equations. To solve these we use the backward scheme and the Rosenbrock algorithm. The standard solution based on the Runge–Kutta method is still used for models with only one dark matter candidate.
Reasons for new version: There are many experiments that are currently searching for the remnants of dark matter annihilation and the relic density is determined precisely from cosmological measurements. In this version we generalize the Boltzmann equations to take into account the possibility of two dark matter candidates. Thus, in solving for the relic density we include interactions between the two dark matter sectors as well as semi-annihilation. The dark matter signals in direct and indirect detection are computed as well.
Summary of revisions:•Generalization of the Boltzmann equations to include two dark matter candidates, their interactions and semi-annihilations, the relative density of the two dark matter components is taken into account when computing direct/indirect detection rates.•Upgrade of the numerical method for solving the Boltzmann equations.•Include sample extensions of the standard model with extra doublet and singlets which contain two stable neutral particles.Unusual features: Depending on the parameters of the model, the program generates additional new code, compiles it and loads it dynamically.
Running time: 4 sec
micrOMEGAs is a code to compute dark matter observables in generic extensions of the standard model. This new version of micrOMEGAs is a major update which includes a generalization of the ...Boltzmann equations to accommodate models with asymmetric dark matter or with semi-annihilation and a first approach to a generalization of the thermodynamics of the Universe in the relic density computation. Furthermore a switch to include virtual vector bosons in the final states in the annihilation cross sections or relic density computations is added. Effective operators to describe loop-induced couplings of Higgses to two-photons or two-gluons are introduced and reduced couplings of the Higgs are provided allowing for a direct comparison with recent LHC results. A module that computes the signature of DM captured in celestial bodies in neutrino telescopes is also provided. Moreover the direct detection module has been improved as concerns the implementation of the strange “content” of the nucleon. New extensions of the standard model are included in the distribution.
Title of program: micrOMEGAs3.
Program obtainable from:http://lapth.cnrs.fr/micromegas
Computers for which the program is designed and others on which it has been tested: PC, Mac
Operating systems under which the program has been tested : UNIX (Linux, Darwin)
Programming language used: C and Fortran
Memory required to execute with typical data: 50 MB depending on the number of processes required.
No. of processors used: 1
Has the code been vectorized or parallelized: no
No. of bytes in distributed program, including test data: 70736 kB
External routines/libraries used: no
CPC Program Library subprograms used: CalcHEP, SuSpect, NMSSMTools, CPSuperH, LoopTools, HiggsBounds
Catalogue identifier of previous version: ADQR_v1_3
Journal reference of previous version: Comput. Phys. Comm. 182 (2011) 842
Does the new version supersede the previous version: yes
Nature of physical problem: Calculation of the relic density and direct and indirect detection rates of the lightest stable particle in a generic new model of particle physics.
Method of solution: 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 as well as some 3-body final states. The cross-sections for all processes are calculated exactly with CalcHEP after definition of a model file. The propagation of the charged cosmic rays is solved within a semi-analytical two-zone model.
Reasons for the new version: There are many experiments that are currently searching for the remnants of dark matter annihilation and the relic density is determined precisely from cosmological measurements. In this version we add the computation of dark matter signals in neutrino telescopes, we generalize the Boltzmann equations so as to take into account a larger class of dark matter models and improve the precision in the prediction of the relic density for DM masses that are below the W mass. We compute the signal strength for Higgs production in different channels to compare with the results of the LHC.
Summary of revisions:•Generalization of the Boltzmann equations to include asymmetric dark matter and semi-annihilations: the DM asymmetry is taken into account when computing direct/indirect detection rates.•Incorporating loop-induced decays of Higgs particles to two-photons and two-gluons, and computing the signal strength for Higgs production in various channels that can be compared to results from LHC searches.•New module for neutrino signature from DM capture in the Sun and the Earth•Annihilation cross sections for some selected 3-body processes in addition to the 2-body tree-level processes. The 3-body option can be included in the computation of the relic density and/or for annihilation of dark matter in the galaxy.•Possibility of using different tables for the effective degrees of freedom in the early Universe•Annihilation cross sections for the loop induced processes γγ and γZ0 in the NMSSM and the CPVMSSM•New function for incorporating DM clumps•New function to define the strange quark content of the nucleon•The LanHEP source code for new models is included•New models with scalar DM are included (Inert doublet model and model with Z3 symmetry)•New implementation of the NMSSM which uses the Higgs self-couplings and the particle spectrum calculated in NMSSMTools_4.1•New versions of spectrum generators used in the MSSM (Suspect_2.4.1) and in the CPVMSSM (CPsuperH2.3)•Extended routines for flavor physics in the MSSM•New facilities to compute DM observables independently of the model•Update in interface tools to read files produced by other codes, this allows easy interface to other codesTypical running time: 4 s
Unusual features of the program: Depending on the parameters of the model, the program generates additional new code, compiles it and loads it dynamically.
We present a new module of
micrOMEGAs devoted to the computation of indirect signals from dark matter annihilation in any new model with a stable weakly interacting particle. The code provides the ...mass spectrum, cross-sections, relic density and exotic fluxes of gamma rays, positrons and antiprotons. The propagation of charged particles in the Galactic halo is handled with a new module that allows to easily modify the propagation parameters.
Program title: micrOMEGAs2.4
Catalogue identifier: ADQR_v2_3
Program summary URL:
http://cpc.cs.qub.ac.uk/summaries/ADQR_v2_3.html
Program obtainable from: CPC Program Library, Queen's University, 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.: 401 126
No. of bytes in distributed program, including test data, etc.: 6 583 596
Distribution format: tar.gz
Programming language: C and Fortran
Computer: PC, Alpha, Mac, Sun
Operating system: UNIX (Linux, OSF1, SunOS, Darwin, Cygwin)
RAM: 50 MB depending on the number of processes required
Classification: 1.9, 11.6
Catalogue identifier of previous version: ADQR_v2_3
Journal reference of previous version: Comput. Phys. Comm. 180 (2009) 747
Does the new version supersede the previous version?: Yes
Nature of problem: Calculation of the relic density and detection rates 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 formulas 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. The propagation of the charged cosmic rays is solved within a semi-analytical two-zone model.
Reasons for new version: There are many experiments that are currently searching for the remnants of dark matter annihilation. In this version we perform the computation of indirect signals from dark matter annihilation in any new model with a stable weakly interacting particle. We include the propagation of charged particles in the Galactic halo.
Summary of revisions:
•
Annihilation cross-sections for all 2-body tree-level processes and for radiative emission of a photon for all models.
•
Annihilation cross-sections into polarised gauge bosons.
•
Annihilation cross-sections for the loop induced processes
γγ and
γ
Z
0
in the MSSM.
•
Modelling of the DM halo with a general parameterization and with the possibility of including DM clumps.
•
Computation of the propagation of charged particles through the Galaxy, including the possibility of modifying the propagation parameters.
•
Effect of solar modulation on the charged particle spectrum.
•
Model independent predictions of the indirect detection signals.
Unusual features: Depending on the parameters of the model, the program generates additional new code, compiles it and loads it dynamically.
Running time: 3 sec
micrOMEGAs: Version 1.3 Bélanger, G.; Boudjema, F.; Pukhov, A. ...
Computer physics communications,
04/2006, Letnik:
174, Številka:
7
Journal Article
Recenzirano
We present the latest version of
micrOMEGAs, a code that calculates the relic density of the lightest supersymmetric particle (LSP) in the minimal supersymmetric standard model (MSSM). All tree-level ...processes for the annihilation of the LSP are included as well as all possible coannihilation processes with neutralinos, charginos, sleptons, squarks and gluinos. The cross-sections extracted from
CalcHEP are calculated exactly using loop-corrected masses and mixings as specified in the
SUSY Les Houches Accord. Relativistic formulae for the thermal average are used and care is taken to handle poles and thresholds by adopting specific integration routines. The input parameters can be either the soft SUSY parameters in a general MSSM or the parameters of a SUGRA model specified at the GUT scale. In the latter case, a link with
Suspect,
SOFTSUSY,
Spheno and
Isajet allows one to calculate the supersymmetric spectrum, Higgs masses, as well as mixing matrices. Higher-order corrections to Higgs couplings to quark pairs including QCD as well as some SUSY corrections (
Δ
m
b
) are implemented. Routines calculating
(
g
−
2
)
μ
,
b
→
s
γ
and
B
s
→
μ
+
μ
−
are also included. In particular the
b
→
s
γ
routine includes an improved NLO for the SM and the charged Higgs while the SUSY large
tan
β
effects beyond leading-order are included. This new version also provides cross-sections for any
2
→
2
process as well as partial decay widths for two-body final states in the MSSM allowing for easy simulation at colliders.
Program title:micrOMEGAs1.3
Catalogue identifier:ADQR_v1_3
Program summary URL:
http://cpc.cs.qub.ac.uk/summaries/ADQR_v1_3
Program obtainable from: CPC Program Library, Queen's University of Belfast, N. Ireland
Licensing provisions:none
Computer:PC, Alpha, Silicon graphics, Sun
Programming language:C and Fortran
Operating system:UNIX (Linux, OSF1, IRIX64, SunOS)
RAM:20 MB depending on the number of processes required
No of lines in distributed program, including test data, etc.:78 314
No. of bytes in distributed program, including test data, etc.:703 112
Distribution format:tar.gz
Number of processors used:1
External routines/libraries:Library of Fortran functions, for example, -lg2c (platform dependent)
Catalogue identifier of previous version:ADQR
Journal reference of previous version:Comput. Phys. Comm. 149 (2002) 103
Does the new version supersede the previous version?:yes
Nature of problem:Calculation of the relic density of the lightest supersymmetric particle in the MSSM.
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 SUSY particles are included. The cross-sections for all processes are calculated exactly with CalcHEP. Higher-order corrections to Higgs masses and Higgs couplings to quark pairs including QCD as well as some SUSY corrections are implemented. The input parameters can be either the soft SUSY parameters in a general MSSM or the parameters of a SUGRA model specified at the GUT scale. In the latter case, a link with
Suspect,
SOFTSUSY,
Spheno and
Isajet allows to calculate the supersymmetric spectrum, Higgs masses, as well as mixing matrices.
Reasons for the new version:This new version contains a more accurate calculation of the relic density of dark matter as well as many new features both for interface with codes that calculate the supersymmetric spectrum as well as for computation of cross-sections and decays relevant for collider physics.
Summary of revisions:
•
Interface with the main codes to calculate the supersymmetric spectrum:
Suspect,
Isajet,
Spheno and
SOFTSUSY in models defined at some high scale.
•
Includes loop corrected sparticle masses and mixing matrices.
•
Includes loop-corrected Higgs masses and widths. QCD corrections to the Higgs couplings to fermion pairs are included as well as, via an effective Lagrangian, the
Δ
m
b
correction relevant at large
tan
β
.
•
Provides exact numerical solution of the Boltzmann equation by Runge–Kutta.
•
Outputs the relative contribution of each channel to
1
/
Ω
.
•
Computes cross-sections for any
2
→
2
process at the parton level.
•
Calculates decay widths for all particles at tree-level including all
1
→
2
decay modes.
•
Calculates constraints on MSSM:
B
s
→
μ
+
μ
−
and NLO corrections to
b
→
s
γ
.
Unusual features:Depending on the parameters of the model, the program generates additional new code for matrix elements, compiles it and loads it dynamically.
Running time:0.2 sec.
We present a full one-loop renormalization of the Higgs sector of the next-to-minimal-supersymmetric-Standard-Model (NMSSM) and its implementation within sloops, a code for the automated computations ...of one-loop processes in theories beyond the Standard Model. The present work is the sequel to the study we performed on the renormalization of the sectors of the NMSSM comprising neutralinos, charginos, and sfermions, thereby completing the full one-loop renormalization of the NMSSM. We have investigated several renormalization schemes based on alternative choices (on-shell or DR¯) of the physical parameters. Special attention is paid to the issue of the mixing between physical fields. To weigh the impact of the different renormalization schemes, the partial widths for the decays of the Higgs bosons into supersymmetric particles are computed at one loop. In many decays large differences between the schemes are found. We discuss the origin of these differences. In particular, we study two contrasting scenarios. The first model is MSSM-like with a small value for the mixing between the doublet and singlet Higgs superfields while the second model has a moderate value for this mixing. We critically discuss the issue of the reconstruction of the underlying parameters and their counterterms in the case of a theory with a large number of parameters, such as the NMSSM, from a set of physical parameters. In the present study this set corresponds to the minimum set of masses for the implementation of the on-shell schemes.
micrOMEGAs is a numerical code to compute dark matter (DM) observables in generic extensions of the Standard Model (SM) of particle physics. We present a new version of micrOMEGAs that includes a ...generalization of the Boltzmann equations governing the DM cosmic abundance evolution which can be solved to compute the relic density of N-component DM. The direct and indirect detection rates in such scenarios take into account the relative contribution of each component such that constraints on the combined signal of all DM components can be imposed. The co-scattering mechanism for DM production is also included, whereas the routines used to compute the relic density of feebly interacting particles have been improved in order to take into account the effect of thermal masses of t-channel particles. Finally, the tables for the DM self-annihilation - induced photon spectra have been extended down to DM masses of 110 MeV, and they now include annihilation channels into light mesons.
Program title: micrOMEGAs6.0
CPC Library link to program files:https://doi.org/10.17632/4ck6jf5vxf.3
Licensing provisions: GNU General Public License 3
Programming language: C and Fortran
Journal reference of previous version: Comput. Phys. Comm. 231 (2018) 173.
Does the new version supersede the previous version?: Yes
Reasons for the new version: Previous versions of micrOMEGAs worked within the assumption that DM is composed of one or two components. The new version allows for more components which can be either weekly or feebly interacting. The possibility of co-scattering is also implemented.
Summary of revisions: This version includes new routines to compute the abundance of multi-component DM that contains either weakly or feebly interacting dark matter particles in generic extensions of the SM of particle physics. The co-scattering mechanism for DM production is also included. The routines to compute the relic density of feebly interacting particles through the freeze-in mechanism have been improved in order to take into account the effect of thermal masses of t-channel particles. The tables for the photon spectra resulting from pair annihilation have been extended down to dark matter masses of 110 MeV and they now include annihilation channels into light mesons.
Nature of problem: DM candidates that satisfy cosmological constraints cover a wide range of masses and interaction strengths. Moreover, the dark sector could contain several stable neutral particles that can all contribute to DM. We provide the first public code to perform a precise computation of the relic density for generic extensions of the SM with more than two component dark matter.
Solution method: We solve N Boltzmann equations treating both the cases where the DM components are in thermal equilibrium with the thermal bath in the early Universe, as well as the case where the DM is too feebly interacting to reach equilibrium. We also include decay terms in the Boltzmann equations. All the signals for DM direct and indirect detection take into account the contribution of each component to the total relic density.
We have completed the one-loop renormalisation of the Next-to-Minimal Supersymmetric Standard Model (NMSSM) allowing for and comparing between different renormalisation schemes. A special attention ...is paid to on-shell schemes. We study a variety of these schemes based on alternative choices of the physical input parameters. In this paper we present our approach to the renormalization of the NMSSM and report on our results for the neutralino-chargino and sfermion sectors. We will borrow some results from our study of the Higgs sector whose full discussion is left for a separate publication. We have implemented the setup for all the sectors of the NMSSM within SloopS, a code for the automatic computation of one-loop corrections initially developed for the standard model and the MSSM. Among the many applications that allows the code, we present here the one-loop corrections to neutralino masses and to partial widths of neutralinos and charginos into final states with one gauge boson. One-loop electroweak and QCD corrections to the partial widths of third generation sfermions into a fermion and a chargino or a neutralino are also computed.
The Supersymmetry Les Houches Accord (SLHA) provides a universal set of conventions for conveying spectral and decay information for supersymmetry analysis problems in high energy physics. Here, we ...propose extensions of the conventions of the first SLHA to include various generalisations: the minimal supersymmetric standard model with violation of CP, R-parity, and flavour, as well as the simplest next-to-minimal model.
The inferred value of the relic density from cosmological observations has reached a precision that is akin to that of the LEP precision measurements. This level of precision calls for the evaluation ...of the annihilation cross sections of dark matter that goes beyond tree-level calculations as currently implemented in all codes for the computation of the relic density. In supersymmetry radiative corrections are known to be large and thus must be implemented. Full one-loop radiative corrections for many annihilation processes have been performed. It is important to investigate whether the bulk of these corrections can be parametrized through an improved Born approximation that can be implemented as a selection of form factors to a tree-level code. This paper is a second in a series that addresses this issue. After having provided these form factors for the annihilation of the neutralinos into fermions, which cover the case of a bino-like lightest supersymmetric particle (LSP), we turn our attention here to a higgsino-like dark matter candidate through its annihilation into ZZ. We also investigate the cases of a mixed LSP. In all cases we compare the performance of the form factor approach with the result of a full one-loop correction. We also study the issue of the renormalization scheme dependence. An illustration of the phenomenon of nondecoupling of the heavy sfermions that takes place for the annihilation of the lightest neutralino into ZZ is also presented.
We describe the main building blocks of a generic automated package for the calculation of Feynman diagrams. These blocks include the generation and creation of a model file, the graph generation, ...the symbolic calculation at an intermediate level of the Dirac and tensor algebra, implementation of the loop integrals, the generation of the matrix elements or helicity amplitudes, methods for the phase space integrations and eventually the event generation. The report focuses on the fully automated systems for the calculation of physical processes based on the experience in developing
GRACE-loop which is a general purpose code applicable to one-loop corrections in the Standard Model. As such, a detailed description of the renormalisation procedure in the Standard Model is given emphasizing the central role played by the non-linear gauge fixing conditions for the construction of such automated codes. These new gauge-fixing conditions are used as a very efficient means to check the results of large scale automated computations in the Standard Model. Their need is better appreciated when it comes to devising efficient and powerful algorithms for the reduction of the tensorial structures of the loop integrals and the reduction of the
N
>
4
point-function to lower rank integrals. A new technique for these reduction algorithms is described. Explicit formulae for all two-point functions in a generalised non-linear gauge are given, together with the complete set of counterterms. We also show how infrared divergences are dealt with in the system. We give a comprehensive presentation of some systematic test-runs which have been performed at the one-loop level for a wide variety of two-to-two processes to show the validity of the gauge check. These cover fermion–fermion scattering, gauge boson scattering into fermions, gauge bosons and Higgs bosons scattering processes. Comparisons with existing results on some one-loop computation in the Standard Model show excellent agreement. These include
e
+
e
-
→
t
t
¯
,
W
+
W
-
,
ZH
;
γ
γ
→
t
t
¯
,
W
+
W
-
;
e
γ
→
eZ
,
ν
W
and
W
+
W
-
→
W
+
W
-
. We also briefly recount some recent development concerning the calculation of one-loop corrections to 3 body final states cross sections in
e
+
e
-
with the help of an automated system.