Inverse problem in the LaMET framework
Pith reviewed 2026-05-22 18:57 UTC · model grok-4.3
The pith
The inverse problem from incomplete lattice harmonics creates large uncertainty in LaMET parton distributions unless strong assumptions are added.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
Using non-perturbative lattice data, the authors show that the uncertainty introduced by this inverse problem in a realistic setup remains significant without very restrictive assumptions, and that the importance of the exact asymptotic behavior is minimal for values of x where the framework is currently applicable. The crux of the inverse problem lies in harmonics of the order of λ = z P_z ∼ 5–15, where the signal in the lattice data is often barely existent or unreliable in current studies.
What carries the argument
The Fourier transform required by LaMET to recover the x-dependence of parton distributions from matrix elements known only over a finite, noisy range of harmonics.
If this is right
- Uncertainties in parton distributions extracted via LaMET stay large with present-day lattice data quality.
- The precise asymptotic decay of missing harmonics has minimal impact on results at accessible x.
- More sophisticated techniques are needed to control the inverse problem in LaMET and in short-distance factorization.
- LaMET does not currently deliver a direct computation of the full x-dependence, contrary to a common misconception.
Where Pith is reading between the lines
- Future lattice work that strengthens the signal specifically in the λ ∼ 5–15 window could shrink the inverse-problem uncertainty without needing extreme assumptions.
- The same truncation issue is likely to appear in any factorization scheme that converts limited lattice matrix elements into momentum-space distributions.
- Cross-checking LaMET results against moment-based extractions on the same ensembles would help quantify how much of the spread comes from the inverse problem.
Load-bearing premise
The non-perturbative lattice data sets employed are representative of typical current calculations and that the signal in the harmonics λ = z P_z ∼ 5–15 is indeed barely existent or unreliable.
What would settle it
A new lattice calculation that delivers reliable data for harmonics well above λ = 15 and shows that enforcing or removing an exponential tail changes the extracted distributions by a large amount at currently accessible x would test the claim.
Figures
read the original abstract
One proposal to compute parton distributions from first principles is the large momentum effective theory (LaMET), which requires the Fourier transform of matrix elements computed non-perturbatively. Lattice quantum chromodynamics (QCD) provides calculations of these matrix elements over a finite range of Fourier harmonics that are often noisy or unreliable in the largest computed harmonics. It has been suggested that enforcing an exponential decay of the missing harmonics helps alleviate this issue. Using non-perturbative data, we show that the uncertainty introduced by this inverse problem in a realistic setup remains significant without very restrictive assumptions, and that the importance of the exact asymptotic behavior is minimal for values of $x$ where the framework is currently applicable. We show that the crux of the inverse problem lies in harmonics of the order of $\lambda=zP_z \sim 5-15$, where the signal in the lattice data is often barely existent in current studies, and the asymptotic behavior is not firmly established. We stress the need for more sophisticated techniques to account for this inverse problem, whether in the LaMET or related frameworks like the short-distance factorization. We also address a misconception that, with available lattice methods, the LaMET framework allows a "direct" computation of the $x$-dependence, whereas the alternative short-distance factorization only gives access to moments or fits of the $x$-dependence.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript examines the inverse problem in the LaMET framework for computing parton distribution functions from lattice QCD matrix elements. Using non-perturbative data, the authors argue that the uncertainty arising from the finite range of Fourier harmonics remains significant without very restrictive assumptions on the missing terms, that the precise asymptotic decay has minimal impact for the x values where LaMET is currently applicable, and that the crux lies in the harmonics with λ = z P_z ∼ 5–15 where the lattice signal is typically weak or unreliable. The paper also addresses a misconception that current LaMET methods permit a direct computation of the x-dependence, in contrast to moment-based or fit-based approaches in short-distance factorization.
Significance. If the findings hold, the work would be significant for the lattice QCD community by quantifying practical limitations of the inverse problem in realistic setups and underscoring the need for more sophisticated techniques. The reliance on actual non-perturbative data to illustrate the points adds practical relevance, though the strength of the conclusions depends on the representativeness of the data sets employed.
major comments (1)
- [Sections discussing non-perturbative data and harmonics λ ∼ 5–15] The central claim that the inverse-problem uncertainty remains significant without restrictive assumptions, and that this affects current LaMET applications in general, rests on the specific non-perturbative lattice data sets used to demonstrate weak signal at λ ∼ 5–15. The manuscript must explicitly identify the ensembles (including references to the original lattice publications), detail parameters such as source-sink separation, smearing, and fitting procedures, and provide evidence that the observed unreliability in this harmonic range is representative of typical contemporary calculations rather than an artifact of the chosen data. Without this justification, the broader conclusion for the field does not fully follow.
minor comments (1)
- [Abstract and introduction] In the abstract and introduction, provide a brief quantitative example or reference to a figure showing the x-range where the framework is currently applicable to make the statement about minimal impact of asymptotic behavior more concrete.
Simulated Author's Rebuttal
We thank the referee for the careful reading of our manuscript and for the constructive feedback on the presentation of our non-perturbative results. We address the single major comment below and will revise the manuscript to incorporate additional details as requested.
read point-by-point responses
-
Referee: [Sections discussing non-perturbative data and harmonics λ ∼ 5–15] The central claim that the inverse-problem uncertainty remains significant without restrictive assumptions, and that this affects current LaMET applications in general, rests on the specific non-perturbative lattice data sets used to demonstrate weak signal at λ ∼ 5–15. The manuscript must explicitly identify the ensembles (including references to the original lattice publications), detail parameters such as source-sink separation, smearing, and fitting procedures, and provide evidence that the observed unreliability in this harmonic range is representative of typical contemporary calculations rather than an artifact of the chosen data. Without this justification, the broader conclusion for the field does not fully follow.
Authors: We agree that explicit identification of the data sets strengthens the manuscript and helps readers evaluate the scope of our conclusions. In the revised version we will add a new subsection (or appendix) that fully specifies the lattice ensembles employed, including complete references to the original publications from which the matrix elements were taken. We will also document the source-sink separations, smearing parameters, and fitting procedures used to obtain the quoted matrix elements. To address representativeness, we will include a brief comparison with other recent LaMET calculations in the literature, noting that the rapid degradation of the signal-to-noise ratio for λ = z P_z ≳ 5–10 is a widely reported feature arising from the exponential fall-off of the correlators and the practical limits on achievable P_z. This discussion will clarify that the behavior we observe is characteristic of current-generation ensembles rather than an isolated artifact. revision: yes
Circularity Check
Analysis relies on external lattice data; no internal circularity in derivation
full rationale
The paper's central claim—that inverse-problem uncertainty remains significant without restrictive assumptions—is supported by analysis of existing non-perturbative lattice matrix elements rather than any quantity defined or fitted inside the paper itself. No self-definitional steps, fitted inputs renamed as predictions, or load-bearing self-citations appear in the provided abstract or described derivation chain. The work draws on external benchmarks (current lattice studies) whose signal quality at λ ∼ 5-15 is treated as representative input, not generated by the present analysis. This yields a minor score consistent with normal self-citation that is not load-bearing for the result.
Axiom & Free-Parameter Ledger
axioms (1)
- domain assumption Lattice QCD provides calculations of matrix elements over a finite range of Fourier harmonics that are often noisy or unreliable in the largest computed harmonics.
Lean theorems connected to this paper
-
IndisputableMonolith/Cost/FunctionalEquation.leanwashburn_uniqueness_aczel unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
Using non-perturbative data, we show that the uncertainty introduced by this inverse problem... remains significant without very restrictive assumptions
What do these tags mean?
- matches
- The paper's claim is directly supported by a theorem in the formal canon.
- supports
- The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
- extends
- The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
- uses
- The paper appears to rely on the theorem as machinery.
- contradicts
- The paper's claim conflicts with a theorem or certificate in the canon.
- unclear
- Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.
Forward citations
Cited by 3 Pith papers
-
Reconstructing the full kinematic dependence of GPDs from pseudo-distributions
Lattice QCD pseudo-distributions at m_π=358 MeV are inverted via multidimensional Gaussian process regression to reconstruct the full kinematic dependence of GPDs H^{u-d} and E^{u-d} while directly extracting double d...
-
Extracting Mellin moments of double parton distributions from lattice data
Skewness dependence in hadronic matrix elements must be accounted for when extracting Mellin moments of double parton distributions from lattice QCD, and model studies show its quantitative impact on current data.
-
Extracting Mellin moments of double parton distributions from lattice data
The skewness dependence of hadronic correlation functions affects Mellin moment extraction for double parton distributions from existing lattice data, as quantified through several models.
Reference graph
Works this paper leans on
-
[1]
AR acknowledges support by U.S. DOE Grant #DE-FG02-97ER41028. KO was supported in part by U.S. DOE Grant #DE-FG02-04ER41302 and would like to acknowledge the hospitality of the American Academy in Rome, where he spent part of his sabbatical. The re- search was conducted in part (DR) under the Laboratory- Directed Research and Development Program at Thomas...
work page 2020
-
[2]
Parton Physics on Euclidean Lattice
X. Ji, Parton Physics on a Euclidean Lattice, Phys. Rev. Lett. 110, 262002 (2013), arXiv:1305.1539 [hep-ph]
work page internal anchor Pith review Pith/arXiv arXiv 2013
-
[3]
Parton Physics from Large-Momentum Effective Field Theory
X. Ji, Parton Physics from Large-Momentum Effective Field Theory, Sci. China Phys. Mech. Astron. 57, 1407 (2014), arXiv:1404.6680 [hep-ph]
work page internal anchor Pith review Pith/arXiv arXiv 2014
-
[4]
Factorization Theorem Relating Euclidean and Light-Cone Parton Distributions
T. Izubuchi, X. Ji, L. Jin, I. W. Stewart, and Y. Zhao, Factorization Theorem Relating Euclidean and Light- Cone Parton Distributions, Phys. Rev. D 98, 056004 (2018), arXiv:1801.03917 [hep-ph]
work page internal anchor Pith review Pith/arXiv arXiv 2018
- [5]
- [6]
- [7]
- [8]
- [9]
-
[10]
On the computation of hadron-to-hadron transition matrix elements in lattice QCD
J. Bulava, M. Donnellan, and R. Sommer, On the com- putation of hadron-to-hadron transition matrix elements in lattice QCD, JHEP 01, 140, arXiv:1108.3774 [hep-lat]
work page internal anchor Pith review Pith/arXiv arXiv
-
[11]
G. S. Bali, B. Lang, B. U. Musch, and A. Sch¨ afer, Novel quark smearing for hadrons with high momenta in lattice QCD, Phys. Rev. D 93, 094515 (2016), arXiv:1602.05525 [hep-lat]
work page internal anchor Pith review Pith/arXiv arXiv 2016
-
[12]
B. U. Musch, P. Hagler, J. W. Negele, and A. Schafer, Exploring quark transverse momentum distributions with lattice QCD, Phys. Rev. D 83, 094507 (2011), arXiv:1011.1213 [hep-lat]
work page internal anchor Pith review Pith/arXiv arXiv 2011
-
[13]
Lattice QCD exploration of pseudo-PDFs
K. Orginos, A. Radyushkin, J. Karpie, and S. Zafeiropou- los, Lattice QCD exploration of parton pseudo- distribution functions, Phys. Rev. D 96, 094503 (2017), arXiv:1706.05373 [hep-ph]
work page internal anchor Pith review Pith/arXiv arXiv 2017
-
[14]
A complete non-perturbative renormalization prescription for quasi-PDFs
C. Alexandrou, K. Cichy, M. Constantinou, K. Had- jiyiannakou, K. Jansen, H. Panagopoulos, and F. Stef- fens, A complete non-perturbative renormalization pre- scription for quasi-PDFs, Nucl. Phys. B 923, 394 (2017), arXiv:1706.00265 [hep-lat]
work page internal anchor Pith review Pith/arXiv arXiv 2017
-
[15]
K. Zhang, Y.-Y. Li, Y.-K. Huo, A. Sch¨ afer, P. Sun, and Y.-B. Yang ( χQCD), RI/MOM renormalization of the parton quasidistribution functions in lattice regulariza- tion, Phys. Rev. D 104, 074501 (2021), arXiv:2012.05448 [hep-lat]. 9
- [16]
- [17]
-
[18]
R. Zhang, H.-W. Lin, and B. Yoon, Probing nucleon strange and charm distributions with lattice QCD, Phys. Rev. D 104, 094511 (2021), arXiv:2005.01124 [hep-lat]
-
[19]
C. Alexandrou, M. Constantinou, K. Hadjiyiannakou, K. Jansen, and F. Manigrasso, Flavor decomposition of the nucleon unpolarized, helicity, and transversity par- ton distribution functions from lattice QCD simulations, Phys. Rev. D 104, 054503 (2021), arXiv:2106.16065 [hep- lat]
- [20]
-
[21]
C. Egerer et al. (HadStruc), Transversity parton dis- tribution function of the nucleon using the pseudodis- tribution approach, Phys. Rev. D 105, 034507 (2022), arXiv:2111.01808 [hep-lat]
- [22]
- [23]
- [24]
- [25]
- [26]
-
[27]
Exclusive processes in position space and the pion distribution amplitude
V. Braun and D. M¨ uller, Exclusive processes in position space and the pion distribution amplitude, Eur. Phys. J. C55, 349 (2008), arXiv:0709.1348 [hep-ph]
work page internal anchor Pith review Pith/arXiv arXiv 2008
-
[28]
A. V. Radyushkin, Quasi-parton distribution func- tions, momentum distributions, and pseudo-parton dis- tribution functions, Phys. Rev. D 96, 034025 (2017), arXiv:1705.01488 [hep-ph]
work page internal anchor Pith review Pith/arXiv arXiv 2017
-
[29]
Exploring hadrons' partonic structure using ab initio lattice QCD calculations
Y.-Q. Ma and J.-W. Qiu, Exploring Partonic Structure of Hadrons Using ab initio Lattice QCD Calculations, Phys. Rev. Lett. 120, 022003 (2018), arXiv:1709.03018 [hep-ph]
work page internal anchor Pith review Pith/arXiv arXiv 2018
-
[30]
J. Karpie, K. Orginos, A. Rothkopf, and S. Zafeiropoulos, Reconstructing parton distribution functions from Ioffe time data: from Bayesian methods to Neural Networks, JHEP 04, 057, arXiv:1901.05408 [hep-lat]
work page internal anchor Pith review Pith/arXiv arXiv 1901
-
[31]
H. Dutrieux, J. Karpie, K. Orginos, and S. Zafeiropou- los, Simple nonparametric reconstruction of parton dis- tributions from limited Fourier information, Phys. Rev. D 111, 034515 (2025), arXiv:2412.05227 [hep-lat]
-
[32]
S. Bhattacharya, K. Cichy, M. Constantinou, J. Dodson, A. Metz, A. Scapellato, and F. Steffens, Chiral-even axial twist-3 GPDs of the proton from lattice QCD, Phys. Rev. D 108, 054501 (2023), arXiv:2306.05533 [hep-lat]
-
[33]
Ji, Large-Momentum Effective Theory vs
X. Ji, Large-Momentum Effective Theory vs. Short- Distance Operator Expansion: Contrast and Comple- mentarity, (2022), arXiv:2209.09332 [hep-lat]
- [34]
- [35]
-
[36]
E. Rietsch, The maximum entropy approach to inverse problems - spectral analysis of short data records and density structure of the Earth., Journal of Geophysics 42, 489 (1976)
work page 1976
-
[37]
A novel Bayesian approach to spectral function reconstruction
Y. Burnier and A. Rothkopf, Bayesian Approach to Spec- tral Function Reconstruction for Euclidean Quantum Field Theories, Phys. Rev. Lett. 111, 182003 (2013), arXiv:1307.6106 [hep-lat]
work page internal anchor Pith review Pith/arXiv arXiv 2013
- [38]
-
[39]
G. Backus and F. Gilbert, The resolving power of gross earth data, Geophysical Journal International 16, 169 (1968), https://academic.oup.com/gji/article- pdf/16/2/169/5891044/16-2-169.pdf
work page 1968
-
[40]
L. Del Debbio, A. Lupo, M. Panero, and N. Tantalo, Bayesian solution to the inverse problem and its rela- tion to Backus–Gilbert methods, Eur. Phys. J. C 85, 185 (2025), arXiv:2409.04413 [hep-lat]
-
[41]
Lattice calculation of hadronic tensor of the nucleon
J. Liang, K.-F. Liu, and Y.-B. Yang, Lattice calculation of hadronic tensor of the nucleon, EPJ Web Conf. 175, 14014 (2018), arXiv:1710.11145 [hep-lat]
work page internal anchor Pith review Pith/arXiv arXiv 2018
- [42]
-
[43]
C. Alexandrou, K. Cichy, M. Constantinou, K. Had- jiyiannakou, K. Jansen, A. Scapellato, and F. Steffens, Transversity GPDs of the proton from lattice QCD, Phys. Rev. D 105, 034501 (2022), arXiv:2108.10789 [hep-lat]
- [44]
-
[45]
S. Bhattacharya et al., Generalized parton distributions from lattice QCD with asymmetric momentum trans- fer: Axial-vector case, Phys. Rev. D 109, 034508 (2024), arXiv:2310.13114 [hep-lat]
-
[46]
On the extraction of spectral densities from lattice correlators
M. Hansen, A. Lupo, and N. Tantalo, Extraction of spec- tral densities from lattice correlators, Phys. Rev. D 99, 094508 (2019), arXiv:1903.06476 [hep-lat]. 10
work page internal anchor Pith review Pith/arXiv arXiv 2019
- [47]
-
[48]
Since the value at λ = 0 is exactly 1 by construction of the dataset, the standard χ2 is infinite without regu- larization. Therefore, we modify the data covariance to attribute an uncorrelated standard deviation of 10 −3 to that point. In contrast, applying a simple SVD cut would effectively exclude the λ = 0 point from contributing to the χ2
-
[49]
J. M. Karpie, On the calculation of parton distributions from Lattice QCD, Ph.D. thesis, William-Mary Coll. (2019)
work page 2019
-
[50]
A. Candido, L. Del Debbio, T. Giani, and G. Petrillo, Bayesian inference with Gaussian processes for the de- termination of parton distribution functions, Eur. Phys. J. C 84, 716 (2024), arXiv:2404.07573 [hep-ph]
-
[51]
We only consider positive values of x since we will only be concerned by the real part of the matrix elements in this numerical exploration, which corresponds to a definite parity in x
-
[52]
For the logarithmical RBF, we use σ = 1, l = 0.7 and a prior mean f(x) = x−0.15
For the RBF kernel, we use σ = 10, l = 0.4 and a prior mean of f(x) = 10. For the logarithmical RBF, we use σ = 1, l = 0.7 and a prior mean f(x) = x−0.15
-
[53]
For the RBF kernel in x-space, we use σ = 50, l = 0 .4 and prior mean of f(x) = 3 while the kernel in λ-space is Σ = 2 and L = 5. For the logarithmical RBF kernel in x-space, we use σ = 4, l = 0 .7 and prior mean of f(x) = x−0.15 while the kernel in λ-space is Σ = 0 .5 and L = 5. The values are the same for both the exponential and Gaussian reconstructions
-
[54]
Ji, Euclidean effective theory for partons in the spirit of Steven Weinberg, Nucl
X. Ji, Euclidean effective theory for partons in the spirit of Steven Weinberg, Nucl. Phys. B 1007, 116670 (2024), arXiv:2408.03378 [hep-ph]
- [55]
-
[56]
Solving Lattice QCD systems of equations using mixed precision solvers on GPUs
M. Clark, R. Babich, K. Barros, R. Brower, and C. Rebbi, Solving Lattice QCD systems of equations using mixed precision solvers on GPUs, Comput. Phys. Commun. 181, 1517 (2010), arXiv:0911.3191 [hep-lat]
work page internal anchor Pith review Pith/arXiv arXiv 2010
-
[57]
Parallelizing the QUDA Library for Multi-GPU Calculations in Lattice Quantum Chromodynamics
R. Babich, M. A. Clark, and B. Joo, Paral- lelizing the QUDA Library for Multi-GPU Cal- culations in Lattice Quantum Chromodynamics, in SC 10 (Supercomputing 2010) (2010) arXiv:1011.0024 [hep-lat]
work page internal anchor Pith review Pith/arXiv arXiv 2010
-
[58]
B. Jo´ o, D. D. Kalamkar, T. Kurth, K. Vaidyanathan, and A. Walden, Optimizing Wilson-Dirac Op- erator and Linear Solvers for Intel ® KNL, in High Performance Computing: ISC High Performance 2016 International Workshops, ExaComm, E-MuCoCoS, HPC-IODC, IXPUG, IWOPH, P 3MA, VHPC, WOPSSS, Frankfurt, Germany, June 19–23, 2016, Revised Selected Papers, edited b...
work page 2016
-
[59]
J. Chen, R. G. Edwards, and W. Mao, Graph Contrac- tions for Calculating Correlation Functions in Lattice QCD, in Platform for Advanced Scientific Computing (2023)
work page 2023
-
[60]
D. Stanzione, J. West, R. T. Evans, T. Minyard, O. Ghat- tas, and D. K. Panda, Frontera: The evolution of lead- ership computing at the national science foundation, in Practice and Experience in Advanced Research Computing, PEARC ’20 (Association for Computing Machinery, New York, NY, USA, 2020) p. 106–111
work page 2020
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.