⚛️ phenomenology

Constraining four-heavy-quark operators with top-quark, Higgs, and electroweak precision data

This paper establishes constraints on dimension-six four-heavy-quark operators in the Standard Model Effective Field Theory by combining LHC top-quark and single-Higgs production data with electroweak precision observables, while highlighting the significant impact of the $\gamma_5$ scheme choice on the resulting fits.

Original authors: Stefano Di Noi, Hesham El Faham, Ramona Gröber, Marco Vitti, Eleni Vryonidou

Published 2026-07-07

📖 6 min read🧠 Deep dive

CC BY 4.0

Original authors: Stefano Di Noi, Hesham El Faham, Ramona Gröber, Marco Vitti, Eleni Vryonidou

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 Standard Model of particle physics as a massive, incredibly detailed instruction manual for how the universe works. For decades, scientists have been checking this manual against reality using the Large Hadron Collider (LHC), which smashes particles together at near-light speeds. So far, they haven't found any "new" particles that aren't in the manual.

Because they haven't found the "smoking gun" (a new particle), scientists are now looking for "fingerprints" left behind. They suspect that if new, heavy physics exists, it might be too heavy to create directly, but it could leave subtle, invisible ripples in the behavior of known particles. This is where SMEFT (Standard Model Effective Field Theory) comes in. Think of SMEFT as a magnifying glass that looks for tiny, weird distortions in the manual's instructions.

This paper focuses on a specific set of instructions involving heavy quarks (like the top quark, the heaviest known particle). Specifically, they are looking at "four-heavy-quark operators." In plain English, these are rules describing how four heavy particles might interact in a way the current manual doesn't fully explain.

Here is a breakdown of what the authors did, using some everyday analogies:

1. The Detective Work: Mixing Different Clues

Usually, to catch a criminal (new physics), you look at the crime scene directly. In particle physics, this means looking at collisions where four top quarks are created at once. However, this is like trying to find a needle in a haystack while wearing thick gloves; the data is messy, and the "needle" is hard to see.

The authors decided to be smarter detectives. Instead of just looking at the messy crime scene, they combined three different types of clues:

The Direct Hit: Looking at four-top-quark collisions (the messy crime scene).
The Side Effect: Looking at how top quarks pair up or pair up with a Higgs boson (a lighter, more stable particle).
The Echo: Looking at precision measurements of the Higgs boson and other electroweak particles (like the Z boson).

By combining these, they hoped to tighten the net around the "criminal." They found that while the direct hits were fuzzy, the "echoes" (precision data) helped pin down the rules much more tightly.

2. The "Left-Handed" vs. "Right-Handed" Problem (The $\gamma_5$ Scheme)

This is the most technical and fascinating part of the paper. In the world of quantum mechanics, particles have a property called "chirality" (think of it as being left-handed or right-handed). To calculate how these particles behave in complex loops (like a particle going around a track and coming back), physicists use a mathematical tool called dimensional regularization.

Imagine you are trying to draw a perfect 3D cube on a 2D piece of paper. You have to make a choice about how to flatten the depth. In particle physics, there are two main ways to do this "flattening" for a specific mathematical object called $\gamma_5$ (which determines left vs. right):

The NDR Scheme: A "naïve" way of doing it.
The BMHV Scheme: A more rigorous, but trickier way.

The authors discovered something surprising: The choice of how you flatten the math changes the result.

It's like measuring a room with two different rulers. If you use Ruler A, the room is 10 feet wide. If you use Ruler B, it's 10.5 feet wide. Both rulers are "correct" in their own system, but if you don't realize you switched rulers, you'll think the room changed size.

The paper shows that if you only look at four-quark interactions and use the "NDR" ruler, you get one set of limits on the new physics. If you use the "BMHV" ruler, you get a different set of limits. In fact, in one scheme, the math creates a "flat direction"—a blind spot where the data can't tell the difference between two different types of interactions. In the other scheme, that blind spot disappears.

The Lesson: You cannot just pick a math trick and run with it. If you want to know the true limits of new physics, you have to be consistent, or you have to include more types of interactions to cancel out the differences between the math tricks.

