Revealing mass-degenerate states in Higgs boson signals

TL;DR

This paper explores the possibility that the observed Higgs signals are due to two nearly identical scalar states, analyzing their couplings and how they could be distinguished from a single Higgs in experiments.

Contribution

It introduces a parameterization of Higgs signal strengths in a quasi-degenerate scenario within the NMSSM, aiding experimental distinction from the Standard Model.

Findings

Combined Higgs signal strengths differ from SM predictions in specific decay channels.

Two parameters effectively describe the quasi-degenerate Higgs scenario.

Experimental strategies can distinguish single from double Higgs states.

Abstract

The observed Higgs boson signals to-date could be due to having two quasi-degenerate 125 GeV scalar states in Nature. This kind of scenario tallies well with the predictions from the Next-to-Minimal Supersymmetric Standard Model (NMSSM). We have analysed the phenomenological NMSSM Higgs boson couplings and derived a parameterization of the signal strengths within the two quasi-degenerate framework. With essentially two parameters, it is shown that the combined strengths of the two quasi-degenerate Higgs states in the leptonic (and b-quark) decay channels depart from the Standard Model values in the opposite direction to those in the vector boson channels. We identify experimental measurements for distinguishing a single from a double Higgs scenarios. The proposed parameterization can be used for benchmarking studies towards establishing the status of quasi-degenerate Higgs scenarios.

Tables

Table 1. Table 1: The 26 pNMSSM parameters and their corresponding flat prior probability density distribution ranges. The SM parameters were varied according to Gaussian distributions with the shown central values and standard deviations. The third column (to be addressed in section.  4 ) shows the 95% Bayesian confidence regions for the posterior sample used in Figure.  4 . For this posterior sample, m h 2 − m h 1 < 3 subscript 𝑚 subscript ℎ 2 subscript 𝑚 subscript ℎ 1 3 m_{h_{2}}-m_{h_{1}}<3  GeV with both m h 1 subscript 𝑚 subscript ℎ 1 m_{h_{1}} and m h 2 subscript 𝑚 subscript ℎ 2 m_{h_{2}} allowed within [ 122 122 122 , 128 128 128 ] GeV.

Parameter	Range	Posterior range
$M_1$	[$-4$ TeV, $4$ TeV]
$M_2$	[$0$ TeV, $4$ TeV]
$M_3,m_{{\tilde{f}}_{Q,U,D,L,E}}^{3rdgen,\mathrm{\hspace{0.17em}\hspace{0.17em}1}st/2ndgen}$	[$100$ GeV, $4$ TeV]
$A_{t,b,\tau }$	[-$8$ TeV, $8$ TeV]
$\tan \beta$	[$2$, $60$]	[$8.8$, $28.3$]
$\lambda$	[${10}^{-4}$, $0.75$]	[$0.17$, $0.52$]
$\kappa$	[$-0.75$, $0.75$]	[$-0.50$, $0.75$]
${\mu }_{\mathrm{eff}}$	[$100$, $400$] GeV	[$111$, $308$] GeV
$A_{\lambda }$	[$50$ GeV, $4$ TeV]	[$1.34$, $4$] TeV
$A_{\kappa }$	[$-2$ TeV, $2$ TeV]	[$-1646$, $846$] GeV
$m_t$	172.6 $\pm$ 1.4 GeV
$m_Z$	91.1876 $\pm$ 0.0021 GeV
$m_b{\left(m_b\right)}^{\overline{MS}}$	4.20 $\pm$ 0.07 GeV
${\alpha }_s{\left(m_Z\right)}^{\overline{MS}}$	0.1172 $\pm$ 0.002

Table 2. Table 2: Summary of the central values and errors for the observables. Theoretical uncertainties have been added in quadrature to the experimental uncertainties quoted.

Observable	Constraint	References
$m_h$	$125.09\pm 3.0$ GeV	[94]
$Br\left(B\to X_s\gamma \right)$	$\left(3.32\pm 0.16\right)\times {10}^4$	[95, 96, 97]
$Br\left(B_s\to {\mu }^{+}{\mu }^{-}\right)$	$\left(3.0\pm 0.6\right)\times {10}^{-9}$	[98, 99, 100]
$\Delta M_{B_s}$	$17.757\pm 0.021$	[100, 101]
$\Delta M_{B_d}$	$0.5064\pm 0.0019$	[100, 101]
$Br\left(B_u\to \tau \nu \right)$	$1.06\pm 0.19$	[102, 103, 104, 105]
$\delta a_{\mu }$	$\left(30.2\pm 8.80\right)\times {10}^{-10}$	[106, 28, 29]
${\Omega }_{CDM}{h}^2$	$0.12\pm 0.02$	[107]
Higgs signal strengths		[78, 2, 79, 80, 81, 82, 83, 84, 85, 69, 86, 67, 74, 75, 76, 87, 88, 89, 90, 91, 92, 93, 68]
CDM direct detection limits		[108, 109, 110, 111, 112, 113, 114]
Constraints in HiggsBounds		[63, 64, 65, 66, 67, 68, 69, 70, 71, 72, 73, 74, 75, 76, 77]
Constraints in SModelS		[41, 42, 43, 44, 45, 46, 47, 48, 49, 50, 51, 52, 53, 54, 55, 56, 57, 58, 59, 60, 61, 62]

Table 3. Table 3: Ten parameter fit of μ F f superscript subscript 𝜇 𝐹 𝑓 \mu_{F}^{f} and μ V f superscript subscript 𝜇 𝑉 𝑓 \mu_{V}^{f} . Table 15 of reference [ 1 ]

Parameter	ATLAS + CMS
${\mu }_{V,\gamma \gamma }$	${1.05}_{-0.41}^{+0.44}$
${\mu }_{V,ZZ}$	${0.47}_{-0.92}^{+1.37}$
${\mu }_{V,WW}$	${1.38}_{-0.37}^{+0.41}$
${\mu }_{V,\tau \tau }$	${1.12}_{-0.35}^{+0.37}$
${\mu }_{V,bb}$	${0.65}_{-0.29}^{+0.31}$
${\mu }_{F,\gamma \gamma }$	${1.16}_{-0.24}^{+0.27}$
${\mu }_{F,ZZ}$	${1.42}_{-0.33}^{+0.37}$
${\mu }_{F,WW}$	${0.98}_{-0.20}^{+0.22}$
${\mu }_{F,\tau \tau }$	${1.06}_{-0.56}^{+0.60}$
${\mu }_{F,bb}$	${1.15}_{-0.94}^{+0.99}$

Equations

$${\varphi }_i=N_{ij}{\psi }_j$$

$$\displaystyle N=\left(\begin{array}[h]{ccc}\cos\beta&\sin\beta&0\\ \sin\beta&-\cos\beta&0\\ 0&0&1\end{array}\right).$$

$$h_i=U_{ij}{\varphi }_j$$

$$g_{h_{i}ZZ}$$

$$g_{h_{i}WW}$$

$$g_{h_{i}tt}$$

$$g_{h_{i}bb}$$

$$\mu =\frac{\sigma (SM\to h_i\to SM){|}_{NMSSM}}{\sigma (SM\to h_{SM}\to SM){|}_{SM}}$$

