Phenomenology of the inflation-inspired NMSSM at the electroweak scale

The concept of Higgs inflation can be elegantly incorporated in the Next-to-Minimal Supersymmetric Standard Model (NMSSM). A linear combination of the two Higgs-doublet fields plays the role of the inflaton which is non-minimally coupled to gravity. This non-minimal coupling appears in the low-energy effective superpotential and changes the phenomenology at the electroweak scale. While the field content of the inflation-inspired model is the same as in the NMSSM, there is another contribution to the μ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mu $$\end{document} term in addition to the vacuum expectation value of the singlet. We explore this extended parameter space and point out scenarios with phenomenological differences compared to the pure NMSSM. A special focus is set on the electroweak vacuum stability and the parameter dependence of the Higgs and neutralino sectors. We highlight regions which yield a SM-like 125GeV\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$125\,\hbox {GeV}$$\end{document} Higgs boson compatible with the experimental observations and are in accordance with the limits from searches for additional Higgs bosons. Finally, we study the impact of the non-minimal coupling to gravity on the Higgs mixing and in turn on the decays of the Higgs bosons in this model.


Introduction
In the history of our universe, there has been a period in which the size of the universe exponentially increased. This short period is known as inflationary epoch, and many models have been developed in order to explain the inflation of the early universe. Unfortunately, most of these models of inflation cannot be tested directly in the laboratory; the observation of the universe is the only discriminator to disfavor or support such models. Therefore, testing the phenomenology of a particle physics model of inflation at the electroweak scale with colliders is of interest both from the point of view of particle physics and cosmology.
One possibility to describe inflation is the extension of a particle physics model by additional scalar fields which drive inflation but are removed from the theory afterwards. A more economical approach is the idea of using the Higgs field of the Standard Model (SM) as inflaton [1][2][3]. The simplest version, however, is under tension as it suffers from a fine-tuning and becomes unnatural [4]. A less minimal version of Higgsportal inflation with an additional complex scalar field can in addition solve further problems of the SM, see Refs. [5,6]. Also the concept of critical Higgs inflation can raise the range of perturbativity to the Planck scale and solve further problems of the SM, see Refs. [7][8][9]. Other solutions are offered by scale-free extensions of the SM. A natural way of such an implementation can be realized in canonical superconformal supergravity (CSS) models as proposed by Refs. [10,11] based on earlier work by Ref. [12].
The Higgs inflation in the supergravity framework is triggered by a non-minimal coupling to Einstein gravity. For the supergravity Lagrangian this can be achieved with an additional term X (ˆ ) R of chiral superfieldsˆ and the curvature multiplet R (the supersymmetrized field version of the Ricci scalar which contains the scalar curvature in the Grassmannian coordinate θ 2 ), following the notation of Ref. [12]. The Lagrangian then reads (1) where X (ˆ ) as well as the Superpotential W(ˆ ) are holomorphic functions of the (left) chiral superfieldsˆ , E is the vierbein multiplet andD a covariant derivative. The ellipses encode further gauge terms. The only possible choice of such a non-minimal coupling suitable for inflation is given by [12] where χ is a dimensionless coupling andĤ d,u contain the two SU (2) L Higgs doublets of the Next-to-Minimal Supersymmetric Standard Model (NMSSM). 1 The extension by an additional scalar singlet like in the NMSSM has been shown to be a viable model for inflation, although this version suffers from a tachyonic instability [13]. In order to avoid this instability, a stabilizer term has been introduced in Refs. [11,13] that is suppressed at low energies. The stabilizer term can be avoided in a model with minimal supergravity couplings where the Kähler potential has a shift symmetry in the doublet fields [14]; however, cosmological phenomenology and observations have meanwhile ruled out this possibility [15].
The simplest implementation of a superconformal model which can accommodate the non-minimal coupling term χĤ u ·Ĥ d is the well-known Z 3 -invariant NMSSM augmented by an additional μ term, which we call μ-extended NMSSM (μNMSSM) in the following. We neglect all additional Z 3 -violating parameters in the superpotential at the tree level (see the discussion below). These terms are not relevant for the physics of inflation: the function X could potentially also contain anŜ 2 term, since it has the same structure asĤ u ·Ĥ d and is allowed by gauge symmetries. However, inflation driven by this term does not lead to the desired properties as pointed out in Ref. [12]. The other term, which is not present in the NMSSM, is a singlet tadpole proportional toŜ that is not quadratic or bilinear in the chiral superfields and thus would need a dimensionful coupling to supergravity instead of the dimensionless χ .
In this work, we are going to study the low-energy electroweak phenomenology of the model outlined in Refs. [10,11] and Ref. [13], where previously the focus was put on the description of inflation and the superconformal embedding of the NMSSM into supergravity. We have generated a model file for FeynArts [16,17], where SARAH [18][19][20][21] has been used to generate the tree-level couplings of the μNMSSM, and we have implemented the one-loop counterterms. The loop calculations have been carried out with the help of FormCalc [22] and LoopTools [22]. In order to predict the Higgs-boson masses, we have performed a one-loop renormalization of the Higgs sector of the μNMSSM which is compatible with the renormalization schemes that have been employed in Refs. [23,24] for the cases of the MSSMand NMSSM, respectively. This allowed us to add the leading MSSM-like two-loop corrections which are implemented in FeynHiggs [25][26][27][28][29][30][31][32] in order to achieve a state-of-the-art prediction for the Higgs masses and mixing. The parameter space is checked for compatibility with the experimental searches for additional Higgs bosons using HiggsBounds version 5.1.0beta [33][34][35][36][37] and with the experimental observation of the SM-like Higgs boson via HiggsSignals version 2.1.0beta [38]. In addition, we check the electroweak vacuum for its stability under quantum tunneling to a non-standard global minimum and for tachyonic Higgs states in the tree-level spectrum. Finally, we investigate some typical scenarios and study their collider phenomenology at the Large Hadron Collider (LHC) and a future electron-positron collider. For this purpose in some analyses we use SusHi [39,40] for the calculation of neutral Higgs-boson production cross-sections. We emphasize the possibility of light CP-even singlets in the spectrum with masses below 100 GeV that could be of interest in view of slight excesses observed in the existing data of the Large Electron-Positron collider (LEP) [41] and the Compact Muon Solenoid (CMS) [42] which are compatible with bounds from A Toroidal LHC Appara-tuS (ATLAS) [43]. For one scenario that differs substantially from the usual NMSSM, we exemplarily discuss the total decay widths and branching ratios of the three lightest Higgs bosons and their dependence on the additional parameters of the μNMSSM.
The paper is organized as follows: we start with a description of our model and the theoretical framework in Sect. 2 by discussing analytically the phenomenological differences of the Higgs potential in the μNMSSM compared to the Z 3invariant NMSSM. We study vacuum stability and the incorporation of higher-order corrections for the Higgs boson masses. Then, we derive the trilinear self-couplings of the Higgs bosons and comment on the remaining sectors of the model which are affected by the additional μ term. In Sect. 3, we focus on the parameter space of interest and investigate the Higgs-boson masses as well as the stability of the electroweak vacuum numerically and also show the neutralino spectrum. Furthermore, we study the effect of the additional μ parameter on Higgs-boson production and decays. Lastly, we conclude in Sect. 4. In the Appendix we present the beta functions for the superpotential and some soft-breaking parameters of the general NMSSM (GNMSSM) [44][45][46] including all Z 3breaking terms.

Theoretical framework
In this section we introduce the model under consideration, the μNMSSM, which differs by an additional μ term from the scale-invariant NMSSM. We derive the Higgs potential and investigate vacuum stability and the prediction for the Higgs-boson masses of the model. Furthermore, we discuss the trilinear self-couplings of the Higgs bosons and comment on the electroweakinos -i.e. charginos and neutralinos -as well as on the sfermion sector. We constrain our analytical investigations in this section mostly to tree-level relations. Higher-order contributions, e.g. for the Higgs-boson masses, are explained generically and are evaluated numerically in the subsequent phenomenological section.

Model description
For the Higgs sector of the NMSSM the superpotential is of the form 2 W Higgs = λŜĤ u ·Ĥ d + 1 3 κŜ 3 .
whereĤ u andĤ d are the well-known SU (2) L doublets of the MSSM, andŜ is the additional SU (2) L singlet. The SU (2) L -invariant productĤ u ·Ĥ d is defined throughĤ u · H d = a,b abĤ a dĤ b u with 21 = 1, 12 = −1 and aa = 0 with a, b ∈ {1, 2}. As outlined in Ref. [11], a Kähler transformation starting from Jordan-frame supergravity introduces a correction in the superpotential, which is of the form The parameter m 3/2 denotes the gravitino mass, and χ is the coupling of Eq. (2). The scalar Higgs fields are denoted by H u , H d and S in the following. During electroweak symmetry breaking, they receive the vacuum expectation values (vevs) v u , v d and v s , respectively. Expanding around the vevs, we decompose the fields as follows: The additional bilinear contribution to the superpotential in Eq. (4) generates a term which is analogous to the μ term of the MSSM, but with μ = 3 2 m 3/2 χ.
When the singlet S acquires its vev, an effective μ eff = λ v s is dynamically generated. Often, the sum (μ + μ eff ) is the phenomenologically more relevant parameter of the model. It takes the form μ + μ eff = 3 2 m 3/2 χ + λ v s (7) and corresponds to the MSSM-like higgsino mass term. In the following, we consider both quantities μ and μ eff as independent input parameters, where μ is linearly dependent on the gravitino mass m 3/2 . In order to be a viable dark-matter candidate, the gravitino mass can range from a few eV to multiple TeV, see e.g. Ref. [47]. The value of χ is a priori not fixed; for cosmological reasons we adopt according to Refs. [11,13]. The additional contribution to the superpotential in the μNMSSM is thus mainly steered by the gravitino mass, whereas v s can be traded for μ eff . If we require a μ parameter above the electroweak scale, μ 1 TeV, and in addition a sizable coupling λ 0.1, the typical gravitino mass turns out to be much below the electroweak scale at m 3/2 10 MeV. However, if we allow for very small values of λ 10 −2 and very large values of μ 1 TeV, the gravitino mass could as well be above the TeV scale. In the latter case, the phenomenology of the μNMSSM is not necessarily similar to the MSSM: the singlets only decouple for λ → 0 with κ ∝ λ and therefore v s → ∞. If the constraint κ ∝ λ is dropped, interesting effects can occur; e.g. we will discuss a scenario with small λ and small μ eff in our numerical studies. In contrast to the NMSSM, the higgsino mass can be generated by μ alone and thus even a vanishing v s is not in conflict with experimental bounds.
In order to avoid the cosmological gravitino problem [48], where the light gravitino dark matter overcloses the universe [49,50], one has to control the reheating temperature in order to keep the production rate of the light gravitinos low [51]. This potential problem may affect the model under consideration for gravitino masses in the range from MeV to GeV; it disappears for much heavier gravitinos ( 10 TeV). In the latter case the inflationary μ term would dominate over the NMSSM-like μ eff and drive the higgsino masses to very high values (unless μ eff is tuned such that the sum (μ + μ eff ) remains small). For gravitino masses m 3/2 > 1 GeV it affects Big Bang Nucleosynthesis via photo-deconstruction of light elements, see Ref. [48]. As discussed in Ref. [11], in the μNMSSM there is no strict constraint on the reheating temperature T R . We note that a reheating temperature below T R 10 8 -10 9 GeV, as advocated in Ref. [52], avoids the gravitino problem. The rough estimate of m 3/2 ∼ 10 MeV even needs T R 10 5 GeV in order to not overclose the universe with thermally produced gravitinos after inflation [53][54][55][56]. Interestingly, such low reheating temperatures preserve high-scale global minima after inflation, see Ref. [57], and disfavor the preparation of the universe in a meta-stable state after the end of inflation [58]. In any case, the reheating temperature at the end of inflation is very model dependent and rather concerns the inflationary physics. A study to estimate the reheating temperature T R is given in Ref. [59]. Therein, a relation is drawn between the decay width of the inflaton and T R . Interestingly, if we naïvely assume that this width at the end of inflation is equal to the SM-like Higgs width h ≈ 4×10 −3 GeV, we can estimate a rather low reheating temperature T R ∼ √ h M Pl ≈ 10 7 GeV with the Planck mass M Pl ≈ 2.4 × 10 18 GeV. For our studies below we assume that a reheating temperature as low as T R 10 9 GeV can be achieved even with large couplings.
Since the bilinear μ term breaks the Z 3 symmetry, additional parameters are allowed compared to the NMSSM. In the general NMSSM (GNMSSM) -including the bilinear singlet mass parameter ν and the singlet tadpole coefficient ξ -the Higgs sector of the superpotential is given by However, we assume that the non-minimal coupling of the Higgs doublets to supergravity is the only source of superconformal and thus Z 3 symmetry breaking -as outlined in Section 5 of Ref. [11]. In this case, all other superpotential parameters that are forbidden by Z 3 symmetry remain exactly zero at all scales: the beta functions for the parameters of the superpotential are proportional to the respective parameter itself and thus they cannot be generated radiatively. Because the Z 3 symmetry is broken (which avoids the typical domain-wall problem of the NMSSM [60]), another symmetry at the high scale is required in order to solve the tadpole problem [61][62][63][64][65][66]: without such a symmetry, Planckscale corrections could possibly induce large contributions to the tadpole term [67]. The superconformal embedding of the μNMSSM, where the μ term is generated from the Kähler potential, serves as this symmetry. As pointed out in Ref. [67], other possibilities consist of discrete or continuous nongauge symmetries, so-called R symmetries. Imposing discrete Z 4 or Z 8 R symmetries as proposed in Refs. [45,68,69] provide a viable solution, since dimensionful linear and bilinear terms are forbidden as long as the symmetry is not broken. 3 Furthermore, each parameter in the superpotential induces a corresponding soft-breaking term; additional mass terms are allowed: It should be noted that the beta functions for soft-breaking parameters are not only proportional to themselves, but also receive contributions from the other soft-breaking parameters. Thus, in contrast to the terms in the superpotential, finite contributions may emerge even if a soft-breaking parameter is set to zero at the tree level. The beta functions for the parameters of the superpotential in Eq. (9) and its corresponding soft-breaking parameters in Eq. (10) can be found in Refs. [44,71,72]; however, since we employ different conventions we list them in Appendix A. Contrary to studies in the GNMSSM (see Refs. [44][45][46]73]), where the MSSM-like μ term can be easily shifted away and absorbed in a redefinition of the other parameters -especially the tadpole contribution -we cannot do so in the inflation-inspired μNMSSM. First of all, the μ term is introduced via the R symmetry-breaking non-minimal coupling to supergravity only. The other parameters in the singlet sector are not supposed to be generated by this breaking. Secondly, by redefining the parameters, we would introduce a tadpole term and shift the effect simply there. Note that the authors of Ref. [45] perform this shift in order to eliminate the linear (i.e. tadpole) term in the superpotential and keep μ, while others (e.g. Ref. [74]) shift the μ term to zero and keep the tadpole and bilinear terms for the singlet in the superpotential. As discussed above, in the μNMSSM considered in this paper due to the superconformal symmetry breaking at the Planck scale solely the Z 3 -breaking μ term is present.

