⚛️ phenomenology

Precise predictions for double Higgs production in association with a vector boson in Effective Field Theory

This paper presents next-to-next-to-leading order QCD predictions for double Higgs production in association with a vector boson ( $pp \to Vhh$ ) within both SMEFT and HEFT frameworks, providing inclusive cross-sections at LHC energies and analyzing how QCD corrections factorize differently for $W^\pm hh$ versus $Zhh$ channels due to the latter's gluon-induced component.

Original authors: Ramona Gröber, Michał Ryczkowski, Gioia Sacchi

Published 2026-07-31

📖 6 min read🧠 Deep dive

CC BY 4.0

Original authors: Ramona Gröber, Michał Ryczkowski, Gioia Sacchi

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

Imagine the universe is a giant, cosmic Lego set, but instead of plastic bricks, it's built from tiny, invisible particles. For decades, scientists have been trying to figure out exactly how these pieces snap together. The most famous piece in the box is the Higgs boson, a particle that gives everything else its "heft" or mass. But here's the tricky part: we don't just want to know how one Higgs works; we want to know how two of them interact with each other. It's like trying to understand how two magnets push or pull when they get close. If we can see how two Higgs bosons behave together, we might finally crack the code on why the universe has the shape and rules it does.

To do this, scientists use the world's biggest particle collider, the Large Hadron Collider (LHC), which smashes protons together at nearly the speed of light. Sometimes, these crashes produce a Higgs boson along with a "vector boson," which is like a messenger particle that carries the weak force (think of it as a delivery truck bringing a package). The big question is: do these collisions happen exactly as our current rulebook (the Standard Model) predicts, or are there hidden rules we haven't found yet? This is where "Effective Field Theory" comes in. Think of it as a magnifying glass that lets scientists look for tiny cracks in the rulebook without needing to know the name of every single new particle that might be hiding in the cracks. It's a way to say, "If there are new, heavy rules, here is how they would subtly change the game."

The Paper: Hunting for Ghosts in the Machine

This paper is like a high-precision weather forecast for a very specific, very rare type of cosmic storm. The authors, Ramona Gröber, Michał Ryczkowski, and Gioia Sacchi, have calculated exactly what should happen when two Higgs bosons are born alongside a messenger particle (either a W or a Z boson) in a proton smash-up. They did this using the most advanced math available, pushing their calculations to "next-to-next-to-leading order" (NNLO). In the world of particle physics, this is like upgrading from a rough sketch to a 4K, high-definition movie. They didn't just guess; they built a super-detailed simulation to see how the numbers change if the universe follows the standard rules or if it follows the "magnifying glass" rules of new physics.

The Main Finding: The "Flat" vs. The "Wobbly"
The team discovered a fascinating split in behavior depending on which messenger particle is involved.

The W-Boson Case (The Steady Hand): When the Higgs pair is produced with a W boson (the charged messenger), the math behaves beautifully. The complex quantum corrections (the "noise" of the calculation) stay separate from the new physics rules. This means the "K-factor"—a number that tells you how much the real-world result differs from the simple guess—stays almost perfectly flat. No matter how much you tweak the new physics rules, the correction remains steady. It's like driving a car on a perfectly paved highway; no matter how fast you go, the road stays smooth.
The Z-Boson Case (The Wobbly Wheel): When the pair is produced with a Z boson (the neutral messenger), things get wobbly. Here, a sneaky, invisible process involving loops of heavy quarks (specifically top quarks) kicks in. This process only shows up at the highest level of calculation. Because of this, the K-factor changes depending on the new physics rules. It's like driving on a road that suddenly gets bumpy if you turn the steering wheel a certain way. This makes the Z-boson channel much more sensitive to detecting new physics, but also much harder to predict perfectly.

What They Found in the Numbers
The authors provided a massive toolkit for other scientists. Instead of just giving one number, they gave a set of "coefficients" (like a recipe) that allows anyone to quickly calculate the result for any combination of new physics rules within a certain range.

They found that for the W-boson channel, the theoretical uncertainty is very low, staying below 4%. This makes it a very clean place to look for deviations.
For the Z-boson channel, the uncertainty is slightly higher (over 5%) because of that tricky loop-induced process.
They showed that if the new physics rules are pushed to their current experimental limits, the production rate could jump significantly—up to 6 times the normal rate for certain scenarios. However, they also noted that current experimental constraints are still roughly 183 to 294 times larger than the Standard Model prediction, meaning we are still in the "searching" phase and haven't seen these jumps yet.

