Kinematically Enhanced Lattice Interpolators

• The paper demonstrates that tailoring lattice operators with momentum-dependent Dirac structures significantly improves overlap with boosted hadrons in lattice QCD.
• Enhanced interpolators yield marked signal-to-noise ratio gains, with improvements up to 40–50 times in boosted pion and nucleon studies.
• This approach extends beyond QCD, informing affine, helicity, and geometric lattice constructions in diverse interpolation frameworks.

Kinematically enhanced lattice interpolation operators are lattice-based operators whose defining structure is adapted to the kinematics or geometry of the target problem rather than being only local, pointlike, or isotropic by construction. In contemporary lattice QCD, the term refers most directly to boosted-hadron interpolators whose Dirac structure is chosen so that the overlap with a fast pion or nucleon carries explicit momentum-dependent factors, producing larger ground-state overlap at large momentum and, because the same enhancement does not apply to the dominant variance channel, improved signal-to-noise ratios (Zhang et al., 1 Jan 2025). A broader, inferential usage includes lattice constructions that preserve affine deformations exactly, encode nonzero-momentum helicity and little-group structure, or engineer spatial profiles through auxiliary three-dimensional fields; in each case, the operator is “enhanced” by embedding physically or geometrically relevant kinematic information into the interpolation or source construction (Zhao, 29 Aug 2025).

1. Terminology and conceptual scope

In the narrow sense established explicitly in recent lattice-QCD work, a kinematically enhanced interpolating operator is a hadron source or sink for which the overlap factor with a boosted state grows with the hadron momentum because the operator has been aligned with the dominant light-cone or Lorentz structure of that state (Zhang et al., 1 Jan 2025). This is distinct from purely spatial source optimization: the enhancement comes from the operator’s spinor or tensor structure, not only from its radial profile.

A broader interpretation is useful but should be treated as inferential. Several earlier literatures do not use the exact phrase “kinematically enhanced,” yet they construct operators that incorporate explicit kinematic consistency conditions. In lattice QCD at nonzero momentum, helicity-based operators encode the reduced symmetry of a cubic box through little-group subduction rather than naive rest-frame spin labels (Thomas et al., 2011). In lattice spectroscopy with three-dimensional auxiliary fermions, the source width is controlled by a tunable 3D mass, and the resulting correlators acquire markedly softer short-distance behavior (Papinutto et al., 2018). In atomistic-to-continuum analysis, Bravais-lattice interpolants are designed to reproduce affine fields exactly, preserve strain norms, and supply smooth nodal reconstructions (Ortner et al., 2012). This suggests that the phrase can function as an umbrella label for lattice operators whose definition is constrained by the relevant kinematic, symmetry, or deformation structure.

The main conceptual divide is therefore between generic interpolation and structure-adapted interpolation. Generic interpolation reconstructs or excites a target from nearby data. Kinematically enhanced constructions add a second requirement: the operator should respect the momentum, helicity, affine, or symmetry content that dominates the target state or field.

2. Boosted-hadron operators in lattice QCD

The clearest explicit realization appears in boosted-hadron lattice QCD. For momentum chosen along the $z$-direction, the quark field is decomposed into light-cone components,

$$\psi=\psi_+ + \psi_-, \qquad \psi_\pm = \frac{1}{2}\gamma_\mp\gamma_\pm \psi,$$

with

$$\gamma_\pm \equiv \frac{1}{\sqrt{2}}(\gamma_t \pm i\gamma_z).$$

The operative idea is that, in the large-momentum limit, the leading dynamics is carried by the plus component $\psi_+$, so interpolators built from that structure should couple more strongly to highly boosted hadrons (Zhang et al., 1 Jan 2025).

For the pion, the standard pseudoscalar interpolator

$O_\pi^{\rm std}(x)=\bar u(x)\gamma_5 d(x)$

is replaced by axial-vector operators

