Pith. Sign in

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

How this review was made: An AI referee (selfhosted/deepseek-v4.1-flash) read the paper on October 29, 2026 and answered a fixed list of questions about it. T0 is the starting level every paper gets. The levels describe how deeply Pith has checked the paper, never the authors.

  1. T0 AI review An AI read the paper and answered a fixed list of questions.
  2. T1 Audit Automatic checks that the references and files the paper relies on exist and match.
  3. T2 Compute The paper's code or calculations re-run with the same results.
  4. T3 Formal statement The claims written out in a formal math language a computer can read.
  5. T4 Proof Those claims proved, and the proof checked by computer.

How Pith checks papers →

· 3 major objections · 10 minor · Ask about this paper · The compiled record · Read the full review →
Challenge this review Comment (0) Read on arXiv ↗
Cite ▾
arXiv IDhep-ph/0405253 (version 1)
Pith IDVF2Z2EMK, Pith's permanent record of this paper and its review
CitationA. Pukhov, A. Semenov, F. Boudjema, G. Belanger. "micrOMEGAs: Version 1.3." arXiv:hep-ph/0405253 (2004). https://pith.science/paper/VF2Z2EMK
Cite this reviewBibTeX for Pith's review of this paper
Claim paper

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

The paper presents version 1.3 of micrOMEGAs, a program that answers a single sharp question: given a point in the minimal supersymmetric standard model, how much of the Universe's energy density today is carried by the lightest supersymmetric particle? The authors' aim is percent-level accuracy in exactly the regions where that number matters, namely where coannihilation with a nearly degenerate partner or annihilation through a Higgs or $Z$ resonance raises the annihilation rate enough to fall under the cosmic-microwave-background upper bound on dark matter. Their strategy is to keep the matrix elements at tree level but to feed them loop-corrected sparticle masses and mixing matrices, so that the sharply sensitive mass differences entering resonances and coannihilation are right. They also solve the abundance evolution equation numerically rather than through the freeze-out approximation, and they add higher-order corrections to the Higgs couplings to quark pairs, including the SUSY-QCD $b$-quark correction that matters at large $tan beta$. The package is meant to be usable, not just readable: it links to four public spectrum calculators, accepts the SUSY Les Houches Accord input format, and exposes cross-sections and two-body decay widths so the same model points can be taken to collider studies.

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

Editorial extensions of the paper, not claims the author makes directly.

  • 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.

Compiled by classifier-1.13-20260917 via OpenRouter, question bank abstract-p1-20260923, 2026-09-23: 1 request, 7,688 input tokens, $0.0003 reported.

What the abstract itself states, as probabilities: claims first 0.03 · states limitation 0.12 · human participants 0.01 · controlled comparison 0.01 · randomized assignment 0.01 · states quantitative result 0.04 · compares with prior results 0.07 · mentions code or data availability 0.32.

The abstract, sentence by sentence

arXiv listing abstract v1. Results and the main claim are underlined; .

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

What the classifier reports for each sentence. The number is the probability it gave the label.
#sentencerolethe authors' own workstates the main conclusionstrength
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 kind of paper the classifier reports this is: method or tool 1.00 · evidence from simulation or computation 1.00.

What was examined, and what was not

sources read
  • arXiv listing abstract v1
not examined
  • full text
  • figures and tables
  • references and citing sentences
  • relations to other papers
levels run
  • P abstract
not run
  • S sections
  • T body sentences
  • C citations
  • R cross-paper relations

Every quoted sentence on this record is the authors'. Every label is from Pith's closed vocabulary, with the probability shown beside it. Every number in the measurements table is the authors', checked by code. Record as JSON

References
How later papers used it

120 later papers cite it

Evidence: some. 110 of the citing papers share no author with it, and 94 of those carry a confident label from Pith's citation classifier.

Pith read 7 of the 135 citing sentences, choosing only those whose citation markers point at this paper alone, and wrote the paraphrases from them. The summary, the counts, the evidence level, the links and the quotes come from code; the summary counts the papers listed below. A later paper is listed under a label when Pith's citation classifier gives one of its sentences that label at probability 0.8 or more (0.95 for support), Pith's reading finds that the citing paper itself did it, and a second check confirms the paraphrase against the sentence. This section lists uses only. It does not list criticism or disagreement, so their absence here says nothing about how the paper has been received. Papers that share an author with this one are listed apart, under "By the same authors". A quote is the citing authors' own sentence, shown when their paper's licence is CC BY, CC BY-SA or CC0. JSON

Uses its method or tool 5 of 99 papers with this label
How the literature uses it

109 arXiv papers cite this in 123 sentences

Each citing sentence is taken from the citing paper's LaTeX source and typed by Pith's classifier. A label appears at probability 0.8 or more (0.95 for support), a contest at 0.7 or more; the quoted words are the citing authors'. How each label was checked · JSON

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
Compares with it 1 paper
Cites it as background 1 paper

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.