GNMSSM: General Next-to-Minimal Supersymmetric Model

• GNMSSM is a supersymmetric extension of the NMSSM that introduces a gauge-singlet superfield and explicit bilinear and tadpole terms to resolve domain wall and tadpole issues.
• It exhibits rich Higgs and neutralino phenomenology, enabling controlled singlet-doublet mixing to reconcile the 125 GeV SM-like Higgs with a 95 GeV singlet excess and muon g-2 anomalies.
• The model offers viable dark matter scenarios through a singlino-dominated LSP while ensuring a natural, stable vacuum and compliance with cosmological and collider constraints.

The General Next-to-Minimal Supersymmetric Standard Model (GNMSSM) is a supersymmetric extension of the Minimal Supersymmetric Standard Model (MSSM) in which the discrete $\mathbb{Z}_3$ symmetry of the standard Next-to-Minimal Supersymmetric Standard Model (NMSSM) is lifted. The GNMSSM augments the theory with a gauge-singlet chiral superfield $\hat S$ and allows all renormalizable, gauge-invariant, and $R$- and $CP$-conserving superpotential terms. This flexible theoretical structure is constructed to remedy cosmological and ultraviolet problems endemic to the $\mathbb{Z}_3$-NMSSM, enables a rich Higgs and neutralino phenomenology, and naturally admits scenarios reconciling experimental anomalies—such as the muon $g-2$ discrepancy and low-mass scalar excesses—with dark matter, Higgs, and collider constraints (Cao et al., 2024, Li et al., 9 Nov 2025, 0910.1785, Cao et al., 2022, Cao et al., 2023, Meng et al., 2024, Cao et al., 2022).

1. Model Structure and Lagrangian

The GNMSSM Lagrangian is defined by extending the MSSM to include a gauge-singlet superfield $\hat S$. The general renormalizable superpotential in the Higgs-singlet sector is

$$W_{\rm GNMSSM} = W_{\rm Yukawa} + \lambda\,\hat S\,\hat H_u\cdot\hat H_d + \frac{\kappa}{3}\,\hat S^3 + \mu\,\hat H_u\cdot\hat H_d + \frac{\mu_S}{2}\hat S^2 + \xi_F\,\hat S\,,$$

where:

• $W_{\rm Yukawa}$: MSSM quark and lepton Yukawa couplings,
• $\lambda,\kappa$: dimensionless singlet-doublet and singlet self-couplings,
• $\hat S$0: supersymmetric Higgsino mass,
• $\hat S$1: supersymmetric singlet mass,
• $\hat S$2: linear singlet (tadpole) term.

The corresponding soft supersymmetry-breaking Lagrangian for the Higgs/singlet sector reads

$\hat S$3

Distinctive to the GNMSSM versus $\hat S$4-NMSSM are the explicit bilinear ($\hat S$5, $\hat S$6) and linear tadpole ($\hat S$7) terms. These parameters control the Higgsino and singlino mass independently and explicitly break the $\hat S$8 symmetry, resolving both tadpole and cosmological domain wall problems for the singlet.

2. Higgs and Neutralino Sectors

Higgs Sector

After electroweak symmetry breaking, the vacuum expectation values are $\hat S$9. The tree-level scalar potential combines $R$0-, $R$1-, and soft terms, ensuring vacuum stability for generic GNMSSM parameter choices (0910.1785, Cao et al., 2022, Cao et al., 2022). The CP-even mass matrix, in the $R$2 basis, is augmented relative to the MSSM via parameters $R$3.

Crucially, the presence of $R$4 and $R$5 enables decoupling of the Higgsino and singlino masses from the singlet scalar vev $R$6:

• The Higgsino mass: $R$7.
• The singlino mass: $R$8.

The flexibility in these parameters allows for larger singlet-doublet mixing in the CP-even sector without requiring large $R$9, increases control over the light singlet-like CP-even Higgs mass $CP$0, and facilitates the simultaneous realization of a SM-like $CP$1 at 125 GeV and a predominantly singlet $CP$2 at $CP$3 GeV (Cao et al., 2024, Cao et al., 2023).

Neutralino Sector

The neutralino mass matrix in the $CP$4 basis becomes

$CP$5

A singlino-dominated lightest neutralino ($CP$6) is achieved for $CP$7, with mixing controlled mainly by the $CP$8 ratio (Meng et al., 2024, Li et al., 9 Nov 2025, Cao et al., 2022). This singlet-dominance is the origin of the "secluded" dark sector phenomenology in the GNMSSM.

3. Solution to Cosmological and UV Problems

The explicit $CP$9-breaking terms in the superpotential and soft Lagrangian address two long-standing issues of the scale-invariant NMSSM:

