REVIEW
micrOMEGAs: Version 1.3
High Energy Physics - Phenomenology · Submitted 25 May 2004 · 50 references
Pith's verdict: A solid, honest tool paper: the new physics is loop-corrected masses inside exact matrix elements plus standalone cross-section and width routines, and the soft spots are accuracy claims asserted rather than measured.
T0 · AI review
Why Pith says this
Highlight any text to comment on it.
Referee report
The paper
The paper claims
micrOMEGAs 1.3 computes the relic density of the lightest supersymmetric particle using exact tree-level annihilation and coannihilation matrix elements evaluated with loop-corrected masses and mixing matrices, and reports agreement with an
The reading
What carries the argument
The load-bearing object is the relativistic thermally averaged annihilation cross-section, a sum over particle pairs weighted by degrees of freedom and by $sqrt{s},K_1(sqrt{s}/T)$ in the integrand and normalised by the equilibrium abundance built from $m_i^2 K_2(m_i/T)$. That average feeds the abundance evolution equation, which the code solves by Runge-Kutta starting from a point determined by requiring the deviation from equilibrium to be a fixed fraction of the equilibrium abundance, and finishing with an analytic integral once the equilibrium abundance has become negligible. The second piece of machinery is the treatment of the Higgs vertices: the $h q q$, $H q q$ and $A q q$ couplings are rewritten using an effective quark mass that reproduces the radiatively corrected Higgs decay widths, with QCD corrections through third order in $alpha_s$ and a SUSY-QCD correction resummed through the factor $1/(1+Delta m_b tan beta)$. Because annihilation near a Higgs resonance depends on the total width, these vertex corrections propagate directly into the relic density.
Extended reading notes
Core claim
The central claim is that a relic-density calculation built on exact tree-level matrix elements, but with loop-corrected masses and mixing matrices inserted throughout, is accurate at the percent level across the parameter regions where supersymmetric dark matter actually saturates the observed density. Two ingredients carry that claim. First, the thermal average uses relativistic formulae and the full sum over annihilation and coannihilation channels, with a Boltzmann-suppression criterion deciding which channels enter and dedicated integration treatment near poles and thresholds. Second, the evolution equation for the abundance is integrated directly by Runge-Kutta from a starting point fixed by an analytic estimate of the departure from equilibrium, then completed by an analytic tail, so no freeze-out approximation is used. The authors check the whole machinery against the independent code DarkSUSY 4.0 on ten benchmark mSUGRA points and find agreement at the 3 percent level, with a residual discrepancy at large $tan beta$ that they trace to a difference in the pseudoscalar Higgs width rather than to the annihilation calculation itself.
Load-bearing premise
The calculation assumes that inserting loop-corrected masses and mixing matrices into otherwise tree-level matrix elements captures the dominant higher-order effect, an approximation the authors themselves describe as only partial and one that they note could produce theoretical inconsistencies such as unitarity problems in some processes.
Editorial extensions
If this is right
- Relic-density predictions in the coannihilation and Higgs-resonance regions, where the measured dark matter density is actually reproduced, inherit the accuracy of the loop-corrected spectrum rather than of tree-level mass relations.
- Users can take a set of high-scale parameters, run any of four spectrum codes, and obtain a relic density computed from the same loop-corrected masses and mixings that the spectrum code produces, so the two calculations are consistent by construction.
- Because the same model file drives $b to s gamma$, $B_s to mu^+ mu^-$, $(g-2)_mu$ and $Delta rho$, a single parameter point can be confronted with several independent precision measurements at once.
- Any two-to-two MSSM process and any two-body decay width can be generated on demand and used outside the relic-density calculation, which turns the package into a general collider-simulation input generator.
- The residual disagreement with DarkSUSY at large $tan beta$ is localised to the pseudoscalar width, so it is a bookkeeping difference in a single input rather than a disagreement about the annihilation dynamics.
Reading between the lines
- If the dominant effect of loops really does enter through the mass spectrum rather than through the matrix elements, then the accuracy quoted here should degrade first in channels with large external-momentum-dependent corrections, which suggests testing the approximation channel by channel rather than only against another code built on similar assumptions.
- The comparison with an independent package is a consistency check between two codes sharing much of the same physics input; a sharper test would be a benchmark where the two codes are forced to use identical masses, widths and couplings, so that any remaining difference is attributable to the matrix elements alone.
- Extending the same mass-insertion logic to the neutralino mixing matrix, as the authors do, means the computed couplings are not derived from a single consistent Lagrangian at any loop order; users pushing into regions with large mixing, such as the stau coannihilation corridor at large $tan beta$, should treat the quoted percent-level precision as indicative rather than guaranteed.
Original abstract
We present the latest version of micromegas, a code that calculates the relic density of the lightest supersymmetric particle in the 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 some high scale (GUT). 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. Higher-order corrections to Higgs couplings to quark pairs including QCD as well as some SUSY corrections (deltaMb) are implemented. Routines calculating (g-2), bsgamma and bsmumu are also included. In particular the bsgamma routine includes an improved NLO for the SM and the charged Higgs while the SUSY large tan(beta) 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.
The compiled record
The paper's own sentences, arranged, kept beside the review and separate from it. Each label is drawn from a closed vocabulary; the number beside it is the probability the classifier gave that label. Numbers in the tables are the authors', checked by code. No sentence in this record was written by a model.
The abstract, sentence by sentence
We present the latest version of micromegas, a code that calculates the relic density of the lightest supersymmetric particle in the 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 some high scale (GUT). 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. Higher-order corrections to Higgs couplings to quark pairs including QCD as well as some SUSY corrections (deltaMb) are implemented. Routines calculating (g-2), bsgamma and bsmumu are also included. In particular the bsgamma routine includes an improved NLO for the SM and the charged Higgs while the SUSY large tan(beta) 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.
background contribution method result interpretation main claim
| # | sentence | role | the authors' own work | states the main conclusion | strength |
|---|---|---|---|---|---|
| s0 | We present the latest version of micromegas, a code that calculates the relic density of the lightest supersymmetric particle in the MSSM. | contribution 1.00 | 0.85 | 0.33 | asserted 1.00 |
| s1 | 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. | method 0.74 | 0.63 | 0.07 | asserted 1.00 |
| s2 | The cross-sections extracted from CalcHEP are calculated exactly using loop-corrected masses and mixings as specified in the SUSY Les Houches Accord. | method 0.99 | 0.59 | 0.05 | asserted 1.00 |
| s3 | Relativistic formulae for the thermal average are used and care is taken to handle poles and thresholds by adopting specific integration routines. | method 0.99 | 0.77 | 0.04 | asserted 1.00 |
| s4 | The input parameters can be either the soft SUSY parameters in a general MSSM or the parameters of a SUGRA model specified at some high scale (GUT). | method 0.89 | 0.37 | 0.04 | asserted 0.98 |
| s5 | 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. | method 0.67 | 0.55 | 0.04 | asserted 0.99 |
| s6 | Higher-order corrections to Higgs couplings to quark pairs including QCD as well as some SUSY corrections (deltaMb) are implemented. | contribution 0.94 | 0.83 | 0.04 | asserted 1.00 |
| s7 | Routines calculating (g-2), bsgamma and bsmumu are also included. | contribution 0.99 | 0.80 | 0.03 | asserted 1.00 |
| s8 | In particular the bsgamma routine includes an improved NLO for the SM and the charged Higgs while the SUSY large tan(beta) effects beyond leading-order are included. | contribution 0.93 | 0.76 | 0.04 | asserted 1.00 |
| s9 | 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. | contribution 1.00 | 0.87 | 0.12 | asserted 1.00 |
What was examined, and what was not
- arXiv listing abstract v1
- full text
- figures and tables
- references and citing sentences
- relations to other papers
- P abstract
- S sections
- T body sentences
- C citations
- R cross-paper relations
References
How later papers used it
120 later papers cite it
5 later papers by other authors are listed below: 5 used its method or tool.
Uses its method or tool 5 of 99 papers with this label
-
Fed its resulting low energy spectrum into the Micromegas program, version 2.0.1, assuming the low energy effective theory is the MSSM with R-parity conservation.
Dark Matter in Gauge Mediation from Emergent Supersymmetry (2007) · 1.00 · paraphrased: the licence does not allow quoting
-
Calculated the relic abundance of the neutralino dark matter using micrOMEGAs 2.4 with the output of SOFTSUSY in the SLHA format.
Discrimination of Supersymmetric Grand Unified Models in Gaugino Mediation (2010) · 1.00 · paraphrased: the licence does not allow quoting
-
Used the micrOMEGAs package version 2.4.5 to calculate the neutralino dark matter relic abundance and other phenomenological constraints with the output of SOFTSUSY.
Positively deflected anomaly mediation in the light of the Higgs boson discovery (2012) · 1.00 · paraphrased: the licence does not allow quoting
-
Linked the MSSM spectrum via the SUSY Les Houches Accord to micrOMEGAs1.3.6 for points not ruled out, which then calculates relic density, branching ratios and the muon anomalous magnetic moment.
Natural Priors, CMSSM Fits and LHC Weather Forecasts (2007) · 1.00 · paraphrased: the licence does not allow quoting
-
Computed some of the associated relic densities with the micrOMEGAs code.
Antimatter cosmic rays from dark matter annihilation: First results from an N-body experiment (2008) · 1.00 · paraphrased: the licence does not allow quoting
How the literature uses it
109 arXiv papers cite this in 123 sentences
Papers citing it · 23 sentences without a confident label
Uses its method or tool 90 papers
-
We also use [ref] to obtain the relic density of the LSP, $ΩLSP$, and its spin-independent cross section to nucleons, $σSI$.
Wino-Higgsino dark matter in the MSSM from the $g-2$ anomaly · Introduction · 1.00
-
This is then interfaced with FeynHiggs 2.14.2 [ref] to calculate the branching ratios of Higgs boson and $ρ$ parameter, Superiso v4.0 [ref] for estimating muon anomalous magnetic moment and B-physics branching ratios and micromegas 5.0.4 [ref] to calculate the relic density of the dark matter which in our case turn out to be lightest neutralino.
An Information Theoretic Exploration of Constrained MSSM · CMSSM and Information Theory · 1.00
-
The relic density and direct detection cross section were calculated numerically with [ref] , using [ref] model files automatically generated with .
Dark Matter in a Constrained $E_6$ Inspired SUSY Model · Scan Procedure · 1.00
Uses its data or values 1 paper
-
The SM value for the branching $$ calculated in [ref] was shifted in order to match the NNLO result of [ref] .
Constraining the MSSM with universal gaugino masses and implication for searches at the LHC · Results · 0.88
Compares with it 1 paper
-
In this way, possible NMSSM DM under-abundance is consistently included in the derivation of the direct detection signal for the NMSSM DM candidate.[We compared our numbers for direct detection with MicrOMEGAs [ref] and found agreement.]
NMSSMScanner: Efficient Scans in the NMSSM Parameter Space Proof of Concept · Experimental Observables · 0.96
Cites it as background 1 paper
-
Several public computer codes that determine the MSSM mass spectrum and other SUSY observables (e.g. SuSpect [ref] , SPheno [ref] , micrOMEGAs [ref] , FeynHiggs [ref] , NMHDECAY [ref] , CPsuperH [ref] and SuperIso [ref] ) contain calculations of $$ in various approximations.
SusyBSG: a fortran code for BR[B -> Xs gamma] in the MSSM with Minimal Flavor Violation · Introduction · 0.86
Discussion
Disagree with the review, or know something it missed? Leave a comment, highlight any passage above to comment on it, or challenge this review formally.