×
近期发现有不法分子冒充我刊与作者联系,借此进行欺诈等不法行为,请广大作者加以鉴别,如遇诈骗行为,请第一时间与我刊编辑部联系确认(《中国物理C》(英文)编辑部电话:010-88235947,010-88236950),并作报警处理。
本刊再次郑重声明:
(1)本刊官方网址为cpc.ihep.ac.cn和https://iopscience.iop.org/journal/1674-1137
(2)本刊采编系统作者中心是投稿的唯一路径,该系统为ScholarOne远程稿件采编系统,仅在本刊投稿网网址(https://mc03.manuscriptcentral.com/cpc)设有登录入口。本刊不接受其他方式的投稿,如打印稿投稿、E-mail信箱投稿等,若以此种方式接收投稿均为假冒。
(3)所有投稿均需经过严格的同行评议、编辑加工后方可发表,本刊不存在所谓的“编辑部内部征稿”。如果有人以“编辑部内部人员”名义帮助作者发稿,并收取发表费用,均为假冒。
                  
《中国物理C》(英文)编辑部
2024年10月30日

MSSM at future Higgs factories

  • In this work, we study the implications of Higgs precision measurements at future Higgs factories for the MSSM parameter space, focusing on the dominant stop sector contributions. We perform a multi-variable fit to both the signal strength for various Higgs decay channels at Higgs factories and the Higgs mass. The χ2 fit results show sensitivity to mA, tan β, stop mass parameter mSUSY , and the stop left-right mixing parameter Xt. We also study the impact of the Higgs mass prediction on the MSSM and compare the sensitivities of different Higgs factories.
  • 加载中
  • [1] G. Aad et al. (ATLAS Collaboration), Phys. Lett. B 716, 1-29 (2012), arXiv:1207.7214 doi: 10.1016/j.physletb.2012.08.020
    [2] S. Chatrchyan et al. (CMS Collaboration), Phys. Lett. B 716, 30-61 (2012), arXiv:1207.7235 doi: 10.1016/j.physletb.2012.08.021
    [3] M. Dong and G. Li (CEPC Study Group Collaboration), CEPC Conceptual Design Report: Volume 2 - Physics & Detector, arXiv: 1811.10545
    [4] A. Abada et al. (FCC Collaboration), Future Circular Collider: Vol. 1 Physics opportunities
    [5] A. Abada et al. (FCC Collaboration), Future Circular Collider: Vol. 2 The Lepton Collider (FCC-ee)
    [6] P. Bambade et al., The International Linear Collider: A Global Project, arXiv: 1903.01629
    [7] S. Chang, P. J. Fox, and N. Weiner, JHEP 08, 068 (2006), arXiv:hep-ph/0511250
    [8] S. P. Martin, A Supersymmetry primer 21, 1-153 (2010), arXiv:hep-ph/9709356
    [9] N. D. Christensen, T. Han, and S. Su, Phys. Rev. D 85, 115018 (2012), arXiv:1203.3207 doi: 10.1103/PhysRevD.85.115018
    [10] J. Fan, M. Reece, and L.-T. Wang, JHEP 08, 152 (2015), arXiv:1412.3107
    [11] J. Fan and M. Reece, JHEP 06, 031 (2014), arXiv:1401.7671
    [12] A. Drozd, J. Ellis, J. Quevillon et al., JHEP 06, 028 (2015), arXiv:1504.02409
    [13] A. Drozd, J. Ellis, J. Quevillon et al., JHEP 03, 180 (2016), arXiv:1512.03003
    [14] M. Carena, J. R. Espinosa, M. Quiros et al., Phys. Lett. B 355, 209-221 (1995), arXiv:hep-ph/9504316
    [15] S. Heinemeyer, W. Hollik, and G. Weiglein, Comput. Phys. Commun. 124, 76-89 (2000), arXiv:hep-ph/9812320 doi: 10.1016/S0010-4655(99)00364-1
    [16] M. Carena, S. Heinemeyer, C. Wagner et al., Eur. Phys. J. C 26, 601-607 (2003), arXiv:hep-ph/0202167 doi: 10.1140/epjc/s2002-01084-3
    [17] A. Dabelstein, Nucl. Phys. B 456, 25-56 (1995), arXiv:hep-ph/9503443
    [18] H. Bahl, T. Hahn, S. Heinemeyer et al., Comput. Phys. Commun. 249, 107099 (2020), arXiv:1811.09073 doi: 10.1016/j.cpc.2019.107099
    [19] H. Bahl, S. Heinemeyer, W. Hollik et al., Eur. Phys. J. C 78(1), 57 (2018), arXiv:1706.00346 doi: 10.1140/epjc/s10052-018-5544-3
    [20] H. Bahl and W. Hollik, Eur. Phys. J. C 76(9), 499 (2016), arXiv:1608.01880
    [21] T. Hahn, S. Heinemeyer, W. Hollik et al., Phys. Rev. Lett. 112(14), 141801 (2014), arXiv:1312.4937 doi: 10.1103/PhysRevLett.112.141801
    [22] M. Frank, T. Hahn, S. Heinemeyer et al., JHEP 02, 047 (2007), arXiv:hep-ph/0611326
    [23] G. Degrassi, S. Heinemeyer, W. Hollik et al., Eur. Phys. J. C 28, 133-143 (2003), arXiv:hep-ph/0212020
    [24] S. Heinemeyer, W. Hollik, and G. Weiglein, Eur. Phys. J. C 9, 343-366 (1999), arXiv:hep-ph/9812472
    [25] CEPC-SPPC Study Group, CEPC-SPPC Preliminary Conceptual Design Report. 1. Physics and Detector, http://cepc.ihep.ac.cn/preCDR/volume.html, 2015
    [26] A. Abada et al. (FCC Collaboration), Eur. Phys. J. C 79(6), 474 (2019
    [27] A. Abada et al. (FCC Collaboration), Eur. Phys. J. ST 228(2), 261-623 (2019 doi: 10.1140/epjst/e2019-900045-4
    [28] M. Bicer et al. (TLEP Design Study Working Group Collaboration), JHEP 01, 164 (2014), arXiv:1308.6176
    [29] N. Alipour Tehrani et al., FCC-ee: Your Questions Answered, in CERN Council Open Symposium on the Update of European Strategy for Particle Physics (A. Blondel and P. Janot, eds.), 6, 2019. arXiv: 1906.02693
    [30] H. Baer, T. Barklow, K. Fujii et al., The International Linear Collider Technical Design Report - Volume 2: Physics, arXiv: 1306.6352
    [31] K. Fujii et al. (LCC Physics Working Group Collaboration), Tests of the Standard Model at the International Linear Collider, arXiv: 1908.11299
    [32] K. Fujii et al., ILC Study Questions for Snowmass 2021, arXiv: 2007.03650
    [33] N. Chen, T. Han, S. Li et al., Type-I 2HDM under the Higgs and Electroweak Precision Measurements, arXiv: 1912.01431
    [34] H. Bahl, S. Heinemeyer, W. Hollik et al., Eur. Phys. J. C 80(6), 497 (2020), arXiv:1912.04199 doi: 10.1140/epjc/s10052-020-8079-3
    [35] T. Han, S. Li, S. Su et al., Comparative Studies of 2HDMs under the Higgs Boson Precision Measurements, arXiv: 2008.05492
    [36] D. Stöckinger and J. Unger, Nucl. Phys. B 935, 1-16 (2018), arXiv:1804.05619 doi: 10.1016/j.nuclphysb.2018.08.005
    [37] R. Harlander, J. Klappert, A. Ochoa Franco et al., Eur. Phys. J. C 78(10), 874 (2018), arXiv:1807.03509 doi: 10.1140/epjc/s10052-018-6351-6
    [38] R. V. Harlander, J. Klappert, and A. Voigt, Higgs mass prediction in the MSSM at three-loop level in a pure DR context, arXiv: 1708.05720
    [39] B. Henning, X. Lu, and H. Murayama, JHEP 01, 023 (2016), arXiv:1412.1837 doi: 10.26524/jap14
    [40] A. Dabelstein, Z. Phys. C 67, 495-512 (1995), arXiv:hep-ph/9409375 doi: 10.1007/BF01624592
    [41] S. Heinemeyer, W. Hollik, and G. Weiglein, Decay widths of the neutral CP even MSSM Higgs bosons in the Feynman diagrammatic approach, hep-ph/0003022
    [42] J. Gu, H. Li, Z. Liu et al., JHEP 12, 153 (2017), arXiv:1709.06103
    [43] M. Carena, D. Garcia, U. Nierste et al., Nucl. Phys. B 577, 88-120 (2000), arXiv:hep-ph/9912516
    [44] S. Heinemeyer, Int. J. Mod. Phys. A 21, 2659-2772 (2006), arXiv:hep-ph/0407244 doi: 10.1142/S0217751X06031028
    [45] M. Carena, S. Heinemeyer, O. Stl et al., Eur. Phys. J. C 73(9), 2552 (2013), arXiv:1302.7033 doi: 10.1140/epjc/s10052-013-2552-1
    [46] E. Bagnaschi et al., Eur. Phys. J. C 79(7), 617 (2019), arXiv:1808.07542 doi: 10.1140/epjc/s10052-019-7114-8
    [47] A. Djouadi, L. Maiani, A. Polosa et al., JHEP 06, 168 (2015), arXiv:1502.05653
    [48] G. Aad et al. (ATLAS Collaboration), Search for heavy Higgs bosons decaying into two tau leptons with the ATLAS detector using pp collisions at √ s = 13 TeV, arXiv: 2002.12223
    [49] G. Aad et al. (ATLAS Collaboration), Search for heavy neutral Higgs bosons produced in association with b-quarks and decaying to b-quarks at √ s = 13 TeV with the ATLAS detector, arXiv: 1907.02749
    [50] A. M. Sirunyan et al. (CMS Collaboration), JHEP 09, 007 (2018), arXiv:1803.06553
    [51] CMS Collaboration, Search for heavy Higgs bosons decaying to a top quark pair in proton-proton collisions at \begin{document}$\sqrt s = 13\;{\rm{TeV}}$\end{document}
    [52] M. Aaboud et al. (ATLAS Collaboration), Eur. Phys. J. C 78(4), 293 (2018), arXiv:1712.06386
    [53] M. Aaboud et al. (ATLAS Collaboration), JHEP 03, 174 (2018), arXiv: 1712.06518, [Erratum: JHEP 11, 051 (2018)]
    [54] M. Aaboud et al. (ATLAS Collaboration), Eur. Phys. J. C 78(1), 24 (2018), arXiv:1710.01123 doi: 10.1140/epjc/s10052-017-5491-4
    [55] G. Aad et al. (ATLAS Collaboration), Phys. Lett. B 800, 135103 (2020), arXiv:1906.02025 doi: 10.1016/j.physletb.2019.135103
    [56] A. M. Sirunyan et al. (CMS Collaboration), JHEP 03, 034 (2020), arXiv:1912.01594
    [57] A. M. Sirunyan et al. (CMS Collaboration), Phys. Lett. B 778, 101-127 (2018), arXiv:1707.02909 doi: 10.1016/j.physletb.2018.01.001
    [58] A. M. Sirunyan et al. (CMS Collaboration), JHEP 03, 065 (2020), arXiv:1910.11634
    [59] M. Aaboud et al. (ATLAS Collaboration), JHEP 09, 139 (2018), arXiv:1807.07915
    [60] A. M. Sirunyan et al. (CMS Collaboration), JHEP 01, 096 (2020), arXiv:1908.09206
    [61] M. Aaboud et al. (ATLAS Collaboration), JHEP 11, 085 (2018), arXiv:1808.03599
    [62] G. Aad et al. (ATLAS Collaboration), Eur. Phys. J. C 80(8), 737 (2020), arXiv:2004.14060 doi: 10.1140/epjc/s10052-020-8102-8
    [63] A. M. Sirunyan et al. (CMS Collaboration), JHEP 05, 032 (2020), arXiv:1912.08887
    [64] A. M. Sirunyan et al. (CMS Collaboration), arXiv: 2008.05936
  • 加载中

Figures(7)

Get Citation
Honglei Li, Huayang Song, Shufang Su, Wei Su and Jin Min Yang. MSSM at future Higgs factories[J]. Chinese Physics C. doi: 10.1088/1674-1137/abe19b
Honglei Li, Huayang Song, Shufang Su, Wei Su and Jin Min Yang. MSSM at future Higgs factories[J]. Chinese Physics C.  doi: 10.1088/1674-1137/abe19b shu
Milestone
Received: 2020-11-22
Article Metric

Article Views(1190)
PDF Downloads(37)
Cited by(0)
Policy on re-use
To reuse of Open Access content published by CPC, for content published under the terms of the Creative Commons Attribution 3.0 license (“CC CY”), the users don’t need to request permission to copy, distribute and display the final published version of the article and to create derivative works, subject to appropriate attribution.
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Email This Article

Title:
Email:

MSSM at future Higgs factories

    Corresponding author: Wei Su, wei.su@adelaide.edu.au, Corresponding author
  • 1. School of Physics and Technology, University of Jinan, Jinan, Shandong 250022, China
  • 2. Department of Physics, University of Arizona, Tucson, Arizona 85721, USA
  • 3. ARC Centre of Excellence for Dark Matter Particle Physics and CSSM, Department of Physics, University of Adelaide, SA 5005, Australia
  • 4. CAS Key Laboratory of Theoretical Physics, Institute of Theoretical Physics, Chinese Academy of Sciences, Beijing 100190, China
  • 5. School of Physical Sciences, University of Chinese Academy of Sciences, Beijing 100049, China

Abstract: In this work, we study the implications of Higgs precision measurements at future Higgs factories for the MSSM parameter space, focusing on the dominant stop sector contributions. We perform a multi-variable fit to both the signal strength for various Higgs decay channels at Higgs factories and the Higgs mass. The χ2 fit results show sensitivity to mA, tan β, stop mass parameter mSUSY , and the stop left-right mixing parameter Xt. We also study the impact of the Higgs mass prediction on the MSSM and compare the sensitivities of different Higgs factories.

    HTML

    I.   INTRODUCTION
    • After the discovery of a 125 GeV SM-like Higgs boson at the LHC [1,2], several proposals for a Higgs factory have been considered for precise measurements of the Higgs mass and couplings, including the Circular Electron Positron Collider (CEPC) in China [3], the electron-positron stage of the Future Circular Collider (FCC-ee) at CERN [4,5], and the International Linear Collider (ILC) in Japan [6]. With about $ 10^6 $ Higgs produced, the Higgs mass can be measured to an accuracy of about 10 MeV. Various Higgs couplings to the Standard Model (SM) particles can be measured at about 1% level, while the $ hZZ $ coupling can be measured at about 0.2%. If no deviation from the SM predictions is observed at future Higgs factories, severe constraints can be imposed on the parameter space of new physics models. On the other hand, if certain deviations are observed, it will provide great insights into new physics models beyond the SM, in particular, on the masses and couplings of new particles, which can be cross checked with direct searches for new particles at future high energy colliders.

      In this work we focus on the implications of Higgs precision measurements for the Minimal Supersymmetric Standard Model (MSSM), which is one of the well-motivated new physics models proposed to solve the hierarchy problem [7,8]. The Higgs sector of the MSSM is the same as that of the Type-II two Higgs doublet model (2HDM), with one Higgs doublet $ H_u $ coupling to up type quarks, and the other Higgs doublet $ H_d $ coupling to down type quarks and charged leptons. After electroweak symmetry breaking, there are five physical Higgses: two CP-even Higgses $ h $ and $ H $, one CP-odd Higgs $ A $ and a pair of charged ones $ H^{\pm} $. In our analyses below, we take the light CP-even Higgs $ h $ to be the observed 125 GeV SM-like Higgs. The other possibility of the heavy CP-even Higgs $ H $ being the 125 GeV Higgs is tightly constrained by the existing experimental searches, as shown in Ref. [9]. At tree level, the masses of MSSM Higgses are completely determined by only two parameters: the CP-odd Higgs mass $ m_A $ and the ratio of two Higgs vacuum expectation values $ \tan\beta $. The mass of the SM-like light CP-even Higgs, however, receives large radiative corrections, with the dominant contribution from the supersymmetric top partners. The Higgs couplings to the SM sector also receive radiative corrections, characterised by the effective mixing angle $\alpha_{\rm eff}$. In addition, the bottom Yukawa coupling receives large vertex corrections. Furthermore, Higgs couplings to a pair of photons or gluons receive loop contributions from the stop sector, which are at about the same order as the SM contributions.

      While in general there are more than 100 parameters in the MSSM, when we focus on the Higgs sector and the dominant loop contributions from the stop sector, only four are the most relevant: $ \tan\beta $, $ m_A $, the soft supersymmetry (SUSY) breaking stop mass parameter $ m_{\rm SUSY} $, and the left-right stop mixing parameter $ X_t $. Other parameters, such as the mass parameters for the sbottom and gluinos, could enter as well. Those effects become important in a particular corner of the parameter space, which is left for future dedicated studies.

      To study the implications of Higgs precision measurements for the parameter space of the MSSM, we perform a multi-variable $ \chi^2 $ fit to both the signal strength $ \mu $ for various Higgs decay channels and the Higgs mass. Earlier works on the implications of Higgs precision measurements mostly focused on the loop-induced channels $ h \rightarrow gg $ and $ h \rightarrow \gamma\gamma $ [10-13], given that both the SM and the MSSM contributions enter at the same order. In our work, we include all the Higgs decay channels measured at Higgs factories, as well as the Higgs mass. The MSSM predictions of those quantities have been studied extensively in the literature [14-16]. For the MSSM corrections to the Higgs couplings to the SM particles, we adopt the $ \alpha_{\rm eff} $ method [14,17]. We also include the additional vertex corrections to the bottom Yukawa and loop-induced couplings of $ hgg $ and $ h\gamma\gamma $. We use the state-of-art program FeynHiggs [15, 18-24] to obtain $ \alpha_{\rm eff} $ and $ m_h $ in the framework of the MSSM.

      In Sec. II, we briefly summarize the Higgs precision measurements at various Higgs factories. We also introduce the $ \chi^2 $ fit formalism used in our analyses. In Sec. III, we discuss the MSSM Higgs sector and stop sector that are needed in our analyses, as well as the SM-like Higgs couplings in the MSSM. In Sec. IV, we summarize the current direct search limits on the mass of the CP-odd Higgs and the stop sector. In Sec. V, we perform detailed analyses of various contributions to the total $ \chi^2 $. In Sec. VI, we present the 95% C.L. allowed region of the MSSM parameter space under the CEPC precisions. In Sec. VII, we compare the reach of the CEPC, the FCC-ee and the ILC. We reserve Sec. VIII for conclusions.

    II.   HIGGS PRECISION MEASUREMENTS AND $ \chi^2 $ FIT
    • Analyses of precision measurements of Higgs decay channels have been performed for the CEPC [3,25], the FCC-ee [26-29], and the ILC [6,30-32] in recent years. A summary of the most up-to-date results on $ \Delta(\sigma\times{\rm Br})/(\sigma\times{\rm Br}) $, as well as the total production cross section $ \Delta \sigma/\sigma $, can be found in Table 3 in Ref. [33], which will be used in our current study. The dominant production channel at 240$ - $250 GeV is associated $ Zh $ production, with the best measured channel being $ h\rightarrow b\bar{b} $, given its large decay branching fraction. A precision of about 0.3% can be achieved for this channel. The precisions for $ h\rightarrow gg,\ WW^*, \tau^+\tau^- $ are about 1%, while $ h\rightarrow c\bar{c} $ is about 2%$ - $3%. The precisions for $ h\rightarrow ZZ^*, \gamma\gamma $ are worse, about 5%$ - $7%, given its suppressed decay branching fractions. The sensitivities for the three Higgs factories are very similar. The weak boson fusion (WBF) process $ e^+e^-\rightarrow \nu \bar\nu h $ becomes more important at higher center of mass energy, with a precision of about 0.23% for the $ h\rightarrow b\bar{b} $ channel at the ILC 500 GeV with 4 $ {\rm ab}^{-1} $ integrated luminosity [6,31].

      To analyze the implications of Higgs precision measurements for the MSSM parameter space, we perform a multi-variable $ \chi^2 $ fit

      $ \begin{aligned}[b] \chi^2_{\rm total} =& \chi^2_{m_h} + \chi^2_{\mu} = \frac{(m_h^{\rm MSSM}-m_h^{\rm obs})^2}{(\Delta m_h)^2}\\& + \sum_{i = f,V..} \frac{(\mu_i^{\rm MSSM}-\mu_i^{\rm obs})^2}{(\Delta \mu_i)^2}, \end{aligned} $

      (1)

      in which $ \mu_i^{\rm{MSSM}} = (\sigma\times \rm{Br}_i)_{\rm{MSSM}}/(\sigma\times \rm{Br}_i)_{\rm{SM}} $ is the signal strength for various Higgs search channels. Here $ \chi^2_{m_h} $ and $ \chi^2_{\mu} $ refer to contributions to the overall $\chi^2_{\rm total}$ from the Higgs mass and signal strength measurements, respectively. For $ \chi^2_{m_h} $, given the small experimental uncertainties and the relatively large theoretical uncertainties in determining $ m_h $ in the MSSM, we set $ \Delta m_h $ to be 3 GeV, taking into account uncertainties coming from higher order radiative corrections [20-23], as well as propagating uncertainties from SM input parameters like $ m_t $. Results with smaller $ \Delta m_h = 1 $ GeV and 2 GeV are also presented in Sec. VI, to show the impact of possible future improvement in $ m_h $ calculation including higher order corrections [34]. For $ \chi^2_\mu $, $ \Delta \mu_i $ is the experimental expected precision in determining the signal strength for a particular Higgs decay channel.

      For future Higgs factories, $ \mu_i^{\rm{obs}} $ are set to be unity in our analyses, assuming no deviations from the SM predictions are observed. If deviations are observed in the future, we can use the same $ \chi^2 $ fit method to determine the constrained parameter space, with $ \mu_i^{\rm{obs}} $ being the observed experimental central value [35]. Usually, the correlations among different search channels at Higgs factories are not provided and are thus assumed to be zero.

      In our analyses, we determine the allowed parameter region at 95% confidence level (C.L.) by a multi-variable fit to the Higgs decay signal strengths of various channels and Higgs masses. For the one-, two- or three-parameter fit, the corresponding $ \Delta\chi^2 = \chi^2-\chi_{\rm{min}}^2 $ at 95% C.L. is 3.84, 5.99 or 7.82, respectively. Note that when we present the results of our three-parameter fit in Sec. VI, we project the three-dimensional space onto a two-dimensional plane for several benchmark points in the third dimension of the parameter space. Most of the results presented below are for the CEPC precisions. We compare the reaches of the CEPC, the FCC-ee and the ILC in Sec. VII.

    III.   HIGGS AND STOP SECTOR OF THE MSSM

      A.   Higgs mass in the MSSM

    • In our analyses, we identify the light CP-even Higgs $ h $ in the MSSM as the observed 125 GeV SM-like Higgs. Its mass receives large radiative corrections, dominantly from the stop sector, as well as the sbottom sector at large $ \tan\beta $. There have been extensive studies of the MSSM loop correction to the Higgs masses up to next-to-next order [14-16], which includes full one-loop contributions as well as the leading two-loop contributions $ \mathcal{O}\left(\alpha_{t} \alpha_{s}, \alpha_{b} \alpha_{s}, \alpha_{t}^{2}, \alpha_{t} \alpha_{b}, \alpha_{b}^{2}\right) $ to the Higgs two-point functions. There are also works considering the three-loop effects at order $ \mathcal{O}(\alpha_{t, b} \alpha_{s}^{2}, \alpha_{t, b}^{2} \alpha_{s}, \alpha_{t, b}^{3}) $ [36], as well as approximate evaluation at order $ \mathcal{O}(\alpha_{t}^2 \alpha_{s}^{2}) $ [37].

      The CP-even Higgs mass matrix is given by

      $ \begin{aligned}[b] \cal{M}_{\rm Higgs} \rm =& \dfrac{ \sin 2 \beta }{2}\left( \begin{array}{*{20}{c}} \cot \beta \ m_Z^2 + \tan \beta \ m_A^2 & - m_Z^2 - m_A^2\\ - m_Z^2 - m_A^2 & \tan \beta \ m_Z^2 + \cot \beta \ m_A^2 \end{array} \right) \\&+ \left( \begin{array}{ll} \Delta_{11} & \Delta_{12} \\ \Delta_{12} & \Delta_{22} \end{array} \right), \end{aligned} $

      (2)

      with the first term being the tree-level contributions, and $ \Delta_{11},\Delta_{12},\Delta_{22} $ in the second term being the loop-induced Higgs mass corrections [14,17,38]. The masses for the CP-even Higgses are obtained by the diagonalization of the mass matrix:

      $ \begin{aligned}[b] M^2_{H,h,\,\rm eff} =& \frac{m_A^2 + m_Z^2 + \Delta_{22} + \Delta_{11}}{2} \pm \Bigg( \frac{ (m_A^2 + m_Z^2)^2 + ( \Delta_{22} - \Delta_{11} )^2}{4}\\& - m_A^2 m_Z^2 \cos^2 2\beta +\frac{( \Delta_{22}- \Delta_{11}) \cos 2\beta}{2} (m_A^2 - m_Z^2) \\&- \frac{\Delta_{12} \sin 2\beta}{2} (m_A^2 + m_Z^2) + \frac{\Delta_{12}^2}{4} \Bigg)^{1/2}. \\[-5pt] \end{aligned} $

      (3)

      The effective mixing angle $\alpha_{\rm eff}$ between CP-even scalars is defined by

      $ \begin{array}{l} \left(\begin{matrix}h\\H\end{matrix}\right) = \left(\begin{matrix} \cos\alpha_{\rm eff} \sin\alpha_{\rm eff}\\ -\sin\alpha_{\rm eff} \cos\alpha_{\rm eff} \end{matrix}\right)\left(\begin{matrix} \text{Re}H_u^0-v_u \\ \text{Re}H_d^0-v_d \end{matrix}\right), \end{array} $

      (4)

      which takes the form of

      $ \tan \alpha_{\rm eff} = \frac{-(m_A^2+m_Z^2)\sin\beta \cos \beta + \Delta_{12}}{m_Z^2 \cos^2 \beta +m_A^2 \sin^2 \beta +\Delta_{11}-m_{h^0,\rm eff}^2}. $

      (5)

      Out of all the supersymmetric particles, the stop sector gives the dominant loop contributions to the Higgs sector. The stop mass matrix depends on the $ H_u -H_d $ mixing parameter $ \mu $ and soft SUSY breaking parameters $ m_{\tilde{Q}} $, $ m_{\tilde{t}_R} $, and trilinear coupling $ A_t $:

      $ \mathcal{M}_{\rm \tilde{t}}^{\rm 2} \rm = \left( \begin{array}{*{20}{c}} m_{\tilde{Q}}^2 + m_t^2 + m_Z^2 \left(\dfrac{1}{2} - \dfrac{2}{3} s_W^2\right) \cos 2 \beta \quad m_t (A_t- \mu \cot \beta ) \\ m_t (A_t - \mu \cot \beta ) \quad m_{\tilde t_R}^2 + m_t^2 +\dfrac{2}{3} m_Z^2 s_W^2 \cos 2 \beta \end{array} \right). $

      (6)

      The stop left-right mixing parameter is defined as $ X_t \equiv A_t - \mu \cot \beta $, which enters the off-diagonal term, and plays an important role in the radiative corrections to the Higgs mass. For our analyses below, we assume mass degeneracy of left- and right-handed top squarks and take the most relevant model parameters as:

      $ \begin{array}{l} \tan \beta, m_A, m_{\rm SUSY}\equiv m_{\tilde Q} = m_{\tilde t_R}, X_t. \end{array} $

      (7)
    • B.   Higgs couplings with $\alpha_{\rm eff}$ method

    • The effective Lagrangian of the Higgs couplings to a pair of fermions and gauge bosons can be written as [39]

      $ \begin{aligned}[b] \mathcal{L} = & \kappa_Z \frac{m_Z^2}{v}Z_{\mu}Z^{\mu}h+\kappa_W \frac{2 m_W^2}{v}W_{\mu}^+ W^{\mu-}h + \kappa_g \frac{\alpha_s}{12 \pi v} G^a_{\mu\nu}G^{a\mu\nu}h \\&+ \kappa_{\gamma} \frac{\alpha}{2\pi v} A_{\mu\nu}A^{\mu\nu}h -\left( \kappa_t \sum_{f = u,c,t} \frac{m_f}{v}f \bar f \right.\\& \left.+ \kappa_b \sum_{f = d,s,b} \frac{m_f}{v}f \bar f+ \kappa_{\tau} \sum_{f = e,\mu,\tau} \frac{m_f}{v}f \bar f \right)h \end{aligned} $

      (8)

      with $\kappa_i = {g_{hii}^{\rm BSM}}/{g_{hii}^{\rm SM}}$ being the Higgs coupling normalized to the SM value. Given that the Yukawa coupling structure of the MSSM is the same as that of the Type-II 2HDM, $ \kappa_u $, $ \kappa_{d,l} $ and $ \kappa_V $ follow the tree-level expression of the Type-II 2HDM, with the mixing angle $ \alpha $ being replaced by the effective mixing angle $\alpha_{\rm eff}$ [14,17], including radiative corrections:

      $ k_u = \frac{\cos \alpha_{\rm eff}}{\sin \beta} ,~ k_{d,l} = -\frac{\sin \alpha_{\rm eff}}{\cos \beta},~ k_V = \sin (\beta-\alpha_{\rm eff}). $

      (9)

      This is the so-called "$\alpha_{\rm eff}$ method" [40], which is used in our analyses to count for the MSSM loop corrections to the SM-like Higgs couplings to the SM particles. This effective method is in good agreement with the full loop results [17,41], under the heavy gluino mass assumption that we adopt in our analyses.

      Given the high experimental precision in the $ h\rightarrow b\bar{b} $ channel, $ \Delta \mu_{b} $ = 0.27% at the CEPC [3,42], and large $ {\rm Br}_{h\rightarrow b\bar{b}} $ = 57.7% , Higgs factories are particularly sensitive to MSSM contributions to $ \kappa_b $. In addition to the loop contributions to $\alpha_{\rm eff}$, which enters $ \kappa_b $ via Eq. (9), additional MSSM loop corrections to $ \kappa_b $ are included in our analyses, which is characterized by $ \Delta m_b $.

      $ \kappa_b = -\frac{\sin \alpha_{\rm eff}}{\cos\beta} \tilde \kappa^{b}_h,\ \ \ \tilde \kappa^{b}_h = \frac{1}{1+\Delta m_b}\left(1-\Delta m_b \frac{1}{\tan \alpha_{\rm eff} \tan \beta}\right). $

      (10)

      Assuming large sbottom and gluino masses, the dominant loop contribution to $ \Delta m_b $ comes from the stop sector [43]:

      $ \Delta m_b^{\rm stop} = \frac{h_t^2}{16 \pi^2}\mu A_t \tan \beta I(m_{\tilde t_1}, m_{\tilde t_2}, \mu). $

      (11)

      The loop-induced Higgs couplings, $ hgg $ and $ h\gamma\gamma $ receive contributions from the SUSY sector as well, which are of the same order as the SM contributions. Therefore, $ hgg $ and $ h\gamma\gamma $ could provide extra sensitivity to the MSSM parameter space [10,12]. In particular, given that the experimental precision for the $ hgg $ channel is about 1% at Higgs factories, this channel is particularly sensitive to stops running in the loop. Contributions from the sbottom sector are typically at least an order of magnitude smaller than those from the stop sector, even in the case of large $ \tan\beta $ [12]. To focus on the dominant effects, we do not include the sbottom effects in our analyses.

      The signal strength $ \mu_i $ that enters the $ \chi^2 $ analyses includes the MSSM contributions to both the Higgs production and decays. We use the state-of-art program FeynHiggs [18,20,34,44] to obtain $ m_h $, $\alpha_{\rm eff}$ and $ \Delta m_b $, calculating the various $ \kappa $s as defined in Eq. (8), which are fed into the evaluation of signal strength $ \mu_i^{\rm MSSM} $.

    IV.   DIRECT SEARCH LIMITS FROM THE LHC
    • As well as the studies of the SM-like Higgs, there have been extensive searches for MSSM heavy Higgses at the LHC. Given the light CP-even Higgs as the observed 125 GeV SM-like Higgs, scenarios such as $m_h^{\rm mod}$ [45], $ M_h^{125} $ [46] and hMSSM [47] are proposed to test the model parameter spaces in the $ m_A-\tan\beta $ plane. Based on the data collected during the LHC Run 2 with an integrated luminosity of 139 $ \,{{\rm{fb}}}^{-1} $ at $ \sqrt{s} = 13 $ TeV, the ATLAS collaboration searched for the heavy neutral Higgs bosons over the mass region 0.2$ - $2.5 TeV with $ A/H \rightarrow \tau^+\tau^- $ decay [48]. In the $ M_h^{125} $ scenario, the data exclude the parameter space of $ \tan\beta>8 $ for $ m_A = 1.0 $ TeV, $ \tan\beta>21 $ for $ m_A = 1.5 $ TeV, and $ \tan\beta>60 $ for $ m_A = 2.0 $ TeV, which are the strongest exclusion limits in the large $ \tan\beta $ region. Exclusion from $ A/H \rightarrow b \bar {b} $ is weaker: for $ \tan\beta $ between $ 20-60 $, $ m_A $ in the mass region of 0.45$ - $0.9 TeV has been excluded with $ bbH/A $ production in the scenarios of hMSSM [49]. Results from CMS are similar [50].

      In the low $ \tan\beta $ region, $ bb $ and $ \tau\tau $ channels are less constraining given the reduced Yukawa couplings. CMS searches with $ A/H \rightarrow t \bar t $ exclude the value of $ m_{A} $ at 400 (700 GeV) for $ \tan\beta $ below 1.5 (1.0) [51]. Decay modes of $ H\rightarrow ZZ $, $ A\rightarrow hZ $, $ H\rightarrow WW $, and $ H\rightarrow hh $ also constrain the parameter space in the low $ \tan\beta $ region. Combining the results from these channels, the mass region of 200$ - $600 GeV is excluded with $ \tan\beta $ value between 1 to 6 at both the ATLAS and CMS experiments [52-57]. In addition, CMS searches of $ A\rightarrow hZ \rightarrow \tau\tau\ell\ell $ exclude $ \tan\beta $ values below 1.6 at $ m_A = 220 $ GeV and 3.7 at $ m_A = 300 $ GeV [58] in the hMSSM scenario.

      Searches for charged Higgses produced either in the top quark decay (for $ m_H^\pm < m_t $) or associated with a top quark (for $ m_H^\pm > m_t $), with the subsequent decay of $ H^+\rightarrow \tau \nu $, are performed in the context of hMSSM at the LHC. For a light charged Higgs with $ m_{H^\pm}<m_t $, $ m_{H+}\leqslant 160 $ GeV is excluded [59]. For a heavy charged Higgs, the region of $ \tan\beta = 20 - 60 $ is excluded with $ m_{H^\pm} $ from 200 to 1100 GeV [59]. $ H^+\rightarrow t b $ decay mode is sensitive to the low $ \tan\beta $ region. Value of $ \tan\beta = 1.5 - 0.4 $ are excluded in the $ m_{H^+} $ range of 200 GeV to 1.5 TeV in the context of $m_h^{\rm mod}$ scenario [60,61].

      For the stop sector, the limits are more complicated, given their dependence on the mass spectrum of charginos and neutralinos, as well as the corresponding decay branching fractions. Several channels of the stop decay to the lighter superparticles have been explored. For $ \tilde{t}_1 \rightarrow t \tilde{\chi} _1^0 /bW\tilde{\chi} _1^0/bff^\prime\tilde{\chi} _1^0 $, the latest results show that the stop mass region of $ m_{\tilde{t}_1}<1.2 $ TeV is excluded for $ \tilde{\chi} _1^0 $ mass below about 500 GeV [62,63]. With a light slepton, the $ \tilde{t}_1 \rightarrow b \tilde{\chi} _1^{+} \rightarrow b\nu \tilde{\ell} \rightarrow b \nu \ell \tilde{\chi} _1^0 $ decay channel can exclude masses up to about 1.4 TeV for $ \tilde{t}_1 $ and 900 GeV for $ \tilde{\chi} _1^0 $ with $ m_{\tilde{\chi} _1^{0}}<m_{\tilde{\ell}}<m_{\tilde{\chi} _1^{+}} $ [64].

    V.   CONTRIBUTIONS TO $ \chi^2 $
    • In this section, we choose several typical sets of MSSM model parameters to study the various contributions to $\chi^2_{\rm total}$:

      $ \begin{aligned}[b] & m_A = 1000/2000 \text{ GeV}, \mu = 500 \text{ GeV}, \tan \beta = 3/30,\\& X_t \in (-5000, 5000) \text{ GeV}, m_{\rm SUSY} \in (200, 3000) \text{ GeV}. \end{aligned} $

      (12)

      To identify the stop contributions, we decouple the masses of other sfermions and gluinos. We scan over the parameter space of $ X_t $ and $ m_{\rm SUSY} $ to explore the various contributions to the overall $\chi^2_{\rm total}$ in Fig. 1, with the colored area being the 95% C.L. allowed region, corresponding to $ \Delta\chi^2 = \chi^2-\chi^2_{\rm min} = 5.99 $ for the two-parameter fit. Different color bands correspond to the $ \chi^2 $ value. The four columns in Fig. 1 are $ \chi^2_{m_h} $ representing the contribution from the Higgs mass, $ \chi^2_{gg+\gamma\gamma} $ representing contributions from the loop-induced processes $ h\rightarrow gg $ and $ h\rightarrow\gamma\gamma $, $ \chi^2_{\mu^\prime} $ representing contributions from tree-level Higgs decays to SM fermions and vector bosons, and $\chi^2_{\rm total}$, from left to right. The three rows are for $ (m_A,\tan\beta) = (1\ {\rm TeV}, 30),\ (2\ {\rm TeV}, 30) $ and $ (2\ {\rm TeV}, 3) $, from top to bottom.

      Figure 1.  (color online) 95% C.L. allowed region in the plane of $m_{\rm SUSY} - X_t$ under CEPC precisions for $(m_A,\tan\beta)=(1\ {\rm TeV}, 30)$ (upper panels), (2 TeV, 30) (middle panels) and (2 TeV, 3) (lower panels). For each row, the panels from left to right show $\chi^2_{m_h}$, $\chi^2_{gg+\gamma\gamma}$, $\chi^2_{\mu^\prime}$, and $\chi^2_{\rm total}$. Different color bands correspond to the $\chi^2$ value. See text for details.

      For $ m_A = 1 $ TeV, $ \tan \beta = 30 $, $ \chi^2_{m_h} $ (top left panel) could push $ m_{\rm SUSY} $ to be at least 1.4 TeV with $ X_t = 0 $, and larger $ m_{\rm SUSY} $ is more preferred for such a zero-mixing case. Two max-mixing branches of $| {X_t}/{m_{\rm SUSY}}|\sim 2$ also appear, given that the radiative correction to the Higgs mass is the largest under such a stop maximal mixing region. Larger values of $ m_{\rm SUSY} $ are disfavored in these two branches. $ \chi^2_{m_h} $ behaviour for $ m_A = 2 $ TeV, $ \tan \beta = 30 $ (middle left panel) is very similar to that of $ m_A = 1 $ TeV, given that the tree-level contributions to $ m_h $ are very similar for large $ m_A $ at large $ \tan\beta $. The lower left panel of $ m_A = 2 $ TeV, $ \tan \beta = 3 $, however, shows very different behaviour: the zero mixing region is completely gone and the max-mixing case is preferred with $ m_{\rm SUSY}\gtrsim 1.2 $ TeV, given the need for large radiative corrections with the reduced tree-level value of $ m_h $.

      For loop-induced contributions, $ \chi^2_{gg+\gamma\gamma} $ (second column), $ m_{\rm SUSY}\leqslant 1 $ TeV for the zero-mixing case of $ X_t = 0 $ GeV, as well as $| {X_t}/{m_{\rm SUSY}}| > 3$ are excluded, which corresponds to too-large radiative corrections to $ h\rightarrow gg,\ \gamma\gamma $. There is, however, a large parameter space which remains viable in the $ m_{\rm SUSY} $ vs. $ X_t $ plane. The dependence of $ \chi^2_{gg+\gamma\gamma} $ on $ m_A $ and $ \tan\beta $ is rather weak.

      There are strong constraints coming from the precision measurements of Higgs Yukawa and gauge couplings, as shown in $ \chi^2_{\mu^\prime} $ plots in the third column. The most constraining channel is $ h\rightarrow b\bar{b} $. As a result, for $ m_A = 1 $ TeV, $ \tan \beta = 30 $, only the positive branch of $ X_t $ survives, as shown in the third panel of the top row. The total $\chi^2_{\rm total}$ including all the contributions gives an even more restricted region of $ m_{\rm SUSY}\ge 1.2 $ TeV and $ X_t/m_{\rm SUSY}\sim 2.6 $. Sensitivity to $ \kappa_b $ is reduced for larger values of $ m_A $. For $ m_A = 2 $ TeV, $ \tan \beta = 30 $, there is a larger allowed parameter region when combining all three $ \chi^2 $s together. For $ m_A = 2 $ TeV, with small $ \tan \beta = 3 $ (bottom row), while the sensitivity to the Higgs precision measurements is similar to that of the large $ \tan\beta $ case, stronger constraints from the Higgs mass lead to the final surviving region to be $ m_{\rm SUSY}>1.5 $ TeV, $| {X_t}/{m_{\rm SUSY}}| \approx 2$ [44].

    VI.   MULTI-VARIABLE $ \chi^2 $ FIT RESULTS
    • In this section, we explore the 95% C.L. allowed region with the Higgs precision measurements at the CEPC in various MSSM parameter spaces. With the four most relevant MSSM parameters ($ m_A $, $ \tan\beta $, $ m_{\rm SUSY} $, $ X_t $), we scan in the range:

      $ \begin{aligned}[b] & m_A\in (200, 3000) \text{ GeV}, \tan \beta\in (1, 50),\\& X_t \in (-5000, 5000) \text{ GeV}, m_{\rm SUSY} \in (200, 3000) \text{ GeV}, \end{aligned} $

      (13)

      with $ \mu = 500 \text{ GeV} $. The fitting results vary little when $ \mu $ varies. For the 3D fit performed in our analyses, we fix one variable to a set of benchmark values. When presenting results in the 2D parameter space, we project the 3D results onto the 2D space for a given set of values of the third parameter.

      In Fig. 2, we show the 95% C.L. allowed region in the $ X_t - m_{\rm SUSY} $ plane for various values of $ \tan\beta $ with the CEPC precisions. The left, middle and right panels are for $ m_A = 1,\ 1.5,\ 2 $ TeV, respectively. The low $ \tan \beta $ case receives strong constraints from the Higgs mass precision, especially for smaller values of $ m_A $, as explained in the last section. For $ m_A = 1 $ TeV (left panel), $ \tan \beta \leqslant 25 $ is completely excluded. The surviving region is around the stop max-mixing section of $ |X_t|\approx 2 m_{\rm SUSY} $. Only $ X_t>0 $ branch survives given the $ \kappa_h^b $ effects, as explained in the last section. For $ m_A = 1.5 $ TeV (middle panel), $ \tan \beta < 10 $ is excluded. For $ \tan\beta = 10 $, a small slice of $ X_t<0 $ survives combining all three contributions to $\chi^2_{\rm total}$. Larger regions open up for larger values of $ \tan\beta $. For $ m_A = 2 $ TeV (right panel), $ \tan \beta $ as small as 3 is still allowed. Precision constraints from both the mass and the couplings are relaxed for larger $ \tan\beta $ and larger $ m_A $, resulting in large survival parameter spaces in $ X_t $ vs. $ m_{\rm SUSY} $.

      Figure 2.  (color online) 95% C.L. allowed region in $X_t$ vs. $m_{\rm SUSY}$ plane with CEPC precisions for $m_A=1$ TeV (left panel), 1.5 TeV (middle panel) and 2 TeV (right panel). For each panel, different colored regions correspond to different values of $\tan\beta$.

      In Fig. 3, we show the 95% C.L. allowed region in the $ m_{\rm SUSY} $ vs. $ m_A $ plane with the CEPC precisions for $ X_t = 0 $ (zero-mixing, left panel) and $ X_t = 2 m_{\rm SUSY} $ (max-mixing, right panel). Regions to the right of the curve are the 95% C.L. allowed regions for different values of $ \tan\beta $. For the zero-mixing case, small $ \tan\beta $ receives the strongest constraints, with $ \tan \beta \leqslant 4 $ excluded totally. For the $ \tan \beta = 50 $, $ m_A\geqslant 1350 $ GeV and $ m_{\rm SUSY}\geqslant 850 $ GeV are still allowed. Note that for the zero-mixing case, the most important constraints for $ m_A $ come from the Higgs gauge and Yukawa couplings, while the most important constraint for $ m_{\rm SUSY} $ comes from the Higgs mass precision, except for the large $ \tan \beta $ case, when the Higgs gauge and Yukawa couplings enter as well.

      Figure 3.  (color online) 95% C.L. allowed region in the $m_{\rm SUSY}$ vs. $m_A$ plane with CEPC precisions for $X_t =0$ (zero-mixing, left panel) and $X_t = 2 m_{\rm SUSY}$ (max-mixing, right panel). For each panel, different colored curves correspond to different values of $\tan\beta$, with the region to the right of the curve allowed.

      In the max-mixing case, values of $ \tan\beta $ as low as 3 could be accommodated. The allowed region is typically larger than that of the zero-mixing case. For small $ \tan\beta $, the strongest constraints for $ m_{\rm SUSY} $ are the Higgs mass precision and loop-induced $ hgg $ and $ h\gamma\gamma $. For $ \tan \beta \geqslant 7 $, the lower limit on $ m_{\rm SUSY} $ mostly comes from the loop-induced $ hgg $ and $ h\gamma\gamma $, which are less sensitive to values of $ \tan\beta $. There are, however, upper limits on $ m_{\rm SUSY} $ from too-large contributions to $ m_h $. Limits on $ m_A $ are mostly determined by the precisions of the Higgs couplings.

      The results of the three-parameter fit for $ m_A, m_{\rm SUSY} $ and $ \tan\beta $ are projected onto the $ m_A $ vs. $ \tan\beta $ plane in Fig. 4. Regions above the curve are the 95% C.L. surviving regions with CEPC precisions for different values of $ m_{\rm SUSY} $. In general, $ m_{\rm SUSY}<900 $ GeV is excluded for both the no-mixing and max-mixing cases. For the no-mixing case, when $ m_{\rm SUSY}<1 $ TeV, $ \tan\beta<40 $ is excluded. Limits on $ \tan\beta $ get lower for larger values of $ m_{\rm SUSY} $, which is sensitive in particular for $ 1\ {\rm TeV}<m_{\rm SUSY}<1.5 $ TeV. For the max-mixing case, limits on $ \tan\beta $ are much lower for $ m_{\rm SUSY} = 1 $ TeV. Those features are mainly due to the Higgs mass constraint. The LHC Run-II direct search limits based on $ A/H\rightarrow\tau\tau $ [48] are shown in the grey shaded region, which is complementary to the indirect limits from Higgs precision measurements.

      Figure 4.  (color online) 95% C.L. allowed region in the $\tan\beta$ vs. $m_A$ plane with the CEPC precisions for $X_t =0$ (zero mixing, left panel) and $X_t = 2 m_{\rm SUSY}$ (max-mixing, right panel). For each panel, different colored curves correspond to different values of $m_{\rm SUSY}$, with the region above the curve allowed. The LHC Run-II direct search limits based on $A/H\to\tau\tau$ [48] are shown in the grey shaded region.

      The results of the three-parameter fit for $ m_A, m_{\rm SUSY} $ and $ \tan\beta $ are projected onto the $ m_{\rm SUSY} $ vs. $ \tan\beta $ plane in Fig. 5. For each panel, different colored curves correspond to different values of $ m_A $, with the region above the curve allowed, except for the $ m_A = 2 $ TeV (blue curves) in the right panel, in which the region between the two curves is allowed. In general, $ m_A<1 $ TeV is excluded for both the zero-mixing and max-mixing cases. The lower limits on $ \tan\beta $ are relaxed for larger values of $ m_A $, and is sensitive to the values of $ m_A $ for $ 1.5\ {\rm TeV}<m_A<2 $ TeV. For the max-mixing case and a given $ m_A $, there are upper limits for $ m_{\rm SUSY} $ at large $ \tan\beta $, as shown in the right panel of Fig. 5. This is due to the too-large contribution to $ m_h $ for larger values of $ m_{\rm SUSY} $. For $ m_A = 3 $ TeV, the upper limit for $ m_{\rm SUSY} $ is larger than 3 TeV, and is therefore not shown in the plot.

      Figure 5.  (color online) 95% C.L. allowed region in the $\tan\beta$ vs. $m_{\rm SUSY}$ plane with the CEPC precisions for $X_t =0$ (zero-mixing, left panel) and $X_t = 2 m_{\rm SUSY}$ (max-mixing, right panel). For each panel, different colored curves correspond to different values of $m_A$, with the region above the curve allowed, except for $m_A=2$ TeV (blue curves) in the right panel, in which the region between the two curves is allowed.

      To illustrate the potential impact of future improvement in the MSSM prediction of $ m_h $, in Fig. 6, we show the 95% C.L. allowed region in $ m_A $ vs. $m_{\rm SUSY}$ plane for $ \Delta m_h = 3 $ GeV (solid curve), 2 GeV (dashed curve), and 1 GeV (dotted curve). The lower limit on $ m_{\rm SUSY} $ for the zero-mixing case, and the upper limit on $ m_{\rm SUSY} $ for the max-mixing case depend sensitively on the values of $ \Delta m_h $. Therefore, it is crucial to improve the precision in the $ m_h $ calculation in the MSSM, which will allow us to obtain tight constraints on the SUSY mass scale, in particular on the stop sector, once Higgs precision measurements are available at future Higgs factories.

      Figure 6.  (color online) 95% C.L. allowed region in the $m_{\rm SUSY}$ vs. $m_A$ plane with CEPC precisions for $X_t =0$ (zero-mixing, left panel) and $X_t = 2 m_{\rm SUSY}$ (max-mixing, right panel). For each panel, different colored curves correspond to different values of $\Delta m_h=1,\ 2, \ 3$ GeV, with the region to the right of the curve allowed.

    VII.   COMPARISON BETWEEN DIFFERENT HIGGS FACTORIES
    • To compare the reach for three different Higgs factories, we show the 95% C.L. allowed region in Fig. 7 in the parameter space of $ m_{\rm SUSY} $ vs. $ m_A $ (upper two panels) and $ \tan\beta $ vs. $ m_A $ (lower two panels) for the zero-mixing and max-mixing cases. While the CEPC and the FCC-ee have similar reach, the reach in $ m_A $ for the ILC is better because of better precisions on Higgs measurements of $ hbb $ and $ hWW $, given the increased rate of the WBF processes at higher center of mass energies. Limits of $ m_{\rm SUSY} $ and $ \tan\beta $ (for large $ m_{\rm SUSY} $) are nearly the same for all three Higgs factories because they are mainly controlled by the precision in Higgs mass, which comes from theoretical uncertainties. Note that for the $ m_A $ vs. $ \tan\beta $ plot, the CEPC results are more constraining than the FCC-ee results for $ m_{\rm SUSY} = 1 $ TeV, while they are nearly the same for $ m_{\rm SUSY} = 2 $ TeV. This is due to the slightly better precision in $ h\rightarrow gg, \gamma\gamma $ channel at the CEPC, which makes it more sensitive for smaller stop mass running in the loop.

      Figure 7.  (color online) 95% C.L. allowed region in the $m_{\rm SUSY}$ vs. $m_A$ plane (top row), and $\tan\beta$ vs. $m_A$ (bottom row) for $X_t =0$ (zero-mixing, left panels) and $X_t = 2 m_{\rm SUSY}$ (max-mixing, right panels), with the CEPC (solid), the FCC-ee (dotted) and the ILC (dashed) precision. The LHC Run-II direct search limits based on $A/H\to\tau\tau$ [48] are shown in the grey shaded region in the bottom panels.

    VIII.   CONCLUSION AND OUTLOOK
    • In this work, we have studied the constraints of Higgs precision measurements from future Higgs factories on the MSSM parameter space. We considered the dominant stop contributions to the Higgs mass, loop-induced $ h\gamma\gamma+hgg $ couplings, and the effective mixing angle $ \alpha_{\rm eff} $, which enters the Higgs couplings to a pair of fermions and gauge bosons, as well as additional loop contributions to the bottom Yukawa coupling $ \kappa_b $. The four relevant parameters under consideration are $ m_A $, $ \tan\beta $, $ m_{\rm SUSY} $ and $ X_t $.

      In the multi-variable $ \chi^2 $ fit, we included all the Higgs decay channels to SM fermions and gauge bosons at Higgs factories, as well as the Higgs mass. We found that $ \chi^2_{m_h} $ dominates for the small $ \tan\beta $ case, while $ \chi^2 $ contributions from the Higgs decays, in particular, $ h\rightarrow b\bar{b} $, dominates for the small to moderate $ m_A $ case. Generally we found that the CP-odd Higgs mass $ m_A $ is sensitive to the precisions of the Higgs decay channels, while $ m_{\rm SUSY} $, $ X_t $ and $ \tan\beta $ are sensitive to the precision of Higgs mass determination. For large $ \tan\beta $, $ m_{\rm SUSY} $ and $ X_t $ are also sensitive to the precisions of fermion and vector gauge boson couplings. For the max-mixing scenario, the loop-induced $ hgg $ and $ h\gamma\gamma $ couplings are the main restrictions on $ m_{\rm SUSY} $ when $ \tan\beta>7 $.

      We obtained the 95% C.L. allowed region given the Higgs factory precisions, and presented the result in the parameter space of $ M_{\rm SUSY} $ vs. $ X_t $, $ M_A $ vs. $ m_{\rm SUSY} $, $ m_A $ vs. $ \tan\beta $ and $ m_{\rm SUSY} $ vs. $ \tan\beta $. We found that small $ \tan\beta $ only survives in the max-mixing case with relatively large $ m_A $, while large regions of $ X_t $ vs. $ m_{\rm SUSY} $ are allowed for large $ \tan\beta $ and large $ m_A $. The lower limits on $ \tan\beta $ depends sensitively on the values of $ m_{\rm SUSY} $ and $ m_A $, in particular, for $ m_{\rm SUSY}<1.5 $ TeV and $ m_A<2 $ TeV. Limits on $ m_{\rm SUSY} $ also depend sensitively on $ \Delta m_h $, indicating the importance of a precise determination of the Higgs mass in the MSSM. For $ \tan\beta = 50 $ of the max-mixing scenario, $ m_{\rm SUSY} \in (0.8,1.2) \,{{\rm{GeV}}} $ when $ \Delta m_h = 1 $ GeV.

      We also compared the reach of the CEPC, the FCC-ee and the ILC. We found that the reach of the CEPC is similar to that of the FCC-ee, while the reach of the ILC is typically better, given the slightly better precision in the Higgs WBF measurements. With the high precision of the Higgs coupling measurements, and the potential improvement of theoretical calculation of $ m_h $ in the MSSM, studying the SM-like Higgs properties at future Higgs factories will offer great insight into the MSSM parameter space, which will be complementary to direct searches for SUSY particles at energy frontier machines.

    ACKNOWLEDGMENTS
    • We thank S. Heinemeyer for insightful discussions on the Higgs sector of the MSSM.

Reference (64)

目录

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return