Pith. Sign in

REVIEW

A Markov Chain Monte Carlo Analysis of the CMSSM

High Energy Physics - Phenomenology · Submitted 2 Feb 2006 · 104 references

Pith's verdict: A careful, genuinely useful Bayesian MCMC scan of the CMSSM whose headline robustness claim holds; the m0 and σ_SI^p numbers are softer than the abstract's precision suggests, and the authors know it.

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 →

· 2 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/0602028 (version 3)
Pith IDT2GPYWOZ, Pith's permanent record of this paper and its review
CitationLeszek Roszkowski, Roberto Ruiz de Austri, Roberto Trotta. "A Markov Chain Monte Carlo Analysis of the CMSSM." arXiv:hep-ph/0602028 (2006). https://pith.science/paper/T2GPYWOZ
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

A Markov-chain fit of the constrained MSSM turns collider, cosmological and rare-decay data into probability intervals for the theory's five parameters, and finds that most of the favored superpartner mass range sits within LHC reach.

The reading

The paper asks where in the Constrained Minimal Supersymmetric Standard Model (CMSSM) — the MSSM with common scalar, gaugino and trilinear masses fixed at the unification scale — the existing data actually point. Instead of scanning a fixed grid and discarding points that fail experimental bounds, the authors run a Markov chain Monte Carlo over eight parameters, four of them supersymmetric ($m_0$, $m_{1/2}$, $A_0$, $tan beta$) and four Standard Model nuisance parameters ($M_t$, $m_b(m_b)$, $alpha_em$, $alpha_s$), weighting each sample by a likelihood that folds theoretical and experimental errors together. The output is a posterior probability density rather than a best-fit point, and its 68% intervals are $0.52$ TeV $< m_{1/2} < 1.26$ TeV, $m_0 < 2.10$ TeV, $-0.34$ TeV $< A_0 < 2.41$ TeV and $38.5 < tan beta < 54.6$. Those numbers matter because they translate directly into reach: the paper puts 78% probability on a gluino below 2.7 TeV, 85% on a right-handed squark below 2.5 TeV and 65% on the lightest chargino below 0.8 TeV, and it locates the spin-independent WIMP-proton cross section mostly between $10^{-10}$ and $10^{-8}$ pb. The paper's own headline conclusion is a statement about robustness: $m_{1/2}$, $A_0$ and $tan beta$ barely move when the prior ranges are doubled, while $m_0$ does, because a broad focus-point region at multi-TeV scalar masses reproduces the measured dark-matter abundance just as well.

What carries the argument

The load-bearing device is the likelihood construction. Experimental upper and lower limits are not applied as step functions; instead each 95% confidence bound is convolved with a Gaussian of width $tau$ representing the theoretical uncertainty of the mapping from parameters to observables, which smears the edge of the allowed region over the scale $tau$ (their Eq. 3.5). Positive measurements are treated as Gaussians with experimental and theoretical errors added in quadrature. Around this likelihood sits a Metropolis sampler that alternates updates of the SUSY block and the nuisance block, proposes steps along directions drawn from the eigenvectors of an estimated posterior covariance matrix, and is monitored with the Gelman-Rubin statistic until $R-1 < 0.05$ across 24 to 32 chains. The spectra come from SOFTSUSY, the relic density from an exact coannihilation calculation, and the rare decays from next-to-leading-order and one-loop treatments with large-$tan beta$ corrections.

Extended reading notes

Core claim

The central claim is that a Bayesian global fit of the CMSSM, with all relevant uncertainties carried explicitly, yields well-defined and mostly prior-independent probability intervals, and that these intervals can be converted into concrete expectations for the LHC and for WIMP direct detection. Concretely, the authors report 68% ranges of $0.52$ TeV $< m_{1/2} < 1.26$ TeV, $m_0 < 2.10$ TeV, $-0.34$ TeV $< A_0 < 2.41$ TeV and $38.5 < tan beta < 54.6$, with $m_{1/2}$, $A_0$ and $tan beta$ stable under a change of prior range while $m_0$ is not. Two secondary findings carry weight. First, the narrow 'WMAP strips' familiar from fixed-grid scans widen into a large region once every parameter is allowed to vary and the others are marginalized, so the relic-density constraint is less restrictive than it had appeared. Second, high posterior probability and good quality of fit are not the same thing: the posterior prefers large $tan beta$, but acceptable fits exist down to $tan beta$ around 15, a discrepancy the authors say only better data can settle. The paper also maps correlations among observables, most usefully between $BR(B_s to mu^+ mu^-)$ and both $BR(bar{B} to X_s gamma)$ and $sigma_p^{SI}$, which it proposes as a handle for distinguishing the CMSSM from other supersymmetric models once a signal is found.

