pith. sign in

arxiv: 2504.17706 · v2 · submitted 2025-04-24 · ✦ hep-lat · hep-ph

Inverse problem in the LaMET framework

Pith reviewed 2026-05-22 18:57 UTC · model grok-4.3

classification ✦ hep-lat hep-ph
keywords LaMETparton distribution functionslattice QCDinverse problemFourier transformhadron structure
0
0 comments X

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.

The paper studies the inverse problem that arises when LaMET is used to extract parton distributions from lattice QCD. Lattice calculations supply matrix elements only up to moderate Fourier harmonics, beyond which the data are noisy or absent. By examining actual non-perturbative results, the authors show that this truncation produces sizable uncertainties in the reconstructed distributions unless very restrictive assumptions are imposed on the missing terms. They further find that the precise form of the asymptotic tail has little practical effect within the x range where LaMET is currently usable. The work calls for more advanced methods to handle the inverse problem in both LaMET and related short-distance factorization approaches.

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

These are editorial extensions of the paper, not claims the author makes directly.

  • 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

Figures reproduced from arXiv: 2504.17706 by Anatoly Radyushkin, Christopher J. Monahan, David Richards, Herv\'e Dutrieux, Joe Karpie, Kostas Orginos, Savvas Zafeiropoulos.

Figure 1
Figure 1. Figure 1: FIG. 1. (top) Various BG reconstructions of the same Fourier [PITH_FULL_IMAGE:figures/full_fig_p004_1.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4. (top) Reconstructions / extrapolations of the same [PITH_FULL_IMAGE:figures/full_fig_p006_4.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3. (top) Reconstructions / extrapolations of the same [PITH_FULL_IMAGE:figures/full_fig_p006_3.png] view at source ↗
Figure 5
Figure 5. Figure 5: FIG. 5. Various [PITH_FULL_IMAGE:figures/full_fig_p007_5.png] view at source ↗
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.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

1 major / 1 minor

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)
  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)
  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

1 responses · 0 unresolved

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
  1. 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

0 steps flagged

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

0 free parameters · 1 axioms · 0 invented entities

The central claim rests on standard lattice QCD domain assumptions about finite-range matrix elements and noisy high harmonics; no new free parameters, invented entities, or ad-hoc axioms are introduced in the abstract.

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.
    Invoked in the opening description of the LaMET setup and the inverse problem.

pith-pipeline@v0.9.0 · 5792 in / 1264 out tokens · 41502 ms · 2026-05-22T18:57:07.446485+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

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

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score.

  1. Reconstructing the full kinematic dependence of GPDs from pseudo-distributions

    hep-lat 2026-04 unverdicted novelty 8.0

    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...

  2. Extracting Mellin moments of double parton distributions from lattice data

    hep-lat 2025-11 unverdicted novelty 5.0

    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.

  3. Extracting Mellin moments of double parton distributions from lattice data

    hep-lat 2025-11 unverdicted novelty 5.0

    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

