REVIEW 3 major objections 5 minor 30 references
This paper establishes three-loop RI'/SMOM-to-MS conversion factors for bilinear lattice QCD operators and two-loop matching plus three-loop anomalous dimensions for three-quark operators, enabling improved lattice studies of baryon distrib
Reviewed by Pith at T0; open to challenge. T0 means a machine referee read the full paper against a public rubric. the ladder, T0–T4 →
T0 review · deepseek-v4-flash
2026-08-02 03:48 UTC pith:62UHK3XN
load-bearing objection A competent review of the authors' own RI/SMOM work; trust the bilinear parts, but the three-quark N=1/2 headline results are only referenced to an unpublished manuscript and cannot be checked. the 3 major comments →
RI/SMOM renormalization for lattice QCD: bilinear and three-quark operators
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
The central claim is that the conversion between RI'/SMOM and MS is now known to three loops for all nonsinglet bilinear operators of phenomenological interest (scalar, vector, tensor currents; twist-two moments n=2,3) and to two loops for the N=0,1 Mellin moments of three-quark operators, with three-loop anomalous dimensions for N=0,1,2. The three-quark results are made possible by renormalizing with open spinor indices: the four-point function is decomposed into 247 independent tensor structures at two loops, renormalized in MS in d dimensions, and only then is d set to 4. This avoids gamma5 ambiguities and makes evanescent operators vanish at d=4. All numerical values are quoted in Landau
What carries the argument
The load-bearing construction is the symmetric subtraction point p_i^2 = -mu^2 with no exceptional momentum, which defines the RI'/SMOM scheme and improves infrared behaviour relative to zero-momentum RI/MOM. Around this point the amputated Green functions are decomposed into scalar form factors; for bilinear operators these are expanded through three loops with IBP-reduced master integrals evaluated by sector decomposition. For three-quark operators, the essential mechanism is the open-spinor-index renormalization: uncontracted spinor structures are renormalized in MS, with a 64- and 247-dimensional form-factor basis at one and two loops in d=4, so evanescent operators are automatically pro
Load-bearing premise
The quoted three-quark matching factors and anomalous dimensions assume that renormalizing with open spinor indices and only setting d=4 at the end defines the same MS scheme as other evanescent-operator prescriptions; if that equivalence fails, the numerical values shift.
What would settle it
Recompute the N=0 three-quark form factor f1 at two loops in an alternative scheme that contracts spinor indices before renormalization and handles gamma5 with a definite d-dimensional prescription; any disagreement with the quoted R_xi-gauge coefficients beyond the stated numerical accuracy would disprove the paper's matching factors.
If this is right
- Lattice determinations of quark masses and low moments of structure functions can include three-loop RI'/SMOM-to-MS conversion, cutting a known systematic uncertainty.
- Baryon distribution amplitudes can now be extracted from lattice correlation functions at two-loop matching accuracy with three-loop evolution for the first two Mellin moments.
- The analytic R_xi-gauge results allow gauge-parameter independence checks of the MS conversion factors.
- The open-spinor-index basis for three-quark operators provides a template for extending the calculation to N=2 matching.
Where Pith is reading between the lines
- One testable extension is to complete the two-loop matching for N=2 and compare the resulting baryon DA moments against sum-rule determinations; the paper flags this as future work.
- Because the bilinear three-loop corrections are numerically large, an analogous size for three-quark three-loop anomalous dimensions suggests that final precision will hinge on the still-missing higher-order matching.
- The same open-spinor-index treatment could be carried over to four-quark operators, where evanescent operators also mix and currently limit MS conversion accuracy.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper is a conference-style review of perturbative RI'/SMOM-to-MS conversion factors relevant for lattice QCD. For bilinear operators, it summarizes three-loop results for local scalar, vector, and tensor currents and for the n=2,3 twist-two moments of nonsinglet structure functions, quoting selected numerical coefficients in Landau gauge and referring to the authors' Refs. [27,28] for the full form factors. For three-quark operators relevant to baryon distribution amplitudes, it summarizes two-loop RI'/SMOM matching for the N=0 Mellin moment, three-loop anomalous dimensions for N=0,1,2, and two-loop matching for N=1, with the N=0 results and a sample form factor shown explicitly in Eqs. (21),(24) and the N=1,2 results attributed to Ref. [30], described as a manuscript in preparation. The stated central claim is that the three-quark results advance the state of the art for lattice studies of baryon distribution amplitudes.
Significance. If the results are correct, the paper summarizes useful perturbative input for lattice QCD: three-loop bilinear conversion factors reduce systematic uncertainties in quark masses and structure-function moments, and the three-quark N=1 matching plus three-loop anomalous dimensions would allow two-loop-accurate conversion and three-loop evolution for baryon-DA operators. The bilinear part is well supported: the three-loop scalar current result agrees with the independent analytic calculation of Ref. [26], the lower-loop coefficients reproduce Refs. [2,3,6,7], and the numerical integrations are performed with established tools. The N=0 three-quark result has a nontrivial check through agreement with Refs. [14,15] after evanescent-operator finite renormalization. However, the paper's headline advance for N=1,2 three-quark operators is not inspectable: no N=1 matching coefficients and no N=1,2 anomalous-dimension coefficients are printed, and the supporting reference is unpublished. This is the decisive weakness for the evaluation of the paper's strongest claim.
major comments (3)
- [§3.4 and §3.5] The central new results — two-loop RI'/SMOM matching for N=1 and three-loop anomalous dimensions for N=1,2 — are attributed entirely to Ref. [30], cited as 'manuscript in preparation.' No coefficient, formula, or derivation for these quantities appears in the paper. A review may cite unpublished work, but the abstract and §4 make an advance-of-the-state-of-the-art claim resting on results the reader cannot inspect. As a referee, I cannot verify correctness, convention choices, or even the existence of the quoted numbers. This is load-bearing because it is the basis for the paper's main conclusion about baryon DAs. The authors should either include the N=1,2 results (or a representative set with a clear derivation/reference to a published or arXiv-published companion), or explicitly mark these parts as announcements of work in progress and soften the abstract and conclusions accordingly.
- [§3.2 and §3.5] The open-spinor-index renormalization procedure is checked for N=0: §3.5 states that after rewriting evanescent structures and finite renormalizations, full agreement with Refs. [14,15] is found. No analogous check is shown for N=1,2. The finite parts of evanescent-operator mixing define a scheme; an MS-like convention for these operators is not uniquely fixed without specifying those finite parts. If Ref. [30] uses a convention different from the one assumed by lattice practitioners (e.g., the operator basis of Refs. [13,17-19]), the quoted N=1 matching and N=1,2 anomalous dimensions would shift at finite order. This is a correctness-risk concern rather than an internal inconsistency, but it is directly relevant to the paper's usability claim. The authors should state the convention explicitly, give the finite renormalization that maps to a standard basis, or provide a nontrivial check
- [§4] The concluding sentence says 'All results discussed here are available in machine-readable form in the respective publications.' This is not currently true for the central three-quark N=1,2 results, which are only in the unpublished Ref. [30]. The paper itself contains no machine-readable supplementary material. For a review whose contribution is to collect results, this statement overstates reproducibility and should be corrected or made conditional on the publication of [30].
minor comments (5)
- [§2.1, Eq. (7)] Typo: 'interprested' should be 'interpreted.'
- [§3.1, Eq. (18)] The third term in the operator definition reads '(D_{\sigma_1}...D_{\sigma_s} s_{b2})'; the color index appears to be b2 but should likely be b3 to match the epsilon tensor. Please check and correct.
- [§3.1, Eq. (16)] The notation mixes quark flavors u,d,s with generic spinor indices; the flavor labels are not carried through the local operators in Eq. (18). A sentence clarifying the flavor structure would help avoid ambiguity.
- [§3.4, Eq. (25)] The decomposition H_{\mu,j} = sum p_{\mu,j} T_n f_{n,j} is written as an equality after 'in the leading twist approximation'; the qualifier should be repeated in the equation or placed directly before it, since Eq. (20) for N=0 was an exact decomposition.
- [§1, Refs. [2-8]] References in the introduction are grouped as [2-8], but some of these refer to three-loop and three-quark results; a more explicit citation breakdown would help readers locate the specific bilinear one- and two-loop papers.
Circularity Check
N=1/2 three-quark matching and anomalous dimensions rest solely on an unpublished self-citation; bilinear and N=0 results are independently checked.
specific steps
-
self citation load bearing
[Sections 3.4–3.5, after Eqs. (25) and (28), and Ref. [30]]
"The amputated four-point functions for N=1 three-quark operators with open spinor indices, H..., is evaluated at two loops in RI'/SMOM kinematics (19) [30]. ... The MS anomalous dimensions ... have been computed through three loops [30]."
The paper's strongest claim — 'three-loop anomalous dimensions for N=0,1,2 together with two-loop matching for N=1 advance the state of the art' — depends on the N=1 matching and N=1,2 anomalous dimensions, but these are attributed only to Ref. [30], a 'manuscript in preparation' by the same two authors. No N=1 matching coefficient or N=1,2 anomalous-dimension coefficient is printed or independently checked; the only agreement quoted (with Refs. [14,15]) is for N=0. The load-bearing step is thus an unverifiable self-citation rather than an external benchmark, and the conclusion that results are 'available in machine-readable form in the respective publications' is not yet true for these central new results.
full rationale
The bilinear section is not circular: it reports parameter-free three-loop computations with lower-loop agreement against Refs. [2,3,6] and analytic three-loop confirmation by Bednyakov–Pikelner (Ref. [26]). The three-quark N=0 content is also anchored externally: the open-spinor-index method is taken from independent Ref. [14], and the N=0 anomalous dimension is stated to agree with Refs. [14,15]. No fitted input is renamed as a prediction anywhere in the paper. The only genuine load-bearing self-citation is Ref. [30], which supplies the N=1 two-loop matching and N=1,2 three-loop anomalous dimensions and exists only as 'manuscript in preparation'. This is a citation burden, not a by-construction equivalence; there is no equation in the paper that reduces to its own input. The possible scheme-convention sensitivity of the open-index MS prescription for N=1,2 is an unvalidated assumption rather than a demonstrated circular step, so it is noted here (Section 3.2, 'Following Ref. [14], we renormalize three-quark operators without contracting spinor indices') but not counted as a second circularity. Score 4 reflects: some load-bearing self-citation for part of the strongest claim, while bilinear and N=0 results retain independent content and external checks.
Axiom & Free-Parameter Ledger
axioms (4)
- standard math Dimensional regularization in d=4-2ε with MS renormalization is a valid scheme for these operators.
- domain assumption The RI'/SMOM subtraction conditions, Eqs. (2) and (19), define a scheme whose difference from MS is computable perturbatively.
- domain assumption Open-spinor-index renormalization following Ref. [14] correctly handles evanescent operators and γ5 ambiguities.
- domain assumption Numerical evaluation of master integrals by sector decomposition and Monte Carlo is accurate to ~1e-4–1e-5 after cancellations.
read the original abstract
We review perturbative matching between the regularization-invariant symmetric MOM (RI/SMOM) and MS schemes at the symmetric subtraction point, which is relevant for lattice QCD simulations. For bilinear operators we summarize three-loop conversion factors for local quark currents and for the n=2,3 twist-two moments of structure functions. For three-quark operators we summarize two-loop RI/SMOM matching for N=0 baryonic operators and three-loop anomalous dimensions together with two-loop matching for the N=1 Mellin moment, enabling improved lattice studies of baryon distribution amplitudes. All numerical results quoted here are given in Landau gauge.
Figures
Reference graph
Works this paper leans on
-
[1]
8 RI′/SMOM renormalization for lattice QCD: bilinear and three-quark operators B.A
G.Martinelli, C.Pittori, C.T.Sachrajda,M.TestaandA.Vladikas, Nucl.Phys.B445(1995) 81 [hep-lat/9411010]. 8 RI′/SMOM renormalization for lattice QCD: bilinear and three-quark operators B.A. Kniehl, O.L. Veretin
Pith/arXiv arXiv 1995
-
[2]
C. Sturmet al., Phys. Rev. D80(2009) 014501 [arXiv:0901.2599 [hep-ph]]
Pith/arXiv arXiv 2009
-
[3]
J. A. Gracey and R. M. Moss, Phys. Lett. B686(2010) 189 [arXiv:1004.1758 [hep-ph]]
Pith/arXiv arXiv 2010
-
[4]
M. Gorbahn and S. Komatsu, JHEP09(2010) 019 [arXiv:1006.3858 [hep-ph]]
Pith/arXiv arXiv 2010
-
[5]
L. G. Almeida and C. Sturm, Phys. Rev. D82(2010) 054017 [arXiv:1004.4613 [hep-ph]]
Pith/arXiv arXiv 2010
-
[6]
J. A. Gracey, Phys. Rev. D83(2011) 116006 [arXiv:1104.2105 [hep-ph]]
Pith/arXiv arXiv 2011
-
[7]
J. A. Gracey, Phys. Rev. D84(2011) 054016 [arXiv:1106.5784 [hep-ph]]
Pith/arXiv arXiv 2011
-
[8]
J. A. Gracey, Phys. Rev. D85(2012) 054003 [arXiv:1112.1784 [hep-ph]]
Pith/arXiv arXiv 2012
-
[9]
A. V. Efremov and A. V. Radyushkin, Theor. Math. Phys.44(1980) 573
1980
-
[10]
G. P. Lepage and S. J. Brodsky, Phys. Rev. D22(1980) 2157
1980
-
[11]
V. L. Chernyak and A. R. Zhitnitsky, Phys. Rept.112(1984) 173
1984
-
[12]
V. M. Braun, R. J. Fries, N. Mahnke and E. Stein, Nucl. Phys. B589(2000) 381 [hep- ph/0007279]
arXiv 2000
-
[13]
V. M. Braun and A. N. Manashov, Nucl. Phys. B542(1999) 91 [hep-ph/9805225]
Pith/arXiv arXiv 1999
-
[14]
S. Krankl, C. A. Lee, A. V. Manohar and M. J. Savage, Nucl. Phys. B856(2012) 83 [arXiv:1107.3718 [hep-ph]]
Pith/arXiv arXiv 2012
-
[15]
J. A. Gracey, Phys. Rev. D86(2012) 034018 [arXiv:1208.5619 [hep-ph]]
Pith/arXiv arXiv 2012
-
[16]
J. A. Gracey, Phys. Rev. D113(2026) 065009 [arXiv:2510.21940 [hep-ph]]
arXiv 2026
-
[17]
G. S. Baliet al., JHEP02(2016) 070 [arXiv:1512.02050 [hep-lat]]
Pith/arXiv arXiv 2016
-
[18]
RQCD Collaboration, Eur. Phys. J. A55(2019) 116 [arXiv:1903.12590 [hep-lat]]
Pith/arXiv arXiv 2019
-
[19]
G. S. Baliet al., Phys. Rev. D111(2025) 094517 [arXiv:2411.19091 [hep-lat]]
Pith/arXiv arXiv 2025
-
[20]
M. Gruber, V. M. Braun, A. V. Manohar and M. J. Savage, Phys. Rev. D97(2018) 034516 [arXiv:1710.07157 [hep-ph]]
Pith/arXiv arXiv 2018
-
[21]
K. G. Chetyrkin and F. V. Tkachov, Nucl. Phys. B192(1981) 159
1981
-
[22]
A. V. Smirnov, Comput. Phys. Commun.189(2015) 182 [arXiv:1408.2372 [hep-ph]]
Pith/arXiv arXiv 2015
-
[23]
K. G. Chetyrkin and M. Steinhauser, Phys. Lett. B639(2006) 263 [hep-ph/0604040]
Pith/arXiv arXiv 2006
-
[24]
A. V. Smirnov, Comput. Phys. Commun.199(2016) 189 [arXiv:1511.03614 [hep-ph]]
Pith/arXiv arXiv 2016
- [25]
-
[26]
9 RI′/SMOM renormalization for lattice QCD: bilinear and three-quark operators B.A
A.BednyakovandA.Pikelner,Phys.Rev.D101(2020)091501[arXiv:2002.12758[hep-ph]]. 9 RI′/SMOM renormalization for lattice QCD: bilinear and three-quark operators B.A. Kniehl, O.L. Veretin
Pith/arXiv arXiv 2020
-
[27]
B.A.KniehlandO.L.Veretin,Phys.Lett.B804(2020)135398[arXiv:2002.10894[hep-ph]]
Pith/arXiv arXiv 2020
-
[28]
B.A.KniehlandO.L.Veretin,Nucl.Phys.B961(2020)115229[arXiv:2009.11325[hep-ph]]
Pith/arXiv arXiv 2020
-
[29]
B.A.KniehlandO.L.Veretin,Nucl.Phys.B992(2023)116210[arXiv:2207.08553[hep-ph]]
Pith/arXiv arXiv 2023
-
[30]
B. A. Kniehl and O. L. Veretin, Renormalization of three-quark operators with up to two derivatives at three loops, manuscript in preparation. 10
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.