Higgs potential
With the superpotential of Eq. (9) and the soft-breaking Lagrangian of Eq. (10), we derive the following Higgs potential, where we stick to real parameters: This potential can be expanded in the components of the Higgs fields in Eq. (5). Defining the vectors in field where the CP-even and CP-odd tadpole coefficients T S and T P , the CP-even, CP-odd and charged squared mass matrices M 2 S , M 2 P and M 2 C are given below, and the trilinear couplings λ i jk andλ i jk are specified in Sect. 2.5, though in a basis where the Goldstone mode corresponds to a mass eigenstate and does not mix with the other states at lowest order. The ellipses denote quadrilinear terms which are immaterial for the following.
We substitute the electroweak vevs v u and v d by their ratio tan β = v u /v d and the sum of their squares v 2 The symbols t β , c β and s β denote tan β, cos β and sin β, respectively. Furthermore, g 1 and g 2 are substituted by the W and Z gauge-boson masses, Using the abbreviations we can write the explicit expressions for the tadpole coefficients T S,P as The minimization of the Higgs potential requires all tadpole coefficients in Eq. (15a) to be equal to zero. With the conditions T S = 0 we choose to eliminate m 2 H d , m 2 H u and m 2 S according to Substituting these expressions in the symmetric mass matrices M S,P,C we find Diagonalizing Eq. (17c) yields zero for the massless charged Goldstone boson, and the charged Higgs-boson mass m H ± at the tree level is given by which we employ as an input parameter. Inserting Eq. (14a) we can then eliminate A λ via Substituting A λ in the abbreviations of Eq. (14) yields (a 2 , a 5 and a 7 are not changed) The tree-level masses of the three neutral CP-even Higgs bosons m 2 h 1,2,3 are determined by diagonalizing Eq. (17a). Analogously, diagonalizing Eq. (17b) yields the masses m 2 a 1,2 of the CP-odd Higgs bosons at the tree level; the third eigenvalue is equal to zero and belongs to the neutral Goldstone boson.

Higgs doublets
The mass-matrix elements of the doublet fields in the upper-left (2 × 2) block matrices of Eqs. (17a)-(17b) contain the abbreviation a 1 . From Eq. (20a) it is apparent that they are determined by SM parameters and m H ± , λ and t β like in the NMSSM.
Neglecting the mixing between the doublet and singlet sector, the mass of the light CP-even doublet state has an upper bound of m 2 Z c 2 2β + λ 2 v 2 s 2 2β . In the limit m H ± m Z , the other two doublet fields decouple and obtain a mass close to m H ± . Smaller values of m H ± increase the mixing of both CP-even doublet fields. Also t β needs to be close to one for large doublet mixing. Higgs singlets The (3, 3) elements of M S and M P in Eqs. (17a) and (17b) set the mass scale of the Higgs singlets. They contain the terms a 4 from Eq. (20c), a 5 from Eq. (14e), and a 7 from Eq. (14g). All Z 3 -violating parameters besides μ and B μ appear in these terms; in our later analysis we set these parameters besides μ and B μ to zero, but for completeness we mention them in the following discussion of this section. The parameter A κ appears only in the term a 5 , whereas B ν only appears in a 7 . Thus it is obvious that the diagonal mass-matrix elements for the singlet fields -and therefore their masses -can be controlled by these two quantities, without changing any other matrix element. If all Z 3violating parameters except μ and B μ were set to zero, we would rediscover the NMSSM-specific feature that A κ is bound from below and above to avoid tachyonic singlet states at the tree level. The ratio κ/λ which appears in both terms, a 5 and a 7 , has sizable impact on the mass scale of the singlets. If κ λ the CP-even singlet entry is purely controlled by a 4 , which in turn is proportional to 1/μ eff ; in the same limit, the CP-odd singlet entry is controlled by a 4 and the remainder of a 7 which is B ν ν. Also note that a 4 contains a term which is linear in μ. In the opposite case κ λ, the term a 5 is likely to dominate the (3, 3) matrix element for the CP-even singlet due to the suppression of a 4 by μ eff if it is of the order of a few 100 GeV. The term a 5 is proportional to (κ/λ) 2 μ 2 eff , such that the CP-even singlet exhibits a strong dependence on μ eff . On the other hand for μ μ eff , the term a 4 can balance the large κenhanced contribution in a 5 ; thus, possible upper bounds on κ as derived in Ref. [75] might be evaded. For the case of the CP-odd singlet, the terms in a 5 and a 7 that are quadratic in μ eff cancel each other. Then the size of the other parameters (especially A κ , μ and μ eff ) determines which contribution is dominant. For moderate values of κ ≈ λ 0.1 together with small A κ the CP-odd singlet develops a dependence on μ/μ eff , as we will discuss later. Lastly, we note that in the case of κ λ and A κ = 0 GeV the CP-even and CPodd singlet masses are controlled through (κ/λ) 2 μ 2 eff and (κ/λ) μ eff A κ , respectively. Later, this will allow us to present a rescaling procedure that keeps both singlet masses constant over a large parameter range. Doublet-singlet mixing The masses of the doublet-like and the singlet-like Higgs states can be significantly shifted by mixing between both sectors. The relevant matrix elements are the ones in the third columns of Eqs. (17a) and (17b). They contain the abbreviations a 2 , a 3 and a 6 , see Eqs. (14b), (20b) and (20d), respectively. The mixing vanishes in the limit λ → 0 with constant κ/λ, and it is enhanced for larger values of λ. For fixed λ it is also strongly enhanced in the limit μ eff → 0 GeV. In the CP-even sector, two terms contribute to the doublet-singlet mixing: a 2 which depends on the sum (μ + μ eff ), and a 3 which does not directly depend on μ, but only on the soft-breaking term B μ μ. In the case of large μ and μ eff of the same sign, a 2 often dominates the mixing with the lighter doublet, eventually yielding a tachyonic singlet or doublet Higgs; this behavior can be avoided by choosing a proper value for B μ (or ξ ) to cancel the large effect in a 2 by a 3 . In the case of similar μ and μ eff of opposite signs, a 3 will always dominate the mixing. Again, the mixing strength can be adjusted by setting B μ (or ξ ). The doublet-singlet mixing in the CP-odd sector contains only one term a 6 which is similar to a 3 with opposite sign. Furthermore, the CP-odd mixing elements can be modified by non-zero ξ and ν. As indicated above, due to the dependences of a 3 and a 6 on 1/μ eff , a small μ eff 100 GeV yields a strong mixing between singlets and doublets.
We subsequently discuss vacuum structure and vacuum stability bounds in the μNMSSM around the electroweak scale. We do not discuss tachyonic instabilities during inflation or the stabilization of the inflationary direction, since they are not of relevance for our study (see e.g. Refs. [11,13]).

