NMSSM interpretations of the observed Higgs signal

While the properties of the signal that was discovered in the Higgs searches at the LHC are consistent so far with the Higgs boson of the Standard Model (SM), it is crucial to investigate to what extent other interpretations that may correspond to very different underlying physics are compatible with the current results. We use the Next-to-Minimal Supersymmetric Standard Model (NMSSM) as a well-motivated theoretical framework with a sufficiently rich Higgs phenomenology to address this question, making use of the public tools HiggsBounds and HiggsSignals in order to take into account comprehensive experimental information on both the observed signal and on the existing limits from Higgs searches at LEP, the TeVatron and the LHC. We find that besides the decoupling limit resulting in a single light state with SM-like properties, several other configurations involving states lighter or quasi-degenerate with the one at about 125 GeV turn out to give a competitive fit to the Higgs data and other existing constraints. We discuss the phenomenology and possible future experimental tests of those scenarios, and compare the features of specific scenarios chosen as examples with those arising from a more global fit.


Introduction
After the discovery of a signal with a mass of about 125 GeV in the Higgs searches at the LHC [1,2], the prime goal is now to identify the underlying nature of the new state and to determine the mechanism of electroweak symmetry breaking. While the properties of the observed state are compatible with the ones predicted for the Higgs boson of the Standard Model (SM) at the current level of precision, also a wide range of alternative interpretations could be possible, corresponding to very different underlying physics. In particular, in models with an extended Higgs sector the observed state would be accompanied by several other Higgs bosons, in contrast to the minimal formulation of the SM where a single SU(2) L -doublet is responsible for electroweak-symmetry breaking.
Supersymmetry (SUSY) [3] is commonly regarded as the most appealing extension of the SM, since it provides a solution for stabilising the huge hierarchy between the Planck scale and the weak scale [4][5][6][7][8] and offers further attractive features such as unification of the gauge couplings and a natural candidate for cold dark matter in the Universe. A crucial JHEP04(2016)095 prediction of supersymmetric extensions of the SM are their extended Higgs sectors: the holomorphicity of the superpotential (as well as the cancellation of gauge anomalies) implies that at least two SU(2) L doublets with opposite hypercharge have to be present, so as to generate mass terms for both up-and down-quarks (in a Type II 2-Higgs-Doublet-Model fashion). The minimal supersymmetric extension of the SM (MSSM) [9,10] is based on the minimal Higgs sector of this kind comprising two Higgs doublets, whereas the Higgs sector of the next-to-minimal supersymmetric extension of the SM (NMSSM) [11,12] contains an additional (complex) gauge-singlet. It has long been recognised that the NMSSM provides an elegant solution [13][14][15] to the "µ-problem" [16] of the MSSM. In the context of the relatively high mass value of about 125 GeV of the observed state, this model has received particular attention lately since the mass of the light doublet-like state receives an additional contribution at lowest order as compared to the MSSM, which dominates at low values of tan β (the ratio of the vacuum expectation values, v.e.v.'s, of the two Higgs doublets). In this case significantly smaller higher-order corrections are required to obtain a Higgs-boson mass in the appropriate range [17] as compared to the MSSM, where the lowest-order prediction for the mass of the light CP-even Higgs boson is bounded from above by the mass of the Z boson, M Z . Furthermore, also the singlet-doublet mixing can give rise to an uplift of the mass of the doublet state, provided that the CP-even singlet state is lighter than the doublet state. It has been argued in this context that the relaxed requirement on the size of the higher-order corrections as compared to the MSSM makes it possible to obtain a Higgs-mass prediction of about 125 GeV in a "more natural" way [18][19][20][21][22][23][24][25].
In the following we will focus on the NMSSM as a theoretically well-motivated alternative to the SM with a potentially rich phenomenology in the Higgs sector. For simplicity, we will restrict to the CP-conserving case, for which the spectrum of physical Higgs states of the NMSSM consists of three CP-even, two CP-odd and a pair of charged Higgs states (while we do not explicitly consider CP-violating effects giving rise to a mixture between the five neutral states, it should be noted that cases where a CP-even and a CP-odd state are nearly mass-degenerate essentially mimick a scenario where a single state is an admixture of CP-even and CP-odd components). Furthermore, while several versions of the NMSSM can be formulated, depending on the form of the singlet and singlet-doublet interaction terms in the superpotential, we will focus on the Z 3 -conserving version only, where the solution to the "µ-problem" is more immediate (on the other hand, this simple model could lead to a domain wall problem [26] but we will not address this question here). The Higgs sector of this version of the NMSSM is characterised by six parameters (at tree-level), in contrast to the two parameters of the MSSM. While we shall borrow most of our notations from [11], we recall the Higgs terms entering the superpotential of the model: where S denotes the singlet (super)field, H u and H d the doublets, while · stands for the SU(2) product. When confronting the predictions of an extended Higgs sector with the observed signal and the limits from the Higgs searches at LEP, the TeVatron and the LHC, the most obvious interpretation of the signal at about 125 GeV is to associate it with the lightest JHEP04(2016)095 CP-even Higgs boson of the considered model. The case where all other Higgs (as well as all new physics) states of a supersymmetric extension of the SM (and the same is true for various other extended Higgs sectors) are significantly heavier corresponds to the "decoupling region" of the model under consideration, where the couplings of the light Higgs boson to gauge bosons and SM fermions are very close to the ones of the SM. Revealing deviations of those couplings from their SM counterparts in such a case will require high-precision measurements, where in many cases the expected deviations do not exceed the level of a few per cent. An additional source of possible deviations from the SM could be decays of the SM-like state into new-physics particles. Such a decay could in particular occur into a pair of dark matter particles, if the mass of the latter is less than half of the mass of the Higgs state, i.e. below about 60 GeV. This would give rise to an invisible decay mode of the observed state, providing a strong motivation for searches of decays of the observed signal into invisible final states.
Besides the interpretation of the observed state as the lightest CP-even Higgs boson of an extended Higgs sector, it is also possible, at least in principle, to identify the observed signal with the second-lightest state of an extended Higgs sector. This interpretation would have the immediate consequence that there should be an (or more generally at least one) additional Higgs state below the one observed at about 125 GeV. The phenomenology of such a scenario is very different from the case of the decoupling limit discussed above, because of the presence of at least one more light state in the spectrum. Within the MSSM this interpretation is in principle possible [27][28][29][30][31][32], but gives rise to a rather exotic Higgs sector where in fact all additional Higgs bosons are light, i.e. in the vicinity of the state at about 125 GeV or below. It is remarkable that a global fit within the MSSM within this interpretation has resulted in an acceptable fit probability [32], but lately this interpretation, which implies in particular a light charged Higgs boson below the mass of the top quark (see in particular ref. [33]), has come under increased pressure from the limits obtained in the charged Higgs searches by ATLAS [34] and, more recently, CMS [35].
The NMSSM provides a well-suited and theoretically well motivated framework for investigating to what extent interpretations that go beyond the obvious case of a single light state in the decoupling limit are compatible with the latest experimental results both with respect to the properties of the observed state and to the limits obtained from Higgs searches (as well as other existing constraints). It is the purpose of the present paper to perform such an analysis. It is obvious that compatibility with the observed signal requires much more than just a Higgs state (or possibly more than one) in the spectrum with a mass of about 125 GeV. In order to properly take into account the latest experimental results from Higgs search limits and from measurements of the properties of the observed state, we make use of the public tools HiggsBounds [36][37][38][39][40] and HiggsSignals [41], which incorporate a comprehensive set of results from ATLAS [42], CMS [43] and the TeVatron [44,45]. We do not explicitly impose limits from the direct searches for SUSY particles at the LHC. While some of the scenarios discussed in this paper could be affected by constraints from SUSY particle searches, we have checked that the qualitative features of the Higgs phenomenology of those scenarios are maintained also for somewhat heavier SUSY particle spectra. Note however that most of the SUSY spectra that we employ (especially for coloured particles) are beyond the mass-range tested in the Run-I of the LHC.

JHEP04(2016)095
As mentioned above, already since the very early hints of a signal at about 125 GeV the Higgs sector of the NMSSM has found a lot of attention in this context. Besides the mass prediction in comparison with the MSSM case [18][19][20][21][22][23][24][25], the possibility of modified rates has been discussed, particularly in the diphoton channel [46][47][48][49][50][51][52]. Furthermore, the case of universal (or semi-universal) SUSY-breaking conditions at the GUT scale [53][54][55][56][57][58], gaugemediation [59,60] and other related scenarios [61][62][63][64] have been considered in this context. The CP-violating version of the NMSSM also received some attention [65,66]. [67][68][69] analysed the fine-tuning in a Z 3 -violating version of the NMSSM and variants. Other groups confronted the presence of a Higgs state at this mass with direct searches for SUSY particles at the LHC or Dark-Matter constraints [70][71][72]. Scenarios with a light singletlike state around ∼ 100 GeV have found considerable interest [73][74][75]. Another possibility involves a singlet and a doublet that are almost mass-degenerate at about 125 GeV and may mix with each other, see ref. [76,77] (and ref. [52]). Several studies also suggested to exploit pair production processes at the LHC in order to distinguish the SM from the NMSSM and/or to look for a light singlet in this fashion [78][79][80][81][82][83]. Scenarios with a very light CP-odd (or CP-even) Higgs boson were addressed with several search proposals in direct production, unconventional light charged-Higgs decays, or cascade decays from SM-like / light singlet states; large Higgs-to-Higgs decays were also considered from the point of view of the SM-compatible nature of the observed state [84][85][86][87]. Recent studies of the properties of a light pseudoscalar in the NMSSM [88,89] have emphasized the relevance of indirect production modes for the investigation of this scenario at the LHC. In a different direction, the authors of [90] focussed on NMSSM Higgs scenarios with a lowscale doublet sector. Furthermore, [91], and more recently [92,93], studied the discovery prospects of NMSSM Higgs states in the LHC run at 13 TeV. In our analysis we go beyond the previous work in several respects: while many of the afore-mentioned analyses discussed scenarios which are compatible with existing limits, our inclusion of a fitting tool allows us to highlight the quality of the various scenarios in view of the available data. Furthermore, we aim at a comprehensive discussion from the point of view of the NMSSM Higgs phenomenology, hence do not confine to a specific scenario (within our assumptions on the model, perturbativity constraints and choices of simplicity with regards to the SUSY spectrum). We also focus on Higgs physics and thus try, without spoiling the physical content, to avoid emphasis on questions of secondary importance with respect to this topic (e.g. the details of the supersymmetric spectrum). Finally, much experimental data has become available in the last few years, narrowing the possibilities in the Higgs sector, and most of the recent developments are included within the tools on which our discussion is based.
The paper is organised as follows: in section 2 we describe the framework that we use for the analyses in this paper, in particular the statistical approach used in our fits, the treatment of external constraints and the tools that we apply. As a first step of our analysis, in section 3 we briefly consider the SM case and the corresponding decoupling limits of the MSSM and the NMSSM. The SM result is used for comparison with χ 2 analyses in different NMSSM scenarios, which we perform in sections 4-8. In section 9, we focus on specific points of the NMSSM parameter space and discuss in more details the

JHEP04(2016)095
Higgs phenomenology and the consequences for future searches, should the corresponding spectrum be realised in nature. In section 10 a more global scan is carried out, and the features observed in the global scan are discussed in view of the results obtained for the specific NMSSM scenarios that we have considered before. Section 11 contains our conclusions.
2 Framework of the analysis: treatment of external constraints and applied tools The NMSSM parameter space is explored with the help of the spectrum generator NMSSMTools 4.4.0 [94][95][96][97][98], computing the Higgs masses up to leading two-loop doublelog order (we will be using the default mode only), in an effective potential approach. This code considers a certain number of phenomenological limits, several among which are kept within our analysis. The first class of such tests are consistency requirements and are (necessarily) included as hard cuts: • stability of the EWSB-vacuum: positivity of the scalar squared-masses, absence of deeper minimum; • absence of Landau poles below the GUT scale: while this requirement is sometimes omitted in order to probe effects associated with large values of the parameter λ and under the assumption that new-physics or specific properties of the non-perturbative regime would smoothen the theoretical difficulty of the Landau poles, we choose to keep this theoretical limit; • requirement for Higgs soft squared-masses at the TeV scale: the potential-minimization procedure in NMSSMTools trades these masses for the Higgs v.e.v.'s, so that the naturalness requirement that soft masses intervene at the TeV scale must be checked explicitly; • requirement for a neutralino LSP (the impact of which, however, is of secondary importance in our discussion).
Another type of constraints are supersymmetric searches at LEP. Given that we are chiefly interested in the Higgs sector, we also keep these limits under the form of a hard cut: • mass lower limits on squarks (mt > 93.2 GeV, mq 1,2 > 100 GeV), gluino (mg > 180 GeV), sleptons (ml > 99.9 GeV), charginos (mg > 103.5 GeV); We remind the reader that LHC limits on SUSY searches are not considered in our analysis. However, the SUSY spectra that we employ are typically beyond the mass-range of the searches in the Run-I. In this context, the inclusion of LEP limits as mentioned above has only a minor impact.

JHEP04(2016)095
NMSSMTools also computes several low-energy observables: • limits from the bottomonium sector: non-observation of a signal in BR(Υ → γA → γl + l − ), excessive contribution to the η b (1S)-A mixing [99][100][101]. Only light CP-odd Higgs below ∼ 10 GeV are concerned by these constraints and the limits are kept as a 95% C.L. cut. 1 • limits from B-factories (under a strong Minimal Flavour Violation hypothesis, i.e. neglecting all possible tree-level flavour-changing neutral currents): [102]. Instead of treating the limits as a hard cut, we combine them in a χ 2 function relying on the central value and error bars computed in NMSSMTools: i are summarized in table 1. The theoretical error estimate is the result of an involved calculation: errors relative to SM-like contributions are taken from the corresponding SM estimate; the uncertainty on new-physics contributions is estimated to 30% (if only leading-order effects are included) / 10% (if next-to-leading α S corrections are present) of the total corresponding contributions and are added linearly to the SM error; additional error sources are mostly CKM matrix elements (taken from tree-level measurements exclusively) and hadronic parameters (decay constants, taken from lattice calculations); to obtain the final theoretical uncertainty range, both higherorder and parametric uncertainties are varied within these previously-discussed limits.
• (g − 2) µ [103]: similarly to B-observables, we add a contribution to the χ 2 with experimental-SM input shown in table 1, where the errors have been added in quadrature. The theoretical uncertainty associated to new-physics contributions and higher orders is calculated as the sum of a fixed error 2.8 · 10 −10 , a 2% error estimate on 1loop contributions (which do not involve coloured particles) and a 30% error estimate on 2-loop effects (involving coloured particles).
Additionally, given that a candidate for the interpretation of the signal observed at the LHC seems necessary, we require that the spectrum produces one CP-even Higgs state in the mass-range [120, 131] GeV. Other limits are deliberately ignored, at least as implemented within NMSSMTools: • dark matter searches: relic-density (via MicrOMEGAs), XENON 100 [104,105]. The reasons for not taking such limits into account come from the observation that they are strongly dependent on the SUSY spectrum, while we want to focus on the Higgs sector: confining to collider constraints allows us to handle simple supersymmetric spectra, which play a secondary part in our analysis, while these would likely have to JHEP04(2016)095 i 3.43 · 10 −4 1.14 · 10 −4 3.2 · 10 −9 1.6 · 10 −6 . . . be finely adjusted if one were to include, e.g., the relic-density bounds. We note also that the dark matter phenomenology may involve mechanisms (e.g. light gravitino LSP) which may alter the conclusions in the dark sector, while all such considerations are not the focus of our discussion.
LEP, TeVatron and LHC limits on Higgs searches are checked through the code HiggsBounds 4.2.0 [36][37][38][39][40], which we interfaced to NMSSMTools via a subroutine. Hig-gsBounds is used with the default settings, using hard cuts on the allowed regions and an individual test for all Higgs states ('full' method). Note that the more sophisticated implementation of limits on the Higgs sector within HiggsBounds leads to divergences with the channel-after-channel checks implemented within NMSSMTools, in particular in the LEP e + e − → hA channels: as a consequence, lighter Higgs (doublet-like) states are accessible, while the corresponding points would be rejected by NMSSMTools. The version of HiggsBounds which we use includes all the released LHC limits till the end of 2014.
Another remark concerns the distinction between 'hard-cut constraints' and 'observables to include in our fit': it is bound to carry some arbitrariness. One could object, for instance, that implementing LEP limits as 'hard cuts', one loses the sensitivity to small deviations in the LEP data so that no benefit is associated in the fit. On the other hand, those may prove to be only statistical fluctuations. In practice, we treat all search limits as hard cuts while observables kept in the fit have been actually measured. HiggsSignals 1.3.1 [41] performs the comparison to the TeVatron+LHC-observed Higgs data, delivering a χ 2 fit to the Higgs-measurement observables. The version we employ collects all released experimental material till the end of 2014. Here follow a few remarks concerning the setting of the options that are offered.

JHEP04(2016)095
• Statistical test method -Two statistical methods, 'peak-centered' (comparing theoretical and experimental signals at masses determined by the experimental signals) and 'mass-centered' (comparing theoretical and experimental signals at the masses that are defined by the theoretical input) are available in HiggsSignals. Only the 'peak-centered' test shall be performed. 81 channels are tested within the analysis.
• Higgs pdf -The probability density function assumed for the Higgs boson mass is modelled as a gaussian.
• Theoretical mass uncertainties -We allow for a ±3 GeV uncertainty on the theoretical Higgs mass predictions delivered by NMSSMTools. This ensures that the phenomenology that we are discussing applies to the NMSSM and not to an arbitrary spectrum. Given also that our study is limited by the density of the scans that we perform, the mass uncertainty smoothens the impact of phenomenological limits. On the other hand, one could argue that, in view of the Higgs-like signals at LHC, the interesting part of the phenomenological study concerns the workings of a given spectrum, irrespectively to the model hiding behind. Moreover, considering that the NMSSM has a large number of parameters in the Higgs sector, it may be possible to absorb a substantial part of the higher-order corrections within a shift in parameter space. In such a context, one could be concerned that large mass uncertainties might blur the phenomenology and provide unlikely spectra with undue attention. Yet this feature does not happen and the scenarios that we propose would appear with qualitatively comparable fit values even though the mass uncertainty was set to be smaller.
This χ 2 test to the Higgs-measurement data, delivered by HiggsSignals, is added to the corresponding tests from B-physics and (g − 2) µ . The resulting quantity, χ 2 (total) , will be at the center of our discussion in the following sections. Let us comment briefly on its interpretation. The χ 2 is the sum of squared deviations between experimental measurements and theoretical inputs, weighed by the corresponding uncertainties. Assuming the corresponding quantities are random and independent (gaussian) variables, the χ 2 would follow a probability law given by the χ 2 -distribution of N th degree P N , where N is the number of variables in the sum: in our case, N = 81 + 7 + 1 = 89, since there are 81 channels involved in the HiggsSignals test, 7 B-observables and (g − 2) µ finally. Correspondingly, one may define the compatibility of the χ 2 test ('P-value') as the probability to fall farther than the obtained χ 2 value (i.e. the probability that the generated deviations be larger). With this definition, the compatibility of a χ 2 value with the test involving N degrees of freedom is obtained as: C(χ 2 ) ≡ ∞ χ 2 P N (x)dx. In this context, for the χ 2 distribution of 89 th degree, the compatibility of the spectrum with the data reaches 90% for χ 2 72, 50% for χ 2 88, 30% for χ 2 95, 5% for χ 2 112 and 1% for χ 2 123 (see figure 1). (Note that HiggsSignals directly provides an approximate P-value for the Higgs-fit.) This test would be the statistically relevant confrontation of one isolated point to the experimental data. This approach raises a few issues, however, as the corresponding statistics is then critically dependent on the list of tested channels and the precise definition of these. Moreover, JHEP04(2016)095 we will not be considering isolated points but scan over various portions of the NMSSM parameter space, thus introducing degrees of freedom which would have to be substracted from the statistical test. Absolute χ 2 values are thus difficult to interpret. Consequently, we will base our phenomenological discussion on relative χ 2 -differences with respect to a best-fit point, which prove to be a more robust interpretation in the given context.

SM and decoupling limits
The Standard-Model limit of the NMSSM is obtained when decoupling the singlet sector -MSSM limit: λ ∼ κ → 0; the singlet states then have vanishing couplings to their doublet counterparts, while an effective µ-term, µ ≡ λ S is generated from the singlet v.e.v. -and pushing the masses of the MSSM non-standard states to very large values. The scale of the heavy doublet sector is controlled at tree-level by the parameter M Athe doublet diagonal entry in the CP-odd Higgs mass matrix at tree-level -which can be used directly as an input (instead of A λ ) within NMSSMTools: the decoupling condition reads M A M Z . Similarly, the following scales enter the supersymmetric spectrum: the sfermion, mf , gaugino, M 1,2,3 , and higgsino, µ, masses, which can be chosen far from M Z . At low energy, one is then left with an effective SM, whose Higgs boson, the remaining light doublet Higgs state, has indeed SM-like couplings and a mass falling within the appropriate range (∼ 125 GeV), provided soft stop terms A t and moderate-to-large tan β > ∼ 10 are chosen accordingly. New physics effects are then suppressed in accordance with the high scales that they involve or to the vanishing couplings of new states to the SM ones, so that this limit is virtually undiscernible (at low energy) from a genuine SM. 2 We consider this trivial scenario in order to 'calibrate' our fit -i.e. set a point of comparison with other NMSSM scenarii -and display our results in figure 2: the χ 2 -value is plotted as a function of the mass of the SM-like Higgs for a multi-TeV heavy supersymmetric and second-Higgs doublet spectrum. The best-fit points receive a χ 2 of about ∼ 87, which statistically places the SM-limit within 1σ compatibility with the considered observables. This fact is not JHEP04(2016)095 While the uncertainty on the Higgs masses, ∆m th , is set to 3 GeV in the rest of the paper, we also include a plot for ∆m th = 1 GeV here, in order to illustrate the impact of this quantity on the fit.
surprising, since the measurements of the Higgs signal at the LHC are grossly consistent with a SM interpretation (within 1σ). Similarly, no tension develops in the B-sector, where the considered observables are also compatible with the SM. On the other hand, the SM-limit is difficult to reconcile with the anomalous magnetic moment of the muon, which generates a χ 2 -pull of ∼ 8-9 (depending on the chosen scale for the sleptons, ml): a limited departure from the strict SM limit -e.g. lowering sleptons and gaugino masseswould remedy this discrepancy. Finally, we wish to comment on the general aspect of the curves in figure 2: the minimum appears there as a broad step, extending from ∼ 123 GeV to ∼ 128 GeV. This appearance is driven by the treatment of the theoretical uncertainty on the Higgs masses, ∆m th , within HiggsSignals, as appears clearly when comparing both cases ∆m th = 1 GeV and ∆m th = 3 GeV.
Without turning to the full SM limit, the most frequent interpretation of the Higgs signal at ∼ 125 GeV within supersymmetric models consists in identifying it with the lightest Higgs state, this prejudice being motivated by the current absence of conclusive experimental signals for a Higgs boson at lower mass values. The corresponding configuration is most naturally achieved in the decoupling limit, that we define as M A M Z , without necessarily requiring that the supersymmetric or the singlet spectra are much heavier: the light doublet state is then largely SM-like, at least at tree-level. The possible presence of new-physics particles at neighbouring scales might then affect the couplings of this light Higgs state at the radiative level, which would allow for tests in precision physics -unless the induced effect is negligible; note that the current LHC results allow for a relatively broad range of coupling strength in the vicinity of the SM values. The major concern in this configuration of the Higgs spectrum actually lies in generating a mass for the light doublet state in the appropriate range in order to identify it with the measured signals. In  the MSSM limit (λ ∼ κ → 0), this can be achieved by saturating the tree-level contribution to the mass ∝ M Z cos 2β -therefore turning to tan β > ∼ 10 -and relying on substantial loop corrections (heavy stops / large trilinear couplings). This solution can be extended to the NMSSM (i.e. departing from λ, κ 1), although specific NMSSM effects can also be employed, as we will show in the following sections. Here, we just illustrate the decoupling limit in figure 3, where we display the results of a scan in the plane {M A , tan β} both for the MSSM limit and a case with non-vanishing λ and κ. Both configurations lead to a fit result where the lowest χ 2 values are obtained in the range of large M A (O(TeV)) and significant tan β, i.e. in the decoupling limit. The best-fit point is indicated in the plots by a flag. Note that the preference for tan β > ∼ 15 is driven only partially by the requirement of a Higgs mass close to ∼ 125 GeV. The discrepancy of the SM with the anomalous magnetic moment of the muon can indeed be cured by supersymmetric contributions, which then favour sizable tan β: this is the main pull in the tan β direction, otherwise the χ 2 distribution is mostly flat as soon as tan β > ∼ 10. It should be noted that, while (g − 2) µ in general favours large values of tan β for a relatively heavy SUSY spectrum, the preferred range does of course depend on the details of the SUSY spectrum, and especially the masses of the sleptons and charginos / neutralinos. We also wish to comment that the main limiting factor at low tan β, in figure 3, rests with the requirement of a Higgs state close to 125 GeV. As such this (LHC) constraint supersedes the LEP lower bound of ∼ 114 GeV for the mass of a SM-like state. Note also that the SUSY spectrum (and especially the masses and mixings in the stop sector) plays a crucial part in the resulting lowest value accessible for tan β (∼ 5-6 in figure 3), as it controls the magnitude of the radiative corrections to the mass of the light Higgs doublet.

Light CP-even singlet
A quite natural NMSSM scenario, already noted for its admissible application to a (speculative) enhancement of the h[∼ 125 GeV] → 2γ decay rate [46][47][48][49][50][51][52], is that of a light CP-even singlet state with mass under ∼ 125 GeV. This setup offers several interesting phenomenological features: the presence of a Higgs state around ∼ 100 GeV could account for a small excess observed in the h → bb LEP data [107]; furthermore, the mixing between singlet and doublet (for non-vanishing λ) offers a lifting mechanism for the mass of the state identified with the LHC-observed signal; finally, the interplay with the singlet allows for an enhanced flexibility in the composition of the state at ∼ 125 GeV -the mixing matrix in the CP-even sector has now three mixing angles, instead of only one in the puredoublet case -so that small deviations from SM-like couplings might be interpreted in this fashion. In contrast to the prejudice according to which light states should already have left tracks in experimental searches, the presence of a light CP-even singlet proves phenomenologically viable, as the large singlet component entails a suppressed production cross-section of this state -via a suppressed coupling to SM-particles, e.g. gauge bosons or fermions -at colliders. Moreover, the singlet induces no major perturbation in the SM fermion and gauge-boson sectors. This scenario can be studied e.g. in the limit of a heavy decoupling SU(2) Higgs doublet.
In figure 4, we show such a region of the NMSSM parameter space, involving a heavy JHEP04(2016)095 Points with λ 2 + κ 2 > ∼ 0.7 2 are discarded by the scan as they would lead to Landau Poles below the GUT scale. Moreover, the region with 'large' λ and moderate κ tends to lead to unstable electroweak symmetry-breaking, as negative Higgs mass-squared are produced via the large singlet-doublet mixing (as soon as tan β > ∼ 5). The best-fit points, with χ 2 down to ∼ 75, involve a light singlet state: this fact is made evident when comparing the plots on both sides of figure 4, as the region including the best-fitting points (left-hand plot) largely coincides with that involving light singlets (right-hand plot). A determining factor for this correlation rests with the uplift of the mass of the light Higgs-doublet via the mixing effect (of only ∼ 1-2 GeV in the particular configuration of figure 4). Note that varying tan β (or the squark spectrum) displaces the favoured region in the {κ, λ} plane: indeed the magnitude of the mass contribution, which originates from the mixing among Higgs-states and shifts the mass of the light doublet state to a value closer to the center of the LHC signals, changes accordingly. Another reason for the improved fit values in the presence of a light singlet is associated with small deviations (at the percent level) from the standard values in the couplings of the light doublet to SM particles: the mixing with the singlet results in an increased flexibility of the doublet-composition of the state, which in turn allows for a possibly improved match with the measured signals.
The composition of the two lightest CP-even states in the scan of figure 4 is displayed in the upper part of figure 5: S ij denotes the orthogonal matrix rotating the CP-even Higgs sector from the gauge eigenstates -second index 'j'; j = 3 stands for the singlet JHEP04(2016)095 component -to the mass eigenbase -first index 'i'; the mass states are ordered with increasing mass. One observes that significant singlet-doublet mixing up to ∼ 20% can be reached in the vicinity of m h 0 1 ∼ 100 GeV, although best-fitting points show a mixing under ∼ 5%. This latter fact is related to the size of the mass-shift bringing the mass of the doublet-like state m h 0 2 in agreement with the window of the LHC signal (larger mixing would lead to m h 0 2 beyond the desirable ∼ 125 GeV range in the present configuration). This mass shift of the doublet state via its mixing with the light singlet, ∆m h 0 2 , is defined in the following fashion: regarding the heavy doublet sector as essentially decoupled, the squared-mass matrix of the singlet and light-doublet CP-even Higgs states may be approximated as the decoupled block: 3 Up to a sign, one can identify cos θ S S 13 , which determines the singlet-doublet mixing angle. In the presence of a lighter singlet-like state, the upward shift of the doublet state is defined as ). This quantity, still for the scan of figure 4, is shown in the lower right-hand portion of figure 5: while the uplift in mass may reach up to ∼ 8 GeV, shifts of only 1-2 GeV are favoured by the fit in this particular scan. Note that the formula that we have just derived only makes sense if h 0 2 is indeed the doublet-like state: for this reason, when displaying ∆m h 0 2 as a function of m h 0 1 , we cut the plot at m h 0 1 < 120 GeV, since, for m h 0 1 > ∼ 120 GeV, the doublet-like state is likely to become h 0 1 in order to match the LHC signals at ∼ 125 GeV. The plot on the lower left-hand side of figure 5 shows the squared-coupling of the singlet-like state h 0 1 to Z bosons -controlling the production cross-section at LEP; it essentially coincides with 1 − S 2 13 here: at m h 0 1 ∼ 100 GeV, this quantity reaches ∼ 5% (for best fits) up to ∼ 20% of its SM value: for memory, the ∼ 2.3σ LEP excess 4 in H → bb observed in this mass-range would be compatible with a Higgs-like state, the squared coupling strength of which is reduced to ∼ 10% of its SM value.
One can observe that (in this particular scan) the case where the lightest state is a doublet -represented by the limit S 2 13 → 0, m h 0 1 → 125 GeV -yields a slightly worse fit than the scenario with a lighter singlet: two factors are at work here. The first one is related to the value of the mass characterising the doublet ('would-be-observed') state in this limit: it typically reaches ∼ 121-123 GeV only -which lies on the margin of the uncertainty-allowed window. Note in particular that the mixing-effect tends to push the mass of the 'visible' state (now h 0 1 ) into the 'wrong direction' (to lower it) when the singlet is heavier. Yet, for some of the points under consideration, m h 0 1 reaches ∼ 125 GeV, hence evades this first argument: in this case, the main penalty with respect to the points involving a lighter singlet originates from the details of the couplings of the 'observed' 3 Note that we derive here an approximate formula, under the assumption that the heavy doublet state has negligible effect. This approach is qualitatively justified, at least at the level of the mass shift, provided MA MZ . A more exact expression, accounting for the heavy doublet state though losing somewhat in clarity, may be derived in a similar fashion however. 4 Local significance without taking into account the 'look-elsewhere'-effect.

