Interpreting Light Scalar Excesses and Heavy Scalar Cascades in the $\mu$-Term Extended NMSSM

This paper demonstrates that the $\mu$-term extended NMSSM can simultaneously accommodate persistent $95~{\rm GeV}$ scalar excesses at LEP and the LHC while satisfying stringent heavy scalar cascade constraints, offering viable parameter regions with distinct coupling patterns and predicting specific long-lived neutralino signatures for positive-$\mu$ scenarios.

The Gist

Imagine the universe is like a giant, complex orchestra. For decades, physicists have been listening to the "Standard Model," which is like a sheet of music that predicts exactly how particles should behave. Most of the time, the music plays perfectly. But recently, the orchestra has started playing a few strange, slightly off-key notes that the sheet music doesn't explain.

This paper is like a detective story where the author, Jingwei Lian, tries to find a new instrument in the orchestra that could explain these weird notes. The instrument is a specific version of a theory called Supersymmetry, specifically a "tuned-up" version called the µNMSSM.

Here is the breakdown of the mystery and the solution, using simple analogies:

The Mystery: Two Strange Noises

The scientists have heard two distinct "glitches" in the data from giant particle colliders (LEP and LHC):

1. The Light Whisper (95 GeV): There is a faint signal of a light, new particle appearing around 95 GeV (a unit of mass). It's showing up in two ways:

   • It's turning into pairs of bottom quarks (heavy particles).
   • It's turning into pairs of light particles (photons).
   • The Problem: The Standard Model says this shouldn't happen, or at least not this strongly.

2. The Heavy Cascade (600–650 GeV): There are hints of a much heavier particle that decays (falls apart) into the famous 125 GeV Higgs boson (the one we already know) plus one of those new light particles.

   • The Twist: Recent searches have gotten stricter. The "heavy" particle isn't showing up exactly where the first hints suggested (650 GeV), but the data is still a bit fuzzy, with some noise appearing around 600 GeV.

The Solution: A New Musical Theory

The author suggests that the "µNMSSM" theory is the right sheet of music to explain these noises. Think of this theory as a house with more rooms than the Standard Model.

• The Standard Model has one main room (the Higgs doublet).
• The µNMSSM adds a secret, hidden room (a "singlet" scalar).

The author argues that the "Light Whisper" (95 GeV) is actually a guest from this hidden room. Because it's mostly a "singlet" (a particle that doesn't interact much with the usual forces), it can hide easily but still leak out enough to be seen as those strange signals.

The Two Ways the Mystery is Solved

The paper finds that this theory works in two distinct "patterns" or "styles," like two different ways to solve a puzzle:

• Pattern 1: The "Quiet" Whisper.
  In this version, the light particle is very shy. It barely talks to bottom quarks. Because it doesn't talk to them much, it doesn't decay into them often. Instead, it turns into photons (light) more frequently. This fits the LHC photon data perfectly but explains the older LEP data poorly.

  • Analogy: Imagine a shy singer who refuses to sing the heavy bass notes (bottom quarks) but is great at singing high-pitched notes (photons).

• Pattern 2: The "Loud" Whisper.
  In this version, the light particle is a bit more social. It mixes with the known particles and talks to bottom quarks more often. This fits the older LEP data well but makes the photon signal weaker.

  • Analogy: This singer loves the bass notes but gets a bit hoarse when trying to hit the high notes.

The author shows that both patterns are mathematically possible and fit the current data within a reasonable margin of error.

The Heavy Cascade: A Domino Effect

The paper also looks at the "Heavy" particle. It suggests that this heavy particle acts like a domino. It falls apart into the known 125 GeV Higgs and the new 95 GeV particle.

• The paper predicts that while the "perfect" signal seen in old data might be gone, there is still a "faint echo" of this cascade happening.
• It also suggests a new type of domino effect involving "CP-odd" particles (a different kind of particle spin) that could explain a different type of noise seen around 600 GeV and 400 GeV in recent data.

The "Gravitino" Twist (The Ghost in the Machine)

There is a special subset of these solutions where the math requires a specific sign for a parameter (called "positive µ").

