Mass-degenerate Higgs bosons near 125 GeV in the NMSSM under current experimental constraints

The observed Higgs signal at the Large Hadron Collider (LHC) may be derived from the two mass-degenerate resonances. We investigate this scenario in the Next-to-Minimal Supersymmetric Standard Model (NMSSM) in which both the two lightest CP-even Higgs bosons have masses around 125 GeV. We perform a comprehensive scan over the parameters in the NMSSM considering the current experimental constraints, especially the constraints from dark matter direct detection. The surviving samples are featured with relatively small $$\mu $$μ and large $$\tan \beta $$tanβ. The samples with large deviation of double ratios have large mixing between doublet field and singlet field.


Introduction
The discovery of Higgs boson with mass about 125 GeV by ATLAS and CMS collaborations [1,2] in 2012 stands for the beginning of a new era for high-energy physics. So far, the properties of the Higgs boson, such as the coupling with other Standard Model (SM) particles, have been measured and the measured properties are in accord with the prediction of SM. Significant improvements are expected from the upcoming runs of the Large Hadron Collider (LHC) or future e + e − colliders. However, the Higgs signal near 125 GeV observed at the LHC may also come from the contributions of two or more nearly mass-degenerate neutral Higgs bosons, and the Higgs resonances cannot be individually resolved by experiment. References [3][4][5] have developed different methods to detect whether two or more resonances exist. Some studies have also been done in the scenario with mass-degenerate Higgs states in the two-Higgsdoublet model (2HDM) [6,7] and Next-to-Minimal Supersymmetric Standard Model (NMSSM) [8][9][10][11][12][13].
Compared with the minimal supersymmetric standard model (MSSM), the NMSSM introduces an extra singlet a e-mail: zxheng@htu.edu.cn (corresponding author) fieldŜ. An effective μ-term is generated dynamically when S acquires the vacuum expectation value (vev), and μ problem in the MSSM can be solved elegantly [14,15]. Meanwhile, the Higgs sector of the NMSSM consists of three CPeven neutral scalars and two CP-odd neutral scalars, which is more complex in comparison to the MSSM. In addition, some experimental constraints on the NMSSM can be weaker due to the mixing between the singlet superfield and the doublet Higgs fields. As a result, some unique phenomenological properties of the Higgs boson may appear in the NMSSM [16,17].
Due to the interaction λŜĤ u ·Ĥ d among the Higgs fields, the squared mass of SM-like Higgs boson in the NMSSM receives a positive contribution at tree-level [18][19][20][21][22], and it can be further enhanced by the mixing between singlet field and doublet field when the next-to-lightest CP-even Higgs boson is SM-like [23][24][25]. For larger values of λ and smaller tan β, the positive contribution is substantial and large radiative correction to the Higgs boson mass is unnecessary and the little hierarchy problem can be avoidable. However, in some cases, for example smaller values of λ and larger tan β, the positive contribution at tree-level may be less than the loop-corrected contributions to achieve a mass around 125 GeV. About the loop corrections to the Higgs boson mass and Higgs boson decays in the NMSSM, a lot of works [26][27][28][29][30][31][32][33][34][35][36][37][38][39][40][41][42] have been done in the recent years. At present there are lots of public codes available for predictions in the NMSSM, such as NMSSMTools [43][44][45], NMSSMCalc [46], SARAH/SPheno [47][48][49][50][51][52][53][54], SoftSUSY [55] and FlexibleSUSY [56][57][58]. Considering the current experimental constraints, especially the constraints from the latest Dark Matter (DM) direct detection, we aim to study this scenario with two nearly mass-degenerate CP even Higgs bosons in the NMSSM, i.e. the observed Higgs signal at the LHC may be due to the superposition of this two resonances.
The paper is organized as follows. In Sect. 2 we will briefly review the basics of the scale invariant NMSSM, and perform a comprehensive scan over the parameter space considering the current experimental constraints. In Sect. 3, we will show the impact of constraints from DM direct detection on the parameter space of NMSSM. In Sect. 4 we show the numerical results of two mass-degenerate states in Higgs boson signals and discuss the main production channels and decay channels relevant for LHC data. Eventually, we give a summary in Sect. 5.