$${\mu }_{\mathrm{VBF}\to h_1\to WW/ZZ}$$

$${\mu }_{\mathrm{VBF}\to h_1\to b\bar{b}}$$

$${\Gamma }_i$$

$$=$$

$${\Gamma }_i/{\Gamma }_{\mathrm{SM}}$$

$$\simeq$$

$${\mu }_{VBF\to h_1\to WW/ZZ}^{\mathrm{an}}$$

$${\mu }_{VBF\to h_2\to WW/ZZ}^{\mathrm{an}}$$

$${\mu }_{VBF\to h_1\to bb}^{\mathrm{an}}$$

$$\text{ and }{\mu }_{VBF\to h_2\to bb}^{\mathrm{an}}$$

$$\underset{m_{H_0}\gg m_{h_0},m_S}{\lim }{\mu }_{VBF\to h_i\to WW/ZZ}^{\mathrm{an}}\simeq (U_{i1}{)}^2,\text{ and }\underset{m_{H_0}\gg m_{h_0},m_S}{\lim }{\mu }_{VBF\to h_i\to bb}^{\mathrm{an}}\simeq (U_{i1}{)}^2.$$

$$\underset{m_{H_0}\gg m_{h_0},m_S}{\lim }{\mu }_{VBF\to h_1\to WW/ZZ}^{\mathrm{an}}+{\mu }_{VBF\to h_2\to WW/ZZ}^{\mathrm{an}}\simeq (U_{11}{)}^2+(U_{21}{)}^2\sim 1\text{ and }$$

$$\underset{m_{H_0}\gg m_{h_0},m_S}{\lim }{\mu }_{VBF\to h_1\to bb}^{\mathrm{an}}+{\mu }_{VBF\to h_2\to bb}^{\mathrm{an}}\simeq (U_{11}{)}^2+(U_{21}{)}^2\sim 1.$$

$$W_{NMSSM}=W_{MSS{M}^{\prime }}-{ϵ}_{ab}\lambda S{H}_1^a{H}_2^b+\frac{1}{3}\kappa S^3,$$

$$W_{MSS{M}^{\prime }}$$

$$L:(1,2,-\frac{1}{2}),\bar{E}:(1,1,1),$$

$$Q:(3,2,\frac{1}{6}),\bar{U}:(\bar{3},1,-\frac{2}{3}),\bar{D}:(\bar{3},1,\frac{1}{3}),$$

$$H_1:(1,2,-\frac{1}{2}),H_2:(1,2,\frac{1}{2}).$$

$$V_{\mathrm{soft}}=V_2+V_3+m_{\mathrm{S}}^2|S{|}^2+(-{ϵ}_{ab}\lambda A_{\lambda }S{H}_1^a{H}_2^b+\frac{1}{3}\kappa A_{\kappa }{S}^3+H.c.),$$

$$V_2$$

$$V_3$$

$$M_{1,2,3};m_{{\tilde{f}}_{Q,U,D,L,E}}^{3rdgen},m_{{\tilde{f}}_{Q,U,D,L,E}}^{1st/2ndgen};A_{t,b,\tau }.$$

$$\tan \beta ,\lambda ,\kappa ,A_{\lambda },A_{\kappa },\lambda v_s.$$

$$\theta =\{M_{1,2,3};m_{{\tilde{f}}_{Q,U,D,L,E}}^{3rdgen},m_{{\tilde{f}}_{Q,U,D,L,E}}^{1st/2ndgen};A_{t,b,\tau ,\lambda ,\kappa };\tan \beta ,\lambda ,\kappa ,{\mu }_{\mathrm{eff}};m_{t,Z,b},{\alpha }_s\}.$$

$$U$$

$$=$$

$$s_{23}=U_{33},$$

$$s_{12}=-\frac{U_{31}}{c_{23}},$$

$$s_{12}\simeq {\theta }_{12},s_{23}\simeq {\theta }_{23},c_{12}\simeq 1,c_{23}\simeq 1$$

$$U$$

Full text

1 Introduction

Higgs boson discovery represents the beginning of a new epoch for fundamental physics. The precise measurements of its couplings is an important aim for particle physics which could possibly give hint to physics beyond the Standard Model. With current data, the Higgs properties are compatible with the prediction of the Standard Model [1, 2]. These same properties could also be due to the combination of effects arising from having two quasi-degenerate scalar states around $125\>\text{GeV}$. Such a tantalizing possibility have been predicted by new physics models such as the Next-to Minimal Supersymmetric Standard Model (NMSSM). The impact of the Higgs properties and precision measurements on the NMSSM scenarios with two quasi-degenerate scalars will contribute towards sharpening our understanding of the Higgs boson data and Nature – it could be that the data might have already contain some indications for new physics.

The current state of findings from the Large Hadron Collider (LHC), i.e. the absence of direct signals of physics beyond the Standard Model (BSM), has been forecasted for the case of supersymmetry (SUSY) by pre-LHC global fits of models to data. For instance, as pointed out in [3, 4, 5] the large mass of the Higgs was already an indication for heavy supersymmetric mass spectra. Within such models, phenomenological studies could be done via two main approaches, namely the simplified models approach [6, 7] and the phenomenological model parameterization [8, 9, 10, 5, 11]. In this article, the latter approach will be used.

Several groups have addressed mass-degenerate Higgs scenarios within the NMSSM. Refs. [12, 13, 14] have considered two quasi-degenerate Higgs states for the real and complex NMSSM, with a mass difference large enough to use the narrow width approximation. Ref. [15] has gone beyond the narrow width approximation and showed that interference effects can account for up to 40% of total cross sections. To be able to conclude that departures from SM prediction are a consequence of the existence of more than one resonance [16, 17] have proposed statistical test based on the analysis of a signal strength matrix, where all the channels are considered independent. A simplified version of their results agrees with what was proposed previously in [12]. In this article, we focus on the possibility of having two mass-degenerate states with different coupling structures that when combined mimic a single Higgs features. The main aim is to derive a set of NMSSM parameters most relevant for quasi-degenerate Higgs studies vis-á-vise collider data. For this, the NMSSM doublet-singlet mixings structure [18, 19, 15] of the Higgs sector will be used.

In section 2 we review the production and decay ratios of the two lightest NMSSM CP-even Higgs states. We focus on the couplings of these to vector bosons and heavy quarks. In section 3 we perform a scan of the parameters of the NMSSM while imposing that the the two lightest CP-even Higgs states reproduce the mass of the standard Higgs measured by the LHC. We describe the allowed parameter space regions and relevant parameter correlations. In section 4 the sample is then used together with analytical relations for the couplings and signal strengths to show that the the quasi-degenerate Higgs properties can be explained approximately by using just two free parameters. We also we show how the superposition of two quasi-degenerate Higgs around 125 GeV could be in agreement with current experimental results. Finally in section 5 we analyse the sample based on signal strength ratios that can discriminate between the single versus double resonance scenarios.

2 Higgs couplings to fermions and vector bosons

Right after the discovery of the Higgs the search for signals of physics beyond Standard Model in the production and decay of the Higgs became a priority. A possible excess in the $\gamma\gamma$ channel motivated a lot of work, some of them within the NMSSM framework [20, 21, 12, 19, 22]. In particular King et. al. [19] pointed out that the signal strengths of the $\gamma\gamma$ channels could be enhanced for large singlet-double mixing. We will take these as a starting point for analysing two quasi-degenerate CP-even Higgs states.