• If this condition is met, the theory predicts that the universe is filled with Gravitinos (ghostly, ultra-light particles) as the main form of Dark Matter.
• In this scenario, the lightest neutral particle (the Neutralino) isn't the Dark Matter; instead, it's a "middleman" that eventually decays into the Gravitino.
• The Catch: Because this decay takes a long time, these particles might travel a few meters or even kilometers inside the detector before disappearing. This makes them very hard to catch with current "instant" detectors, but future experiments looking for "delayed" signals could find them.

The Bottom Line

The author concludes that this specific version of Supersymmetry (µNMSSM) is a viable candidate to explain the strange 95 GeV signals and the heavy particle searches.

• It successfully fits the data without breaking other known laws of physics.
• It predicts that future experiments should look for specific patterns: a mix of light particles turning into photons and bottom quarks, and heavy particles falling apart in specific ways.
• It also suggests that if we look for "delayed" signals (particles that hang around a bit before vanishing), we might find evidence of a Gravitino-based Dark Matter.

In short, the paper says: "The Standard Model is missing a few notes. Our theory adds a hidden room and a secret instrument that explains the noise, and we have two different ways the music could be playing."

Technical Summary

Technical Summary: Interpreting Light Scalar Excesses and Heavy Scalar Cascades in the µ-Term Extended NMSSM

Problem Statement
The paper addresses two persistent, yet conflicting, experimental anomalies in the Higgs sector. First, there are hints of a light scalar resonance near 95 GeV, observed as a mild excess in the LEP $b\bar{b}$ channel ($\mu_{b\bar{b}} \approx 0.117$) and a more significant excess in the LHC diphoton channel ($\mu_{\gamma\gamma} \approx 0.24$). Second, searches for a heavier scalar resonance decaying into the 125 GeV Higgs and an additional scalar have evolved. While an earlier CMS analysis suggested an excess near $(M_X, M_Y) \approx (650, 95)$ GeV in the $\gamma\gamma b\bar{b}$ channel, subsequent analyses (including the latest CMS $b\bar{b}b\bar{b}$ search) have not confirmed this specific point, instead showing the largest local deviations near $(600, 400)$ GeV. The challenge is to construct a Beyond the Standard Model (BSM) framework that can simultaneously accommodate the 95 GeV light scalar signals and the emerging heavy scalar structures without violating current Higgs, flavor, and direct SUSY search limits.

Methodology
The authors investigate these anomalies within the $\mu$-term extended Next-to-Minimal Supersymmetric Standard Model ($\mu$NMSSM). This model is a restricted version of the General NMSSM (GNMSSM), retaining an explicit MSSM-like $\mu$ term while omitting singlet bilinear and tadpole terms. This setup is motivated by inflation-inspired supergravity constructions where the $\mu$ term is induced by non-minimal Higgs-gravity coupling.

The analysis proceeds via a comprehensive parameter space scan using the MultiNest algorithm. The scan focuses on nine key parameters: the singlet-doublet couplings $\lambda$ and $\kappa$, the mixing parameter $\delta$, the doublet and singlet mass scales $m_A$ and $m_B$, $\tan\beta$, and the Higgsino mass parameters $\mu_{tot}$ and $\mu_{eff}$. Other SUSY parameters (e.g., gluino mass fixed at 3 TeV) are set to decouple heavy colored states.

The likelihood function combines:

1. Anomaly Fit: A $\chi^2$ term targeting the central values of the 95 GeV diphoton and $b\bar{b}$ excesses, as well as the earlier 650 GeV $\gamma\gamma b\bar{b}$ excess.
2. Constraints: Hard cuts and Gaussian penalties for:
   • The 125 GeV Higgs signal rates (using HiggsSignals).
   • Exclusion limits on extra Higgs bosons from LEP, Tevatron, and LHC (using HiggsBounds).
   • Flavor observables ($B_s \to \mu^+\mu^-$, $B \to X_s\gamma$).
   • Electroweak precision observables ($S, T, U$).
   • Direct SUSY searches (using SModelS).
   • Vacuum stability checks (using Vevacious) applied as a posterior diagnostic.

Key Contributions and Results
The study identifies viable regions in the $\mu$NMSSM parameter space that accommodate the 95 GeV excesses within $2\sigma$. The viable points naturally split into two distinct phenomenological scenarios based on the reduced couplings of the light scalar ($h_s$):

