The N$^3$LO Twist-2 Matching of Linearly Polarized Gluon TMDs

This paper computes the twist-2 matching of transverse-momentum-dependent linearly polarized gluon parton distribution and fragmentation functions at next-to-next-to-next-to-leading order (N$^3$LO) in QCD, supplemented by next-to-next-to-leading logarithmic (NNLL) small-$x$ resummation, to provide high-precision theoretical inputs for future Electron-Ion Collider studies of hadron spin structure and three-dimensional tomography.

The Gist

Imagine the universe is built out of tiny, invisible Lego bricks. The most fundamental of these bricks are quarks and gluons. Quarks are the bricks that stick together to form protons and neutrons (the stuff inside your body), while gluons are the "glue" that holds them together.

For a long time, physicists have been trying to understand exactly how these bricks are arranged inside a proton. It's not just a static pile of Lego; it's a swirling, chaotic storm of particles moving at nearly the speed of light.

This paper is a massive, high-precision update to the "instruction manual" that physicists use to map out this storm. Here is the breakdown in simple terms:

1. The Problem: The "Glue" is Lopsided

Most people think of the glue (gluons) inside a proton as just being there, holding things together. But gluons have a weird property: they can be linearly polarized.

• The Analogy: Imagine a spinning top. Usually, we just care that it's spinning. But imagine if the top was also wobbling in a specific, flat direction, like a coin spinning on a table. That "wobble" is the linear polarization.
• Why it matters: This wobble changes how the proton behaves when it smashes into other things. If we ignore this wobble, our predictions about particle collisions are slightly off.

2. The Goal: A Better Map

Physicists use a tool called TMD (Transverse Momentum Dependent) functions to map the proton. Think of this as a 3D GPS for the particles inside.

• The Old Map: Previous versions of this map were good, but they were like a low-resolution photo. You could see the general shape, but the details were blurry.
• The New Map: This paper provides a 4K Ultra-HD version of the map. The authors have calculated the "matching coefficients" (the mathematical rules that translate the blurry photo into a sharp image) up to a level of precision called N3LO.

What does N3LO mean?
In physics, calculations are done in steps of complexity.

• LO (Leading Order): A rough sketch.
• NLO (Next-to-LO): A better drawing.
• NNLO: A detailed painting.
• N3LO (Next-to-Next-to-Next-to-LO): This is the "God-tier" calculation. It includes tiny, almost invisible corrections that only show up when you look extremely closely. It's the difference between guessing the weather and knowing exactly when the rain will start, down to the second.

3. The Challenge: The "Small x" Fog

There is a tricky part of the proton map where things get very messy. This is the "small x" region, where the gluons are carrying very little energy but there are millions of them.

• The Analogy: Imagine trying to count individual raindrops in a hurricane. At the edge of the storm (small x), the rain is so dense and chaotic that standard counting methods fail.
• The Solution: The authors didn't just calculate the map; they also added a special "fog-clearing lens" called resummation. This mathematical trick allows them to predict what happens in that chaotic, high-energy zone without the numbers blowing up.

4. Why Should You Care? The Future of the "Super Microscope"

You might ask, "Why do we need a map this detailed?"
The answer is the Electron-Ion Collider (EIC), a giant particle accelerator currently being built. It will act like a super-microscope, taking 3D movies of protons and neutrons.

• The Connection: To make sense of the pictures the EIC takes, scientists need a perfect theoretical map to compare against. If the map is blurry, the EIC data will look confusing. If the map is this new, ultra-precise one, scientists can finally see the "wobble" of the gluons clearly.
• The Payoff: This helps us understand:
  • How the proton gets its spin (why it spins).
  • How matter is built from nothing but energy.
  • The fundamental forces that hold the universe together.

Summary

Think of this paper as the ultimate software update for the physics of the proton.

1. They calculated the rules for how "wobbly" gluons behave with unprecedented precision (N3LO).
2. They fixed the math for the chaotic, high-energy zones (Small-x resummation).
3. They double-checked their work using a theoretical "supersymmetry" test (like a spell-check for physics equations) to make sure it's perfect.

This work ensures that when the Electron-Ion Collider turns on, we won't just be taking pictures; we'll be able to read the story the proton is telling us.

Technical Summary

1. Problem Statement

The paper addresses a critical gap in the theoretical framework of Transverse Momentum Dependent (TMD) factorization in Quantum Chromodynamics (QCD). Specifically, it focuses on the linearly polarized gluon distribution, denoted as $h_1^{\perp g}$.

• Context: While the unpolarized gluon distribution ($f_1^g$) matching coefficients are known up to Next-to-Next-to-Next-to-Leading Order (N3LO), the matching coefficients for the linearly polarized gluon TMDs were previously only known up to Next-to-Next-to-Leading Order (NNLO).
• Significance: The linearly polarized gluon distribution contributes to azimuthal asymmetries in various high-energy processes (e.g., diphoton production, Higgs transverse momentum, dijets). Precision predictions for future facilities like the Electron-Ion Collider (EIC) require N3LO accuracy to match experimental precision and to properly constrain the 3D tomography of nucleons.
• Gap: The N3LO matching coefficient for $h_1^{\perp g}$ was unknown, limiting the precision of fixed-order calculations and the reliability of small-$x$ resummation.

2. Methodology

The author employs a rigorous perturbative QCD approach within the Soft-Collinear Effective Theory (SCET) framework.

