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. You need that distribution to calculate 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, perturbative 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 use a mathematical framework where such paths don’t exist.

The HOPE Collaboration, with members at 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). The OPE 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 can be written as a sum over Gegenbauer moments, which decompose the LCDA’s shape, each weighted by a perturbatively calculable coefficient.

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. They use the Generalized Eigenvalue Problem (GEVP), a linear-algebra technique for disentangling contributions from multiple quantum states, to construct an optimized operator that suppresses excited-state contamination.

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 tensor 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. These are two independent observables with complementary sensitivity to the distribution’s moments. The Gegenbauer moments convert to Mellin moments $\langle\xi^n\rangle_M$: 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 makes a good test case because it contains a strange quark paired with an up or down antiquark, producing a flavor asymmetry. Its LCDA is not symmetric around $\xi = 0$, so the odd Mellin moments are nonzero. That asymmetry is a direct signature of SU(3) flavor-symmetry breaking, and it has real phenomenological consequences.

Systematic uncertainties are tracked throughout: 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, reactions where specific particles go in and specific particles come out. The pion electromagnetic form factor, B-meson decays, and transition form factors relevant to Standard Model tests all depend on them. 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 their discriminating power, and experimental signals can’t be confidently read as evidence for new physics.

The HOPE method is one of several recent approaches to the light-cone problem in lattice QCD, alongside quasi-LCDAs, pseudo-LCDAs, and the Compton amplitude method. HOPE is distinguished by its systematic use of heavy-quark expansion technology, which gives it clearly defined power corrections and a controlled approximation scheme. As lattice simulations improve and computing resources grow, 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, pushing precision QCD predictions for exclusive processes closer to reality.