Basics of the NMSSM
Compared with the MSSM, the Higgs sector of the NMSSM is rather complicated by adding an extra singlet fieldŜ. The superpotential and the tree-level Higgs potential of the model are given by [14]: where W F represents the MSSM superpotential without the μ-term, the fieldsĤ u ,Ĥ d andŜ are Higgs superfields and their scalar components are H u , H d and S respectively, the dimensionless parameters λ and κ are the coupling strength in the Higgs sector, g 1 and g 2 are the U (1) Y and SU (2) L gauge couplings, and the remaining quantities m 2 H u ,H d ,S and A λ,κ are soft breaking parameters.
The tree-level Higgs potential of the NMSSM is made up of F-term and D-term of the superfields, and also the soft breaking terms. Using the minimization condition of the scalar potential, the soft breaking masses m 2 H u , m 2 H d and m 2 S can be replaced by v u , v d and v s , and v u , v d , v s are the vevs of the fields H u , H d and S, respectively. As a result, in the Higgs sector, six independent parameters in the following are needed, We define h 0 = cos β H u + ε sin β H * d and H 0 = sin β H u −ε cos β H * d with ε being two-dimensional antisymmetric tensor and ε 12 = −ε 21 = 1 and ε 11 = ε 22 = 0, then the Higgs fields in the CP-conserving NMSSM can be written as [59]: where G + and G 0 stand for Goldstone bosons and v 2 = v 2 u + v 2 d . The charged fields H ± are already physical mass eigenstates, and their masses at tree-level can be given by By diagonalizing the squared mass matrix M 2 with rotation matrix V , we can get the physical mass eigenstates . Similarly, one can also get the CP-odd mass eigenstates A 1 and A 2 . In general, we assume m H 1 < m H 2 < m H 3 and m A 1 < m A 2 , and call H i the SMlike Higgs boson if it is dominated by the S 2 field. Without the mixing among S i fields, the squared mass of SM-like Higgs boson receives an additional contribution λ 2 v 2 sin 2 2β compared with the case in MSSM. Furthermore, the mass of SM-like boson can also be enhanced by the (S 2 , Therefore, the mass of the SMlike Higgs boson can get its measured value without large radiative corrections [61][62][63].
For large tan β and large M A , the tree-level physical masses of the two lightest CP-even Higgs bosons may be approximately expressed as follows [59], Hence the mass difference between the tree-level masses should be as small as possible in order to obtain two massdegenerate Higgs bosons.

Scan strategies and constraints on the parameter space of NMSSM
From Eq. (6), we can see that the parameter A λ determines the size of M H ± . If the mixing among the doublet fields should be negligible, then M H ± should be large. So we fix A λ to be 2 TeV. 1 Considering the loop-corrected contributions to the Higgs boson mass, the two nearly mass-degenerate CP-even Higgs states strongly depend on the parameters λ, κ, tan β, μ, A κ and A t (A t is stop trilinear coupling), so we perform a random scan over these parameters in the following ranges: For other unimportant parameters including the soft breaking masses of gauginos, the soft breaking parameters in the slepton sector and squark sector except A t , we fix them to be 2 TeV.
In the calculation, we use the package NMSSMTools [43,44] to generate the particle mass spectrum, the relevant couplings and the decay branching ratios of Higgs bosons. We also consider the constraints on direct search for Higgs bosons from LEP, Tevatron and LHC with the package HiggsBounds [64,65], and perform fits for the 125 GeV Higgs data with the package HiggsSignals [66][67][68]. Moreover, we use the package micrOMEGAs [69,70] to compute the DM relic density and also the spin-independent (SI) and spin-dependent (SD) cross sections. We take the limits from LUX-2017 [71] for SD cross section and XENON1T-2018 [72] for SI cross section, respectively. We only require the DM relic density is less than the measured central value, i.e. h 2 < 0.1187. Since we assume DM in the NMSSM is only one of DM candidates [73], the SI and SD cross sections for DM-nucleon scattering should also be scaled by a factor ( h 2 )/0.1187.

Impact of constraints from DM direct detection on the parameter space of NMSSM
In contrast to limits from direct searches for Higgs bosons and DM relic density, DM direct detection experiments has a strong constraint on the parameter space of NMSSM when we require there exist two nearly mass-degenerate Higgs bosons. At present, the strong limits on the SI and SD cross sections of DM-nucleon scattering come from the XENON-1T experiment in 2018 and the LUX measurement in 2017, respectively. Considering the constraints from the DM relic density, Higgs data from LEP, Tevatron and LHC experiments, we calculate the SI and SD cross sections separately and project the samples on the σ − mχ0 1 plane in Fig. 1. In the left plane, samples with green color above (below) the red dashed line have been excluded by the XENON-1T experiment (LUX experiment). In the right plane, samples with cyan color above (below) the magenta dash-dotted line have been excluded by the LUX experiment (XENON-1T experiment). The samples with red color in the left plane and samples with magenta color in the right plane are still experimentally allowed. We can see that the two experiments are complementary to limit the parameter space of NMSSM. In the following, we will discuss the features of the parameter space in the NMSSM, which is limited tightly by the DM direct detection experiments.
In Fig. 2, we project the samples on λ − κ plane, λ − A t plane and μ − tan β plane. The samples marked by green color have been excluded by DM direct detection experiments and the red ones are still experimentally allowed. We see that the samples with large λ and small κ shown in Fig. 2a and also samples with large μ shown in Fig. 2c are excluded by the DM direct searches. This is because that DM corresponding to these samples are composed mainly of higgsinos and the couplings of Higgs-DM-DM are relatively large. The figure also shows that the parameter space of NMSSM shrinks somewhat after considering the constraints of DM direct detection. The surviving parameter space of NMSSM is featured by 100 GeV μ 250 GeV and 6 tan β 25. For most of the surviving samples, λ is not too large, tan β is relatively large and the value of parameter A t is relatively large. Hence, the tree-level contributions for the Higgs boson mass from λ 2 v 2 sin 2 2β is not large and the contributions from the stop-loop radiative corrections are important to obtain two mass-degenerate Higgs boson near 125 GeV.

Searching for two mass-degenerate Higgs bosons near 125 GeV
At the LHC, the main production channels for Higgs boson are gluon-gluon fusion and vector boson fusion (VBF). We focus on these two main production modes (denoted by Y = gg or Y = VBF) and consider the Higgs decay channels (denoted by X ) to γ γ, Z Z * , bb and τ + τ − . We define the following ratio where To test the existence of two mass-degenerated Higgs states, Ref. [3] defined the following double ratios, Each of these double ratios should be unity if only a single Higgs boson contributes to the 125 GeV SM-like Higgs signal. However, these ratios may be deviated from 1 if two Higgs bosons contribute to the observed Higgs signal.
After implementing all the constraints mentioned above, we perform a global χ 2 fit to a total of 102 most up-todate Higgs boson observables at the LHC by the package HiggsSignals-2.2.3beta. Then in Fig. 3 we display the double ratios D 1 , D 2 , D 3 as functions of the mass difference m H 2 − m H 1 with the color-bar representing the total χ 2 . From the figure, we can see that the three ratios are very close to unity, and the maximum deviation can reach to 10%. Both D 1 and D 3 can be larger or smaller than 1 because they depend on the same decay channel H i → bb, while D 2 is not bigger than 1 because it relies on only the bosonic signal strengths. Meanwhile, the small values of χ 2 correspond to samples with small deviation of each ratios and the samples with large deviation correspond to relatively larger χ 2 , such as for samples with deviation reach to 10% (4%), the value of χ 2 is about 125 (100) in Fig. 3a, c, and samples with deviation of 3% (1%), the value of χ 2 is about 120 (100) in Fig. 3b. Figure 4 shows the field components of H 1 and H 2 , respectively. For the left (right) upper corner samples, H 1 (H 2 ) is dominated by singlet field S 3 and H 2 (H 1 ) consists mainly of SM-like Higgs field S 2 . Their corresponding Higgs signal characters D i (i = 1, 2, 3) are around one, which stands for the singlet dominated Higgs make little contribution to the Higgs signal. But for the middle samples, the mixing between S 2 and S 3 fields is large, both H 1 and H 2 have relatively large components of SM-like Higgs filed S 2 and their double ratios D i (i = 1, 2, 3) slightly deviate from one. That is to say, the observed Higgs signal may be due to the contributions from the nearly mass-degenerate Higgs bosons near 125 GeV.

Summary
Considering the current experimental limits from LHC Higgs data and DM relic density, especially limits from DM direct detection experiments, there still exist samples predicting two nearly mass-degenerate neutral CP-even Higgs bosons lie in the mass range between 122 and 128 GeV. The surviving samples are featured with relatively small μ and large tan β. We can use the double ratios defined in Eq. (10) to distinguish the existence of two mass-degenerate Higgs bosons. The maximum deviation of the double ratio from 1 can reach to about 10%. The samples with large deviation of each ratios correspond to large mixing between the S 2 and S 3 fields.

Note added
We note that varying the bino parameter M 1 would affect the properties of DM somewhat and we will talk this scenario in detail in our next work.

Data Availability Statement
This manuscript has no associated data or the data will not be deposited. [Authors' comment: This is a theoretical study and there is no experimental data associated with it.] Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permit-ted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecomm ons.org/licenses/by/4.0/. Funded by SCOAP 3 .