Light Higgs bosons in the general NMSSM

Physics beyond the Standard Model (SM) may manifest itself as small deviations from the SM predictions for Higgs signal strengths at 125 GeV. Then, a plausible and interesting possibility is that the Higgs sector is extended and at the weak scale there appears an additional Higgs boson weakly coupled to the SM sector. Combined with the LEP excess in $e^+e^-\to Z(h\to b\bar b)$, the diphoton excess around 96 GeV recently reported by CMS may suggest such a possibility. We examine if those LEP and CMS excesses can be explained simultaneously by a singlet-like Higgs boson in the general next-to-minimal supersymmetric Standard Model (NMSSM). Higgs mixing in the NMSSM relies on the singlet coupling to the MSSM Higgs doublets and the higgsino mass parameter, and thus is subject to the constraints on these supersymmetric parameters. We find that the NMSSM can account for both the LEP and CMS excesses at 96 GeV while accommodating the observed 125 GeV SM-like Higgs boson. Interestingly, the required mixing angles constrain the heavy doublet Higgs boson to be heavier than about 500 GeV. We also show that the viable region of mixing parameter space is considerably modified if the higgsino mass parameter is around the weak scale, mainly because of the Higgs coupling to photons induced by the charged higgsinos.


Introduction
Although the Standard Model (SM) successfully describes the observed particle physics up to energy scales around TeV, it is clear that a more fundamental theory is needed to provide a complete description of nature. So far the LHC has seen no clear signal for physics beyond the SM, and the discovered 125 GeV Higgs boson has properties compatible with the SM predictions [1,2]. Yet an interesting possibility is that the Higgs sector is extended to include an additional light Higgs boson which is accessible to collider experiments and result in some deviations of the 125 GeV Higgs boson from the SM predictions. The CMS has recently announced that Higgs searches in the diphoton final state show a local excess of about 3σ at 96 GeV [3]. The results from the ATLAS do not show a relevant excess, but are well compatible with the CMS limit [4]. Combined with the 2.3σ local excess observed in the LEP searches for e + e − → Z(h → bb) [5,6], the CMS results provide a motivation to consider the possibility that the Higgs sector involves an additional scalar boson at 96 GeV, which has been studied recently in refs. [7][8][9][10][11][12][13][14][15][16][17].
In this paper we explore if the next-to-minimal supersymmetric SM (NMSSM) can explain the LEP and CMS excesses around 96 GeV while accommodating the observed 125 GeV Higgs boson. The NMSSM extends the Higgs sector to include a gauge singlet scalar which generates the higgsino mass parameter µ via its coupling λ to the MSSM Higgs doublets [18,19]. As noticed in refs. [20,21], there are intriguing relations between Higgs mixings and the model parameters λ and µ that hold for the general NMSSM. Those relations are quite useful when examining how much Higgs mixings, which determine how the neutral Higgs bosons couple to SM particles [22][23][24][25][26][27][28][29], are constrained by the requirements on the model such as radiative corrections to the Higgs masses, the perturbativity bound on λ, and the chargino mass limit on µ. The viable region of mixing parameter space would be further constrained if one specifies singlet self-interactions. For instance, there are no tadpole and quadratic terms for the singlet in the Z 3 -symmetric NMSSM, 1 for which the mixing between the neutral singlet and doublet Higgs bosons has a dependence on the mass of the CP-odd singlet scalar.
Our analysis is based on the relations between Higgs mixings and the model parameters, and is performed for the general NMSSM without specifying singlet self-interactions. We first examine if a singlet-like Higgs boson at 96 GeV can be responsible for the LEP and CMS excesses within the range of mixing angles allowed by the current LHC data on the 125 GeV Higgs boson, under the assumption that the gauginos, squarks and sleptons are heavy enough, above TeV as indicated by the LHC searches for supersymmetry (SUSY), while the higgsinos can be significantly lighter. We then impose the constraints on λ and µ to find the viable mixing angles. It turns out that the general NMSSM can accommodate the SM-like 125 GeV Higgs boson compatible with the current LHC data, and also a singlet-like 96 GeV Higgs boson explaining both the LEP and CMS excesses. The allowed range of mixing angles is considerably modified if µ is around the weak scale because the charged higgsinos enhance the Higgs coupling to photons. Interestingly, if the excesses around 96 GeV are due to the singlet-like Higgs boson, the heavy doublet Higgs boson should be heavier than about 500 GeV. This paper is organized as follows. In section 2, we briefly discuss the effects of the neutral Higgs boson mixings on Higgs phenomenology and examine the relations between the mixing angles and the NMSSM parameters. The region of the mixing parameter space compatible with the current LHC data on the SM-like Higgs boson is presented in section 3. Section 4 is devoted to our main results, which show the mixing angles required to explain the LEP and CMS excesses, while satisfying the various constraints on the NMSSM parameters. It is also shown that the allowed mixing angles constrain the heavy Higgs boson to have a mass in a certain range. The final section is for the summary and comments.

