Revealing mass-degenerate states in Higgs boson signals

The observed Higgs boson signals to-date could be due to having two quasi-degenerate 125GeV\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$125 \,\text {GeV}$$\end{document} 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.


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 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. a e-mail: abdussalam@sbu.ac.ir b e-mail: meugenia@ecfm.usac.edu.gt 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 [5,[8][9][10][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 [15,18,19] of the Higgs sector will be used.
In Sect. 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 Sect. 3 we perform a scan of the parameters of the NMSSM while imposing that 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 Sect. 4 the sample is then used together with analytical relations for the couplings and signal strengths to show that the quasi-degenerate Higgs properties can be explained approximately by using just two free parameters. We also show how the superposition of two quasi-degenerate Higgs around 125 GeV could be in agreement with current experimental results. Finally in Sect. 5 we analyse the sample based on signal strength ratios that can discriminate between the single versus double resonance scenarios.

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 γ γ channel motivated a lot of work, some of them within the NMSSM framework [12,[19][20][21][22]. In particular King et. al. [19] pointed out that the signal doublet of the γ γ channels could be enhanced for large singlet-doublet mixing. We will take these as a starting point for analysing two quasidegenerate 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 ψ = (H d , H u , s) and φ = (h 0 , H 0 , s) in such a way that h 0 = v and H 0 = 0: The Higgs states h = (h 1 , h 2 , h 3 ) are related to ψ and φ in the following way, where U i j 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: 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, 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: (I) h 1 or h 2 is the Higgs state detected at the LHC, and (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: (I) Small singlet-doublet mixing, and (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 ∼ 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 → bb and therefore the variation of the width is controlled by the square of the Higgs coupling to bottom quarks, g h 1 bb . Using the couplings described in Eq. (4) the signal strengths of the vector-boson fusion production of h 1 and further decay to W W/Z Z and bb are approximately, whereĝ = g NMSSM /g SM , the couplings g NMSSM are those in Eq. (4), and g SM are the Standard Model (SM) couplings. The enhancement or suppression of the first signal strengths depends on tan β 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 μ VBF→h 1 →WW/ZZ with respect to μ VBF→h 1 →bb . 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 ∼ 1 is not valid any more. The assumption that the width of h 1 is almost totally controlled by h 1 → bb is no longer a good approximation. The size of tan β 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 where h i → SM rest represents the rest of the decay channels.
For large singlet-doublet mixing the widths of h 1 and h 2 could be much smaller than 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. (10) 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 W W/Z Z and bb can be written as, and μ an V BF→h2→bb Note that for a large singlet-doublet mixing the relative size of tan β 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 0 and U 22 0), 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, lim Notice that the last (approximate)equalities require U 31 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 β U 12 /U 11 for h 1 and on tan β 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 .

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].

The phenomenological NMSSM (pNMSSM)
We shall consider an R-parity conserving NMSSM with superpotential, where The chiral superfields have the following Q : 3, 2, The corresponding soft SUSY-breaking terms are with A tilde-sign over the superfield symbol represents the scalar component. However, an asterisk over the superfields as in, for example,ũ * R represents the scalar component ofŪ . The SU (2) L fundamental representation indices are donated by a, b = 1, 2 while the generation indices by i, j = 1, 2, 3. 12 = 12 = 1 is a totally antisymmetric tensor.
In an approach similar to that of the pMSSM [5,[8][9][10], the pNMSSM parameters are defined at the weak scale with the non-Higgs sector set, Here, M 1,2,3 and mf are respectively the gaugino and the sfermion mass parameters. A t,b,τ represent the trilinear scalar couplings. With electroweak symmetry breaking an effective μ-term, μ eff = λ v s is developed. The μ-term, the ratio of the MSSM-like Higgs doublets' vevs tan β = H 2 / H 1 and the Z-boson mass, m Z lead to the tree-level Higgs sector parameters Next, including four SM nuisance parameters, namely, the top and bottom quarks m t,b , m Z and the strong coupling constant, α s , makes the pNMSSM parameters: 3.2 The scanning procedure . tan β is allowed between 2 and 60. For minimising fine-tuning, we subjectively let μ eff = λ 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.

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 ± 3 GeV, where the ±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. 1 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. 1 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 h2 − m h1 < 3 GeV can be considered as a mass degeneracy condition for which the two Higgs cannot be resolved by CMS run-2.
It is useful to have an explicit form for the Higgs mixing matrix U . We parameterise this using three angles θ 13 , θ 12 , and θ 23 such that Here c i j = cos θ i j and s i j = sin θ i j . Given the mixing matrix, obtained numerically by the SUSY spectra calculator NMSSMtools, then the mixing angles can be extracted as: 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 μ V B F→h 1,2 →W W/Z Z and μ V B F→h 1,2 →bb close to one, either the value of μ V B F→h i →W W/Z Z /μ V B F→h i →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 case. 2 From Eqs. (12)- (15) 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. (29). On an other hand, Eq. (16) implies that the superposition of h 1 and h 2 can reproduce the standard Higgs signal for U 31 ∼ 0 (i.e. large values of m H 0 ). For this to happen either θ 12 has to be very small or θ 23 has to be close to ±π/2. In summary, θ 12 ∼ 0 and θ 23 ∼ 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 θ 12 and θ 23 , where we have neglected O(θ 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 standardlike 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 θ 13 , θ 23 and θ 12 to the fundamental parameters of the Higgs sector. Using Eq. (29) we relate the terms of the mass matrix with the physical masses by introducing two new parameters: m 2 h , the central value of the two lightest CP-even Higgs states, and δm 2 h , half of the squared mass difference, To simplify the expressions obtained from Eq. (33) we factorise c 12 and c 23 to write U in terms of t kl ≡ tan θ kl and use the approximations: where kl = 12, 23 and mn = 12, 23. Finally, we will focus on the relations in terms of the mass matrix elements M 22 and M 23 since M tree 22 and M tree 23 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.
We can further simplify Eq. (35) where in the last equation we have further considered that δm 2 h cos(2θ) m 2 h . Using the approximation of large tan β and large M A from reference [115]: We have checked numerically that m h 3 ≈ M A is a good approximation for the pNMSSM points considered. Now, let us take M 23 from reference [115] 3 where v s = √ 2 μ/λ and ξ = μ/M A . Left panel of Fig. 1 shows in the x-axes the value of θ 23 computed by NMSSMtools and in the y-axes the analytical approximation described in Eq. (41), 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 θ 23 comes from the first term of Eq. (41) (blue points). Right panel of Fig. 1 shows the relation between θ 23 and m h 3 for constant values of λ. There is a trend: larger values of |θ 23 | correspond to smaller values of m h 3 , except for very small values of |θ 23 | where the two parameters seem to be uncorrelated. Still, Eq. (41) shows that the value of tan θ 23 is not directly related to the scale of the heaviest Higgs, but instead it is related to the value of λ, μ and tan β. 4 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 β and M 2 A , where v S = √ 2μ/λ (see Eq. (32) of [115]). In order to get a constrain for the initial parameters from the condition of 4 Let us remember that in the decoupling limit of H 0 , small mass difference between the two lightest Higgs states, we require a small mass difference between the tree level masses showed in Eq. (42). 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 2 Z , meaning that both terms inside the square root should be smaller than M 4 Z . Let us focus on the first term, for A κ M Z there should be a correlation between A κ and κv s such that there is a cancellation that leads to an order M 2 Z 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 κ | M Z we expect, Figure 2 shows the relation between A κ and κv s , as manifested in the figure for |A κ | 600 GeV the approximation of Eq. (43) works within an error smaller than 5%.
Furthermore, using Eq. (43) 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, Putting Eq. (43) into Eq. (44) we write the mass of the lightest pseudoscalar in terms of κ and v s , Figure 3 shows the comparison between Eq. (45) and the value computed by NMSSMtools. It can be seen that for m A 1 > 500 GeV Eq. (45) 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. (42) is suppressed by a factor v −2 s , as such we do not expect to get any good correlation of parameters from there.  Comparison between the mass of the lightest pseudoscalar computed by NMSSMtools and the approximate analytical value described in Eq. (45). The colour code shows the value of A κ , which as described in Eq. (43) it is related with the value of κv s 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 massdegenerate Higgs region(s).

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. 5 Figure 4 shows the comparison between the signal strengths 5 To 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 computed by NMSSMtools, μ num , and the analytic approximations showed in Eqs. (12)- (15), μ an , for VBF→ h i → W W/Z Z (left panel) and VBF→ h i → 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 θ 12 and θ 23 , see Eqs. (4) and (32), Using Eqs. (10) and (46) we get, Finally, Eqs. (12)-(13) can be written in terms of the mixing angles as   (54) From Eqs. (51)-(54) we see that the signal strengths depend on four parameters: θ 13 , θ 23 , θ 12 and tan β. However, in the limit where θ 12 tan β θ 13 , which is the case for the pNMSSM posterior sample analysed, the number of parameters reduces to two: From Eqs. (51)-(54), one can see that the dependence on θ 12 always appears as a factor in the expression cos θ 13 (1 + θ 12 tan β) or sin θ 13 (1+θ 12 tan β). Therefore for θ 12 tan β 1 the contribution of θ 12 is negligible.
To understand the dependence of the signal strengths with respect to θ 13 and θ 23 tan β let us start analysing the relation between the signal strengths for a given Higgs state. The top row of Fig. 5 shows the correlations between μ V B F→h i →W W/Z Z and μ V B F→h i →bb for h 1 (top left) and h 2 (top right); for h 1 we can see that the difference between the bb and W W/Z Z 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 θ 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. (49) and (50), show that the term proportional to θ 23 tan β 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 → bb while increasing(decreasing) the decay rate of h 2 → bb as |θ 23 | increases its value.
The bottom row of Fig. 5 shows the width of h 1 and h 2 as function of θ 13  Let us analyse the global signal strengths. Figure 6 shows the sum of the signal strengths of vector-boson fusion production and decay to W W/Z Z (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. 6 There are several points we would like to comment from Fig. 6, the departure of the signal strength increases with the size of θ 23 tan β 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 There are two regions that seem to be in full agreement with the SM (the signal strength is 1): the region where θ 23 0 and the region where θ 13 0, as we expected. There is a third region where θ 13 is between 0.2 and 0.4, where for a very precise value of θ 23 the signal strength is very close to one. On the other hand, for small values of θ 23 , let's say θ 23 tan β −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 Figs. 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: θ 13 and θ 23 tan β, instead of three, indicating that θ 12 tan β 1 for the set of successful scanned points.
So far we have focused our study to two channels: VBF→ h i → W W/Z Z and VBF→ h i → bb, but the current measurements of the Higgs couplings constrain several more channels. Let us comment about the most relevant production and decays: (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 Eq. (4), which describe the couplings of the Higgs states to top quarks, Comparingĝ g i tt withĝ g i bb we see that the contribution from θ 23 is cot 2 β times smaller forĝ h i tt than forĝ h i bb , therefore we expect the contribution of θ 23 to be very tiny and the production processes of GGF and ttH to behave as vector-boson fusion for given values of θ 13 and θ 23 tan β.
(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. (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 Fig. 7 we show the signal strengths for GGF→ h 1,2 → W W/Z Z (left panel) and GGF→ h 1,2 → γ γ (right panel). As we expected, the gluongluon 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 θ 13 and θ 23 tan β. In the limit where θ 12 0, we could write a biunivocal function to determine one (of these parameters) in terms of the other. An approximate relation between θ 13 and θ 23 tan β might be useful to study the region around 0.2 θ 13 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: By taking μ V B F→h 1 →W W/Z Z and μ V B F→h 2 →W W/Z Z from Eqs. (51) and (52), neglecting the terms proportional θ 12 , and rewriting the sin θ 13 and cos θ 13 in terms of sin(2 θ 13 ) and cos(2 θ 13 ) we can simplify Eq. (55) to get a quadratic equation in cot(2 θ 13 ). So, there are two solutions for θ 13 : where BR bb = 1 − BR bb . For δ = 0 the solution simplifies to With Eq. (56) we are able to determine θ 13 in terms of θ 23 tan β and δ. Figure 8 shows the comparison between the semi-analytical relation in Eq. (56) and the numerical results from our scans. Although it is not a precise relation, Eq. (56) gives a very good approximation to the correlation between θ 13 and θ 23 for a fixed value of δ.

Searching for mass-degenerate Higgses
As commented in references [12,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: μ i j = μ i→ j , where i represents the production mode and j the decay channel. Considering two square matrices, the condition for the determinant to be non-zero can be written in terms of the ratios 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 Sect. 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. (59). The upper (lower) panel shows all the points that are within three (one) sigma of the val-   Table 3), it doesn't include the SM value at one sigma. The left panel of Fig. 9 shows that the ratios between W W and γ γ 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 τ τ and γ γ 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 μ V B F,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 gluongluon fusion.
Let us consider the matrices, 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: 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 Figs. 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. Equations (4) show thatĝ tth i has an extra factor −U i2 / tan β with respect to the coupling to vector bosons, using the approximation of small θ 23 and negligible θ 12 , the extra factor simplify to θ 23 / tan β times cos θ 13 (sin θ 13 ) for h 1 (h 2 ), a factor suppressed by tan β. Therefore, unless tan β is close to one, or θ 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, Eq. (4) show thatĝ bbh i has an extra factor U i2 tan β with respect to vector boson coupling, the factor is tan 2 β larger than in the case ofĝ tth i . For non-negligible values of θ 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 β. For small values of tan β and large enough values of θ 23 the determinant of R A and R B will also show a departure from unity. Figure 11 shows the values of θ 12 , θ 23 , tan β and m h 3 for the pNMSSM posterior sample with m h 3 larger that 1 TeV and values of tan β larger than 10. As we expected the value of θ 23 / tan β is tiny, which explains why the determinant of R A and R B is very close to zero. The large values of tan β 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 β, to complete our analysis we analyse a new set of points with smaller values of tan β relative to the first sample set. We perform another small scan giving more preference to the region of small tan β and small m h 3 , covering tan β 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 Fig. 12. The top row of the figure shows the values of θ 12 and θ 23 with respect to m h 3 and tan β. To analyse these two plots in comparison with Fig. 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 Fig. 12. Note that the range of values for |θ 23 tan β| 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 θ 12 , still |θ 23 | is one order of magnitude larger than |θ 12 |, meaning that the approximation of θ 12 ∼ 0 is still valid. The top-right plot of Figs. 11 and 12 compare the values of θ 23 tan β with θ 23 / tan β that illustrate the contribution of θ 23 to the Higgs production associated to bottom quarks (x-axis) and gluon-gluon fusion production (y-axis).
The bottom row of Fig. 12 show the values of R B and R D for the new set of scanned points. Here, points with θ 23 tan β ∼ 0.7 correspond to |θ 23 / tan β| 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 β values and lighter h 3 . The bottom-right plot of Fig. 12 shows the relevant ratios to compute the determinant of R D . There is a discrepancy in the region with |θ 23 tan β| larger than ∼ 0.65. According to the top-row plots of Fig. 12, points with |θ 23 tan β| > 0.7 correspond to m h 3 smaller that 1 TeV and tan β smaller than 10. Getting relatively larger values for |θ 23 tan β| in the new set of points scanned compared to the first pNMSSM posterior sample is in accord with the fact that |θ 23 | increases as m h 3 decreases for a fixed value of λ (as discussed in Sect. 3.3). So in the new scan by exploring m h 3 < 1 TeV, we expand the range of exploration for |θ 23 tan β|.

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 θ 23 tan β and θ 13 , where θ 23 is the mixing between the singlet and the heaviest neutral Higgs of the Higgs doublet H 0 and θ 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 κ , λ, κ and μ. 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 θ 23 can be written in terms of μ/λ and tan β. The allowed range for |θ 23 tan β| is between 0.0 and 0.7. Greater values can be obtained if m h s 1 TeV and tan β 8 are imposed. There are no direct constraints on the mass spectra from specific values of θ 23 but it is possible to reproduce various values of m h 3 for a fixed value of θ 23 and different values of λ.
• 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 |θ 23 tan β|. • With respect to expectations due to previous studies, it was surprising to find that for medium to large values of tan β, it is rather difficult to distinguish the two degener-ate 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 → bbh can be an important channel in searches for multiple Higgs states degenerate around 125 GeV.
Acknowledgements 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.

Data Availability Statement
This manuscript has no associated data or the data will not be deposited. [Authors' comment: The scanned points used in this analysis can be obtained using the free source codes as described in Sect.

3.2.]
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecomm ons.org/licenses/by/4.