JHEP04(2016)095
state to SM particles, hence of its production and decay rates at the LHC. While both configurations -with a lighter singlet or a lighter doublet -provide (doublet) couplings within a few percent of each other, and of those values that a SM Higgs boson at this mass would take, small deviations can provide a closer match to the LHC data. In our particular scan, for instance, the γγ rate is slightly enhanced when the lighter state is dominantly singlet, resulting in an improved agreement with the ATLAS measurement. More generally, the effects that the presence of a light singlet state may have on the couplings of the light doublet are related to the increased flexibility inherent in the 3 × 3 Higgs-mixing matrix S ij when compared to the case of a pure doublet 2 × 2 matrix. While one degree of freedom controls the singlet-composition of the Higgs state at ∼ 125 GeV -i.e. its 'invisible', for phenomenological reasons subdominant, component -the other two modulate the relative proportions amongst the two doublet components H u and H d -which are fixed in the case of pure doublets: in the limit M A M Z , the corresponding H d /H u ratio would be ∼ tan −1 β -therefore granting room for small deviations (or not, if the relative proportions are left unchanged) at the level of the couplings, with respect to the naive SM-like case. While this mechanism would offer an interpretation for slightly non-SM couplings, should this case be motivated by precision measurements of the Higgs properties, note that similar effects can also be obtained e.g. via loop-effects involving the supersymmetric spectrum. We shall come back to the question of non-standard couplings of the 'observed' state in the following section (5).
We complete this discussion with figure 6, whose scan differs from the previous one only by a lower value of the trilinear stop coupling |A t | -which thus tends to decrease the magnitude of the corrections to the mass of the light doublet originating in radiative effects. The situation is essentially comparable to the previous case, except for the fact that larger uplifts of the mass of the doublet-like h 0 2 -∆m h 0 2 ∼ 2-4 GeV -are now favoured. Larger singlet-doublet mixings, hence larger squared couplings of the light singlet to Z bosons, 1 − S 2 13 15-20%, are correspondingly preferred, in the vicinity of m h 0 1 ∼ 100 GeV. With slightly heavier singlets m h 0 1 110-115 GeV, we observe that large mixings, up to ∼ 25% may appear. Note that the best-fit point lies in an isolated region with mass close to ∼ 60 GeV: this isolated position results both from the limited scan density and from the marginal situation with respect to the LEP limits.
Finally, we note that, in the plots of figure 5 and 6, the mass of the singlet may reach values as low as ∼ 62 GeV without spoiling the quality of the fit. The case of states under ∼ 62 GeV opens the possibility of h 0 2 → 2h 0 1 decays and will be treated in a separate section (6). While precision tests at the level of the couplings of the Higgs state at ∼ 125 GeV could provide arguments in favour of this scenario involving a light singlet, as we already mentioned, the most convincing evidence would lie in the detection of the light singlet itself. The latter is likely to appear as a 'miniature' Higgs boson, i.e. with decay rates grossly comparable to those of a SM Higgs boson at the same mass but a reduced production crosssection (and a smaller width) -indeed the singlet component decouples from SM-particles and production is thus only achieved via a small doublet component. The observability of the singlet state thus critically rests with the magnitude of its doublet component. In the JHEP04(2016)095 discussion above, we stressed that the doublet composition (1 − S 2 13 ) of the light CP-even singlet may reach O(10%-20%) and it is likely that the corresponding signal would be large enough to allow detection -at least in the form of a local excess, while discovery at the JHEP04(2016)095 5σ level could remain challenging, provided the LHC searches are extrapolated to the lowmass region. Yet, apart from LHCb searches in the τ τ channel [108], which are still quite far from the necessary sensitivity in order to probe such a scenario, only the recent ATLAS results in the γγ channel [109] consider masses below ∼ 110 GeV. (Note that, in the scans of figure 5 and 6, the typical cross-section at 8 TeV in the diphoton channel lies below 1 fb; see also the discussion of figure 11 below.) On the other hand, smaller singlet-doublet mixings, at the percent level or below, would also fit the picture adequately, allowing for a mass-uplift of the doublet and / or small variations in its couplings to SM particles -such are actually the best-fit points of the scan in figure 4, which can be found in the appendix. The visibility of the light singlet in direct production should become increasingly difficult as its doublet component becomes small.
Alternative search strategies have been suggested, as the Higgs-to-Higgs couplings need not follow the same pattern as couplings to SM particles: singlet-doublet Higgs couplings could in principle allow for singlet production from e.g. the observed state, via Higgs-pair productions. It was stressed, however, that even such channels did not ensure the visibility of the light singlet, as the presence of this particle does not necessarily entail significantly larger inclusive rates than those of a single doublet state [78][79][80][81][82][83] (although an enhancement by a factor up to 2-3 has been reported for certain points). Let us highlight the fact that trilinear singlet-doublet couplings, while possibly as large as SM Higgs-to-Higgs couplings, 192 GeV, where we took m SM H = 125.6 GeV, v 174 GeV, may also remain much smaller without contradicting the light-singlet scenario. In figure 7, we display the strength of the triple Higgs couplings involving the light singlet and the light doublet in the scan of figure 4: while the trilinear coupling of the state h 2 with mass ∼ 125 GeV remains SM-like, the couplings involving both singlet and doublet-like states reach only up to ∼ 30% of the SM-strength. Assuming that the cross section for pair production follows a similar pattern -which is only justified at high center of mass energy -we see that a discovery in such channels would be challenging experimentally. On the other hand, pair production close to threshold is very sensitive to a small imbalance among triangle and box contributions so that the presence of a light singlet may affect this observable. An estimate of such effects goes beyond the scope of this paper, and we refer the reader to the discussions in [78][79][80][81][82][83]. As a summary, let us stress that, if Higgs-pair production is a viable search channel for light singlet states, should this state be present in nature, it does not automatically ensure the discovery of the singlet at the LHC. Production associated to Higgs-gauge couplings was also discussed in [91] but again depends critically on the magnitude of the doublet component of the light singlet-like state.
Other production modes would involve the supersymmetric spectrum or the heavy Higgs states. In particular, the decays h 0 In this respect, NMSSM effects may intervene at several levels: • The doublet-to-doublet Higgs couplings differ from their MSSM equivalent due to NMSSM-specific terms in the Higgs potential (∝ λ, κ) so that the associated width may differ significantly.
• Singlet-doublet couplings induce a decay of the heavy Higgs into the singlet-like state. 192 GeV (see text).

