Kinematic Riffs and Interference Effects in Triple Higgs Production in the N2HDM

This paper investigates resonant triple Higgs production within the Next-to-minimal Two Higgs Doublet Model (N2HDM), demonstrating that interference effects and additional decay channels significantly alter kinematic distributions and necessitate fully differential studies over simplified approximations to accurately probe extended Higgs sectors at the LHC.

The Gist

Imagine the Large Hadron Collider (LHC) as a giant, high-speed particle smasher. Its main job is to crash protons together to see what tiny pieces fly out. For a long time, scientists have been looking at the "Higgs boson," a particle that gives other particles mass. Usually, they look for just one Higgs boson appearing after a crash. But now, they are trying to catch a much rarer event: three Higgs bosons appearing at once.

This paper is like a detailed investigation into what happens when we try to catch these "triple Higgs" events, specifically looking at a theoretical model called the N2HDM (Next-to-minimal Two Higgs Doublet Model). Think of this model as a slightly more complex version of the standard rules of physics, where there are extra, heavier "sibling" Higgs particles hiding in the mix.

Here is the breakdown of their findings using simple analogies:

1. The "Double-Resonance" Shortcut vs. The Full Reality

In the past, scientists often tried to understand these complex crashes by looking for a specific, simple pattern. They imagined a "domino effect":