$$O_{\pi,\mu}(x)=\bar u(x)\gamma_\mu\gamma_5 d(x), \qquad \mu=t,z,$$

and especially by the light-cone combination

$O_{\pi,+}(x)=\bar u(x)\gamma_+\gamma_5 d(x).$

The key overlap relation is the axial-current matrix element

$$\langle \Omega | \bar u\gamma_\mu\gamma_5 d | \pi(\vec P)\rangle = i f_\pi P_\mu.$$

Because the overlap contains $P_\mu$, the two-point signal is enhanced parametrically like $\mathcal O(P_\mu^2/m_\pi^2)$, or in the temporal channel like $$\psi=\psi_+ + \psi_-, \qquad \psi_\pm = \frac{1}{2}\gamma_\mp\gamma_\pm \psi,$$0, relative to a conventional pseudoscalar construction (Zhang et al., 1 Jan 2025).

For the nucleon, the standard interpolator

$$\psi=\psi_+ + \psi_-, \qquad \psi_\pm = \frac{1}{2}\gamma_\mp\gamma_\pm \psi,$$1

with

$$\psi=\psi_+ + \psi_-, \qquad \psi_\pm = \frac{1}{2}\gamma_\mp\gamma_\pm \psi,$$2

is supplemented by boosted-state-optimized operators such as

$$\psi=\psi_+ + \psi_-, \qquad \psi_\pm = \frac{1}{2}\gamma_\mp\gamma_\pm \psi,$$3

and

$$\psi=\psi_+ + \psi_-, \qquad \psi_\pm = \frac{1}{2}\gamma_\mp\gamma_\pm \psi,$$4

Their overlaps are parameterized by

$$\psi=\psi_+ + \psi_-, \qquad \psi_\pm = \frac{1}{2}\gamma_\mp\gamma_\pm \psi,$$5

and

$$\psi=\psi_+ + \psi_-, \qquad \psi_\pm = \frac{1}{2}\gamma_\mp\gamma_\pm \psi,$$6

By contrast, the standard $$\psi=\psi_+ + \psi_-, \qquad \psi_\pm = \frac{1}{2}\gamma_\mp\gamma_\pm \psi,$$7 operator has no explicit $$\psi=\psi_+ + \psi_-, \qquad \psi_\pm = \frac{1}{2}\gamma_\mp\gamma_\pm \psi,$$8 term in its leading overlap decomposition. The nucleon signal therefore acquires an enhancement of order $$\psi=\psi_+ + \psi_-, \qquad \psi_\pm = \frac{1}{2}\gamma_\mp\gamma_\pm \psi,$$9 (Zhang et al., 1 Jan 2025).

This boosted-hadron construction is complementary to momentum smearing rather than a replacement for it. Momentum smearing adjusts the spatial wave packet. Kinematic enhancement changes the Dirac structure so that the same boosted hadron is created with a larger overlap factor.

3. Signal-to-noise mechanism and empirical status

The underlying spectral logic is standard: $$\gamma_\pm \equiv \frac{1}{\sqrt{2}}(\gamma_t \pm i\gamma_z).$$0 where $$\gamma_\pm \equiv \frac{1}{\sqrt{2}}(\gamma_t \pm i\gamma_z).$$1. Kinematic enhancement acts by increasing $$\gamma_\pm \equiv \frac{1}{\sqrt{2}}(\gamma_t \pm i\gamma_z).$$2 at large momentum. Its distinctive claim is stronger than a mere rescaling of the correlator: the dominant variance channel is governed by different states, so the same boost factors do not propagate into the noise (Zhang et al., 1 Jan 2025).

For the pion,

$$\gamma_\pm \equiv \frac{1}{\sqrt{2}}(\gamma_t \pm i\gamma_z).$$3

The late-time variance is dominated by two pions at rest, so

$$\gamma_\pm \equiv \frac{1}{\sqrt{2}}(\gamma_t \pm i\gamma_z).$$4

