Extracting the distribution amplitude of pseudoscalar mesons using the HOPE method

The Big Picture

Imagine a pion, one of the most common particles in the universe, zipping along at nearly the speed of light. Inside it, a quark and an antiquark are locked together by the strong nuclear force, sharing momentum in some proportion. If you could freeze that particle mid-flight and ask “who’s carrying how much of the momentum right now?”, the answer would follow a probability distribution called the light-cone distribution amplitude, or LCDA. That distribution is an essential ingredient in calculating the outcomes of high-energy collisions at facilities like the Large Hadron Collider.

The problem: this distribution is fiendishly hard to compute from first principles. The strong force, described by quantum chromodynamics (QCD), behaves very differently at different energy scales. At high energies, step-by-step approximation techniques work well. At the lower energies governing a pion’s internal structure, they break down entirely.

Physicists have long relied on lattice QCD, which discretizes spacetime onto a grid for numerical computation on supercomputers. But there’s a fundamental mismatch. LCDAs are defined along a light-like direction, a path traveling at the speed of light, and lattice calculations work in a mathematical framework where such paths don’t exist.

The HOPE Collaboration, spanning MIT, Fermilab, RIKEN, Argonne, and institutions in Taiwan and Spain, has developed a workaround. They report progress computing the low Mellin moments (statistical averages that quantify the shape of the distribution) of the LCDAs for pions and kaons, using the Heavy-Quark Operator Product Expansion (HOPE) method.

Key Insight: By introducing a fictitious heavy quark as a mathematical probe, the HOPE method converts an incalculable light-like quantity into ordinary lattice QCD correlation functions.

How It Works

The HOPE method sidesteps the light-like operator entirely. Instead of computing the LCDA directly, researchers calculate a quantity involving two axial-vector currents built from light quarks and a fictional heavy quark that doesn’t exist in nature. This fictitious particle is a mathematical scaffold, not a real constituent. It lets the team probe a meson’s internal structure using currents that can be placed at standard lattice spacetime separations.

The connection back to the LCDA comes through the Operator Product Expansion (OPE), a technique that simplifies the product of two field operators when they are close together in spacetime, expressing it as a sum of simpler terms with calculable coefficients. When the heavy quark is sufficiently massive, the two-current matrix element (a quantity encoding the meson’s internal structure as probed by those currents) can be written as a sum over Gegenbauer moments, a mathematical decomposition of the LCDA’s shape, weighted by perturbatively calculable coefficients.

In practice, the calculation proceeds in three steps:

1. Compute correlation functions on the lattice. The team calculates two-point and three-point correlation functions using lattice QCD, employing the Generalized Eigenvalue Problem (GEVP) to construct an optimized operator that suppresses excited-state contamination. The GEVP is a linear-algebra technique for disentangling contributions from multiple quantum states.

2. Extract the hadronic matrix element. By forming a ratio of three-point to two-point correlators at large time separations, the researchers isolate the hadronic tensor $R_M^{\mu\nu}(t^-, p, q)$ as a function of momentum and current separation. This structured object encodes how the meson responds to the two current insertions.

3. Fit to the HOPE formula. The extracted matrix element is fit to the one-loop HOPE expression, which predicts its dependence on kinematic variables in terms of the Gegenbauer moments $\phi_{M,n}$.

The fitting procedure uses both the $t^-$-even and $t^-$-odd components of the hadronic matrix element, two independent observables providing complementary sensitivity to the distribution’s moments. These Gegenbauer moments convert to Mellin moments $\langle\xi^n\rangle_M$, representing the average value of $\xi^n$ (where $\xi$ is the momentum fraction imbalance between quark and antiquark) weighted by the LCDA.

The team applied this procedure to both the pion and the kaon using lattice ensembles with physical or near-physical quark masses. The kaon is particularly interesting because it contains a strange quark paired with an up or down antiquark, creating a flavor asymmetry. Its LCDA is not symmetric around $\xi = 0$, so the odd Mellin moments are nonzero. This is a direct signature of SU(3) flavor-symmetry breaking with real phenomenological consequences.

Systematic uncertainties are carefully tracked throughout, including higher-twist contamination (corrections from subleading OPE terms), renormalization scale dependence (sensitivity to the energy scale at which coupling constants are defined), and discretization artifacts from the finite lattice spacing.

Why It Matters

LCDAs enter directly into predictions for exclusive processes, precisely measured reactions where specific particles go in and specific particles come out. Examples include the pion electromagnetic form factor, B-meson decays, and transition form factors relevant to Standard Model tests. Many tensions between theory and experiment in flavor physics hinge on the precise values of these non-perturbative quantities. Without accurate distribution amplitudes, theoretical predictions lose interpretive power, and experimental signals can’t be confidently read as evidence for new physics.

The HOPE method is one of several recently developed approaches to circumventing the light-cone problem in lattice QCD, a family that includes quasi-LCDAs, pseudo-LCDAs, and the Compton amplitude method. What sets HOPE apart is its systematic use of heavy-quark expansion technology, providing a rigorous theoretical framework with clearly defined power corrections. As lattice simulations improve and computing resources scale up, methods like HOPE will supply the non-perturbative inputs that high-energy experiments need.

Bottom Line: The HOPE Collaboration has a working, systematically controlled method for extracting pion and kaon distribution amplitudes from lattice QCD, bringing precision QCD predictions for exclusive processes closer to reality.