Vacuum structure and vacuum stability bounds
The space of model parameters can be constrained using experimental exclusion limits and theoretical bounds. Those constraints can be applied to rule out certain parts of the parameter space. In this context, constraints from the stability of the electroweak vacuum appear to be very robust and theoretically well motivated. It has already been noticed in the early times of supersymmetry that constraints from the electroweak vacuum stability on the trilinear soft SUSY-breaking parameters can be important [76][77][78][79][80][81][82][83][84]. Recently they have been rediscussed in light of the Higgs discovery [85][86][87][88][89]. These constraints are usually associated with non-vanishing vacuum expectation values of sfermion fields (e.g. staus or stops) and thus known under the phrase "charge-and colorbreaking minima". Such minima can invalidate the electroweak vacuum and therefore lead to unphysical parameter configurations (see below). However, the existence of charge-and color-breaking minima is only a necessary condition for the destabilization of the electroweak vacuum. Clearly one has to compare the value of the potential at this new minimum with the desired electroweak one, and only if the non-standard vacuum is deeper the corresponding scenario is potentially excluded. In fact, some of the points with a deeper non-standard vacuum may be valid when accepting meta-stable vacua under the condition that the transition time from the local electroweak vacuum to the global true vacuum appears to be longer than the age of the universe [90]. However, the possibility of the existence of meta-stable vacua is of limited practical relevance for our analysis: typically only parameter points in close neighborhood to the stable region are affected by such considerations; well-beyond the boundary region, the false vacua become rather short-lived and thus are strictly excluded. In addition, there are thermal corrections in the early universe which give a sizable and positive contribution to the effective potential as the one-loop corrections are proportional to m 2 (φ) T 2 for the field-dependent masses m(φ). For finite temperature, they shift the ground state to the symmetric phase around φ = 0 GeV [91,92]. We presume, however, that our inflationary scenario preselects a vacuum at field values different from zero and, thanks to the relatively low reheating temperatures in our scenario, gets caught in it, see Ref. [57]. Following the inflationary scenario of Ref. [11], the trajectory in field space lies at β = π/4 with h 2 u = h 2 d = h 2 and s = 0 GeV; the presence of the singlet field S is needed for the stabilization of the inflationary trajectory in order to not fall into the tachyonic direction as pointed out by Refs. [11,13]. Inflation ends at field values h = O(0.01) in units of the Planck mass. For small λ ∼ 10 −2 , the D-flat trajectory remains stable after inflation ends according to Ref. [11], and will change to β = π/4 and s = 0 GeV when the SUSY-breaking terms become important. NMSSM-specific effects like the relevance of singlet Higgs bosons and the additional contribution to the 125 GeV Higgs boson are usually connected to a large value of λ. This is not necessarily the case in the μNMSSM, where striking differences also appear for small values of μ eff . Moreover, we will take it as a working assumption that after inflation ends, even for larger values of λ the universe will remain in the state with the inflationary field direction until it settles down in a minimum closest to this direction. If it is the global minimum of the zero-temperature potential, reheating may not be sufficient to overcome the barrier and to select a false (and maybe metastable) vacuum. The thermal history of the universe plays then no role for the choice of the vacuum, and in this case the universe would remain in the global minimum. Accordingly, we adopt the prescription to exclude all points with a global minimum that does not coincide with the electroweak vacuum. This means that we do not consider meta-stable electroweak vacua as they are excluded by the selection rule. A similar discussion and argument has been given in Ref. [93], where a selection of the vacuum with the largest expectation values was promoted, irrespective whether or not it is the global minimum of the theory.
We will see that actually in most cases scenarios are excluded because of a tachyonic Higgs mass. Tachyonic masses are related to the fact that the electroweak pointaround which the potential is expanded -is not a local minimum in the scalar potential, but rather resembles a saddle point or even local maximum, and the true vacuum lies at a deeper point along this tachyonic direction. Thus, the true vacuum has vevs different from the input values, and the electroweak breaking condition T S = 0 in Eq. (15a) does not select a minimum.
We briefly sketch how to get constraints on the relevant model parameters in the (neutral) Higgs sector of the μNMSSM. Similar observations for the NMSSM have been intensively discussed in the literature [94,95]. Already the presence of an additional Higgs singlet (see e.g. Refs. [96][97][98]) invalidates the well-known results that no chargebreaking Higgs vevs exist at lowest order in the MSSM (see e.g. Refs. [82,99]) and in two-Higgs-doublet models (see e.g. Refs. [100,101]). On the other hand, in the NMSSM the inclusion of such charge-breaking minima has rather little impact on the overall vacuum stability and gives no further information, see Ref. [102]. In a similar manner, we neglect non-vanishing squark vevs (see discussion below) and therefore we only have to deal with the following potential: where we just presented the real fields as we do not consider spontaneous CP violation. 4 Notice also that we do not consider the shifted theory with all fields In our case for the stability analysis, the potential vanishes at the origin, and the electroweak minimum is one of the minima not located at the origin. It is not necessarily the global minimum. Furthermore, compared to Eq. (11), we neglect all additional Z 3 -breaking terms besides the contri-butions of μ and B μ μ of the μNMSSM (see the discussion above). The "desired" electroweak vacuum can be constructed by fulfilling the minimization conditions at the tree level, T S = 0, with T S given by Eq. (15a). The vevs of the doublet fields are taken as fixed input parameters, whereas the value of μ eff is treated as variable similar to μ. These equations can be solved for the soft-breaking masses m 2 H u , m 2 H d and m 2 S according to Eq. (16). The masses of the Higgs sector are determined in such a way that the desired vacuum with h u = v u , h d = v d and s = μ eff /λ is a viable vacuum of the potential V in Eq. (21). However, one has to ensure that there is no deeper minimum of V . This can only be achieved reasonably-well through a numerical evaluation. For that purpose, we determine the stationary points of the potential V and then compare the corresponding values of V at these points with the desired minimum given by From the expression in Eq. (22), one can derive a few general results: (a) for small values of λ the desired minimum gets deeper and -as the singlet contribution decouples from the rest of the potential -it becomes more difficult for a non-standard vacuum to appear and to be deeper than the desired minimum; (b) the (μ)NMSSM potential at the desired minimum is usually deeper than in the case of the MSSM 5 and is mainly driven by μ eff ; (c) the contribution of A λ plays a subdominant role compared to A κ whose impact is strongly influenced by μ eff and λ; (d) parameter points with V des min > 0 have to be excluded because the trivial minimum at h u = h d = s = 0 GeV is obviously deeper.
In our analysis, we focus for clarity on constraints from the tree-level potential, considering the appearance of global non-standard minima and, as discussed above, disregarding the possibility of meta-stable false vacua. Employing higher-order (i.e. one-loop) corrections does not necessarily give more accurate predictions of vacuum stability, see Ref. [103]. An approach to include one-loop effects using a certain numerical procedure has been implemented in the public code collection of Vevacious, see Ref. [104], including a tunneling calculation also at finite temperature using CosmoTransitions [105]. The tree-level evaluation is much faster and numerically more stable; moreover, it has been argued that the one-loop effective potential is problematic for tunneling rate calculations [106].
Constraints on the NMSSM parameters There are two main constraints known for the trilinear soft SUSY-breaking parameters A κ and A λ . The first constraint relies on the existence of a non-vanishing singlet vev to generate μ eff = 0 GeV. This can be easily derived from the Higgs potential with only s = 0 GeV and is given by the requirement [75] This lower bound on A κ is inappropriate for the μNMSSM, as there always exists a non-vanishing higgsino mass term from μ = 3 2 m 3/2 χ . As shown in Sect. 3, this constraint has hardly any impact on our analyses. We simply keep it for illustrative reasons.
The second constraint, on A λ , follows from a nontachyonic charged Higgs mass, since a tachyonic mass (m 2 < 0 GeV 2 ) means that the potential has negative curvature at this stationary point derived by the minimization conditions. Thus, the true vacuum would have some non-zero vev for a charged Higgs component. Configurations like this are possible in the NMSSM, whereas they do not exist as global or local minima in the MSSM [82]. From the (tree-level) charged Higgs mass in Eq. (18), we get an indirect bound on A λ . Taking m H ± as input value, we can eliminate A λ as free parameter, see Eq. (19). Hence, we can ensure that m 2 H ± is always positive. Still, it is worth noticing that by this procedure A λ gets strongly enhanced for small μ eff (compared to m H ± ) and thus drives tachyonic neutral Higgs bosons.
Charge and color breaking There exist quite strong constraints in the MSSM from the formation of non-standard minima which break the electric and color charges, known as charge-and color-breaking (CCB) minima. The famous "A-parameter bounds" read traditionally [76,80,82,107] where m 2Q and m 2 t,b are the soft SUSY-breaking masses for the superpartners of the left-handed SU (2) L quark doublet,Q, and of the right-handed quark singlets,t andb. Several modifications and improvements of Eq. (24) are present in the literature, see e.g. Refs. [82,84,90]. These constraints follow from the "D-flat" directions in the scalar potential of the MSSM, i.e. h u =t L =t R and h d =b L =b R , respec-tively. Thus the quartic terms associated with squared gauge couplings vanish. In addition, one has to be reminded that Eq. (24) are only necessary conditions for the formation of a non-trivial minimum with non-vanishing squark vevs in that specific direction. In the case of a violation of Eq. (24), one has to check that the generated CCB vacuum is actually deeper than the electroweak minimum. In the MSSM the desired minimum takes on a comparably small numerical value, only depending on c 2β (and the B μ term which can be replaced by the CP-odd Higgs mass M A ): In principle, the A-parameter bounds (24) can be simply transferred to the μNMSSM, where μ has to be replaced by (μ+μ eff ), as they can be transferred to the NMSSM [108]. The net effect is roughly the same in the MSSM, NMSSM and μNMSSM; if A t fulfills Eq. (24a), no CCB will appear. Constraints on μ eff alone may get weakened, because the desired minimum also gets deeper for larger μ eff . Moreover, the additional singlet direction stabilizes the potential with respect to CCB minima since the μ eff term originates from a quadrilinear scalar coupling, and the vacuum with nonvanishing μ eff or v s is typically deeper than a CCB vacuum. Generically, constraints from the coupling to the wrong Higgs doublet relating down-type sfermion vevs to the uptype Higgs and vice versa, see Refs. [109,110], are expected to be valid for (μ + μ eff ) and not weakened if the singlet is fixed at its vev. Similarly, there are bounds on A t,b not related to D-flat directions as discussed in Ref. [111]. These can be reasonably-well determined only numerically. Generically speaking, for the μNMSSM the risk of generating a CCB vacuum is reduced because (a) the dependence of the desired minimum on μ eff drives the electroweak vevs to be more stable, and (b) not as large values of A t are needed to raise the SM-like Higgs mass because of the additional NMSSM-specific tree-level contribution. Constraints from CCB minima as given in Eq. (24), are less important in comparison to the MSSM for both, the NMSSM and the μNMSSM, even if large stop corrections are needed to shift the SM-like Higgs mass (as in the case for small λ). If the singlet-field direction were neglected and the stop Dflat directiont R =t L =t defined, one could directly apply Eq. (24) for the μNMSSM, keeping v s = 0 GeV and replacing μ → μ + μ eff . However, with the singlet as dynamical degree of freedom, the stability of the electroweak vacuum is improved as the only singlet-stop contribution is actually a quadrilinear term λ h d st 2 and the occurrence of a true vacuum with h u,d = v u,d , s = v s and t = 0 GeV is disfavored.
Meta-stability and tunneling rates Lastly, we comment on vacuum-to-vacuum transitions in case of a local electroweak vacuum. It is in general of interest to see how long such a meta-stable state could survive compared with the life-time of the universe. We have outlined some arguments why -in view of the inflationary history of the universe -we disregard meta-stable long-lived vacua. We will see in Sect. 3.3 that totally stable points survive in a wide range of the parameter space.
For an estimate of the bounce action of the unstable configuration [112], we define an effectively single-field scalar potential linearly interpolating between the electroweak local minimum and the true vacuum found by the numerical minimization of the scalar potential at different field values and apply an exact solution of the quartic potential given by Ref. [113]. See also Ref. [114] for the application of this method to the μNMSSM.