while the signal decays as $$\gamma_\pm \equiv \frac{1}{\sqrt{2}}(\gamma_t \pm i\gamma_z).$$5. The standard asymptotic behavior is therefore

$$\gamma_\pm \equiv \frac{1}{\sqrt{2}}(\gamma_t \pm i\gamma_z).$$6

Because the noise-dominating state is not boosted, the $$\gamma_\pm \equiv \frac{1}{\sqrt{2}}(\gamma_t \pm i\gamma_z).$$7-dependent enhancement affects the signal but not the leading variance, and the signal-to-noise ratio inherits the $$\gamma_\pm \equiv \frac{1}{\sqrt{2}}(\gamma_t \pm i\gamma_z).$$8 gain. For nucleons the corresponding Parisi–Lepage channel is dominated by three zero-momentum pions, leading to the analogous $$\gamma_\pm \equiv \frac{1}{\sqrt{2}}(\gamma_t \pm i\gamma_z).$$9 expectation (Zhang et al., 1 Jan 2025).

Proof-of-principle calculations were performed on a $\psi_+$0 ensemble with $\psi_+$1, a valence pion mass near $\psi_+$2, momentum smearing with $\psi_+$3, pion momenta from $\psi_+$4 to $\psi_+$5, and nucleon momenta up to about $\psi_+$6 (Zhang et al., 1 Jan 2025). At $\psi_+$7, the standard $\psi_+$8 correlator becomes noise-dominated beyond about $\psi_+$9, whereas the enhanced pion operators remain usable to about $O_\pi^{\rm std}(x)=\bar u(x)\gamma_5 d(x)$0. The observed ordering at large momentum is

$O_\pi^{\rm std}(x)=\bar u(x)\gamma_5 d(x)$1

The Lanczos ground-state analysis gives pion SNR improvements of about $O_\pi^{\rm std}(x)=\bar u(x)\gamma_5 d(x)$2 for $O_\pi^{\rm std}(x)=\bar u(x)\gamma_5 d(x)$3, and the measured gain reaches $O_\pi^{\rm std}(x)=\bar u(x)\gamma_5 d(x)$4 at $O_\pi^{\rm std}(x)=\bar u(x)\gamma_5 d(x)$5, corresponding to an effective statistics increase of order $O_\pi^{\rm std}(x)=\bar u(x)\gamma_5 d(x)$6 (Zhang et al., 1 Jan 2025). For nucleons, the enhanced operators $O_\pi^{\rm std}(x)=\bar u(x)\gamma_5 d(x)$7 show visibly better precision than the standard $O_\pi^{\rm std}(x)=\bar u(x)\gamma_5 d(x)$8, with Lanczos SNR improvement in the ground-state energy of roughly $O_\pi^{\rm std}(x)=\bar u(x)\gamma_5 d(x)$9. Three-point tests for the bare unpolarized quasi-PDF matrix element at $$O_{\pi,\mu}(x)=\bar u(x)\gamma_\mu\gamma_5 d(x), \qquad \mu=t,z,$$0 also show improved correlator-level precision at $$O_{\pi,\mu}(x)=\bar u(x)\gamma_\mu\gamma_5 d(x), \qquad \mu=t,z,$$1, particularly for $$O_{\pi,\mu}(x)=\bar u(x)\gamma_\mu\gamma_5 d(x), \qquad \mu=t,z,$$2 and $$O_{\pi,\mu}(x)=\bar u(x)\gamma_\mu\gamma_5 d(x), \qquad \mu=t,z,$$3 (Zhang et al., 1 Jan 2025).

The method is not uniformly superior. At rest, the standard pion pseudoscalar operator remains best for the ground state, and enhanced pion operators can couple more strongly to heavier states. The construction is therefore fundamentally a large-boost technique.

4. Earlier lattice-QCD lineages