For the discussion of the following sections it is important to have a clear picture of how the widths and therefore the Higgs branching ratios depend on the singlet-doublet mixing. Let us start introducing some notation, we define $\psi=(H_{d},H_{u},s)$ and $\phi=(h_{0},H_{0},s)$ in such a way that $\langle h_{0}\rangle=v$ and $\langle H_{0}\rangle=0$:

[TABLE]

where

[TABLE]

The Higgs states $h=(h_{1},h_{2},h_{3})$ are related to $\psi$ and $\phi$ in the following way,

[TABLE]

where $U_{ij}$ are the elements of the mixing matrix, $U$. We consider it convenient to use the elements of $U$ to parameterise the couplings; for example $U_{i1}$ and $U_{i2}$ are respectively the $h_{0}$-component and $H_{0}$-component of $h_{i}$. In this way it is easier to make the comparison to the standard Higgs.

Using the above notation we write the tree-level Higgs couplings to vector bosons and heavy quarks as:

[TABLE]

In the $H_{0}$ decoupling limit (i.e. $U_{12}=U_{22}=0$) all the couplings are proportional to $U_{11}$, the $h_{0}$-component of $h_{1}$. We are interested in the departure of the production and decay signals of $h_{1}$ in the $Z_{3}$-invariant NMSSM with respect to the one of the standard Higgs. To weight this we will use the signal strength,

[TABLE]

Because of the small width of the Higgs states we assume they are produced on-shell, therefore the total cross sections are evaluated as the production cross section times the branching ratio.

Now, in order to obtain the required properties for the Higgs states to reproduce ATLAS and CMS measurements we consider two possibilities:

1. I)

$h_{1}$ or $h_{2}$ is the Higgs state detected at the LHC, and 2. II)

$h_{1}$ and $h_{2}$ are the Higgs states measured by the LHC, where $h_{1}$ and $h_{2}$ are mass degenerate.

We will show that these two possibilities correspond, respectively, to:

1. I)

Small singlet-doublet mixing, and 2. II)

Large singlet-doublet mixing.

Let us analyse the case with small singlet-doublet mixing where $h_{1}$ is mainly $h_{0}$, in other words $U_{11}\sim 1$. For this case it is a good approximation to consider that the width of $h_{1}$ is dominated by the decay rate of $h_{1}\rightarrow b\bar{b}$ and therefore the variation of the width is controlled by the square of the Higgs coupling to bottom quarks, $g_{h_{1}b\bar{b}}$. Using the couplings described in eq. (2) the signal strengths of the vector-boson fusion production of $h_{1}$ and further decay to $WW/ZZ$ and $b\bar{b}$ are approximately,

[TABLE]

where $\hat{g}=g_{\mathrm{NMSSM}}/g_{\mathrm{SM}}$, the couplings $g_{\mathrm{NMSSM}}$ are those in eq. (2), and $g_{\mathrm{SM}}$ are the Standard Model (SM) couplings. The enhancement or suppression of the first signal strengths depends on $\tan\beta\,\,U_{12}/U_{11}$. As such, the absolute value and sign of this factor determines respectively the magnitude of the ratio between the signal strengths and whether there is an enhancement or suppression of $\mu_{\mathrm{VBF}\rightarrow h_{1}\rightarrow\mathrm{WW/ZZ}}$ with respect to $\mu_{\mathrm{VBF}\rightarrow h_{1}\rightarrow b\bar{b}}$. A similar analysis holds when $h_{2}$ is considered the Higgs state measured at the LHC.

Next, let us examine the case with large singlet-doublet mixing where $h_{1}$ has non-negligible S content. In this case, the approximation $U_{11}\sim 1$ is not valid any more. The assumption that the width of $h_{1}$ is almost totally controlled by $h_{1}\rightarrow b\bar{b}$ is no longer a good approximation. The size of $\tan\beta\,U_{12}/U_{11}$ may take very large values and therefore the branching ratio could significantly differ with respect to the standard Higgs. So, we would like to have a simple expression for the widths appropriate for all values of $U_{i1}$. In terms of the standard Higgs decay rates, one can write

[TABLE]

where $h_{i}\rightarrow\mathrm{SM}_{\mathrm{rest}}$ represents the rest of the decay channels. The dominant contribution for the rest of decay channels is the decay to gluons through a top loop. For simplicity we are going to consider that the rest of the decay modes behave as the ones of the standard Higgs. For this reason we took the corresponding decay rate proportional to the square of $h_{i}$’s $h_{0}$ content, $U_{i1}$. By writing the decay rates in terms of the SM branching ratios we get

[TABLE]

For large singlet-doublet mixing the widths of $h_{1}$ and $h_{2}$ could be much smaller than $\Gamma_{\mathrm{SM}}$, producing large departures of the branching ratios with respect to the ones of the standard Higgs, unless the widths and the decay rates of each Higgs state change at the same proportion. From now on we will use eq. (13) as the enhancement(suppression) rate of the width with respect to the SM value.

The analytic expressions for the signal strengths for vector-boson fusion production and decay to $WW/ZZ$ and $b\bar{b}$ can be written as,

[TABLE]

Note that for a large singlet-doublet mixing the relative size of $\tan\beta\,U_{12}/U_{11}$ has a larger range of variation than in the case of small singlet-doublet mixing, as consequence there might be larger enhancement(suppression) to the signals. Moreover, since the $H_{0}$-component of the Higgs states is the one responsible for large variations of the branching ratios, it is interesting to see that in the $H_{0}$ decoupling limit ($U_{12}\simeq 0$ and $U_{22}\simeq 0$),

[TABLE]

Hence for large singlet-doublet mixing it is not possible to reproduce the experimental data with a single Higgs state. But, if $h_{1}$ and $h_{2}$ are mass quasi-degenerate, assumed to be unresolved away from each other by experiments, the superposition of the two states could show up in signals as single standard Higgs with,

[TABLE]

Notice that the last (approximate)equalities require $U_{31}\simeq 0$ to fulfill the unitarity condition for U.

It is interesting to compare the departure of the signal strengths for different channels of the same Higgs state. As described earlier, the ratio between signal strengths depends on $\tan\beta\,U_{12}/U_{11}$ for $h_{1}$ and on $\tan\beta\,U_{22}/U_{21}$ for $h_{2}$. As such, the departure of the global signal strength will depend on the relation between $U_{12}$ and $U_{22}$.

In the following sections we analyse the scenario with large singlet-doublet mixing. We will assume that the Higgs signal measured by ATLAS and CMS is a superposition of the production and decay of two Higgs states. To get the global enhancement(suppression) we will sum the contribution of the two Higgs states. Notice that for this approximation to be valid the widths should be much smaller that the mass difference between $h_{2}$ and $h_{1}$.

3 The phenomenological NMSSM Parameters scan

Let us consider the case where the Higgs signal measured by ATLAS and CMS is a superposition of the production and decay of $h_{1}$ and $h_{2}$, meaning that the Higgs states are close enough not to be resolved by the experiments, but with large enough separation to have negligible interference effects. To study the region of the parameter space of the NMSSM where this condition is fulfilled we perform a parameter scan as done in [23].

