⚛️ phenomenology

Resonant Di-Higgs Searches in $b\bar{b}\tau^+\tau^-$ at HL-LHC: Supersymmetry versus Compositeness Benchmarks

This paper investigates the High-Luminosity LHC's potential to discover resonant di-Higgs production in the $b\bar{b}\tau^+\tau^-$ channel, finding that while Supersymmetry benchmarks face significant background challenges, Compositeness models offer substantial discovery potential in a background-free environment.

Original authors: Atri Dey, Carl Johan Konigsson, Stefano Moretti, Luigi Delle Rose, Luca Panizzi, Stefania De Curtis

Published 2026-09-18

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: Atri Dey, Carl Johan Konigsson, Stefano Moretti, Luigi Delle Rose, Luca Panizzi, Stefania De Curtis

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

The universe is held together by a delicate balance of forces, and at the heart of this balance sits a tiny, invisible particle known as the Higgs boson. Discovered in 2012, this particle is the physical proof of a field that gives mass to other fundamental particles, acting like a cosmic molasses that slows them down. While the discovery was a triumph, it also left physicists with a nagging question: is the Higgs boson a solitary figure, or is it part of a larger family? The standard model of physics, our best current map of reality, suggests there is only one. However, many theories propose that there are heavier, hidden cousins waiting to be found. If these heavier versions exist, they might occasionally decay into pairs of the familiar Higgs bosons, creating a unique signature that could reveal the deeper structure of the universe.

To find these hidden signals, scientists are looking toward the High-Luminosity Large Hadron Collider, a massive machine in Switzerland designed to smash protons together with unprecedented energy and frequency. In a recent study, researchers simulated what would happen if these collisions produced pairs of Higgs bosons that then decayed into specific, detectable particles: two bottom quarks and two tau leptons. They focused on two leading theories that attempt to solve the biggest mysteries in physics: one based on supersymmetry, which suggests every known particle has a heavier partner, and another based on compositeness, which proposes that the Higgs boson is actually a composite object made of even smaller, tightly bound pieces. By running detailed computer simulations of billions of collisions, the team investigated whether the upcoming collider could distinguish between these two theories and, more importantly, whether it could actually see the new particles they predict.

The researchers set up a virtual experiment to mimic the conditions of the future collider, generating millions of simulated collision events. They looked for a specific pattern where a heavy, unseen particle would briefly appear and then split into two lighter Higgs bosons. These Higgs bosons would then immediately break apart into the two bottom quarks and two tau leptons the team was tracking. The challenge is immense because the background noise from ordinary particle interactions is overwhelming. It is like trying to hear a single violin in a stadium full of cheering fans. The team had to design a sophisticated filter to separate the rare, interesting events from the sea of common ones. They focused on the energy and direction of the particles, looking for signs that the Higgs bosons were moving very fast, a clue that they had been born from the decay of a much heavier parent particle.

When they applied their filters to the simulated data, the results for the two theories diverged sharply. For the supersymmetry scenario, the signal was faint and difficult to separate from the background noise. Even with the massive amount of data expected from the future collider, the researchers found that the difference between the new physics and the standard background was too subtle to be certain. The heavy particles predicted by this theory were not heavy enough to create a distinct enough signature to stand out clearly against the noise. The study suggests that while the experiment might see a hint of something unusual, it would likely fall just short of a definitive discovery, leaving the question of supersymmetry in this specific context unresolved.

In stark contrast, the compositeness scenario painted a much brighter picture. In this model, the heavy particles involved are significantly more massive, and the physics governing their creation creates a much cleaner signal. When the researchers applied the same filters to this scenario, the background noise vanished almost entirely, leaving a clear, distinct group of events that could not be explained by the standard model. The simulation showed that the future collider would have a very strong chance of discovering these new particles, with the signal rising far above the threshold needed to claim a discovery. This difference arises because the heavy partners in the compositeness theory are so massive that they impart a much stronger "kick" to the resulting particles, making them easier to spot against the backdrop of ordinary collisions.