60 extracted references · 60 canonical work pages · cited by 2 Pith papers · 17 internal anchors

  1. [1]

    DOE Grant #DE-FG02-97ER41028

    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...

  2. [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]

  3. [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]

  4. [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]

  5. [5]

    L.-B. Chen, W. Wang, and R. Zhu, Next-to-Next-to- Leading Order Calculation of Quasiparton Distribu- tion Functions, Phys. Rev. Lett. 126, 072002 (2021), arXiv:2006.14825 [hep-ph]

  6. [6]

    Li, Y.-Q

    Z.-Y. Li, Y.-Q. Ma, and J.-W. Qiu, Extraction of Next- to-Next-to-Leading-Order Parton Distribution Functions from Lattice QCD Calculations, Phys. Rev. Lett. 126, 072001 (2021), arXiv:2006.12370 [hep-ph]

  7. [7]

    V. M. Braun, A. Vladimirov, and J.-H. Zhang, Power corrections and renormalons in parton quasidistributions, Phys. Rev. D 99, 014013 (2019), arXiv:1810.00048 [hep- ph]

  8. [8]

    Zhang, J

    R. Zhang, J. Holligan, X. Ji, and Y. Su, Leading power accuracy in lattice calculations of parton distributions, Phys. Lett. B 844, 138081 (2023), arXiv:2305.05212 [hep- lat]

  9. [9]

    Zhang, A

    R. Zhang, A. V. Grebe, D. C. Hackett, M. L. Wagman, and Y. Zhao, Kinematically-enhanced interpolating op- erators for boosted hadrons, (2025), arXiv:2501.00729 [hep-lat]

  10. [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]

  11. [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]

  12. [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]

  13. [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]

  14. [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]

  15. [15]

    Zhang, Y.-Y

    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. [16]

    X. Ji, Y. Liu, A. Sch¨ afer, W. Wang, Y.-B. Yang, J.-H. Zhang, and Y. Zhao, A Hybrid Renormalization Scheme for Quasi Light-Front Correlations in Large-Momentum Effective Theory, Nucl. Phys. B 964, 115311 (2021), arXiv:2008.03886 [hep-ph]

  17. [17]

    Huo et al

    Y.-K. Huo et al. (Lattice Parton (LPC)), Self- renormalization of quasi-light-front correlators on the lat- tice, Nucl. Phys. B 969, 115443 (2021), arXiv:2103.02965 [hep-lat]

  18. [18]

    Zhang, H.-W

    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. [19]

    Alexandrou, M

    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. [20]

    X. Gao, A. D. Hanlon, S. Mukherjee, P. Petreczky, P. Scior, S. Syritsyn, and Y. Zhao, Lattice QCD Determi- nation of the Bjorken-x Dependence of Parton Distribu- tion Functions at Next-to-Next-to-Leading Order, Phys. Rev. Lett. 128, 142003 (2022), arXiv:2112.02208 [hep- lat]

  21. [21]

    Egerer et al

    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. [22]

    Z. Fan, W. Good, and H.-W. Lin, Gluon parton distribu- tion of the nucleon from (2+1+1)-flavor lattice QCD in the physical-continuum limit, Phys. Rev. D 108, 014508 (2023), arXiv:2210.09985 [hep-lat]

  23. [23]

    X. Gao, A. D. Hanlon, S. Mukherjee, P. Petreczky, Q. Shi, S. Syritsyn, and Y. Zhao, Transversity PDFs of the proton from lattice QCD with physical quark masses, Phys. Rev. D 109, 054506 (2024), arXiv:2310.19047 [hep- lat]

  24. [24]

    Delmar, C

    J. Delmar, C. Alexandrou, K. Cichy, M. Constantinou, and K. Hadjiyiannakou, Gluon PDF of the proton using twisted mass fermions, Phys. Rev. D 108, 094515 (2023), arXiv:2310.01389 [hep-lat]

  25. [25]

    H.-T. Ding, X. Gao, S. Mukherjee, P. Petreczky, Q. Shi, S. Syritsyn, and Y. Zhao, Three-dimensional imaging of pion using lattice QCD: generalized parton distributions, JHEP 02, 056, arXiv:2407.03516 [hep-lat]

  26. [26]

    W. Good, K. Hasan, and H.-W. Lin, Toward the first gluon parton distribution from the LaMET, J. Phys. G 52, 035105 (2025), arXiv:2409.02750 [hep-lat]

  27. [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]

  28. [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]

  29. [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]

  30. [30]

    Reconstructing parton distribution functions from Ioffe time data: from Bayesian methods to Neural Networks

    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]

  31. [31]

    Dutrieux, J

    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. [32]

    Bhattacharya, K

    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. [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. [34]

    Egerer, R

    C. Egerer, R. G. Edwards, C. Kallidonis, K. Orginos, A. V. Radyushkin, D. G. Richards, E. Romero, and S. Zafeiropoulos (HadStruc), Towards high-precision par- ton distributions from lattice QCD via distillation, JHEP 11, 148, arXiv:2107.05199 [hep-lat]

  35. [35]

    V. M. Braun, M. Koller, and J. Schoenleber, Renormalons and power corrections in pseudo- and quasi-GPDs, Phys. Rev. D 109, 074510 (2024), arXiv:2401.08012 [hep-ph]

  36. [36]

    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)

    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)

  37. [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]

  38. [38]

    Liang, T

    J. Liang, T. Draper, K.-F. Liu, A. Rothkopf, and Y.- B. Yang (XQCD), Towards the nucleon hadronic ten- sor from lattice QCD, Phys. Rev. D 101, 114503 (2020), arXiv:1906.05312 [hep-ph]

  39. [39]

    Backus and F

    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

  40. [40]

    Del Debbio, A

    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. [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]

  42. [42]

    M. Bhat, K. Cichy, M. Constantinou, and A. Scapellato, Flavor nonsinglet parton distribution functions from lat- tice QCD at physical quark masses via the pseudodis- tribution approach, Phys. Rev. D 103, 034510 (2021), arXiv:2005.02102 [hep-lat]

  43. [43]

    Alexandrou, K

    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. [44]

    M. Bhat, W. Chomicki, K. Cichy, M. Constantinou, J. R. Green, and A. Scapellato, Continuum limit of parton distribution functions from the pseudodistribution ap- proach on the lattice, Phys. Rev. D 106, 054504 (2022), arXiv:2205.07585 [hep-lat]

  45. [45]

    Bhattacharya et al., Generalized parton distributions from lattice QCD with asymmetric momentum trans- fer: Axial-vector case, Phys

    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. [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

  47. [47]

    X. Gao, A. D. Hanlon, J. Holligan, N. Karthik, S. Mukherjee, P. Petreczky, S. Syritsyn, and Y. Zhao, Un- polarized proton PDF at NNLO from lattice QCD with physical quark masses, Phys. Rev. D 107, 074509 (2023), arXiv:2212.12569 [hep-lat]

  48. [48]

    Therefore, we modify the data covariance to attribute an uncorrelated standard deviation of 10 −3 to that point

    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. [49]

    J. M. Karpie, On the calculation of parton distributions from Lattice QCD, Ph.D. thesis, William-Mary Coll. (2019)

  50. [50]

    Candido, L

    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. [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. [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. [53]

    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

    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. [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. [55]

    R. G. Edwards and B. Joo (SciDAC, LHPC, UKQCD), The Chroma software system for lattice QCD, Nucl. Phys. B Proc. Suppl. 140, 832 (2005), arXiv:hep- lat/0409003

  56. [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]

  57. [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]

  58. [58]

    Jo´ o, D

    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...

  59. [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)

  60. [60]

    Stanzione, J

    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