3.1 The phenomenological NMSSM (pNMSSM)

We shall consider an R-parity conserving NMSSM with superpotential,

[TABLE]

where

[TABLE]

The chiral superfields have the following $SU(3)_{C}\otimes SU(2)_{L}\otimes U(1)_{Y}$ quantum numbers,

[TABLE]

The corresponding soft SUSY-breaking terms are

[TABLE]

with

[TABLE]

A tilde-sign over the superfield symbol represents the scalar component. However, an asterisk over the superfields as in, for example, $\tilde{u}_{R}^{*}$ represents the scalar component of $\bar{U}$. The $SU(2)_{L}$ fundamental representation indices are donated by $a,b=1,2$ while the generation indices by $i,j=1,2,3$. $\epsilon_{12}=\epsilon^{12}=1$ is a totally antisymmetric tensor.

In an approach similar to that of the pMSSM [8, 9, 10, 5], the pNMSSM parameters are defined at the weak scale with the non-Higgs sector set,

[TABLE]

Here, $M_{1,2,3}$ and $m_{\tilde{f}}$ are respectively the gaugino and the sfermion mass parameters. $A_{t,b,\tau}$ represent the trilinear scalar couplings. With electroweak symmetry breaking,an effective $\mu$-term, $\mu_{\mathrm{eff}}=\lambda\,v_{s}$ is developed. The $\mu$-term, the ratio of the MSSM-like Higgs doublets’ vevs $\tan\beta=\left<H_{2}\right>/\left<H_{1}\right>$ and the Z-boson mass, $m_{Z}$ lead to the tree-level Higgs sector parameters

[TABLE]

Next, including four SM nuisance parameters, namely, the top and bottom quarks $m_{t,b}$, $m_{Z}$ and the strong coupling constant, $\alpha_{s}$, makes the pNMSSM parameters:

[TABLE]

3.2 The scanning procedure

$M_{1,2}$ affects the gaugino masses for which a wide range, $\mathcal{O}(\mathrm{GeV})$ to $\mathcal{O}(\mathrm{TeV})$, is possible. We let $M_{1}\in[-4,4]\>\text{TeV}$ and same for $M_{2}>0$. With the LHC in mind, we let the gluino and squark mass parameters be within $[100\>\text{GeV},4\>\text{TeV}]$, and the trilinear scalar couplings allowed in $[-8\>\text{TeV},8\>\text{TeV}]$. $\tan\beta$ is allowed between 2 and 60. For minimising fine-tuning, we subjectively let $\mu_{\mathrm{eff}}=\lambda\,v_{s}$ to vary within 100 to 400 GeV not too far away from the Z-boson mass. The remaining Higgs-sector parameters were set within the ranges shown in Table 1.

The selected pNMSSM points were required pass all the constraints summarised in Tab.2. These are: the Higgs boson mass $m_{h}$, the neutralino cold dark matter (CDM) relic density $\Omega_{CDM}h^{2}$, anomalous magnetic moment of the muon $\delta a_{\mu}$, and the B-physics related limits summarised in the upper part of Table 2. The experimental constraints used were those implemented in NMSSMTools [24, 25, 26, 27, 28, 29], Lilith [30], MicrOMEGAs [31, 32, 33, 34, 35, 36, 37, 38, 39, 40], SModelS’[41, 42, 43, 44, 45, 46, 47, 48, 49, 50, 51] implementation of ATLAS and CMS limits[52, 53, 54, 55, 56, 57, 58, 59, 60, 61, 62], and HiggsBounds [63, 64, 65, 66, 67, 68, 69, 70, 71, 72, 73, 74, 75, 76, 77]. The Higgs boson signal strength measurements from Tevatron [78], ATLAS [2, 79, 80, 81, 82, 83, 84, 85, 69, 86, 67] and CMS [74, 75, 76, 87, 88, 89, 90, 91, 92, 93, 68] as implemented in Lilith v1.1 (with data version 15.09) [30] were also included.

3.3 Constraints on the parameters of the Higgs

sector

From the pNMSSM parameter scan, we use a sample with two quasi-degenerate lightest CP-even Higgs bosons. It was required that $h_{1}$ and $h_{2}$ have mass equal to $125\pm 3$ GeV, where the $\pm 3$ GeV accounts to the theoretical errors associated to the values of the masses computed by NMSSMtools. In addition it was required that the mass difference, $m_{h_{2}}-m_{h_{1}}<3$ GeV 111 The CMS resolutions for Higgs bosons are channel dependent and typically around 2.5 to 4 GeV [74, 75] for bosonic channels. As such $m_{h_{2}}-m_{h_{1}}<3$ GeV can be considered as a mass degeneracy condition for which the two Higgs cannot be resolved by CMS run-2.. We focus on the regions of the Higgs sector parameters for studying the correlations within those parameters and for relating them to other parameters which are directly connected with the signals measured at the LHC such as the CP-even Higgs mixing matrices.

It is useful to have an explicit form for the Higgs mixing matrix $U$. We parameterise this using three angles $\theta_{13}$, $\theta_{12}$, and $\theta_{23}$ such that

[TABLE]

Here $c_{ij}=\cos\theta_{ij}$ and $s_{ij}=\sin\theta_{ij}$. Given the mixing matrix, obtained numerically by the SUSY spectra calculator NMSSMtools, then the mixing angles can be extracted as:

[TABLE]

Now, considering that we want to reproduce a standard Higgs signal, we determine the expected ranges for the mixing angles. In order to get the ratio between $\mu_{VBF\rightarrow h_{1,2}\rightarrow WW/ZZ}$ and $\mu_{VBF\rightarrow h_{1,2}\rightarrow bb}$ close to one, either the value of $\mu_{VBF\rightarrow h_{i}\rightarrow WW/ZZ}/\mu_{VBF\rightarrow h_{i}\rightarrow bb}$ for each Higgs state has to be close to one, or a fine cancellation should take place. In this work we focus on the first case222In other words, this means that we restrict our analyses to the scenario where $H_{0}$ is much heavier than $h_{0}$ and $S$.. From eqs. (15)-(18) one can see that this condition is possible when $U_{12}$ and $U_{22}$ are very small and as a result $s_{12}$ and $s_{23}$ should also be very small according to eq. (41). On an other hand, eq. (19) implies that the superposition of $h_{1}$ and $h_{2}$ can reproduce the standard Higgs signal for $U_{31}\sim 0$ (i.e. large values of $m_{H_{0}}$). For this to happen either $\theta_{12}$ has to be very small or $\theta_{23}$ has to be close to $\pm\pi/2$. In summary, $\theta_{12}\sim 0$ and $\theta_{23}\sim 0$ will guarantee that we are working in the regime where the superposition of the two Higgs states agrees with experimental measurements.

In the limit of small $\theta_{12}$ and $\theta_{23}$,

[TABLE]

and the mixing matrix eq. (41) reduces to

[TABLE]

where we have neglected $\mathcal{O}(\theta^{2})$ terms. For the results of our scan this approximation works with a $0.5\%$ error.

