Pith. sign in

REVIEW 3 major objections 3 minor 41 references

Four loop Green's functions involving the $n$ $=$ $2$ moment of the Wilson operator

T0 review · 3 major / 3 minor · reviewed 2026-08-12 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗ Review JSON

Pith's one-line read This paper claims the first four-loop Green's functions for the second-moment twist-2 Wilson operator in a quark two-point function, in all three single-scale momentum configurations, plus four-loop correlation functions of the operator…

desk verdict A solid four-loop extension of the n=2 Wilson operator matrix elements, but the new finite parts rest on an unshipped Forcer pipeline and the reported checks do not directly constrain them. read the letter →

arxiv 2411.10106 v2 pith:RL6F4YTS submitted 2024-11-15 hep-ph hep-lat

J.A. Gracey This is my paper · ORCID
classification hep-phhep-lat PACS 12.38.Bx11.10.Gh12.38.-t
keywords four-loopQCDtwist-2WilsonoperatormatrixelementsrenormalizationlatticematchingmasslessFeynmanintegralsMSbarschemeanomalousdimensions
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

The paper's aim is to obtain the first four-loop expressions, in the MSbar scheme and the massless limit, for the Green's functions obtained by inserting the second moment of the twist-2 flavour non-singlet Wilson operator into a quark two-point function, covering all three single-scale momentum configurations: zero momentum at the insertion, and the two asymmetric routings where a non-zero momentum flows through the operator and one quark leg is nullified. It also computes the four-loop correlation functions of the two gauge-invariant operators in the same scheme. These matrix elements are the continuum objects that lattice QCD determinations of parton distributions must match onto at high energy, so one more perturbative order reduces the systematic uncertainty in that matching. The author deliberately keeps a non-zero gauge parameter as a consistency check and provides the full Lorentz decomposition, which is enough to construct any lattice RI' matching scheme.

What carries the argument

The argument is carried by the automatic evaluation of massless four-loop two-point Feynman integrals in dimensional regularization, implemented in the Forcer package and driven by Qgraf for graph generation. The operator insertions are projected onto a complete basis of two symmetric traceless tensors, $T_1$ and $T_2$, so each Green's function reduces to scalar form factors; the asymmetric configurations require the upper-triangular mixing matrix between $O_1$ and $O_2$, including the total derivative operator that is physical with zero anomalous dimension. The correlation functions are decomposed into a basis of three Lorentz tensors and require additional contact renormalization constants, whose renormalization group functions are extracted from the four-loop expressions.

What would settle it

Recompute one of the new four-loop finite parts, say the SU(3) form factor $\Sigma^{(1)}_1$ at $N_f = 3$, with an independent integral-reduction algorithm and require the resulting anomalous dimension to match $\gamma_{11}$; any disagreement in the finite part would falsify the paper's four-loop claim.

Watch

Extended reading notes

Core claim

The central claim is that the complete set of scalar form factors for the operator pair $O_1 = S(\bar\psi\gamma^\mu D^\nu\psi)$ and $O_2 = S\partial^\mu(\bar\psi\gamma^\nu\psi)$ is determined to four loops in the MSbar scheme: the zero-momentum insertion reproduces the known four-loop anomalous dimension, the two asymmetric insertions determine the off-diagonal mixing renormalization constant, and the operator correlation functions fix the contact renormalization functions to the same order. These are the first four-loop results for these Green's functions, extending the three-loop forward matrix elements and three-loop correlation functions by one order. The results are analytic in the linear covariant gauge parameter and in the colour factors, with SU(3) numerical values provided.

Load-bearing premise

Everything new rests on the assumption that the automatic four-loop integral engine evaluates correctly all of the more than 32,000 Feynman graphs listed in the tables, since only the anomalous dimensions and the three-loop limits are independently verified.

Editorial extensions

If this is right

  • The zero-momentum form factors provide an independent four-loop confirmation of the non-singlet operator anomalous dimension.
  • The asymmetric form factors give the four-loop off-diagonal mixing constant $Z_{12}$ and verify $\gamma_{12} = -\tfrac{1}{2}\gamma_{11}$.
  • The gauge-invariant operator correlation functions supply four-loop matching data that do not require a gauge-fixed lattice action.
  • Setting $\alpha = 0$ in the analytic expressions yields Landau-gauge matrix elements for direct lattice comparison.
  • The full four-loop form factors allow lattice collaborations to define any RI' projection scheme and convert the result to MSbar.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • If confirmed, the same approach can be aimed at transversity and other spin-dependent twist-2 operators, which the paper names as the natural next target.
  • The four-loop relation $\gamma_{12} = -\tfrac{1}{2}\gamma_{11}$ supplies a strong consistency condition that any future five-loop computation of these operators should satisfy.
  • A lattice collaboration could combine the four-loop correlation functions with the forward matrix elements to cross-check RI' to MSbar conversion factors without re-deriving the operator renormalization.
  • Because the gauge parameter is kept non-zero throughout, one could test gauge invariance of physical combinations at four loops, a check that was previously only available at three loops.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