Kinematic enhancement in the boosted-hadron sense has important antecedents in lattice-QCD operator design. One lineage arises from helicity operators for mesons in flight. There the natural organizing principle at nonzero momentum is not the rest-frame spin projection $$O_{\pi,\mu}(x)=\bar u(x)\gamma_\mu\gamma_5 d(x), \qquad \mu=t,z,$$4 but helicity $$O_{\pi,\mu}(x)=\bar u(x)\gamma_\mu\gamma_5 d(x), \qquad \mu=t,z,$$5, followed by subduction into irreducible representations of the relevant little group. Continuum-like operators $$O_{\pi,\mu}(x)=\bar u(x)\gamma_\mu\gamma_5 d(x), \qquad \mu=t,z,$$6 are rotated to definite-helicity operators,

$$O_{\pi,\mu}(x)=\bar u(x)\gamma_\mu\gamma_5 d(x), \qquad \mu=t,z,$$7

and then subduced via

$$O_{\pi,\mu}(x)=\bar u(x)\gamma_\mu\gamma_5 d(x), \qquad \mu=t,z,$$8

with $$O_{\pi,\mu}(x)=\bar u(x)\gamma_\mu\gamma_5 d(x), \qquad \mu=t,z,$$9. The resulting correlator matrices are approximately block-diagonal in helicity magnitude, and the variationally optimized operators give rapid relaxation in Euclidean time to the targeted state (Thomas et al., 2011). Although the phrase “kinematically enhanced” is not used there, the construction is explicitly momentum-resolved, helicity-resolved, and symmetry-adapted.

A second lineage comes from extended interpolating fields built from quenched three-dimensional fermions. In this approach, auxiliary 3D fermion fields live on a spatial time slice and are coupled to physical 4D quarks through pseudoscalar bilinears such as $O_{\pi,+}(x)=\bar u(x)\gamma_+\gamma_5 d(x).$0. The lattice action is

$O_{\pi,+}(x)=\bar u(x)\gamma_+\gamma_5 d(x).$1

and the extended baryon operator has the form

$O_{\pi,+}(x)=\bar u(x)\gamma_+\gamma_5 d(x).$2

The corresponding extended quark propagator is a 3D–4D–3D sandwich,

$O_{\pi,+}(x)=\bar u(x)\gamma_+\gamma_5 d(x).$3

Here the tunable 3D mass $O_{\pi,+}(x)=\bar u(x)\gamma_+\gamma_5 d(x).$4 controls the source width; lighter 3D fermions generate broader sources, and the broadened source exhibits an exponential long-distance decay rather than the Gaussian-like profile of Jacobi smearing (Papinutto et al., 2018).

This 3D-fermion construction carries two properties that are highly relevant to the broader idea of kinematic enhancement. First, it gives a field-theoretic control parameter for source extent and overlap. Second, it changes the ultraviolet structure of the correlator itself: for baryons, the short-distance divergence softens from

$O_{\pi,+}(x)=\bar u(x)\gamma_+\gamma_5 d(x).$5

for a local operator to

$O_{\pi,+}(x)=\bar u(x)\gamma_+\gamma_5 d(x).$6

for the 3D-extended operator, while an extended meson correlator becomes finite at short distance (Papinutto et al., 2018). The numerical studies on CLS $O_{\pi,+}(x)=\bar u(x)\gamma_+\gamma_5 d(x).$7 configurations show much better short-distance regularity than point or Jacobi sources, similar noise to Jacobi smearing, and a source family that is potentially useful in GEVP bases (Scardino et al., 2016).

5. Generalized meanings in lattice interpolation theory

Outside lattice QCD, the same phrase can be generalized only inferentially, but several mathematically precise constructions fit the pattern of embedding kinematic structure into lattice interpolation.

In Bravais-lattice interpolation for atomistic and atomistic/continuum multiscale methods, the foundational requirement is exact affine reproduction: $O_{\pi,+}(x)=\bar u(x)\gamma_+\gamma_5 d(x).$8 Starting from the nodal interpolant