What They Rule Out (and What They Don't)
The paper explicitly clarifies that the QCD corrections (the messy math of the strong force) do not behave the same way for both W and Z bosons. They prove that for the W boson, the corrections largely factorise (stay separate) from the new physics effects, but for the Z boson, the loop-induced part changes the game entirely, creating a stronger dependence on the new physics rules.
They also clarify what they didn't do. They didn't calculate the next level of complexity (N3LO) for the gluon-induced part of the Z-boson process because it's too mathematically heavy for this specific study. They also didn't include every single possible new physics operator, focusing only on the most relevant ones for this specific collision.

How Sure Are They?
The authors are very confident in their calculations. They have simulated the process with extreme precision, including scale uncertainties, PDF (parton distribution function) uncertainties, and strong coupling constant ( $\alpha_s$ ) uncertainties. They are not claiming to have found new physics; they are saying, "Here is exactly what we expect to see if the Standard Model is right, and here is exactly how the numbers will wiggle if new physics is hiding in the cracks." Their predictions are ready for the High-Luminosity LHC (HL-LHC), the next upgrade of the collider, which will smash protons with much higher intensity.

In short, this paper builds a super-accurate map for a treasure hunt. The treasure is "new physics," and the map tells us exactly where to look (the Z-boson channel) and what to expect (a wobbly K-factor) when we get there. While the current map shows no treasure yet, it gives the explorers the best possible tools to find it when the HL-LHC turns up the volume.

Technical Summary: Precise predictions for double Higgs production in association with a vector boson in Effective Field Theory

Problem and Motivation
Higgs boson pair production ($hh$) is a flagship measurement for the High-Luminosity Large Hadron Collider (HL-LHC), offering direct access to the trilinear Higgs self-coupling and the shape of the Higgs potential. While gluon fusion is the dominant production mode, its cross-section is significantly suppressed compared to single Higgs production. Consequently, subleading channels, such as associated production with a weak gauge boson ( $pp \to Vhh$ with $V=W^\pm, Z$ ), are of considerable interest. These channels can probe whether electroweak symmetry breaking is realized linearly or non-linearly by testing the coupling patterns of one or two Higgs bosons to vector bosons.

Current experimental constraints on $Vhh$ production are still orders of magnitude above the Standard Model (SM) predictions. To exploit the discriminatory power of the HL-LHC, precise theoretical predictions are required. This work addresses the need for Next-to-Next-to-Leading Order (NNLO) Quantum Chromodynamics (QCD) predictions for $Vhh$ production within Effective Field Theory (EFT) frameworks, specifically the Standard Model Effective Field Theory (SMEFT) and the Higgs Effective Field Theory (HEFT).

Methodology
The authors calculate NNLO QCD corrections for $pp \to Vhh$ in both SMEFT (truncated at dimension six) and HEFT (leading order in the chiral expansion). The calculation relies on the factorization of the quark-induced production process into the Drell-Yan production of an off-shell vector boson ( $V^*$ ) and its subsequent decay into $Vhh$.

Factorization: The QCD corrections affect only the $pp \to V^*$ production vertex, which is a Drell-Yan process. These corrections, computed at NNLO in Refs. [23, 24], are applied to the EFT-modified decay amplitude.
SMEFT Implementation: The authors utilize the Warsaw basis for dimension-six operators. They distinguish between operators affecting electroweak vertices (modifying $V V h$ , $V V hh$, and $qqV$ couplings) and those inducing new Lorentz structures. The calculation includes contributions from operators $O_\phi$ , $O_{\phi\Box}$ , $O_{\phi D}$ , $O_{\phi W}$ , $O_{\phi B}$ , $O_{\phi WB}$ , and fermionic operators $O_{\phi q}^{(1,3)}$ , $O_{\phi u}$ , $O_{\phi d}$ , $O_{\phi t}$ , $O_{\phi Q}^{(1,3)}$ , and $O_{t\phi}$ . The trilinear Higgs coupling is modified by $O_\phi$ , $O_{\phi\Box}$ , and $O_{\phi D}$ .
HEFT Implementation: The calculation employs a multiplicative $\kappa$ -framework approach, modifying the SM couplings via factors $\kappa_V$ , $\kappa_\lambda$ , $\kappa_{2V}$ , and $\kappa_t$ . The HEFT power counting is separated from the perturbative QCD expansion, allowing the authors to remain at Leading Order (LO) in the HEFT expansion while including NNLO QCD corrections.
Gluon-Induced Contributions: A critical component of the $Zhh$ calculation is the inclusion of the loop-induced $gg \to Zhh$ subprocess. This contribution enters formally at NNLO QCD (via heavy-quark loops: triangle, box, and pentagon topologies) and does not interfere with the quark-induced process at this order. This channel introduces sensitivity to top-quark Yukawa couplings and $Ztt$ interactions.
Parametrization: Instead of providing raw cross-sections for every parameter point, the authors express the inclusive cross-sections as linear combinations of numerical coefficients ( $A_i$ for HEFT, $B_i$ for SMEFT) multiplied by the relevant EFT parameters. This allows for fast reconstruction of predictions for arbitrary EFT parameters within the considered ranges.

Key Results
The study provides inclusive cross-section predictions at $\sqrt{s} = 13.6$ and $14.0$ TeV, including uncertainties from renormalization/factorization scales, Parton Distribution Functions (PDFs), and $\alpha_s$ .

$W^\pm hh$ Production:
- The NNLO QCD corrections largely factorize from the EFT dependence.
- Consequently, the K-factors ( $K = \sigma_{NNLO}/\sigma_{LO}$ ) are almost flat across the parameter space, showing only mild dependence on EFT coefficients.
- Theoretical uncertainties remain small (below 4%) and relatively constant regardless of the EFT parameter values.
- Maximal cross-section enhancements can reach $\mathcal{O}(6)$ for extreme values of $\kappa_\lambda$ (e.g., $\kappa_\lambda = 6.1$ ).
$Zhh$ Production:
- The behavior differs significantly due to the gluon-induced $gg \to Zhh$ component.
- The K-factor exhibits a strong dependence on EFT coefficients, particularly $\kappa_\lambda$ (varying between 1.25 and 1.6) and top-Yukawa related parameters.
- Theoretical uncertainties are generally higher (up to 5% or more) compared to $W^\pm hh$ , driven by the gluon-initiated contribution and PDF uncertainties in the high- $x$ region for specific Wilson coefficients (e.g., $C_{\phi q}^{(1)}$ ).
- Maximal enhancements can reach $\mathcal{O}(6)$ for $\kappa_\lambda = 6.1$ , and up to $10 \times \sigma_{SM}$ for specific SMEFT Wilson coefficients like $C_{\phi q}^{(1)} = -3/\text{TeV}^2$ .
SMEFT vs. HEFT:
- In SMEFT, the $Zhh$ channel shows significant sensitivity to operators affecting the top Yukawa coupling ( $C_{t\phi}$ ) and the $Ztt$ interaction, which are absent in the $W^\pm hh$ channel at this order.
- The paper provides tables of numerical coefficients ( $A_i$ and $B_i$ ) for both frameworks at LO and NNLO, enabling efficient phenomenological studies.

Significance and Claims
The authors claim that their results constitute a necessary step toward a systematic EFT program for sub-leading Higgs pair production channels at the HL-LHC.

Precision: The inclusion of NNLO QCD corrections reduces theoretical uncertainties to a level sufficient for the expected precision of future HL-LHC measurements.
Discrimination Power: The distinct behavior of K-factors between $W^\pm hh$ (flat) and $Zhh$ (dependent on EFT parameters due to the gluon component) offers a tool to disentangle different EFT contributions.
Current Status: The authors note that current experimental constraints ( $183 \times \sigma_{SM}$ and $294 \times \sigma_{SM}$ ) are still an order of magnitude larger than the maximal theoretical enhancements predicted in their EFT scans. However, they assert that the gap between experimental sensitivity and theoretical predictions is expected to narrow substantially with the HL-LHC upgrade.
Utility: The provided numerical coefficients allow for the rapid evaluation of cross-sections for arbitrary EFT parameters, facilitating immediate use in experimental analyses and global fits.

The paper concludes that while current data cannot yet probe the predicted EFT enhancements, the theoretical framework established here is robust and ready for the high-precision era of the HL-LHC. Future extensions mentioned include fully differential EFT predictions and a more complete treatment of higher-order operators in the HEFT expansion.

The Paper: Hunting for Ghosts in the Machine

More like this