Skewness-dependent moments of the pion GPD from nonlocal quark-bilinear correlators

This paper presents lattice QCD calculations of the odd Mellin moments of the pion valence-quark generalized parton distribution up to fifth order across a range of skewness values, utilizing boosted pion states and advanced renormalization techniques to extract skewness-dependent moments through polynomiality-constrained fits.

The Gist

Imagine the pion not as a solid marble, but as a bustling, fuzzy cloud of tiny particles called quarks and gluons. For decades, physicists have tried to map this cloud to understand how the universe's most fundamental forces hold matter together. Usually, they've only been able to take a "flat" snapshot of this cloud, seeing how the particles move forward. But this paper takes a giant leap forward by creating a 3D movie of the pion, showing how the cloud deforms and shifts when you look at it from different angles.

Here is a simple breakdown of what the researchers did and found:

1. The Challenge: Seeing the Invisible

Think of the pion's internal structure as a secret recipe. Scientists know the ingredients (quarks), but they can't see how they are arranged.

• The Old Way: Previous experiments were like looking at a shadow puppet. You could see the outline, but you couldn't tell if the puppet was leaning left or right, or how its arms were positioned. This is called "zero skewness"—looking straight on.
• The New Goal: The researchers wanted to see the "skewness." Imagine taking a photo of a spinning dancer. If you take the picture when they are facing you, it looks one way. If you take it when they are leaning sideways, the shape looks different. This paper is the first to successfully calculate what the pion looks like when it's "leaning" (non-zero skewness).

2. The Tool: A Super-Computer Microscope

To see these tiny particles, you can't use a regular microscope. You need Lattice QCD, which is like building a giant, digital grid (a lattice) of space and time.

• The Simulation: The team ran massive simulations on supercomputers. They created a virtual pion and "boosted" it to incredibly high speeds (up to 2.4 GeV).
• The Analogy: Imagine trying to study the wind inside a hurricane. If the hurricane is stationary, it's hard to see the details. But if you fly a plane through it at high speed, the wind patterns become clearer. By boosting the pion, the researchers could "freeze" the quantum fuzziness enough to take a clear picture of its internal structure.

3. The Method: Piecing Together the Puzzle

The researchers didn't just take one picture; they took thousands of snapshots from different angles and distances.

• The "Moments": They calculated specific mathematical "moments." Think of these like the average weight of the cloud at different distances from the center. They calculated up to the fifth "moment," which is like checking the cloud's shape very far out from the center.
• The "Polynomial" Rule: Nature has a rulebook. The shape of the pion must follow a specific mathematical pattern (called polynomiality). The researchers used this rule like a puzzle guide. Even though their data was a bit noisy, they knew the pieces had to fit a specific curve, which helped them solve the puzzle accurately.

4. The Results: What They Found

• The "Lean" Matters: They confirmed that as the pion "leans" more (higher skewness), the distribution of its internal particles changes. The particles don't just stay in a neat circle; the cloud stretches and shifts.
• Fading Out: They found that as you look further away from the center of the pion (higher momentum transfer) or as the pion leans more, the "weight" of the higher-order moments gets smaller. It's like the edges of the cloud become thinner and less significant.
• A New Contrast: Interestingly, they found that the pion behaves differently than a proton (the particle in the center of an atom). While a proton's internal structure shifts one way when it leans, the pion shifts the opposite way. It's as if the proton and pion are mirror images in how they react to being pushed.

5. Why It Matters (According to the Paper)

This work is a "first-principles" calculation, meaning they didn't guess; they calculated it directly from the laws of Quantum Chromodynamics (QCD).

• The Map: They have created the first reliable map of the pion's 3D structure that includes these "leaning" angles.
• The Future Guide: While the paper doesn't claim to cure diseases or build new engines, it provides a crucial "ground truth" for future experiments. Upcoming facilities like the Electron-Ion Collider will try to measure these same things in the real world. This paper gives those experimentalists a theoretical map to check their results against.

In short: The team used a super-computer to simulate a speeding pion, figured out how to measure its shape from different angles, and discovered that the pion's internal cloud deforms in a specific, predictable way that is opposite to how a proton deforms. They successfully mapped the first few layers of this 3D structure, setting a new standard for understanding the building blocks of matter.

Technical Summary

Technical Summary: Skewness-Dependent Moments of the Pion GPD from Nonlocal Quark-Bilinear Correlators

Problem Statement
Generalized Parton Distributions (GPDs) provide a comprehensive framework for understanding the three-dimensional structure of hadrons, encoding information on both longitudinal momentum and transverse spatial distributions. While GPDs are accessible experimentally through exclusive processes like Deeply Virtual Compton Scattering (DVCS), their extraction remains challenging due to convolution integrals. Lattice Quantum Chromodynamics (QCD) offers a first-principles approach to calculating GPDs. However, previous lattice studies have been largely restricted to the zero-skewness ($\xi=0$) limit. Calculations at non-zero skewness are significantly more complex due to the enlargement of the kinematic configuration space, increased computational costs, and the mixing of moments under Efremov–Radyushkin–Brodsky–Lepage (ERBL) scale evolution. Furthermore, traditional approaches based on local operators suffer from power-divergent operator mixing at higher dimensions, limiting the extraction of higher Mellin moments.