The study concludes that the path to understanding the Higgs boson's true nature depends heavily on which of these theories describes reality. If nature follows the compositeness model, the upcoming collider will likely provide a clear window into a new world of physics, revealing a heavy particle that decays into pairs of Higgs bosons with high confidence. If nature follows the supersymmetry model, the same experiment will likely struggle to find a clear answer, as the signal remains too buried in the noise to be certain. The researchers emphasize that their work provides a roadmap for experimentalists, showing exactly which types of data to look for and how to analyze it. They have demonstrated that the specific combination of particles they studied offers a powerful way to test these theories, but the ultimate success will depend on the actual data the machine produces and the ability to control the uncertainties in the background noise.

Technical Summary: Resonant Di-Higgs Searches in $b\bar{b}\tau^+\tau^-$ at HL-LHC

Problem Statement
The discovery of the 125 GeV Higgs boson confirmed the mechanism of Electroweak Symmetry Breaking (EWSB) within the Standard Model (SM). However, the Higgs sector remains poorly constrained, and open questions regarding the hierarchy problem, dark matter, and baryon asymmetry suggest the existence of Beyond the SM (BSM) physics. A primary probe for extended Higgs sectors is SM di-Higgs ($hh$) production, which provides direct access to the Higgs trilinear self-coupling. In the SM, this process is highly suppressed due to destructive interference between triangle and box diagrams in gluon-gluon fusion (ggF). Consequently, observing SM di-Higgs production is challenging even at the High-Luminosity LHC (HL-LHC).

Many BSM scenarios, specifically Supersymmetry (SUSY) and Composite Higgs Models (CHM), predict resonant di-Higgs production via a heavier scalar state ( $H$ ) decaying into a pair of SM-like Higgs bosons ( $gg \to H \to hh$ ). This paper investigates the discovery potential of the resonant process $gg \to H \to hh \to b\bar{b}\tau^+\tau^-$ at the HL-LHC, comparing two specific theoretical frameworks: the Next-to-Minimal Supersymmetric SM (NMSSM) and the Composite 2-Higgs Doublet Model (C2HDM). The study aims to determine if characteristic differences between these models survive realistic detector reconstruction and to assess the sensitivity of the $2b2\tau$ final state, which offers a balance between large branching ratios and manageable backgrounds.

Methodology
The authors perform a comprehensive collider analysis using Monte Carlo (MC) simulations at the detector level, assuming HL-LHC operating conditions ( $\sqrt{s} = 14$ TeV, integrated luminosity $L = 4000$ fb $^{-1}$ ).

Theoretical Frameworks and Benchmarks:
- NMSSM: A benchmark point (BP) is selected featuring a heavy doublet Higgs ( $m_H \approx 809$ GeV) and a singlet ( $m_S \approx 1193$ GeV). The model includes contributions from top squarks (stops) in the loop, leading to constructive interference in specific mass intervals.
- C2HDM: A representative BP is chosen with a heavy Higgs mass ( $m_H \approx 1182$ GeV) and eight Heavy Top Partners (HTPs). The model features effective vertices involving HTPs that modify the production amplitude and kinematic spectrum.
- Both benchmarks are tuned to yield an inclusive di-Higgs cross section approximately 50% higher than the SM irreducible background, despite having different mass spectra and underlying dynamics.
Simulation and Reconstruction:
- Signal and background events are generated using MadGraph5_aMC@NLO, showered and hadronized with Pythia8, and passed through Delphes3 with an HL-LHC detector configuration.
- The analysis focuses on the $b\bar{b}\tau^+\tau^-$ final state, requiring at least two $b$ -tagged jets and two hadronically decaying $\tau$ -tagged jets.
- Backgrounds: The dominant reducible backgrounds are $t\bar{t}$ (both fully leptonic and semileptonic) and $Z$ +jets. QCD multi-jet production is also considered, though its impact is mitigated by high-mass cuts. The irreducible background is SM di-Higgs production.
Analysis Strategy:
- The study utilizes a suite of kinematic observables to distinguish signal from background and to differentiate between the two BSM scenarios. Key variables include:
  - Invariant masses: $M(bb)$, $M(\tau\tau)$ , and the four-body mass $M(bb\tau\tau)$ .
  - Transverse variables: Total hadronic activity ( $H_T$ ) and leading object transverse momenta ( $p_T$ ).
  - Angular separations: $\Delta R(bb)$ and $\Delta R(\tau\tau)$ to probe boosted topologies.
  - Polarization-sensitive observable: The leading pion momentum fraction ( $Z_g = p_T(\pi_1)/p_T(\tau_1)$ ) within the $\tau$ -jet, which exploits the different helicity structures of $\tau$ leptons from scalar Higgs decays versus $W$ boson decays.