JHEP04(2016)095
• Kinematically accessible additional final states (e.g. singlinos) also affect the branching ratios of h 0 3 . More generally, the branching ratios are strongly dependent on the details of the supersymmetric spectrum.
• The decay rates (into e.g. bb, γγ) of the decay products h 0 1 , h 0 2 may vary with respect to the naive standard rates.
Unless the decay products h 0 can be observed separately (using kinematical cuts), it would thus be difficult to infer the presence of a singlet from the inclusive rates, as several factors could explain a deviation from the predictions of a type II 2HDM.
Another remark on the decay rates of the light singlet is in order: as mentioned above, one may naively expect them to coincide (coarsely) with those of a SM Higgs boson at the same mass. Yet, it has also been stressed that these rates could show unconventional behaviours in specific cases, at low [110] or large tan β [73][74][75]: couplings to down-type quarks could be suppressed indeed, which would lead to an apparent enhancement of decay channels such as cc or γγ. Extreme cases with an up-to-seven-times enhanced diphoton branching fraction of the singlet, allowing for cross-sections at the level of their SM equivalent -despite the reduced production cross-section -have received much attention in view of their remarkable consequences on the possible discovery of such a singlet. We illustrate this possibility of non-conventional singlet rates with figure 8: with an intermediate value of tan β = 12, we observe that the bb branching fraction may be strongly suppressed, while the other rates (here cc and γγ) are enhanced, together with acceptable fit valuesnote however that the best-fit points lie in a region of more SM-like behaviour.