3. The Results: Tightening the Net

The authors ran a massive statistical fit (a giant puzzle solver) using data from the LHC and previous experiments.

Linear vs. Quadratic: They looked at the data in two ways.
- Linear: Assuming the new physics is a small, gentle nudge.
- Quadratic: Assuming the new physics might be a stronger push, where the effects square themselves (like how a small increase in speed makes a car crash much more dangerous).
The Finding: When they included the "quadratic" effects (the stronger pushes), the constraints became much tighter. The "four-top-quark" data was the strongest driver here.
The Synergy: The most important finding is that combining the messy four-top-quark data with the precise Higgs and electroweak data gave the best results. The precise data filled in the gaps left by the messy data, and vice versa.

4. The Bottom Line

This paper is a warning and a guide for future physicists:

Don't ignore the math tricks: The way you handle the "left-handed/right-handed" math ( $\gamma_5$ ) actually changes the physical limits you set on new physics. You can't just pick one and ignore the other.
Combine your clues: You can't rely on just one type of experiment. You need the messy, high-energy collisions and the precise, low-energy measurements to get a clear picture.
The "Four-Heavy" operators are still a mystery: While they have tightened the rules, these specific interactions involving four heavy quarks are still among the least constrained parts of the Standard Model. There is still plenty of room for new physics to hide there.

In short, the authors built a better net to catch new physics, but they also realized that the shape of the net depends heavily on which mathematical scissors you use to cut the string. To get the right shape, you need to use all the tools in the toolbox.

Technical Summary: Constraining Four-Heavy-Quark Operators with Top-Quark, Higgs, and Electroweak Precision Data

Problem Statement
The absence of direct evidence for new light particles at the Large Hadron Collider (LHC) has necessitated a shift toward indirect searches within the Standard Model Effective Field Theory (SMEFT) framework. While global SMEFT analyses have placed constraints on many operator classes, operators involving four heavy-quark fields (specifically four-top-quark $t\bar{t}t\bar{t}$ , four-bottom-quark $b\bar{b}b\bar{b}$ , and mixed $t\bar{t}b\bar{b}$ interactions) remain among the least constrained. Current bounds on these operators are weak because they are primarily probed by inclusive multi-top production, which suffers from large experimental uncertainties and "flat directions" in the parameter space where different color and chirality structures cannot be distinguished. Furthermore, theoretical predictions for these operators often involve loop computations where the treatment of the chiral matrix $\gamma_5$ in dimensional regularization introduces scheme dependence that has not been fully quantified in the context of global fits.

Methodology
The authors perform a global fit to constrain dimension-six four-heavy-quark operators by synthesizing data from three distinct sectors:

Direct Probes: Leading-order (LO) predictions for four-top-quark ( $pp \to t\bar{t}t\bar{t}$ ) and top-bottom pair ( $pp \to t\bar{t}b\bar{b}$ ) production.
Indirect Loop-Induced Probes: Next-to-leading-order (NLO) predictions for top-pair ( $pp \to t\bar{t}$ ) and top-pair-plus-Higgs ( $pp \to t\bar{t}H$ ) production, which receive contributions from four-heavy operators at the one-loop level.
High-Precision Observables: Two-loop contributions to single-Higgs production via gluon fusion ( $gg \to H$ ) and Higgs decays ( $H \to gg, \gamma\gamma$ ), as well as two-loop corrections to Electroweak Precision Observables (EWPO).

The analysis utilizes the SMEFT framework with a specific flavor assumption ( $U(2)^2 \times U(3)^3$ ) focusing on third-generation quarks. Computations are performed using MadGraph5_aMC@NLO and the SMEFT@NLO package, with analytic expressions derived for Higgs processes and EWPO. A critical methodological component is the explicit comparison of two distinct $\gamma_5$ continuation schemes in dimensional regularization: the Naïve Dimensional Regularization (NDR) scheme and the Breitenlohner–Maison–'t Hooft–Veltman (BMHV) scheme. The authors demonstrate that while physical observables should be scheme-independent, restricting the fit to a subset of operators (four-heavy quarks only) without including the necessary counter-terms from lower-loop operators leads to scheme-dependent bounds.