We have been able to constrain the parameters of the mixing matrix requiring conditions that will give us a standard-like Higgs signal. This conditions will affect the masses or couplings of the heaviest and pseudoscalar Higgs bosons. To see this, it will be useful to relate the mixing angles $\theta_{13}$, $\theta_{23}$ and $\theta_{12}$ to the fundamental parameters of the Higgs sector. Using eq. (41) we relate the terms of the mass matrix with the physical masses by introducing two new parameters: $m_{h}^{2}$, the central value of the two lightest CP-even Higgs states, and $\delta m_{h}^{2}$, half of the squared mass difference,

[TABLE]

To simplify the expressions obtained from eq. (52) we factorise $c_{12}$ and $c_{23}$ to write U in terms of $t_{kl}\equiv\tan\theta_{kl}$ and use the approximations:

[TABLE]

where $kl=12,23$ and $mn=12,23$. Finally, we will focus on the relations in terms of the mass matrix elements $\mathcal{M}_{22}$ and $\mathcal{M}_{23}$ since $\mathcal{M}_{22}^{\mathrm{tree}}$ and $\mathcal{M}_{23}^{\mathrm{tree}}$ reproduce pretty well the values computed by NMSSMtools, and because we wish to get simple relations between the Higgs sector parameters, masses and mixing angles. We have checked numerically that for the rest of mass matrix elements the tree level expression are not precise enough.

[TABLE]

We can further simplify eq. (54) taking into account that $\delta m_{h}^{2}$ and $m_{h}^{2}$ are smaller than $m_{h_{3}}$ and $\mathcal{M}_{22}$. Using the last approximation of eq. (53) we get that terms proportional to $t_{12}^{2}$, $t_{23}^{2}$ and $t_{12}t_{23}$ in the right hand of eq. (54) are negligible. Regarding eq. (56), using the approximation eq. (53) and eq. (55) one gets $\mathcal{M}_{23}+\delta m_{h}^{2}\,t_{12}\sin(2\theta_{13})\simeq\mathcal{M}_{23}$, allowing us to neglect the term proportional to $t_{12}$ in eq. (56) (besides that, for the sample of pNMSSM points described in section 3 the values of $\theta_{12}$ are much smaller than the values of $\theta_{23}$). Hence eqs. (54)-(56) can be rearranged to get,

[TABLE]

where in the last equation we have further considered that $\delta m_{h}^{2}\cos(2\theta)\ll m_{h}^{2}$.

Using the approximation of large $\tan\beta$ and large $M_{A}$ from reference [115]:

[TABLE]

We have checked numerically that $m_{h_{3}}\approx M_{A}$ is a good approximation for the pNMSSM points considered. Now, let us take $\mathcal{M}_{23}$ from reference [115]333Since they perform a different rotation, written in eq. 16 of [115], we transform the mass matrix as follow:

\displaystyle\mathcal{M}=\left(\begin{array}[h]{ccc}M_{22}&-M_{21}&M_{23}\\ -M_{12}&M_{11}&-M_{13}\\ M_{32}&-M_{31}&M_{33}\end{array}\right)

(63)

[TABLE]

and replace it in eq. (59), considering that $M_{A}$ is much heavy than $m_{h}^{2}$ one can write $t_{23}$ as,

[TABLE]

where $v_{s}=\sqrt{2}\,\mu/\lambda$ and $\xi=\mu/M_{A}$.

Left panel of Figure 1 shows in the x-axes the value of $\theta_{23}$ computed by NMSSMtools and in the y-axes the analytical approximation described in eq. (64), as one can see in the Figure there is a good agreement between the analytical expression and the numerical value (green points), and it is clear that the main contribution to $\theta_{23}$ comes from the first term of eq. (64) (blue points). Right panel of Figure 1 shows the relation between $\theta_{23}$ and $m_{h_{3}}$ for constant values of $\lambda$. There is a trend: larger values of $|\theta_{23}|$ correspond to smaller values of $m_{h_{3}}$, except for very small values of $|\theta_{23}|$ where the two parameters seem to be uncorrelated. Still, eq. (64) shows that the value of $\tan{\theta_{23}}$ is not directly related to the scale of the heaviest Higgs, but instead it is related to the value of $\lambda$, $\mu$ and $\tan\beta$ 444Let us remember that in the decoupling limit of $H_{0}$,

$m_{h_{3}}^{2}\simeq M_{A}^{2}=\frac{2\mu}{\sin{2\beta}}\left(A_{\lambda}+\frac{\kappa}{\lambda}\mu\right)$

.

Although the Higgs boson masses get important contributions from loop corrections, it is possible to get some information from the tree level expressions for $m_{h_{1}}$ and $m_{h_{2}}$. For large values of $\tan\beta$ and $M_{A}^{2}$,

[TABLE]

where $v_{S}=\sqrt{2}\mu/\lambda$ (see Eq. (32) of [115]). In order to get a constrain for the initial parameters from the condition of small mass difference between the two lightest Higgs states, we require a small mass difference between the tree level masses showed in eq. (65). But, since the tree level expression do not precisely reproduce the masses of the Higgs states we request the mass square difference at tree level to be smaller than $M_{Z}^{2}$, meaning that both terms inside the square root should be smaller than $M_{Z}^{4}$.

Let us focus on the first term, for $A_{\kappa}\gg M_{Z}$ there should be a correlation between $A_{\kappa}$ and $\kappa v_{s}$ such that there is a cancellation that leads to an order $M_{Z}^{2}$ value. Note that the average of the tree-level squared masses also requires this cancellation to occur in order to get the masses of the Higgs states in the desired range.

For $|A_{\kappa}|\gg M_{Z}$ we expect,

[TABLE]

Figure 2 shows the relation between $A_{\kappa}$ and $\kappa v_{s}$, as manifested in the figure for $|A_{\kappa}|\gtrsim 600$ GeV the approximation of eq. (66) works within an error smaller than $5\%$.

Furthermore, using eq. (66) it is possible to simplify other parameters relevant in the Higgs sector, eq. (30) of [115] gives a simplified expression for the mass of the light pseudoscalar,

[TABLE]

Putting eq. (66) into eq. (67) we write the mass of the lightest pseudoscalar in terms of $\kappa$ and $v_{s}$,

[TABLE]

Figure 3 shows the comparison between eq. (68) and the value computed by NMSSMtools. It can be seen that for $m_{A_{1}}>500$ GeV eq. (68) is a pretty good representation for the light pseudoscalar mass.

For completeness, it is worth mentioning that the second term inside the squared root of eq. (65) is suppressed by a factor $v_{s}^{\,-2}$, as such we do not expect to get any good correlation of parameters from there.

All the information, presented above, are useful for determining an optimal range of parameters in order to perform a specialised parameters scan dedicated for studying mass-degenerate Higgs region(s).

4 The two lightest CP-even Higgses at the LHC

In this section we will use the results of the scan and the analytical relations for the couplings and signal strengths to study the parameter space where the two lightest CP-even Higgs states mimic the SM-Higgs signals.

