Light Higgs bosons in the general NMSSM

TL;DR

This paper investigates whether the general NMSSM can simultaneously explain the LEP and CMS excesses around 96 GeV while maintaining the observed 125 GeV Higgs, highlighting the role of Higgs mixing and supersymmetric parameters.

Contribution

It demonstrates that the NMSSM can account for both excesses with specific Higgs mixing, constraining heavy Higgs masses and emphasizing the impact of higgsino mass parameters.

Findings

NMSSM can explain LEP and CMS excesses at 96 GeV.

Heavy doublet Higgs mass constrained to > 500 GeV.

Higgs coupling to photons affected by higgsino mass parameter.

Abstract

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.

Equations

$$W=\lambda S{H}_u{H}_d+f(S)+(\text{MSSM Yukawa terms}),$$

$$\left(\begin{array}{c}\hat{h}\\ \hat{H}\\ \hat{s}\end{array}\right)=\sqrt{2}\left(\begin{array}{ccc}[r]\sin \beta & \cos \beta & 0\\ -\cos \beta & \sin \beta & 0\\ 0& 0& 1\end{array}\right)\left(\begin{array}{c}\mathrm{Re}{H}_u^0-v\sin \beta \\ \mathrm{Re}{H}_d^0-v\cos \beta \\ \mathrm{Re}S-⟨S⟩\end{array}\right)$$

$$V_{\rm mix}=\lambda^{2}|S|^{2}(|H_{u}|^{2}+|H_{d}|^{2})+\Big{(}A_{\lambda}\lambda SH_{u}H_{d}+(\partial_{S}f)^{*}\lambda H_{u}H_{d}+{\rm h.c.}\Big{)},$$

$$\mu$$

$$B\mu$$

$${\hat{M}}_{33}^2=({\partial }_S^{2}f{)}^2+\left({\partial }_{S}f-\frac{1}{2}\lambda v^2\sin 2\beta \right)\left({\partial }_S^{3}f-\frac{{\partial }_S^{2}f}{S}\right)+\frac{1}{2}\lambda v^2\frac{A_{\lambda }}{S}\sin 2\beta +\left({\partial }_S^{2}F-\frac{{\partial }_{S}F}{S}\right),$$

$${\hat{M}}_{11}^2$$

$${\hat{M}}_{12}^2$$

$${\hat{M}}_{13}^2$$

$${\hat{M}}_{22}^2$$

$${\hat{M}}_{23}^2$$

$$\Lambda \equiv A_{\lambda }+⟨{\partial }_S^{2}f⟩.$$

$$m_0^2$$

$$\Delta m_{12}^2$$

$$\Delta m_{22}^2$$

$$U{\hat{M}}^2{U}^{\mathrm{T}}=\mathrm{diag}(m_h^2,m_H^2,m_s^2),$$

$$U=\left(\begin{array}{ccc}{c}_{{\theta }_1}{c}_{{\theta }_2}& -s_{{\theta }_1}& -c_{{\theta }_1}{s}_{{\theta }_2}\\ s_{{\theta }_1}{c}_{{\theta }_2}{c}_{{\theta }_3}-s_{{\theta }_2}{s}_{{\theta }_3}& c_{{\theta }_1}{c}_{{\theta }_3}& -c_{{\theta }_2}{s}_{{\theta }_3}-s_{{\theta }_1}{s}_{{\theta }_2}{c}_{{\theta }_3}\\ s_{{\theta }_1}{c}_{{\theta }_2}{s}_{{\theta }_3}+s_{{\theta }_2}{c}_{{\theta }_3}& c_{{\theta }_1}{s}_{{\theta }_3}& c_{{\theta }_2}{c}_{{\theta }_3}-s_{{\theta }_1}{s}_{{\theta }_2}{s}_{{\theta }_3}\end{array}\right)$$

$$m_0^2$$

$$\lambda v\mu$$

$${\lambda }^2{v}^2$$

$$\Delta m_{12}^2=-\frac{m_0^2-m_Z^2}{\tan \beta }+ϵv^2,$$

$$ϵ=-\frac{3y_t^4}{8{\pi }^2}\frac{\mu }{M_S}\left(1+\frac{1}{{\tan }^2\beta }\right)\left(\frac{X_t}{M_S}-\frac{X_t^3}{6M_S^3}\right).$$

$${\mathcal{L}}_{\mathrm{eff}}$$

