Kinematic enhancement for nucleon interpolators

Motivated by future Electron-Ion Collider physics, this paper demonstrates that kinematically enhanced interpolators significantly improve the precision of renormalized nucleon matrix elements at high momenta while showing no dependence on lattice spacing, thereby establishing them as a promising standard for modern lattice QCD parton physics calculations.

The Gist

Imagine you are trying to take a photograph of a hummingbird in flight. If you use a standard camera with a slow shutter speed, the bird will look like a blurry mess. To get a sharp picture, you need a very fast shutter speed and a lot of light. In the world of particle physics, scientists are trying to take "photos" of protons (nucleons) to understand what they are made of. But instead of light, they use complex math simulations on supercomputers, and instead of a hummingbird, they are looking at particles moving at incredibly high speeds.

Here is the simple story of what this paper does, using everyday analogies.

The Problem: The "Blurry Photo" of Fast Particles

Scientists use a method called Lattice QCD (Quantum Chromodynamics) to simulate how particles like protons behave. To understand how protons are built from smaller parts called "quarks" (which is crucial for future particle colliders), they need to simulate protons moving very fast.

However, there is a major problem: The Signal-to-Noise Ratio.

• The Signal: The actual data about the fast-moving proton.
• The Noise: Random mathematical "static" that gets louder and louder as the proton moves faster.

Think of it like trying to hear a whisper (the signal) in a room where a jet engine is revving up (the noise). As the proton gets faster, the jet engine gets louder, and the whisper becomes impossible to hear. This makes it very hard to get accurate results for fast-moving protons.

The Solution: A "Kinematic Booster"

The authors of this paper tested a new tool, which they call "kinematically enhanced interpolators."

Imagine you are trying to catch a specific type of fish in a river.

• The Old Way: You use a standard net that catches everything—fish, leaves, rocks, and mud. You have to sift through a huge pile of junk to find the one fish you want. The more water (momentum) flows, the more junk you catch, making it harder to find your fish.
• The New Way: The authors designed a "smart net" that is shaped exactly like the fish they are looking for. It only catches the fish and lets the leaves and rocks pass right through.

In physics terms, they changed the mathematical "net" (the interpolator) used to create the proton in the simulation. By tuning this net to match the specific shape of a fast-moving proton, they filtered out the "junk" (noise) before it even started.

What They Found

The team ran these simulations on three different supercomputing setups (called "ensembles") to make sure their results were real and not just a fluke. Here is what happened:

1. A Massive Boost in Clarity: When they used the new "smart net," the quality of their data improved by ten times (an order of magnitude). It's like going from a grainy, black-and-white photo to a crystal-clear, high-definition 4K image.
2. No New Distortions: Sometimes, when you fix one problem, you create another. They worried that this new method might introduce "excited state contamination" (a fancy way of saying the simulation might get confused about which proton state it is looking at). They checked this carefully and found no new confusion. The new method is just as clean as the old one, but much sharper.
3. Consistency Across Scales: They tested this on three different "grid sizes" (lattice spacings). Even though the grids were different, the results were the same. This proves the method is robust and reliable, not just a trick that works on one specific setting.

The "Secret Sauce": The Gamma-Plus Trick

The paper highlights a specific mathematical trick they used, involving a symbol called $\gamma_+$.
Think of this as a special filter that cuts the work in half.

• Normally, the computer has to calculate information in all directions (up, down, left, right, forward, backward).
• The $\gamma_+$ filter realizes that for a fast-moving proton, only the "forward" information matters. It tells the computer, "Ignore everything else."
• This not only makes the data cleaner but also cuts the computing time and cost in half because the computer doesn't have to do unnecessary calculations.

The Bottom Line

This paper proves that by using these new, smarter mathematical "nets," scientists can finally get clear, high-quality pictures of fast-moving protons without needing to wait for even bigger supercomputers.

This is a big deal because it opens the door to studying the internal structure of protons with much higher precision. This is essential for understanding the physics that future particle colliders (like the Electron-Ion Collider) will explore. The authors conclude that this method should become a standard tool for anyone doing this kind of high-speed particle physics.

Technical Summary

Technical Summary: Kinematic enhancement for nucleon interpolators