Higgs bosons in the general NMSSM
In this section we describe how the neutral Higgs boson mixings depend on the NMSSM parameters, in particular, on the singlet coupling λ to the MSSM Higgs doublets and the higgsino mass parameter µ. Such relations should be taken into account when examining the Higgs mixings consistent with the experimental constraints. We also discuss the properties of Higgs bosons within the low energy effective theory constructed by integrating out heavy superparticles under the assumption that the higgsinos can be light. Note that our approach is applicable to a general form of NMSSM.

Dependence of Higgs mixing on NMSSM parameters
Taking an appropriate redefinition of superfields, one can always write the superpotential of the general NMSSM as with a canonical Kähler potential. Here H u and H d are the Higgs doublet superfields, and S is the gauge singlet superfield. There are various types of NMSSM, depending on the form of the singlet superpotential f (S). Our subsequent discussion applies for a general form of f (S), but for simplicity we will assume no CP violation in the Higgs sector. After the electroweak symmetry breaking, the CP-even neutral Higgs bosons mix with each other due to the potential terms where H 0 u = v sin β and H 0 d = v cos β with v = 174 GeV, and A λ is the SUSY breaking trilinear coupling. With the above Higgs potential terms, the effective Higgs µ and Bµ parameters are given by The supersymmetric parameters λ and µ, on which the Higgs mixing depends, are subject to the perturbativity bound and the chargino mass bound, respectively. Imposing the conditions for the electroweak symmetry breaking, the mass squared matrix for (ĥ,Ĥ,ŝ) reads 2M HereM 2 11 ,M 2 12 , andM 2 22 include radiative corrections, which can be sizable as arising from top and stop loops [30]: The value ofM 2 33 is determined by the singlet superpotential f (S) and the associated soft SUSY breaking terms F (S) according tô where all the terms are evaluated at the vacuum.
is the geometric mean of the two stop mass-eigenvalues, and A t is the SUSY breaking trilinear coupling associated with the top quark Yukawa y t = m t /v. Note that m 0 corresponds to the SM-like Higgs boson mass at large tan β in the decoupling limit of the MSSM. The LHC results constrain the stops to be heavier than TeV, and thus m 0 cannot be lower than about 115 GeV as long as stop mixing has X 2 t 10M 2 S as is the case in the conventional mediation schemes of SUSY breaking. The stop loop corrections maximizes m 0 at X t = ± √ 6M S , i.e. for maximal stop mixing. On the other hand, the charged Higgs boson has a mass around the square-root ofM 2 22 , and it should be heavier than about 350 GeV to avoid the experimental constraint associated with b → sγ [31].
To find the mass eigenstates, one needs to diagonalize the mass squared matrix as where the orthogonal mixing matrix U can be parametrized as with s θ ≡ sin θ and c θ ≡ cos θ, for which the angles θ 1 , θ 2 and θ 3 representĥ-Ĥ,ŝ-ĥ and s-Ĥ mixing, respectively. Obviously each matrix element ofM 2 can be expressed in terms of the mass eigenvalues {m h , m H , m s } and the mixing angles {θ 1 , θ 2 , θ 3 }. Among such relations, the following ones are particularly relevant for our subsequent discussion: because the constraints on the parameters λ, µ and m 0 can be converted into those on the Higgs boson masses and mixing angles, and also vice versa. Here one should note that ∆m 2 12 in eq. (2.11) is written with given by The above shows that vanishes at X t = 0 and X t = ± √ 6M S , and therefore ∆m 2 12 is tightly correlated with m 0 near the regions of minimal and maximal stop mixing. The correlation of the stop corrections in large regions of parameter space has been noted in ref. [30]. It is also easy to see that | | is smaller than about 0.1 for X 2 t 10M 2 S . For stop mixing with X t far from 0 or ± √ 6M S , the -contribution to ∆m 2 12 can be sizable only when |µ| is large, close to M S , and tan β is low. Keeping this feature in mind, we shall neglect the -contribution in our analysis unless stated otherwise.