• Domain wall problem: The accidental discrete symmetry leads to degenerate vacua and late-time domain walls, which are cosmologically problematic. The explicit breaking terms lift vacuum degeneracy, collapsing walls before nucleosynthesis (0910.1785, Cao et al., 2022).
• Tadpole problem: Planck-suppressed operators in supergravity can generate large singlet tadpoles, destabilizing the weak scale. The GNMSSM allows for appropriate tuning of tadpole and bilinear terms to avoid destabilization and maintain naturalness over a broad parameter region.

Unlike the $\mathbb{Z}_3$0-NMSSM, where the effective $\mathbb{Z}_3$1-term is $\mathbb{Z}_3$2, the GNMSSM's $\mathbb{Z}_3$3 and $\mathbb{Z}_3$4 ensure that neither fine-tuning nor cosmologically dangerous consequences are forced by discrete symmetries (0910.1785).

4. Collider and Low-Energy Phenomenology

Anomalies and Excesses

The GNMSSM provides unified explanations for:

• Muon anomalous magnetic moment ($\mathbb{Z}_3$5): Light electroweakinos and smuons, with $\mathbb{Z}_3$6 enhanced, yield $\mathbb{Z}_3$7 predominantly via wino–Higgsino–smuon (WHL) loops (Cao et al., 2024, Cao et al., 2022). Analytic expressions for all leading diagrams—including Bino–Higgsino–(L,R)-slepton and Bino-LR mixing contributions—are given by

$\mathbb{Z}_3$8

where $\mathbb{Z}_3$9 are loop functions, and relating parameters of the GNMSSM directly to the measured $g-2$0.

• Low-mass Higgs signals: Observed diphoton and $g-2$1 excesses near 95 GeV (LHC, LEP) are naturally interpreted as resonant production of the singlet-dominated CP-even Higgs $g-2$2. The couplings to SM states are suppressed but non-negligible due to controlled doublet admixture: $g-2$3 with $g-2$4 as singlet/doublet mixing. Required mixing to match observed strengths: $g-2$5, $g-2$6 (Cao et al., 2024, Cao et al., 2023).

Parameter Space and Experimental Constraints

Global parameter scans with flat priors over $g-2$7 show compatibility with:

• 125 GeV SM-like Higgs mass and couplings (HiggsBounds/HiggsSignals)
• Planck relic density, LZ spin-independent/direct detection bounds,
• B-physics ($g-2$8, $g-2$9),
• Vacuum stability and perturbative unitarity (Vevacious, SARAH),
• LHC SUSY and extra Higgs searches (CheckMATE, SModelS), requiring, for viable points:
  • $\hat S$0 GeV,
  • $\hat S$1 GeV,
  • $\hat S$2 GeV,
  • $\hat S$3 GeV,
  • $\hat S$4.

5. Dark Matter Phenomenology

The GNMSSM realizes a "secluded" WIMP dark matter scenario via a singlino-dominated $\hat S$5 annihilating into singlet-sector scalars: $\hat S$6 with $\hat S$7, $\hat S$8 singlet-dominated CP-even/odd Higgses. The annihilation cross sections are approximately

$\hat S$9

The relic density is achieved for $$W_{\rm GNMSSM} = W_{\rm Yukawa} + \lambda\,\hat S\,\hat H_u\cdot\hat H_d + \frac{\kappa}{3}\,\hat S^3 + \mu\,\hat H_u\cdot\hat H_d + \frac{\mu_S}{2}\hat S^2 + \xi_F\,\hat S\,,$$0–$$W_{\rm GNMSSM} = W_{\rm Yukawa} + \lambda\,\hat S\,\hat H_u\cdot\hat H_d + \frac{\kappa}{3}\,\hat S^3 + \mu\,\hat H_u\cdot\hat H_d + \frac{\mu_S}{2}\hat S^2 + \xi_F\,\hat S\,,$$1 at $$W_{\rm GNMSSM} = W_{\rm Yukawa} + \lambda\,\hat S\,\hat H_u\cdot\hat H_d + \frac{\kappa}{3}\,\hat S^3 + \mu\,\hat H_u\cdot\hat H_d + \frac{\mu_S}{2}\hat S^2 + \xi_F\,\hat S\,,$$2–$$W_{\rm GNMSSM} = W_{\rm Yukawa} + \lambda\,\hat S\,\hat H_u\cdot\hat H_d + \frac{\kappa}{3}\,\hat S^3 + \mu\,\hat H_u\cdot\hat H_d + \frac{\mu_S}{2}\hat S^2 + \xi_F\,\hat S\,,$$3 GeV (Meng et al., 2024, Li et al., 9 Nov 2025, Cao et al., 2022).