$$C_V^i=U_{i1},C_t^i=U_{i1}-U_{i2}\cot \beta ,C_b^i=C_{\tau }^i=U_{1i}+U_{i2}\tan \beta .$$

$$\Delta C_g^i$$

$$\Delta C_{\gamma }^i$$

$$A_{1/2}(\tau )$$

$$A_1(\tau )$$

$$f(\tau )=\{\begin{array}{ll}{\arcsin }^2\sqrt{\tau }& \mathrm{for}\tau \le 1\\ -\frac{1}{4}{\left(\ln \frac{1+\sqrt{1-{\tau }^{-1}}}{1-\sqrt{1-{\tau }^{-1}}}-i\pi \right)}^2& \mathrm{for}\tau >1\end{array}.$$

$$\frac{\lambda }{\sqrt{2}}\sum _i{U}_{i3}{\varphi }_i{\tilde{H}}_u^{+}{\tilde{H}}_d^{-}.$$

$$\delta C_{\gamma }^i\simeq \frac{\lambda v}{6|\mu |}\left(1+\frac{7}{30}\frac{m_{{\varphi }_i}^2}{4|\mu {|}^2}\right)U_{i3}$$

$$C_V^h=c_{{\theta }_1}{c}_{{\theta }_2},C_t^h=c_{{\theta }_1}{c}_{{\theta }_2}+s_{{\theta }_1}\cot \beta ,C_b^h=C_{\tau }^h=c_{{\theta }_1}{c}_{{\theta }_2}-s_{{\theta }_1}\tan \beta ,$$

$$\Delta C_g^h$$

$$\Delta C_{\gamma }^h$$

$$r\equiv \frac{\lambda v}{|\mu |},$$

$$\frac{\Gamma (h\to bb)}{\Gamma ({\varphi }_{125}\to bb)}$$

Full text

1 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 ATLAS-CONF-2019-005 ; Sirunyan:2018koj . 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\sigma$ at $96$ GeV Sirunyan:2018aui . The results from the ATLAS do not show a relevant excess, but are well compatible with the CMS limit ATLAS-CONF-2018-025 . Combined with the $2.3\sigma$ local excess observed in the LEP searches for $e^{+}e^{-}\to Z(h\to b\bar{b})$ Barate:2003sz ; Schael:2006cr , 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. Fox:2017uwr ; Richard:2017kot ; Haisch:2017gql ; Biekotter:2017xmf ; Tao:2018zkx ; Liu:2018xsw ; Domingo:2018uim ; Hollik:2018yek ; Wang:2018vxp ; Biekotter:2019kde ; Cline:2019okt .

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 $\mu$ via its coupling $\lambda$ to the MSSM Higgs doublets Maniatis:2009re ; Ellwanger:2009dp . As noticed in refs. Choi:2012he ; Choi:2013lda , there are intriguing relations between Higgs mixings and the model parameters $\lambda$ and $\mu$ 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 Ellwanger:2010nf ; Jeong:2012ma ; Badziak:2013bda ; Jeong:2014xaa ; Ellwanger:2015uaz ; Beskidt:2017dil ; Jeong:2017jqp ; LiuLiJia:2019kye , are constrained by the requirements on the model such as radiative corrections to the Higgs masses, the perturbativity bound on $\lambda$, and the chargino mass limit on $\mu$. 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,111The possibility of accommodating both LEP and CMS excesses in the $Z_{3}$-symmetric NMSSM was firstly explored in ref. Domingo:2018uim . Our study considers the general NMSSM without any additional symmetry or matter. 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 $\lambda$ and $\mu$ 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 $\mu$ 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.

2 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 $\lambda$ to the MSSM Higgs doublets and the higgsino mass parameter $\mu$. 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.

2.1 Dependence of Higgs mixing on NMSSM parameters

Taking an appropriate redefinition of superfields, one can always write the superpotential of the general NMSSM as

[TABLE]

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

[TABLE]

mix with each other due to the potential terms

[TABLE]

where $\langle H^{0}_{u}\rangle=v\sin\beta$ and $\langle H^{0}_{d}\rangle=v\cos\beta$ with $v=174$ GeV, and $A_{\lambda}$ is the SUSY breaking trilinear coupling. With the above Higgs potential terms, the effective Higgs $\mu$ and $B\mu$ parameters are given by

[TABLE]

The supersymmetric parameters $\lambda$ and $\mu$, 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 $(\hat{h},\hat{H},\hat{s})$ reads222 The value of $\hat{M}^{2}_{33}$ is determined by the singlet superpotential $f(S)$ and the associated soft SUSY breaking terms $F(S)$ according to

$\hat{M}^{2}_{33}=(\partial_{S}^{2}f)^{2}+\left(\partial_{S}f-\frac{1}{2}\lambda v^{2}\sin 2\beta\right)\left(\partial_{S}^{3}f-\frac{\partial_{S}^{2}f}{S}\right)+\frac{1}{2}\lambda v^{2}\frac{A_{\lambda}}{S}\sin 2\beta+\left(\partial_{S}^{2}F-\frac{\partial_{S}F}{S}\right),$

(5)

where all the terms are evaluated at the vacuum.

[TABLE]

with $\Lambda$ defined by

[TABLE]

Here $\hat{M}^{2}_{11}$, $\hat{M}^{2}_{12}$, and $\hat{M}^{2}_{22}$ include radiative corrections, which can be sizable as arising from top and stop loops Carena:2015moc :

[TABLE]

with $X_{t}=A_{t}-\mu\cot\beta$ and $Y_{t}=A_{t}+\mu\tan\beta$, where $M_{S}=\sqrt{m_{\tilde{t}_{1}}m_{\tilde{t}_{2}}}$ 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\beta$ 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}\lesssim 10M_{S}^{2}$ as is the case in the conventional mediation schemes of SUSY breaking. The stop loop corrections maximizes $m_{0}$ at $X_{t}=\pm\sqrt{6}M_{S}$, i.e. for maximal stop mixing. On the other hand, the charged Higgs boson has a mass around the square-root of $\hat{M}^{2}_{22}$, and it should be heavier than about $350$ GeV to avoid the experimental constraint associated with $b\to s\gamma$ Gambino:2001ew .