JHEP04(2016)095
Note that points involving light singlets are quite common in the NMSSM parameter space. The only difficulty consists in stabilizing the low singlet mass and keeping the singlet-doublet mixing under control: too strong a mixing would push the squared mass of the lightest state towards negative values. The typical scale entering the singlet mass is κ λ µ, so that light singlets favour low ratios κ/λ. As tan β increases, however, the balance among terms entering the mixing of the light doublet and singlet CP-even states is disturbed, such that the region with large λ and low κ becomes increasingly unstable. For larger values of κ/λ, one observes that µ tends to be driven to low values ∼ 100 GeV -in order to keep the singlet mass at the electroweak scale without relying too much on accidental cancellations. This issue becomes even more severe when λ ∼ κ becomes large: one then relies exclusively on the accidental cancellation in the singlet diagonal and the singlet-doublet mixing massmatrix entries. It is therefore most natural to consider the scenario involving a light singlet in the low tan β regime, allowing for small κ/λ: this is the focus of the next section.

Low tan β and large λ
The region with low tan β (∼ 2) and large λ (∼ 0.6-0.7) is particularly interesting in the NMSSM parameter space: as we just mentioned, light CP-even singlets under ∼ 125 GeV appear most naturally there; furthermore, the squared mass of the light-doublet Higgs state receives an F-term contribution of the form λ 2 v 2 sin 2 2β at tree-level (see e.g. eq. (2.23)

JHEP04(2016)095
of [11]), which is maximized in this regime, so that significantly smaller radiative corrections than in the MSSM case are needed to reach the experimental value of the observed signal; finally, the region of low tan β is 'specific' to the NMSSM, in the sense that current phenomenological requirements on the Higgs sector forbid it in the MSSM. One may typically consider this setup in the Peccei-Quinn (PQ) limit κ/λ 1 -note that, in view of the Landau Pole constraints, turning to relatively large values of λ automatically implies moderate values of κ anyway.
Given that, in this context, one does not have to rely on large corrections from the stop sector to uplift the mass of the light doublet Higgs state into the mass range ∼ 125 GeV where it can be identified with the LHC-observed signal, we will consider moderate stop masses (∼ 0.5 TeV) and trilinear couplings (∼ 0.1 TeV). While we do not check the compatibility of this choice with LHC limits on supersymmetric searches, note that this is not a binding requirement but purely an illustration of the fact that radiative corrections to the Higgs mass are of lesser importance in this region of the parameter space: the specific tree-level contribution and / or the mixing with a light singlet are mechanisms enough to generate the mass of the light doublet Higgs state in the ∼ 125 GeV range. In fact, if in this regime the SUSY sector is such that it generates large radiative corrections to the mass of the light Higgs doublet state (e.g. via a sizable mixing in the stop sector), this would lead to predictions for the mass of the light doublet state in the region of low tan β and large λ that tend to be higher than the mass value detected at the LHC -up to ∼ 140 GeV if it is the lightest CP-even state and beyond if it is the second lightest: see e.g. [17]. Yet, a Higgs mass that is compatible with the observed value can still be obtained even in the presence of large contributions from both the NMSSM tree-level and SUSY radiative corrections through the effect of the mixing between the doublet and the singlet. This mixing effect lowers the mass of the light CP-even doublet, provided the singlet is heavier. Note finally that B-physics is of almost no concern in this scenario since both charged-Higgs and supersymmetric effects remain small, due to a heavy H ± (decoupling limit) and low tan β (absence of an enhancing-factor in radiative effects). On the other hand, low values of tan β tend to suppress supersymmetric contributions to (g − 2) µ , translating into a typical pull of a few units in χ 2 . The presence of light sleptons / charginos / neutralinos (with mass close to 100 GeV) could balance this effect, however, so that we will assume low masses for these states in the following. Figure 9 shows the χ 2 distribution in a scan where tan β = 2, λ = 0.7, κ = 0.1, and the supersymmetric spectrum is relatively light (at least for the third generation of sfermions, with small trilinear couplings). The plots illustrate how the fit distributes in terms of the mass and singlet composition -i.e. singlet component squared S 2 i3 -of the first and second lightest CP-even Higgs states. Points where the lightest Higgs state has a dominantly singlet nature (S 2 13 > 0.5) are shown in yellow-red shades. Bluish tones correspond to points where the lightest state is dominantly doublet (S 2 13 < 0.5). Both configurations give an excellent agreement with the measurements reported by the LHC and the TeVatron, and their fit values improve somewhat on the SM limit. One observes that the best fit points (with χ 2 ∼ 76-79) tend to cluster in the vicinity of S 2 i3 = 0 or 1, that is for moderate singlet-doublet mixing. This seems reasonable since one expects a JHEP04(2016)095 'full' doublet state at ∼ 125 GeV: the corresponding experimental signals would have been suppressed in proportion to the singlet composition otherwise. Still, fairly large values of the mixing (S 2 i3 ∼ 0.5) turn out to be a possible scenario, giving acceptable χ 2 ∼ 77-80: this situation occurs only when both states are almost degenerate and within a few GeV of ∼ 125 GeV. This last configuration will be studied in more detail in section 7.
In figure 10, we display some information concerning the coupling properties of the light states to SM particles obtained with the scan of figure 9: the plots on the left-hand side illustrate the possible effects that the three-state mixing could have on the couplings of the doublet-like state -the one which is identified with the LHC-observed signal. The plot on the upper part shows the proportion of H d /H u components in this state as a function of the mixing: deviations of up to ∼ 25% appear relative to the naive ratio 0.5 = tan −1 β (expected in the decoupling limit), and the maximal effects are obtained for maximal mixing S 2 13 ∼ 0.5. Deviations at the level of the branching ratios reach ∼ 25% for bb and ∼ 50% for γγ. Note that best fit regions favour more moderate variations however. Note also that, even for large mixings, the deviations in doublet proportions or in branching ratio may remain negligible.
The plots on the right-hand side of figure 10 provide similar information but concerning the singlet-like state. We observe that the corresponding coupling properties may depart very significantly from the naive behaviour when the state is almost purely singlet -that is, in the vicinity of S 2 13 ∼ 0 and ∼ 1 -and that the corresponding points offer excellent fit to the Higgs data. The bb channel may then become subdominant while the diphoton channel is enhanced by a factor of up to ∼ 7. In such a spectacular case, the light singlet JHEP04(2016)095 Figure 10. Same scan as in figure 9: consequences of the three-state mixing on the couplings of the light states are shown. The plots on the left concern the mostly-doublet state (identified with the observed signal at ∼ 125 GeV); those on the right give information concerning the lighter (yellow-red points) or heavier (bluish points) singlet. The first pair of plots displays the proportion of H d /H u components (in comparison to the value tan −1 β = 0.5 expected in the decoupling limit). The branching ratios into bb and γγ are provided in the lower part of the figure.
JHEP04(2016)095 Figure 11. On the left: gluon-gluon-fusion cross-section for the mostly-singlet state, then decaying into a pair of photons, for a center of mass energy of 8 TeV, resulting from the scan of figure 9; the corresponding value for a SM Higgs boson is given by the green curve. On the right, a reproduction of the ATLAS limit on the fiducial cross-section for a light Higgs state (in the presence of the ∼ 125 GeV one) decaying into photons.
could be more easily observed in direct production at the LHC. On the other hand, the fit tends to associate the large diphoton branching fraction tightly with the limit of a pure singlet state, that is with vanishing production cross-sections. Note also that the naive scenario of a singlet-like state with dominant decays towards down-type fermions is also represented and actually provides the best-fit points of the scan. Unconventional decay rates also appear as a possibility when the singlets are beyond ∼ 125 GeV (blue points), even though the maximal diphoton rates remain below ∼ 1%.
In figure 11, we study how the Higgs production cross-section at 8 TeV compares to the ATLAS limits on the fiducial cross-section for the diphoton decay channel [109]. We estimated the cross-section for the light Higgs states of the scan of figure 9 in the following way: we multiplied the SM gluon-gluon-fusion cross-section delivered by SusHi [111,112] by the squared effective coupling of h 0 1 to gluons, relative to its SM value at the same mass, and the diphoton branching ratio of h 0 1 . We observe that the cross-section may almost reach the order of magnitude probed experimentally, both when the singlet is heavier or lighter than 125 GeV (note that in the immediate vicinity of 125 GeV, comparing the cross-section of the mostly-singlet state with the ATLAS limit has limited sense, due to the possibly large mixing between singlet and doublet states), although the best-fitting points tend to cluster around much smaller values -at or below the 1 fb range. Further searches in the low-mass region, in the diphoton but also in the fermionic channels, would be an interesting probe and place limits on the light-singlet scenario.
In figure 12, we vary tan β and λ somewhat so as to modulate the strength of the Fterm contribution to the tree-level doublet Higgs mass. As a result, larger singlet-doublet mixings are favoured: the two-state mixing uplift can indeed compensate the decreased tree-level contribution and thus help maintain the mass of the light doublet state in the vicinity of ∼ 125 GeV. In agreement with our discussion in section 4, we observe that large singlet-doublet mixing, up to ∼ 25%, may be achieved for a singlet mass in the range [90][91][92][93][94][95][96][97][98][99][100] GeV, with excellent fit-values to the Higgs measurement data. Therefore, this low JHEP04(2016)095 Concerning the prospects of discovery of the light state in pair production, the Higgsto-Higgs couplings in the scan of figure 12 are displayed on the right-hand side of this figure.
The typical magnitude is about 10-40% of g SM simply estimated as the latter depends on several interfering diagrams. We see however that the typical couplings reach ∼ 30% of the pure-doublet value.
Although all these observations are essentially similar to our discussion in section 4, the crucial point rests upon the fact that such a Higgs phenomenology is also achievable in this low tan β / large λ regime, without relying on large radiative corrections to the Higgs masses. This provides a motivation for relatively-light supersymmetric spectra (at least, as far as the third generation is concerned). In the case where the search for stops at the LHC would provide experimental support for such a configuration, deviations of the Higgs couplings from the SM expectations could be generated at the loop level and be addressed in precision tests.