Direct detection cross sections scale as $$W_{\rm GNMSSM} = W_{\rm Yukawa} + \lambda\,\hat S\,\hat H_u\cdot\hat H_d + \frac{\kappa}{3}\,\hat S^3 + \mu\,\hat H_u\cdot\hat H_d + \frac{\mu_S}{2}\hat S^2 + \xi_F\,\hat S\,,$$4 (with moderate singlet-doublet Higgs mixing) and $$W_{\rm GNMSSM} = W_{\rm Yukawa} + \lambda\,\hat S\,\hat H_u\cdot\hat H_d + \frac{\kappa}{3}\,\hat S^3 + \mu\,\hat H_u\cdot\hat H_d + \frac{\mu_S}{2}\hat S^2 + \xi_F\,\hat S\,,$$5 if $$W_{\rm GNMSSM} = W_{\rm Yukawa} + \lambda\,\hat S\,\hat H_u\cdot\hat H_d + \frac{\kappa}{3}\,\hat S^3 + \mu\,\hat H_u\cdot\hat H_d + \frac{\mu_S}{2}\hat S^2 + \xi_F\,\hat S\,,$$6 is heavy, ensuring compliance with the LZ bound for $$W_{\rm GNMSSM} = W_{\rm Yukawa} + \lambda\,\hat S\,\hat H_u\cdot\hat H_d + \frac{\kappa}{3}\,\hat S^3 + \mu\,\hat H_u\cdot\hat H_d + \frac{\mu_S}{2}\hat S^2 + \xi_F\,\hat S\,,$$7–$$W_{\rm GNMSSM} = W_{\rm Yukawa} + \lambda\,\hat S\,\hat H_u\cdot\hat H_d + \frac{\kappa}{3}\,\hat S^3 + \mu\,\hat H_u\cdot\hat H_d + \frac{\mu_S}{2}\hat S^2 + \xi_F\,\hat S\,,$$8.

Nested-sampling and Bayesian analyses favor a singlino-dominated LSP in $$W_{\rm GNMSSM} = W_{\rm Yukawa} + \lambda\,\hat S\,\hat H_u\cdot\hat H_d + \frac{\kappa}{3}\,\hat S^3 + \mu\,\hat H_u\cdot\hat H_d + \frac{\mu_S}{2}\hat S^2 + \xi_F\,\hat S\,,$$9–$W_{\rm Yukawa}$0 of the parameter space, with annihilation typically dominated by $W_{\rm Yukawa}$1 ($W_{\rm Yukawa}$2), $W_{\rm Yukawa}$3 ($W_{\rm Yukawa}$4), and $W_{\rm Yukawa}$5 ($W_{\rm Yukawa}$6) (Meng et al., 2024, Li et al., 9 Nov 2025, Cao et al., 2022).

Characteristic mass hierarchies:

• $W_{\rm Yukawa}$7: light Bino, $W_{\rm Yukawa}$8 GeV, annihilation via $W_{\rm Yukawa}$9 or $\lambda,\kappa$0, mild tuning.
• $\lambda,\kappa$1: heavy Bino, $\lambda,\kappa$2 GeV, requires small $\lambda,\kappa$3 for direct detection, large tuning.

Benchmarks in the literature exemplify points yielding correct $\lambda,\kappa$4, $\lambda,\kappa$5, $\lambda,\kappa$6, $\lambda,\kappa$7, and direct detection rates, for both Bino- and singlino-dominated scenarios (Cao et al., 2024, Cao et al., 2023).

6. Experimental and Future Tests

A comprehensive suite of collider, dark matter, and low-energy measurements constrain the GNMSSM, but large portions of parameter space remain viable:

• High-Luminosity LHC ($\lambda,\kappa$8): Can probe compressed electroweakino spectra ($\lambda,\kappa$9 down to $\hat S$00 GeV) and direct slepton production ($\hat S$01 TeV) (Cao et al., 2024).
• Future $\hat S$02 colliders (ILC, CLIC, FCC-ee): Expected sensitivity to $\hat S$03 couplings at a few percent and improved $\hat S$04 GeV Higgs mass resolution, facilitating precision studies of singlet-like Higgs states.
• Direct detection: LZ 2024 and future multi-ton experiments will test $\hat S$05 down to $\hat S$06 cm$\hat S$07; future improvements by a factor of 5 would strongly impact the allowed parameter space (Meng et al., 2024, Li et al., 9 Nov 2025).
• Muon $\hat S$08 (FNAL/J-PARC): Ongoing improvements will further challenge or confirm the surviving corners of GNMSSM parameter space.
• Higgs property measurements: Precision determinations of the 125 GeV Higgs couplings to the percent level will critically test the singlet-doublet mixing structure required for low-mass excesses (Cao et al., 2023).
• Dedicated LHC searches: Targeted analyses for extended decay chains with soft leptons and multi-step cascades will be essential for probing the fully-realized GNMSSM scenario (Li et al., 9 Nov 2025, Cao et al., 2022).

The broad decoupling and flexible parameter structure of the GNMSSM ensure its continued empirical testability and its capacity to synthesize diverse anomalies within a natural, UV-complete, and cosmologically-viable supersymmetric framework.