To find the mass eigenstates, one needs to diagonalize the mass squared matrix as

[TABLE]

where the orthogonal mixing matrix $U$ can be parametrized as

[TABLE]

with $s_{\theta}\equiv\sin\theta$ and $c_{\theta}\equiv\cos\theta$, for which the angles $\theta_{1}$, $\theta_{2}$ and $\theta_{3}$ represent $\hat{h}$–$\hat{H}$, $\hat{s}$–$\hat{h}$ and $\hat{s}$–$\hat{H}$ mixing, respectively. Obviously each matrix element of $\hat{M}^{2}$ can be expressed in terms of the mass eigenvalues $\{m_{h},\,m_{H},\,m_{s}\}$ and the mixing angles $\{\theta_{1},\,\theta_{2},\,\theta_{3}\}$. Among such relations, the following ones are particularly relevant for our subsequent discussion:

[TABLE]

because the constraints on the parameters $\lambda$, $\mu$ 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 $\Delta m^{2}_{12}$ in eq. (2.1) is written

[TABLE]

with $\epsilon$ given by

[TABLE]

The above shows that $\epsilon$ vanishes at $X_{t}=0$ and $X_{t}=\pm\sqrt{6}M_{S}$, and therefore $\Delta 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. Carena:2015moc . It is also easy to see that $|\epsilon|$ is smaller than about $0.1$ for $X^{2}_{t}\lesssim 10M_{S}^{2}$. For stop mixing with $X_{t}$ far from [math] or $\pm\sqrt{6}M_{S}$, the $\epsilon$-contribution to $\Delta m^{2}_{12}$ can be sizable only when $|\mu|$ is large, close to $M_{S}$, and $\tan\beta$ is low. Keeping this feature in mind, we shall neglect the $\epsilon$-contribution in our analysis unless stated otherwise.

2.2 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 Carmi:2012in

[TABLE]

where $(\phi_{1},\,\phi_{2},\,\phi_{3})=(h,\,H,\,s)$ and $f$ denotes the SM fermions.

At tree level, the Higgs couplings to massive SM particles are given by

[TABLE]

The Higgs couplings to gluons and photons are radiatively generated, which results in

[TABLE]

where $\delta C^{i}_{g}$ and $\delta C^{i}_{\gamma}$ are additional contributions from sparticle loops, and the loop functions are given by

[TABLE]

where $\tau^{i}_{j}\equiv m^{2}_{\phi}/(4m^{2}_{j})$ and

[TABLE]

For the case when the superpartners of SM particles are heavier than TeV, sparticles give negligible contributions to $\Delta C^{i}_{g}$. However, if $\mu$ is small, the charged higgsinos are light and can give a sizable contribution to $\Delta C^{i}_{\gamma}$ through the following Higgs-higgsino couplings

