Parton physics from a heavy-quark operator product expansion: Dynamical lattice QCD calculation of moments of the pion and kaon light-cone distribution amplitudes

This paper reports progress on calculating the first three nontrivial Mellin moments of kaon light-cone distribution amplitudes and summarizes recent continuum-limit results for the pion fourth moment using dynamical lattice QCD within the heavy-quark operator product expansion framework, demonstrating the method's feasibility for accessing higher moments.

The Gist

Imagine you are trying to understand the shape of a cloud. You can't just look at it from far away; you need to know how the water droplets are distributed inside it. In the world of particle physics, protons, neutrons, and particles like pions and kaons are like those clouds. They aren't solid balls; they are swirling bags of even smaller particles called quarks and gluons.

Physicists want to know the "recipe" of these particles: How is the momentum (the "oomph") shared between the quarks inside? This recipe is called the Light-Cone Distribution Amplitude (LCDA).

The problem is, this recipe is written in a language that is incredibly hard to read directly. It involves "light-like" separations, which are like trying to take a photo of a moving car with a camera that only works when the car is standing still. You can't do it directly on the supercomputers (called Lattice QCD) that physicists use to simulate the universe.

The Solution: The "Heavy Quark" Trick (HOPE)

This paper introduces a clever workaround called HOPE (Heavy-Quark Operator Product Expansion). Here is how it works, using an analogy:

The Analogy: The Heavy Anchor
Imagine you are trying to measure the shape of a tiny, wobbly boat (the pion or kaon) floating on a rough sea. It's too unstable to measure directly.

• The Trick: You drop a giant, heavy anchor (a fictitious heavy quark) onto the boat.
• The Result: The heavy anchor stabilizes the boat and creates a clear, measurable "shadow" or pattern in the water. Because the anchor is so heavy, the physics around it becomes simpler and easier to calculate, like a calm pond.
• The Math: Physicists use a mathematical tool called an Operator Product Expansion (OPE). Think of this as a translator. It takes the simple, calm "shadow" created by the heavy anchor and translates it back into the complex language of the original, wobbly boat.

By studying how this heavy anchor interacts with the boat, the team can figure out the boat's true shape without ever having to measure the boat directly in its chaotic state.

What Did They Actually Do?

The team, led by researchers from MIT, Taiwan, and Japan, used this "Heavy Anchor" method to calculate specific numbers called Mellin Moments.

• What is a Moment? Imagine the LCDA is a curve showing where the quarks are likely to be. A "moment" is like a summary statistic of that curve.
  • The 1st moment tells you the average position.
  • The 2nd moment tells you how wide or "spread out" the curve is.
  • The 3rd moment tells you if the curve is lopsided.
  • The 4th moment gives even finer details about the shape.

The Achievements:

1. The Kaon: They successfully calculated the first three non-trivial moments for the Kaon (a particle made of a strange quark and a down quark). This is like taking the first three detailed measurements of the Kaon's internal shape.
2. The Pion: They confirmed a previous result for the Pion (made of up and down quarks), specifically the 4th moment. This is a very high-level detail, like measuring the tiny ripples on the surface of the cloud.

Why Does This Matter?

Think of the Pion and Kaon as the "standard bricks" of the universe's matter. If we want to understand how they interact with light or other particles in high-energy experiments (like at the Large Hadron Collider), we need to know their exact internal shape.

• Old Way: Trying to guess the shape based on models or indirect measurements was like trying to guess the recipe of a cake by tasting the frosting.
• New Way (HOPE): This method allows physicists to calculate the recipe directly from the fundamental laws of physics (Quantum Chromodynamics) with high precision.

The Bottom Line

This paper is a progress report. The team has built a new, powerful tool (the Heavy Quark trick) and has successfully used it to take the first few precise measurements of the internal structure of pions and kaons.

They are currently refining their calculations (adding more data, fixing small errors) to get the final, perfect numbers. Once finished, this will help physicists better understand the fundamental building blocks of our universe, proving that even the most chaotic, invisible particles can be measured if you know the right trick to stabilize them.

Technical Summary

Here is a detailed technical summary of the paper "Parton physics from a heavy-quark operator product expansion: Dynamical lattice QCD calculation of moments of the pion and kaon light-cone distribution amplitudes."

1. Problem Statement

The Light-Cone Distribution Amplitude (LCDA), denoted as $\phi_M(\xi, \mu)$, is a fundamental non-perturbative quantity describing the momentum fraction distribution of valence quarks within a meson (such as the pion or kaon). It is essential for calculating exclusive scattering processes in Quantum Chromodynamics (QCD) via factorization theorems.

However, extracting the LCDA directly from Lattice QCD is challenging because:

• Euclidean vs. Minkowski: The definition of the LCDA involves a gauge-invariant operator with a light-like separation ($z^2=0$), which cannot be directly simulated in Euclidean lattice QCD.
• Operator Mixing: Traditional methods to compute Mellin moments (integrals of $\xi^n \phi_M$) rely on local twist-two operators. For higher moments ($n > 1$), these suffer from severe power-divergent mixing and operator proliferation, making calculations increasingly difficult and unreliable.
• Endpoint Sensitivity: Approaches like Large-Momentum Effective Theory (LaMET) struggle with finite momentum effects and power corrections near the endpoints of the distribution ($\xi \to \pm 1$).

