⚛️ phenomenology

Subleading Effects in Soft-Gluon Emission at One-Loop in Massive QCD

This paper completes the one-loop soft-gluon approximation for massive QCD amplitudes by deriving a new soft operator that preserves on-shell conditions and momentum conservation, while also providing the necessary subleading collinear expansion for massless quark-antiquark pairs.

Original authors: Michał Czakon, Kilian Erhard Minguez, Felix Eschment

Published 2026-02-13

📖 4 min read🧠 Deep dive

CC BY 4.0

Original authors: Michał Czakon, Kilian Erhard Minguez, Felix Eschment

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 you are trying to predict the outcome of a chaotic dance party where particles are the dancers. In the world of Quantum Chromodynamics (QCD), which is the physics of how quarks and gluons interact, these dances are incredibly complex.

This paper is about understanding what happens when one specific dancer—a gluon—decides to leave the party very quietly and very slowly. In physics terms, this is called a "soft" emission.

Here is the breakdown of the paper's story, using simple analogies:

1. The Problem: The "Whispering" Gluon

Usually, when physicists calculate how particles scatter, they look at the main action. But sometimes, a gluon is emitted with almost zero energy. It's like a dancer whispering a secret before sneaking out the door.

The Old Way: In the past, scientists had a good rulebook for predicting what happens when these "whispering" gluons leave, but only if all the other dancers were massless (like light beams).
The Missing Piece: Real-world particles, like the top quark, are heavy (they have mass). The old rulebook didn't work for them. If you tried to use the old rules for heavy dancers, the math would break. This paper provides the missing instruction manual for how heavy particles behave when a soft gluon leaves the party.

2. The Solution: The "Magic Remote Control"

The authors created a new mathematical tool, which they call a "Soft Operator."

Think of the entire particle collision as a complex machine. The "Soft Operator" is like a magic remote control that you can point at the machine.

What it does: Instead of recalculating the entire chaotic dance from scratch, you just press a button on the remote.
The Result: The remote instantly tells you how the machine changes when that one quiet gluon leaves. It adjusts the "color" (a type of charge in particle physics) and the "spin" (how the particle is rotating) of the remaining heavy dancers, while making sure they don't break the laws of physics (like momentum conservation).

3. The Twist: Heavy vs. Light Dancers

In the world of massless particles, the math is surprisingly simple because the dancers don't care about their weight. But heavy particles are different.

The Analogy: Imagine a feather (massless) and a bowling ball (massive) both trying to whisper. The feather's whisper barely affects its movement. The bowling ball, however, has to shift its weight significantly to whisper.
The Discovery: The authors found that for heavy particles, there are extra "hidden" effects coming from the hard, violent parts of the collision that we usually ignore. These effects are like the recoil of a cannon when it fires a tiny bullet. Even though the bullet is small, the heavy cannon moves back. The paper calculates exactly how much that "cannon" (the heavy quark) moves back.

4. The "Collinear" Puzzle (The Side Quest)

To build their main remote control, the authors had to solve a side puzzle. They needed to know what happens when a massless particle splits into two others that are moving in almost the exact same direction (collinear).

The Metaphor: Imagine a single stream of water splitting into two streams that are so close together they look like one.
The Breakthrough: They found a missing piece of the puzzle for when this split involves a quark and an anti-quark. Interestingly, to solve this, they had to look at the "high-energy limit," which is like zooming out so far that the details blur, revealing a simple, elegant pattern.

5. Why Should You Care?

You might ask, "Why do we need to know how a heavy quark whispers?"

Precision is Key: Modern particle physics (like at the Large Hadron Collider) is so precise that we can't afford to ignore these "whispers." If we want to find new physics or understand the Higgs boson better, we need to subtract the background noise perfectly.
The Result: This paper gives physicists the exact formula to subtract that noise for heavy particles. It's like giving a sound engineer a perfect noise-canceling headphone specifically tuned for heavy bass.

Summary

In short, this paper is the final piece of the puzzle for predicting how heavy particles behave when they emit a very low-energy gluon. The authors built a "magic remote control" (the soft operator) that simplifies incredibly complex math, allowing scientists to make more accurate predictions about the fundamental building blocks of our universe. They also fixed a side problem regarding how particles split when moving in the same direction, ensuring the whole theory holds together.

1. Problem Statement

The paper addresses a critical gap in the perturbative calculation of Quantum Chromodynamics (QCD) scattering amplitudes involving massive quarks. Specifically, it aims to provide the missing ingredient required to approximate one-loop amplitudes in the limit where a single gluon becomes soft (vanishing energy, $q \to 0$ ) up to terms suppressed by the gluon energy (subleading power).

While the subleading soft behavior for massless QCD was established in a previous work by the same authors [1], the presence of massive quarks introduces new complexities:

Mass Singularities: Massive quarks do not exhibit the same collinear singularities as massless ones, altering the factorization structure.
Hard Virtual Contributions: In the massive case, contributions from hard virtual momenta appear in the soft expansion that are not contained in the one-loop amplitude of the process without the gluon.
Incomplete Factorization: Previous attempts to describe the splitting $g \to q\bar{q}$ at tree-level in the subleading collinear limit lacked a universal expression for the term involving the subleading soft asymptotics of an (anti-)quark.

The goal is to derive a complete one-loop soft operator ( $S^{(1)}$ ) that acts on color and spin space, manipulating hard parton momenta while keeping them on-shell and conserving momentum, specifically for processes involving massive quarks.

2. Methodology

The authors employ a combination of expansion-by-regions, Ward identities, and diagrammatic analysis to derive the soft operator.