$O_{\pi,+}(x)=\bar u(x)\gamma_+\gamma_5 d(x).$9

the theory constructs the smoothed quasi-interpolant

$$\langle \Omega | \bar u\gamma_\mu\gamma_5 d | \pi(\vec P)\rangle = i f_\pi P_\mu.$$0

and the smooth nodal interpolant

$$\langle \Omega | \bar u\gamma_\mu\gamma_5 d | \pi(\vec P)\rangle = i f_\pi P_\mu.$$1

The exact preservation of affine fields is the decisive kinematic consistency condition, while $$\langle \Omega | \bar u\gamma_\mu\gamma_5 d | \pi(\vec P)\rangle = i f_\pi P_\mu.$$2 combines nodal exactness with additional smoothness, and $$\langle \Omega | \bar u\gamma_\mu\gamma_5 d | \pi(\vec P)\rangle = i f_\pi P_\mu.$$3 supplies a local quasi-interpolant with approximation estimates up to $$\langle \Omega | \bar u\gamma_\mu\gamma_5 d | \pi(\vec P)\rangle = i f_\pi P_\mu.$$4 without requiring pointwise nodal traces (Ortner et al., 2012). In that literature, “enhancement” refers not to large momentum but to exact low-order deformation kinematics, smoothing, and stable discrete–continuum norm equivalences.

A related but distinct nonlinear notion appears in the theory of lattice Lipschitz operators. There a map is diagonal with respect to a basis $$\langle \Omega | \bar u\gamma_\mu\gamma_5 d | \pi(\vec P)\rangle = i f_\pi P_\mu.$$5 if

$$\langle \Omega | \bar u\gamma_\mu\gamma_5 d | \pi(\vec P)\rangle = i f_\pi P_\mu.$$6

and it is lattice Lipschitz when each output coordinate is controlled coordinatewise by the corresponding input coordinate. For almost diagonal maps, coordinatewise McShane and Whitney extensions are modified to

$$\langle \Omega | \bar u\gamma_\mu\gamma_5 d | \pi(\vec P)\rangle = i f_\pi P_\mu.$$7

$$\langle \Omega | \bar u\gamma_\mu\gamma_5 d | \pi(\vec P)\rangle = i f_\pi P_\mu.$$8

after first estimating approximate eigenvector directions. This is again a structure-adapted interpolation theory, now enhanced by basis selection and approximate diagonal dynamics rather than by momentum (Arnau et al., 2023).

In zonotopal interpolation, another exact lattice-delta construction appears. For a totally unimodular vector configuration $$\langle \Omega | \bar u\gamma_\mu\gamma_5 d | \pi(\vec P)\rangle = i f_\pi P_\mu.$$9, the projected Todd operator

$P_\mu$0

acts on the box spline $P_\mu$1 so that

$P_\mu$2

for each interior lattice point $P_\mu$3. The operators $P_\mu$4 form a canonical interpolation basis attached to the lattice geometry of the zonotope (Lenz, 2013). This is not “kinematic” in the QCD sense, but it is a paradigmatic example of interpolation enriched by combinatorial and geometric structure.

Domain	Representative construction	Structural feature
Atomistic/continuum	$P_\mu$5, $P_\mu$6 (Ortner et al., 2012)	affine reproduction, smooth nodal exactness
Almost diagonal maps	lattice McShane/Whitney extensions (Arnau et al., 2023)	basis- and eigenvector-adapted coordinatewise interpolation
Zonotopal interpolation	$P_\mu$7 (Lenz, 2013)	canonical delta interpolation on interior lattice points

6. Algorithmic realizations and broader analogues

A recurring practical issue is that kinematically or geometrically enriched lattice operators are often more expensive than simple nearest-neighbor or tensor-product schemes. One response is automatic code generation for spline evaluation on general lattices. A general framework based on the convolution sum

$P_\mu$8

