REVIEW 4 major objections 4 minor 1 cited by
Unpolarized gluon PDF of the nucleon from lattice QCD at physical point in the continuum limit
T0 review · 4 major / 4 minor · reviewed 2026-08-04 · deepseek-v4-flash
Pith's one-line read Lattice QCD, extrapolated to the continuum and infinite-momentum limits, reproduces the nucleon's unpolarized gluon PDF within global-fit uncertainty.
desk verdict First LaMET gluon-PDF extraction with three lattice spacings plus a momentum extrapolation, but the abstract overclaims (five spacings, physical pion mass, 3 GeV are not in the body), and the final result is 'consistent with zero,' so agreement with global fits is weak. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The argument rests on a specific chain of techniques. First, a gluonic operator O(z) that is multiplicatively renormalizable, computed with distillation to improve the signal. Second, hybrid renormalization with a self-renormalization factor Z_R that removes linear and logarithmic ultraviolet divergences. Third, the large-λ extrapolation h_R(λ)=l1 λ^{-a1} e^{-λ/λ0}, used to extend the renormalized matrix elements from the measured range out to all distances before the Fourier transform. Fourth, the NLO hybrid-scheme matching kernel that converts the quasi-PDF to the light-cone PDF. Fifth, the joint continuum/infinite-momentum extrapolation ansatz of Eq. (8), whose intercept xg0(x) is the fin
What would settle it
Compute the renormalized matrix elements at higher nucleon momentum (P_z ≳ 3 GeV) on a finer lattice and check whether the large-λ data deviate significantly from the fitted power-times-exponential tail used here. Alternatively, independently compute the gluon momentum fraction <x>_g and test whether ∫ dx xg0(x) from the extrapolated curve is consistent with it within the combined uncertainty; a violation would falsify the extrapolated PDF.
Extended reading notes
Core claim
On its own terms, the paper establishes that the nucleon's unpolarized gluon PDF can be obtained from lattice QCD through the LaMET procedure. The renormalized quasi-PDF matrix elements are Fourier transformed to momentum space using a model-based tail extension, matched at one loop to the light-cone PDF, and finally extrapolated via xg(x,P_z,a) = xg0(x) + a^2 f(x) + a^2 P_z^2 h(x) + d(x)/P_z^2 to the continuum and infinite-momentum limits. The resulting xg0(x)/<x> is consistent with the current global fits within errors, with the large-x region close to zero (and slightly negative central values), indicating suppression of the gluon at high momentum fraction. The paper explicitly states tha
Load-bearing premise
The load-bearing premise is that the true large-distance (large λ = zP_z) behavior of the renormalized gluon quasi-PDF has the functional form l1 λ^{-a1} e^{-λ/λ0}, because the Fourier transform integrates this model over all λ while the lattice data only constrain it out to λ ≈ 10–14.
Editorial extensions
If this is right
- If the central claim holds, lattice QCD can supply direct, model-independent input for the gluon PDF in regions where global fits are weakly constrained, especially at large x.
- The consistency with global fits suggests that the LaMET machinery—hybrid renormalization, NLO matching, and extrapolation—is under control for gluonic operators, opening the way to other gluonic observables such as the gluon momentum fraction.
- The observed large-x suppression of the gluon distribution is a quantitative confirmation of a known qualitative feature of the nucleon's gluonic structure.
- The uncertainty budget identifies the combined continuum/infinite-momentum extrapolation and the tail model as the dominant systematics, indicating where future calculations (finer lattices, higher momenta, more spacings) will pay off.
Reading between the lines
- Editorial inference: The paper's title says 'physical point,' but the body uses a pion mass of roughly 300 MeV, not the physical 136 MeV; a chiral extrapolation over pion mass is not performed here.
- Editorial inference: The large-λ tail form is only constrained out to λ ≈ 10–14; the systematic error varies the fitting range but not the functional form, so a different tail shape would propagate directly into the small- and mid-x PDF.
- Editorial inference: The slightly negative central values in the extrapolated band are unphysical; a stronger test would check the momentum sum rule against an independent lattice determination of the gluon momentum fraction.
- Editorial inference: The continuum and infinite-momentum extrapolation uses three lattice spacings at one pion mass; adding finer spacings and higher momenta would test the stability of the a^2 and d(x)/P_z^2 terms in Eq. (8).
Signed reviews
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper reports a LaMET lattice QCD calculation of the unpolarized gluon PDF of the nucleon. Bare gluonic matrix elements are computed on three CLQCD ensembles (a = 0.105, 0.0897, 0.0775 fm; m_pi ~ 287-300 MeV) with distillation, renormalized in a hybrid scheme with parameters fitted to zero-momentum data, extended to large Ioffe time using a power-times-exponential tail, matched at NLO to the light-cone PDF, and then extrapolated to the continuum and infinite momentum via Eq. (8). The resulting xg(x)/<x> is compared to CT18, NNPDF, and JAM24 and reported as consistent within errors; the central value is consistent with zero at most x. The paper claims to be the most complete lattice gluon PDF determination to date.
Significance. If the result is correct, it would be one of the first continuum-extrapolated determinations of the nucleon gluon PDF from LaMET and would demonstrate that the large-x gluon is suppressed as in global fits. The work has clear strengths: multiple lattice spacings, high statistics, distillation, hybrid renormalization, NLO matching, and a breakdown of statistical and systematic uncertainties. The comparison with three global fits is a useful sanity check. However, the final uncertainty is very large, and the central value being consistent with zero means the agreement with CT18/NNPDF/JAM24 is a weak compatibility statement rather than a precision test. The main technical risk is the model dependence of the large-distance tail and the underconstrained continuum/infinite-momentum extrapolation; these need to be addressed before the central claim can be accepted.
major comments (4)
- [Light-cone PDF, Eq. (7), Eq. (14), Table IV] The Fourier transform in Eq. (14) requires h_R(z,P_z) at all lambda, but the lattice data end at lambda ~ 10-14, and for F32P30 only P_z <= 1.5 GeV is used. The unmeasured tail is modeled by Eq. (7), h_R = l1 lambda^{-a1} e^{-lambda/lambda0}, fitted per ensemble and momentum. The systematic uncertainty is estimated only by moving the fit range (Table IV), never by changing the functional form. Since a1 is only 'associated with' the endpoint power law and not derived from it, Eq. (7) is an ansatz; if the true tail decays differently, the integral in Eq. (14) shifts xg(x) across the whole x range, including the comparison region 0.2<x<0.8. This is a load-bearing model assumption and needs a derivation, alternative-form scan, or an insensitivity demonstration.
- [Eq. (8), Table I, Fig. 3] Equation (8) contains four x-dependent functions (xg0, f, h, d) for the continuum and infinite-momentum extrapolation. Per x, the available data are at most three momenta for C24P29/E32P29 and two for F32P30, so the fit has very few degrees of freedom; the separation between a^2 f(x), a^2 P_z^2 h(x), and d(x)/P_z^2 is weakly determined. No chi^2 or stability test for Eq. (8) is reported (e.g., dropping the a^2 P_z^2 term, or fitting only the two largest momenta). The large uncertainty of the gray band in Figs. 3-4 is a direct consequence. To support the central claim, the extrapolation fit quality and its robustness to ansatz variations must be presented.
- [Title, abstract, Table I] The title and abstract state the calculation is at the physical point with five lattice spacings, pion masses 136-317 MeV, and momenta up to 3 GeV. The body and Table I report three lattice spacings, m_pi ~ 287-300 MeV, and P_z up to 1.97 GeV. This is not a minor wording issue: the claimed parameter coverage does not match the actual data, and no chiral extrapolation is performed. The authors must correct the abstract/title to match the body, or provide the missing ensembles and analysis.
- [Results, after Eq. (8)] The ensembles have m_pi = 287-300 MeV, yet the text states 'we do not include pion mass dependence in this work, as its effect is expected to be weak based on previous studies [31].' The available data do not test this assumption because the pion-mass spread is only ~13 MeV. For a paper titled 'at physical point,' this is a central limitation; either a dedicated pion-mass study or a clear caveat in the abstract and conclusions is required.
minor comments (4)
- [Throughout] Several typos and grammatical errors should be corrected: 'nulceon' in the abstract, 'in in the supplemental material' near Eq. (8), 'resloved' in the supplementary, and 'anaylsis' in Ref. [46].
- [Fig. 12 caption] The caption says 'the width of the (nonoverlapping) colored bands denotes the size of each uncertainty,' but the bands in the figure appear to overlap. Please clarify how the widths are to be read.
- [Fig. 4] The gray bands marking x<0.2 and x>0.8 are introduced without justification. A brief explanation of why LaMET is unreliable in these regions would help the reader.
- [Conclusion] The claimed agreement with CT18/NNPDF/JAM24 is based on visual overlap. Given the very large final uncertainty, a quantitative measure (e.g., chi^2 per data point or a confidence level) would be more informative.
Circularity Check
No significant circularity: the CT18/NNPDF/JAM24 comparisons are external anchors, and the renormalization/tail fits are to lattice data rather than to the target PDF.
full rationale
The claimed prediction is the continuum/infinite-momentum xg0(x) from Eq. (8), obtained by Fourier transforming measured renormalized matrix elements (Eq. 14), NLO matching (Eq. 11), and a joint extrapolation. The global-fit PDFs (CT18, NNPDF, JAM24) enter only in the final comparison (Fig. 4); nothing is fitted to them, so the agreement is not forced. The self-renormalization parameters (k, Lambda_QCD, m0, d, g(z)) are fitted to zero-momentum bare matrix elements via Eq. (6) and then used in Eq. (3) to define renormalized long-distance matrix elements for all momenta; this is a standard renormalization condition, and the target xg0 is obtained from independently measured nonzero-momentum matrix elements, not from the fit function. The large-lambda tail (Eq. 7) is an explicit extrapolation ansatz fitted to the measured h_R in a restricted range (Table IV) and varied for systematics; while it is a genuine model-form risk (the functional form is never varied, and the Fourier integral extends beyond the measured lambda), the result is not equivalent to the tail model by construction - it is a data-driven extrapolation with a stated assumption. The citations to Refs. [13,40] for the hybrid scheme and tail form are published methodology with stated assumptions, not a self-citation uniqueness theorem. No step exhibits a parameter fitted to the target quantity and then renamed as a prediction.
Assumptions & free parameters
free parameters (8)
- k (linear divergence coefficient in Z_R) =
1.92(1.18)
- Lambda_QCD in Z_R =
0.59(0.15) GeV
- m0 (finite-mass renormalization ambiguity) =
2.78(1.51) GeV
- d (subleading-log coefficient) =
-0.079(0.403)
- g(z) (long-distance residual function) =
free function, not tabulated
- l1, a1, lambda0 (lambda-tail extrapolation) =
per (ensemble, Pz), not tabulated
- f(x), h(x), d(x) in Eq. (8) =
per-x functions, not tabulated
- Boundaries z0, z1, zs, z_max =
0.15, 0.3, 0.3, 1.0 fm
assumptions (7)
- domain assumption LaMET factorization (Eq. 11): quasi-PDF differs from the light-cone PDF by the NLO kernel plus power corrections Lambda_QCD^2/(xPz)^2 and Lambda_QCD^2/((1-x)Pz)^2.
- domain assumption NLO hybrid matching kernel (Eqs. 12-13) is adequate without RG, threshold, or renormalon resummation at mu = 2 GeV.
- domain assumption Hybrid renormalization: Z_R of Eq. (4) removes all UV and linear divergences, with the residual renormalization ambiguity absorbed by m0 z and g(z).
- ad hoc to paper Large-lambda tail has the form l1 lambda^{-a1} exp(-lambda/lambda0) (Eq. 7).
- ad hoc to paper Pion-mass dependence is negligible (m_pi = 287-300 MeV vs the physical point).
- domain assumption Discretization and higher-twist terms have the form a^2 f(x) + a^2 Pz^2 h(x) + d(x)/Pz^2 (Eq. 8).
- domain assumption The operator combination Eq. (1) is multiplicatively renormalizable and free of quark-operator mixing under this smearing.
Cite this review
Pith. "Pith review of Unpolarized gluon PDF of the nucleon from lattice QCD at physical point in the continuum limit." pith.science (2026). https://pith.science/paper/RCUSQAJQ
@misc{pith2026251026425,
author = {Pith},
title = {Pith review of: Unpolarized gluon PDF of the nucleon from lattice QCD at physical point in the continuum limit},
year = {2026},
howpublished = {\url{https://pith.science/paper/RCUSQAJQ}},
note = {Machine review of arXiv:2510.26425}
}
read the original abstract
We report a state-of-the-art lattice QCD calculation of the nucleon unpolarized gluon parton distribution function employing large-momentum effective theory. The calculation is carried out on the 2+1 flavor CLQCD ensembles with five lattice spacings a={0.105,0.0897,0.0775, 0.0688, 0.0519} fm and various pion masses ranging from 136 MeV to 317 MeV, covering nulceon momenta up to 3 GeV. Distillation technique is applied to improve the signal of two-point correlators. We then apply the state-of-the-art hybrid renormalization and one-loop perturbative matching, and extrapolate the result to the continuum limit, infinite momentum limit and physical pion mass.
Figures
Figures from the paper (8 more)
Forward citations
Cited by 1 Pith paper
-
Parton distribution functions from lattice QCD
A review of lattice-QCD approaches to parton distribution functions concludes that the field is moving from feasibility studies to quantitatively controlled calculations.
Reference graph
Works this paper leans on
-
[31]
Z. Fan, W. Good, and H.-W. Lin, Phys. Rev. D108, 014508 (2023), 2210.09985
arXiv 2023
- [1]
- [2]
-
[3]
R. D. Ball et al. (NNPDF), Eur. Phys. J. C77, 663 (2017), 1706.00428
arXiv 2017
-
[4]
R. D. Ball et al. (NNPDF), Eur. Phys. J. C82, 428 (2022), 2109.02653
arXiv 2022
-
[5]
A. Accardi, L. T. Brady, W. Melnitchouk, J. F. Owens, and N. Sato, Phys. Rev. D93, 114017 (2016), 1602.03154
arXiv 2016
-
[6]
S. Dulat, T.-J. Hou, J. Gao, M. Guzzi, J. Huston, P. Nadolsky, J. Pumplin, C. Schmidt, D. Stump, and C. P. Yuan, Phys. Rev. D93, 033006 (2016), 1506.07443
arXiv 2016
- [7]
Show all 57 references
-
[8]
X. Ji, Sci. China Phys. Mech. Astron.57, 1407 (2014), 1404.6680
2014 arXiv
- [9]
-
[10]
Lin, J.-W
H.-W. Lin, J.-W. Chen, X. Ji, L. Jin, R. Li, Y.-S. Liu, Y.-B. Yang, J.-H. Zhang, and Y. Zhao, Phys. Rev. Lett. 121, 242003 (2018), 1807.07431
2018 arXiv
- [11]
-
[12]
Zhang et al
Q.-A. Zhang et al. (Lattice Parton), Phys. Rev. Lett. 125, 192001 (2020), 2005.14572
2020 arXiv
-
[13]
Huo et al
Y.-K. Huo et al. (Lattice Parton (LPC)), Nucl. Phys. B 969, 115443 (2021), 2103.02965
2021 arXiv
-
[14]
Hua, M.-H
J. Hua, M.-H. Chu, P. Sun, W. Wang, J. Xu, Y.-B. Yang, J.-H. Zhang, and Q.-A. Zhang (Lattice Parton), Phys. Rev. Lett.127, 062002 (2021), 2011.09788
2021 arXiv
-
[15]
Chu et al
M.-H. Chu et al. (Lattice Parton (LPC)), Phys. Rev. D 106, 034509 (2022), 2204.00200
2022 arXiv
-
[16]
Zhang, X
K. Zhang, X. Ji, Y.-B. Yang, F. Yao, and J.-H. Zhang (Lattice Parton (LPC)), Phys. Rev. Lett.129, 082002 (2022), 2205.13402
2022 arXiv
- [17]
- [18]
- [19]
-
[20]
He, M.-H
J.-C. He, M.-H. Chu, J. Hua, X. Ji, A. Sch¨ afer, Y. Su, W. Wang, Y.-B. Yang, J.-H. Zhang, and Q.-A. Zhang (Lattice Parton Collaboration (LPC)), Phys. Rev. D 109, 114513 (2024), 2211.02340
2024 arXiv
- [21]
-
[22]
Zhang, Y.-K
K. Zhang, Y.-K. Huo, X. Ji, A. Schaefer, C.-J. Shi, P. Sun, W. Wang, Y.-B. Yang, and J.-H. Zhang (Lattice Parton), Phys. Rev. D110, 074505 (2024), 2405.14097
2024 arXiv
-
[23]
C. Chen, Y. Geng, L. Liu, P. Sun, Y.-B. Yang, F. Yao, J.-H. Zhang, and K. Zhang (CLQCD, Lattice Parton), Phys. Rev. D111, 074506 (2025), 2408.12819
2025 arXiv
- [24]
- [25]
-
[26]
Zhang, J.-C
Q.-A. Zhang, J.-C. He, M.-H. Chu, J. Hua, X. Ji, A. Sch¨ afer, Y. Su, W. Wang, Y.-B. Yang, and J.-H. Zhang, PoSQNP2024, 052 (2025)
2025
-
[27]
W. Good, K. Hasan, and H.-W. Lin, J. Phys. G52, 035105 (2025), 2409.02750
2025 arXiv
-
[28]
Fan, Y.-B
Z.-Y. Fan, Y.-B. Yang, A. Anthony, H.-W. Lin, and K.-F. Liu, Phys. Rev. Lett.121, 242001 (2018), 1808.02077
2018 arXiv
-
[29]
Z. Fan, R. Zhang, and H.-W. Lin, Int. J. Mod. Phys. A 36, 2150080 (2021), 2007.16113
2021 arXiv
- [30]
-
[32]
W. Good, Z. Fan, and H.-W. Lin, in30th International Workshop on Deep-Inelastic Scattering and Related Sub- jects(2023), 2307.06916
2023 arXiv
-
[33]
Delmar, C
J. Delmar, C. Alexandrou, K. Cichy, M. Constanti- nou, and K. Hadjiyiannakou, Phys. Rev. D108, 094515 (2023), 2310.01389
2023 arXiv
- [34]
-
[35]
W. Good, K. Hasan, A. Chevis, and H.-W. Lin, Phys. Rev. D109, 114509 (2024), 2310.12034
2024 arXiv
-
[36]
Salas-Chavira, Z
A. Salas-Chavira, Z. Fan, and H.-W. Lin, Phys. Rev. D 106, 094510 (2022), 2112.03124
2022 arXiv
-
[37]
W. Good, F. Yao, and H.-W. Lin (2025), 2505.13321
2025
-
[38]
NieMiera, W
A. NieMiera, W. Good, H.-W. Lin, and F. Yao (2025), 2510.17758
2025
-
[39]
Peardon, J
M. Peardon, J. Bulava, J. Foley, C. Morningstar, J. Dudek, R. G. Edwards, B. Joo, H.-W. Lin, D. G. Richards, and K. J. Juge (Hadron Spectrum), Phys. Rev. D80, 054506 (2009), 0905.2160
2009 arXiv
-
[40]
X. Ji, Y. Liu, A. Sch¨ afer, W. Wang, Y.-B. Yang, J.-H. Zhang, and Y. Zhao, Nucl. Phys. B964, 115311 (2021), 2008.03886
2021 arXiv
- [41]
- [42]
-
[43]
Zhang, X
J.-H. Zhang, X. Ji, A. Sch¨ afer, W. Wang, and S. Zhao, Phys. Rev. Lett.122, 142001 (2019), 1808.10824
2019 arXiv
-
[44]
Balitsky, W
I. Balitsky, W. Morris, and A. Radyushkin, Phys. Lett. B808, 135621 (2020), 1910.13963
2020 arXiv
-
[45]
Zhang, A
R. Zhang, A. V. Grebe, D. C. Hackett, M. L. Wagman, and Y. Zhao (2025), 2501.00729
2025
-
[46]
See the supplemental material for further details on the theoretical framework and data anaylsis. (2025)
2025
-
[47]
J.-W. Chen, X. Ji, and J.-H. Zhang, Nucl. Phys. B915, 1 (2017), 1609.08102
2017 arXiv
-
[48]
Izubuchi, X
T. Izubuchi, X. Ji, L. Jin, I. W. Stewart, and Y. Zhao, Phys. Rev. D98, 056004 (2018), 1801.03917
2018 arXiv
-
[49]
Alexandrou, K
C. Alexandrou, K. Cichy, M. Constantinou, K. Had- jiyiannakou, K. Jansen, H. Panagopoulos, and F. Stef- fens, Nucl. Phys. B923, 394 (2017), 1706.00265
2017 arXiv
- [50]
-
[51]
V. M. Braun, A. Vladimirov, and J.-H. Zhang, Phys. Rev. D99, 014013 (2019), 1810.00048
2019 arXiv
-
[52]
Li, Y.-Q
Z.-Y. Li, Y.-Q. Ma, and J.-W. Qiu, Phys. Rev. Lett.122, 062002 (2019), 1809.01836. 7
2019 arXiv
-
[53]
Anderson, W
T. Anderson, W. Melnitchouk, and N. Sato (2024), 2501.00665
2024
-
[54]
Zhang, J
R. Zhang, J. Holligan, X. Ji, and Y. Su, Phys. Lett. B 844, 138081 (2023), 2305.05212
2023 arXiv
-
[55]
Y. Su, J. Holligan, X. Ji, F. Yao, J.-H. Zhang, and R. Zhang, Nucl. Phys. B991, 116201 (2023), 2209.01236
2023 arXiv
-
[56]
Balitsky, W
I. Balitsky, W. Morris, and A. Radyushkin, Phys. Rev. D105, 014008 (2022), 2111.06797
2022 arXiv
-
[57]
F. Yao, Y. Ji, and J.-H. Zhang, JHEP11, 021 (2023), 2212.14415. 8 SUPPLEMENT AL MA TERIAL Theoretical framework The light-cone PDF of gluon under adjoint represen- tation is defined as xg(x, µ) = Z dz− 2πP + e−ixP +z− × ⟨P|Fa +µ(z−)U ab(z−,0)F bµ +(0)(µ)|P⟩,(9) wherez ± = 1 2 ...
2023 arXiv
Reviewed August 4, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.