First, we have to verify the validity of the analytic expressions for the signal strengths comparing these expressions with the numerical values computed by NMSSMtools.555To perform this comparison we flip the order of the mass eigenstates computed by NMSSMtools, in such a way that $h_{1}$ has the largest component of $h_{0}$, and it is not necessary the lightest mass eigenstate. The need of this transformation is due to the convention used for the Higgs mixing matrix in NMSSMtools. The determinant of this matrix could be positive or negative depending on $h_{0}$-fraction of the lightest eigenstate. It is positive if $h_{1}$ is $h_{0}$-dominated and negative if it is $S$-dominated.

The reason why we perform the flip of states is because we want to make a comparison of the analytic relations as function of the mixing angles, for this we need to assume a specific form of the mixing matrix U.

Figure 4 shows the comparison between the signal strengths computed by NMSSMtools, $\mu^{\mathrm{num}}$, and the analytic approximations showed in eqs. (15)-(18), $\mu^{\mathrm{an}}$, for VBF$\rightarrow h_{i}\rightarrow WW/ZZ$ (left panel) and VBF$\rightarrow h_{i}\rightarrow bb$ (right panel). From the figure we see that there is a good agreement between the analytical approximation and the numerical computation.

Now, let us identify the relevant parameters that produce deviation from experimental measurement. Writing the couplings, widths and signal strengths in terms of the mixing angles, for small values of $\theta_{12}$ and $\theta_{23}$, see eqs. (2) and (51),

[TABLE]

Using eq. (77) and eq. (13) we get,

[TABLE]

Finally, eq. (15)-(16) can be written in terms of the mixing angles as

[TABLE]

From eqs.(80)-(83) we see that the signal strengths depend on four parameters: $\theta_{13}$, $\theta_{23}$, $\theta_{12}$ and $\tan\beta$. However, in the limit where $\theta_{12}\tan\beta\ll\theta_{13}$, which is the case for the pNMSSM posterior sample analysed, the number of parameters reduces to two:

[TABLE]

From eqs. (80)-(83), one can see that the dependence on $\theta_{12}$ always appears as a factor in the expression $\cos{\theta_{13}}(1+\theta_{12}\tan\beta)$ or $\sin{\theta_{13}}(1+\theta_{12}\tan\beta)$. Therefore for $\theta_{12}\tan\beta\ll 1$ the contribution of $\theta_{12}$ is negligible.

To understand the dependence of the signal strengths with respect to $\theta_{13}$ and $\theta_{23}\tan\beta$ let us start analysing the relation between the signal strengths for a given Higgs state. The top row of Figure 5 shows the correlations between $\mu_{VBF\rightarrow h_{i}\rightarrow WW/ZZ}$ and $\mu_{VBF\rightarrow h_{i}\rightarrow bb}$ for $h_{1}$ (top left) and $h_{2}$ (top right); for $h_{1}$ we can see that the difference between the $b\bar{b}$ and $WW/ZZ$ channel signal strengths is not small. In fact, this could be taken to imply that it is not possible to reproduce the experimental results with such differences. However, looking at the right panel of the Figure and using the colour code to select regions with constant values of $\theta_{13}$, it is possible to compare the rates of the signal strengths for both Higgs bosons. The plots show that the enhancement(suppression) of one channels of $h_{1}$ is more or less compensated with a suppression(enhancement) in the same channel of $h_{2}$.

The analytic expressions for the widths of the Higgs states, eqs. (78) and (79), show that the term proportional to $\theta_{23}\tan\beta$ has a minus sign in the width of $h_{1}$ and plus sign in the width of $h_{2}$, decreasing(increasing) the decay rate of $h_{1}\rightarrow bb$ while increasing(decreasing) the decay rate of $h_{2}\rightarrow bb$ as $|\theta_{23}|$ increases its value.

The bottom row of Figure 5 shows the width of $h_{1}$ and $h_{2}$ as function of $\theta_{13}$ and $\theta_{23}\tan\beta$. The figure agrees with what we expected from the approximate expressions, eqs.(78) and (79), a function dominated by $\cos^{2}{\theta_{13}}$ for $h_{1}$ and $\sin^{2}{\theta_{13}}$ for $h_{2}$, the phase of the distributions varies with the values of $|\theta_{23}\tan\beta|$.

Let us analyse the global signal strengths. Figure 6 shows the sum of the signal strengths of vector-boson fusion production and decay to $WW/ZZ$ (left panel) and to $bb$ (right panel), these factors represent the global enhancement or suppression of the superposition of the two signals respect to the signal of the standard Higgs. It is important to keep in mind that to get the global signal strengths we sum the contributions of the individual signal strengths, which is allowed since we require the mass difference of the two lightest CP-even Higgs states to be small enough not to be resolved by current experiments, but much larger than the width of the particles to neglect interference effects.666With this approach we are not considering the shape of the signal distribution. The analysis of the shape of the distribution goes beyond the scope of this work

There are several points we would like to comment from Figure 6, the departure of the signal strength increases with the size of $\theta_{23}\tan\beta$ as in the case of the individual signal strengths. The modification of the signal strengths for $h_{1}$ is “compensated” by the modifications of the signal strengths for $h_{2}$ and therefore the total effect is smaller than the one for the individual rates but still not negligible. Regarding the relation between the two global signals strengths it is clear from Figure 6 that $\mu_{VBF\rightarrow h_{1,2}\rightarrow WW/ZZ}$ has opposite behaviour and larger range with respect to $\mu_{VBF\rightarrow h_{1,2}bb}$.

There are two regions that seem to be in full agreement with the SM (the signal strength is $\simeq 1$): the region where $\theta_{23}\simeq 0$ and the region where $\theta_{13}\simeq 0$, as we expected. There is a third region where $\theta_{13}$ is between 0.2 and 0.4, where for a very precise value of $\theta_{23}$ the signal strength is very close to one. On the other hand, for small values of $\theta_{23}$, let’s say $\theta_{23}\tan\beta\gtrsim-0.25$, the deviation from one of the signal strength is very small, very precise measurements will be necessary to resolve it.

There is one last comment about Figures 5 and 6. We are able to fully describe the rates and the widths of $h_{1}$ and $h_{2}$ in terms of two parameters: $\theta_{13}$ and $\theta_{23}\tan\beta$, instead of three, indicating that $\theta_{12}\tan\beta\ll 1$ for the set of successful scanned points.

So far we have focused our study to two channels: VBF$\rightarrow h_{i}\rightarrow WW/ZZ$ and VBF$\rightarrow h_{i}\rightarrow bb$, but the current measurements of the Higgs couplings constrain several more channels. Let us comment about the most relevant production and decays:

1. a)

Production processes like gluon-gluon fusion (GGF) and Higgs production associated to top quarks (ttH) are very important. To analyse these let us go back to eqs.(2), which describe the couplings of the Higgs states to top quarks,

[TABLE]

Comparing $\hat{g}_{g_{i}tt}$ with $\hat{g}_{g_{i}bb}$ we see that the contribution from $\theta_{23}$ is $\cot^{2}\beta$ times smaller for $\hat{g}_{h_{i}tt}$ than for $\hat{g}_{h_{i}bb}$, therefore we expect the contribution of $\theta_{23}$ to be very tiny and the production processes of GGF and ttH to behave as vector-boson fusion for given values of $\theta_{13}$ and $\theta_{23}\tan\beta$. 2. b)

