How alive is constrained SUSY really?

How alive is constrained SUSY really?

Philip Bechtle Klaus Desch Herbert K. Dreiner Matthias Hamer Michael Krämer Ben O’Leary Werner Porod Björn Sarrazin Tim Stefaniak Mathias Uhlenbrock Peter Wienemann Physikalisches Institut, University of Bonn, Germany Bethe Center for Theoretical Physics, University of Bonn, Germany Centro Brasileiro de Pesquisas Fisicas, Rio de Janeiro, Brazil Institute for Theoretical Particle Physics and Cosmology, RWTH Aachen, Germany SLAC National Accelerator Laboratory, Stanford University, USA Institut für Theoretische Physik und Astrophysik, University of Würzburg, Germany Santa Cruz Institute for Particle Physics, University of California, Santa Cruz, USA

Constrained supersymmetric models like the CMSSM might look less attractive nowadays because of fine tuning arguments. They also might look less probable in terms of Bayesian statistics. The question how well the model under study describes the data, however, is answered by frequentist -values. Thus, for the first time, we calculate a -value for a supersymmetric model by performing dedicated global toy fits. We combine constraints from low-energy and astrophysical observables, Higgs boson mass and rate measurements as well as the non-observation of new physics in searches for supersymmetry at the LHC. Using the framework Fittino, we perform global fits of the CMSSM to the toy data and find that this model is excluded at more than confidence level.

supersymmetry phenomenology, LHC

1 Introduction

As a fully renormalizable theory beyond the Standard Model (SM) supersymmetry Wess:1974tw () makes a wide range of predictions, allowing to probe it with many different measurements. Indirect constraints from low energy measurements, the discovery of an SM-like Higgs boson, direct searches for supersymmetric particles and additional Higgs bosons as well as astrophysical observations constrain the parameter space.

As one of the simplest supersymmetric models the Constrained Minimal Supersymmetric Standard Model (CMSSM) Nilles:1983ge () has been studied extensively in light of this data over the last years. It has five free parameters: The ratio of the vacuum expectation values of the two Higgs doublets (), the sign of the Higgs mixing parameter (), and universal soft-breaking parameters defined at the GUT scale for scalar masses (), gaugino masses (), and trilinear couplings (). Global fits are used to determine the parameter region which is in best agreement with these measurements.

It has become clear from these studies that the non-observation of convincing hints of new physics at the LHC becomes challenging for this model once constraints from low energy observables are taken into account. However, it has so far been an open question what this means quantitatively for the validity of the model. In this work, we address this question for the first time systematically using global toy fits to calculate a -value for the CMSSM.

The remainder of this article is organized as follows: Section 2 gives a brief review of the status of the CMSSM. Section 3 introduces the -value. Section 4 details which measurements we use for our analysis, while Section 5 summarizes the codes used to obtain the corresponding model predictions. The technique to compare the measurements and their predictions is described in Section 6. Profile likelihood based results are shown in Section 7, while results from toy fits are given in Section 8. We conclude in Section 9.

2 What has happened so far

Before the start of the LHC, global fits of the CMSSM preferred relatively low values of the mass parameters. This was mainly driven by the anomalous magnetic moment of the muon and one of the reasons why it was anticipated to find supersymmetry early at the LHC. Since this has not happened yet, some tension arises between low energy observables and exclusions from direct searches. While still giving an excellent fit, this effect could already be seen after the first searches for supersymmetry by ATLAS or CMS with an integrated luminosity of  pb at  TeV were included in the fit Bechtle:2011dm (). With increasing luminosity up to at 8 TeV this tension increased significantly Bechtle:2012zk (); Bechtle:2013mda ().

While there is also no convincing hint for new physics from direct and indirect astrophysical searches for dark matter, bounds from these experiments are currently easily fulfilled by the CMSSM. Similarly, while no extended Higgs sector has been found, the discovery of an SM-like Higgs boson is in agreement with CMSSM predictions. The inclusion of numerous SM-like Higgs rate measurements even improves the ratio of the fit Bechtle:2013mda (). This is due to the fact that the CMSSM, when forced to large supersymmetric masses by direct searches, enters the so-called decoupling limit where an SM-like Higgs is predicted.

3 The missing piece

The status described in Section 2 seems to leave us in an unsatisfactory situation after the end of the first run of the LHC. On the one hand, there is some tension between low energy measurements and direct searches, on the other hand, the model describes many measurements very well and there is no single clear piece of evidence against it. This hampers a final judgement. Here, we escape from this dead end by calculating the -value of the CMSSM. In high energy physics, the -value is the standard measure for the agreement between model and data. It is defined as the probability to obtain an agreement at least as bad as the one observed assuming the model is true.

