⚛️ phenomenology

Unitarity limits on HH production from VBF mediated channels

This paper derives joint unitarity constraints on anomalous quartic $VVHH$ couplings arising from four-derivative effective operators in both SMEFT and HEFT frameworks, specifically analyzing their impact on off-shell Vector Boson Fusion mediated double Higgs production at high-energy colliders.

Original authors: A. Dobado, D. Domenech, M. Herrero, R. A. Morales

Published 2026-10-07

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: A. Dobado, D. Domenech, M. Herrero, R. A. Morales

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

In the vast, high-stakes arena of particle physics, scientists are constantly searching for cracks in the Standard Model, the prevailing theory that describes how the universe's fundamental particles interact. One of the most promising places to look for new physics is the behavior of the Higgs boson, the particle responsible for giving mass to other particles. While the Standard Model predicts how Higgs bosons should behave, it does not account for every possible scenario. If the model is incomplete, there might be hidden forces or new particles influencing how Higgs bosons are created, particularly when they appear in pairs. To test this, researchers use powerful accelerators to smash particles together at incredible speeds, hoping to catch a glimpse of these rare double-Higgs events. However, as the energy of these collisions increases, the mathematical descriptions used to predict what happens can sometimes break down, producing impossible results that violate the basic laws of probability. This creates a puzzle: how can scientists use these theories to search for new physics without the math itself becoming nonsensical?

A team of physicists has tackled this problem by focusing on a specific way Higgs pairs are produced, known as vector boson fusion. In this process, two incoming particles, such as quarks or electrons, exchange force-carrying particles called vector bosons, which then merge to create two Higgs bosons. The researchers were interested in how this process behaves when the theory includes "anomalous" interactions—hypothetical deviations from the Standard Model that would signal new physics. These deviations are often described using mathematical tools called effective field theories, which act like a map for exploring energy levels beyond what current experiments can directly reach. The team realized that previous studies had tried to set safety limits on these theories by looking at the intermediate step where the force-carrying particles collide as if they were real, stable particles. But in the actual collisions happening inside accelerators, these intermediate particles are fleeting and unstable, never existing as free, stable objects. By treating them as if they were stable, earlier calculations were missing a crucial piece of the physical reality.

To correct this, the researchers performed a detailed simulation of the entire process, from the initial collision of the fermions to the final appearance of the two Higgs bosons and the two outgoing fermions. They used advanced computational techniques to break down the complex interactions into simpler components based on angular momentum, a property that describes how the particles spin and move relative to each other. By analyzing these components, they could determine exactly where the mathematical predictions would start to fail and violate the laws of physics. They found that the limits on how strong these new interactions can be depend heavily on the total energy of the collision and the specific type of particles involved. For instance, collisions involving up and down quarks, which are the building blocks of protons, impose much stricter limits than those involving electrons. The study revealed that the constraints on these hypothetical interactions are significantly different when calculated for the full, realistic process compared to the simplified, intermediate steps used in the past.

The team discovered that the most restrictive limits come from collisions where the incoming particles have specific spin orientations, and these limits become tighter as the energy of the collision increases. At an energy of 10 teraelectronvolts, a scale relevant for future high-energy colliders, the strength of the interaction involving one of the key parameters, which they call delta, must be less than about 0.04 to keep the theory consistent. The other parameter, eta, is allowed to be slightly larger, but still within strict bounds. These findings are crucial because they define the "safe zone" for future experiments. If scientists ignore these limits and try to interpret data using values outside this zone, their conclusions about new physics could be based on mathematical errors rather than real discoveries. The researchers also compared two different methods for applying these safety limits to experimental data. One method, which looks at the mass of the Higgs pair alone, tends to cut off valid data too early, making the search for new physics less sensitive. The other method, which considers the mass of the entire final group of particles, preserves more of the data while still ensuring the math remains valid.

This work provides a clearer roadmap for the next generation of particle colliders, such as the High-Luminosity Large Hadron Collider and proposed future machines like the Future Circular Collider. By establishing precise boundaries for where these theories remain valid, the researchers ensure that when experimentalists eventually find a deviation from the Standard Model, it will be a genuine discovery and not an artifact of a broken calculation. The study emphasizes that for the most accurate and reliable results, future analyses must account for the fact that the intermediate particles in these fusion processes are never truly on their own, but are always part of a larger, complex interaction. This shift in perspective allows for a more honest and robust search for the hidden layers of the universe, ensuring that the tools used to explore the unknown are as solid as the theories they test.

Technical Summary: Unitarity Limits on HH Production from VBF Mediated Channels

Problem Statement
The study of double Higgs boson production via Vector Boson Fusion (VBF) channels ( $f_1f_2 \to f_3f_4H_5H_6$ ) is a primary window for probing Beyond Standard Model (BSM) physics, specifically anomalous quartic couplings between two electroweak gauge bosons and two Higgs bosons ( $VVHH$ ). These interactions are often parameterized within Effective Field Theories (EFTs), such as the Standard Model EFT (SMEFT) at dimension-8 and the Higgs EFT (HEFT) at chiral dimension-4, using operators containing four derivatives.

