Deeply virtual pion production through two-loop order

This paper presents the first calculation of next-to-next-to-leading order (NNLO) QCD radiative corrections for deeply virtual pion production, demonstrating that these two-loop corrections substantially improve the agreement between perturbative QCD predictions and experimental data from JLab while also refining theoretical descriptions of transverse single-spin asymmetries for future facilities like the EIC and EicC.

The Gist

The Big Picture: Taking a 3D X-Ray of a Proton

Imagine a proton (a tiny particle inside an atom) not as a solid marble, but as a bustling city filled with smaller, invisible residents called quarks. For a long time, scientists have only had a "flat map" of this city, showing how many residents live there and how fast they move. But they wanted a 3D hologram to see exactly where the residents are located in space and how they move together.

To build this hologram, scientists use a process called Deeply Virtual Meson Production (DVMP). Think of this as firing a high-speed, virtual "flashbulb" (a photon) at the proton city. The flash hits a resident, who then pops out of the city as a new particle (a pion), leaving a "scuff mark" on the city's structure. By studying these scuff marks, scientists can reconstruct the 3D map of the proton.

The Problem: The Blueprint Was Outdated

To interpret these scuff marks, scientists need a mathematical "blueprint" (theory) to predict what should happen.

• The Old Blueprint: For about 20 years, the best blueprint scientists had was like a sketch drawn with a pencil. It was good, but it missed a lot of fine details. In physics terms, this was the "Next-to-Leading Order" (NLO) calculation.
• The Reality Check: When scientists compared this old sketch to real data from the Jefferson Lab (JLab), the lines didn't quite match up. The prediction was off.

The Solution: A Super-Computer Upgrade (NNLO)

The authors of this paper decided to upgrade the blueprint. They performed a massive calculation called Next-to-Next-to-Leading Order (NNLO).

• The Analogy: If the old calculation was like a sketch, the new NNLO calculation is like a high-definition, 3D architectural rendering that includes every tiny bolt, wire, and shadow.
• The Work: They had to calculate the interactions of particles through "two loops." Imagine a particle traveling a path, but instead of going straight, it takes a detour, loops back, interacts with itself, and then continues. Doing this math for two loops is incredibly complex—like trying to solve a puzzle where every piece is moving and changing shape.

The Key Discovery: The "Pure Singlet" Puzzle Piece

One of the hardest parts of this job was a specific type of interaction called the "Pure Singlet" contribution.

• The Metaphor: Imagine trying to hear a whisper in a noisy room. Most of the noise (the "Non-Singlet" part) is loud and easy to hear. But the "Pure Singlet" part is a very quiet, specific frequency that gets drowned out by the noise and the rules of quantum mechanics (specifically a tricky math problem involving a symbol called $\gamma_5$).
• The Breakthrough: The team developed a clever new method to isolate this quiet whisper without getting confused by the noise. They successfully calculated this piece for the first time.

The Results: The Map Finally Fits

When they added these new, high-definition corrections to their predictions, something amazing happened:

1. The Fit Improved: The new predictions lined up much better with the actual data collected at JLab. It was like taking a blurry photo and suddenly sharpening the focus until the details were crystal clear.
2. The Correction Was Huge: The new math didn't just add a tiny tweak; it added a substantial boost. In some cases, the correction was so large it doubled the predicted signal. This proves that to get an accurate map of the proton, you must include these complex, two-loop details.
3. Future Proofing: The authors show that this high-precision blueprint is essential for future experiments at big facilities like the Electron-Ion Collider (EIC). Without this new level of detail, future experiments would be trying to navigate with an outdated map.

What About the "Spin"?

The paper also looked at something called Transverse Single-Spin Asymmetry (TSSA).

• The Analogy: Imagine spinning a top. If you hit it from the side, does it wobble left or right? This asymmetry tells us about the "spin" of the proton's residents.
• The Finding: The new, complex math didn't change the size of this wobble much (it was already stable), but it confirmed that the direction and shape of the wobble depend heavily on how we model the proton's internal structure. It acts like a sensitive test to see which model of the proton is correct.

Summary

In short, this paper is about upgrading the math used to understand the inner structure of protons. The authors built a much more precise, "two-loop" version of the theory. When they used this new version, their predictions matched real-world experiments much better than before. This means we are finally getting a clear, high-resolution 3D picture of how the building blocks of our universe are arranged.

Technical Summary

1. Problem Statement

Deeply Virtual Meson Production (DVMP), specifically the production of pions ($\pi^+, \pi^0$) in electron-proton scattering ($\gamma^* p \to \pi N$), is a "gold-plated" channel for extracting Generalized Parton Distributions (GPDs). GPDs are essential for understanding the 3D internal structure of the nucleon, its spin decomposition, and mechanical properties.

While experimental facilities like JLab, the future Electron-Ion Collider (EIC), and EicC are generating high-precision data, the theoretical predictions have lagged.

• Current Status: Phenomenological analyses of DVMP have historically relied on Next-to-Leading Order (NLO) QCD calculations, which were established two decades ago.
• The Gap: Recent advancements in Deeply Virtual Compton Scattering (DVCS) have reached Next-to-Next-to-Leading Order (NNLO). However, DVMP calculations remained at NLO. Given that NLO corrections for DVMP are sizable, there is a critical uncertainty regarding whether NNLO corrections are necessary for the precise extraction of GPDs from upcoming high-luminosity data.
• Specific Challenge: Calculating the NNLO corrections for DVMP involves complex two-loop diagrams. A specific technical hurdle exists for neutral pion ($\pi^0$) production: the "pure singlet" (PS) contribution, which first appears at the two-loop level, involves $\gamma_5$ matrices that are ill-defined in dimensional regularization when using standard covariant projector methods.