The paper addresses the need for a robust, complementary method to determine the Mellin moments $\langle \xi^n \rangle_M$ of the pion and kaon LCDAs without relying on resolving the full $\xi$-dependence or suffering from the mixing issues of local operators.

2. Methodology: The HOPE Framework

The authors utilize the Heavy-Quark Operator Product Expansion (HOPE) method. This approach introduces a fictitious valence heavy quark ($\Psi$) with mass $m_\Psi$ to provide a hard scale, enabling a short-distance expansion.

Key Technical Steps:

1. Heavy-Light Currents: The method defines heavy-light axial currents $J_\mu^f(x) = \bar{\Psi}(x)\gamma_\mu\gamma_5 q_f(x) + \text{h.c.}$, where $q_f$ is a light or strange quark.
2. Euclidean Correlators: Instead of light-like separations, the calculation focuses on a Euclidean current-current correlator with space-like separation. The hadronic tensor is parametrized in terms of kinematic variables $Q^2 = -q^2$ and $\tilde{Q}^2 = Q^2 + m_\Psi^2$.
3. OPE Matching: The correlator is expanded using the Operator Product Expansion (OPE). The expansion coefficients (Wilson coefficients) are known at one-loop order in perturbation theory. The expansion is expressed in terms of Gegenbauer moments $\phi_{n,M}(\mu)$, which are directly related to the Mellin moments.
   $$\mathcal{V}_M(\tilde{\omega}, \tilde{Q}^2, \mu) = \sum_{n=0}^{\infty} F_n(\tilde{\omega}, \tilde{Q}^2, \mu, m_\Psi) \phi_{n,M}(\mu)$$
4. Lattice Implementation:
   • Ensembles: Calculations were performed on $N_f = 2+1$ dynamical gauge ensembles generated by the CLS collaboration (using $O(a)$-improved Wilson fermions).
   • Ratio Construction: To isolate the ground state and remove excited-state contamination, the authors construct ratios of three-point functions (involving the heavy-light currents) to two-point functions.
   • Excited-State Control: They employ an extrapolation strategy varying the Euclidean time extent ($\tau_e$) at fixed current separation ($\tau$) to extrapolate to the $\tau_e \to \infty$ limit, effectively removing excited-state contributions.
   • Channel Separation: By analyzing the real and imaginary parts of the correlator and combining data at $q$ and $-q$, they separate even and odd moments. The even-real channel is sensitive to the second moment, while odd channels provide sensitivity to the first and third moments.

3. Key Contributions

• Dynamical Kaon Calculation: This work presents the first dynamical ($N_f=2+1$) lattice QCD determination of the first three nontrivial Mellin moments ($\langle \xi \rangle, \langle \xi^2 \rangle, \langle \xi^3 \rangle$) for the kaon LCDA using the HOPE framework.
• Continuum Limit Benchmark for Pion: The paper summarizes a recently published result for the fourth Mellin moment of the pion LCDA in the quenched approximation ($N_f=0$), which represents the first continuum-limit determination of this specific moment using the HOPE method.
• Validation of HOPE: The study demonstrates the feasibility of the HOPE method for accessing higher moments (up to $n=4$) without the operator mixing problems associated with local operators.

4. Results

• Kaon Moments (Preliminary):
  • The authors performed fits to the lattice data for the first three nontrivial moments at a renormalization scale of $\mu = 2$ GeV.
  • Preliminary results show sensitivity to the first (odd), second (even), and third (odd) moments.
  • The analysis successfully separated even and odd contributions, confirming a non-zero signal for the odd channels (which vanish for the isospin-symmetric pion but are non-zero for the kaon due to $u/s$ quark mass differences).
  • Data was generated across multiple lattice spacings ($a \approx 0.049$ fm to $0.075$fm) and pion masses ($287$MeV to$428$ MeV) to allow for future continuum and physical mass extrapolations.
• Pion Fourth Moment (Published):
  • In the quenched approximation, the fourth moment was determined as:
    $$\langle \xi^4 \rangle_\pi = 0.039(28)(11)$$
    (First error: statistical + systematic; Second error: truncation of Wilson coefficients).
  • The second moment was found to be $\langle \xi^2 \rangle_\pi = 0.202(8)(9)$.
  • These results are consistent with previous HOPE calculations and other literature methods (LaMET, sum rules) after accounting for quenching and mass effects.

5. Significance

• Overcoming Lattice Limitations: The HOPE method offers a viable alternative to the "local operator" approach for higher moments, bypassing the power-divergent mixing that plagues standard lattice calculations for $n \ge 2$.
• Non-Perturbative Constraints: Providing precise values for the Gegenbauer moments of the kaon and pion is crucial for reducing uncertainties in QCD predictions for exclusive processes (e.g., $B \to \pi \ell \nu$, meson-photon transitions).
• Feasibility Demonstration: By successfully extracting the fourth moment of the pion and the first three of the kaon, the authors prove that the HOPE framework can be extended to higher orders and applied to flavor-asymmetric mesons (kaons) in full QCD.
• Future Outlook: The work lays the groundwork for a complete, continuum-limit determination of the kaon LCDA moments in full QCD, which will significantly improve the precision of Standard Model tests involving heavy-flavor decays.

In summary, this paper validates the HOPE method as a powerful tool for extracting parton distribution information from lattice QCD, successfully applying it to the challenging case of the kaon in a dynamical setting and providing a benchmark continuum result for the pion.