Next-to-Minimal Supersymmetric Standard Model

• NMSSM is a supersymmetric extension of the Standard Model that introduces a gauge-singlet chiral superfield to dynamically solve the μ-problem and modify the Higgs and neutralino sectors.
• It enhances the tree-level Higgs mass via a λ-dependent quartic coupling and induces rich mixing among three CP-even Higgs states, two CP-odd Higgs states, and five neutralinos.
• The model yields distinct experimental signatures including non-standard Higgs decays and singlino/mixed dark matter candidates, offering new avenues for collider and cosmological studies.

The Next-to-Minimal Supersymmetric Standard Model (NMSSM) is a supersymmetric extension of the Standard Model in which the Higgs sector is enlarged relative to the Minimal Supersymmetric Standard Model (MSSM) by one gauge-singlet chiral superfield $\hat S$. In its simplest $Z_3$-invariant or “scale-invariant” form, the explicit MSSM $\mu$-term is absent and is replaced by the trilinear coupling $\lambda \hat S \hat H_u \hat H_d$, so that the singlet vacuum expectation value generates an effective Higgsino mass parameter $\mu_{\text{eff}} = \lambda \langle S\rangle$. This simultaneously modifies the Higgs and neutralino sectors, yielding three CP-even Higgs bosons, two CP-odd Higgs bosons, one charged Higgs pair, and five neutralinos including a singlino. More general $R$-parity-conserving singlet extensions allow explicit $\mu$, $\mu'$, and tadpole terms, but the standard NMSSM is usually identified with the $Z_3$-symmetric realization (Teixeira, 2011, Allanach et al., 2013).

1. Defining structure and theoretical motivation

In the MSSM, the renormalizable, $R$-parity-conserving superpotential contains the supersymmetric Higgs mass term $Z_3$0. Phenomenology requires $Z_3$1, but there is no symmetry reason for $Z_3$2 to be tied to the SUSY-breaking scale. The NMSSM addresses this by adding a singlet and replacing the explicit $Z_3$3-term with dimensionless couplings. In the simplest scale-invariant form, the superpotential is

$Z_3$4

while the Higgs-sector soft terms contain

$Z_3$5

When the scalar component of $Z_3$6 develops a vacuum expectation value $Z_3$7, one obtains $Z_3$8, so the only mass scale in the superpotential is replaced by a quantity generated after SUSY breaking (Teixeira, 2011).

A more general $Z_3$9-parity-conserving singlet extension permits

$\mu$0

and the Next-to-Minimal SOFTSUSY implementation likewise treats both the $\mu$1-invariant case and the case with general $\mu$2-violating terms in the superpotential and soft sector (Allanach et al., 2013). This distinction is conceptually important. In the $\mu$3-invariant NMSSM, the $\mu$4-problem is addressed dynamically, whereas in the S-MSSM, which retains an explicit $\mu$5-term and a supersymmetric singlet mass $\mu$6, the singlet is used primarily to alleviate the little hierarchy problem rather than to solve the $\mu$7-problem (Delgado et al., 2010).

A recurrent conceptual issue is that the $\mu$8 symmetry often imposed to obtain the scale-invariant superpotential can lead to domain walls; realistic models must break this symmetry, usually by tiny non-renormalizable terms. This has motivated systematic study of both $\mu$9-symmetric and $\lambda \hat S \hat H_u \hat H_d$0-violating variants (Teixeira, 2011, Allanach et al., 2013).

2. Vacuum structure, electroweak symmetry breaking, and particle content

After electroweak symmetry breaking,

$\lambda \hat S \hat H_u \hat H_d$1

with $\lambda \hat S \hat H_u \hat H_d$2 and $\lambda \hat S \hat H_u \hat H_d$3. In the standard $\lambda \hat S \hat H_u \hat H_d$4-invariant model, the singlet vev generates $\lambda \hat S \hat H_u \hat H_d$5, while in the more general formulation one defines $\lambda \hat S \hat H_u \hat H_d$6 and related effective bilinear quantities for the Higgs sector (Teixeira, 2011, Allanach et al., 2013).

The physical Higgs spectrum consists of three CP-even states, two CP-odd states, and one charged Higgs pair. In the CP-even sector, a convenient basis is

$\lambda \hat S \hat H_u \hat H_d$7

related to the physical states $\lambda \hat S \hat H_u \hat H_d$8 by a real orthogonal rotation

$\lambda \hat S \hat H_u \hat H_d$9

