The N$^3$LO Twist-2 Matching of Helicity TMDs and SIDIS $q_st$ Spectrum

We compute the twist-2 matching of transverse momentum dependent (TMD) helicity parton distribution and fragmentation functions at next-to-next-to-next-to-leading order (N$^3$LO) in QCD. This calculation entails the complete set of next-to-next-to-leading order (NNLO) Dokshitzer-Gribov-Lipatov-Altarelli-Parisi (DGLAP) splitting functions govering the evolution of helicity-dependent parton distribution functions (PDFs) and fragmentation functions (FFs). Within TMD factorization framework, we quantify the impact of radiative corrections by completing the next-to-next-to-next-to-leading logarithmic (N$^3$LL) prediction for lepton-hadron transverse momentum imbalance in semi-inclusive deep inelastic scattering (SIDIS). Our results provide the most precise theoretical input for probing the helicity structure and confined motion of quarks and gluons at future electron-ion collider (EIC).

💡 Research Summary

The paper presents a landmark calculation in perturbative QCD: the twist‑2 matching of helicity (spin‑dependent) transverse‑momentum‑dependent (TMD) parton distribution functions (PDFs) and fragmentation functions (FFs) at next‑to‑next‑to‑next‑to‑leading order (N³LO). The authors first derive the complete set of next‑to‑next‑to‑leading order (NNLO) polarized DGLAP splitting functions, which govern the scale evolution of helicity PDFs and FFs. These splitting functions, previously known only up to NLO, are obtained analytically for both space‑like and time‑like processes in the MS scheme, with careful treatment of γ₅ and scheme‑dependent terms.

Using Soft‑Collinear Effective Theory (SCET), the paper constructs a factorization theorem for a novel observable in semi‑inclusive deep‑inelastic scattering (SIDIS): the transverse‑momentum imbalance q* between the scattered lepton and the observed hadron. Unlike the traditional q_T spectrum, q* is defined in the lepton‑target frame and probes the component of the lepton’s transverse momentum orthogonal to the event plane. The authors introduce an exponential rapidity regulator for both soft and collinear sectors, perform zero‑bin subtractions, and demonstrate that the soft function reduces to the standard TMD soft factor within this scheme.

The factorized cross‑section reads schematically as a product of a hard matching coefficient C(Q²,μ), a soft function S(b,μ,ν), a beam function B_f/N(x,b,μ,ν) (the TMD PDF), and a fragmentation function D_h/f(z,b,μ,ν) (the TMD FF). Spin dependence is isolated via a Fierz decomposition, yielding separate contributions proportional to the unpolarized and helicity beam functions, ΔB_f/N. The hard coefficient is computed to three loops, and the rapidity evolution is resummed to next‑to‑next‑to‑next‑to‑leading logarithmic (N³LL) accuracy.

The core technical achievement is the N³LO matching of helicity TMDs. Starting from the operator definitions of the helicity TMD PDF and FF, the authors perform collinear mass factorization, renormalization‑group (RG) evolution, and scheme transformations. They evaluate three‑loop diagrams using modern algebraic tools (FORM, FIRE, LiteRed) and extract the finite matching kernels. The resulting kernels contain polylogarithmic structures up to weight five, ζ‑values, and logarithms of the impact‑parameter b. Small‑x expansions are provided, showing that the new terms can modify the low‑x behavior of helicity distributions, a region of particular interest for the proton spin puzzle.

Armed with these analytic results, the authors present a phenomenological study of the q* spectrum in polarized SIDIS at an Electron‑Ion Collider (EIC). They implement the N³LO matching coefficients together with the NNLO polarized splitting functions, evolve PDFs and FFs using the five‑loop QCD beta function, and perform N³LL resummation of rapidity logarithms. Numerical results indicate that the N³LO corrections shift the q* distribution by roughly 5–10 % relative to the N²LO prediction, with the most pronounced effect at low q* (≲ 1 GeV). The helicity‑dependent part of the cross‑section, σ_L, exhibits a noticeable enhancement, implying that future high‑precision spin asymmetry measurements at the EIC will be sensitive to the newly computed higher‑order terms.

In the conclusion, the authors emphasize that this work places helicity‑dependent TMD observables on the same theoretical footing as their unpolarized counterparts, providing the most accurate perturbative input available. They outline several immediate applications: global fits of polarized TMDs, improved extractions of the gluon helicity distribution at small x, and refined Monte‑Carlo implementations of polarized parton showers. Appendices contain the full flavor decomposition of the NNLO splitting functions, checks of RG consistency, the explicit form of the QCD beta function up to five loops, anomalous dimensions, and the renormalization constants used throughout the calculation.

Overall, the paper delivers a comprehensive, state‑of‑the‑art calculation that will become a cornerstone for spin‑dependent QCD phenomenology, especially in the era of the EIC where unprecedented precision on helicity observables is expected.