Load-bearing premise

The quoted intervals are only as wide as the hand-chosen theoretical error bars fed into the likelihood, above all a single 10% error on the dark-matter relic abundance applied everywhere and a flat 5% error on every superpartner mass, even though the authors say the relic-density error in the focus-point region is almost certainly much larger.

Editorial extensions

If this is right

  • The LHC should cover a large share of the favored region: the paper quotes 78% probability for $m_{gluino} < 2.7$ TeV, 85% for $m_{squark_R} < 2.5$ TeV and 65% for $m_{chargino} < 0.8$ TeV.
  • Direct WIMP detection has a 68% probability of landing between $10^{-10}$ and $10^{-8}$ pb, with a nearly $tan beta$-independent band at about $1.6 times 10^{-8}$ pb produced by the focus-point region and expected to be probed first.
  • A Tevatron measurement of $BR(B_s to mu^+ mu^-)$ above roughly $2 times 10^{-8}$ would strongly disfavor large $m_0$, including the focus-point region.
  • Removing the $(g-2)_mu$ measurement from the likelihood relaxes the constraint on $m_0$ almost entirely and shifts the predicted $BR(B_s to mu^+ mu^-)$ and $sigma_p^{SI}$ downward, lowering the chance that sleptons and squarks fall inside LHC reach.
  • Because the posterior and the quality of fit disagree, a Bayesian preference for large $tan beta$ should not be read as a statement that the best-fitting points live there.

Reading between the lines

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

  • The two-branch structure of the $BR(B_s to mu^+ mu^-)$ versus $sigma_p^{SI}$ plane gives a clean experimental discriminator: a joint signal (or joint null) in a rare-decay measurement and a direct-detection experiment would test whether the focus-point branch is really populated, which no single measurement can do.
  • If the focus-point reading is correct, a null result near $10^{-8}$ pb in the next direct-detection round would not kill the model but would push the preferred region toward lower $m_0$, where sleptons are more accessible to the LHC, so the two searches cover each other's blind spots in a way that could be planned for deliberately.
  • The gap between posterior preference and best-fit quality for $tan beta$ is a warning about how global-fit results get quoted generally: a parameter can look pinned by a Bayesian scan purely because the allowed volume shrinks, and the check is to report the profile likelihood alongside the marginal posterior.
  • The authors' own admission that the relic-density error in the focus-point region may reach order unity suggests the honest test is not a better central value but a propagated error band on $Omega_chi h^2$ computed region by region; that would sharpen or dissolve the large-$m_0$ tail that drives the high-$sigma_p^{SI}$ prediction.
Original abstract

We perform a comprehensive exploration of the Constrained MSSM parameters using a Markov Chain Monte Carlo technique and a Bayesian analysis. We compute superpartner masses and other collider observables as well as a cold dark matter abundance, and compare them with experimental data. We include uncertainties arising from theoretical approximations as well as from residual experimental errors on relevant SM parameters. We delineate probability distributions of the CMSSM parameters, the collider and cosmological observables as well as a dark matter direct detection cross section. The 68% probability intervals of the CMSSM parameters are: 0.52 TeV < m_{1/2} < 1.26 TeV, m_0 <2.10 TeV, -0.34 TeV < A_0 < 2.41 TeV and 38.5< tan(beta) <54.6. Generally, large fractions of high probability ranges of the superpartner masses will be probed at the LHC. We highlight a complementarity between LHC and WIMP dark matter searches in exploring the CMSSM parameter space. We further expose a number of correlations among the observables, in particular between BR(B_s \to \mu^+ \mu^-) and BR({\bar B}\to X_s\gamma) or sigma_p^{SI}. Once SUSY is discovered, this and other correlations may prove helpful in distinguishing the CMSSM from other supersymmetric models. The robustness of our results is investigated in terms of the assumed ranges of CMSSM parameters and the effect of the (g-2)_mu anomaly which shows some tension with the other observables. We find that the results for m_0, and the observables which strongly depend on it, are sensitive to our assumptions, while our conclusions for the other variables are robust.

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, 8,447 input tokens, $0.0004 reported.

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

The abstract, sentence by sentence

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