The state $\mu_{\text{eff}} = \lambda \langle S\rangle$0 is frequently identified with the observed Higgs boson near $\mu_{\text{eff}} = \lambda \langle S\rangle$1 GeV in coupling-fit analyses (Buttazzo, 2014). In the low-$\mu_{\text{eff}} = \lambda \langle S\rangle$2 regime, however, the SM-like Higgs can be either the lightest or the second-lightest CP-even state, whereas an SM-like heaviest CP-even scalar is described as unlikely (Christensen et al., 2013).

The CP-odd sector contains a doublet pseudoscalar and a singlet pseudoscalar. In a commonly used approximation,

$\mu_{\text{eff}} = \lambda \langle S\rangle$3

so $\mu_{\text{eff}} = \lambda \langle S\rangle$4 measures the doublet fraction of the lightest CP-odd state. For small $\mu_{\text{eff}} = \lambda \langle S\rangle$5, the lightest pseudoscalar can be predominantly singlet, with suppressed couplings to Standard Model matter (Almarashi et al., 2010).

The neutralino sector is enlarged from four to five states. In the basis $\mu_{\text{eff}} = \lambda \langle S\rangle$6, the mass matrix contains the usual bino, wino, and higgsino entries, together with the singlino mass term $\mu_{\text{eff}} = \lambda \langle S\rangle$7 and higgsino–singlino mixings proportional to $\mu_{\text{eff}} = \lambda \langle S\rangle$8. This makes possible bino-like, higgsino-like, singlino-like, or mixed lightest neutralinos, depending on $\mu_{\text{eff}} = \lambda \langle S\rangle$9, $R$0, $R$1, and the gaugino masses (Teixeira, 2011, Cheng et al., 2013).

3. Higgs-mass enhancement, mixing regimes, and naturalness

A central structural difference from the MSSM is the tree-level Higgs-mass enhancement. In the NMSSM, the “doublet-only” SM-like Higgs mass parameter is

$R$2

where the first term is the usual MSSM tree-level contribution, the second is the NMSSM tree-level enhancement from $R$3, and $R$4 encodes the leading top–stop loop correction (Buttazzo, 2014). The familiar tree-level bound is therefore modified to

$R$5

up to mixing effects and loop corrections (Teixeira, 2011, Ellwanger, 2011).

This modification underlies most naturalness arguments in the NMSSM literature. A sizable $R$6 and moderate $R$7 permit a $R$8–$R$9 GeV Higgs with smaller $\mu$0, corresponding to lighter stops than in the MSSM. In the “natural scalars” analysis, $\mu$1 GeV is taken as a representative value, corresponding to average stop masses around $\mu$2 GeV, and the “most natural” region is characterized by $\mu$3, moderate $\mu$4–$\mu$5, stop masses around or below $\mu$6 GeV, and Higgs mixing angles small enough to respect LHC Higgs coupling measurements (Buttazzo, 2014).

Two limiting regimes are especially useful. In the “singlet decoupled” case, the extra CP-even state is doublet-like and the phenomenology resembles that of the heavier CP-even MSSM Higgs. In the “doublet decoupled” case, the extra CP-even state is singlet-like, its couplings to Standard Model fields are uniformly scaled by the singlet–doublet mixing angle, and the decay $\mu$7 can be dominant when kinematically open (Buttazzo, 2014). Complementarily, low-$\mu$8 scans show that all Higgs bosons can remain relatively light, near or below the electroweak scale, and that Higgs-to-Higgs decays such as $\mu$9 or $\mu'$0 can strongly alter the pattern of visible branching fractions (Christensen et al., 2013).

The relation between singlet mixing and Higgs-mass enhancement is subtle. A mostly SM-like Higgs above $\mu'$1 GeV is possible in the NMSSM, but in the scenarios reviewed in “Higgs Bosons in the Next-to-Minimal Supersymmetric Standard Model at the LHC” it is necessarily accompanied by a lighter state with a large gauge-singlet component (Ellwanger, 2011). By contrast, in the S-MSSM the singlet may raise the light Higgs mass to $\mu'$2–$\mu'$3 GeV with sub-TeV stops while the singlet itself decouples and the low-energy phenomenology remains close to the MSSM; this illustrates that not all singlet extensions address the same structural problem (Delgado et al., 2010).

4. Collider phenomenology and characteristic search channels

The NMSSM Higgs sector permits production and decay patterns not available in the MSSM. Higgs coupling fits constrain the mixing of the observed $\mu'$4 GeV state with both the second doublet and the singlet. Using