Determining the -value in global fits is challenging because the -function used in the fit is not necessarily -distributed due to non-Gaussian uncertainties and non-linear dependences of the observables on the model parameters. Therefore, it is not possible to calculate the -value by simply integrating over the -distribution for the appropriate number of degrees of freedom. Instead, it is necessary to determine the shape of the distribution by performing toy fits as will be described in Section 6. Using the framework Fittino Bechtle:2004pc (), this makes it possible to calculate a confidence level for an exclusion of the CMSSM for the first time.

4 Measurements

We use the same set of low energy observables as in Ref. Bechtle:2012zk (), but with updated measurements as given in Table 1. They include the anomalous magnetic moment of the muon , the effective weak mixing angle , the masses of the top quark and boson, the oscillation frequency , as well as the branching ratios , , and .

In order to have a good dark matter candidate, we require the lightest supersymmetric particle to be a neutralino. We use the dark matter relic density as measured by the Planck collaboration Ade:2013zuv () and bounds on the spin-independent WIMP-nucleon scattering cross section as measured by the LUX experiment Akerib:2013tjd ().

Concerning direct searches at the LHC, we use the final state with jets and missing transverse momentum Aad:2014wea () because this channel is by far the most dominant one in the region preferred by the low energy observables. In addition, we make sure that the LEP bound on the chargino mass,  GeV, is fulfilled LEPSUSYWG (). We apply the limits from Higgs boson searches as implemented in HiggsBounds 4.1.1 Bechtle:2013wla (). Using the program HiggsSignals 1.2.0, we include Higgs signal rate measurements in the channels , , , , and from both ATLAS and CMS. Higgs boson mass measurements in the channels and are used as well ( see Ref. Bechtle:2013xfa () and references therein).

Bennett:2006fi (); Davier:2010nc ()
ALEPH:2005ab ()
 GeV ATLAS:2014wva ()
 GeV Group:2012gb ()
Beringer:1900zz ()
CMSandLHCbCollaborations:2013pla ()
Amhis:2012bh ()
Beringer:1900zz ()
Table 1: Low energy observables used in the fit.

5 Model predictions

1 GeV
Table 2: Theoretical uncertainties of the low energy observables used in the fit.

To calculate the corresponding model predictions of the observables, we employ the following public codes: The spectrum calculation is done using SPheno 3.2.4 Porod:2003um (); Porod:2011nf (). Properties of the Higgs bosons as well as , , and are taken from FeynHiggs 2.10.1 Heinemeyer:1998yj (), which – compared to FeynHiggs 2.9.5 and earlier versions – contains a significantly improved calculation of the Higgs boson mass Hahn:2013ria (). B-physics branching ratios are calculated by SuperIso 3.3 Mahmoudi:2008tp (). Concerning astrophysical observables we use MicrOMEGAs 3.6.9 Belanger:2001fz (); Belanger:2004yn () for the dark matter relic density and DarkSUSY 5.0.5 Gondolo:2004sc () via the interface program AstroFit AstroFit () for the direct detection cross section. For the calculation of expected event numbers at the LHC, we use the Monte Carlo event generator Herwig++ Bahr:2008pv () and a carefully tuned and validated version of the fast parametric detector simulation Delphes Ovyn:2009tx (). We reweight the events depending on their production channel according to NLO cross sections obtained from Prospino Beenakker:1996ed (). Renormalization and factorization scales have been chosen such that NLO+NLL results are reproduced. For all predictions we take theoretical uncertainties into account, most of them are parameter dependent. For the low energy observables, they are given in Table 2. For the dark matter relic density we assume a theoretical uncertainty of , for the neutralino-nucleon cross section , for the Higgs boson mass prediction , and for Higgs rates we use the uncertainties given by the LHC Higgs Cross Section Working Group Heinemeyer:2013tqa ().

6 Techniques

The framework Fittino compares measurements and model predictions using a -function given by

where is the vector of measurements, the vector of corresponding model predictions, and the covariance matrix containing all experimental and theoretical uncertainties and their correlations. The term contains the contributions from Higgs limits, LUX and the direct searches. This function is used in a Markov chain Monte Carlo in order to sample the CMSSM parameter space efficiently. We vary the continuous CMSSM parameters , , and and fix in order to get positive supersymmetric contributions to the anomalous magnetic moment of the muon. We fit the top quark mass as an additional free parameter.

In order to determine the -value, we perform toy fits. As described in section 3 the -value is given by

where is the minimal observed and is the probability to get a minimal between and if the model is true. Toy fits are needed in order to determine the a priori unknown distribution . We smear the observable values within their uncertainties around their predictions at the best fit point and fit the model to each of these toy measurements. In the limit of large statistics, the -value is then by definition given as the fraction of toy fits that have a minimal at least as bad as the one observed.

7 Profile likelihood based results

Figure 1: Preferred parameter region in the (, ) plane.
Figure 2: Preferred parameter region in the (, ) plane.