We perform a comprehensive exploration of the Constrained MSSM parameters using a Markov Chain Monte Carlo technique and a Bayesian analysis. We compute superpartner masses and other collider observables as well as a cold dark matter abundance, and compare them with experimental data. We include uncertainties arising from theoretical approximations as well as from residual experimental errors on relevant SM parameters. We delineate probability distributions of the CMSSM parameters, the collider and cosmological observables as well as a dark matter direct detection cross section. The 68% probability intervals of the CMSSM parameters are: 0.52 TeV < m_{1/2} < 1.26 TeV, m_0 <2.10 TeV, -0.34 TeV < A_0 < 2.41 TeV and 38.5< tan(beta) <54.6. Generally, large fractions of high probability ranges of the superpartner masses will be probed at the LHC. We highlight a complementarity between LHC and WIMP dark matter searches in exploring the CMSSM parameter space. We further expose a number of correlations among the observables, in particular between BR(B_s \to \mu^+ \mu^-) and BR({\bar B}\to X_s\gamma) or sigma_p^{SI}. Once SUSY is discovered, this and other correlations may prove helpful in distinguishing the CMSSM from other supersymmetric models. The robustness of our results is investigated in terms of the assumed ranges of CMSSM parameters and the effect of the (g-2)_mu anomaly which shows some tension with the other observables. We find that the results for m_0, and the observables which strongly depend on it, are sensitive to our assumptions, while our conclusions for the other variables are robust.

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 perform a comprehensive exploration of the Constrained MSSM parameters using a Markov Chain Monte Carlo technique and a Bayesian analysis. method 0.51 0.97 0.04 asserted 0.99
s1 We compute superpartner masses and other collider observables as well as a cold dark matter abundance, and compare them with experimental data. contribution 0.68 0.96 0.05 asserted 0.99
s2 We include uncertainties arising from theoretical approximations as well as from residual experimental errors on relevant SM parameters. method 0.90 0.88 0.02 asserted 0.99
s3 We delineate probability distributions of the CMSSM parameters, the collider and cosmological observables as well as a dark matter direct detection cross section. contribution 0.58 0.95 0.08 asserted 1.00
s4 The 68% probability intervals of the CMSSM parameters are: 0.52 TeV < m_{1/2} < 1.26 TeV, m_0 <2.10 TeV, -0.34 TeV < A_0 < 2.41 TeV and 38.5< tan(beta) <54.6. result 1.00 0.95 0.21 asserted 1.00
s5 Generally, large fractions of high probability ranges of the superpartner masses will be probed at the LHC. interpretation 0.85 0.78 0.13 hedged possibility 0.65
s6 We highlight a complementarity between LHC and WIMP dark matter searches in exploring the CMSSM parameter space. interpretation 0.55 0.86 0.18 asserted 1.00
s7 We further expose a number of correlations among the observables, in particular between BR(B_s \to \mu^+ \mu^-) and BR({\bar B}\to X_s\gamma) or sigma_p^{SI}. result 0.98 0.93 0.07 asserted 1.00
s8 Once SUSY is discovered, this and other correlations may prove helpful in distinguishing the CMSSM from other supersymmetric models. interpretation 0.91 0.24 0.07 hedged possibility 1.00
s9 The robustness of our results is investigated in terms of the assumed ranges of CMSSM parameters and the effect of the (g-2)_mu anomaly which shows some tension with the other observables. limitation 0.33 0.91 0.06 asserted 0.91
s10 We find that the results for m_0, and the observables which strongly depend on it, are sensitive to our assumptions, while our conclusions for the other variables are robust. result 0.69 0.96 0.24 asserted 1.00

What kind of paper the classifier reports this is: simulation study 1.00 · evidence from simulation or computation 0.99.

What was examined, and what was not

sources read
  • arXiv listing abstract v3
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