1. Scenario I (Diphoton Favored):

   • Characterized by $\mu_{\gamma\gamma} > 0.11$.
   • The light scalar $h_s$ is predominantly singlet-like with a suppressed $h_s b\bar{b}$ coupling.
   • This suppression reduces the total width of $h_s$, enhancing the branching ratio to photons ($Br(h_s \to \gamma\gamma)$), thereby fitting the LHC diphoton excess.
   • Consequently, the LEP $b\bar{b}$ rate is small ($\mu_{b\bar{b}} \sim \mathcal{O}(10^{-2})$), consistent with the lack of a strong LEP signal in this channel.
   • This scenario requires lighter Higgsino-like charginos to enhance the loop-induced $h_s \gamma\gamma$ coupling.

2. Scenario II (LEP Favored):

   • Characterized by $\mu_{\gamma\gamma} < 0.11$.
   • The light scalar $h_s$ has a larger doublet component, leading to stronger couplings to the $Z$ boson and $b$-quarks.
   • This enhances the LEP $e^+e^- \to Z h_s$ rate, fitting the LEP excess ($\mu_{b\bar{b}} \approx 0.15-0.20$).
   • However, the increased $h_s \to b\bar{b}$ partial width suppresses the diphoton branching fraction, resulting in a lower LHC diphoton signal.

Heavy Scalar Cascades:

• CP-Even Sector: The heavy doublet-like scalar $H$ (mass $\sim 625-670$ GeV) decays via cascades $H \to h h_s$ and $H \to h_s h_s$. The predicted rate for the $\gamma\gamma b\bar{b}$ final state is generally below the original CMS best-fit value (typically $10^{-2}$ to $10^{-1}$ fb) but remains within reach of future searches. The model also predicts observable rates for $H \to h_s h_s \to b\bar{b}b\bar{b}$ and $H \to h_s h_s \to \gamma\gamma b\bar{b}$.
• CP-Odd Sector: The paper highlights a new topology involving the CP-odd sector: $A_2 \to h A_1 \to 4b$. This cascade naturally populates the mass region $(M_X, M_Y) \approx (600, 400)$ GeV, where the latest CMS $4b$ analysis observes its largest local deviation. The predicted rates are below current limits but offer a specific target for dedicated fits.

Gravitino Dark Matter (Positive-$\mu$ Subset):
For the subset of points with a positive explicit $\mu$ term (motivated by inflationary supergravity), the model admits a gravitino as the Lightest Supersymmetric Particle (LSP).

• The gravitino mass ranges from the keV scale to tens of MeV.
• The lightest neutralino ($\tilde{\chi}^0_1$) becomes the Next-to-Lightest Supersymmetric Particle (NLSP).
• Due to the small decay widths into $\gamma \tilde{G}$, $Z \tilde{G}$, and $h \tilde{G}$, the neutralino NLSP is typically long-lived at collider scales (decaying outside the main detectors or with displaced vertices).
• This long lifetime renders standard prompt SUSY searches ineffective, while the gravitino relic abundance is determined by the reheating temperature ($T_R$).

Significance and Claims
The paper claims that the $\mu$NMSSM provides a viable and testable framework for the 95 GeV anomalies, despite being more constrained than the general GNMSSM.

• It demonstrates that the model can reconcile the light scalar hints with the heavy scalar search landscape, including the shift from the $(650, 95)$ GeV to the $(600, 400)$ GeV region.
• It identifies two distinct mechanisms (suppressed $b\bar{b}$ coupling vs. enhanced doublet mixing) to explain the tension between LEP and LHC light scalar signals.
• It connects the collider phenomenology to a specific cosmological scenario (inflation-induced $\mu$ term) where the gravitino is the LSP and the neutralino is a long-lived NLSP.
• The authors modestly conclude that while the model is consistent with current data, the preferred regions are tightly constrained. They emphasize that sharper statements regarding the viability of these scenarios require dedicated future analyses of $4b$, $\gamma\gamma b\bar{b}$, $b\bar{b}\tau\bar{\tau}$, and long-lived particle signatures. Vacuum stability checks indicate that some benchmark points may be short-lived, suggesting that strict theoretical consistency could further reduce the viable parameter space.