3 major / 3 minor

Summary. This paper presents a four-loop computation, in the MS scheme and the massless limit, of the Green's functions obtained by inserting the second moment of the twist-2 flavour nonsinglet Wilson operator into a quark two-point function, in all three single-scale momentum routings: the symmetric insertion G, and two asymmetric routings bG and eG in which momentum flows out through the operator itself. Operator mixing with the total-derivative operator O_2^{μν} is included for the asymmetric configurations. The paper also computes the correlation functions Π_ij of the two gauge-invariant operators to four loops, including contact-term renormalization. Results are presented as analytic expressions for general colour group and linear-covariant-gauge parameter, with Landau-gauge form factors in Appendix A, and as numerical SU(3) expressions. The author reproduces the four-loop anomalous dimension γ_11, verifies γ_12 = −γ_11/2 and γ_22 = 0 at four loops, reproduces the three-loop correlation functions of [34], and checks the gauge independence of the operator correlation functions. The new four-loop O(ε^0) finite parts are claimed as the first results of this order for these objects.

Significance. The computation is a direct Feynman-graph evaluation with no free parameters and no reverse-engineered constants: every quoted check anchors the poles or lower orders externally, and the manuscript is explicit about what has been verified. If correct, the four-loop finite parts are of genuine reference value for precision matching between continuum MS-scheme results and lattice determinations, and the decomposition into the Lorentz bases of (3.2) and (4.3) is appropriate for RI'-scheme conversions. Strengths include the reproduction of γ_11 against [5,6], the non-trivial verification of the mixing relations γ_12 = −γ_11/2 and γ_22 = 0, the three-loop comparison for the correlation functions, and the α-independence check. The central weakness, however, is that the genuinely new a^4 finite parts are not independently checked: the reported tests constrain pole terms and lower orders only, so a systematic error in the new master integrals or in the operator Feynman rules would pass every check listed in Sections 3 and 4. Because the asymmetric configurations lack the three-loop comparison that the correlators receive, this verification gap lands precisely on the paper's headline claim.

major comments (3)
  1. [Section 3, Eqs. (3.4)-(3.7); Section 4, Eqs. (4.8)-(4.15)] The verification reported for the new four-loop results constrains only the pole structure and lower orders, not the new O(ε^0) finite parts. Reproducing γ_11 from the 1/ε poles, verifying γ_12 = −γ_11/2 and γ_22 = 0, and reproducing the three-loop correlators of [34] do not exercise the a^4 finite terms that constitute the paper's central claim; the gauge-independence check is likewise insensitive to colour- and α-independent errors in those terms. I ask for an additional check that directly constrains the new finite parts. Concretely, one of the following would suffice: (i) an independent evaluation of a subset of the new four-loop master integrals by a second method (e.g., an alternative reduction or high-precision numerical integration); (ii) a second, independent code path for at least one of the new form factors at a^4; or (iii) a detailed three-loop comparison of the b/e form factors with refs. [26-28] together with the release of the reduction inputs. As it stands, the central claim rests on the correctness of a single reduction pipeline.
  2. [Section 3, Eqs. (3.10), (3.11); Appendix A, Eqs. (A.1)-(A.3)] No comparison with the existing three-loop non-forward results [26-28] is reported for the asymmetric (b and e) configurations, even though the introduction cites these references as the three-loop state of the art for non-forward Wilson-operator matrix elements. The b/e form factors are precisely the non-forward generalizations of [23] in which the operator Feynman rules at non-zero operator momentum and the asymmetric routing enter; a three-loop comparison would therefore test exactly the ingredients that are new at four loops. The paper explicitly reports the analogous check for the correlation functions ('We have reproduced the previous three loop expressions, [34]'), and the eΣ and bΣ expressions in the Landau gauge should similarly be checked against [26-28] where the momentum configurations overlap. This is load-bearing because the b/e form factors appear only as bare output in the main text and Appendix A.
  3. [Data Availability Statement; Section 3] The ancillary file contains the final expressions but not the computational workflow: the manuscript does not state the versions of Qgraf, Form, and Forcer used, nor does it provide the diagram-generation scripts, the operator insertion rules for non-zero operator momentum, or the reduction inputs. A bug in the user-supplied operator Feynman rules or in the reduction setup would corrupt exactly the new finite parts while leaving every reported check unchanged, since those checks are insensitive to the finite parts of the new four-loop integrals. Please either include the run scripts and input files in the ancillary material, or specify the exact software versions and enough configuration detail to allow independent reproduction of the new a^4 coefficients.