Key Contributions

Comprehensive Global Fit: The paper presents a fit combining direct multi-top measurements, top-pair production, Higgs production/decays, and EWPO to constrain five specific four-heavy-quark Wilson coefficients ( $c_{QQ}^{(1,8)}, c_{Qt}^{(1,8)}, c_{tt}^{(1)}$ ).
Two-Loop EWPO and Higgs Calculations: The authors incorporate two-loop SMEFT corrections to EWPO and single-Higgs processes mediated by four-heavy operators, highlighting that these corrections are non-negligible, particularly for bottom-quark asymmetry observables ( $A_b, A_{b,FB}$ ).
Scheme Dependence Analysis: A dedicated study reveals that the choice of the $\gamma_5$ scheme (NDR vs. BMHV) significantly alters the extracted bounds on Wilson coefficients when only four-quark operators are considered. The authors show that in the BMHV scheme, a flat direction emerges in the $(c_{Qt}^{(1)}, c_{Qt}^{(8)})$ plane due to the vanishing of specific diagram topologies, whereas the NDR scheme lifts this degeneracy.
Operator Basis Clarification: The work clarifies the mapping between the Warsaw basis and the "top basis" (used in SMEFT@NLO), noting that evanescent operators can induce numerical shifts in the results depending on the definition of the color-octet operator.

Results

Linear vs. Quadratic Fits: In the linear fit (considering only $\mathcal{O}(\Lambda^{-2})$ terms), the combination of EWPO and top-quark production data is essential to reduce the allowed parameter space, as individual channels often exhibit flat directions. In the quadratic fit (including $\mathcal{O}(\Lambda^{-4})$ terms), the four-top-quark production channel provides the strongest constraints, significantly tightening the bounds compared to the linear case.
Complementarity: The study confirms that indirect probes (Higgs and EWPO) provide complementary information to direct multi-top production. While inclusive Higgs production probes only a subset of coefficients, the inclusion of EWPO is crucial for lifting degeneracies, particularly in the $(c_{QQ}^{(8)}, c_{QQ}^{(1)})$ plane where four-top data alone is blind.
Scheme Sensitivity: The fit results for single-Higgs production differ between NDR and BMHV schemes. In BMHV, the fit finds a flat direction for the linear combination $c_{Qt}^{(1)} + c_F c_{Qt}^{(8)}$ , whereas NDR distinguishes these coefficients. The authors attribute this to the fact that the BMHV scheme requires specific symmetry-restoring counterterms (involving other operators) to cancel scheme-dependent terms, which are absent when the fit is restricted to four-quark operators.
Numerical Bounds: The paper provides 95% Confidence Level (CL) bounds for all five Wilson coefficients. For instance, in the quadratic fit, the combination of all data yields tighter constraints than individual channels, with the four-top channel dominating the quadratic limits.

Significance
The paper establishes that constraining four-heavy-quark operators requires a global approach that leverages the complementarity between tree-level multi-top production and loop-induced precision observables. A primary finding is the demonstration that the interpretation of SMEFT limits on these operators is sensitive to the $\gamma_5$ regularization scheme if the analysis is not performed within a complete operator basis. The authors argue that to obtain scheme-independent results, one must include operators that enter at lower loop orders to compensate for the scheme dependence of the four-quark contributions. Consequently, the paper concludes that a comprehensive assessment of these operators necessitates their inclusion in a global fit alongside the complete set of relevant SMEFT operators, rather than in isolation. The work highlights that current experimental data, when combined with state-of-the-art two-loop theoretical predictions, can significantly improve constraints on these elusive interactions, though care must be taken regarding theoretical scheme choices.

1. The Detective Work: Mixing Different Clues

2. The "Left-Handed" vs. "Right-Handed" Problem (The γ5\gamma_5γ5​ Scheme)

3. The Results: Tightening the Net

4. The Bottom Line

More like this

2. The "Left-Handed" vs. "Right-Handed" Problem (The $\gamma_5$ Scheme)