$\mu'$5

one finds that current data already give strong constraints on doublet mixing, while constraints on singlet admixture are weaker (Buttazzo, 2014).

Direct searches bifurcate according to whether the relevant state is doublet-like or singlet-like. In the singlet-decoupled case, the extra CP-even scalar has couplings analogous to the MSSM heavy Higgs and decays mainly into $\mu'$6 or $\mu'$7, with vector-boson decays suppressed by the small mixing angle required by Higgs coupling data (Buttazzo, 2014). In the doublet-decoupled case, the singlet-like scalar has Standard-Model branching ratios rescaled by a universal mixing factor, and the decay $\mu'$8 can dominate. The most promising collider signatures are then di-Higgs final states such as $\mu'$9, $Z_3$0, and possibly $Z_3$1 (Buttazzo, 2014).

Light CP-odd Higgs phenomenology is especially distinctive. In the analysis of low-mass Higgs signals at the LHC, the lightest CP-odd Higgs $Z_3$2 can be singlet-like with a tiny doublet component in large regions of parameter space, yet still be produced in association with $Z_3$3 at large $Z_3$4. The paper finds that an $Z_3$5 with mass $Z_3$6 can be extracted from Standard Model backgrounds by using the $Z_3$7 decay channel, a possibility precluded to the MSSM, whereas the $Z_3$8 decay mode is overwhelmed by backgrounds despite the fact that the branching ratio of this mode can reach unity when $Z_3$9 is a pure singlet (Almarashi et al., 2010). This tension between production and branching ratio recurs throughout NMSSM light-Higgs searches.

Broader low-mass Higgs scans show that the SM-like Higgs may be either the lightest or second-lightest CP-even scalar, that the correlations of $R$0 and $R$1 can be substantially altered, and that non-standard decays into Higgs pairs often dominate when kinematically accessible (Christensen et al., 2013). Related reviews emphasize that Higgs-to-Higgs decays can spoil simple “no-lose theorem” expectations familiar from the MSSM, while also creating channels such as $R$2, $R$3, or $R$4 that are specific probes of the NMSSM (Teixeira, 2011, Ellwanger, 2011).

Charged-Higgs phenomenology supplies an additional discriminator. In the scale-invariant NMSSM, a relatively light charged Higgs can decay to a singlet-like pseudoscalar and a $R$5 boson,

$R$6

and the PYTHIA-FastJet simulation of “Non-standard charged Higgs decay at the LHC in Next-to-Minimal Supersymmetric Standard Model” finds that such a scenario can be probed with early data of $R$7 fb$R$8 at $R$9 and $Z_3$00 TeV. Because this decay evades the recent bounds on charged Higgs from the LHC, it exemplifies how MSSM-motivated search strategies can miss viable NMSSM Higgs sectors (Bandyopadhyay et al., 2015). Complementarily, the study “Light Higgs Bosons in NMSSM at the LHC” emphasizes that, although $Z_3$01 and $Z_3$02 decay modes appear to be the most promising, for a substantial region of parameter space the two-photon decay mode has a remarkably large rate, which can be exploited to find the NMSSM Higgs signal and can also be a potential avenue to distinguish the NMSSM from the MSSM (Guchait et al., 2015).

5. Dark matter, low-energy observables, and early-universe applications

The singlino extends neutralino dark matter beyond MSSM bino–higgsino–wino scenarios. The lightest neutralino can be bino-like, higgsino-like, singlino-like, or a mixture, and the extra Higgs states open new annihilation channels and resonances such as $Z_3$03, $Z_3$04, or $Z_3$05. This permits viable relic densities in regions where the MSSM would overclose the Universe and, in some constrained realizations, yields an almost pure singlino LSP with extremely small direct and indirect detection cross sections (Teixeira, 2011).

The electroweak-supersymmetry analysis of the NMSSM studies three scenarios—$Z_3$06-parity conservation with one dark matter candidate, $Z_3$07-parity conservation with multi-component dark matter, and $Z_3$08-parity violation—and obtains the minimal $Z_3$09 of $Z_3$10, $Z_3$11, and $Z_3$12, respectively. After LHC neutralino/chargino and slepton searches are imposed, the majority of viable parameter space preferred by the muon anomalous magnetic moment is excluded except for the parameter space with moderate to large $Z_3$13 ($Z_3$14). The most favorable parameter space has relatively large $Z_3$15, moderate $Z_3$16, small $Z_3$17, heavy squarks/gluino, and the second lightest CP-even neutral Higgs boson with mass around $Z_3$18 GeV (Cheng et al., 2013).