2. Methodology

The authors performed the first NNLO QCD calculation for the longitudinal cross-section and transverse single-spin asymmetries (TSSA) of $\gamma^*_L p \to \pi^+ n$ and $\gamma^*_L p \to \pi^0 p$ within the collinear factorization framework.

Key Technical Steps:

• Factorization Framework: The amplitude is factorized into a perturbative hard-scattering kernel ($T$), a twist-2 pion distribution amplitude ($\phi_\pi$), and nucleon axial-vector GPDs ($\tilde{H}, \tilde{E}$).
• Non-Singlet (NS) Channel ($\pi^+$ and part of $\pi^0$):
  • Instead of calculating the two-loop diagrams from scratch, the authors utilized a "shortcut." They established a mapping between the hard-scattering kernel of DVMP and the electromagnetic form factor (EMFF) of the pion.
  • By replacing the kinematic variables in the known NNLO coefficient function of the $\pi^+$ EMFF, they derived the NNLO coefficient function for DV$\pi^+$P.
• Pure Singlet (PS) Channel ($\pi^0$ only):
  • This channel receives no contribution at tree-level or one-loop due to color and C-parity conservation; it first appears at two-loops.
  • Handling $\gamma_5$: To avoid inconsistencies with $\gamma_5$ in dimensional regularization, the authors avoided taking Dirac traces. Instead, they kept external spinor indices open, expanded the amplitude in a basis of totally anti-symmetric gamma matrices, and enforced the Ward-Takahashi identity (current conservation).
  • This method allowed them to isolate UV and IR finite results for the PS contribution.
• Numerical Implementation:
  • Tools: Used QGRAF and HepLib for diagram generation; Apart and FIRE for integral reduction; differential equation methods for master integrals; and AmpRed for analytic expression generation.
  • Inputs: Adopted the RQCD lattice QCD prediction for the pion distribution amplitude ($a_2$) and two popular GPD parametrizations: the Goloskokov-Kroll (GK) model and the GUMP (Universal Moment Parametrization) model.
  • Scales: Renormalization and factorization scales were varied to estimate theoretical uncertainties.

3. Key Contributions

1. First NNLO Calculation: This work presents the first complete NNLO QCD calculation for deeply virtual pion production, covering both charged ($\pi^+$) and neutral ($\pi^0$) channels.
2. Pure Singlet Derivation: The authors successfully derived the analytic expression for the two-loop pure singlet hard-scattering kernel for $\pi^0$ production, overcoming the technical difficulties associated with $\gamma_5$ in multi-loop calculations.
3. Methodological Shortcut: Demonstrated the validity of inferring DVMP coefficients from pion EMFF corrections via analytic continuation, significantly reducing computational complexity for the non-singlet sector.

4. Results

The study analyzed differential longitudinal cross-sections ($d\sigma_L/dt$) and Transverse Single-Spin Asymmetries (TSSA, $A_{UT}$) at kinematics relevant to JLab, EicC, and EIC.

• Magnitude of Corrections:
  • The NNLO corrections are positive and substantial.
  • For JLab kinematics ($Q^2 \approx 4$ GeV$^2$), the NNLO correction is as significant as the NLO correction, enhancing the cross-section by over 100% in some intermediate $Q^2$ regions.
  • The relative importance of NNLO corrections decreases as $Q^2$ increases, but they remain non-negligible for precision physics.
• Comparison with Data:
  • $\pi^+$ Production: Including NNLO corrections significantly improves the agreement between perturbative QCD predictions and existing JLab experimental data, particularly when using the GK model.
  • $\pi^0$ Production: The NNLO prediction based on the GUMP parametrization aligns well with unpolarized cross-section data. However, the GK model (dominated by the pion pole term) still underestimates the data even with NNLO corrections, suggesting potential limitations in the model's truncation or the pion pole dominance assumption.
• Pure Singlet Impact:
  • The PS contribution to the cross-section is modest (reducing the NNLO prediction by only $\sim 1\%$).
  • However, the PS contribution has a more noticeable effect on TSSA, altering predictions by up to 10% at large $|t|$.
• TSSA Behavior:
  • The magnitude of TSSA is robust against NNLO corrections (similar to NLO findings).
  • However, the $t$-dependence of TSSA is highly sensitive to the GPD parametrization. The GK model predicts a large asymmetry at low $t$ due to the pion pole, whereas GUMP does not. This divergence offers a crucial test for the pion-pole dominance hypothesis.

5. Significance

• Essential for GPD Extraction: The results demonstrate that NNLO accuracy is indispensable for the reliable extraction of nucleon GPDs from DVMP data. Relying on NLO predictions would lead to significant systematic errors, especially in the kinematic regimes of JLab and the upcoming EicC/EIC.
• Future Facilities: The improved theoretical precision is critical for interpreting data from the JLab 22 GeV upgrade, the US EIC, and China's EicC, where high-precision tomography of the nucleon is a primary goal.
• Model Discrimination: The sensitivity of TSSA to GPD models (specifically the pion pole contribution) provides a new observable to distinguish between different GPD parametrizations and validate theoretical models of nucleon structure.
• Theoretical Benchmark: The successful calculation of the pure singlet two-loop contribution sets a new standard for handling complex multi-loop QCD processes involving $\gamma_5$ and singlet channels.