The Higgs decay to photons was one of the most important channels for the discovery of a new particle, where the main contribution to the decay of the standard Higgs to photons is through a loop of W bosons. We expect that the decay of the Higgs states to photons with respect to the value of the standard Higgs scale as the decay to WW/ZZ. 3. c)

The decay of the Higgs states to taus with respect to the value of the standard Higgs will scale as the decay of the Higgs states to bottom quarks.

To complete the description of the signals of the two lightest CP-even Higgs states, in Figure 7 we show the signal strengths for GGF$\rightarrow h_{1,2}\rightarrow WW/ZZ$ (left panel) and GGF$\rightarrow h_{1,2}\rightarrow\gamma\gamma$ (right panel). As we expected, the gluon-gluon fusion production of the Higgses and decay to WW/ZZ is pretty similar to the vector-boson fusion production, on the other hand, the decay to photons shows a larger departure.

So far we have seen that the leading behaviour of the signal strengths is given by $\theta_{13}$ and $\theta_{23}\tan\beta$. In the limit where $\theta_{12}\simeq 0$, we write a biunivocal function to determine one (of these parameters) in terms of the other. An approximate relation between $\theta_{13}$ and $\theta_{23}\tan\beta$ might be useful to study the region around $0.2\lesssim\theta_{13}\lesssim 0.4$ where it seems possible to mimic the signal of the standard Higgs and make it indistinguishable even for very precise experimental measurements. To determine the relation between the parameters we choose the to solve the equation:

[TABLE]

By taking $\mu_{VBF\rightarrow h_{1}\rightarrow WW/ZZ}$ and $\mu_{VBF\rightarrow h_{2}\rightarrow WW/ZZ}$ from eqs. (80) and (81), neglecting the terms proportional $\theta_{12}$, and rewriting the $\sin\theta_{13}$ and $\cos(\theta_{13})$ in terms of $\sin(2\,\theta_{13})$ and $\cos(2\,\theta_{13})$ we can simplify eq. (84) to get a quadratic equation in $\cot(2\,\theta_{13})$. So, there are two solutions for $\theta_{13}$:

[TABLE]

where $\overline{\mathrm{BR}}_{bb}=1-\mathrm{BR}_{bb}$. For $\delta=0$ the solution simplifies to

[TABLE]

With eq. (85) we are able to determine $\theta_{13}$ in terms of $\theta_{23}\tan\beta$ and $\delta$. Figure 8 shows the comparison between the semi-analytical relation in eq. (85) and the numerical results from our scans. Although it is not a precise relation, eq. (85) gives a very good approximation to the correlation between $\theta_{13}$ and $\theta_{23}$ for a fixed value of $\delta$.

5 Searching for mass-degenerate Higgses

As commented in references [12] and [16] there are ways to test the existence of mass-degenerate states. The determinant of a signal strengths square matrix could give information about the number of resonances. If the determinant of the square matrix is equal to zero then the existence of a single Higgs resonance will be enough to reproduce the signal strengths.

For simplicity we will use a compact notation: $\mu_{ij}=\mu_{i\rightarrow j}$, where $i$ represents the production mode and $j$ the decay channel. Considering two square matrices,

[TABLE]

the condition for the determinant to be non-zero can be written in terms of the ratios

[TABLE]

To check if it is possible to establish the existence of two resonances in the NMSSM we consider the set of pNMSSM posterior sample described in section 3 and check for points which are within one and three sigma of the particular signal strengths listed in Table 3.

Figure 9 shows the comparison between the ratios of the signal strengths in eq. (92). The upper (lower) panel shows all the points that are within three (one) sigma of the values of the individual rates. The points are ordered in such a way that smaller values of $|\theta_{23}\tan\beta|$ are on top. Notice that in the lower panel the one sigma region do not contain the point $\{1,1\}$, which is what we expect from a standard Higgs, this is because the experimental value of $\mu_{VBF,bb}$ is $0.65^{+0.31}_{-0.29}$ (see Table 3), it doesn’t include the SM value at one sigma. The left panel of Figure 9 shows that the ratios between $WW$ and $\gamma\gamma$ signal strength are basically the same, meaning that the determinant of $R_{A}$ is approximately zero and therefore in agreement with a single resonance hypothesis. On the other hand the ratios between $\tau\tau$ and $\gamma\gamma$ signal strength are slightly separated from the dotted line, the determinant of $R_{B}$ is different from zero. In general we would expect that if there is more than one Higgs state the ratio between two signal strengths with the same production process and different decay product is not going to be equal to one. However, we get that this ratio is almost the same for the rate between gluon-gluon fusion and for vector-boson fusion production processes, which indicates that both production cross-sections are very similar for a given Higgs state. Therefore, it doesn’t seem possible to distinguish between single and double resonances from those measurements for this set of scanned points.

Is there any observable that could be used to distinguish between single and double resonance signals? From the discussion of the previous sections we have learned that $\mu_{VBF,bb}$ have an opposite behaviour with respect to the other signal strength we have considered, therefore we may suspect that the production of Higgs states associated to bottom quarks compared to the production associated to vector bosons would give a larger departure from the SM signal than the comparison between vector-boson fusion and gluon-gluon fusion.

Let us consider the matrices,

[TABLE]

where BBF represents the Higgs productions associated to bottom quarks. To obtain a determinant different from zero requires that ratios of the signal strengths follow:

[TABLE]

To compute the signal strength of Higgs production associated to bottom quarks we use the reduced couplings to bottom quarks computed by NMSSMtools.

Figure 10 shows the comparison of the ratios described in eq. (10) for points that fulfill the experimental signal strength listed in Table 3 within three sigma. The figure shows that the determinant of the $R^{C}$ and $R^{D}$ is different from zero for a large part of the points, and therefore it gives a clear signature for the existence of more than one Higgs resonance.

It may be surprising to see such a large deviation from zero in the determinant of $R^{C}$ and $R^{D}$ and not in the determinant of $R^{A}$ and $R^{B}$, the main reason lies in the difference between the production processes. Although it does not seem straight forward from the analytic expressions of the full signal strength to single out this differences and directly relate them with the value of the determinants, one can always compare the production cross-sections for each Higgs state separately. If they are approximately the same, then the ratios shown in Figures 9 and 10 will be the same – and the determinant of the matrix $R$ will be approximately equal to zero.

For simplicity let us consider that the gluon-gluon fusion cross section is dominated by the coupling of the Higgs to top quarks, this consideration will allow us to have more insights of the source of discrepancy between the determinants. Eqs.(2) show that $\hat{g}_{tth_{i}}$ has an extra factor $-U_{i2}/\tan\beta$ with respect to the coupling to vector bosons, using the approximation of small $\theta_{23}$ and negligible $\theta_{12}$, the extra factor simplify to $\theta_{23}/\tan\beta$ times $\cos\theta_{13}$($\sin\theta_{13}$) for $h_{1}$($h_{2}$), a factor suppressed by $\tan\beta$. Therefore, unless $\tan\beta$ is close to one, or $\theta_{23}$ is large, we would expect very similar signal strengths for gluon-gluon fusion and vector-boson fusion for each Higgs state, in consequence the total signal strengths for the same final state will be also very similar, and the determinant of $R^{A}$ and $R^{B}$ will be close to zero.