[TABLE]

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

[TABLE]

for a Higgs boson with $m^{2}_{\phi_{i}}\ll 4|\mu|^{2}$.

3 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

[TABLE]

while the couplings to gluons and photons read

[TABLE]

Here $r$ is defined by

[TABLE]

and measures the chargino contribution, which has been approximated by using the fact that it is non-negligible only when $|\mu|$ is not far above the electroweak scale for $\lambda$ below the perturbative bound, and the chargino search at LEP requires $|\mu|>104$ GeV. One should note that the couplings of the SM-like Higgs boson are determined by four parameters: $\theta_{1}$, $\theta_{2}$, $r$ and $\tan\beta$.

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 $\phi_{125}$ of the minimal SM with mass 125 GeV:

[TABLE]

Assuming that $h$ does not decay to non-SM particles, one also finds its total decay rate to be

[TABLE]

with $\Gamma_{\rm tot}(\phi_{125})$ being the total decay width of $\phi_{125}$. Here we have used the branching ratios of the SM Higgs boson listed in ref. Heinemeyer:2013tqa . 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

[TABLE]

for the inclusive $WW/ZZ$ channel, where ${\rm Br}(h\to ii)$ is the branching ratio of the indicated mode. For other channels, one obtains

[TABLE]

It is obvious that one should have $\mu^{ii}_{h}=1$ if $\theta_{1}=\theta_{2}=0$ and $r=0$.

The ATLAS collaboration has recently updated the measurements on the Higgs signal strengths using the $13$ TeV data ATLAS-CONF-2019-005 :333 Both ATLAS and CMS collaborations have recently reported their analyses results on the Higgs coupling measurements using the LHC Run 2 data ATLAS-CONF-2019-005 ; Sirunyan:2018koj . The ATLAS analysis used the larger amount of the Run 2 data, up to the integrated luminosity of 80 fb*-1*. As there is no combined global fit for the full Run 2 data yet, we adopted only the ATLAS result in our study.

[TABLE]

and for the other channels

[TABLE]

Note that the NMSSM leads to $\mu^{WW}_{h}=\mu^{ZZ}_{h}$, which is within the $1\sigma$ 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 mixing angles $\theta_{1}$ and $\theta_{2}$ for given values of $r$ and $\tan\beta$. For instance, $\mu^{VV}_{h}=1$ is obtained if $\theta_{1}$ and $\theta_{2}$ satisfy444 Although we will not pursue in this paper, a large value of $\theta_{1}$ satisfying

$\displaystyle\theta_{1}\approx\frac{2.1-1.9s^{2}_{\theta_{2}}}{\tan\beta}$

can also lead to $\mu^{VV}_{h}=1$. In this case, $C^{h}_{b}$ is negative and thus leads to wrong sign Yukawa couplings Choi:2012he , and one needs a large $\lambda$ beyond the perturbativity bound Coyle:2018ydo .

[TABLE]

for which the branching ratio for $h\to VV$ 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 $\lambda$ should be smaller than about 0.7 in order for the NMSSM to remain perturbative up to the GUT scale, and $|\mu|$ should be larger than $104$ GeV to satisfy the LEP bound on the chargino mass.