Higher-order corrections to Higgs-boson masses and mixing
It is well-known that perturbative corrections beyond the tree level alter the Higgs masses and mixing significantly in supersymmetric models. For instance, in the MSSM such large corrections are needed to lift the lightest CP-even Higgs mass beyond the Z -boson mass. On the other hand, in the NMSSM and similarly the μNMSSM there are scenarios where an additional tree-level term lowers the tension between the tree-level SM-like Higgs mass and the measured value of the SM-like Higgs boson at 125 GeV. Still, since loop corrections to the Higgs spectrum have a large impact, in our phenomenological analysis we take into account contributions of higher order as described in the following. The masses of the Higgs bosons are obtained from the complex poles of the full propagator matrix. The inverse propagator matrix is a (6 × 6) matrix that readŝ −1 Hereˆ S andˆ P denote the matrices of the renormalized self-energy corrections to the neutral CP-even and CP-odd Higgs fields. In the CP-conserving limit there are no transition elements between CP-even and CP-odd degrees of freedom, which is why Eq. (26) is block diagonal.
In principle, contributions from mixing with the longitudinal Z boson have to be considered as well. However, these contributions as well as those from mixing with the Goldstone mode enter the mass predictions only at sublead-ing two-loop level [115,116]. Since these contributions are numerically small [117] we neglect them in the following and use a (5 × 5) propagator matrix. The (5 × 5) matrices are denoted by the symbolsˆ hh for the propagators andˆ hh for the renormalized self-energies in the following. The complex poles of the propagator are given by the values of the squared external momentum k 2 for which the determinant The real part, M 2 h i , of each pole yields the loop-corrected mass of the corresponding Higgs boson h i .
In this work, a model file for FeynArts [16,17] of the GNMSSM at the tree level has been generated with the help of SARAH [18][19][20][21]. In addition, the one-loop counterterms for all vertices and propagators have been implemented, and a renormalization scheme which is consistent with Refs. [23,24] for the cases of the MSSM and NMSSM has been set up. All Z 3 -violating parameters are renormalized in the DR scheme, see Appendix A for a list of the respective beta functions. The numerical input values of all DR-renormalized parameters are understood to be given at a renormalization scale which equals the top-quark pole mass. The renormalized self-energies of the Higgs bosonsˆ hh are evaluated with the help of FormCalc [22] and LoopTools [22] by taking into account the full contributions from the GNMSSM at the one-loop order. For other variations of the NMSSM, similar calculations of Higgsmass contributions up to the two-loop order have been performed in Refs. [118][119][120][121][122][123][124][125][126]. A comparison of results from public codes using different renormalization schemes can be found in Refs. [127,128].
As an approximation, we have added the leading two-loop contributions in the MSSM of O(α t α s ) [129] and O α 2 t [130,131] at vanishing external momentum to their MSSM-like counterparts in the μNMSSM (for a discussion of this approximation in the NMSSM see Ref. [124]). They are taken from their current implementation in FeynHiggs [25][26][27][28][29][30][31][32]. 6 We thus havê We note that the two-loop contributions of O(α b α s ) to the MSSM-like Higgs self-energies are not included in our calculation. However, in the definition of the bottom-Yukawa coupling we employ a running DR bottom mass at the scale m t [116] which entersˆ (1L) hh k 2 GNMSSM , and we take into account large t β -enhanced contributions to the bottom mass as discussed in Refs. [116,[135][136][137][138][139][140]. We expect that the missing two-loop piece of O(α b α s ) is numerically subleading (for a discussion in the MSSM see [141,142]).
Higher-order propagator-type corrections are not only needed for predicting the Higgs-boson masses, but also for the correct normalization of S-matrix elements involving Higgs bosons as external particles. The wave-function normalization factors incorporating the effects of the mixing between the different Higgs bosons can be written as a nonunitary matrix Z mix i j . It is constructed from the Higgs selfenergies and their derivatives with respect to k 2 , evaluated at the various physical poles; for details we refer the reader to Refs. [24,[143][144][145][146]. A recent application in the framework of the NMSSM can be found in Ref. [147]. Here, we follow the setup outlined in Section 2.6 of Ref. [24] and determine the matrix elements of Z mix i j from the eigenvalue equation The normalization of each eigenvector is fixed by ⎡ In our numerical analysis we denote the three CP-even mass eigenstates h i as h 0 , H 0 and s 0 , and the two CP-odd mass eigenstates a i as A 0 and a s . These assignments become ambiguous as soon as loop corrections are included. In our analysis we use the largest admixture to a loop-corrected mass state in order to define the assignment. For this purpose we employ the previously discussed loop-corrected mixing matrix Z mix i j . In this way s 0 denotes the dominantly singletlike state. The light doublet-like state is named h 0 and the heavy doublet-like state is H 0 . The CP-odd Higgs bosons are the predominantly singlet-like state a s and the doublet-like state A 0 .

Trilinear Higgs-boson self-couplings
In order to discuss possible distinctions between the NMSSM and the μNMSSM, the Higgs-boson self-couplings are particularly relevant. Experimentally these self-couplings can be probed through Higgs pair production or through decays of a heavier Higgs boson to two lighter ones. Through electroweak symmetry breaking there is also a strong correlation with Higgs-boson decays into Higgs bosons and gauge bosons, e.g. A 0 → Zh 0 or H 0 → Za s . For both, the Higgs mixing between singlets and doublets is essential. We take both types of decays into account when checking against experimental limits from Higgs boson searches, but only exemplify the parameter dependence for the decays involving only Higgs bosons in our numerical analysis below.
The Higgs self-couplings are introduced in Eq. (12). In order to simplify their presentation in the neutral sec- where in the CP-odd sector the Goldstone boson is in a basis where it does not mix with the other Higgs bosons at lowest order (see discussion in Sect. 2.2). 7 We denote the couplings as λ i jk for the interactions among three Higgs bosons φ i φ j φ k in the basis . For the couplings among the CP-even components -expressed in gauge couplings (see Eq. (13) for the relation to the gauge-boson masses) -we obtain at the tree level The couplings of CP-even components to CP-odd components are given by Similarly we can write down the couplingsλ i for the interaction φ i H + H − of the neutral Higgs bosons in the basis to the physical charged Higgs bosons (the Goldstone bosons are again in a basis where they do not mix) as follows: The remaining couplings which are not present above are equal to zero. Again s x and c x are defined as s x = sin(x) and c x = cos(x). In most of the cases when μ or μ eff appear, the coupling depends on the sum (μ + μ eff ). For the interactions of the neutral Higgs bosons, only a few couplings carry an (additional) proportionality to μ eff itself, see λ 123 , λ 345 and λ 333 which all involve the singlet state. This dependence manifests itself for the former two couplings in the Higgsto-Higgs decays In the charged Higgs sector, the decay s 0 → H + H − has a direct dependence on μ eff at the tree level in addition to (μ + μ eff ) for a dominantly singlet-like state s 0 , as can be seen inλ 3 . For both cases a very pronounced mixing of the singlet states with the Higgs doublets, and an individual dependence on μ eff and on the sum (μ + μ eff ) can also occur in other Higgs-to-Higgs decays. We will emphasize later that Higgs mixing is crucial for the observed dependences on μ eff and μ. We consider the decays at the tree level, however, including the external corrections to Higgs-boson masses and mixing as discussed in Sect. 2.4. Though, we emphasize that higher-order contributions to Higgs-boson self-couplings and Higgs-boson decays can be large, see Refs. [147][148][149][150] for corresponding calculations in the NMSSM.

Neutralino and chargino masses
We write the neutralino and chargino sector in the gaugeeigenstate bases which includes the bino componentB 0 , the neutral and charged wino componentsW 0 3 andW ± , the neutral and charged higgsino componentsh 0 u,d andh ± u,d , and the singlino components 0 in the form of Weyl spinors. Their mass terms in the Lagrangian can be written in the form The symmetric mass matrix of the neutralinos and the mass matrix of the charginos are given by The abbreviations s w = g 2 / g 2 1 + g 2 2 and c w = g 1 / g 2 1 + g 2 2 denote the sine and cosine of the weak-mixing angle, respectively. We see that the mass scale of the MSSMlike higgsinos is given by the sum (μ + μ eff ), and the mass scale of the singlino is controlled by (2 κ/λ μ eff + ν). If only the electroweakinos were taken into account at the tree level, it is apparent that the μNMSSM would be indistinguishable from the NMSSM, since any shift in masses and mixing induced through μ could be compensated through shifts in μ eff . However, such shifts will induce differences in the Higgs sector.
Including the singlino elements (with ν = 0 GeV as discussed in Sect. 2.1), an NMSSM-like neutralino spectrum can be generated, where (μ + μ eff ) serves as the NMSSMlike μ eff term and κ is rescaled as This rescaling on the other hand affects the Higgs spectrum, thus giving a possible handle to distinguish the μNMSSM from the NMSSM. For the case where κ and λ are kept fixed, an interesting behavior can be observed for light higgsinos. For small (μ + μ eff ) huge cancellations may occur between the two contributions with large μ > 0 GeV and μ eff of the same size but opposite sign. As a consequence, the singlino state becomes much heavier compared to the case of the NMSSM (of the order of μ eff ). Such a scenario is displayed in Fig. 1 where the neutralino-chargino spectrum is shown for the cases μ ∈ {0, 200, 1000} GeV (ν is set equal to zero). The left column with μ = 0 GeV corresponds to the case of the NMSSM. The masses are obtained by diagonalizing the tree-level mass matrices in Eq. (35). With respect to the Z 3invariant NMSSM, the most significant alteration is visible in the singlino component (blue): the mass shows an aboutlinear increase with μ since the sum (μ + μ eff ) is kept fixed. Due to the varying mixing, some influence on the masses of the other two neutral higgsino states (orange) can be seen despite a constant higgsino mass parameter (μ + μ eff ); the impact on the gaugino states (red and purple) remains negligible. The chargino masses (rose) are not influenced by the different choices.
In a scenario as discussed above, with light higgsinos as well as large μ and μ eff of opposite signs, the lightest neutralino is typically not the singlino state as the singlino mass is pushed up, see Fig. 1. The lightest supersymmetric particle (LSP), however, tends to be the gravitino, which is at risk to overclose the universe as dark matter candidate. In this case, the inflationary scenario has to be such that the reheating temperature stays below a certain value and gravitinos are not overproduced in the early universe, see our discussion in Sect. 2.1.

Sfermion masses
The mass term for each charged sfermion -for which we distinguish the superpartners of the left-and right-handed components by the notationf L andf R , respectively -takes the following form in the Lagrangian where the squared mass matrix reads Therein we denote the fermion mass by m f , the bilinear soft-breaking parameters by mf L,R , the trilinear soft-breaking parameter by A f , and the electric and weak charges by Q f and T (3) f . In this sector we encounter the sum (μ + μ eff ) in the off-diagonal elements of the sfermion mass matrices as the only difference compared to the NMSSM or MSSM. If this sum becomes large, A f /θ f needs to be adjusted in order to avoid tachyonic sfermions in particular for the third genera- The mostly bino-and wino-like statesB 0 (purple) andW 0 (red) as well as the charginosχ ± (rose) have (nearly) constant masses. The masses of the two mostly higgsino-like statesH 0 (orange) and the mostly singlino-like stateS 0 (blue) vary visibly tion squarks. In that case, bounds from vacuum stability (see e.g. Eq. (24)) can also constrain the viable size of (μ + μ eff ).

Phenomenological analysis
In this section we investigate various scenarios of the μNMSSM with a particular focus on the μ parameter. We will point out differences between the μNMSSM and the ordinary Z 3 -preserving NMSSM, where the latter corresponds to the limit μ = 0 GeV of the μNMSSM. At first we qualitatively define the investigated scenarios, before we numerically analyze them.

Viable parameter space compatible with theoretical and experimental bounds
In the previous sections we have analytically discussed the relevant sectors of the μNMSSM with respect to effects of the inflation-inspired μ parameter. Before we provide a phenomenological analysis -including the higher-order effects specified in Sect. 2.4 -we discuss the viability of various parameter regions. As discussed in Sect. 2.1 we focus on scenarios with non-zero μ and B μ , but set all other Z 3 -violating parameters in the superpotential (9) and soft-breaking Lagrangian (10), i.e. ξ , C ξ , ν and B ν , equal to zero.
The μ parameter of the model is positive by construction in the inflation-inspired model, see Eqs. (6) and (8). Furthermore, we only investigate scenarios with μ 2 TeV to stay in the phenomenologically interesting region for the collider studies. Still, we point out that also much larger scales are viable from the inflationary point of view. As discussed in Sect. small μ 1 GeV in the case of small μ also the softbreaking term B μ μ becomes small. Since in addition we set all other Z 3 -violating parameters to zero, we recover the standard NMSSM in this limit (see the discussion in Fig. 1). Thus, differences between the NMSSM and the μNMSSM can directly be deduced by comparing scenarios with zero and non-zero μ parameter. large μ 1 TeV with μ eff −μ as discussed in Sect. 2.6, the higgsino masses depend only on the sum (μ + μ eff ) at the tree level. The same combination contributes to the sfermion mixing in combination with the trilinear soft SUSY-breaking terms. In order to keep these quantities small at a large value of μ, one can assign the same value with opposite sign to μ eff ; note, however, that the region |μ + μ eff | 100 GeV is experimentally excluded by direct searches for charginos [151,152]. An immediate consequence of large, opposite sign μ eff and μ is that the singlino and the singlet-like Higgs states receive large masses of the order of |μ eff | [see the (5,5) (24) with μ promoted to (μ + μ eff )), have a much smaller impact in the μNMSSM as compared to the MSSM [108]. A special case of this scenario is the possibility of having μ at the electroweak scale in combination with an almost vanishing |μ eff | μ. This implies that (μ + μ eff ) remains at the electroweak scale. In contrast to the standard NMSSM this scenario allows the occurrence of both, κ λ and a light singlet sector. As discussed in Sect. 2.2, the mixing between singlets and doublets is in this case dominated by terms proportional to μ −1 eff . We will explicitly discuss such a scenario in Sect. 3.5.
There are more parameters that are relevant for the following phenomenological studies. We keep those fixed which behave similarly as in the MSSM and NMSSM. The choice of our constant input values is given in Table 1. Furthermore, we specify the values of t β , κ, λ, and A κ directly at the respective places. Besides the analyses where we explicitly study the dependence on B μ , we use B μ = 0 GeV as default value.
As our analysis is focused on the impact of the inflation model, we are not going to discuss the influence of the sfermion parameters. If not mentioned otherwise, we use mf ≡ mf L = mf R and A f 3 /mf = 2, which maximizes the prediction for the SM-like Higgs-boson mass at μ + μ eff = 0 GeV. The gluino mass parameter M 3 is set well above the squark masses of the third generation which is in accordance with the existing LHC bounds. For completeness, we also give the parameters of the SM which are most relevant for our numerical study in Table 1.
The gaugino-mass parameters M 1 and M 2 do not play a big role in the following analysis, but are necessary input parameters for the mass matrices of the charginos and neutralinos in Eqs. (35). We set M 2 = 500 GeV and fix M 1 via the usual GUT relation, see Table 1. Our phenomenological analysis is most sensitive to the neutralino and chargino spectrum if a Higgs boson can decay into them. This is in particular the case if the particle spectrum contains light higgsinos, whose masses are controlled through (μ + μ eff ). For a scenario with light higgsinos and a light singlino we will later also discuss the electroweakino phenomenology at a linear collider, see Sect. 3.4.
As we use μ 3 2 m 3/2 10 5 λ and focus on μ 2 TeV, we are considering scenarios where the gravitino typically is the LSP. We do not specify the mediator mechanism of SUSY breaking; however, we assume that such a light gravitino is always possible. Although the gravitino is the Dark Matter candidate, traditional collider searches for a neutralino LSP do apply in our case: for instance, if the nextto LSP (NLSP) is gaugino-like, it can decay into a photon and the gravitino, where the NLSP lifetime is typically so large that it can escape the detector [153]. We roughly estimate the NLSP phenomenology via the approximate partial decay width of the neutralino NLSP into a photon or Z boson and gravitino ψ 3/2 according to Refs. [154][155][156] where we expanded in a small gravitino mass m 3/2 and use s w and c w for the sine and cosine of the weak mixing angle, respectively. The neutralino mixing matrix elements N i j follow from the diagonalization of Eq. (35a). As an example for the decay of the NLSP with mχ0 1 100 GeV and m 3/2 10 MeV, we find a lifetime of τ ≡ 1/ = O(1 s). Thus, the NLSP decays outside of the detector and is counted as missing energy. Nevertheless, such decays might be of certain interest with respect to future experimental searches for long-lived particles like the MATHUSLA experiment [157]. Note that for a higgsino-like NLSP the decay into a Z boson and the gravitino is obtained by replacing the mixing factor in Eq. (39) by |−N 13 c β + N 14 s β | 2 . If kinematically open, also the decay into a (singlet-like) CP-even or CP-odd Higgs boson and the gravitino can occur (see Ref. [155]), but this decay mode does not change the qualitative features described above.
We have chosen m H ± as an input parameter and adjust A λ according to Eq. (19). If not denoted otherwise, we set m H ± = 800 GeV. We use HiggsBounds version 5.1.0beta [33][34][35][36][37] in order to implement the constraints on the parameter space of each of our scenarios resulting from the search limits for additional Higgs bosons. In this context, the exclusion limits from H, A → τ τ decays are particularly important. For relatively low values of tan β the choice of m H ± = 800 GeV is well compatible with these bounds. The code HiggsBounds determines for each parameter point the most sensitive channel and evaluates whether the parameter point is excluded at the 95% confidence level (C.L.). We use those exclusion bounds as a hard cut in the parameter spaces of our analyses.
We also indicate the regions of the parameter space which provide a Higgs boson that is compatible with the observed state at 125 GeV. These regions are obtained with the help of HiggsSignals version 2.1.0beta [38]. The code HiggsSignals evaluates a total χ 2 value, obtained as a sum of the χ 2 values for each of the 85 implemented observables. Four more observables are added, which test the compatibility of the predicted Higgs-boson mass with the observed value of 125 GeV. This latter test includes a theoretical uncertainty on the predicted Higgs-boson mass of about 3 GeV, such that a certain deviation from the four measured mass values (from the two channels with either a γ γ or a Z Z ( * ) final state from both experiments ATLAS and CMS) is acceptable. Thus, in total HiggsSignals tests 89 observables.
Since all our two-dimensional figures include a region with a SM-like Higgs boson, 9 we classify the compatibility 9 The minimal χ 2 value obtained in our numerical analysis is χ 2 m = 74.6. All subsequently discussed benchmark planes include a parameter region with χ 2 m < 80. Further details are provided below.  Table 1 with the observed state as follows: we determine the minimal value of χ 2 , denoted by χ 2 m , in the two-dimensional plane and then calculate the deviation χ 2 = χ 2 − χ 2 m from the minimal value in each parameter point. We allow for a maximal deviation of χ 2 < 5.99, which corresponds to the 95% C.L. region in the Gaussian limit. All parameter points that fall in this region χ 2 < 5.99 are considered to successfully describe the observed SM-like Higgs boson.
Lastly, we note that HiggsBounds and HiggsSignals are operated through an effective-coupling input. We will comment on the results of the two codes where appropriate.
For our implementation of the constraints from the electroweak vacuum stability we refer to Sect. 2.3. For informative reasons, we distinguish long-lived vacua from short-lived ones in the numerical analysis. We do not explicitly enforce a perturbativity bound on κ and λ, but discuss this issue below.

Higgs-boson and neutralino mass spectra
In this section, we point out the differences of the Higgsboson and neutralino mass spectra in the μNMSSM with respect to the NMSSM. Similar to the case of the MSSM, the charged and the CP-even heavy doublet as well as the MSSM-like CP-odd Higgs bosons are (for sufficiently large m H ± M Z ) quasi-degenerate. In Fig. 2, we show the masses of the Higgs bosons for vanishing A κ in the left, A κ = 100 GeV in the middle frame, and the masses of the neutralinos in the right frame. Each frame contains three different scenarios which are characterized by the three values μ ∈ {0, 200, 1000} GeV while keeping all other parameters fixed: μ + μ eff = −200 GeV, t β = 3.5, λ = 0.2, κ = 0.2 λ, and the other parameters as given in Table 1. The additional μ term has the biggest influence on the singlet-like states s 0 and a s , as well as the singlino-like stateS 0 . In analogy to the discussion in Fig. 1, the reason for this behavior is the fixed sum (μ + μ eff ): an increase in μ causes a larger negative μ eff which primarily drives the singlet-mass terms in the (3, 3) elements of Eqs. (17a) and (17b), and the singlino-mass term in the (5, 5) element of Eq. (35a) to large values. In the investigated parameter region, the mass of the CP-odd singlet is also very sensitive to A κ : in order to avoid a tachyonic state a s over a large fraction of the parameter space, it is essential to keep A κ sufficiently large. However, in the left frame a scenario is shown where even a vanishing A κ is possible. It generates a rather light CP-odd singlet-like state, whereas a sizable A κ = 100 GeV (middle) lifts this mass up. There is thus the potential for a distinction between the NMSSM-limit for μ = 0 GeV and the μNMSSM with a large μ = 1 TeV. Note that in the middle frame for μ = 200 GeV, the purple and blue lines are on top of each other.
The masses of the neutralino sector do not depend on A κ at the tree level. Concerning the Higgs sector, only the two cases in Fig. 2 with μ = 0 GeV and A κ ∈ {0, 100} GeV yield a SM-like Higgs boson that is compatible with the experimental data with χ 2 values of maximal 77. These two cases are also compatible with searches for additional Higgs bosons probed by HiggsBounds. The two cases with μ = 1 TeV and A κ ∈ {0, 100} GeV yield minimal χ 2 Fig. 3 In a similar manner as in Fig. 2, the spectra of Higgs bosons and neutralinos are shown in the μNMSSM. The neutralino masses are invariant under changes in μ by identifying the sum (μ + μ eff ) of the μNMSSM with the μ eff term of the NMSSM, and by rescaling κ according to Eq. (36). We set κ = 0.8 λ, and for μ = 0 GeV we assign μ eff = −200 GeV. The Higgs mass spectra are slightly affected by the rescaling values of 82.6 and 84.0, respectively. The larger values of χ 2 mainly arise because the SM-like Higgs-boson mass is slightly below 122 GeV. The large variation with μ for the mass prediction of the mostly SM-like Higgs boson is mainly induced by a large mixing with the CP-even singlet. The mixing for μ = 200 GeV in this scenario becomes very large for both values of A κ such that these cases are outside the parameter region that is compatible with the constraints by HiggsSignals. Note that the apparent preference for μ = 0 GeV over μ ∈ {200, 1000} GeV in this scenario is purely accidental and could be reversed by a slight shift in the input parameters, see the discussion below.
As already mentioned in Sect. 2.6, the electroweakino sector alone, at least at the tree level, does not allow one to distinguish the μNMSSM from the NMSSM: one can keep the neutralino-chargino spectrum at the tree level invariant by identifying the sum (μ + μ eff ) with the μ eff term of the NMSSM, and rescaling κ according to Eq. (36). However, as pointed out above, the rescaling does have an impact on the Higgs spectrum. We show in Fig. 3 spectra for μ ∈ {0, 200, 1000} GeV and A κ ∈ {0, 100} GeV with fixed μ + μ eff = −200 GeV. The neutralino spectrum is shown in only one column in the very right frame. In analogy to Fig. 2, the left and middle frames show the Higgs-boson masses for the two values of A κ where one still can see the effect of a varying μ term. While contributions to the mass matrices in Eqs. (17) which are proportional to (μ + μ eff ) or κ μ eff are kept constant, other terms ∝ μ −1 eff , μ −2 eff induce variations. Accordingly, the singlet-like Higgs masses in Fig. 3 are only slightly sensitive to μ, much less than the changes observed in Fig. 2. A rising μ slightly increases the mass splitting between the singlet-like and the SM-like Higgs state.
Still, while the Higgs masses remain almost constant for not too small μ eff , the doublet-singlet mixing can be strongly affected by varying μ and μ eff (but keeping their sum constant), in particular if the doublet-singlet mixing almost vanishes at a certain choice of μ and μ eff . In general, the mixing between the singlet and doublet states is affected by a large |μ eff |. However, by rescaling κ according to Eq. (36) all contributions linear in μ eff are absorbed, while the contributions ∝ μ −1 eff depend on the values of t β , M H ± and B μ , see Eqs. (20a) and (20b). 10 In Sect. 3.5 we will further investigate scenarios with very small μ eff and enhanced Higgs-boson mixing.
In Fig. 3 only the case A κ = 100 GeV in combination with μ = 0 GeV is allowed by HiggsBounds and HiggsSignals (χ 2 = 80.1), since the other scenarios are either ruled out by the decay of the SM-like Higgs into a pair of light CP-odd singlets or by a too large deviation of the SM-like Higgs-boson mass from 125 GeV. In addition to our discussion above, we emphasize that in particular the latter exclusion can be easily avoided through a slight adjustment of the input parameters.

Parameter scan
We have discussed above the dependence of the Higgs masses and of the condition for the stability of the electroweak vacuum on the model parameters. Apart from the fixed parameters in Table 1, we choose seven "free" parameters that we vary in the following regimes for our analyses: where the largest values of λ and κ in the specified range of (40) violate the approximate perturbativity bound λ 2 + κ 2 0.5. 11 For the results presented in the following, this bound is always fulfilled and lies outside the plot ranges. Values of tan β 4 push the model into the MSSM-like regime and are of less interest for studying the μNMSSM effects.
We have performed a scan over the parameter space defined in (40) and identified regions which are compatible with current observations concerning the properties of the SM-like Higgs boson at 125 GeV and the limits from searches for additional Higgs bosons with HiggsBounds 10 In the GNMSSM, there are further possibilities of absorbing shifts in μ eff through a redefinition of other Z 3 -violating parameters. 11 This perturbativity bound was explicitly derived for the NMSSM in Ref. [158]. According to the beta functions for λ and κ (see appendix A) no additional scale-dependent contribution is introduced by the μNMSSM at the one-loop order. and HiggsSignals as described above. In the following, we present a selection of results from this scan; different regions of vacuum stability are illustrated, and the experimental constraints from Higgs physics are indicated. While we display some typical examples, it should be noted that similar observations hold for other regions in the parameter space as well.
In Figs. 4, 5, 6, 7, we present a selection of parameter regions. Before we discuss them individually, their common features are explained. The colored dots in the background display different states of the electroweak vacuum: we distinguish stable (blue), long-lived meta-stable (purple), short-lived meta-stable (red), and tachyonic (rose). As discussed above, we regard not only tachyonic but also metastable regions as excluded in the context of this inflationary scenario, but nevertheless display long-and short-lived meta-stable regions for illustration. Furthermore, we indi-cate those points that do not fulfill Eq. (23) and thus have no singlet vev (orange), although, as explained in Sect. 2.3, this constraint is not relevant for the μNMSSM. We overlay mass contours for the SM-like Higgs h 0 (black), the CPeven singlet-like Higgs s 0 (blue), and the CP-odd singletlike Higgs a s (red). The spectrum is calculated taking into account the full one-loop and the known MSSM-like twoloop contributions as described in Sect. 2.4. The assignment of the labels h 0 , s 0 and a s to the loop-corrected states is determined by the largest respective contribution in the mixing matrix Z mix i j . We emphasize again that the parameters of the stop sector specified in Table 1 for the given scale of SUSY masses maximize the SM-like Higgs mass for μ + μ eff = 0 GeV; therefore, lower values for the SM-like Higgs mass could easily be obtained by reducing the mixing in the stop sector. Finally, we also indicate a naïve exclusion bound from direct searches for charginos Fig. 4 Contours for the SM-like Higgs mass (black) and the masses of the two singlet-like states (CP-even in blue and CP-odd in red) in the plane κ/λ versus (μ + μ eff ), where λ = 0.6 and μ = 500 GeV are kept fixed and κ and μ eff vary. In the left plot A κ = 0 GeV is used; in the right one A κ = 100 GeV. Furthermore, tan β = 2.5 is set in both plots. The other relevant parameters are listed in Table 1. The few red and purple points have a short-and long-lived meta-stable electroweak vacuum, respectively, whereas blue points have a stable electroweak vacuum. Rose points are excluded because of tachyonic tree-level masses. The orange points cannot reproduce a non-vanishing μ eff at the electroweak vacuum via the constraint of Eq. (23). With the gray vertical band we mark a naïve direct experimental exclusion bound from the chargino mass m χ ± > 94 GeV. Green areas are allowed by HiggsBounds and HiggsSignals (indicated as "HBHS" in the legend) Fig. 5 The same as Fig. 4, except that A κ = 100 GeV is used in both plots, and the parameter μ is set to μ = 1000 GeV (left) and 1500 GeV (right) Fig. 6 The same as Fig.4 but for μ = 1000 GeV, tan β = 3.5 and λ = 0.3 Fig. 7 Dependence of mass contours and vacuum stability, see Fig. 4 for an explanation of the color code, on the Z 3 -breaking soft SUSYbreaking B μ term and (μ + μ eff ) for λ = 0.5. On the left-hand side, the value A κ = 0 GeV was chosen, while on the right A κ = 100 GeV by the gray-shaded band: Ref. [152] reports a lower bound on the chargino mass of 94 GeV which translates into the requirement that |μ + μ eff | must be above that value in the μNMSSM. Lastly, all Figs. 4, 5, 6, 7 show the region of parameter points that successfully passed HiggsBounds and HiggsSignals and thus, in particular, yield a SM-like Higgs boson compatible with the observed state at 125 GeV. This region is represented through the larger, light green dots in the background. We refer to Sect. 3.1 for our statistical interpretation of the results obtained from the two codes.
A large part of the parameter region that is consistent with the measured SM-like Higgs mass is also in concordance with an absolutely stable electroweak vacuum. Small intersections between stable regions and regimes with tachyonic Higgs states exist, where there are meta-stable non-standard vacua. The strongest constraints arise from the existence of tachyonic masses for one of the physical Higgs states at the tree level. In the remaining region only a small fraction of points has a global minimum which does not coincide with the electroweak vacuum whereas the majority has a true electroweak vacuum. For the short-lived meta-stable regions, the vacuum lifetime is longer than the age of the universe.
In Fig. 4 we indicate the Higgs-mass contours and the constraints from vacuum stability in the plane of (μ + μ eff ) and κ/λ with fixed μ and λ. Note that for this choice of variables the tree-level doublet sector in Eqs.  Table 1; in the left plot A κ = 0 GeV is used, while in the right plot A κ = 100 GeV. The value of A κ has an impact in particular on the mass scale of the CP-odd singlet-like Higgs which is much lighter on the left-hand side. In fact, for a light CP-odd singlet-like Higgs a parameter region opens up where decays of the SM-like Higgs into a pair of them become kinematically allowed. The CP-even singlet-like Higgs is also somewhat lighter for A κ = 0 GeV, while the SM-like Higgs is scarcely affected. The contour lines of the Higgs masses stop when one Higgs becomes tachyonic. The reason why this does not exactly coincide with the border between the blue and pink dotted regions are the loop corrections to the Higgs spectrum while the constraints from vacuum stability were investigated at the tree level. It can be seen that the boundaries at the left of the stable region are parallel to one of the displayed Higgsmass contours -the corresponding particle becomes tachyonic at this boundary. The boundary to the right of the stable region can be understood when comparing the right plots of Figs. 4 and 5, which differ from each other by the value of μ: in the right plot of Fig. 5 a contour for the SM-like Higgs mass which is parallel to the tachyonic border appears around μ + μ eff = 250 GeV and κ/λ = 0.5. In Fig. 4 such a contour is not visible as this particular parameter region is excluded by a tachyonic SM-like state at the tree level. Note that the NMSSM-and μNMSSM-specific one-loop contributions to the Higgs spectrum are particularly large in that region (about 60 GeV additional shift compared to the same scenario in the MSSM-limit with λ → 0 and κ/λ constant), see also Ref. [124]; a dedicated analysis taking into account two-loop effects beyond the MSSM-limit might be necessary for a robust prediction of the Higgs mass close to the right border of the stable region, see e.g. Ref. [123]. It should be noted that in Fig. 4 the region where the Higgs mass is close to the right border of the stable region is disfavored by the limits from chargino searches at LEP.
As expected, the region allowed by HiggsBounds and HiggsSignals is a subset of the region where the SMlike Higgs has a mass in the vicinity of 125 GeV. In the green-marked region, χ 2 is at maximum 5.99. The minimal value χ 2 m from HiggsSignals is 74.6 in both figures. One can see on the left-hand side of Fig. 4 that this region is split into two: in between the two regions the SM-like Higgs can decay into a pair of CP-odd singlet-like Higgs bosons h 0 → a s a s with a branching ratio of up to 90 %; this behavior is not compatible with the observed signal strengths implying a limit on decays of the state at 125 GeV into non-SM particles. For a very light CP-odd singlet, the admixture between the SM-like Higgs and the CP-even singlet component is reduced, since the latter becomes heavier in this region. In the scenario under consideration, the decay h 0 → a s a s is dominated by the coupling among the two singlet states, λ 355 in Eq. (31q), such that a reduced admixture between h 0 and s 0 also closes the decay h 0 → a s a s . This is why -despite the very light CP-odd Higgs a s -the region at μ + μ eff −300 GeV and κ/λ 0.4 is allowed by the constraints from both HiggsSignals and HiggsBounds.
In Fig. 5 we present scenarios similar to the right-hand side of Fig. 4 with A κ = 100 GeV, but with different values of μ (note the larger scale at the x-axis). Thus, the influence of this parameter that distinguishes the μNMSSM from the NMSSM can be seen directly. Obviously, the parameter region with a stable vacuum is enlarged: for a given value (μ + μ eff ) the tachyonic border moves to smaller ratios of κ/λ as μ increases. Concerning the Higgs spectrum, the most notable difference is seen for the SM-like Higgs mass: for μ = 1 TeV a turning point at about μ + μ eff = −800 GeV is visible, which moves to smaller values of κ/λ for μ = 1.5 TeV. For the larger value of μ one can see that the possibility emerges for scenarios with the correct SM-like Higgs mass but positive (μ + μ eff ). Again all tested points which yield a SM-like Higgs boson close to 125 GeV successfully pass the constraints implemented in HiggsBounds and HiggsSignals. The minimal values of χ 2 m from HiggsSignals are 74.9 and 74.6 on the left-hand and on the right-hand side of Fig. 5, respectively. Figure 6 shows scenarios with larger tan β and smaller λ compared to the previous figures. Like in Fig. 4 we set A κ = 0 GeV on the left, and A κ = 100 GeV on the right-hand side, but μ = 1 TeV is used. We observe again that a larger value of A κ widens the allowed parameter region, because the mass of the CP-odd singlet is lifted up, giving rise to a drastic effect in this case. In fact, for A κ = 0 GeV only a rather small area in the plane of (μ + μ eff ) and κ/λ is allowed, while the allowed region is very significantly enhanced for A κ = 100 GeV. In the plot on the righthand side one can see a (nearly) closed 125 GeV contour for the mass of the SM-like Higgs with even larger values in the enclosed area. Adjusting the parameters of the stop sector in order to obtain a smaller contribution to the SMlike Higgs mass can render a SM-like Higgs with a mass of about 125 GeV in the whole enclosed region. Close to the tachyonic borders we find larger regions with a long-lived meta-stable vacuum (purple) than in Figs. 4 and 5. However, in this part of the plot the prediction for the mass of the SMlike Higgs is below the experimental value. On the right-hand side of Fig. 6 a large region is allowed by the constraints from HiggsBounds and HiggsSignals. Only low values of |μ + μ eff | < m h /2 are excluded by HiggsSignals due to the decay of the SM-like Higgs boson into a pair of higgsinos. However, this region is anyhow not compatible with the LEP bound on light charginos. The minimal values of χ 2 m from HiggsSignals are 74.7 in both plots. In Fig. 7 we change the parameter on the y-axis: B μ is varied and κ is kept fixed. We set A κ = 0 GeV on the lefthand side, and A κ = 100 GeV on the right-hand side. One can see that non-zero values for B μ can have a significant impact on the predicted Higgs masses and might determine whether or not a scenario is excluded. For larger negative values of B μ , one can see an area where the electroweak vacuum is meta-stable and long-lived, while the area in the lower left corner of the plots indicates that the electroweak vacuum is unstable and short-lived. The effect of a larger A κ mainly lifts the tachyonic boundary at the top so that values of B μ = 1 TeV are allowed for A κ = 100 GeV and leaves the other regions invariant. However, towards the upper limit of B μ , there is a small short-lived area. As a new feature, we find large regions with a meta-stable vacuum but a SMlike Higgs with a mass of 125 GeV for both values of A κ . Accordingly, scenarios with too large negative values of B μ are excluded due to a rapidly decaying vacuum despite providing a SM-like Higgs boson close to the observed mass. The constraints from HiggsBounds and HiggsSignals indicate that a large part of the region with the correct Higgs mass is compatible with the experimental data. For both plots HiggsSignals yields a minimal value of χ 2 m = 74.9. Only in those scenarios where the decay channel h 0 → a s a s is kinematically allowed -which happens in the plot for A κ = 0 GeV for μ + μ eff −300 GeV and μ + μ eff −700 GeV -the parameter region is incompatible with the data on the detected Higgs boson.
We briefly summarize the observed features and give an outlook for the phenomenological studies in the following. The allowed parameter region is mainly constrained by configurations where one Higgs field is tachyonic at the tree level. It can be seen that the tachyonic boundaries follow the Higgs mass contours in the Figs. 4, 5, 6, 7; in addition, there are effects from μ −1 eff terms as discussed in Sect. 2.2 which enhance the doublet-singlet mixing and eventually cause tachyons. This feature can be observed towards the right end of the Figs. 4, 5, 6. The experimental limits and constraints confine the allowed regions further around the region where the SM-like Higgs has a mass of about 125 GeV and exclude parameter regions where for instance the decay of the SM-like Higgs into a pair of light CP-odd singlets has a large branching ratio. In this context, the singlet sector has a significant impact on the features discussed in Figs. 4, 5, 6, 7. In the NMSSM, one usually expects to find the phenomenologically most interesting regions (accommodating a 125 GeV Higgs) for rather large values of λ 0.1, since the NMSSM contribution to the SM-like Higgs mass at the tree level is enhanced. In addition, large λ enhances the doublet-singlet mixing. However, in the μNMSSM, there is another way to obtain a large doublet-singlet mixing also for small values of λ: this is the region of low μ eff where terms proportional to μ −1 eff become large, as discussed in Sect. 3.1. We will investigate this class of scenarios, which are not possible in the NMSSM but generic to the μNMSSM, in Sect. 3.5 in more detail.
Similar to the NMSSM, the chosen value of A κ has a strong influence on the singlet-like Higgs masses, which is relevant for the tachyonic regions. In a large part of the viable parameter space the relation sign (A κ ) = − sign (μ eff ) applies, where for A κ = 0 GeV both signs of μ eff are allowed in general. This dependence on the relative signs of A κ and μ eff can be derived from the discussion in Sect. 2.2 about the Higgs singlets and especially the functional dependence of a 5 in Eq. (14e) versus a 4 in Eq. (20c): large negative values of the sum (a 4 + a 5 ) drive the CP-even singlet tachyonic. In the investigated scenarios above, which have either A κ = 0 GeV or A κ = 100 GeV, the sign of μ eff is negative in most of the viable parameter space. Accordingly, there is a preference for negative (μ + μ eff ). The allowed region with small positive values occurs where the negative value of μ eff is overcompensated by the positive value of μ. In Sect. 3.5 we will investigate a scenario where we keep (μ + μ eff ) fixed at a positive value, while for A κ small negative and small positive values are used for μ eff > 0 GeV and μ eff < 0 GeV, respectively. There we will also discuss the dependence of the singlet masses on μ and μ eff in more detail.

Higgs-boson and electroweakino production
In this and the next section we discuss phenomenological features of Higgs-boson mixing and thus consequences on Higgs-boson production and decays due to the μ parameter of the μNMSSM. For vanishing μ the phenomenology of the Higgs bosons equals the one of the NMSSM, for which typical benchmark scenarios can be found in Ref. [159] (see also Ref. [160]). Naturally they differ from MSSM-type benchmark scenarios through singlet states modifying the phenomenology: since the singlet states s 0 and a s neither directly couple to fermions nor to gauge bosons, but only through their admixture with the doublet states, their direct production -both at a hadron collider and a lepton collideris negligible in many scenarios. However, besides their direct production light singlet states can also be potentially observable via their production in cascade decays of heavier Higgs bosons, as we will discuss in the following.
In most parts of our numerical study, we make use of the approximation of SM-normalized effective couplings of a Higgs boson to gluons -calculated at leading orderwhich we insert into HiggsBounds for the evaluation of the Higgs-production cross-sections for the neutral Higgs bosons at the LHC. This treatment should be sufficiently accurate for determining the allowed regions in our scans over the parameter space. In the following, however, we will investigate to what extent the μNMSSM can accommodate the slight excesses in the data over the background expectation at a mass around 95-98 GeV that have been reported recently by CMS [42] in the γ γ channel 12 and earlier at LEP [41] in the bb channel. For this purpose we use more sophisticated predictions for the Higgs-production cross-sections in order to compare with the experimental results. We obtain those predictions from SusHi [39,40,[161][162][163][164][165][166][167], for which a dedicated version for the NMSSM exists [168]. The predictions include N 3 LO QCD corrections for the top-quark contribution of the light CP-even Higgs bosons, while we have neglected contributions from heavy squarks and gluinos beyond the resummed contributions in the bottom-Yukawa coupling.
In the NMSSM, the observed excesses in the data around 95-98 GeV can be interpreted in terms of a singletlike state s 0 , see Ref. [74] for a discussion of the LEP result, and Ref. [169] for a discussion of the CMS data. At first sight it seems to be non-trivial to describe both excesses simultaneously, since accommodating the LEP excess would require a rather large rate s 0 → bb, which in turn would suppress the channel s 0 → γ γ that is employed in the interpretation of the CMS excess. As it was pointed out in Ref. [147] based on a detailed analysis of the Higgs mixing properties, this is nevertheless possible -albeit in a relatively narrow region of the parameter space, which is somewhat enlarged if the possibility of non-vanishing phases giving rise to CP-violating effects is taken into account. We investigate in the following to which extent the additional freedom that is present in the μNMSSM with respect to the possible values of the masses in combination with the mixing properties has an impact regarding a possible interpretation of the observed excesses. In Table 2 we present four scenarios with s 0 masses in the range 95-98 GeV that have a phenomenology addressing the excesses observed both at LEP and CMS. Scenarios 1 and 3 have a small value of μ and are NMSSM-like (inspired by the scenarios investigated in Ref. [147]), while Scenarios 2 and 4 both have μ values that significantly differ from zero, and Scenario 4 furthermore has a non-zero value of B μ . These two μNMSSM scenarios are intrinsically different from the NMSSM. Similar scenarios could also be obtained by changing the signs of (μ + μ eff ) and A κ simultaneously.
Interpreting the LEP excess as the contribution of a singlet-like state s 0 in the considered mass range yields a "signal strength" of σ e + e − → Zs 0 · BR s 0 → bb while a "signal rate" of σ ( pp → s 0 → γ γ ) 0.1 pb would be compatible with the CMS observation. As mentioned above, the cross-section gg → s 0 in our analysis is obtained from SusHi [161,162] for the 13 TeV LHC at N 3 LO QCD. The renormalization-and factorization-scale uncertainties amount to about ±5%. Sizable values for the cross-sections gg → s 0 and e + e − → Zs 0 as well as the branching ratio BR(s 0 → bb) arise if the admixture of s 0 with the SM-like Higgs boson is sufficiently large. A sizable BR(s 0 → γ γ ) can occur as a consequence of a significant H u component of the singlet state s 0 , whereas a small H d component suppresses the decay into bb. In all the listed scenarios the CP-odd singlet-like Higgs boson a s has a mass below 300 GeV. It should be noted that the occurrence of the state s 0 at low masses in combination with a very heavy a s state through a large value of A κ would usually yield a meta-stable (long-lived) vacuum. The listed scenarios involve a certain amount of tuning in the choice of A κ since an increase in A κ by a few GeV yields a tachyonic s 0 state. It is well-known from the NMSSM that a too large A κ yields a Table 3 Cross-sections for electroweakinos at an electron-positron collider for Scenario 1 defined in tachyonic CP-even singlet-like Higgs boson s 0 , see Eq. (37) in Ref. [158] or Eq. (26) in Ref. [170] for lower and upper bounds on A κ . Similarly, we have noted a very pronounced dependence of the masses of both states, s 0 and a s , on A κ for the μNMSSM scenarios investigated here. Of course, a large admixture of s 0 with the SM-like Higgs boson in turn has an impact on the SM-like Higgs properties, visible through the increase in χ 2 returned by HiggsSignals. In fact, from the listed scenarios only Scenario 3 with χ 2 = 82 is compatible with the SM-like Higgs boson at the 95% C.L. The other scenarios have χ 2 values outside of the 95% C.L. region, as they have a slightly larger mixing of the singlet state with the SM-like Higgs boson. The enhanced mixing increases the s 0 cross-sections, but on the other hand yields reduced relative couplings to fermions and gauge bosons for the SM-like Higgs boson h 0 . It is thus apparent that explaining the excesses through a singlet state that only couples to SM particles through its admixture with the SM-like Higgs boson is under a certain tension from the measured SM-like Higgs-boson properties for both the μNMSSM and the NMSSM, if one requires signal rates that fully saturate the amount of deviation from the SM indicated by the excesses observed by LEP and CMS.
Scenarios with light singlet-like Higgs bosons tend to have a light singlino. For Scenario 1 we provide the light electroweakino spectrum, i.e. the masses ofχ 0 1,2,3 andχ ± 1 , in Table 3. Due to μ + μ eff = −140 GeV the scenario has light higgsino-like states, whereas the gauginos are close in mass to M 1 = 239 GeV and M 2 = 500 GeV. The higgsinolike states are strongly admixed with the singlino, e.g. the singlino-fraction ofχ 0 2 is 59%, the singlino-fraction ofχ 0 3 is 40%. It is apparent that the three lightest neutralinos and the light chargino are very close to each other in mass. At the LHC, ATLAS and CMS have only recently started to probe such compressed mass spectra by dedicated analyses, see e.g. Refs. [171,172]. In fact, an electron-positron collider may be required to ultimately probe scenarios of this kind, see for instance Ref. [173] tackling such compressed higgsino-like scenarios at the International Linear Collider (ILC). For Scenario 1 we provide the cross-sections for the two centerof-mass energies √ s = 350 GeV and √ s = 500 GeV, which are considered for Higgs-boson and top-quark precision studies at the ILC [174], in Table 3. Although in this scenario the LSP is the gravitino, the lightest neutralinoχ 0 1 has a lifetime of a few milliseconds such that it only gives rise to a missing-energy signature. Besides the possibility to tag e + e − →χ 0 1χ 0 1 through initial-state radiation (ISR), the production of one or more heavier neutralinos or charginos results in detectable SM particles. The possibility to polarize the initial state is an important tool to enhance the signalto-background ratio, and allows one to minimize systematic uncertainties. This capability is mandatory for performing precision measurements. In Table 3 we provide results for three different polarizations: an unpolarized initial state (as reference only), and polarizations of ±80% and ∓30% for the initial-state electron and positron, respectively. Such polarizations are foreseen in the current baseline design of the ILC. As one can see from Table 3 u . As discussed above, the electroweakino spectrum of the μNMSSM is a priori indistinguishable from the NMSSM if one restricts the analysis to information from the electroweakino sector and employs tree-level predictions, see Sect. 3.2. Previous studies of the electroweakino sector, see e.g. Ref. [173] and Refs. [176,177], discussed the ILC capabilities for distinguishing the MSSM from the NMSSM electroweakino sector. From such studies one can infer that a determination of the parameters of the electroweakino sector with an accuracy at the percent level is possible using the measurements of the light electroweakino masses and the corresponding production cross-sections for different polarizations, see e.g. Refs. [178,179] and references therein. This holds even if only the lightest electroweakinos are accessible. Based on earlier comprehensive studies where similar rates as in the scenarios of Table 3 were considered, the input parameters of the corresponding sector can be extracted: as an example, the values of M 1 and M 2 can be determined from the measurement of light gauginos, or the value of (μ + μ eff ) from the measurement of light higgsinos. In this regard beam polarization plays a crucial role: it allows one to even resolve scenarios where only a few light particles are kinematically accessible. Furthermore, the clean environment at an electron-positron collider allows the application of an ISR method [173] to detect and precisely measure scenarios where the light spectrum is close together in mass, as it is the case for instance for the compressed electroweakino spectrum in Scenario 1 leading to very soft decay characteristics. Complementing the particle spectrum via measuring additional heavier electroweakino masses and parts of the scalar and colored sector at the LHC would allow global fits of the model parameters, so that a model distinction between the μNMSSM, NMSSM and the MSSM might be feasible.

Higgs-boson mixing and decays
We now extend our previous discussion on Higgs-boson mixing and consider Higgs-boson decays. In this context we investigate in particular the influence of μ and μ eff on the masses of the two light singlets. For our discussion of the possibilities for distinguishing the μNMSSM from the NMSSM we assume that the sum (μ + μ eff ) is identified with the μ eff term of the NMSSM, and κ is rescaled according to Eq. (36). As discussed in Sect. 3.2, the rescaling of κ significantly reduces the dependence of all Higgs masses on μ and μ eff over a large region of the parameter space. Light CP-even singlets and decays of the SM-like Higgs boson into them were already part of the discussion in Figs. 4, 5, 6, 7, but therein we focused on constraints from vacuum stability and general features in the parameter space. Now we investigate Higgs-boson decays in more detail.
For the discussion of the dependence of the masses of the light singlet states on μ and μ eff we choose a scenario based on Table 1 and fix in addition λ = 1/4, κ = 1/5, tan β = 4, A κ = 7 GeV. We vary (μ+μ eff ) between −600 and 450 GeV and μ between 1 and 1000 GeV. The lower end of the range of μ values corresponds to the NMSSM-limit, μ → 0 GeV. The results are depicted in Fig. 8, where we show the masses of the three lightest neutral Higgs bosons as a function of μ and (μ + μ eff ). The background colors indicate the constraints from vacuum stability and from the experimental results on the Higgs sector using the same color coding as in Figs. 4, 5, 6, 7. It is apparent that with increasing μ the  Table 1 range in (μ + μ eff ) that is allowed by the constraints from vacuum stability is also increasing, which is in accordance to our observations in the previous section. Since the mixing between the light CP-even doublet state and the CP-even singlet-like state is large in the parts of the displayed parameter plane where the masses of the two lightest CP-even states are close to each other, we label the CP-even states in Fig. 8  We now focus on the region of large μ where both states can be clearly separated: the region with a SM-like Higgs mass of m h 0 1 ∼ 125 GeV is strongly affected by the values of λ and/or tan β through their impact on the NMSSMlike tree-level contribution to the doublet states, see the quantity a 1 in Eq. (20). In addition, it is well-known that this state receives large radiative corrections that depend on the mass splitting in the stop sector, which is proportional to X t = A t − (μ + μ eff )/ tan β. As discussed above, we have chosen A t in such a way that the contribution is maximized at μ + μ eff = 0 GeV and thus decreases to both directions. This behavior is visible in Fig. 8 for the contours displaying the mass of h 0 1 at values of μ 300 GeV.
The behavior of the singlet-like states is important for the phenomenology: as explained in Sect. 2.2, for a scenario with κ ∼ λ the mass of the CP-even singlet s 0 is mainly controlled by a 5 , see Eq. (14e), and therefore proportional to μ eff . As in our scenario A κ is small, the mass of the CPodd singlet a s is dominated by a 4 , see Eq. (20c), and therefore proportional to √ |μ/μ eff |. Those mass dependences of the singlet-like m h 0 2 150 GeV (blue) and m a s (red) can be clearly identified in Fig. 8: for m h 0 2 , the lines are roughly diagonal, and thus the mass contours follow lines with constant μ eff . For the mass of a s , the dependence on the square root √ |μ/μ eff | leads to the shape of the contours. We emphasize that the behavior displayed in Fig. 8 is specific to a small value of A κ μ or μ eff for κ ∼ λ. As above, the parameter range allowed by constraints from HiggsBounds and HiggsSignals is indicated by light-green dots in the background. For values of μ 300 GeV, the light CP-even Higgs h 0 1 corresponds to the SMlike state. In this region, the decay h 0 1 into higgsinos forbids low values of |μ + μ eff |, while the decay h 0 1 → a s a s is kinematically closed. For μ 500 GeV the mixing between the CP-even doublet state and the singlet state becomes larger which is not compatible with the observation of the properties of the SM-like Higgs. The minimal value of χ 2 in this figure is χ 2 m = 76.4 and thus slightly worse than the scenarios studied in Sect. 3.3. There is also a small allowed region with χ 2 < 5.99 at low values of μ, where a SM-like Higgs boson is present. In this region the mass of the singlet state s 0 has crossed the mass of the doublet state h 0 , and the states h 0 1 and h 0 2 have changed their character as discussed previously. The doublet-singlet mixing in this case yields a positive contribution to the mass of the SM-like state h 0 2 , lifting the treelevel value towards the experimentally allowed mass window (in the allowed region at low values of μ the mass of h 0 2 is about 126 GeV). In this region the decay h 0 2 → a s a s is kinematically open, but sufficiently suppressed to be in accordance with experimental observations. We conclude that the additional μ term of the μNMSSM lifts up the CP-odd Higgs mass and enlarges the allowed parameter space compared to the NMSSM. Still, in particular due to the large admixture of the singlet and doublet states, such a scenario is difficult to distinguish from the standard NMSSM, if not all Higgs states are fully determined. As a consequence of the strong admixture of the Higgs bosons and the influence of their masses on the kinematics, all decay modes show a non-trivial dependence on the coupling structure. The decay rates of the heavy Higgs bosons H 0 and A 0 into any combination of the three light Higgs bosons remain small throughout the parameter plane, i.e. the branching ratios are below 3%. The maximal branching ratios for A 0 → h 0 a s and A 0 → s 0 a s are reached at large μ and |μ eff |, i.e. in the lower right corner of Fig. 8. The two decays show a different dependence on μ and μ eff , which is in accordance with our discussion in Sect. 2.5. Whereas h 0 2 → a s a s is kinematically only allowed for very low μ in this scenario, the decay h 0 2 → h 0 1 h 0 1 is -when kinematically open -strongly dependent on μ eff . We will demonstrate below the dependence of the different decay modes on μ and μ eff in a scenario with essentially fixed Higgs-boson masses.
We now discuss a scenario that is intrinsically different from the NMSSM and shows a peculiar dependence of Higgs mixing and thus Higgs-boson decays on μ and μ eff . As indicated in the third item of Sect. 3.1 the μNMSSM allows large values of κ in combination with low values of μ eff and λ without being constrained by higgsino-like states. In the following we vary μ from 0 to 240 GeV and fix μ + μ eff = 160 GeV and thus simultaneously reduce μ eff from 160 to −80 GeV. We choose a very small value of λ = 0.02 and a value of κ = 0.02, which we rescale as κ →κ according to Eq. (36). The Higgs bosons therefore stay almost constant in mass, such that differences in Higgs-boson decays are solely induced by differences in the mixing among the Higgs bosons and not by kinematics. Note that in the limit μ eff → 0 GeV the rescaled parameterκ gets pushed beyond the perturbativity limit. This and the fact that at the tree-level one scalar mass becomes tachyonic are the reasons why the region around μ eff = 0 GeV is omitted for the lines in Figs. 9, 10, 11. Besides the parameters in Table 1 we set tan β = 4, A κ = − sign (μ effκ ) 1.3 GeV, which is like in Fig. 8 small compared to μ and |μ eff |, and B μ = 0 GeV. Note that the rescaling procedure for κ according to Eq. (36) turnsκ negative in this scenario, when μ + μ eff > 0 is fixed and μ takes on values larger than (μ + μ eff ). This case might be unsuitable for inflation, see Ref. [10]. The opposite case μ eff > 0 and μ + μ eff < 0 cannot appear in our model, since μ is always positive. The green-shaded area, which we show in the figures depicting SM-like Higgs properties, indicates the region that is compatible with the constraints from HiggsSignals, where as before we demand χ 2 < 5.99 with a minimum of χ 2 m = 77.5. For μ eff ∈ [−22, 22] GeV the decay h 0 → a s a s is enhanced such that this region is not compatible with the allowed fraction of non-SM decays of the SM-like Higgs boson. It should be noted that all of the shown area is allowed by the constraints from HiggsBounds. We emphasize that the effects that will be discussed in the following are related to a small value of μ eff : i.e. in the scenario discussed in Fig. 3, despite the same rescaling procedure for κ with respect to μ eff , the mixing among the Higgs states is much less influenced by the choice of μ, since a large value of μ results in an even larger negative value of μ eff .
In the standard NMSSM a measurement of the masses of the whole neutralino and neutral Higgs spectrum would fix all free parameters, in particular μ eff , λ, κ and A κ . With these parameters also the Higgs mixing is completely Fig. 9 The mass spectrum (left) and total decay widths (right) for the three lightest Higgs bosons in a scenario with small λ = 0.02 and fixed μ + μ eff = 160 GeV are shown. The value of κ is rescaled as κ → κ according to Eq. (36) with κ = 0.02. Furthermore, |A κ | = 1.3 GeV with the opposite sign of (μ effκ ). The other parameters are given in Table 1. The green area indicates compatibility with the constraints from HiggsSignals (right) Branching ratios of s 0 into non-SM particles and Higgs bosons; both for the same scenario as in Fig. 9. The green area indicates com-patibility with the constraints from HiggsSignals. The branching ratio BR(s 0 →χ iχ j ) includes all branching ratios of s 0 into pairs of neutralinos and charginos determined (at the tree level). A small value of λ in any case implies a small mixing between the singlet and doublet states of the Higgs sector. This is not the case in the μNMSSM: we show our results in Figs. 9, 10, 11. As explained above, in the considered parameter region the Higgs-boson masses are almost constant, see Fig. 9 on the left-hand side. The two heavy Higgs bosons H 0 and A 0 both have a mass very close to 800 GeV within a range of 3 GeV. The neutralino masses are constant, in detail mχ0 i = {134.7, 163.9, 252.1, 320.0, 516.1} GeV, where the particle with mass mχ0 3 = ms = 320 GeV corresponds to the singlino-like state with a purity of 99.9%. The two lightest neutralinos are higgsino-like states. Though the mixing in the neutralino sector remains constant, the mixing between the light CP-even Higgs boson h 0 and the singlet component s 0 is strongly enhanced for μ eff → 0 GeV. We depict the total widths for the three lightest Higgs bosons on the right-hand side of Fig. 9. The enhancement of the total width of h 0 for μ eff → 0 GeV is due to the decay h 0 → a s a s . This is also apparent in the left plot of Fig. 10, where the branching ratio for the decay of the SM-like state h 0 into a pair of light CP-odd singlets is displayed. For s 0 both the decays into h 0 h 0 and a s a s are of relevance, whereas other non-standard decay modes -e.g. into a pair of higgsinos -have a small rate, see the right-hand sides of Figs. 9 and 10. Apart from decays into Higgs bosons, s 0 decays Fig. 11 Branching ratios of A 0 (left) and H 0 (right) into pairs of lighter Higgs bosons; both for the same scenario as in Fig. 9. The branching ratios BR(A 0 →χ iχ j ) and BR(H 0 →χ iχ j ) include all kinematically allowed channels into pairs of neutralinos and charginos into massive SM gauge bosons. As mentioned above, all of the shown area in Fig. 10 is allowed by the constraints from HiggsBounds, i.e. both a s and s 0 are compatible with searches for additional Higgs bosons. However, for the state s 0 the region around μ eff = ±10 GeV is close to the boundary of the region that is excluded by the limits from Higgs searches, see below. As a result, we conclude that in this scenario with small μ eff the singlet s 0 can again be directly produced at a hadron collider through its admixture with the two CP-even doublets, see the discussion in Sect. 3.4. For μ = {150, 170} GeV the mass of the singlet is m s 0 = {323.9, 324.8} GeV, and the gluon-fusion production cross-section is σ (gg → s 0 ) = {270, 274} fb. The production rates through bottom-quark annihilation is negligible. Given the large branching ratios BR(s 0 → a s a s ) ∼ 57%, BR(s 0 → h 0 h 0 ) ∼ 19%, BR(s 0 → W + W − ) ∼ 15% and BR(s 0 → Z Z) ∼ 7%, the most sensitive searches are those with a decay into a pair of SM-like Higgs or gauge bosons. As an example, for m X ∼ 320 GeV the upper limits σ ( pp → X → h 0 h 0 ) 500 fb [180] and σ ( pp → X → Z Z) 200 fb [181] are already within a factor 10 of the signal rates that can be obtained at μ = {150, 170} GeV. Lastly, also the decays of the heavy Higgs bosons -whose total decay widths only vary within 10% for the considered scenario -show potentially observable branching ratios into pairs of lighter Higgs bosons in the limit μ eff → 0 GeV, see Fig. 11. At such low values of tan β both heavy Higgs bosons H 0 and A 0 are not predominantly decaying into bottom quarks or tau leptons, but decay into a pair of top quarks with a branching ratio of about 30%. Thus, decay modes into Higgs bosons could actually serve as discovery modes. However, note that our scenario includes light electroweakinos, into which heavy Higgs bosons tend to decay with large branching fractions. The branching ratios BR(A 0 →χ iχ j ) and BR(H 0 →χ iχ j ), shown in Fig. 11, both exceed 60% except for small values of μ eff . Both branching ratios include all kinematically allowed decays into pairs of neutralinos and charginos. This adds to the motivation for dedicated searches for heavy Higgs bosons decaying either into a pair of lighter Higgs bosons or into supersymmetric particles, see also the discussion in Refs. [182,183].
We conclude that a small value of μ eff in the discussed scenario strongly enhances the mixing among the Higgs bosons despite a low value of λ, which makes both singlet states potentially accessible at colliders. We have demonstrated that the Higgs-boson decays are not only controlled through the self-coupling dependences given in Sect. 2.5 for gauge eigenstates, but are also strongly dependent on the mixing of the Higgs bosons. In the standard NMSSM, light singlet states are usually associated with κ < λ, since the limits from chargino searches at LEP imply |μ eff | 120 GeV, and therefore v s 120 GeV. Accordingly, only small κ λ < 1 results in two light singlet states. However, in the μNMSSM scenario that we have considered κ/λ and μ can be large in combination with a small μ eff . Whereas the CP-odd singlet a s can be as light as a few GeV, the CPeven singlet s 0 is usually in the ballpark of a few hundred GeV in such scenarios. This scenario is intrinsically different from the behavior of Higgs masses and mixing known in the NMSSM.
We are left with a discussion of vacuum stability in this scenario of large values of κ together with small values of λ and μ eff , for which we consider a wider range of parameters, i.e. we allow for different values of μ + μ eff . Our results are shown in Fig. 12, where we varied both μ and μ+μ eff similar to Fig. 8 but for small λ = 0.02 and initial κ = λ. In contrast to Fig. 8 we rescale κ according to Eq. (36) in order to achieve a flat neutralino spectrum. Except at vanishing values of μ eff , which correspond to a diagonal line from (μ, μ + μ eff ) = (0, 0) GeV to (μ, μ + μ eff ) = (500, 500) GeV, also the Fig. 12 Vacuum stability analysis for a scenario with small value of λ, but potentially large values ofκ along the diagonal line of μ eff ≈ 0 GeV. The masses of the SM-like Higgs boson and the singlet-like CP-even and -odd states are indicated in black, blue and red, respectively. The green area is allowed by HiggsBounds and HiggsSignals; with the gray band, the direct exclusion bound on light charginos is shown. Superimposed with orange dashes is the line of constant μ + μ eff = 160 GeV, along which Figs. 9, 10, 11 are defined Higgs spectrum is almost constant throughout the plane due to the rescaling of κ. The behaviour discussed in Figs. 9, 10, 11 appears along the horizontal (orange, dashed) line at μ + μ eff = 160 GeV indicated in Fig. 12. Very similar results to the ones described in Figs. 9, 10, 11 are obtained for smaller or larger values of μ + μ eff > 0 GeV close to μ eff ≈ 0 GeV. Fig. 12 demonstrates that a large fraction of the (μ, μ + μ eff )-plane yields a stable vacuum. The vacuum appears to be unstable, but long-lived, in a valley around the diagonal line of μ eff ≈ 0 GeV discussed above as well as the lower left corner of the plot and along μ + μ eff 500 GeV. Along the same line and for small values of μ + μ eff tachyonic states rule out the parameter space. Those two bands, i.e. μ eff ≈ 0 GeV and μ + μ eff ≈ 0 GeV, are also the only regions in parameter space, which do not include a SM-like Higgs boson compatible with experimental measurements, indicated by the green area obtained with HiggsBounds and HiggsSignals. Two more comments are on order: In the upper right triangle of Fig. 12, i.e. μ + μ eff > 0 GeV and μ eff < 0 GeV,κ is negative, which might not be compatible with inflation. Moreover we want to point out that a very similar scenario can be found for even smaller λ, as long as the initial value κ ≈ λ is kept.

Conclusions
We have analyzed the phenomenology at the electroweak scale of an inflation-inspired extension of the Next-to-Minimal Supersymmetric Standard Model (NMSSM). We have put special emphasis on the spectra of additional, non-SM-like Higgs bosons and the branching ratios of their decays. This model has the same field content as the NMSSM, but at early times in the universe the D-flat direction of the Higgs doublet plays the role of the inflaton. Such a model can successfully describe inflation without the need of introducing a dedicated inflaton field. The singlet super-fieldŜ of the NMSSM is needed to stabilize the inflationary direction at the origin ofŜ = 0. Inflation occurs due to a non-minimal coupling of the doublet Higgs fields to gravity ∼ χ H u · H d , where the proportionality factor involves the gravitino mass m 3/2 at low energies. Thus, this model is characterized by an MSSM-like μ term, which is generated from the coupling χ and involves m 3/2 , in addition to the usual effective μ eff term of the NMSSM. The latter arises since the scalar component of the singlet superfield acquires a vacuum expectation value as in the NMSSM. At low energies, i.e. the electroweak scale, this model differs from the NMSSM by the additional μ term which breaks the accidental Z 3 symmetry of the NMSSM; we denote this model as the μNMSSM. The higgsino-mass term in the μNMSSM is composed of the sum (μ + μ eff ). We have classified and discussed various scenarios regarding the prospects to distinguish the μNMSSM from the NMSSM, where the latter corresponds to the limit μ = 0 GeV of the μNMSSM. We have derived constraints on the model parameters from theoretical and phenomenological considerations. For that purpose, we have computed the SM-like Higgs mass at the one-loop order in the μNMSSM and added as approximation at the two-loop level the known two-loop results from the MSSM which are implemented in FeynHiggs. We have probed our scenarios against the rate measurements of the SM-like Higgs boson and the limits from searches for additional Higgs bosons at colliders with the codes HiggsBounds and HiggsSignals. Furthermore, we have checked whether the electroweak ground state of the Higgs potential corresponds to the absolute minimum of the theory, i.e. the true vacuum, or whether the Higgs potential has a deeper non-standard minimum such that the electroweak vacuum eventually decays. In the inflationary scenario considered here, configurations with a meta-stable electroweak vacuum in general do not yield a viable phenomenology. In fact, the most stringent constraints arise from the possible appearance of tachyonic Higgs states at the tree level.
The additional freedom of varying μ and B μ in the μNMSSM allows one to choose values for the parameters of the NMSSM which would otherwise be excluded. In this extended parameter space, we have focused on relatively small values of tan β, since in this case -like in the NMSSMthe light doublet-like Higgs mass squared is increased by a shift ∝ λ 2 v 2 ; in this way the loop corrections which are required in order to acquire a SM-like Higgs at 125 GeV can be smaller. As expected, in particular the requirement of a SM-like Higgs boson at about 125 GeV yields impor-tant constraints on the parameter space. Concerning the constraints from vacuum stability, we find that the region with a phenomenologically viable Higgs spectrum is strongly correlated with the region of a stable electroweak vacuum, where the electroweak ground state corresponds to the true vacuum at the electroweak scale. An exception is the case where the soft SUSY-breaking B μ μ term is large. We have demonstrated that large negative values of B μ μ destabilize the vacuum.
For most of the numerical analyses in this paper we have fixed the sum (μ + μ eff ), since μ enters at the tree level only in this combination in the mass matrices for the charginos and sfermions as well as in the MSSM-like part of the neutralino mass matrix. Accordingly, the particle spectrum of the μNMSSM in those sectors resembles the one of the NMSSM if the sum (μ + μ eff ) in the μNMSSM is identified with the μ eff term of the NMSSM. Moreover, we have pointed out the possibility to further reduce the influence of the non-minimal coupling to supergravity ∼ μ on the neutralino sector by a rescaling of the parameter κ. This rescaling compensates the dependence of the singlino component of the neutralino mass matrix on μ eff , so that the neutralino, chargino and sfermion sectors of the μNMSSM and the NMSSM become indistinguishable from each other at tree level. We have demonstrated that the dependence of the Higgs masses on μ is significantly weakened after this transformation, but the individual dependences on μ and μ eff still have a large impact on the Higgs mixing and thus the branching ratios of Higgs decays. The modified value of κ resulting from the rescaling can also have an important influence on Higgs phenomenology.
Since with the above parameter settings the neutralino sector of the μNMSSM is NMSSM-like, we have not performed a detailed numerical analysis of the neutralino sector -besides our discussion of Higgs decays into electroweakinos. In general, the gravitino is found to be the LSP since it is tightly connected to the size of μ. For phenomenological reasons in our scenarios it typically has a mass of O(10 MeV). The NLSP, which is either singlino-or bino-like, tends to be sufficiently long-lived such that it only gives rise to missingenergy signatures in collider searches. Accordingly, typical constraints from SUSY searches including missing energy apply without large modifications. The character of the NLSP is influenced by a variation of the corresponding parameters, i.e. (μ+μ eff ) for the higgsino mass and M 1 or M 2 for the bino or wino mass, respectively. Our choices for M 1,2 are rather arbitrary in this context. Their impact could be scrutinized with a dedicated study of the neutralino phenomenology in the μNMSSM.
In some of our analyses we have kept λ large in order to lift up the mass of the SM-like Higgs boson at the tree level through genuine NMSSM effects, and in order to allow for sizable doublet-singlet mixing. However, we emphasize that large mixing between the doublet and singlet fields can also be achieved through small λ in combination with nearly vanishing μ eff . Such a scenario is viable in the μNMSSM and gives rise to a phenomenology that significantly differs from the NMSSM.
A phenomenologically very interesting set of scenarios includes light singlet states. The direct production of these states at colliders suffers from their nature as gauge singlets: couplings to SM particles only emerge through the admixture with doublet-like Higgs states. Similarly, Higgsto-Higgs decays involving doublet and singlet fields are strongly correlated with Higgs mixing. We have shown this effect exemplarily for decays of the SM-like Higgs boson into a pair of light CP-odd singlets, which depends on the fraction of the CP-even singlet component in the SM-like Higgs boson h 0 . In the μNMSSM, this mixing is not only controlled through λ, but also depends sensitively on the values of μ and μ eff . We conclude that in order to distinguish the Higgs sectors of the μNMSSM and the NMSSM, the detection of singlet states in the Higgs spectrum and their couplings to other Higgs bosons and the SM particles will be crucial. We have discussed four scenarios that yield a light CPeven singlet-like Higgs around 97 GeV, motivated by slight excesses in experimental searches performed with CMS and at LEP. These scenarios are associated with a compressed spectrum of light electroweakinos. We have pointed out that searches at a future electron-positron collider would provide complementary information to the results achievable at the LHC in scenarios of this kind.