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

We provide the last missing ingredient necessary to approximate one-loop amplitudes in QCD with massive quarks in the limit of vanishing energy of a single gluon up to terms suppressed by this energy. Our main result is a soft operator acting in color and spin space that manipulates the momenta of the hard partons while keeping them on-shell and respecting momentum conservation. Additionally, we provide a complete expression for the subleading term of the expansion of an arbitrary tree-level amplitude in the limit where the momenta of a massless quark and a massless anti-quark of the same flavor become collinear. This limit is necessary to obtain the one-loop soft approximation whenever the process involves such a quark-anti-quark pair. Interestingly, the result involves a high-energy limit.

💡 Research Summary

The paper addresses a long‑standing gap in the theoretical description of QCD amplitudes involving massive quarks when a single gluon becomes soft. Building on the authors’ earlier work on subleading soft‑gluon behavior in massless QCD, the present study derives a universal soft operator that is valid at one‑loop order for processes with any number of massive (and possibly massless) partons, provided the emitted gluon carries vanishing energy. The operator acts in the combined color‑and‑spin space, reshuffles the momenta of the hard partons while keeping them on‑shell, and respects overall momentum conservation.

The soft operator, denoted S₁^{(q)}, is decomposed into a non‑abelian (NA) piece and an abelian (A) piece. The NA contribution generalizes the familiar eikonal structure by incorporating explicit mass dependence through scalar products s_{ij}=2p_i·p_j, s_{iq}=2p_i·q, etc. The A contribution is entirely new for massive QCD; it reproduces the known QED subleading soft factor and vanishes in the massless limit. Both pieces are expressed as a double sum over unordered pairs of external partons (i,j). Color factors appear as products of generators T_i^a T_j^b, while spin effects are encoded in little‑group rotation operators K_i^{μν} acting on the polarization or spinor of parton i.

A crucial technical ingredient is a set of six master integrals that capture all loop‑momentum dependence of the operator. One is the standard massive two‑point function A₀(p,m²); the remaining five are specific scalar four‑point integrals I_{α₁α₂α₃α₄} with various powers of propagators. Their ε‑expansions (to O(ε⁰)) are reproduced from Ref.