-
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. -
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) -
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} $ . -
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]. -
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.
MSSM at future Higgs factories
- Received Date: 2020-11-22
- Available Online: 2021-04-15
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.