We find a minimal of with degrees of freedom. The best fit point is shown in Figures 1 and 2 as a black star, together with the projections of the and regions in red and blue, respectively. Figure 1 shows two disjoint areas in the (, ) plane, which are mainly characterized by different dominant dark matter annihilation mechanisms. At lower masses, stau coannihilation and the Higgs funnel are present. At higher masses, the focus point region and chargino coannihilation are relevant. While the high mass region has been disfavored in our previous fits Bechtle:2013mda (), it is now even part of the region. This is due to the new Higgs boson mass calculation in FeynHiggs, which allows larger Higgs boson masses when going to large supersymmetric mass scales. We did not scan the parameter space for values of larger than  TeV. Figure 2 shows the (, ) plane. In the low mass region negative values of are preferred and a wide range of values of is allowed, whereas the high mass region prefers positive values of and larger values of above 40.

Figure 3: Predicted mass spectrum.

The corresponding mass spectrum is shown in Figure 3. Due to the high mass region, squark masses of about 10 TeV are now allowed at . At the best fit point, squarks and gluinos have masses of about 2 TeV, the heavy Higgs boson masses are about 1.5  TeV. Charginos, neutralinos and sleptons stay lighter.

Figure 4: Distribution of best fit values of obtained from toy fits.

8 Toy fit results

We repeat the fit described above about 1000 times with toy measurements. The obtained distribution of best fit points is shown in in Figure 4. Most toy fits have their minimum in the low mass region, but some of them end up in the high mass region.

Figure 5: Distribution of predicted values of at the best fit points obtained from toy fits.

As an example, Figure 5 shows the corresponding distribution of predicted values of at the best fit points. While the measurement of has been smeared according to a Gaussian distribution, this distribution clearly shows a non-Gaussian behavior. Its long tail shows that there is more space in the relevant CMSSM parameter space to go to larger values of than to lower values, which would be closer to the actual measurement. Nevertheless, the toy values are only in slight conflict with this measurement.

Figure 6: Distribution of predicted values of at the best fit points obtained from toy fits.

The situation is different for as shown in Figure 6. There is some variation in the best fit predictions of the toy fits, showing that in the relevant CMSSM parameter space still has some influence on the best fit point. However, this variation is small compared to the disagreement with the actual measurement. As a result, contributes to the much less in the toy fits than in the fit to the actual measurements. Because of this, a minimal as bad as the one observed is found only in a small fraction of the toy fits.

This can be seen in Figure 7, which shows the distribution of minimal values and in comparison a -distribution for 22 degrees of freedom. Since the latter is shifted to the right, it overestimates the goodness of fit. Shown in blue is the fraction of toy fits giving a as bad as the one observed, resulting in a -value of , excluding the CMSSM with confidence.

Figure 7: Distribution of minimal -values obtained from toy fits. A -distribution for 22 degrees of freedom is shown for comparison.

9 Conclusion

The first run of the LHC has improved our knowledge significantly: Among other results, there are the discovery of an SM-like Higgs boson, the first measurement of , an improved top quark mass world average, and stronger constraints on new physics. In this work, all this information has been combined with the results from low energy and astrophysical experiments. In this way, we could show that the CMSSM as one of the simplest and most popular models is excluded with more than confidence. More details of this analysis will be given in an upcoming publication.

Of course this is by no means a statement about supersymmetry in general. However, it sets the scene for applying the described method of global toy fits to more general models, where, for instance, the strong and electroweak sectors are disentangled. In this way, it will be possible to make full use of the potential of the upcoming results from the second run of the LHC.


We are grateful to the organizers of the ICHEP 2014 conference for the opportunity to present this work. The work of MK was supported by the Deutsche Forschungsgemeinschaft through the collaborative research centre SFB-TR9 ”Computational Particle Physics”, by the German Federal Ministry of Education and Research BMBF, and by the U.S. Department of Energy under contract DE-AC02-76SF00515. The work of TS was partially supported by a Feodor-Lynen research fellowship sponsored by the Alexander von Humboldt Foundation.


Comments 0
Request Comment
You are adding the first comment!
How to quickly get a good reply:
  • Give credit where it’s due by listing out the positive aspects of a paper before getting into which changes should be made.
  • Be specific in your critique, and provide supporting evidence with appropriate references to substantiate general statements.
  • Your comment should inspire ideas to flow and help the author improves the paper.

The better we are at sharing our knowledge with each other, the faster we move forward.
The feedback must be of minimum 40 characters and the title a minimum of 5 characters
Add comment
Loading ...
This is a comment super asjknd jkasnjk adsnkj
The feedback must be of minumum 40 characters
The feedback must be of minumum 40 characters

You are asking your first question!
How to quickly get a good answer:
  • Keep your question short and to the point
  • Check for grammar or spelling errors.
  • Phrase it like a question
Test description