minor comments (3)
  1. [Section 4, text before Eq. (4.11)] The sentence 'γij(a) are given by (3.4), (3.6) and 3.7)' has a missing parenthesis before the last reference; it should read '(3.7)'.
  2. [Eq. (3.11)] The SU(3) a^4 coefficient of bΣ_2^(2) is printed as 2014630.937750, roughly an order of magnitude larger than the neighbouring a^4 coefficients (for comparison, bΣ_1^(2) gives 145064.413688); please confirm that this value, and the varying decimal precision of the printed coefficients (six to eight digits), coincide with the analytic expressions in the ancillary file.
  3. [Abstract and Section 1] The abstract and the body use both 'MSbar' and 'MS' for the renormalization scheme; since the finite parts of operator Green's functions are scheme-dependent, a single notation should be fixed (the body uses MS throughout).

Circularity Check

0 steps flagged · score 0.0 of 10

No significant circularity: the four-loop results are direct Forcer evaluations, benchmarked against external anomalous dimensions and three-loop checks.

full rationale

The paper's central claim is the four-loop MSbar evaluation of operator matrix elements and correlation functions for the n = 2 twist-2 Wilson operator. These are obtained by applying the external Forcer package to Qgraf-generated diagrams, with no fitted parameters and no quantity defined in terms of the result it is said to derive. The verification steps cited in the manuscript reproduce the four-loop anomalous dimension γ11 from the independent results of refs. [1-6], check the relations γ12 = -1/2 γ11 and γ22 = 0, recover the three-loop correlation functions of [34], and confirm gauge-parameter independence of the correlators. These checks anchor the computation to published, independent data. The author's self-citations ([23], [26], [27], [28], [34]) are used as benchmarks, notation conventions, or earlier three-loop context, not as the provenance of the new a^4 finite parts. The skeptic's concern that a bug in Forcer or in the implementation could corrupt the new four-loop finite parts while passing all checks is a verification and correctness risk, not a circularity: the derivation is not equivalent to its inputs by construction. Accordingly, no circular step is identified.

Assumptions & free parameters 0 free parameters · 3 assumptions · 0 invented entities

The central claim relies on the correctness of the automated Forcer pipeline, on the assumed completeness of the two-operator basis, and on the standard dimensional regularization and MSbar framework. No free parameters are fitted and no new entities are introduced.

assumptions (3)
  • domain assumption Forcer correctly evaluates the four-loop massless two-point integrals in the generated graphs.
    Section 3 assigns all integral evaluation to Forcer; no independent check of the new four-loop finite parts is provided.
  • domain assumption The two-operator basis O1, O2 is complete and the renormalization mixing matrix is upper triangular with Z22 = 1.
    Section 2 asserts completeness from integration by parts and verifies gamma12 = -gamma11/2 and gamma22 = 0; a hidden additional mixing operator would change the extracted renormalization constants.
  • standard math Dimensional regularization in d = 4 - 2 epsilon with MSbar subtraction is the correct framework for these Green's functions.
    Standard tool of the field; all expressions are series in epsilon with MSbar counterterms.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Four loop Green's functions involving the $n$ $=$ $2$ moment of the Wilson operator." pith.science (2026). https://pith.science/paper/RL6F4YTS