- A sequential cut-based analysis is applied to optimize the signal-to-background ratio.

Key Contributions

Comparative Study: The paper provides a direct comparison of resonant di-Higgs production in NMSSM and C2HDM under identical experimental conditions, highlighting how distinct theoretical mechanisms (stops vs. HTPs) manifest in detector-level observables.
Polarization Sensitivity: It demonstrates the utility of the $Z_g$ observable in suppressing the dominant $t\bar{t}$ background by exploiting the helicity differences between signal and background $\tau$ leptons.
Realistic Detector Analysis: Unlike parton-level studies, this work incorporates full detector effects, jet reconstruction uncertainties, and neutrino smearing, providing a realistic assessment of the resonance reconstruction capabilities at the HL-LHC.

Results
The analysis yields significantly different sensitivities for the two models:

NMSSM Case:
- The resonant signal exhibits a moderate enhancement over the SM background.
- While the cut strategy effectively suppresses reducible backgrounds ( $t\bar{t}$ , $Z$ +jets), the residual dileptonic $t\bar{t}$ background and the irreducible SM di-Higgs contribution remain substantial.
- The final statistical significance for the NMSSM signal (combined ATLAS+CMS) reaches approximately 1.3 $\sigma$ .
- The sensitivity is heavily limited by systematic uncertainties; even with 1% background uncertainty, the significance does not reach the discovery threshold. The kinematic differences between the NMSSM signal and the SM background are not sufficiently distinct to allow for clear discrimination after detector effects.
C2HDM Case:
- Due to the higher mass scale of the heavy scalar and the HTPs, the C2HDM signal exhibits a much harder kinematic spectrum (higher $H_T$ , $p_T$ , and $M(bb\tau\tau)$ ) compared to both the SM and the NMSSM.
- The optimized selection strategy, particularly the high-mass cut ( $M(bb\tau\tau) > 500$ GeV) and the polarization cut ( $Z_g$ ), successfully eliminates all simulated reducible backgrounds.
- The analysis results in a "background-free" environment (within MC statistics) where only the irreducible SM di-Higgs and the C2HDM signal remain.
- The final statistical significance for the C2HDM excess over the SM is 4.12 $\sigma$ (single experiment) and 5.83 $\sigma$ (combined ATLAS+CMS), well within the discovery range.
- This result remains robust even under pessimistic assumptions of 50% background systematic uncertainty.

Significance and Claims
The paper concludes that the $b\bar{b}\tau^+\tau^-$ channel is a promising probe for BSM physics at the HL-LHC, but its effectiveness is highly model-dependent.

Modest Claims for NMSSM: The authors state that for the specific NMSSM benchmark chosen, the discovery potential is limited. The kinematic features (peaks, thresholds) visible at the parton level are washed out by detector effects, making it difficult to distinguish the NMSSM signal from the SM continuum or to map the signal back to the specific underlying theory.
Strong Claims for C2HDM: In contrast, the C2HDM scenario offers substantial scope for discovery. The distinct kinematic hardness induced by the heavy top partners allows for efficient background rejection, leading to a clear excess over the SM prediction.
Methodological Value: The authors advocate for the specific selection strategies developed (utilizing $Z_g$ and invariant mass ratios) as robust tools for experimental collaborations. They emphasize that while the $2b2\tau$ channel can probe BSM scenarios beyond the electroweak scale, the ability to disentangle the specific BSM mechanism (SUSY vs. Compositeness) relies heavily on the mass scale and the resulting kinematic hardness of the signal, which may be obscured in lighter scenarios like the NMSSM benchmark studied.

Technical Summary: Resonant Di-Higgs Searches in bbˉτ+τ−b\bar{b}\tau^+\tau^-bbˉτ+τ− at HL-LHC

More like this

Technical Summary: Resonant Di-Higgs Searches in $b\bar{b}\tau^+\tau^-$ at HL-LHC