Could two NMSSM Higgs bosons be present near 125 GeV?

We examine next-to-minimal supersymmetric model scenarios with partial parameter unification at the grand unification scale in which both $h_1$ and $h_2$ lie in the 123–128 GeV mass range. Very substantially enhanced $\gamma \gamma$ and other rates are possible. Broadened mass peaks are natural.

Figure 1

Correlation of $gg\to (h_1,h_2)\to \gamma \gamma$ signal strengths when both $h_1$ and $h_2$ lie in the 123–128 GeV mass range. The circular points have $\Omega h^2<0.094$, while diamond points have $0.094\le \Omega h^2\le 0.136$. Points are color coded according to $m_{h_2}-m_{h_1}$ as indicated on the figure.

Figure 2

The first three plots show $R_{gg}^h(X)$ for $X=\gamma \gamma$, $VV$, $b\overline{b}$ versus the appropriate effective Higgs mass $m_h$ as defined in text. The final plot gives $R_{gg}^h(\gamma \gamma )$ versus $R_{W^{*}\to Wh}^h(b\overline{b})=R_{\mathrm{VBF}}^h(b\overline{b})$. The color code here and in the following figures is green for points with $2\mathrm{GeV}<m_{h_2}-m_{h_1}\le 3\mathrm{GeV}$, blue for $1\mathrm{GeV}<m_{h_2}-m_{h_1}\le 2\mathrm{GeV}$, and red for $m_{h_2}-m_{h_1}\le 1\mathrm{GeV}$.

Figure 3

Average light-flavor squark mass, $m_{\tilde{q}}$, versus gluino mass, $m_{\tilde{g}}$, for the points plotted in the previous figures.

Figure 4

Top row: $\Omega h^2$ and spin-independent cross section on protons versus LSP mass for the points plotted in previous figures. Bottom row: $\Omega h^2$ versus LSP higgsino (left) and singlino (right) components.