A critical issue arises because these 4-derivative operators induce effective vertices that grow with energy ( $O(p^2)$ ), leading to a violation of perturbative unitarity in scattering amplitudes at high energies. Previous studies on unitarity constraints for the relevant coefficients ( $\delta$ and $\eta$ ) have largely focused on the Vector Boson Scattering (VBS) subprocess ( $VV \to HH$ ), assuming the initial vector bosons are on-shell. However, in actual collider processes (both hadron-hadron and lepton-lepton), the intermediate $VV$ pairs in VBF channels propagate off-shell. Consequently, applying unitarity limits derived from on-shell $VV \to HH$ subprocesses to the full $f_1f_2 \to f_3f_4H_5H_6$ process may not accurately reflect the physical constraints at high-energy colliders.

Methodology
The authors perform a detailed analysis of unitarity constraints directly on the on-shell scattering process $f_1f_2 \to f_3f_4H_5H_6$ , where $f$ represents fermions (quarks or leptons). The methodology involves:

Theoretical Framework: The study utilizes both HEFT and SMEFT, focusing on 4-derivative operators. The relevant BSM contributions are parametrized by two coefficients, $\delta$ and $\eta$ , which are related to specific Wilson coefficients in both EFT frameworks (e.g., $c^{(1,2,3)}_{\Phi4D4}$ in SMEFT and $add_{VV1,2}$ in HEFT).
Partial Wave Analysis for 2 $\to$ 4 Scattering: Unlike standard 2 $\to$ 2 analyses, the authors address the complexity of a four-particle final state. They employ the massless spinor-helicity formalism to project helicity amplitudes onto a basis of kinematic polynomials with definite total angular momentum $J$ .
Numerical Integration: The partial wave amplitudes ( $T^J$ ) are computed by integrating over the 4-body phase space. The authors use a spherical parametrization and the VEGAS integrator to handle the non-trivial integrals arising from the off-shell propagators of the intermediate $VV$ bosons.
Unitarity Condition: The perturbative unitarity condition is imposed as $|T^J| < 1$ . The analysis considers projections for $J=0, 1, 2$ and various helicity configurations for processes including $ud \to duHH$ , $uu \to uuHH$ , $dd \to ddHH$ , and $e^-e^+ \to e^-e^+HH$ .
Clipping Method Comparison: The paper compares two methods for ensuring unitarity-safe predictions at colliders:
- Clipping based on the invariant mass of the final state jets and Higgs bosons ( $m_{HHjj} = \sqrt{s_{12}}$ ), derived from the current study.
- Clipping based on the invariant mass of the Higgs pair ( $m_{HH} = \sqrt{s_{56}}$ ), derived from the traditional on-shell $VV \to HH$ analysis.

Key Contributions and Results

Derivation of On-Shell Constraints: The paper provides the first derivation of unitarity limits specifically for the full VBF-mediated process $f_1f_2 \to f_3f_4H_5H_6$ , acknowledging the off-shell nature of the intermediate vector bosons.
Dominant Channels: The analysis identifies that the tightest unitarity constraints arise from the $J=0$ partial wave with same-sign negative helicities ( $(-,-; -,-, 0, 0)$ ). Among the processes studied, $ud \to duHH$ (involving both WBF and ZBF contributions) yields the most restrictive bounds due to the dominance of the WBF contribution.
Energy and Coefficient Dependence: The unitarity limits are presented as functions of the total initial energy $\sqrt{s_{12}}$ and the coefficients $\delta$ and $\eta$ . The results show that constraints on $\delta$ are generally tighter than those on $\eta$ . Furthermore, the constraints for hadron-hadron collisions (driven by valence quarks) are significantly stronger than for lepton-lepton collisions at comparable energies.
Comparison with VBS Limits: The authors demonstrate that the unitarity limits derived from the full on-shell process $f_1f_2 \to f_3f_4H_5H_6$ are looser than those derived from the on-shell subprocess $VV \to HH$ . This suggests that using $VV \to HH$ limits for collider analyses may be overly conservative or physically inappropriate for VBF topologies.
Implications for Clipping: In the context of LHC predictions (14 TeV), the paper shows that applying a clipping method based on $m_{HHjj}$ (the correct variable for the full process) preserves unitarity while allowing for a broader range of coefficient values compared to the traditional $m_{HH}$ clipping. The $m_{HH}$ method cuts off the cross-section at lower coefficient values, potentially underestimating the sensitivity to BSM physics.

Significance and Claims
The paper claims that for a realistic and safe analysis of BSM effective interactions in double Higgs production at future high-energy colliders (such as HL-LHC, HE-LHC, FCC-hh, and muon colliders), unitarity constraints must be applied to the full on-shell process rather than the intermediate subprocess.

The authors argue that the standard "clipping" method, which typically relies on the invariant mass of the vector boson pair or Higgs pair derived from on-shell VBS studies, is insufficient for VBF channels where the intermediate bosons are off-shell. They propose that using the invariant mass of the entire final state ( $m_{f_3f_4H_5H_6}$ ) as the clipping variable is the physically correct approach. This method ensures perturbative unitarity is preserved without unnecessarily restricting the parameter space of the EFT coefficients $\delta$ and $\eta$ , thereby providing a more accurate framework for interpreting future experimental data. The work emphasizes that while unitarity limits reduce the sensitivity to these parameters compared to uncut analyses, using the correct kinematic variable ( $m_{HHjj}$ ) is essential for deriving reliable exclusion limits.

More like this