• A heavy particle (let's call it H3) is created.
• It instantly breaks into a medium-heavy particle (H2) and a normal Higgs.
• The medium-heavy particle (H2) then instantly breaks into two more normal Higgses.

This is called the "Double-Resonance" scenario. It's like watching a magician pull a rabbit out of a hat, and then the rabbit pulls two more rabbits out of its own hat. It's a clean, easy-to-follow story.

The Paper's Discovery: The authors found that relying only on this simple "domino" story is dangerous. While it happens, it's not the whole story. The real crash is like a chaotic traffic jam where cars (particles) are swerving, merging, and crashing into each other in ways that don't follow a straight line.

2. The "Interference" Effect (The Noise in the Signal)

The most important finding in this paper is about interference. In physics, when different ways of creating the same result happen at the same time, they can either boost each other up or cancel each other out.

• The Analogy: Imagine two people singing the same note. If they sing in perfect sync, the sound gets louder (constructive interference). If one sings slightly out of phase, they might cancel each other out, and you hear silence (destructive interference).
• The Result: The authors found that in these triple Higgs crashes, the "simple domino" path often gets cancelled out by other messy paths happening at the same time. Sometimes, the messy paths cancel the simple path so much that the total number of events is actually lower than if you just looked at the simple path alone.

This means that if you only look for the "clean domino" pattern, you might miss the event entirely, or you might think you see more events than are actually there.

3. Why Mass Matters (The Weight of the Particles)

The paper tested different "weights" (masses) for these heavy sibling particles.

• Lighter Weights: When the heavy particles are just heavy enough to break apart into the lighter ones, the "domino" story works pretty well. It's like a heavy box that easily splits into two smaller boxes.
• Heavier Weights: When the particles get much heavier, the "domino" story falls apart. The particles can break apart in many different, messy ways at once. The paper shows that even if the "domino" path is the most common single path, the messy, non-domino paths are still doing a lot of work, changing the shape of the data.

4. The "Fingerprint" of the Crash

How do scientists tell the difference between the simple story and the messy reality? The paper suggests looking at specific "fingerprints" left behind in the data:

• Invariant Mass: This is like weighing the total debris from the crash. The simple story predicts specific weights (peaks) where the debris should pile up. The messy reality shows extra piles of debris in unexpected places.
• Transverse Momentum ($p_T$): This is like measuring how hard the debris is flying sideways. The simple story predicts the debris flies in a certain way. The messy reality shows the debris flying much harder or softer than expected, creating a "tail" in the data that the simple story can't explain.

The Bottom Line

The main message of this paper is a warning to physicists: Don't oversimplify.

If you try to understand the complex world of triple Higgs production by only looking for the clean, step-by-step "domino" effect, you will get the wrong answer. The real world is full of "interference" and messy, off-path events that change the numbers and the shapes of the data.

To truly understand what's happening in the universe (and to find new physics beyond our current understanding), scientists need to look at the entire chaotic picture, not just the clean parts. They need to account for all the swerving, cancelling, and messy interactions, or they might miss the discovery entirely.

Technical Summary

Technical Summary: Kinematic Riffs and Interference Effects in Triple Higgs Production in the N2HDM

Problem Statement
The production of three Higgs bosons ($gg \to hhh$) at the Large Hadron Collider (LHC) represents one of the most complex final states in hadron collider physics. In the Standard Model (SM), this process is heavily suppressed by small Higgs self-couplings and loop-induced nature, rendering it effectively unobservable. However, in models with extended scalar sectors, such as the Next-to-minimal Two-Higgs-Doublet Model (N2HDM), production rates can be significantly enhanced.

A common strategy in theoretical and experimental analyses involves simplifying the process by focusing on "resonant topologies," specifically double-resonant (DR) contributions where intermediate heavy scalars decay on-shell (e.g., $gg \to H_3 \to H_2 h \to hhh$). While this factorized approach offers a tractable framework, the authors posit that it may fail to capture the full complexity of the underlying dynamics. Specifically, the interplay between resonant and non-resonant diagrams, off-shell effects, and quantum interference can significantly alter kinematic observables and total rates. The paper investigates whether simplified DR approximations are sufficient to describe triple Higgs production within the rich scalar spectra of the N2HDM, which features three CP-even neutral scalars ($H_1, H_2, H_3$).

Methodology
The study employs the N2HDM as a representative framework for an extended scalar sector with three neutral CP-even states. The methodology proceeds through the following steps:

1. Parameter Scanning: Using the ScannerS tool, the authors scan the N2HDM parameter space, constrained by theoretical consistency (perturbative unitarity, vacuum stability) and experimental bounds (Higgs signal strengths, direct searches, electroweak precision data) via HiggsTools.
2. Benchmark Selection: Specific mass configurations for the heavier scalars ($m_{H_2}, m_{H_3}$) are selected to maximize the double-resonant cross-section ($\sigma_{DR}$). These configurations range from near-threshold regions ($m_{H_2} \sim 2m_h, m_{H_3} \sim 3m_h$) to heavier mass hierarchies where additional decay channels (e.g., $H_2 \to hhh$) become kinematically accessible.
3. Event Generation: For selected benchmark points, events are generated using MadGraph5_aMC@NLO at the one-loop level in QCD. To isolate specific contributions, the authors manipulate coupling orders and manually set specific diagrams to zero. The generated samples include:
   • Full Process: All relevant topologies.
   • Double-Resonant (DR): Only the sequential on-shell decay chain ($H_3 \to H_2 h \to hhh$).
   • Non-DR: All processes excluding the DR topology.
   • Box Diagrams: Specific contributions involving box diagrams with $H_2$ or $H_3$ propagators.
4. Kinematic Analysis: The study analyzes differential distributions for the invariant masses of the di-Higgs ($M_{hh}$) and tri-Higgs ($M_{hhh}$) systems, as well as the transverse momentum ($p_T$) of the Higgs bosons, comparing the full process against the DR approximation.

Key Contributions and Results

• Breakdown of Factorization: The analysis demonstrates that even in regions where the double-resonant cross-section is maximized and numerically close to the total cross-section, the DR approximation often fails to reproduce the full kinematic distributions. The ratio $\sigma_{DR}/\sigma_{full}$ is shown to be an indicative but insufficient measure of the DR contribution's dominance due to interference effects.
• Significant Interference Effects: Destructive interference between resonant and non-resonant amplitudes is observed across all benchmark points. In several cases, this interference reduces the total cross-section below the DR contribution alone, leading to scenarios where $\sigma_{DR} > \sigma_{full}$. This highlights that the total rate cannot be understood as a simple sum of independent topologies.
• Kinematic Discrepancies:
  • Invariant Masses: While DR contributions produce clear peaks at $m_{H_2}$ and $m_{H_3}$, the full process exhibits significant tails and additional structures arising from off-shell decays and continuum contributions. For instance, in heavier mass configurations, the full process shows peaks corresponding to three-body decays ($H_2 \to hhh$) and box diagrams that are absent in the pure DR approximation.
  • Transverse Momentum ($p_T$): The $p_T$ distributions reveal that the DR approximation falls off rapidly at high momenta, whereas the full process maintains a longer tail due to continuum and non-resonant contributions.
• Mass Hierarchy Dependence: The study explores various mass hierarchies. In scenarios where $m_{H_2} > 3m_h$, the opening of the $H_2 \to hhh$ channel introduces additional resonant structures and interference patterns that further complicate the separation of contributions.
• Angular Variables: The authors note that pseudorapidity ($\eta$) and azimuthal angle ($\phi$) distributions show negligible structural differences between the full process and the DR approximation, suggesting these variables are less effective for disentangling the underlying production mechanisms compared to invariant mass and $p_T$.

Significance and Claims
The paper argues that simplified descriptions based solely on factorized double-resonant topologies are generally insufficient to capture the dynamics of $gg \to hhh$ in extended scalar sectors. The authors claim that:

1. Interference is Ubiquitous: Interference effects and off-shell contributions play a critical role in shaping both total rates and differential distributions, even in regimes where resonant enhancements are kinematically favored.
2. Need for Full Differential Studies: Relying on simplified approximations can lead to misleading conclusions regarding signal rates and kinematic shapes. A faithful interpretation of potential triple Higgs signals requires fully differential studies that account for the complete kinematic structure, including interference.
3. Experimental Implications: The distinct features in invariant mass and transverse momentum distributions suggest that these observables should be explicitly included in multivariate analysis strategies (e.g., boosted decision trees, graph neural networks) to effectively disentangle the underlying production mechanisms.
4. Model Generality: While the analysis is specific to the N2HDM, the authors assert that the observed features—correlations between masses and couplings, multiple interfering topologies, and deviations from sequential on-shell decay descriptions—are expected to be generic to renormalizable scalar extensions of the SM with three or more degrees of freedom (e.g., C2HDM, Georgi-Machacek model).

In conclusion, the work underscores the necessity of moving beyond simplified resonant pictures to fully exploit the potential of triple Higgs production as a probe of the scalar sector and Beyond the Standard Model physics at the HL-LHC.