analyzes lattice cosets, support regions, piecewise-polynomial structure, and symmetries to derive fast evaluation schemes for box splines and Voronoi splines on BCC, FCC, Cartesian, and even $P_\mu$9 lattices (Horacsek et al., 2021). This does not introduce the phrase “kinematically enhanced,” but it supplies an implementation template for sophisticated lattice interpolants whose quality depends on non-Cartesian geometry and support structure.

A second analogue arises in scalable Gaussian processes. There the standard SKI factorization

$\mathcal O(P_\mu^2/m_\pi^2)$0

is reinterpreted on the permutohedral lattice, with $\mathcal O(P_\mu^2/m_\pi^2)$1 given by barycentric interpolation inside the enclosing simplex. Each input has only $\mathcal O(P_\mu^2/m_\pi^2)$2 neighbors rather than the $\mathcal O(P_\mu^2/m_\pi^2)$3 neighbors of a cubic lattice, and the inducing lattice points are stored sparsely (Kapoor et al., 2021). This suggests a different type of enhancement: the operator is adapted to lattice geometry in a way that improves high-dimensional scaling rather than physical-state overlap.

A third analogue comes from discrete hydrodynamics. Using discrete velocities $\mathcal O(P_\mu^2/m_\pi^2)$4 and weights $\mathcal O(P_\mu^2/m_\pi^2)$5 satisfying isotropic moment conditions, one obtains lattice differential operators such as

$\mathcal O(P_\mu^2/m_\pi^2)$6

and

$\mathcal O(P_\mu^2/m_\pi^2)$7

These are not interpolation operators strictly speaking, but they are derived from lattice kinematics and built-in isotropy, and they illustrate how discrete velocity structure can organize operator design (Ramadugu et al., 2012).

7. Limitations, misconceptions, and open questions

A common misconception is to treat kinematic enhancement as a synonym for generic smearing. In the explicit boosted-hadron literature, that is incorrect: the main novelty is the momentum-dependent Dirac structure, while momentum smearing remains in use and is explicitly complementary (Zhang et al., 1 Jan 2025). Likewise, the earlier 3D-fermion source construction is not merely another heuristic smearing recipe; its claims concern renormalization, tunable source width, and improved ultraviolet behavior (Papinutto et al., 2018).

A second misconception is that enhancement at large momentum automatically solves state isolation. The precision review of LaMET is explicit that the maximum attainable momentum remains a limiting factor, finer lattices remain desirable, and excited-state contamination still requires separate treatment through interpolator optimization, GEVP methods, or Lanczos-based extraction (Zhao, 29 Aug 2025). The new boosted-hadron operators improve the overlap and the signal-to-noise ratio, but they do not remove discretization effects, finite-volume effects, or excited-state systematics.

The earlier lattice-QCD constructions also retain open questions. The 3D-fermion extended operators were tested only on an exploratory $\mathcal O(P_\mu^2/m_\pi^2)$8, $\mathcal O(P_\mu^2/m_\pi^2)$9-symmetric ensemble at relatively heavy pion mass, and their authors explicitly stated that a satisfactory understanding of the interplay between short-distance behavior and excited states requires a scaling study (Scardino et al., 2016). They did not present boosted-hadron, nonzero-momentum, or matrix-element applications. Helicity-subduced operators for mesons in flight establish a clean group-theoretic framework, but they also show that reduced rotational symmetry and the use of spatial rather than four-vector operator structures complicate continuum spin identification at finite momentum (Thomas et al., 2011).

More broadly, the phrase itself remains domain-dependent. In lattice QCD it now has a concrete and technically specific meaning tied to large-momentum hadrons. In atomistic, spline-theoretic, and nonlinear approximation settings, any broader use is interpretive rather than terminologically standard. A plausible long-term implication is that the phrase may continue to bifurcate: one branch centered on boosted-hadron source design and large-momentum QCD, and another referring more generally to lattice operators that encode the relevant affine, helicity, or geometric kinematics directly into the interpolation mechanism.