Factorization theorem for quasi-TMD distributions with kinematic power corrections

This paper derives a frame-invariant factorization theorem for quasi-TMD distributions that incorporates all-order kinematic power corrections, revealing a non-multiplicative TMD evolution structure and demonstrating that these corrections significantly improve the agreement between lattice simulations and phenomenological extractions of the Collins-Soper kernel.

The Gist

Imagine you are trying to take a high-resolution photograph of a tiny, fast-moving particle inside a proton. In the world of particle physics, this is called a "Tomography" (like a CT scan, but for subatomic particles).

To do this, scientists use two main tools:

1. Real-world experiments: Smashing particles together in giant colliders (like the LHC).
2. Computer simulations: Using "Lattice QCD," which is like a giant 3D grid where they simulate the laws of physics on a supercomputer.

The Problem: The "Blurry" Snapshot
The computer simulations are powerful, but they have a limitation. To get a clear picture, the particle needs to be moving very, very fast. However, current supercomputers can only simulate particles moving at "moderate" speeds.

Because the particles aren't moving fast enough, the standard mathematical formulas used to interpret the data are slightly "blurry." They ignore a subtle effect: Kinematic Power Corrections (KPCs).

Think of it like this:

• The Standard Formula (Leading Power): Imagine you are estimating the speed of a car by looking at its speedometer. You assume the car is driving on a perfectly flat, straight highway. This works great if the car is going 200 mph.
• The Reality (KPCs): But what if the car is going slower, or the road has bumps? The speedometer reading isn't wrong, but your interpretation of it is slightly off because you ignored the bumps and the car's actual weight. In physics, these "bumps" are the particle's sideways wobble (transverse momentum) and the fact that the proton has mass.

The Breakthrough: A Sharper Lens
This paper, written by physicists from the University of Madrid, derives a new, sharper mathematical lens (a Factorization Theorem) that accounts for these "bumps" and "wobbles" all the way up to the highest possible precision.

Here is the core discovery, explained simply:

1. The "Recipe" Changed

Previously, scientists used a recipe that said: "Take the computer data, multiply it by this one magic number (the Soft Factor), and you get the real physics." It was a simple multiplication.

The new paper says: "Wait, that's too simple. Because the particle is wobbling, you can't just multiply. You have to mix the data with a complex recipe."

• Analogy: Instead of just adding a pinch of salt to a soup (multiplication), you now have to blend the salt into the broth in a specific way that depends on how much broth you have (convolution). The math is more complex, but it's much more accurate.

2. Why It Matters (The "Aha!" Moment)

The authors tested their new formula against existing computer data. They found that the old, blurry formulas were underestimating the true values by about 10% to 20%.

In the world of science, a 20% error is huge! It's like a weather forecast saying "20% chance of rain" when it's actually going to pour.

3. Solving the Mystery

For a while, there was a disagreement between:

• The Computer Simulations: Which gave one set of numbers.
• The Real-World Experiments: Which gave a slightly different set of numbers.

Scientists were confused. "Why don't they match?"

The authors of this paper showed that if you use their new, sharper lens (including the Kinematic Power Corrections), the computer simulation numbers suddenly match perfectly with the real-world experiments.

The Takeaway

This paper is like finding the missing piece of a puzzle.

• Before: We thought the computer simulations and real experiments disagreed because the simulations were "wrong."
• Now: We realize the simulations were actually right all along; we just needed a better mathematical tool to read them correctly.

By fixing the "blur" caused by the particles' wobbles, the authors have reconciled the gap between theory and experiment, giving us a much clearer, more accurate map of how the tiny building blocks of our universe are arranged.

Technical Summary

1. Problem Statement

Lattice QCD simulations of nonlocal operators offer a unique pathway to determine Transverse-Momentum-Dependent (TMD) parton distribution functions (TMDPDFs) and the Collins–Soper (CS) kernel. These calculations rely on quasi-TMD correlators, which are defined in Euclidean space with a finite hadron momentum $P_z$.

The standard theoretical framework for extracting TMDs from these correlators relies on a Leading-Power (LP) factorization theorem, which assumes the limit $P_z \to \infty$. However, current lattice simulations are limited to moderate momenta ($P_z \simeq 3$ GeV), which is only slightly larger than the nonperturbative QCD scale ($\Lambda_{QCD}$). Consequently, power corrections (terms suppressed by powers of $1/P_z$) are expected to be significant.

While various types of power corrections exist (genuine higher-twist, target-mass, $q_T$-corrections), the authors argue that Kinematic Power Corrections (KPCs) are the most critical. KPCs arise from the non-zero transverse momentum of partons ($k_T$) and are necessary to restore frame invariance (reparameterization invariance) and other symmetries broken at the LP level. Previous studies only addressed KPCs up to Next-to-Leading Power (NLP); this work aims to derive the factorization theorem including KPCs to all orders.