• Operator Definition & Calculation:

  • The calculation starts from the SCET definition of the gluon beam function (TMD PDF) and fragmentation function (TMD FF) using gauge-invariant collinear fields.
  • The matching is performed by computing the bare gluon beam function via splitting amplitudes integrated over the collinear phase space.
  • Regulator: An exponential rapidity regulator ($e^{-b_0 \tau P \cdot K / P^+}$) is utilized to handle rapidity divergences, a method previously established in the author's work.
  • Dimensional Regularization: Calculations are performed in $d = 4 - 2\epsilon$ dimensions.

• Renormalization & Factorization:

  • The bare coefficient functions are renormalized using multiplicative operator renormalization constants ($Z_B$) and zero-bin subtractions ($S_{0b}$) to remove soft and collinear overlaps.
  • The Renormalization Group (RG) equations for the coefficient functions are solved to separate scale-dependent logarithms ($L_\perp, L_Q, L_\nu$) from the scale-independent matching coefficients.

• Numerical Fitting:

  • Due to the complexity of the analytic expressions (involving harmonic polylogarithms and multiple scales), the author provides numerical fits for the scale-independent parts of the N3LO coefficients.
  • The fitting functions use elementary logarithmic terms ($\ln x, \ln(1-x)$) and polynomial expansions, achieving an accuracy better than $10^{-3}$ across the kinematic range.

• Small-$x/z$ Resummation:

  • For TMD Fragmentation Functions (FFs), the author performs NNLL (Next-to-Next-to-Leading Logarithmic) resummation of small-$z$ (or small-$x$) logarithms.
  • This is achieved using Mellin space techniques, where small-$z$ logarithms map to poles in the Mellin variable $N$. An ansatz for the asymptotic behavior of unrenormalized collinear functions is used to extract the resummation coefficients.

3. Key Contributions

1. First N3LO Calculation: This work presents the first analytic computation of the N3LO twist-2 matching coefficients for linearly polarized gluon TMDs ($h_1^{\perp g}$) for both Parton Distribution Functions (PDFs) and Fragmentation Functions (FFs).
2. Numerical Implementation: The full analytic results are provided in ancillary files, while the paper presents high-precision numerical fits for immediate phenomenological use.
3. NNLL Resummation for FFs: The paper derives the NNLL resummation of small-$z$ logarithms specifically for the linearly polarized gluon TMD FFs, extending previous work that focused mainly on unpolarized distributions.
4. Supersymmetry Sum Rule Check: The author verifies the N=1 supersymmetry momentum conservation sum rule at the three-loop level. By setting $C_A = C_F = N_f$ and $T_F = 1/2$, the integral of the difference between gluon-gluon and gluon-quark matching coefficients vanishes, providing a non-trivial consistency check for the complex N3LO calculation.

4. Key Results

• Matching Coefficients ($I'_{gi}$ and $C'_{ig}$):
  • Explicit expressions for the N3LO coefficients for the $q \to g$ and $g \to g$ channels are provided.
  • Threshold Limit ($z \to 1$): The behavior near the threshold is analyzed, showing specific logarithmic structures ($\ln(1-z)$) and constant terms.
  • Small-$x/z$ Limit: The leading power behavior in the small-$x$ limit is confirmed to agree with previous Leading Logarithmic (LL) predictions. For FFs, the singlet sector exhibits double logarithmic divergences ($\ln^2 z$) at each order, which are resummed.
• Perturbative Convergence:
  • Analysis of the first moment of the linearly polarized gluon TMD PDF shows that perturbative uncertainties are well-controlled for $x > 10^{-3}$ once N3LO corrections are included.
  • For very small $x$ ($x < 10^{-3}$), resummation effects become dominant, necessitating the inclusion of the derived resummation terms.
• Phenomenological Application (Higgs Production):
  • The results are applied to the transverse momentum ($p_T$) distribution of Higgs bosons at the LHC ($13$ TeV).
  • The study demonstrates that the linearly polarized gluon contribution, while suppressed by $\alpha_s$, induces measurable azimuthal modifications. The N3LO corrections stabilize the prediction for the cumulative cross-section at low $p_T$.
• Resummation Impact:
  • For TMD FFs, the resummed NNLL predictions are compared with fixed-order results. The resummation remains crucial for $z < 10^{-2}$, significantly altering the shape of the coefficient functions compared to fixed-order truncations.

5. Significance

• Precision Physics: These results provide the essential theoretical inputs required for high-precision TMD phenomenology. They allow for global fits of TMDs at N3LO accuracy, reducing theoretical uncertainties in extracting parton distribution functions.
• EIC Readiness: The Electron-Ion Collider (EIC) aims to perform unprecedented measurements of gluon polarization and 3D nucleon structure. This work equips the community with the necessary fixed-order and resummed tools to interpret future EIC data on semi-inclusive DIS and dijet production.
• Theoretical Consistency: The verification of the N=1 SUSY sum rule at three loops serves as a robust validation of the calculation method, increasing confidence in the N3LO results for other TMD observables.
• Small-$x$ Dynamics: The derivation of NNLL resummation for linearly polarized gluon FFs opens new avenues for studying gluon saturation and high-energy QCD dynamics where small-$x$ effects are dominant.

In summary, this paper represents a major advancement in perturbative QCD, bridging the gap between NNLO and N3LO precision for a crucial spin-dependent observable, thereby enabling the next generation of high-energy physics experiments to probe the internal spin structure of the nucleon with greater fidelity.