Figures
Figure 1. An illustration of the likelihood function p(σ, ξlim|θ, τ) used for a quantity for which only a lower bound is available, including the theoretical uncertainty τ and setting the experimental error σ = 0. The dashed green line is the sharp 95% CL bound (ξlim = 1), the solid red curve includes the theoretical uncertainty of τ = 0.05 which smears out the limit. confidence limit only, without the possibility of deriving a form of the likelihood analogous to (3.5). In the absence of fuller information about the experimental value of σ, we can take in (3.5) the limit σ ≪ τ and simply use…
Figure 2. The 2–dimensional probability densities in the planes spanned by the CMSSM pa￾rameters: m1/2, m0, A0 and tan β for the “4 TeV range” analysis (see
Figure 3. The 2–dimensional relative probability densities as in
Figure 4. The 1–dimensional relative probability densities p(θi |d) for each of the CMSSM param￾eter, θi = m0, m1/2, A0 and tan β. All other parameters have been marginalized over. The two curves compare the results for two different prior ranges (see
Figure 6. The 2–dim relative probability density p(Ωχh 2 , θi |d), where θi = m0, m1/2 and tan β. Note that the measured value of Ωχh 2 , Eq. (3.10), has been included in computing the relative probability density. m0 (TeV) δ aSUSY µ × 1010 1 2 3 4 0 10 20 30 m1/2 (TeV) 0.5 1 1.5 2 0 10 20 30 tanβ 20 40 60 0 10 20 30 Relative probability density 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1
Figure 10. An illustration of the tensions between different observables. While for mb(mb)MS and Ωχh 2 the likelihood and the posterior pdf agree very well (top panels), for sin2 θeff and δaSUSY µ (bottom panel) they exhibit a substantial discrepancy. The latter case is a clear sign that the other constraints are strongly pulling away the posterior pdf from the value preferred by the anomalous magnetic moment measurement. rameters in the problem after we have seen the data and given our choice of priors, as we have explained in some detail in section 2.1. It is important to remark that regions
Figure 11. The quality of fit in the planes spanned by the CMSSM parameters: m1/2, m0, A0 and tan β for the “4 TeV range” scan. This figure should be compared with
Figure 13. Top panel: the 2–dimensional relative probability density p(mχ, σSI p |d), with the contours containing 68% and 95% probability also marked. Bottom panel: the mean quality of fit (likelihood) with the same 68% and 95% probability contours as in the top panel to facilitate the comparison. Current 90% experimental upper limits are also shown. be the first to be probed by DM search experiments. Finally, after marginalizing over all other parameters, we obtain the following 1– – 33 –
Figure 14. The 2–dimensional relative probability density p(ξi , ξj |d) for the pairs of observables marked along horizontal and vertical axes. dimensional regions encompassing 68% and 95% of the total probability: 1.0 × 10−10 pb < σSI p < 1.0 × 10−8 pb (68% region), 0.5 × 10−10 pb < σSI p < 3.2 × 10−8 pb (95% region). (4.4) Currently running experiments (most notably CDMS–II but also Edelweiss–II and ZEPLIN–II) should be able to reach down to a few ×10−8 pb, on the edge of exploring the FP region. A future generation of “one–tonne” detectors are expected to reach down to σ SI p ∼ > 10−10 pb,…
Figure 15. Values of mh (upper left), Ωχh 2 (upper right), σ SI p (lower left) and δaSUSY µ (lower right) in the (m1/2, m0) plane for samples uniformly selected from our MC chains. While σ SI p and δaSUSY µ show a clear correlation with the values of (m1/2, m0), the CDM abundance Ωχh 2 can take any value within the 2σ WMAP range in the plane. Narrow “WMAP strips” in the (m1/2, m0) plane for fixed values of tan β and A0 thus disappear when these parameters are allowed to vary and all other parameters are correctly marginalized over. The values of σ SI p (bottom left panel of
Figure 16. The 1–dimensional relative posterior pdf’s as in
Figure 18. The probability for the superpartner masses to lie below a given mass. All other parameters have been marginalized over. marized in

Showing 12 of the paper's 18 figures. See all figures in the PDF ↗ Figures reproduced from arXiv:hep-ph/0602028 by the authors.

References
How later papers used it

74 later papers cite it

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

Pith read 18 of the 129 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 8 of 29 papers with this label · 4 by the same authors

By the same authors

Uses its data or values 1 paper
How the literature uses it

73 arXiv papers cite this in 127 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 · 49 sentences without a confident label

Uses its method or tool 28 papers
Uses its data or values 1 paper
Compares with it 8 papers
  • The sparticle mass distributions for priors that are flat in $tanβ$ were displayed in Refs. [ref] and show a spread up to much higher values of the masses.

    Natural Priors, CMSSM Fits and LHC Weather Forecasts · CMSSM Fits With the New Priors · 0.88

  • [ref] However, the precise value of this upper bound turns out to be irrelevant, as the region of very large $tanβ$ is strongly suppressed (see the discussion after eq.

    MSSM Forecast for the LHC · High-energy vs Low-energy regions and the EW breaking · 0.97

  • Despite our updated constraints and additional observables, the flat $tanβ$ posterior looks indistinguishable to the eye to previous determinations [ref] , where the connection between a preference for the focus point and good direct detection prospects were pointed out.

    Panglossian Prospects for Detecting Neutralino Dark Matter in Light of Natural Priors · Direct Detection · 0.99

Cites it as background 25 papers

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.