Expansion-by-Regions: The calculation separates the loop integrals into "soft" and "hard" regions.
- Soft Region: The loop momentum $l$ and soft gluon momentum $q$ are both treated as small ( $l, q \sim \lambda$ ). The integrands are expanded in these variables.
- Hard Region: The loop momentum is large, while the soft gluon momentum is expanded. The authors identify a "mismatch" where the hard-region contribution does not simply factorize into the tree-level soft operator acting on the one-loop amplitude. This mismatch is resolved by identifying specific hard-region contributions that must be included in the soft operator itself.
Ward Identities and Gauge Invariance: The derivation heavily relies on the Ward identity to relate internal-line gluon emissions to derivatives of the non-radiative amplitude. This allows the replacement of complex internal diagrams with differential operators acting on the external states.
Master Integrals: The result is expressed in terms of six master integrals ( $A_0, M_{1i}, M_{1j}, M_{2i}, M_{2j}, M_3$ ). The authors provide the explicit Laurent expansions of these integrals in dimensional regularization ( $d=4-2\epsilon$ ) for various mass configurations (both massive, one massless, both massless).
Tree-Level Collinear Limit: To close the system of equations for the one-loop soft operator, the authors derive the subleading collinear limit for the tree-level splitting $g \to q\bar{q}$ . This involves analyzing the high-energy limit where the quark and anti-quark become collinear.

3. Key Contributions

A. The One-Loop Soft Operator ( $S^{(1)}$ )

The primary result is the explicit formula for the one-loop soft operator acting on massive quarks. It is decomposed into non-abelian (NA) and abelian (A) contributions:
$P_{g} S^{(1)} = \sum_{i \neq j} \left( \dots \right)^{\epsilon} \left[ 2if^{abc}T^a_i T^b_j S^{(1)NA}_{ij} + T^c_i T_i \cdot T_j S^{(1)A}_{ij} \right] + \text{Single-line terms}$

Non-Abelian Part: Contains the standard color structure and coefficients ( $g^{(1)}_{ij}, a_i, b_{ij}, c_{ij}$ ) dependent on master integrals and kinematic variables ( $\alpha_i, \alpha_j, v$ ).
Abelian Part: A new feature for massive QCD (absent in massless QCD at this order) that agrees with QED results. It involves terms proportional to the field-strength tensor $F_{\mu\nu}$ and spin operators.
Hard-Region Contributions: The authors explicitly include terms proportional to the master integral $A_0(m^2)$ (which vanishes for massless particles). These terms arise from the hard region of the loop expansion and are necessary to maintain gauge invariance and correct the factorization formula. They correspond to corrections to the chromomagnetic operator in Heavy Quark Effective Theory (HQET).

B. Subleading Collinear Limit for $g \to q\bar{q}$

The paper provides the complete expression for the subleading behavior of tree-level amplitudes when a gluon splits into a massless quark-antiquark pair ( $g \to q\bar{q}$ ).

The result is expressed in terms of universal factors and process-dependent amplitudes.
Crucially, it identifies a term corresponding to the limit $|x| \to \infty$ (where $x$ is the momentum fraction), which effectively "eikonalizes" the collinear quarks. This term is represented by a straight eikonal line with multiple gluon attachments, a high-energy limit feature.

C. Spurious Pole Cancellation

The authors demonstrate that the pole terms ( $1/\epsilon^2, 1/\epsilon$ ) of their derived soft operator correctly reproduce the known infrared singularities of massive one-loop gauge-theory amplitudes. They show that when massless partons are involved, the inclusion of jet operators and flavor-off-diagonal soft contributions ensures the cancellation of spurious poles, matching the universal singularity structure.

4. Results and Validation

Analytic Expressions: The paper provides closed-form analytic expressions for the coefficients of the soft operator for all mass configurations ( $m_i, m_j \neq 0$ ; one massless; both massless).
Numerical Verification: The authors performed a numerical study for two processes:
1. $e^+e^- \to t\bar{t}g$ (massive top quarks).
2. $e^+e^- \to b\bar{b}t\bar{t}g$ (mixed masses).
  Using the Recola library for amplitudes and OneLOop for integrals, they calculated the relative error ( $\Delta$ ) between the exact one-loop amplitude and the soft approximation.
- Finding: The subleading power (NLP) approximation shows a relative error scaling as $O(q_0/\sqrt{s})$ , confirming the validity of the derived formula. The leading power (LP) approximation shows a larger error, as expected.

5. Significance and Future Directions

Completeness of QCD Factorization: This work completes the set of ingredients necessary to approximate one-loop QCD amplitudes with massive quarks up to the first non-singular term in the soft expansion. This is essential for high-precision phenomenology at the LHC, particularly for processes involving top quarks or heavy BSM particles.
Bridge to Effective Theories: The identification of hard-region contributions links the soft expansion to HQET and the matching coefficients of the chromomagnetic operator, providing a deeper theoretical understanding of factorization at subleading power.
Resummation Frameworks: The authors note that while their factorization formula is compact, it is currently written in a form distinct from the SCET operator formalism. They suggest that their results could serve as a starting point for improving existing resummation frameworks or developing new ones that handle massive particles more naturally.
Future Work: The authors identify the need to extend these results to:
- Emission of two or more soft gluons at one-loop.
- Single gluon emission beyond one-loop.
- The regime where the gluon energy is comparable to the quark mass ( $q_0 \sim m$ ), requiring non-trivial jet operators dependent on the energy-to-mass ratio.

In summary, this paper provides the rigorous mathematical framework and explicit formulas required to handle subleading soft-gluon effects in massive QCD at the one-loop level, resolving previous inconsistencies and enabling more accurate theoretical predictions for collider physics.