Constrained Supersymmetry after the Higgs Boson Discovery: A global analysis with Fittino
Philip Bechtle, Klaus Desch, Herbert K. Dreiner, Matthias Hamer, Michael Krämer, Ben O'Leary, Werner Porod, Xavier Prudent, Björn Sarrazin, Tim Stefaniak, Mathias Uhlenbrock, Peter Wienemann
BBONN-TH-2013-19
Constrained Supersymmetry after the Higgs BosonDiscovery: A global analysis with F
ITTINO
Philip Bechtle, Klaus Desch, Björn Sarrazin, Mathias Uhlenbrock, PeterWienemann
Physikalisches Institut, Bonn University, GermanyE-mail: [email protected] , [email protected] , [email protected] , [email protected] , [email protected] Herbert K. Dreiner, Tim Stefaniak ∗ Physikalisches Institut and Bethe Center for Theoretical Physics, Bonn University, GermanyE-mail: [email protected] , [email protected] Matthias Hamer
II. Physikalisches Institut, University of Göttingen, Göttingen, GermanyE-mail: [email protected]
Michael Krämer
Institute for Theoretical Particle Physics and Cosmology, RWTH Aachen, GermanyE-mail: [email protected]
Werner Porod, Ben O’Leary
Institut für Theoretische Physik und Astrophysik, University of Würzburg, GermanyE-mail: [email protected] , [email protected] Xavier Prudent
Institut für Kern- und Teilchenphysik, TU Dresden, Dresden, GermanyE-mail: [email protected]
We present preliminary results from the latest global fit analysis of the constrained minimal super-symmetric standard model (CMSSM) performed within the F
ITTINO framework. The fit includeslow-energy and astrophysical observables as well as collider constraints from the non-observationof new physics in supersymmetric searches at the LHC. Furthermore, the Higgs boson mass andsignal rate measurements from both the LHC and Tevatron experiments are included via the pro-gram H
IGGS S IGNALS . Although the LHC exclusion limits and the Higgs mass measurements puttight constraints on the viable parameter space, we find an acceptable fit quality once the Higgssignal rates are included.
The European Physical Society Conference on High Energy Physics -EPS-HEP201318-24 July 2013Stockholm, Sweden ∗ Speaker. c (cid:13) Copyright owned by the author(s) under the terms of the Creative Commons Attribution-NonCommercial-ShareAlike Licence. http://pos.sissa.it/ a r X i v : . [ h e p - ph ] O c t onstrained SUSY after the Higgs Boson Discovery Tim Stefaniak
1. Introduction
Supersymmetry (SUSY) provides an elegant solution to the fine tuning problem of the Stan-dard Model (SM) if supersymmetric particles are realized around the TeV scale. Then, the SUSYparameter space underlies various constraints from SM and astrophysical observables as well asdirect sparticle searches at colliders. Moreover, since the Higgs boson masses and couplings arepredictions of the theory, the mass and signal strength measurements of the recently discoveredHiggs boson [1] severely constrain the parameter space. We perform a global fit of the constrainedMinimal Supersymmetric Standard Model (CMSSM) to these observables in order to answer thefollowing questions: ( i ) What is the allowed SUSY model parameter space after including allavailable and relevant observables and constraints? ( ii ) To what extend are the observables andconstraints in mutual agreement?The CMSSM is defined by a few parameters at the grand unification (GUT) scale ∼ GeV:Universal soft-breaking mass parameters for the scalars and gauge fermions, M and M / , respec-tively, a universal soft-breaking trilinear coupling, A , the ratio of the vevs of the two Higgs dou-blets, tan β , and the sign of the Higgs mixing parameter, sgn ( µ ) , which we fix to be positive in ourstudy. In addition to the four free CMSSM fit parameters, we allow for the top quark mass m t as anadditional nuisance parameter.In a previous F ITTINO analysis [2] we found rather grim prospects for a discovery of super-symmetric effects within the CMSSM for the near future: The sparticles and the remaining Higgsspectrum was most likely beyond the LHC 8 TeV reach. The light Higgs boson signal rates andthe branching ratio for B s → µ µ were predicted to be close to their SM prediction. Furthermore,no dark matter (DM) signal was expected in current direct and indirect searches. This presenta-tion updates and extends the previous analysis [2] by an implementation of the Higgs boson signalrate measurements. Here, we show first preliminary results of the on-going work, which will bepresented in more detail in a future publication.
2. The F
ITTINO framework and technical implementation
The F
ITTINO framework [3, 4] incorporates an auto-adaptive Markov Chain Monte Carlo(MCMC) algorithm, allowing to efficiently sample the multi-dimensional SUSY parameter space.The results presented here are based on a high statistics sample with ∼ scan points and apurely frequentist interpretation. For each scan point, the SUSY particle spectrum is calculatedwith SP HENO -3.1.11 [5], followed by a dedicated evaluation of the Higgs masses and couplingsas well as the SUSY contribution to the anomalous magnetic moment of the muon, a SUSY µ , withF EYN H IGGS -2.9.4 [6]. We use S
UPER I SO -3.3 [7] for predictions of the heavy flavor observ-ables and M ICR
OMEGA S -2.4.5 [8] for the DM relic density calculation. A STRO F IT [9] andD ARK
SUSY-5.0.5 [10] are employed to include direct detection limits from DM searches. We in-clude exclusion limits from Higgs boson searches at LEP, Tevatron and LHC using H
IGGS B OUNDS -3.8.1 [11, 12]. The LEP constraints on Higgs boson production are incorporated as a reconstructed χ likelihood [2, 12]. The χ contribution from the Higgs boson signal rate and mass measurementsfrom the Tevatron and LHC experiments is evaluated with H IGGS S IGNALS -1.0.0 [13].2 onstrained SUSY after the Higgs Boson Discovery
Tim Stefaniak χ L HC M500 1000 1500 2000 2500 ( G e V ) / M = 8 TeVs L dt = 20.3 / fb, ∫ Atlas median expected limitAtlas observed limitFittino implementation (a) Comparison of our reconstructed likelihoodwith the ATLAS 95% C.L. exclusion contours. (b) Relative difference of the grid-interpolated χ andthe true χ contribution in the ( A , tan β ) plane, exam-plarily shown for ( M , M / ) = ( , ) GeV.
Figure 1:
Implementation of the exclusion limit from the full hadronic ATLAS SUSY search [17].
3. Experimental constraints
A detailed description of the experimental measurements and constraints (including refer-ences) can be found in Ref. [2]. Here, we briefly discuss the new or updated observables in thecurrent fit and their implementation.We include the LHCb measurement of BR ( B s → µ µ ) = ( . ± . ± . ) · − [14], theupdated Belle measurement of BR ( B → τν ) = ( . ± . ± . ± . ) · − [15] and the newrelic density abundance determination from the Planck collaboration, Ω CDM h = . ± . ± . ( B → τν ) measurement is in better agreementwith both the SM and CMSSM predictions than the old value used in Ref. [2].We update our implementation of the LHC constraints from searches for direct sparticle pairproduction (see Ref. [2] for a description) to accommodate the latest ATLAS results in the fullhadronic channel using 20 . − of data at √ s = ( M , M / ) plane. Wefind that in the parameter region around the exclusion limit, the contribution from light ˜ t pairsto the signal yield is non-negligible. This introduces a dependence on the trilinear soft-breakingparameter A , which strongly influences the ˜ t mass in the renormalization group (RG) evolution.We found that the uncorrected χ contribution (using only the ( M , M / ) acceptance grid) can besmaller than the true χ contribution (obtained from full Monte-Carlo simulation) by up to 80%for | A | (cid:38) − ( A , tan β ) planefor ( M , M / ) values along the exclusion contour. After this correction, the differences betweenthe acceptance grid based χ and the true χ are typically (cid:46) IGGS S IGNALS -1.0.0 (see Ref. [13] for a description and references).Moreover, we include the four Higgs mass measurements from the ATLAS and CMS H → γγ and H → ZZ ( ∗ ) → (cid:96) analyses. Note, that H IGGS S IGNALS treats the uncertainties for the signal rateand Higgs mass predictions as well as the luminosity uncertainty as fully correlated Gaussian un-3 onstrained SUSY after the Higgs Boson Discovery
Tim Stefaniak (GeV) M0 200 400 600 800 1000120014001600180020002200 ( G e V ) / M
2D 95% CL 1D 68% CL (a) ( M , M / ) plane. (GeV) A 5000 4000 3000 2000 1000 0 1000 β t a n
2D 95% CL 1D 68% CL (b) ( A , tan β ) plane. Figure 2:
Preferred two-dimensional CMSSM parameter regions. We profile over the remaining fit param-eters. The red and blue areas correspond to the one-dimensional 68% C.L. ( ∆ χ =
1) and two-dimensional95% C.L. ( ∆ χ = .
18) parameter region, respectively. The black star indicates the best-fit point. certainties. We use the same estimates for the Higgs boson production cross section and branchingratio uncertainties for the SM and the CMSSM and assume a theoretical (Gaussian) uncertainty of3 GeV on the light Higgs boson mass prediction.
4. Results
We show the preferred parameter space in the ( M , M / ) and ( A , tan β ) plane in Fig. 2. Afterincluding the Higgs boson observables, we find that the stau co-annihilation region is preferredover the focus-point region since it can accommodate slightly better the correct Higgs boson mass.Governed by the observed Higgs boson mass, the fit prefers large negative values of A , where thestop mixing is (nearly) maximized, and a vanishing trilinear coupling, A = β is rather unconstrained. We find a weak linear correlation between tan β and M for the preferred parameter space. The best-fit point is found at M =
504 GeV , M / = , A = − , tan β = , m t = .
74 GeV , (4.1)with a good fit quality of χ / ndf = . /
59. The tension observed in previous studies [2] betweenthe direct LHC limits and the Higgs boson mass measurements (preferring a heavy colored SUSYparticle spectrum) versus the anomalous magnetic moment of the muon (preferring a light uncol-ored SUSY particle spectrum) becomes more severe with the updated LHC exclusion limit. In theremaining viable CMSSM parameter space the SUSY contribution to the anomalous magnetic mo-ment is rather negligible. This tension motivates the consideration of more general models whichabrogate the strong connection of the colored and uncolored sparticle masses of the CMSSM.With the inclusion of the Higgs boson signal rate measurements — which generally are in verygood agreement with the predictions for the SM Higgs boson — the fit quality improves signifi-cantly. This is because the CMSSM parameter regions, which are not excluded by the LHC SUSYsearches, naturally feature a decoupled heavy Higgs spectrum and therefore a SM-like lightestHiggs boson. 4 onstrained SUSY after the Higgs Boson Discovery
Tim Stefaniak ) SM ZZ → h Γ / CMSSM ZZ → h Γ ) / ( SM Γ / CMSSM Γ (0.9 0.95 1 1.05 1.1 1.15 1.2 1.25 Best fit point σ
1D 1 σ
1D 2 WW → h γγ → h ττ → h bb → h gg → h all observables = 3 GeV htheo m ∆ (a) (SM normalized) partial decay widths of the lightestHiggs boson, normalized to the h → ZZ decay mode. h A H + H χ χ χ χ χ χ R l~ L l~ τ∼ τ∼ R q~ L q~ b~ b~ t~ t~ g~ P a r t i c l e M a ss ( G e V ) Environment σ
1 Environment σ (b) Mass spectrum of the Higgs bosons and sparticles. Figure 3:
Predictions for the Higgs and SUSY particle spectrum in the allowed parameter region.
In Fig. 3(a) we show the (SM normalized) Higgs boson partial decay widths, normalized tothe h → ZZ decay mode, for the preferred parameter regions. Deviations from the SM predictionare at most ∼ O ( ) , making the CMSSM extremely difficult to probe via Higgs boson signalrate measurements even at a future linear collider. Note also, that there are remaining theoreticaluncertainties of the F EYN H IGGS calculation of these rates.The sparticle and Higgs boson mass spectrum predicted by the fit is shown in Fig. 3(b). Theheavy Higgs bosons with masses (cid:38) ∼ t may be as light as (cid:38)
750 GeV. The lightest neutralino and stau arenearly mass degenerate with masses around 350 −
600 GeV.
5. Conclusions and outlook
We presented preliminary results from an ongoing F
ITTINO global fit analysis of the CMSSM,including various up-to-date observables and constraints from low-energy and flavor physics, as-trophysics, direct LHC SUSY searches and Higgs boson searches. For the lightest Higgs bosonwe include the Higgs mass and signal rate measurements from the Tevatron and LHC experimentsusing the program H
IGGS S IGNALS . The direct LHC SUSY limit and the Higgs boson mass mea-surements drive the fit to regions with very heavy sparticles, making the CMSSM incapable ofexplaining the discrepancies observed in the anomalous magnetic moment of the muon. The inclu-sion of the Higgs signal rates improves the fit quality, since the already preferred region naturallyfeatures a SM-like light Higgs boson. However, it does not change significantly the general picture.The prospects for a future direct or indirect discovery of the CMSSM are rather grim. Phenomeno-logically, the CMSSM looks like the SM with a viable explanation of dark matter.The next steps within this project comprise a dedicated p -value calculation with repeated fitsto randomly generated pseudo-measurements in order to provide a quantitative statement about thefit quality. Furthermore, we plan to investigate the fit outcome and p -value dependence using amore inclusive Higgs boson signal rate observable set with ∼ O ( ) measurements.5 onstrained SUSY after the Higgs Boson Discovery Tim Stefaniak
Acknowledgments
We thank the organizers of the EPS HEP 2013 conference for the opportunity to present thiswork. TS is grateful for the warm hospitality and support of N
ORDITA during his extended stayin Stockholm and thanks the participants and organizers of the N
ORDITA
Workshop “Beyond theLHC” for a stimulating workshop in a very pleasant atmosphere.
References [1] G. Aad et al. [ATLAS Collaboration], Phys. Lett. B (2012) 1 [arXiv:1207.7214 [hep-ex]];S. Chatrchyan et al. [CMS Collaboration], Phys. Lett. B (2012) 30 [arXiv:1207.7235 [hep-ex]].[2] P. Bechtle, T. Bringmann, K. Desch, H. Dreiner, M. Hamer, C. Hensel, M. Kramer and N. Nguyen etal. , JHEP (2012) 098 [arXiv:1204.4199 [hep-ph]].[3] P. Bechtle, K. Desch and P. Wienemann, Comput. Phys. Commun. (2006) 47 [hep-ph/0412012].[4] P. Bechtle, K. Desch, M. Uhlenbrock and P. Wienemann, Eur. Phys. J. C (2010) 215[arXiv:0907.2589 [hep-ph]].[5] W. Porod, Comput. Phys. Commun. (2003) 275 [hep-ph/0301101]; W. Porod and F. Staub,Comput. Phys. Commun. (2012) 2458 [arXiv:1104.1573 [hep-ph]].[6] S. Heinemeyer, W. Hollik and G. Weiglein, Comput. Phys. Commun. (2000) 76[hep-ph/9812320]; S. Heinemeyer, W. Hollik and G. Weiglein, Eur. Phys. J. C (1999) 343[hep-ph/9812472]; G. Degrassi, S. Heinemeyer, W. Hollik, P. Slavich and G. Weiglein, Eur. Phys. J. C (2003) 133 [hep-ph/0212020].[7] F. Mahmoudi, Comput. Phys. Commun. (2008) 745 [arXiv:0710.2067 [hep-ph]]; F. Mahmoudi,Comput. Phys. Commun. (2009) 1579 [arXiv:0808.3144 [hep-ph]].[8] G. Belanger, F. Boudjema, A. Pukhov and A. Semenov, Comput. Phys. Commun. (2002) 103[hep-ph/0112278]; G. Belanger, F. Boudjema, A. Pukhov and A. Semenov, arXiv:1305.0237 [hep-ph].[9] N. Nguyen, D. Horns and T. Bringmann, arXiv:1202.1385 [astro-ph.HE].[10] P. Gondolo, J. Edsjö, P. Ullio, L. Bergström, M. Schelke and E. A. Baltz, JCAP (2004) 008[astro-ph/0406204].[11] P. Bechtle, O. Brein, S. Heinemeyer, G. Weiglein and K. E. Williams, Comput. Phys. Commun. (2010) 138 [arXiv:0811.4169 [hep-ph]]; P. Bechtle, O. Brein, S. Heinemeyer, G. Weiglein andK. E. Williams, Comput. Phys. Commun. (2011) 2605 [arXiv:1102.1898 [hep-ph]].[12] P. Bechtle, O. Brein, S. Heinemeyer, O. Stål, T. Stefaniak, G. Weiglein and K. Williams, PoSCHARGED (2012) 024 [arXiv:1301.2345 [hep-ph]]. See also new manual available on http://higgsbounds.hepforge.org .[13] P. Bechtle, S. Heinemeyer, O. Stål, T. Stefaniak and G. Weiglein, arXiv:1305.1933 [hep-ph]; O. Ståland T. Stefaniak, in proceedings of The European Physical Society Conference on High EnergyPhysics , PoS(EPS-HEP 2013)314 .[14] R. Aaij et al. [LHCb Collaboration], Phys. Rev. Lett. (2013) 021801 [arXiv:1211.2674 [hep-ex]].[15] I. Adachi et al. [Belle Collaboration], Phys. Rev. Lett. (2013) 131801 [arXiv:1208.4678 [hep-ex]].[16] P. A. R. Ade et al. [Planck Collaboration], arXiv:1303.5076 [astro-ph.CO].[17] The ATLAS Collaboration, ATLAS-CONF-2013-047.[Planck Collaboration], arXiv:1303.5076 [astro-ph.CO].[17] The ATLAS Collaboration, ATLAS-CONF-2013-047.