Effective Higgs couplings to the SM sector
At energy scales around the electroweak scale, the properties of the Higgs bosons can be examined within an effective theory constructed by integrating out heavy sparticles. The LHC results on the Higgs sector and the searches for new physics indicate that SUSY, if exists, would be broken at a scale above TeV. Taking this into account, we assume that all gauginos, squarks and sleptons have masses above TeV, while the higgsinos and additional Higgs bosons can be lighter than TeV. Then the effective lagrangian describing how the neutral Higgs bosons interact with the SM fermions and gauge bosons is written as [32] L At tree level, the Higgs couplings to massive SM particles are given by The Higgs couplings to gluons and photons are radiatively generated, which results in where δC i g and δC i γ are additional contributions from sparticle loops, and the loop functions are given by For the case when the superpartners of SM particles are heavier than TeV, sparticles give negligible contributions to ∆C i g . However, if µ is small, the charged higgsinos are light and can give a sizable contribution to ∆C i γ through the following Higgs-higgsino couplings Under the assumption that the gauginos are significantly heavier than the higgsinos, the Higgs coupling to photons induced by charginos can be approximated to be for a Higgs boson with m 2

Mixing consistent with the 125 GeV Higgs boson
For small scalar mixing, h has properties close to those of the SM Higgs boson. In this paper, we identify h with the SM-like Higgs boson observed at the LHC and examine how the scalar mixing is constrained by the measured signal strengths. The SM-like Higgs boson has m h = 125 GeV, and its couplings to the massive SM particles are given by while the couplings to gluons and photons read ∆C h g 0.97c θ 1 c θ 2 + (1.03 cot β + 0.06 tan β)s θ 1 , Here r is defined by and measures the chargino contribution, which has been approximated by using the fact that it is non-negligible only when |µ| is not far above the electroweak scale for λ below the perturbative bound, and the chargino search at LEP requires |µ| > 104 GeV. One should note that the couplings of the SM-like Higgs boson are determined by four parameters: θ 1 , θ 2 , r and tan β. The partial decay rates of the SM-like Higgs boson h can be easily estimated by using the well-known decay properties of the hypothetical Higgs boson φ 125 of the minimal SM with mass 125 GeV: (3.4) Assuming that h does not decay to non-SM particles, one also finds its total decay rate to be with Γ tot (φ 125 ) being the total decay width of φ 125 . Here we have used the branching ratios of the SM Higgs boson listed in ref. [33]. The production of the SM-like Higgs boson is dominated by the gluon fusion process, and the signal strength normalized by the SM value is given by (3.7) It is obvious that one should have µ ii h = 1 if θ 1 = θ 2 = 0 and r = 0. The ATLAS collaboration has recently updated the measurements on the Higgs signal strengths using the 13 TeV data [ Note that the NMSSM leads to µ W W h = µ ZZ h , which is within the 1σ range. Let us now examine how severely the Higgs mixing is constrained by the LHC experimental results on the Higgs boson at 125 GeV. The signal rates of h are determined by two for which the branching ratio for h → V V is suppressed compared to the case of the SM Higgs boson, but such effect is compensated by the enhancement of production rate via the gluon fusion process. Note that r, which is relevant for the diphoton signal strength, is below about 1.2 because λ should be smaller than about 0.7 in order for the NMSSM to remain perturbative up to the GUT scale, and |µ| should be larger than 104 GeV to satisfy the LEP bound on the chargino mass. Figure 1 shows the 2σ range of (θ 1 , θ 2 ) allowed by the current LHC data on the Higgs boson. The left panel is for r = 0.1, for which the Higgs coupling to photons is rarely affected by the charged higgsinos, and the right panel is for r = 1. Here we have taken 1.5 ≤ tan β ≤ 15 to see how the allowed region changes with tan β. From the figure, one can see that the allowed region of θ 1 gets smaller if tan β increases, but a broad range of θ 2 is allowed insensitively to tan β. This is because the Higgs coupling C h b gets sensitive to θ 1 at large tan β while the couplings to the top quark and gauge bosons do not. If r is around unity, the charged higgsinos can significantly enhance the diphoton signal rate, 4 Although we will not pursue in this paper, a large value of θ1 satisfying can also lead to µ V V h = 1. In this case, C h b is negative and thus leads to wrong sign Yukawa couplings [20], and one needs a large λ beyond the perturbativity bound [34]. excluding θ 2 in the range between about 0.4 and 0.8. We can understand this feature from the fact that the charged higgsinos induce a Higgs coupling to photons, δC h γ ≈ −0.17rθ 2 for small mixing angles, whereas the Higgs couplings to other SM particles only quadratically depend on θ 2 .

LEP and CMS excesses around 96 GeV
The CMS collaboration has recently reported a local excess of 2.8σ in the diphoton channel around m γγ = 96 GeV [3]. The signal strength amounts to where φ 96 denotes the hypothetical SM Higgs boson with mass 96 GeV, and Br(ϕ → γγ) denotes the branching ratio for the diphoton channel [13]. Intriguingly, there is another 2.3σ local excess at the similar mass region from the Higgs searches in the Z-boson associated Higgs production (e + e − → Zϕ) at LEP [5]. The signal strength is [35] µ LEP = σ(e + e − → Zϕ) σ(e + e − → Zφ 96 ) It has long been known that the LEP excess can be explained by a light singlet-like Higgs boson in the NMSSM. At this stage, a naturally occurring question is whether this singletlike Higgs boson can explain the CMS diphoton excess as well.
If both excesses were arisen due to the light singlet-like Higgs boson s, the signal strengths can be expressed by the effective couplings in eq. (2.14) as follows, assuming that the CP-odd singlet scalar and the singlino are heavy enough so that s decays only into the SM particles. Here we have used HDECAY [36,37] to calculate the decay widths of the 96 GeV Higgs boson with the SM couplings. The effective couplings of s are written in terms of the mixing angles as and those to gluons and photons read including the contribution from the loops of charged higgsinos.
Before performing a numerical analysis, we present approximate analytic relations between mixing angles holding if s is responsible for the LEP and CMS excesses. The LEP signal rate given in eq. (4.4) is approximated by for small θ 2 , and thus the LEP excess is explained if s 2 θ 2 ∼ 0.1. On the other hand, the ratio between the CMS and LEP signal rates is given by where k 32 ≡ s θ 3 /s θ 2 , and it should be around 6 to account for the LEP and CMS excesses.
Here the first factor of the right-hand side represents the Higgs production dominated by gluon fusion, while the second one concerns the branching ratio for the diphoton mode. Note that both effects are enhanced if k 32 is negative, and the charged higgsinos can further increase the branching ratio into photons for s θ 2 < 0. These features help to understand the numerical analysis given below.
Let us explain in detail how to search the viable region of mixing angles in the parameter scan. The signal rates of h and s are functions of the mixing angles, a combination of λ and µ, and tan β: where we have taken = 0 in eq. (2.12) since it is negligibly small in most of the parameter space of our interest. It is obvious that m H is determined by θ 1,2,3 and tan β once we fix r. The above relations allow us to analyze the viable mixing angles as follows. We first examine the (θ 1 , θ 2 ) space to see in which region µ ii h are consistent with the current LHC data, and then continue to check if it is further possible to explain both µ LEP and µ CMS .
In figure 2, the Higgs signal rates µ ii h are consistent with the measurements in the shaded region, and the excesses at 96 GeV are explained in the red shaded region. Figure 3 shows the region of θ 3 for the excesses at 96 GeV. The next thing that one has to examine is if the viable region above is consistent also with the various constraints on λ, µ, and m 0 . Here we have imposed the conditions λ < 0.7, |µ| > 104 GeV, m 0 ≥ 115 GeV, (4.11)  as required by the perturbativity up to the GUT scale, the LEP limit on the chargino mass, and the radiative contribution to the Higgs mass from stops above TeV, respectively, provided that stop mixing is not too large as would be the case in the conventional medi- ation models of SUSY breaking. 5 Then it follows r ≤ 1.1. Note that r parameterizes the radiative effect of the charged higgsinos on Higgs decays. As benchmark points, we have taken r = 0.1, 1 for 1.5 ≤ tan β ≤ 15. 6 Each color in figure 3 represents how much the above constraints reduce the viable region. We note that the bound on m 0 gets important when tan β is small and r is around 1 or above.
In the parameter region for the LEP and CMS excesses, the main effect of stop loop corrections ∆m 2 12 is to increase (decrease) the coupling λ if ∆m 2 12 is negative (positive), as can be deduced from the last relation in eq. (2.11). This implies that the parameter space compatible with both excesses shrinks for larger negative ∆m 2 12 due to the perturbativity bound on λ. In the parameter region with m 0 ≥ 115 GeV, however, the dependence on ∆m 2 12 becomes quite weak because λ should be small in order to get m h = 125 GeV. On the other hand, the LEP limit on the chargino mass given in eq. (4.11) implies λ > 0.6r, following from r ≡ λv/|µ|. Taking this together with the perturbativity bound on λ, one can find that λ would be more severely constrained at larger r when the stop correction ∆m 2 12 is negative. We have checked these features by taking analysis for nonzero values of between −0.05 and 0.05.
We close this section by pointing out that the LEP and CMS excesses can constrain the masses of the heavy Higgs boson and higgsinos, if they are due to the singlet-like Higgs boson. Eqs. (2.11) enable us to extract the information on the region of µ and m H 5 The lower bound on m0 is set by considering the numerical result that m0 is around 113 GeV for stops at 1 TeV in the limit of vanishing stop mixing [38][39][40]. 6 We have cross checked the results of our analysis for some parameter points by using NMSSMTools [38][39][40].
compatible with the Higgs signal strengths, µ ii h , µ LEP and µ CMS . Figure 4 shows the allowed region of (µ, m H ), where we have taken tan β = 2 (left) and 5 (right) with 0 < r < 1.1. As discussed already, the m 0 cut is relevant for small tan β. It is important to note that the CMS and LEP excesses put a lower and upper bound on m H . The lower bound turns out to be m H 500 GeV, nearly irrespectively of the values of µ and tan β, while the upper bound depends on those parameters and is found to increase with tan β.

Summary
Extended with an additional gauge singlet scalar, the Higgs sector of the NMSSM offers a rich phenomenology to be explored at collider experiments. In particular, as experimentally allowed to be light, a singlet-like Higgs boson could be observable in the searches for e + e − → Z(h → bb) and pp → h → γγ if it couples to the SM sector via the Higgs mixing. It is thus interesting to examine if the excesses reported by LEP and CMS in those channels can be interpreted as signals of a singlet-like Higgs boson with mass around 96 GeV within the NMSSM.
For the case that the gauginos, squarks and sleptons have masses above TeV, while the Higgsinos can be significantly lighter, which is perfectly consistent with the null results for SUSY searches at LHC so far, we have found that the general NMSSM can successfully accommodate such a light singlet-like Higgs boson explaining the LEP and CMS excesses simultaneously, as well as the 125 GeV Higgs boson compatible with the current LHC data. The range of mixing angles required to explain the 96 GeV excesses can be considerably modified if the higgsinos are around the weak scale, because the singlet-like Higgs coupling to photons is enhanced.
To examine a viable region of mixing parameter space, it should be taken into account that Higgs mixing is subject to various constraints on the NMSSM parameters. We have shown that, if a singlet-like Higgs boson is responsible for the LEP and CMS excesses, Higgs mixing is strongly constrained by the LEP bound on the charged higgsino mass and the perturbativity bound on the singlet coupling to the Higgs doublets. Interestingly, in the viable mixing space, the heavy doublet Higgs boson is found to be heavier than about 500 GeV.
The physics underlying the electroweak symmetry breaking may manifest itself as slight deviations from the SM predictions for the Higgs signal strengths at 125 GeV. It is then a plausible possibility that there exist additional light Higgs bosons weakly coupled to the SM sector, which would provide crucial information on how the Higgs sector is extended. The excesses reported by LEP and CMS, both of which are interestingly around 96 GeV, would thus deserve more attention.