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
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
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
- 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.
The abstract, sentence by sentence
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
| # | sentence | role | the authors' own work | states the main conclusion | strength |
|---|---|---|---|---|---|
| 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 was examined, and what was not
- arXiv listing abstract v3
- 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
Figures
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
5 later papers by other authors are listed below: 4 used its method or tool and 1 used its data or values. 4 more papers by the same authors are listed: 4 used its method or tool.
Uses its method or tool 8 of 29 papers with this label · 4 by the same authors
-
Used the Markov Chain Monte Carlo routines of Superbayes from this paper, which generate Markov chains by means of the Metropolis algorithm, in its FORTRAN code.
A novel determination of the local dark matter density (2009) · 1.00 · paraphrased: the licence does not allow quoting
-
Accounted for positive measurements with a Gaussian likelihood and implemented upper limits as described in this paper.
Explaining the Fermi Galactic Centre Excess in the CMSSM (2015) · 1.00 · paraphrased: the licence does not allow quoting
-
Followed this paper's approach for 95% CL exclusion limits, formulating a likelihood that is a step function at the observed limit smeared by a 5% gaussian uncertainty.
Disfavouring Electroweak Baryogenesis and a hidden Higgs in a CP-violating Two-Higgs-Doublet Model (2016) · 1.00 · paraphrased: the licence does not allow quoting
-
Constructed smeared step-function likelihoods at 95% CL for generic limits implemented in another work, as in this paper.
The Power of Genetic Algorithms: what remains of the pMSSM? (2018) · 0.95 · paraphrased: the licence does not allow quoting
By the same authors
-
Referred the reader to this paper for a detailed description of the likelihood function.
Dark Matter Searches: The Nightmare Scenario (2011) (same authors) · 0.82 · paraphrased: the licence does not allow quoting
-
Included the uncertainty in the likelihood following this paper's procedure for the treatment of upper bounds.
Bayesian Implications of Current LHC Supersymmetry and Dark Matter Detection Searches for the Constrained MSSM (2012) (same authors) · 1.00 · paraphrased: the licence does not allow quoting
-
Following this paper's procedure, implemented positive measurements through a Gaussian likelihood and upper or lower limits through an error function smeared with theory and, when available, experimental error.
The Constrained NMSSM with a 126 GeV Higgs boson: A global analysis (2012) (same authors) · 1.00 · paraphrased: the licence does not allow quoting
-
Following this paper's procedure, implemented positive measurements through a Gaussian likelihood with experimental and theoretical uncertainties added in quadrature.
Two ultimate tests of constrained supersymmetry (2013) (same authors) · 1.00 · paraphrased: the licence does not allow quoting
Uses its data or values 1 paper
-
Based its updated bounds upon the sparticle spectrum bounds collected in this paper.
Natural Priors, CMSSM Fits and LHC Weather Forecasts (2007) · 0.95 · paraphrased: the licence does not allow quoting
How the literature uses it
73 arXiv papers cite this in 127 sentences
Papers citing it · 49 sentences without a confident label
Uses its method or tool 28 papers
-
We use our updated and modified version of the [ref] computer program to perform four scans of the CMSSM:
Bayesian Implications of Current LHC and XENON100 Search Limits for the Constrained MSSM · The statistical framework · 0.98
-
To this end, we follow the strategy outlined in Refs. [ref] .
Constrained MSSM favoring new territories: The impact of new LHC limits and a 125 GeV Higgs boson · The framework · 1.00
-
We refer the reader to [ref] for a detailed description of the likelihood function.
Dark Matter Searches: The Nightmare Scenario · Cosmological and collider constraints · 0.82
Uses its data or values 1 paper
-
Bounds upon the sparticle spectrum have been updated and are based upon the bounds collected in Ref. [ref] .
Natural Priors, CMSSM Fits and LHC Weather Forecasts · The Likelihood · 0.95
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
-
In [ref] , Bayesian statements were made about the posterior probability density of the CMSSM, after indirect data had been taken into account.
Natural Priors, CMSSM Fits and LHC Weather Forecasts · Introduction · 0.90
-
More recently, many studies have been performed to extract the preferred values for the CMSSM parameters using low-energy precision data, bounds from astrophysical observables and flavour related observables [ref] .
Prediction for the Lightest Higgs Boson Mass in the CMSSM using Indirect Experimental Constraints · Introduction · 0.85
-
Recently, Bayesian parameter estimation techniques using the Markov Chain Monte Carlo (MCMC) sampling have been applied to the study of mSUGRA, performing a multi–dimensional Bayesian fit to indirect constraints [ref] .
Bayesian Selection of sign(mu) within mSUGRA in Global Fits Including WMAP5 Results · Introduction · 0.93
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.