2. Methodology

The authors derive the factorization theorem using the background field method and the concept of TMD-twist.

• TMD-Twist Decomposition: They classify operators by their TMD-twist. KPCs correspond specifically to the series of terms involving only twist-two TMD distributions. This allows the authors to isolate KPCs from genuine higher-twist effects (which involve multiparton correlations) by computing the factorization theorem with vanishing background gluon fields and using the QCD equations of motion to re-express "bad" quark field components in terms of "good" ones.
• Derivation Strategy:
  1. Tree-Level: They perform algebraic manipulations on the quasi-TMD matrix element, expressing the result in terms of standard TMD correlators and a reduced soft factor.
  2. Loop-Level: They compute the one-loop coefficient function in general kinematics (without expanding in $1/P_z$). They demonstrate that the coefficient function remains the same as the LP result but with a modified Lorentz-invariant argument.
  3. Rapidity Renormalization: They address the rapidity renormalization scales ($\zeta, \bar{\zeta}$). In the KPC formalism, the condition $\zeta \bar{\zeta} = \mu^2 (2v \cdot k)^2$ replaces the standard LP condition. To simplify the analysis, they set $\zeta = \bar{\zeta}$, which eliminates boost-invariant logarithmic complications while preserving the essential kinematic structure.

3. Key Contributions

The paper presents a TMD-with-KPC factorization theorem (Eq. 23 and 25) with the following distinct features:

• Frame Invariance: The derived expression is explicitly frame-invariant, correcting the symmetry violations present in the LP approximation.
• Non-Multiplicative Evolution: Unlike the LP case where the evolution factor and soft factor factorize multiplicatively, the KPC expression involves a convolution. The parton's longitudinal momentum fraction $\xi$ becomes dependent on the transverse momentum $k_T$:
  $$\xi = \frac{x}{2} \left( 1 + \sqrt{1 + \frac{k_T^2}{x^2 P_v^2}} \right)$$
  This kinematic constraint prevents the simple separation of the evolution factor from the TMD distribution.
• All-Orders Validity: The derivation captures KPCs to all orders in the power expansion, restricted to twist-two distributions.
• Reduced Soft Factor: The reduced soft factor still factorizes multiplicatively, allowing it to be computed and removed separately, preserving the utility of quasi-TMDs as clean observables.

4. Results and Numerical Estimates

The authors perform a numerical analysis using the ART25 phenomenological TMD extraction as input to estimate the magnitude of KPCs in current lattice setups.

• Magnitude of Corrections:
  • For typical lattice momenta ($P_v \sim 1.7$ GeV), the KPCs induce a 10–20% correction to the reduced quasi-TMD correlator compared to the LP approximation.
  • The corrections are negative, implying that LP-based lattice extractions underestimate the absolute value of the TMD distributions.
  • The corrections are highly dependent on the transverse distance $b$ and the momentum fraction $x$. They are most significant at small and large $x$, but exhibit a stable plateau in the region $0.3 \lesssim x \lesssim 0.7$.
• Impact on the Collins–Soper (CS) Kernel:
  • The CS kernel is typically extracted from the ratio of quasi-TMDs at different momenta. The authors show that using the LP formula on data generated with KPCs leads to a mismatch between the extracted kernel and the true input kernel.
  • However, when the KPC-corrected formalism is applied to the same data, the extracted CS kernel aligns perfectly with the ART25 input.
  • This suggests that previous discrepancies between lattice results and phenomenological extractions may be resolved by including KPCs.
• Impact on TMDPDFs:
  • Similar to the CS kernel, the extracted TMDPDFs using the LP formula underestimate the true values by 10–20% in the relevant kinematic regions.

5. Significance

• Theoretical Advancement: This work provides the first complete, all-orders factorization theorem for quasi-TMDs including kinematic power corrections, establishing a rigorous framework for high-precision lattice QCD.
• Resolution of Discrepancies: The study demonstrates that the inclusion of KPCs significantly improves the agreement between lattice simulations and phenomenological extractions (like ART25). It suggests that current "disagreements" in the literature may be artifacts of neglecting these kinematic corrections.
• Guidance for Future Simulations: The results indicate that while increasing $P_z$ reduces KPCs, current lattice momenta ($P_z \sim 1-2$ GeV) are insufficient to neglect them. Future analyses must employ the TMD-with-KPC formalism to avoid systematic errors of 10–20% in the determination of nucleon structure.
• Methodological Shift: The paper highlights that the evolution factor cannot be simply factored out in position space when KPCs are included; instead, a convolution integral over $k_T$ is required, fundamentally changing the data analysis strategy for lattice TMDs.