Methodology
This work presents a lattice QCD calculation of the odd Mellin moments of the pion valence-quark GPD up to the fifth order ($\langle x^4 \rangle$) for skewness values in the range $\xi \in [-0.33, 0]$. The study employs the following methodological framework:

• Lattice Setup: Calculations are performed on a $2+1$ flavor Highly Improved Staggered Quark (HISQ) gauge ensemble generated by the HotQCD collaboration with a lattice spacing of $a = 0.04$ fm. The valence sector utilizes the Wilson-Clover action with HYP smearing, tuned to a valence pion mass of 300 MeV. To suppress frame-dependent effects and power corrections, boosted pion states are employed with momenta up to $P_z \approx 2.428$ GeV and momentum transfers reaching $-t \approx 2.748$ GeV$^2$.
• Operator Formalism: Instead of local operators, the study utilizes nonlocal quark-bilinear operators within the Large-Momentum Effective Theory (LaMET) and Short-Distance Factorization (SDF) frameworks. The matrix elements of these operators are related to light-cone GPD moments via an Operator Product Expansion (OPE).
• Renormalization and Matching: Bare matrix elements are renormalized using the ratio scheme to remove linear and logarithmic ultraviolet divergences. The renormalized quasi-GPDs are matched to light-cone moments using perturbative matching coefficients.
  • For the zero-skewness case, the analysis employs Next-to-Next-to-Leading Logarithmic (NNLL) resummation based on Next-to-Next-to-Leading Order (NNLO) matching coefficients.
  • For the non-zero skewness case, the analysis utilizes Next-to-Leading Logarithmic (NLL) resummation based on Next-to-Leading Order (NLO) matching coefficients, including off-diagonal elements that account for moment mixing.
• Extraction Strategy: The extraction of moments involves simultaneous polynomiality-constrained fits to the lattice data across multiple values of $\xi$ and momentum transfer $t$. The polynomiality property, a fundamental consequence of Lorentz covariance, is enforced by parameterizing the Generalized Form Factors (GFFs) $A_{n,k}(t)$ using either a monopole-like form or a $z$-expansion. The $z$-expansion is adopted for the final analysis to provide more conservative uncertainty estimates.

Key Contributions and Results

• First Non-Zero Skewness Calculation: This work presents the first lattice QCD determination of pion GPD moments at non-zero skewness. The study successfully extracts the odd Mellin moments up to $\langle x^4 \rangle$ across the skewness range $[-0.33, 0]$.
• Validation at Zero Skewness: At $\xi=0$, the results for the lowest moments are consistent with previous lattice studies. The lowest moment $A_{1,0}$ reproduces the pion electromagnetic form factor, showing agreement with lattice results at the same pion mass (300 MeV) and experimental data, though discrepancies with physical mass experiments are attributed to the heavier valence pion mass used. Higher moments ($A_{3,0}$ and $A_{5,0}$) exhibit the expected monotonic decrease with increasing momentum transfer $-t$.
• Skewness Dependence: By combining matrix elements at multiple $\xi$ and $t$, the study extracts skewness-dependent GFFs ($A_{n,k}$ with $k>0$).
  • The third Mellin moment $H_3(\xi, t)$ is found to decrease with increasing $|\xi|$. The associated form factor $A_{3,2}(t)$ is negative, a feature distinct from the nucleon case where positive values have been reported.
  • The fifth Mellin moment $H_5(\xi, t)$ similarly shows suppression with increasing $|\xi|$ and $-t$.
  • The combined fits demonstrate that the extracted moments satisfy the polynomiality constraints required by Lorentz covariance, serving as a stringent internal consistency check.
• Perturbative Stability: The inclusion of NLL and NNLL resummation significantly improves the stability of the results compared to fixed-order calculations. The study identifies a stable short-distance fitting window ($0.08 \text{ fm} \leq z \leq 0.24 \text{ fm}$) where fixed-order, NLL, and NNLL results agree within statistical uncertainties.

Significance
The paper establishes that skewness-dependent pion GPD moments can be reliably determined on the lattice using nonlocal quark-bilinear correlators and OPE. By successfully extracting moments at non-zero skewness, this work provides crucial first-principles input for global analyses of experimental data, as most experimental measurements (e.g., at Jefferson Lab, AMBER, and the future EIC) are performed at non-zero skewness. The results offer new insights into the dynamical origin of skewness effects in the pion's three-dimensional structure, specifically demonstrating how partonic momentum and spatial distributions evolve beyond the forward limit. The study confirms that the lattice methodology can handle the complexities of moment mixing and frame dependence, paving the way for future calculations involving lighter pion masses, singlet and gluon contributions, and even moments.