Figure 1 shows the $2\sigma$ range of $(\theta_{1},\theta_{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\leq\tan\beta\leq 15$ to see how the allowed region changes with $\tan\beta$. From the figure, one can see that the allowed region of $\theta_{1}$ gets smaller if $\tan\beta$ increases, but a broad range of $\theta_{2}$ is allowed insensitively to $\tan\beta$. This is because the Higgs coupling $C^{h}_{b}$ gets sensitive to $\theta_{1}$ at large $\tan\beta$ 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, excluding $\theta_{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, $\delta C^{h}_{\gamma}\approx-0.17r\theta_{2}$ for small mixing angles, whereas the Higgs couplings to other SM particles only quadratically depend on $\theta_{2}$.

4 LEP and CMS excesses around 96 GeV

The CMS collaboration has recently reported a local excess of 2.8$\sigma$ in the diphoton channel around $m_{\gamma\gamma}=96$ GeV Sirunyan:2018aui . The signal strength amounts to

[TABLE]

where $\phi_{96}$ denotes the hypothetical SM Higgs boson with mass $96$ GeV, and ${\rm Br}(\varphi\to\gamma\gamma)$ denotes the branching ratio for the diphoton channel Domingo:2018uim . Intriguingly, there is another 2.3$\sigma$ local excess at the similar mass region from the Higgs searches in the $Z$-boson associated Higgs production ($e^{+}e^{-}\to Z\varphi$) at LEP Barate:2003sz . The signal strength is Cao:2016uwt

[TABLE]

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 singlet-like 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. (14) as follows,

[TABLE]

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 Djouadi:1997yw ; Djouadi:2018xqq 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

[TABLE]

and those to gluons and photons read

[TABLE]

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. (34) is approximated by

[TABLE]

for small $\theta_{2}$, and thus the LEP excess is explained if $s^{2}_{\theta_{2}}\sim 0.1$. On the other hand, the ratio between the CMS and LEP signal rates is given by

[TABLE]

where $k_{32}\equiv s_{\theta_{3}}/s_{\theta_{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_{\theta_{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 $\lambda$ and $\mu$, and $\tan\beta$:

[TABLE]

for $r\equiv\lambda v/|\mu|$, with $m_{h}=125$ GeV and $m_{s}=96$ GeV. On the other hand, the relations (2.1) exhibit how the parameters $\lambda$, $\mu$ and $m_{0}$ change with the mixing angles

[TABLE]

where we have taken $\epsilon=0$ in eq. (12) since it is negligibly small in most of the parameter space of our interest. It is obvious that $m_{H}$ is determined by $\theta_{1,2,3}$ and $\tan\beta$ once we fix $r$. The above relations allow us to analyze the viable mixing angles as follows. We first examine the $(\theta_{1},\,\theta_{2})$ space to see in which region $\mu^{ii}_{h}$ are consistent with the current LHC data, and then continue to check if it is further possible to explain both $\mu_{\rm LEP}$ and $\mu_{\rm CMS}$.

In figure 2, the Higgs signal rates $\mu^{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 $\theta_{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 $\lambda$, $\mu$, and $m_{0}$. Here we have imposed the conditions

[TABLE]

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 mediation models of SUSY breaking.555 The lower bound on $m_{0}$ is set by considering the numerical result that $m_{0}$ is around $113$ GeV for stops at $1$ TeV in the limit of vanishing stop mixing Ellwanger:2004xm ; Ellwanger:2005dv ; Ellwanger:2006rn . Then it follows $r\leq 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\leq\tan\beta\leq 15$.666 We have cross checked the results of our analysis for some parameter points by using NMSSMTools Ellwanger:2004xm ; Ellwanger:2005dv ; Ellwanger:2006rn .

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\beta$ 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 $\Delta m_{12}^{2}$ is to increase (decrease) the coupling $\lambda$ if $\Delta m_{12}^{2}$ is negative (positive), as can be deduced from the last relation in eq. (2.1). This implies that the parameter space compatible with both excesses shrinks for larger negative $\Delta m_{12}^{2}$ due to the perturbativity bound on $\lambda$. In the parameter region with $m_{0}\geq 115$ GeV, however, the dependence on $\Delta m_{12}^{2}$ becomes quite weak because $\lambda$ 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. (41) implies $\lambda>0.6r$, following from $r\equiv\lambda v/|\mu|$. Taking this together with the perturbativity bound on $\lambda$, one can find that $\lambda$ would be more severely constrained at larger $r$ when the stop correction $\Delta m^{2}_{12}$ is negative. We have checked these features by taking analysis for nonzero values of $\epsilon$ 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.1) enable us to extract the information on the region of $\mu$ and $m_{H}$ compatible with the Higgs signal strengths, $\mu^{ii}_{h}$, $\mu_{\rm LEP}$ and $\mu_{\rm CMS}$. Figure 4 shows the allowed region of $(\mu,m_{H})$, where we have taken $\tan\beta=2$ (left) and $5$ (right) with $0<r<1.1$. As discussed already, the $m_{0}$ cut is relevant for small $\tan\beta$. 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}\gtrsim 500$ GeV, nearly irrespectively of the values of $\mu$ and $\tan\beta$, while the upper bound depends on those parameters and is found to increase with $\tan\beta$.

5 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^{-}\to Z(h\to b\bar{b})$ and $pp\to h\to\gamma\gamma$ 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.

Acknowledgements.

We would like to thank S. Heinemeyer for useful comments on the manuscript. This work was supported by IBS under the project code, IBS-R018-D1 (KC and CBP), and by the National Research Foundation of Korea (NRF) grant funded by the Korea government (MSIP) (NRF-2018R1C1B6006061) (SHI and KSJ).