@misc{pith2026241110106,
  author       = {Pith},
  title        = {Pith review of: Four loop Green's functions involving the $n$ $=$ $2$ moment of the Wilson operator},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/RL6F4YTS}},
  note         = {Machine review of arXiv:2411.10106}
}
abstract

We evaluate the Green's function for the insertion of the second moment of the twist-$2$ flavour nonsinglet Wilson operator in a quark $2$-point function in all three different single scale external momentum configurations at four loops in the MSbar scheme and the chiral limit. One configuration is where the operator is inserted at zero momentum while the other two are where a non-zero momentum flows out through the operator itself with one external quark momentum nullified. In the latter two configurations mixing of the operator with a total derivative twist-$2$ operator is included for renormalization group consistency. In addition we compute the correlation functions of both gauge invariant operators to four loops in the same scheme.

Figures

Figures reproduced from arXiv: 2411.10106 by the authors.

Figure 1
Figure 1. Operator correlation function Πµ1µ2ν1ν2 ij (q). 4 Operator correlation functions. The second class of Green’s function that we examine is that concerning the correlation function of the two gauge invariant operators O µν i defined as Π µ1µ2ν1ν2 ij (q) = (4π) 2 i Z d dx eiqx⟨0|Oµ1µ2 i (x)O ν1ν2 j (0)|0⟩ (4.1) where q is the momentum with q 2 = − Q2 . The correlation function is formally illustrated in [PITH_FULL_IMA… view at source ↗

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

41 extracted references · 33 canonical work pages

  1. [34]

    Gracey, Eur

    J.A. Gracey, Eur. Phys. J. C83 (2023), 181

    work page 2023

  2. [23]

    G¨ ockeler, R

    M. G¨ ockeler, R. Horsley, Y. Nakamura, H. Perlt, D. Pleiter, P.E.L. Rakow, A. Sch¨ afer, G. Schierholz, A. Schiller, H. St¨ uben & J. Zanotti, Phys. Rev. D82 (2010), 114511; Phys. Rev. D86 (2012), 099903

    work page 2010

  3. [1]

    This basis choice is complete and sufficient for the lattice to construct a specific RI ′ scheme for the operator matrix element measurement itself

    As Forcer requires Lorentz scalar Feynman graphs the initial step is to decompose (2.10) into a complete tensor basis and to fix our notation we set Gµν i (p) = 2X j=1 Σ(j) i (p)T µν j (p) 5 bGµν i (p) = 2X j=1 bΣ(j) i (p)T µν j (p) eGµν i (p) = 2X j=1 eΣ(j) i (p)T µν j (p) (3.1) for i = 1 and 2 where the two basis tensors are, [23], T µν 1 (p) = γµpν + γ...

    work page arXiv 1984

  4. [2]

    Gross & F.J

    D.J. Gross & F.J. Wilczek, Phys. Rev. D9 (1974), 980

    work page 1974

  5. [3]

    Floratos, D.A

    E.G. Floratos, D.A. Ross & C.T. Sachrajda, Nucl. Phys. B129 (1977), 66; B139 (1978), 545(E)

    work page 1977

  6. [4]

    Moch, J.A.M

    S. Moch, J.A.M. Vermaseren & A. Vogt, Nucl. Phys. B688 (2004), 101

    work page 2004

  7. [5]

    Bl¨ umlein, P

    J. Bl¨ umlein, P. Marquard, C. Schneider & K. Sch¨ onwald, Nucl. Phys.B971 (2021), 115542

    work page 2021

  8. [6]

    Baikov & K.G

    P.A. Baikov & K.G. Chetyrkin, Nucl. Phys. Proc. Suppl. 160 (2006), 76

    work page 2006

Show all 41 references
  1. [7]

    Velizhanin, Nucl

    V.N. Velizhanin, Nucl. Phys. B860 (2012), 288

    2012

  2. [8]

    Moch, J.A.M

    S. Moch, J.A.M. Vermaseren & A. Vogt, Phys. Lett. B825 (2022), 136853

    2022

  3. [9]

    S. Moch, B. Ruijl, T. Ueda, J.A.M. Vermaseren & A. Vogt, Phys. Lett. B849 (2024), 138468

    2024

  4. [10]

    Davies, A

    J. Davies, A. Vogt, B. Ruijl, T. Ueda & J.A.M. Vermaseren, Nucl. Phys. B915 (2017), 335

    2017

  5. [11]

    Gracey, Phys

    J.A. Gracey, Phys. Lett. B322 (1994), 141

    1994

  6. [12]

    Herzog, S

    F. Herzog, S. Moch, B. Ruijl, T. Ueda, J.A.M. Vermaseren & A. Vogt, Phys. Lett. B790 (2019), 436

    2019

  7. [13]

    T. Ueda, B. Ruijl & J.A.M. Vermaseren, PoS LL2016 (2016), 070

    2016

  8. [14]

    T. Ueda, B. Ruijl & J.A.M. Vermaseren, Comput. Phys. Commun. 253 (2020), 107198

    2020

  9. [15]

    Vermaseren, math-ph/0010025

    J.A.M. Vermaseren, math-ph/0010025

    arXiv

  10. [16]

    Tentyukov & J.A.M

    M. Tentyukov & J.A.M. Vermaseren, Comput. Phys. Commun. 181 (2010), 1419

    2010

  11. [17]

    G¨ ockeler, R

    M. G¨ ockeler, R. Horsley, D. Pleiter, P.E.L. Rakow & G. Schierholz, Phys. Rev. D71 (2005), 114511

    2005

  12. [18]

    Chambers, R

    A.J. Chambers, R. Horsley, Y. Nakamura, H. Perlt, P.E.L. Rakow, G. Schierholz, A. Schiller, K. Somfleth, R.D. Young & J.M. Zanotti, Phys. Rev. Lett. 118 (2017), 242000. 25

    2017

  13. [19]

    Harris, G

    T. Harris, G. von Hippel, P. Junnarkar, H.B. Meyer, K. Ottnad, J. Wilhelm, H. Wittig & L. Wrang, Phys. Rev. D100 (2019), 034513

    2019

  14. [20]

    Mondal, R

    S. Mondal, R. Gupta, S. Park, B. Yoon, T. Bhattacharya & H.-W. Lin, Phys. Rev. D102 (2020), 054512

    2020

  15. [21]

    Gao, A.D

    X. Gao, A.D. Hanlon, J. Holligan, N. Karthik, S. Mukherjee, P. Petreczky, S. Syritsyn & Y. Zhao, Phys. Rev. D107 (2023), 074509

    2023

  16. [22]

    Rodekamp, M

    M. Rodekamp, M. Engelhardt, J.R. Green, S. Krieg, S. Liuti, S. Meinel, J.W. Negele, A. Pochinsky & S. Syritsyn, Phys. Rev. D109 (2024), 074508

    2024

  17. [24]

    Gracey, Nucl

    J.A. Gracey, Nucl. Phys. B667 (2003), 242

    2003

  18. [25]

    Gracey, J

    J.A. Gracey, J. High Energy Phys. 10 (2006), 040

    2006

  19. [26]

    Gorishny, S.A

    S.G. Gorishny, S.A. Larin, L.R. Surguladze & F.K. Tkachov, Comput. Phys. Commun. 55 (1989), 381

    1989

  20. [27]

    Gracey, J

    J.A. Gracey, J. High Energy Phys. 03 (2011), 109

    2011

  21. [28]

    Gracey, Phys

    J.A. Gracey, Phys. Rev. D84 (2011), 016002

    2011

  22. [29]

    Kniehl & O.L

    B.A. Kniehl & O.L. Veretin, Phys. Lett. B804 (2020), 13598

    2020

  23. [30]

    Martinelli, C

    G. Martinelli, C. Pittori, C.T. Sachrajda, M. Testa & A. Vladikas, Nucl. Phys. B445 (1995), 81

    1995

  24. [31]

    Franco & V

    E. Franco & V. Lubicz, Nucl. Phys. B531 (1998), 641

    1998

  25. [32]

    Chetyrkin & A

    K.G. Chetyrkin & A. R´ etey, Nucl. Phys.B583 (2000), 3

    2000

  26. [33]

    Chetyrkin & A

    K.G. Chetyrkin & A. R´ etey, hep-ph/0007088

    arXiv

  27. [35]

    Gracey, J

    J.A. Gracey, J. High Energy Phys. 04 (2009), 127

    2009

  28. [36]

    ’t Hooft, Nucl

    G. ’t Hooft, Nucl. Phys. B61 (1973), 455

    1973

  29. [37]

    Nogueira, J

    P. Nogueira, J. Comput. Phys. 105 (1993), 279

    1993

  30. [38]

    van Ritbergen, A.N

    T. van Ritbergen, A.N. Schellekens & J.A.M. Vermaseren, Int. J. Mod. Phys. A14 (1999), 41

    1999

  31. [39]

    Moch & S

    S. Moch & S. van Thurenhout, Nucl. Phys. B971 (2021), 115536

    2021

  32. [40]

    Chetyrkin, Phys

    K.G. Chetyrkin, Phys. Lett. B390 (1997), 309

    1997

  33. [41]

    Chetyrkin, J.H

    K.G. Chetyrkin, J.H. K¨ uhn & A. Kwiatkowski, Phys. Rept. 277 (1996), 189. 26

    1996

Pith tools

Reviewed August 12, 2026 · model on record in the stance chip above.