On the contrary, if instead of gluon-gluon fusion production process we consider Higgs production associated to bottom quarks, eqs. (2) show that $\hat{g}_{bbh_{i}}$ has an extra factor $U_{i2}\tan\beta$ with respect to vector boson coupling, the factor is $\tan^{2}\beta$ larger than in the case of $\hat{g}_{tth_{i}}$. For non-negligible values of $\theta_{23}$ there will be a significant departure of signal strength of the Higgs production associated to bottom quarks with respect to the vector-boson fusion for the same final state. When computing the ratio of the total signal strength for different final states we would expect a larger deviation, in consequence the determinant of $R^{C}$ and $R^{D}$ will be different from zero.

These arguments describe very well a set of points with medium to large values of $\tan\beta$. For small values of $\tan\beta$ and large enough values of $\theta_{23}$ the determinant of $R^{A}$ and $R^{B}$ will also show a departure from unity. Figure 11 shows the values of $\theta_{12}$, $\theta_{23}$, $\tan\beta$ and $m_{h_{3}}$ for the pNMSSM posterior sample with $m_{h_{3}}$ larger that 1 TeV and values of $\tan\beta$ larger than 10. As we expected the value of $\theta_{23}/\tan\beta$ is tiny, which explains why the determinant of $R^{A}$ and $R^{B}$ is very close to zero. The large values of $\tan\beta$ also explain the large departure from one for the determinant of $R^{C}$ and $R^{D}$.

Our scan focused on the region of the parameter space with medium to large values of $\tan\beta$, to complete our analysis we analyse a new set of points with smaller values of $\tan\beta$ relative to the first sample set. We perform another small scan giving more preference to the region of small $\tan\beta$ and small $m_{h_{3}}$, covering $\tan\beta$ in the range of $[2.5,21]$ and $m_{h_{3}}$ in the range of [435 GeV to 2 TeV], the results are summarized in figure 12. The top row of the figure shows the values of $\theta_{12}$ and $\theta_{23}$ with respect to $m_{h_{3}}$ and $\tan\beta$. To analyse these two plots in comparison with Figure 11 we have used the same range for the variables plotted in the colour bar to make easier the comparison. First let us focus on the top-left plot of Figure 12. Note that the range of values for $|\theta_{23}\tan\beta|$ is almost the same for both samples suggesting that this parameter is directly constrained by the experimental measurements of the Higgs couplings. Smaller values of $m_{h_{3}}$ are correlated with larger values of $\theta_{12}$, still $|\theta_{23}|$ is one order of magnitude larger than $|\theta_{12}|$, meaning that the approximation of $\theta_{12}\sim 0$ is still valid. The top-right plot of Figures 11 and 12 compare the values of $\theta_{23}\tan\beta$ with $\theta_{23}/\tan\beta$ that illustrate the contribution of $\theta_{23}$ to the Higgs production associated to bottom quarks (x-axis) and gluon-gluon fusion production (y-axis).

The bottom row of Figure 12 show the values of $R^{B}$ and $R^{D}$ for the new set of scanned points. Here, points with $\theta_{23}\tan\beta\sim 0.7$ correspond to $|\theta_{23}/\tan\beta|$ up to 0.030, which is around fifty times larger than our first scan. This increment will be reflected in the value of $R^{B}$, which involves the rate plotted in the left panel of the figure. Previous studies, like [12, 13, 14] pointed out that the determinant of $R^{A}$ and $R^{B}$ will be useful to determine the existence of more than one resonance. Our analyses indicate that this is indeed the case but mostly for pNMSSM regions with relatively smaller $\tan\beta$ values and lighter $h_{3}$. The botton-right plot of Figure 12 shows the relevant ratios to compute the determinant of $R^{D}$. There is a discrepancy in the region with $|\theta_{23}\tan\beta|$ larger than $\sim 0.65$. According to the top-row plots of Figure 12, points with $|\theta_{23}\tan\beta|>0.7$ correspond to $m_{h_{3}}$ smaller that 1 TeV and $\tan\beta$ smaller than 10. Getting relatively larger values for $|\theta_{23}\tan\beta|$ in the new set of points scanned compared to the first pNMSSM posterior sample is in accord with the fact that $|\theta_{23}|$ increases as $m_{h_{3}}$ decreases for a fixed value of $\lambda$ (as discussed in section 3.3). So in the new scan by exploring $m_{h_{3}}<$ 1 TeV, we expand the range of exploration for $|\theta_{23}\tan\beta|$.

6 Conclusions

We studied the phenomenology of the two mass degenerate CP-even Higgs bosons in the NMSSM using a sample set from the parameter scan of the pNMSSM. In this scenario it is possible to reproduce the experimental signal measured by ATLAS and CMS. We parameterised the Higgs boson signal strengths using three angles and found that it is possible to write approximate expressions in terms of two parameters $\theta_{23}\tan\beta$ and $\theta_{13}$, where $\theta_{23}$ is the mixing between the singlet and the heaviest neutral Higgs of the Higgs doublet $H_{0}$ and $\theta_{13}$ the mixing between the lightest neutral scalar of the Higgs doublet and the singlet. We have focused our analysis into observables that could help to determine the existence of more that one Higgs state, leading to the following conclusions.

• •

To obtain two mass degenerate CP-even Higgs bosons there is required tuning associated to large values of $A_{\kappa}$, $\lambda$, $\kappa$ and $\mu$. An approximate relation between those parameters could be obtained from the tree level mass relations, although this relation simplifies the expression for the mass of the lightest pseudoscalar it does not point out to specific mass relations.

• •

An approximate expression for $\theta_{23}$ can be written in terms of $\mu/\lambda$ and $\tan\beta$. The allowed range for $|\theta_{23}\tan\beta|$ is between 0.0 and 0.7. Greater values can be obtained if $m_{h_{s}}\lesssim 1$ TeV and $\tan\beta\lesssim 8$ are imposed. There are no direct constraints on the mass spectra from specific values of $\theta_{23}$ but it is possible to reproduce various values of $m_{h_{3}}$ for a fixed value of $\theta_{23}$ and different values of $\lambda$.

• •

Analysing the Higgs bosons couplings to fermions and vector bosons, and the signal strengths, we found that the signal of the superposition of the Higgs bosons decaying to leptons (and bottom quarks) depart from the SM signal in an opposite direction with respect to vector boson final states. This is proportional to $|\theta_{23}\tan\beta|$.

• •

With respect to expectations due to previous studies, it was surprising to find that for medium to large values of $\tan\beta$, it is rather difficult to distinguish the two degenerate Higgs from the single Higgs scenario when the matrix of signal strengths are for vector-boson and gluon-gluon fusion Higgs productions (with the Higgs decaying to vector boson).

• •

By including Higgs production in association with bottom quarks in the signal strengths square matrix we found that the matrix determinant departs significantly large from the single resonance value.Therefore the process $pp\rightarrow bbh$ can be an important channel in searches for multiple Higgs states degenerate around $125\>\text{GeV}$.

Acknowledgment

Thanks to Alberto Casas for very useful comments and discussions, and to Fernando Quevedo for encouragements towards the NMSSM project. Maria Cabrera thanks ICTP and CERN Theory Division for hosting and supporting her as short-term visitor.