Light CP-odd (even) Higgs states under m h[125] /2
A durably-considered NMSSM scenario (see for instance [113,114]) is that involving light neutral states, e.g. a light CP-even singlet or a light CP-odd state, allowing for unconventional Higgs decays (h → 2A 1 /2h 1 ). As the width of a SM-like Higgs is quite narrow -∼ 4 MeV at ∼ 125 GeV -Higgs-to-Higgs decays can easily dominate the width of a Higgs state as soon as they are kinematically allowed. This configuration could have explained the invisibility of a doublet Higgs state h at LEP -provided m h > ∼ 85 GeV (i.e. above the limit from the decay-mode independent search [115]) and A 1 → τ + τ − would represent the dominant decay channel of the pseudoscalar. Furthermore, this mechanism suggested an interpretation of the 2.3 σ (local) excess in e + e − → Zh → bb originating from suppressed branching ratios into SM particles of a CP-even state at ∼ 100 GeV, instead of a suppressed production as in the case of a light singlet, and would have led to a reduced visibility of the conventional channels at the LHC -even if the light CP-even doublet had been heavier than 100 GeV. Considering the success of these conventional searches and the approximately SM-like behaviour of the LHC-observed state, Higgs-to-Higgs decays of the SM-like state no longer appear as a favoured option, unless one would correspondingly enhance the SM-like modes, which is not easily achieved in the context of the NMSSM. Yet, we already noted in the previous sections that the existence of light states -under the kinematic threshold at ∼ 62 GeV -was still possible (see e.g. figure 5); we will discuss here under what conditions. In brief, if the invisible decay channel is not kinematically closed, compatibility with the measured rates of the observed signal demands the decoupling of the light state -suggesting a singlet-like nature with reduced Higgs-to-Higgs couplings [84][85][86][87]. Figure 13 shows how the fit distributes in the presence of light CP-odd Higgs states. Note that the SUSY sector, in particular the mixing in the stop sector, is set to provide radiative corrections to the mass of the light doublet Higgs of the appropriate magnitude in view of the signal at ∼ 125 GeV, as already discussed in section 3. While tan β is left free, larger values > ∼ 10 are accordingly favoured by the requirement of a SM-like Higgs state at ∼ 125 GeV, but also by (g − 2) µ (for which the exact prefered range is determined by our choice of spectrum in the slepton, neutralino and chargino sectors). One observes in the plot on the top left-hand quadrant that, when the CP-odd Higgs is lighter than half of JHEP04(2016)095 the mass of the 'observed' Higgs, its doublet composition 5 must fall under ∼ 1% to provide an acceptable fit to the data, while fairly larger singlet-doublet mixings are allowed beyond this mass. This condition can be interpreted as follows: if the channel h[∼ 125 GeV] → 2A 1 is not kinematically forbidden and A 1 has a significant doublet component, doublet Higgs couplings lead to a large (dominant) decay of the state h[∼ 125 GeV] -which is identified with the LHC Higgs discovery -into a pair of pseudoscalars; yet, this is in contradiction with the roughly standard behaviour of the signal observed at the LHC, hence a disfavoured possibility; it thus follows that the light pseudoscalar with mass under the threshold must JHEP04(2016)095 be essentially singlet in nature, for its presence not to spoil the fit to LHC Higgs data. Note, on the other hand, that for masses beyond 62 GeV, the doublet composition of the light CP-odd state can reach up to ∼ 100%, as we will discuss in the light-doublet section (8). However, apart from this rather isolated possibility, LHC limits -in particular searches in the τ τ channel -do not leave much room for heavy Higgs doublets under ∼ 400 GeV (see e.g.  In view of the early LHC results, the maximal branching ratio into lighter Higgs states of the state at ∼ 125 GeV which remains compatible with the observed signals had been estimated at ∼ 20% [84][85][86][87]116]. In the bottom left-hand quadrant of figure 13, the unconventional branching ratio BR(h[∼ 125 GeV] → 2A 1 ) is displayed as a function of the Higgs mass. Best-fit (yellow-to-red) points cluster within BR(h[∼ 125 GeV] → 2A 1 ) < 20% indeed. Interestingly, a light CP-odd Higgs may coexist with a light CP-even singlet-like state at reduced κ/λ: the CP-even singlet may then decay dominantly into 2A 1 , as shown in the plot at the bottom right-hand corner of figure 13. In this configuration, the observability of such a CP-even singlet state is hindered both by its reduced production cross-section and its unconventional decay. However, this setup typically requires that BR(h[∼ 125 GeV] → 2A 1 ) be suppressed through an accidental cancellation of the h[∼ 125 GeV] − A 1 − A 1 coupling.
The case of light CP-even states with mass under the decay threshold of the observed state essentially follows the same principles: from figures 5, 6, 9 and 12, we observe that light CP-even singlets under ∼ 62 GeV do not seem to offer particularly good fits to the Higgs-measurement data, as points with acceptable fit values drastically disappear from the low-mass tail under ∼ 62 GeV. Note that in the corresponding plots, λ tended to be large, so that the decay h 2 → 2h 1 would typically retain a significant branching fraction (which is disfavoured by the fit), even though h 1 is purely singlet: this explains why no red points persist at S 2 13 1 below 62 GeV in these plots, contrarily to the case of P 2 13 1 in figure 13. An example for this scenario can be found in table 4 (Point 9). In figure 8, on the contrary, we observe points with mass down to 40 GeV, some of them having an excellent fit value. The fact that these points are dominantly singlet was already a requirement from LEP limits. However, in figure 15 the threshold.
From this discussion, it follows that the success of the conventional Higgs searches prohibits 6 possible light Higgs states below ∼ 62 GeV to play a significant part in the phenomenology of the state which was discovered by the LHC, but also to intervene in connection to other standard particles: we have observed indeed that a condition for suppressed Higgs-to-Higgs decays laid in an almost-pure singlet nature of the light Higgs state. Therefore, the latter does not possess relevant couplings to SM-fermions or gauge bosonsthese develop only in proportion to the doublet components. As a consequence, direct production of the light state, e.g. via its coupling to a quark line (in associated production with tops or b's), proves even less promising than searches in Higgs-to-Higgs processes: the corresponding cross-sections would indeed receive suppression factors of the form ∼ (1−S 2 13 ) / (1 − P 2 13 ) for a CP-even / CP-odd state. 7 On the other hand, the existence of such avoidant states remains a phenomenological possibility. The most straightforward search channel remains a possible decay (with branching ratio of a few percent) of the state observed at ∼ 125 GeV towards the light (pseudo)scalars, which may result in pairs of bb's, ττ 's or γγ's in the final state. Note that these remarks are consistent with the results of the recent analysis [88,89]. Other search channels would involve new Higgs or supersymmetric particles, depending critically on the characteristics of these states.
As a side-remark, let us mention that another type of unconventional Higgs decays would involve light supersymmetric particles, e.g. a light singlino. This scenario has raised interest in dark matter studies [104,105]. While our framework is not suited to discuss dark matter phenomenology, we may still comment on the tentative presence of a light neutralino with mass under ∼ 62 GeV from the point of view of Higgs physics. Decays 6 Note that this conclusion holds for the NMSSM only. Although this is likely to prove a common trend in most models where such light (pseudo)scalars are directly associated with the Higgs sector, there also exist theoretical frameworks where the light state has naturally suppressed couplings to the SM-like Higgs boson, without necessarily decoupling from SM-particles. 7 Extreme values of tan β may partially compensate this suppression in the couplings to down-type quarks and leptons, however. of the Higgs state which is identified with the signal observed at the LHC into a pair of light neutralinos (invisible decay) would indeed lead to a suppression of the rates in the conventional search channels, which contradicts the LHC / TeVatron results (unless this effect would be compensated by an enhanced production mode). In the case of a light singlino, the corresponding width is related to λ 2 , so that large values of λ are again disfavoured as long as this decay channel is kinematically open. Figure 16 in the low tan β / large λ regime illustrates this comment: the region below the threshold for h[125 GeV] → 2χ 0 1 gives rise to rather large χ 2 values. A comparable limit would apply in the presence of light binos: the couplings to a light doublet Higgs would then be of electroweak magnitude, which could lead to a disfavoured decay of the Higgs state towards neutralinos. However, contrarily to the case of light scalar states, the constraints on light new fermions are typically milder since such decays are of a similar type to the standard channels, hence less likely to dominate the branching ratio as soon as the kinematical threshold opens.

Two states in the vicinity of ∼ 125 GeV
In [76,77] the possible presence of two CP-even Higgs states in the vicinity of ∼ 125 GeV was highlighted, the rates of both states adding while the experimental resolution remained too broad to distinguish between them. Note that, while individual channels -e.g. searches in the diphoton final state -provide an excellent precision on the mass where the Higgs signals are centered, it will be challenging for ATLAS and CMS, because of the limited detector resolution, to resolve two Higgs signals separated by less than 2-3 GeV. On the theoretical side, a typical separation scale is the SM Higgs width of 4 MeV. The aim in [76,77] originally consisted in exploiting the presence of two states in the signal region, so as to enhance the diphoton and ZZ rates in the context of a NUHM version of the NMSSM -thence explaining tentative deviations of the apparent Higgs rates from the standard values via this approximate degeneracy of Higgs states. In this section, we shall discuss such configurations in more detail.
While the CP-even singlet and the light doublet may be close in mass, requiring them to lie within a few 100 MeV / 1 GeV demands that the mixing entry in the mass matrix JHEP04(2016)095 approximately vanishes. This is a tuning requirement, in general, unless one moves to the MSSM limit of the NMSSM (and it is then still necessary to adjust the diagonal entries of the mass matrix so that they are approximately equal). As a consequence of the small size of this off-diagonal mass entry, both states can be quasi degenerate while their composition, hence their coupling properties to SM particles, may vary from that of pure singlet / doublet to strong admixtures. Note that the typical widths involved ( < ∼ 4 MeV) are still very small compared to the mass differences which we are considering, so that interference effects such as those discussed in [117] should remain small. Therefore, as long as the two single states are not resolved, the overall behaviour in interactions with SM particles would naively coincide with that of the doublet-like state taken alone. Only the apparent Higgs-to-Higgs couplings may show deviations -to which even a HL-LHC is unlikely to be sensitive to, due to the limited experimental precision, and which would still have to be disentangled from other sources of non-standard Higgs-to-Higgs couplings for the state at ∼ 125 GeV. On the other hand, the composition of the apparent light doublet may differ slightly from that of a pure doublet state due to the three-state mixing effects (similarly to what we discussed in sections 4 and 5): here again, this effect could account for small deviations from the standard rates. On the side of the singlet, the simplest case consists in having all particles to which it directly couples (CP-odd singlet, higgsinos) heavier than half its mass: then, the only possible decay products are SM particles, the corresponding widths depending exclusively on the doublet component. In this context, the singlet does not contribute to the total width. If, on the contrary, significant decays of the singlet towards other new-physics states are allowed, then we are back to a case similar to what we discussed in section 6: indeed, the singlet-doublet mixing would then dilute the branching ratios into standard particles, and could hence contradict the signals observed at LHC -large singlet-doublet mixings would then be disfavoured. This second configuration will not emerge from the scans.
We have already come across points involving quasi-degenerate CP-even Higgs states, e.g. in section 5: some examples have been recorded in table 5 below. As a first illustration of this scenario, we consider the MSSM limit in figure 17: the off-diagonal mass term between singlet and light doublet state is then naturally suppressed. One observes, however, that significant singlet-doublet mixing, up to 50% can develop for favourable fit values, provided the two states are almost degenerate. On the other hand, no effect at the level of the rates is expected in the MSSM limit (because of the decoupling of the singlet). We display the apparent branching ratio into bb and γγ as a function of the singlet-doublet mixing on the right-hand side of figure 17, for points where m h 0 2 − m h 0 1 < 1 GeV: while some variations are present, these are entirely independent of the mixing or the presence of the singlet, and their origin can actually be traced back to the scale of the heavy doublet sector (M A ) and the value of tan β, hence appears as a pure doublet effect. Similarly, consequences on Higgs pair production are minimal since the singlet-doublet interactions are suppressed in this limit. While singlet-singlet interactions formally enter the relevant Higgs-to-Higgs couplings, the associated effects are projected onto the doublet component -due to the coupling to SM particles in the process of the Higgs production -and vanish when summed over both degenerate states. In the lower left-hand quadrant of figure 17 JHEP04(2016)095 one observes indeed that the apparent trilinear Higgs coupling accounting for both states in the vicinity of ∼ 125 GeV remains SM-like. Therefore, in this scenario, only the presence of two separate peaks in the Higgs searches -as the result of a large mixing between the singlet and the doublet state -with widths adding to ∼ 4 MeV, would be observable and document the quasi-degeneracy of the singlet with the light doublet. On the other hand, the singlet could also be almost degenerate but remain (almost) unmixed, in which case, its 'peak' would remain invisible. The same configuration is considered in figure 18, but in a different regime: we turn to the low tan β / large λ region again. Only points where the two lightest CP-even states have masses within 1 GeV of each other are kept in this scan. Deviations at the percent level can be observed in the apparent rates into SM-particles. The corresponding effect is obviously associated with the three-state mixing -this is most visible in the case of the γγ channel -in contrast to the case of the MSSM limit. Note that pure doublet effects are negligible here as M A M Z and tan β is fixed. Large singlet-doublet couplings affecting the apparent triple-Higgs coupling at ∼ 125 GeV may develop as well, leading to a ∼ 10% increase (at most) of the apparent trilinear Higgs coupling (see the lower left-hand quadrant). Note that while this kind of scenario may occur in different parameter regimes, the existence of a doublet state at ∼ 125 GeV always relies on the mechanisms described in the previous sections, i.e. substantial radiative corrections driven by the SUSY spectrum and larger tan β or the specific tree-level contribution associated with large λ and low tan β. On the other hand, singlet-doublet mixing does not provide a mass uplift for the doublet in the scenario under consideration.

JHEP04(2016)095
Another possibility would consist in the presence of a CP-odd state in mass-proximity to the SM-like state. Such a scenario might even be considered as a CP-conserving approximation of the CP-violating case, where the observed state would appear as a superposition of CP-even and CP-odd components. Considering LHC limits associated with Higgs searches in the τ τ channel as well as (indirectly) with top decays to a charged Higgs, it proves difficult to admit a (pure) CP-odd doublet state in the desired mass range. The possibility of a CP-odd Higgs close to ∼ 125 GeV thus rests with singlet-like states, which may however carry a significant doublet component. In figure 13, one observes that a doublet composition of 5-10% still receives an acceptable fit to the data in the vicinity of ∼ 125 GeV. The fermion rates at ∼ 125 GeV would be the observables that are most significantly affected in this setup, with an apparent increase of a few percent, while also the diphoton channel may receive a small subsidiary contribution.
Note that in the CP-even as in the CP-odd case, the additional state may well be purely singlet so that it would remain undetected in direct production due to vanishing couplings to SM particles. Only at the level of pair productions, making use of Higgs-

JHEP04(2016)095
to-Higgs couplings, would there be a deviation from the naive one-particle case -which would be extremely challenging to detect with sufficient precision experimentally. Note also that, e.g. in the MSSM limit, a singlet that is mass degenerate with the state at ∼ 125 GeV could even have no impact at all on the rates or on multi-Higgs production.

Two-light-doublet scenario
Identifying the LHC signal with a heavy doublet state -which suggests that a lighter Higgs doublet state has eluded searches so far -is a scenario which was originally considered in the MSSM [27][28][29][30][31][32]. It was recently put under further pressure by the publication of charged Higgs searches in top decays [34] -see the newest reference in [36][37][38][39][40] for a discussion in the context of the MSSM -and [35], however. In this section, we will discuss the situation in the context of the NMSSM.
The presence of a light CP-even Higgs -with mass below 125 GeV -in the spectrum does not entail major phenomenological difficulties by itself. Indeed we have already outlined two possible strategies which allow such a scenario to evade LEP limits: • the production cross-section at LEP could be suppressed by a small coupling of the light Higgs state to electroweak gauge bosons -we have discussed above how a singlet Higgs naturally fulfills this requirement; • the decays in the standard search channels could be blurred by large non-conventional decays -in this case, the difficulty lies in explaining why these unconventional decays do not affect the observed state at ∼ 125 GeV.
While the success of the LHC Higgs searches makes the second approach difficult to implement in a realistic model, the first strategy remains viable, even in the case of a doublet state -even though the cancellation of the couplings of the light Higgs to electroweak gauge bosons is then largely accidental. Turning to direct limits from the LHC, it is to be noted that searches in the τ τ channel for an additional Higgs boson essentially disfavour neutral states with mass above ∼ 110 GeV, while the light doublet can be well below 100 GeV. In practice, we could indeed find NMSSM parameter points where a light doublet Higgs with vanishing couplings to electroweak gauge bosons -typically in the 70 GeV mass range -escapes all direct limits. Yet, if all the CP-even doublet states are at or below 125 GeV, the correlations of the masses in the doublet sector will force the other doublet masses -those of the CP-odd and charged states -to be close. Indeed, as a good approximation, we have: In this list of (approximate) masses for the doublet states, M A,eff does not exactly correspond to the NMSSMTools input M A , but represents a corrected value absorbing radiative corrections. There are, of course, additional deviations at the loop level, and δ rad only denotes the bulk of the large shift due to top/stop effects -hence associated to the H u flavour.
The CP-odd doublet Higgs causes limited concern as its mass typically falls under 100 GeV and its presence was tested at LEP only via pair production processes (since a CP-odd Higgs has vanishing couplings to gauge bosons). Only when its mass falls below the threshold ∼ 125/2 GeV does it entail indirect limits from the LHC measurement -in accordance with our discussion in section 6.
On the other hand, the existence of a light charged Higgs is problematic in view of the existing limits. This possibility is already severely constrained in view of the tensions that it produces in the B-sector -e.g. in the channels B → X s γ, B + → τ ν τ or B s → µ + µ − . Yet, to the price of a χ 2 pull of at least ∼ 4, these flavour limits could be satisfied. On the other hand, a light charged Higgs opens the possibility of top decays: these channels have been investigated by ATLAS [34] and CMS [35], e.g. in the bτ ν τ final state, and the current constraints exclude all the points that were preferred by the fit performed in figure 19.
In order to illustrate the workings of the recent LHC constraints on t → bH + , we turn to older versions of HiggsBounds results of the summer 2014 had not yet been included. To set a reference, the best-fit point in the SM-limit then receives a χ 2 of 8 ∼ 94. We then consider the region in the NMSSM parameter space where M A = 130 GeV is set as a fixed input. The results are shown in the {κ, λ} plane in figure 19. The best-fitting points are obtained for λ ∼ 0.2. They have a somewhat better χ 2 than the SM limit -χ 2 ∼ 88.0 -and provide the following spectrum: m h 2 125.0 GeV, m H ± 107 GeV and light H d -like doublet states (CP-even and odd) around 70-75 GeV (i.e. M A,eff ∼ 70 GeV); the singlet states are much heavier (beyond 1 TeV) and play no role in the electroweak Higgs phenomenology. The h 2 -rates are essentially SM-like, which explains the quality of the fit. On the right part of figure 19, one observes that the squared coupling of the light CP-even doublet state to Z-bosons is reduced (of the order 10 −2 ), while its typical production cross-section in the 8 TeV run of the LHC, accompanied by a diphoton decay, would amount to about a few fb's only.
Version 4.1.3 of HiggsBounds included the ATLAS limits on top decays [34], and we show on the upper part of figure 20 how the points of the scan in figure 19 compare with these constraints -note that the corresponding experimental bounds have been obtained under the assumption of a 100% H + → τ + ν τ decay, so that we rescale our points by a factor BR(H + → τ + ν τ ): while sitting on the edge of the exclusion limit, the light doublet scenario appears compatible with these constraints. On the other hand, the more recent CMS limits on t → H + b [35] had not been included within HiggsBounds, and we display them in the lower part of figure 20. All the best-fitting region is excluded while the few remaining points at m H ± 130 GeV -disconnected from the best-fit region due to the interplay of constraints -do not offer interesting fit qualities. Therefore, while the light doublet scenario might still be realized in regions with reduced BR(t → H + b), the current limits discard most of the associated parameter space.

Highlight of specific points
In this section, we focus on specific points in the parameter space, found in the vicinity of the best-fit points of the various plots presented thus far (see table 5, in the appendix), and aim at discussing the associated Higgs phenomenology in more detail. Note again that, although the supersymmetric spectra have some impact on the phenomenology -Higgs mass, (g − 2) µ , B-physics, etc., they are not strictly tied to the respective scenarios -i.e. analogous Higgs properties should be accessible with different, e.g. heavier SUSY spectra: the corresponding characteristics are thus purely illustrative and we shall not discuss the prospects of discovering such states at the LHC. Table 3 provides the NMSSMTools inputs and the Higgs spectra of the considered points, while table 2 gives the Higgs couplings to SM particles. While most of the qualitative features that we discussed in the previous sections enter the characteristics of the few points presented below, we would like to caution the reader against reducing the discussed phenomenology to those specific points.
The first point is characteristic of the decoupling limit: the heavy-doublet states are at ∼ 1 TeV, so as to decouple from the SM-like Higgs, while tan β 1 to maximize the JHEP04(2016)095  figure 19; the limits from ATLAS [34] (upper plot) and CMS [35] (lower plot) on this branching ratio BR(t → H + b) -assuming a 100% decay H + → τ + ν τ -are reproduced for comparison.
tree-level contribution to the mass of this light state (we stress that the rather large value tan β = 22.5 is actually driven by the anomalous magnetic moment of the muon). While λ and κ do not vanish, which means that one is far from the MSSM limit, the impact of the singlet states on the light doublet Higgs is negligible, and the CP-even singlet is actually quite heavy (∼ 1200 GeV). Correspondingly the light doublet state, at 125.0 GeV, is SM-like, with couplings to SM-particles within a few percent of their values at the same mass for a genuine SM Higgs boson: the corresponding rates at the LHC are thus comparable (still within a few percent). Distinguishing this point from the SM via the properties of the 'observed Higgs' would thus require a high precision on the couplings, achievable most likely at a Linear Collider only -note also that one may get even closer to a SM Higgs scenario via further decoupling of the heavy states. The heavy doublet JHEP04(2016)095 h 3 W W 6 · 10 −5 6 · 10 −6 3 · 10 −5 3 · 10 −7 8 · 10 −7 2 · 10 −6 h 3 ZZ 6 · 10 −5 6 · 10 −6 3 · 10 −5 3 · 10 −7 8 · 10 −7 2 · 10 −6

JHEP04(2016)095
states are essentially H d in nature, hence give very large couplings (enhanced by the large value of tan β) to bottoms and taus (for the neutral states). All other couplings to SM particles are suppressed (which is characteristic of the large tan β regime), so that decays into bottoms and taus dominate the branching ratio of the heavy neutral states. Typical production modes would be gluon-gluon fusion or production in association with b's. Note that, even though LHC searches at 13/14 TeV could likely probe such heavy states, e.g. in the τ + τ − channel, even larger masses of the heavy doublet states would be easily possible in this scenario. Moreover, should the heavy doublet states be discovered, the question of distinguishing such a point, in the Higgs sector, from a MSSM scenario remains open: the heavy CP-even singlet, lightly coupled to SM-particles, is likely to evade detection. As an additional feature for this point -additional but not binding: light CP-odd singlets are not a generic feature of this decoupling scenario -however, one observes that the CP-odd singlet is quite light, ∼ 110 GeV. Although its large singlet component of ∼ 99.7% leads to a significant suppression of the couplings of this state to SM particles (hence of its production cross-section), the couplings to down-type fermions are again enhanced by the large value of tan β, so that one may look for a signal in the τ + τ − and bb channels -in associated production with b's. The signal in the γγ channel would be strongly suppressed so that the current ATLAS limits [109] have no impact for this point. Remember also from figure 13 that CP-odd singlets between 63 and ∼ 150 GeV with slightly larger doublet components are also possible, which would improve their observability, although this feature would typically be associated with a lowered mass for the heavy Higgs-doublet states. Another possibility lies in exploiting triple-Higgs couplings for a production of A 1 in pairs. Furthermore, the associated production of A 1 with a Z-boson only proceeds via the tiny doublet component of A 1 , hence is also suppressed. The second point is representative of the light-singlet scenario of the NMSSM. Much that has been said in connection with the previous point, especially concerning the heavy doublet states, remains valid -note however that the lower value of tan β decreases the importance of the branching ratios of the heavy doublet states into down-type fermions in favour of cascade decays, via light Higgs or SUSY states, which hence affects the search strategy . The light CP-even doublet again shows SM-like couplings (within a few percent), except for the couplings to down-type fermions, which are slightly reduced: this results from the perturbation among H u and H d components, which is itself related to the triple-state mixing in the presence of a singlet component. As a consequence, LHC rates into gauge bosons are slightly enhanced (while the bb and τ + τ − channels are slightly suppressed): increased precision on the Higgs decays would thus prove interesting in such a configuration. The central test of this scenario however rests with the observability of the light CP-even singlet at ∼ 106 GeV: its reduced doublet component ≤ 3% indeed entails suppressed rates (due to the reduced production cross-section), at the percent level of those of a SM Higgs boson at the same mass. The corresponding cross-section at 8 TeV for the diphoton final state is at the level of 10 −3 fb. Direct searches thus carry little chance of success on short timescales. Again, one could then search the light state in Higgs pair production.
The third point is also an example of a light singlet Higgs, although it differs from the previous one in two respects: it involves very low tan β = 2, and the light singlet has a JHEP04(2016)095 significant doublet component ∼ 9%. The first feature, low tan β, plays a key role, together with λ ∼ 0.7, in increasing the mass of the 'observed' state in the presence of a light squark spectrum. We remind the reader that this effect is not bound to the presence of a light CP-even singlet (see section 5). The second aspect, S 2 13 ∼ 9%, has consequences on the -however still challenging -observability of the light singlet at ∼ 102 GeV, as apparent LHC rates into bb and τ + τ − would now reach a somewhat larger fraction (∼ 6-8%) of the corresponding values for a SM Higgs boson at this mass -remember from sections 4, 5 that the doublet components of the light singlet may take phenomenologically viable values as large as ∼ 20%. The corresponding signal would also fit the LEP excess in the bb channel adequately enough. Note, though, that the h 1 → γγ rate is not particularly large (∼ 25% of its SM value at this mass), and the corresponding cross-section at 8 TeV is under 1 fb for this particular example. On the other hand, as for the previous point, the rates of the SM-like state would give rise to slight excesses in the vector channels / a slight suppression in the down-type fermion channels.
Several (a priori separable) features are present for the fourth point. First, as for the previous point, the low tan β = 2 and large λ = 0.7 ensure that the mass of the light doublet state is generated in the correct range ∼ 125 GeV, despite a light sfermion spectrum and negligible trilinear couplings. Furthermore, the heavy doublet sector is essentially decoupled again. The most remarkable aspect of this point however rests with the fact that the light CP-even doublet and singlet states are almost mass-degenerate (within ∼ 0.1 GeV). Moreover, both quasi-degenerate states are mixed at almost ∼ 50%, which means that they carry 'half a SM-like Higgs boson' each (in their contribution to the LHC rates). Ideally, one should try and resolve experimentally the two states (their width is about 2 MeV). Another strategy consists in looking for deviations in the Higgs pair production rate, as singlet-doublet contributions would modify the apparent trilinear Higgs couplings at ∼ 125 GeV. The various couplings read (in units of g SM The apparent triple-Higgs coupling accounting for both states (in connection to external SM particles) gives ∼ 205 GeV, which represents a ∼ 7% increase with respect to the SM value (very difficult to resolve).
The fifth point is another example of the decoupling limit. Its aspect that we wish to discuss, however, is the presence of a CP-odd state at ∼ 44 GeV, that is below the h 1 → 2A 1 threshold. Yet, the corresponding width is strongly suppressed, in accordance with the LHC data, resulting in a branching ratio of only BR(h 1 → 2A 1 ) ∼ 10 −4 . The couplings of the CP-even state at ∼ 125.1 GeV and its rates at the LHC are thus essentially SM-like. On the other hand, the light CP-odd Higgs is almost a pure singlet so that it will evade detection in direct production. Production of A 1 from a heavier state is also problematic, as only the CP-even singlet h 2 has a sizable decay into this state. Therefore, the most serious hope for the observation of such a A 1 would still lie in increased precision in the characteristics of h 1 at future colliders. This point illustrates the extreme possibility of a light CP-odd state devoid of any phenomenological impact on the standard sector. As we stressed in section 6, however, it still makes sense to test a BR(h 1 → 2A 1 ) in the few percent range.

JHEP04(2016)095
In the case of the sixth point, the whole doublet sector lies below ∼ 125 GeV, with very light neutral states at ∼ 63 GeV while the charged Higgs has a mass of 101.5 GeV. Expectedly, the couplings of h 1 to SM gauge bosons are suppressed, which accounts for its invisibility at LEP. Its diphoton cross-section at 8 TeV is at the level of 1 fb, as that of the light A 1 state. The couplings of the heavy doublet state, at ∼ 125.6 GeV, are reasonably close to their SM value, resulting in SM-like cross-sections for the corresponding state (within a few percent): it thus offers an acceptable interpretation of the observed signals as the competitive χ 2 proves. Yet this point is discarded by Higgsbounds: indeed, BR(t → H + b) ∼ 1.12 · 10 −2 (together with BR(H − → τ ν τ ) 0.97) is beyond the recent limits established by CMS. Much tension is also present in B-physics observables where BR(B → τ ν τ ) and BR(B s → µ + µ − ) lie beyond 2 standard deviations.
While electroweak precision observables were not included within our fit, predictions for M W are displayed in table 3, using the tools presented in [118]; the values given in brackets (in MeV) are estimates of the theoretical uncertainty. Comparing these values with the current experimental measurement (M exp. W = 80.385 ± 0.015 GeV), we observe that all points remain within 1σ deviation of the experimental average, with possible tensions at low tan β or in the presence of light doublet states.

Attempt at a 'global' scan
Performing a global scan over the NMSSM parameter space proves to be both a challenging and potentially even misleading task. The NMSSM Higgs sector involves 6 degrees of freedom (λ, κ, tan β, M A , µ, A κ ) at lowest order, to which one should also add the stop masses and trilinear coupling (A t ), which have an important impact via loop corrections. The rest of the supersymmetric spectrum cannot be altogether left out from a scan if one keeps e.g. (g − 2) µ and the requirement for a neutralino LSP in the process. Moreover, some scenarii that we described in the previous sections and which provide a reasonable or even an excellent fit to LHC/TeVatron data involve specific parameter configurations (so that one obtains the appropriate spectrum/couplings), hence disappear or become uncompetitive if the scanning density is too loose. On the other hand, keeping a tight scanning density over so many variables would require both long scan durations and huge storage capacities.
Under these conditions, we settled for the following procedure. Since new NMSSM effects, with respect to the MSSM, essentially involve the tree-level parameters, we decided to freeze the supersymmetric input to a 'good' MSSM value (large A t and relatively massive squarks, along with light sleptons) and scanned over the 6 tree-level Higgs input quantities:  tan β is not spoilt, however (the MSSM tree-level mass would still be too small to provide a phenomenologically viable Higgs boson at low tan β).
We show the outcome of this scan in the {κ, λ} and the {M A , tan β} planes in figure 21: one observes that the best-fitting points (χ 2 < ∼ 77.5) lay in the λ ∼ κ ∼ 0.3-0.5 region and draw a band from (M A ∼ 500 GeV, tan β ∼ 5) to (M A ∼ 2000 GeV, tan β ∼ 15): they involve a light CP-even singlet. Next follow points fitting very well (χ 2 ∼ 77.5-79) in the decoupling limit with large tan β (M A > ∼ 1 TeV, tan β > ∼ 20-25). Points in the low tan β limit do not give competitive fits as they suffer from a large pull from (g − 2) µ (they still receive χ 2 ∼ 85, improving on the SM-limit). Figure 22 shows the singlet composition of the lightest Higgs state. In addition to the best-fitting scenarios that we just discussed -JHEP04(2016)095 → 125 GeV -one observes that light singlet states under ∼ 62 GeV may come with acceptable fit values, although their doublet components have to be strongly suppressed. Points involving a large singlet-doublet mixing are also found in the vicinity of m h 0 1 = 125 GeV. Furthermore, one observes a few points close to the x-axis for m h 1 < ∼ 100 GeV: they correspond to the light doublet scenario. While they result in a poor fit to the observed Higgs data χ 2 > ∼ 135, they escape exclusion from the unobserved top decay to a charged Higgs (with mass ∼ 110 GeV here) due to a suppressed branching ratio at large values of tan β ( > ∼ 20). Finally, note that the light CP-odd Higgs scenario is only represented by light CP-odd doublets (hence offers a poor fit to the Higgs data due to suppressed conventional channels for the state at ∼ 125 GeV) in our scan: light CP-odd singlets are barely represented as the scan density allows for very few points with low |A κ |. Figure 23 shows the properties of the heavier Higgs states. The plot on the top lefthand corner shows that the points which give the best fit tend to come with a light secondlightest Higgs state: this can be either the SM-like doublet state at ∼ 125 GeV or a slightly heavier singlet. We also find a good fit in a region with mass values of > ∼ 1 TeV: there, the second-lightest Higgs can be either singlet or heavy-doublet in nature. Interestingly, sizable mixings of the singlet with the heavy doublet state lead to fit values of secondary quality (orange vertical lines, where the 'line' shape is an artefact of the scan using discrete values JHEP04(2016)095  of M A ). In the top right-hand corner, the plot shows the composition of the heaviest CPeven state. Unsurprisingly, this state is the heavy doublet state for the best fit pointswhich involve light singlet states, as follows from our discussion above. However, a good fit is also obtained for the case where the singlet is very heavy, as indicated by the red points forming a horizontal line at S 2 23 = 1 (note that the apparent interruption in red points at 2.9 TeV is an artefact of the scan and that the corresponding 'red line' actually continues for heavier masses, not displayed in the figure). The plot in the lower row shows the χ 2 distribution on the plane defined by the masses of the two heaviest CP-even Higgs states. Two favoured regions appear: one, for m h 2 ∼ 125 GeV and m h 3 > ∼ 0.7 TeV, corresponds to the scenario with a light singlet state; on the contrary, in the other preferred region both m h 2 and m h 3 are heavier than 1 TeV, corresponding to the scenario where one doublet state is light while both the singlet state and the other doublet state are very heavy. Note that both these preferred regions share the common feature to involve a heavy doublet state with a mass close to or above 1 TeV: this is characteristic of the decoupling limit and gives rise to the fact that the light doublet state has SM-like properties.
Let us now investigate the outcome of this global scan in terms of the various scenarii that we considered. Representatives of all the types of spectra that we discussed are found:
• Point A is the best-fit point of the whole scan. It exhibits a light CP-even singlet at ∼ 101 GeV. As for the doublet sector, it is characteristic of the decoupling limit. That this best-fitting point involves a light CP-even singlet demonstrates clearly that this scenario is currently the most interesting and motivated one arising in the context of the NMSSM. The doublet component of the light state is of the order of 2%, so that the associated signals in direct production should be expected at the percent level of those of a SM Higgs boson.
• Point B involves a light CP-even doublet state below 100 GeV. The associated χ 2 is quite high, illustrating the limited pertinence of this interpretation of the LHC data.
• Points C and D are two representatives of the NMSSM at low tan β and large λ: in the first case, the CP-even singlet state is lighter than the doublet at ∼ 125 GeV, while it is heavier in the second case. Note that (g − 2) µ is a major pull in the present case, due to both low tan β and moderately light sleptons/neutralinos/charginos. In both cases, the mostly-singlet state has a doublet component at the percent level only, leading to suppressed production rates.
• Point E is a representative of the decoupling limit at 'large' tan β, although not in the MSSM limit. The large value of tan β ∼ 18 is essentially driven by (g − 2) µ . The heavy doublet states have large masses (∼ 2 TeV), while a CP-even singlet with a ∼ 2% doublet component appears at ∼ 200 GeV.
• Point F gives the best fit value for the MSSM limit: this is again a 'Decoupling Limit'-like point with large tan β, heavy H ± , H 0 , A 0 doublet states and even heavier singlets. The presence of relatively light sleptons and a bino at ∼ 250 GeV together with large tan β allow for a satisfactory fit to (g − 2) µ , while, for all other aspects, this point could be regarded as belonging to the SM limit.
• Two CP-even states, mixing singlet and light doublet almost at the level of 25%, intervene close to ∼ 125 GeV in Point G. While ATLAS and CMS may be able to resolve their mass gap of about 1.5 GeV, note that the degeneracy could be at the level of 100 MeV. While points with smaller mass gaps were present in the global scan, none showed a mixing between singlet and doublet components as large as Point G: the latter would allow for the observability of two separate peaks of comparable width in the spectrum.
• Point H involves, in addition to a CP-even doublet, one singlet-like CP-odd state close to ∼ 125 GeV. As the corresponding state is singlet at 99.8%, the effect in direct production is completely negligible. Effects in pair production could develop however, due to the h 1 − A 1 − A 1 coupling.
• For the Points I and J, light states are present below m h SM /2, hence potentially opening unconventional Higgs decays. For Point I, the light state is a CP-even singlet at ∼ 50 GeV: even though λ and κ are quite large, the unconventional Higgs decay is suppressed due to an accidental cancellation of the various terms entering the h 2 − h 1 − h 1 coupling: this avoids conflict with the LHC/TeVatron data. The doublet component of the light state is under 1%, explaining its compatibility with the LEP constraints. For Point J, a light doublet-like CP-odd state is present, leading to a large h 1 → 2A 1 branching ratio ∼ 0.7: the quality of the fit is correspondingly quite poor, despite the presence of a CP-even doublet at ∼ 124 GeV. Note that the case of light CP-odd Higgs states is not really probed in this scan since the low scan density does not allow for many points in the low |A κ | region (which provides the light CP-odd singlet states).

Conclusions
This discussion of the NMSSM Higgs scenarii in view of the Higgs-measurement data shows that, even though a conventional SM-like Higgs is expected at ∼ 125 GeV, many scenarios that differ very significantly from the SM case remain possible and even give competitive fits to the LHC/TeVatron data. In most cases, these configurations involve light additional degrees of freedom, CP-even and/or CP-odd, singlet or doublet-like, hiding invisibly at lower masses or combining their signals with those of the SM-like state. The most prominent NMSSM Higgs scenario which turned out to be favoured by the fits carried out in this paper involves a CP-even mostly singlet state with mass below that of the light doublet state (which is identified with the signal observed at LHC) and includes a small mixing of the singlet and doublet components. From the point of view of the model, this setup has the desirable feature of uplifting somewhat the mass of the doublet-like state, arguably rendering a mass of ∼ 125 GeV more natural. Small deviations, at the percent level, from the SM rates, e.g. slightly suppressed decays into down-type fermions or slightly enhanced vector channels, can also be generated for the state at ∼ 125 GeV via a disturbed proportion in H u and H d components, which is associated to the appearance of a small singlet component (three-state mixing). Such features are shared with many other configurations however (e.g. perturbation of the couplings via heavier particles contributing at the loop level), and the crucial test of this scenario would lie in the detection of the light CP-even singlet itself. In the most favourable cases, the latter may acquire doublet components of the order of 10% and could thus appear as a 'miniature' Higgs boson, with a production cross-section reduced to ∼ 10% of that of a SM state at the same mass. Its dominant decays would likely be bb and τ + τ − , although unconventional rates, e.g. enhanced γγ or even cc channels, could also appear in certain configurations of the doublet composition of the lightest state. As an additional feature, the 2.3 σ (local) excess noted at LEP in e + e − → Zh → bb could be accomodated with such a singlet carrying a sizable doublet component in the vicinity of ∼ 100 GeV. Direct searches, e.g. in the γγ or the τ + τ − channel, may thus eventually detect such a particle, provided that those searches are carried out in the low-mass region. On the other hand, reduced doublet components, in the percent range or below, are also compatible with the data and may equally well account for an uplift of the mass of the second lightest CP-even state. Discovery of the light state JHEP04(2016)095 could then prove difficult and, while effects in Higgs pair production are possible, improved searches would be required to ensure that the singlet does not remain unnoticed.
In contrast, the popular NMSSM scenario where the SM-like Higgs boson decays unconventionally and sizably into lighter Higgs states (CP-even or CP-odd) has been essentially emptied of its phenomenological content, due to the success of the Higgs searches in the standard channels. Yet, light states below ∼ 63 GeV are still an open possibility provided the decay of the 'observed' Higgs towards them is suppressed. From the point of view of the model, this amounts to a constraint on the Higgs sector. The main consequence for the hypothetical light states is that their coupling to the observed state should remain small, which implies, in particular, that they should be almost purely singlets. This further reduces the prospects for detecting the light Higgs state via channels alternative to the decays from the doublet states, as all couplings to SM particles (hence the cross-section in production channels involving fermions or vector bosons) are correspondingly suppressed.
A more exotic possibility is a scenario with additional quasi-degenerate states at ∼ 125 GeV. The NMSSM indeed allows for CP-even singlet states at this mass. While the offdiagonal singlet-doublet mass-entry should then remain small (so that the quasi-degeneracy is not lifted by this off-diagonal term), the mixing of such a singlet with the SM-like state could be as large as 50% (but also as small as 0%). Such a possible scenario motivates experimental efforts to improve the resolution in the mass measurement of the state at ∼ 125 GeV, as two separate peaks may eventually become distinguishable. Other possible, but not guaranted, effects would arise from modified apparent triple-Higgs couplings at 125 GeV. The situation is somewhat different if the additional state is CP-odd. While current limits then favour a mostly singlet state, a sizable doublet component may also be present (hence allow the state to be produced and contribute to the LHC rates). As the CP-odd Higgs does not couple to electroweak gauge bosons, its decay products, essentially fermionic, would add to those of the SM-like state and displace the apparent proportion among the various channels. Note that this latter case mimicks a scenario with CP-violation in the Higgs sector, where the observed state would have a CP-odd component.
A noteworthy aspect of the NMSSM, as compared to the MSSM, lies in its resilience at low values of tan β, where specific tree-level contributions associated with a large λ become important. Dependence on the supersymmetric spectrum to generate a Higgs mass close to ∼ 125 GeV is of secondary importance in this case, so that the hypothetical presence of light stops and small trilinear couplings can be compatible in the NMSSM with the discovered signal. In addition to the most commonly considered case where the SM-like state is the lightest CP-even Higgs, the previous scenarios, involving e.g. a lighter CP-even singlet or quasi-degenerate states, can be accomodated in this context as well.
Finally, the prospect of a light doublet Higgs sector, with light neutral states below 100 GeV, decoupling from electroweak gauge bosons, cannot be discarded altogether yet. While a significant adjustment of the parameters is necessary to avoid the relevant experimental limits, this scenario can serve as a motivation to cover remaining loopholes in the experimental searches. Further searches for light neutral but also charged Higgs states should soon provide more information about the validity of this scenario.
In summary, while the properties of the detected signal are so far compatible with the Higgs boson of the SM, the analyses of this paper clearly demonstrate that the current JHEP04(2016)095 results do not necessarily imply a minimalistic phenomenology but, on the contrary, call for a comprehensive investigation and precision tests of the Higgs properties.

JHEP04(2016)095
Open Access. This article is distributed under the terms of the Creative Commons Attribution License (CC-BY 4.0), which permits any use, distribution and reproduction in any medium, provided the original author(s) and source are credited.