Problem Statement
The reliable extraction of partonic structure information from lattice QCD simulations, particularly within the Large Momentum Effective Theory (LaMET) framework, requires the study of hadrons at high momenta ($P_z$). However, as the momentum of the hadron increases, the signal-to-noise ratio (SNR) of hadronic observables deteriorates rapidly. For the nucleon, this precision bottleneck typically limits accessible momenta to $P_z \lesssim 3$ GeV on contemporary ensembles, hindering the approach to the infinite momentum limit required for reliable matching to light-cone parton distribution functions (PDFs). While kinematically enhanced interpolators were previously proposed to improve SNRs for mesons and two-point functions, their efficacy for nucleon quasi-distribution matrix elements and their impact on excited-state contamination (ESC) and discretization effects had not been systematically benchmarked.

Methodology
The authors investigate the application of kinematically enhanced interpolating operators to the calculation of unpolarized isovector nucleon quark matrix elements. The study utilizes three CLS ensembles (H102, S400, N203) with $N_f=2+1$ O($a$) improved Sheikholeslami-Wohlert fermions, covering lattice spacings $a \in [0.064, 0.085]$ fm and a pion mass of $m_\pi \approx 350$ MeV.

The core methodology involves comparing three distinct choices of interpolator kernels $(T, \Gamma)$:

1. Conventional: $(P_+, 1)$, using standard parity projection.
2. $\gamma_t$ enhanced: $(\gamma_t, \gamma_t)$.
3. $\gamma_+$ enhanced: $(\gamma_+, \gamma_+)$.

The kinematic enhancement relies on the fact that at large momentum, the "+" component of the quark field dominates. By constructing interpolators with projectors $T$ and diquark kernels $\Gamma$ that select these "+" components (specifically using $T=\gamma_+$ and $\Gamma=\gamma_+$), the overlap with the ground state is kinematically enhanced by a factor proportional to $P_+^{n/2}$, where $n$ is the number of valence quarks. Crucially, the $\gamma_+$ projector eliminates "−" components of the propagators, reducing the computational cost of inversions by half and providing an additional variance reduction factor.

The authors compute two-point and three-point correlation functions for nucleon momenta $P_z \approx 2.5$ GeV. They employ a two-state fit to extract ground-state matrix elements, thereby controlling excited-state contamination. The matrix elements are renormalized using a hybrid scheme (ratio scheme for short distances, self-renormalization for long distances) to remove linear and logarithmic divergences associated with the Wilson line.

Key Contributions and Results

• Variance Reduction: The study demonstrates that kinematically enhanced interpolators, particularly the $(\gamma_+, \gamma_+)$ configuration, improve the statistical precision of nucleon matrix elements by approximately one order of magnitude at $P_z \sim 2.5$ GeV compared to conventional interpolators. This improvement is observed across all three lattice ensembles.
• Excited-State Contamination: Contrary to concerns that kinematic enhancement might increase overlap with excited states, the two-state fits show no statistically significant increase in ESC when using enhanced interpolators. The curvature of the data remains consistent with the conventional case, indicating that the ground-state matrix elements are reliably extracted.
• Discretization Effects: By comparing results at nearly identical momenta across three different lattice spacings, the authors observe no statistically significant dependence on the lattice spacing $a$ in the renormalized matrix elements. This suggests that discretization effects remain small even at momenta as high as $P_z \sim 2.5$ GeV.
• Computational Efficiency: The use of the $\gamma_+$ projector reduces the number of required inversions by half because "−" components of the propagators cancel out in the contraction with the sink. This yields a factor-of-four gain in variance (two from combining $\gamma_t$ and $\gamma_z$ information, and two from eliminating "−" components).

Significance
The paper concludes that kinematically improved operators are highly advantageous for parton physics calculations involving boosted nucleons. The results suggest that:

1. The statistical reach of LaMET calculations can be significantly extended, potentially allowing for simulations at higher momenta ($P_z \gtrsim 3$ GeV) where power corrections are smaller.
2. The method is robust, as it does not introduce new systematic uncertainties regarding excited states or discretization effects within the tested range.
3. The technique is a promising candidate for becoming a standard component of modern lattice QCD calculations for baryon observables, supporting the physics programs of future facilities like the Electron-Ion Collider (EIC) and the Electron-Ion collider in China (EIcC).

The authors note that while the current results are encouraging, further work is needed to address $O((aP_z)^n)$ discretization effects at even higher momenta, potentially by combining these interpolators with methods like the Coulomb-gauge approach or using finer lattices.