Cosmological applications go beyond dark matter. In the $Z_3$19-invariant, scale-invariant NMSSM, the conditions for a strongly first-order electroweak phase transition can be analyzed semi-analytically, and parameter regions can be found that simultaneously yield successful baryogenesis criteria, the observed relic density for the neutralino LSP, and inflation driven by a gauge invariant MSSM flat direction made up of right-handed squarks, while remaining constrained by the recent Higgs mass bound, branching ratios of rare flavour-violating decays, and the invisible $Z_3$20 decay width (Balazs et al., 2013). This places the NMSSM in a broader class of supersymmetric frameworks in which the same weak-scale parameters enter collider observables, dark matter phenomenology, and early-universe dynamics.

High-scale boundary conditions also matter. In the $Z_3$21-invariant NMSSM with universal scalar masses and trilinear couplings, non-universal gaugino masses can significantly relax the restrictive vacuum conditions that arise under universality at the unification scale. In that setup, a higgsino can be the lightest SUSY particle and its thermal relic abundance can reproduce the observed dark matter density in a wide range of parameter space in which the $Z_3$22 GeV Higgs-boson mass is obtained. The resulting higgsino-like dark matter may be probed in direct detection experiments, and many model points predict colored particles such as the gluino to be within the reach of a future $Z_3$23 TeV collider (Kawamura et al., 2018).

6. Constrained realizations, computational implementations, and enduring issues

A particularly economical realization is the fully constrained NMSSM (cNMSSM), in which universal high-scale soft terms are assumed: $Z_3$24 Electroweak symmetry breaking and phenomenological constraints reduce the effective freedom dramatically. In the cNMSSM, $Z_3$25 must be small so that the singlet soft mass runs little and a non-zero singlet vev can be obtained; the resulting LSP is almost pure singlino, nearly degenerate with a stau NLSP, and relic density typically requires singlino–stau co-annihilation. LEP Higgs constraints require $Z_3$26, and the viable parameter space becomes essentially one-dimensional, largely controlled by $Z_3$27, with $Z_3$28 favoured by collider bounds and $Z_3$29 (Teixeira, 2011).

Because NMSSM phenomenology depends sensitively on extended Higgs and neutralino sectors, dedicated spectrum and decay tools are required. NMSDECAY computes widths and branching ratios of sparticle decays in the NMSSM, generalizing SDECAY to include the extended Higgs and neutralino sectors and adding slepton 3-body decays, possibly relevant in the case of a singlino-like lightest supersymmetric particle. It is part of the NMSSMTools package, which computes Higgs and sparticle masses and Higgs decays in the NMSSM (Das et al., 2011). Next-to-Minimal SOFTSUSY extends SOFTSUSY to calculate the sparticle spectrum in both the $Z_3$30-symmetric and general $Z_3$31-violating NMSSM, solving the renormalisation group equations numerically between the weak scale and a high energy scale using a nested iterative algorithm (Allanach et al., 2013).

Several enduring issues remain intrinsic to the model. One is practical rather than formal: if the singlet decouples, the NMSSM can mimic the MSSM closely, making experimental discrimination difficult outside the Higgs sector (Teixeira, 2011). Another is structural: if Higgs-to-Higgs decays dominate and singlet admixtures are significant, a universal “no-lose theorem” for Higgs discovery does not hold in the same simple form as in the MSSM (Teixeira, 2011, Ellwanger, 2011). A further issue is conceptual: the $Z_3$32-invariant superpotential solves the $Z_3$33-problem elegantly but introduces the domain-wall problem, motivating $Z_3$34-breaking completions (Teixeira, 2011, Allanach et al., 2013).

Taken together, these studies establish the NMSSM as the simplest singlet extension of the MSSM that dynamically generates $Z_3$35, relaxes the MSSM Higgs-mass tension through the $Z_3$36-dependent quartic interaction, enlarges the Higgs and neutralino sectors in phenomenologically consequential ways, and supports a broad range of collider, dark-matter, and cosmological applications. Its distinctive signatures are not a single prediction but a family of correlated possibilities: light singlet-like scalars or pseudoscalars, singlino-dominated or higgsino-dominated dark matter, non-standard Higgs decays, altered Higgs coupling patterns, and constrained or non-universal high-scale realizations with sharply different low-energy spectra (Teixeira, 2011, Buttazzo, 2014).