pith. sign in

arxiv: 2606.28035 · v1 · pith:GSLHUAQ3new · submitted 2026-06-26 · ✦ hep-lat · astro-ph.CO· hep-ph· hep-th

The QCD energy-momentum tensor on the lattice: non-perturbative renormalization with N_f=3

Matteo Bresciani , Mattia Dalla Brida , Leonardo Giusti , Mitsuaki Hirasawa , Michele Pepe , Luca Virz\`i This is my paper

Pith reviewed 2026-06-29 02:00 UTC · model grok-4.3

classification ✦ hep-lat astro-ph.COhep-phhep-th
keywords lattice QCDenergy-momentum tensornon-perturbative renormalizationWard identitiesshifted boundary conditionsNf=3clover discretization
0
0 comments X

The pith

The QCD energy-momentum tensor on the lattice is renormalized non-perturbatively for three flavors by fixing constants via discretized Ward identities.

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

The paper constructs the traceless components of the energy-momentum tensor on the lattice for QCD with three quark flavors so that their correlation functions satisfy the appropriate continuum Ward identities. The hypercubic lattice symmetry splits the SO(4) nonet into triplet and sextet representations, each requiring independent gluonic and fermionic renormalization constants. These constants are fixed non-perturbatively through Monte Carlo simulations that impose discretized versions of the Ward identities on one-point functions using shifted boundary conditions and an imaginary chemical potential. The resulting constants achieve a few percent accuracy for bare couplings in the range 0 to 0.96, enabling consistent continuum extrapolations for physical observables built from the tensor.

Core claim

The traceless components of the energy-momentum tensor are constructed on the lattice for Nf=3 QCD using the Wilson plaquette action for gluons and O(a)-improved Wilson action for quarks. Gluonic parts employ the clover discretization of the field strength tensor while fermionic parts use symmetric covariant derivatives. For each of the resulting triplet and sextet representations under the hypercubic group, the two independent renormalization constants are determined non-perturbatively by imposing discretized continuum Ward identities on one-point correlation functions with shifted boundary conditions and imaginary chemical potential, with Monte Carlo simulations yielding final accuracy of

What carries the argument

Discretized continuum Ward identities imposed on one-point functions with shifted boundary conditions and imaginary chemical potential to fix the two independent renormalization constants per multiplet.

Load-bearing premise

The chosen clover discretization, symmetric covariant derivatives, and boundary conditions suffice to determine the renormalization constants without residual lattice artifacts that would prevent the operators from satisfying continuum Ward identities.

What would settle it

Renormalized one-point functions that fail to match the expected continuum Ward identity values after taking the lattice spacing to zero at fixed physical volume would show the constants are not correctly fixed.

read the original abstract

We construct the traceless components of the energy-momentum tensor on the lattice for QCD with $N_f=3$ flavours, such that their correlation functions satisfy the appropriate Ward identities in the continuum limit. To carry out this program, we define the theory on the lattice by the Wilson-plaquette and the $O(a)$-improved Wilson actions for gluons and quarks respectively. The discretization of the space-time entails that (i) the irreducible nonet representation of the SO($4$) group splits into a triplet and a sextet irreducible representations of the hypercubic group, and (ii) for each multiplet non-perturbative determinations of the the gluonic and fermionic renormalization constants are required. The bare gluonic components of the energy-momentum tensor are defined via the clover discretization of the field strength tensor, while the fermionic ones are discretized by appropriate combinations of symmetric covariant derivatives. Either for the triplet or the sextet representations, the two independent renormalization constants are then fixed non-perturbatively by imposing discretized versions of continuum Ward identities for one-point correlation functions in the presence of shifted boundary conditions and an imaginary chemical potential. The non-perturbative calculation is then carried out by Monte Carlo simulations, and the resulting renormalization constants are determined with a final accuracy of a few percent for values of the bare coupling constant squared in the range $0 \leq g_0^2\leq 0.96$.

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 paper constructs the traceless components of the QCD energy-momentum tensor on the lattice for N_f=3 using the Wilson plaquette action for gluons and O(a)-improved Wilson fermions. It splits the nonet into triplet and sextet representations under the hypercubic group, defines bare gluonic operators via the clover discretization and fermionic ones via symmetric covariant derivatives, and determines the two independent renormalization constants per multiplet non-perturbatively by imposing discretized continuum Ward identities on one-point functions with shifted boundary conditions and imaginary chemical potential. Monte Carlo simulations yield the constants to a few-percent accuracy for 0 ≤ g₀² ≤ 0.96.

Significance. If the results hold, the work supplies non-perturbative renormalization constants essential for lattice studies of the EMT in three-flavor QCD, supporting calculations of thermodynamic observables, transport coefficients, and matrix elements without perturbative matching uncertainties. The non-perturbative fixing via Ward identities on one-point functions with shifted boundaries is a technically sound approach that avoids circularity.

major comments (1)
  1. [renormalization procedure (abstract and method outline)] The central claim of few-percent accuracy for the renormalization constants rests on the assumption that the clover discretization for gluonic components and symmetric covariant derivatives for fermionic components, together with the chosen shifted boundary conditions and imaginary chemical potential, suffice to eliminate residual O(a) artifacts that would affect the continuum matching. This assumption is load-bearing but receives limited explicit validation in the provided description (e.g., no detailed continuum extrapolations or comparisons at weak coupling are referenced).
minor comments (1)
  1. [Monte Carlo simulations] The range 0 ≤ g₀² ≤ 0.96 includes relatively strong couplings; it would be useful to clarify how the Monte Carlo statistics and autocorrelation times behave in this regime to support the quoted accuracy.

Simulated Author's Rebuttal

1 responses · 0 unresolved

We thank the referee for their thorough review and positive evaluation of the significance of our work. We provide a point-by-point response to the major comment below.

read point-by-point responses

  1. Referee: [renormalization procedure (abstract and method outline)] The central claim of few-percent accuracy for the renormalization constants rests on the assumption that the clover discretization for gluonic components and symmetric covariant derivatives for fermionic components, together with the chosen shifted boundary conditions and imaginary chemical potential, suffice to eliminate residual O(a) artifacts that would affect the continuum matching. This assumption is load-bearing but receives limited explicit validation in the provided description (e.g., no detailed continuum extrapolations or comparisons at weak coupling are referenced).

    Authors: The non-perturbative renormalization is performed by imposing discretized Ward identities that become exact in the continuum limit. The clover discretization of the field strength and the use of symmetric covariant derivatives are chosen to be consistent with the O(a)-improved action, thereby reducing leading discretization errors. The shifted boundary conditions with imaginary chemical potential provide a practical way to compute the one-point functions without introducing additional systematic errors beyond those controlled by the lattice spacing. The few-percent accuracy quoted refers to the statistical uncertainty from the Monte Carlo simulations. While the manuscript focuses on the method and results, we acknowledge that more explicit discussion of residual O(a) effects would strengthen the presentation. We will add a dedicated subsection discussing the expected size of artifacts and comparisons to perturbative expectations at weak coupling in the revised version. revision: yes

Circularity Check

0 steps flagged

No significant circularity; renormalization fixed by external continuum Ward identities

full rationale

The paper determines the two renormalization constants per multiplet by imposing discretized versions of continuum Ward identities on lattice one-point functions (with shifted boundaries and imaginary chemical potential). This is a standard non-perturbative renormalization procedure in which the continuum symmetries serve as external input to fix lattice Z factors; the output is not a fit or redefinition of the same data used to define the constants. No self-citation is load-bearing, no ansatz is smuggled, and no uniqueness theorem from prior author work is invoked to force the result. The derivation chain is self-contained against the continuum Ward identities and does not reduce by construction to its inputs.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The central claim rests on the assumption that the lattice theory with the chosen actions and discretizations approaches continuum QCD in the limit a→0, and that the discretized Ward identities correctly isolate the renormalization constants. No new entities are postulated; the renormalization constants themselves are the computed outputs rather than free parameters.

axioms (2)
  • domain assumption The continuum limit of the lattice-regularized theory with Wilson plaquette and O(a)-improved Wilson actions reproduces continuum QCD.
    Invoked when stating that the correlation functions must satisfy continuum Ward identities after renormalization.
  • domain assumption The hypercubic symmetry splitting of the SO(4) nonet into triplet and sextet representations is correctly accounted for by separate renormalization constants for each.
    Stated explicitly in the abstract as a consequence of the lattice discretization.

pith-pipeline@v0.9.1-grok · 5822 in / 1547 out tokens · 30358 ms · 2026-06-29T02:00:33.208900+00:00 · methodology

discussion (0)

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

Reference graph

Works this paper leans on

44 extracted references · 1 canonical work pages

  1. [1]

    Weinberg,Gravitation and cosmology: principles and applications of the general theory of relativity, John Wiley & Sons (1972)

    S. Weinberg,Gravitation and cosmology: principles and applications of the general theory of relativity, John Wiley & Sons (1972)

    1972

  2. [2]

    Polyakov and P

    M.V. Polyakov and P. Schweitzer,Forces inside hadrons: pressure, surface tension, mechanical radius, and all that,Int. J. Mod. Phys. A33(2018) 1830025 [1805.06596]

    Pith/arXiv arXiv 2018

  3. [3]

    Landau and E

    L. Landau and E. Lifshitz,Course of Theoretical Physics, Vol. 6: Fluid Mechanics(1987)

    1987

  4. [4]

    Kubo,Statistical mechanical theory of irreversible processes

    R. Kubo,Statistical mechanical theory of irreversible processes. 1. General theory and simple applications in magnetic and conduction problems,J. Phys. Soc. Jap. 12(1957) 570

    1957

  5. [5]

    R. Kubo, M. Yokota and S. Nakajima,Statistical-Mechanical Theory of Irreversible Processes. II. Response to Thermal Disturbance,J. Phys. Soc. Jap. 12(1957) 1203

    1957

  6. [6]

    Caracciolo, G

    S. Caracciolo, G. Curci, P. Menotti and A. Pelissetto,The Energy Momentum Tensor for Lattice Gauge Theories,Annals Phys.197(1990) 119

    1990

  7. [7]

    Caracciolo, G

    S. Caracciolo, G. Curci, P. Menotti and A. Pelissetto,Renormalization of the Energy Momentum Tensor and the Trace Anomaly in Lattice QED,Phys. Lett. B228(1989) 375

    1989

  8. [8]

    Giusti and M

    L. Giusti and M. Pepe,Energy-momentum tensor on the lattice: Nonperturbative renormalization in Yang-Mills theory,Phys. Rev. D91(2015) 114504 [1503.07042]

    Pith/arXiv arXiv 2015

  9. [9]

    Dalla Brida, L

    M. Dalla Brida, L. Giusti and M. Pepe,Non-perturbative definition of the QCD energy-momentum tensor on the lattice,JHEP04(2020) 043 [2002.06897]

    arXiv 2020

  10. [10]

    Giusti and H.B

    L. Giusti and H.B. Meyer,Thermal momentum distribution from path integrals with shifted boundary conditions,Phys. Rev. Lett.106(2011) 131601 [1011.2727]

    Pith/arXiv arXiv 2011

  11. [11]

    Giusti and H.B

    L. Giusti and H.B. Meyer,Thermodynamic potentials from shifted boundary conditions: the scalar-field theory case,JHEP11(2011) 087 [1110.3136]

    Pith/arXiv arXiv 2011

  12. [12]

    Giusti and H.B

    L. Giusti and H.B. Meyer,Implications of Poincare symmetry for thermal field theories in finite-volume,JHEP01(2013) 140 [1211.6669]

    Pith/arXiv arXiv 2013

  13. [13]

    L¨ uscher,Properties and uses of the Wilson flow in lattice QCD,JHEP08 (2010) 071 [1006.4518]

    M. L¨ uscher,Properties and uses of the Wilson flow in lattice QCD,JHEP08 (2010) 071 [1006.4518]

    Pith/arXiv arXiv 2010

  14. [14]

    Suzuki,Energy–momentum tensor from the Yang–Mills gradient flow,PTEP 2013(2013) 083B03 [1304.0533]

    H. Suzuki,Energy–momentum tensor from the Yang–Mills gradient flow,PTEP 2013(2013) 083B03 [1304.0533]. 34

    Pith/arXiv arXiv 2013

  15. [15]

    Del Debbio, A

    L. Del Debbio, A. Patella and A. Rago,Space-time symmetries and the Yang-Mills gradient flow,JHEP1311(2013) 212 [1306.1173]

    Pith/arXiv arXiv 2013

  16. [16]

    Makino and H

    H. Makino and H. Suzuki,Lattice energy-momentum tensor from the Yang–Mills gradient flow—inclusion of fermion fields,PTEP2014(2014) 063B02 [1403.4772]

    Pith/arXiv arXiv 2014

  17. [17]

    Capponi, L

    F. Capponi, L. Del Debbio, A. Patella and A. Rago,Renormalization constants of the lattice energy momentum tensor using the gradient flow,PoS LATTICE2015(2016) 302 [1512.04374]. [18]WHOT-QCDcollaboration,ExploringN f = 2+1 QCD thermodynamics from the gradient flow,Phys. Rev. D96(2017) 014509 [1609.01417]

    Pith/arXiv arXiv 2016

  18. [18]

    Harlander, Y

    R.V. Harlander, Y. Kluth and F. Lange,The two-loop energy–momentum tensor within the gradient-flow formalism,Eur. Phys. J.C78(2018) 944 [1808.09837]

    arXiv 2018

  19. [19]

    Pavan, O

    P. Pavan, O. Kaczmarek, G.D. Moore and C. Schmidt,Shear viscosity from quenched to full lattice QCD,PoSLATTICE2024(2025) 199 [2503.11395]

    arXiv 2025

  20. [20]

    Roberge and N

    A. Roberge and N. Weiss,Gauge Theories With Imaginary Chemical Potential and the Phases of QCD,Nucl. Phys. B275(1986) 734

    1986

  21. [21]

    Dalla Brida, R

    M. Dalla Brida, R. H¨ ollwieser, F. Knechtli, T. Korzec, A. Ramos, S. Sint et al., High-precision calculation of the quark–gluon coupling from lattice QCD,Nature 652(2026) 328

    2026

  22. [22]

    Dalla Brida, L

    M. Dalla Brida, L. Giusti, T. Harris, D. Laudicina and M. Pepe,Non-perturbative thermal QCD at all temperatures: the case of mesonic screening masses,JHEP 04(2022) 034 [2112.05427]

    arXiv 2022

  23. [23]

    Bresciani, M

    M. Bresciani, M. Dalla Brida, L. Giusti and M. Pepe,QCD Equation of State with Nf=3 Flavors up to the Electroweak Scale,Phys. Rev. Lett.134(2025) 201904 [2501.11603]

    arXiv 2025

  24. [24]

    Bresciani, M

    M. Bresciani, M. Dalla Brida, L. Giusti and M. Pepe,QCD equation of state at very high temperature: Computational strategy, simulations, and data analysis, Phys. Rev. D113(2026) 034506 [2511.09160]. [26]ALPHAcollaboration,Monte Carlo errors with less errors,Comput. Phys. Commun.156(2004) 143 [hep-lat/0306017]

    arXiv 2026

  25. [25]

    Ramos,Automatic differentiation for error analysis of Monte Carlo data, Comput

    A. Ramos,Automatic differentiation for error analysis of Monte Carlo data, Comput. Phys. Commun.238(2019) 19 [1809.01289]

    Pith/arXiv arXiv 2019

  26. [26]

    Joswig, S

    F. Joswig, S. Kuberski, J.T. Kuhlmann and J. Neuendorf,pyerrors: A python framework for error analysis of Monte Carlo data,Comput. Phys. Commun.288 (2023) 108750 [2209.14371]. 35

    arXiv 2023

  27. [27]

    Laine and M

    M. Laine and M. Vepsalainen,On the smallest screening masses in hot QCD, JHEP09(2009) 023 [0906.4450]

    Pith/arXiv arXiv 2009

  28. [28]

    Giusti and M

    L. Giusti and M. L¨ uscher,Topological susceptibility atT > T c from master-field simulations of the SU(3) gauge theory,Eur. Phys. J. C79(2019) 207 [1812.02062]

    arXiv 2019

  29. [29]

    Bonati, M

    C. Bonati, M. D’Elia, M. Mariti, G. Martinelli, M. Mesiti, F. Negro et al.,Axion phenomenology andθ-dependence fromN f = 2 + 1lattice QCD,JHEP03(2016) 155 [1512.06746]

    Pith/arXiv arXiv 2016

  30. [30]

    Borsanyi et al.,Calculation of the axion mass based on high-temperature lattice quantum chromodynamics,Nature539(2016) 69 [1606.07494]

    S. Borsanyi et al.,Calculation of the axion mass based on high-temperature lattice quantum chromodynamics,Nature539(2016) 69 [1606.07494]

    Pith/arXiv arXiv 2016

  31. [31]

    Caracciolo, G

    S. Caracciolo, G. Curci, P. Menotti and A. Pelissetto,The Energy Momentum Tensor on the Lattice: The Scalar Case,Nucl. Phys.B309(1988) 612

    1988

  32. [32]

    Caracciolo, P

    S. Caracciolo, P. Menotti and A. Pelissetto,One loop analytic computation of the energy momentum tensor for lattice gauge theories,Nucl. Phys.B375(1992) 195

    1992

  33. [33]

    Caracciolo, P

    S. Caracciolo, P. Menotti and A. Pelissetto,Analytic determination at one loop of the energy momentum tensor for lattice QCD,Phys. Lett.B260(1991) 401

    1991

  34. [34]

    Bruno and R

    M. Bruno and R. Sommer,On fits to correlated and auto-correlated data,Comput. Phys. Commun.285(2023) 108643 [2209.14188]

    arXiv 2023

  35. [35]

    Giusti and M

    L. Giusti and M. Pepe,Equation of state of the SU(3) Yang–Mills theory: A precise determination from a moving frame,Phys. Lett.B769(2017) 385 [1612.00265]

    Pith/arXiv arXiv 2017

  36. [36]

    L¨ uscher,openQCD: Simulation programs for lattice QCD, https://luscher.web.cern.ch/luscher/openQCD/

    M. L¨ uscher,openQCD: Simulation programs for lattice QCD, https://luscher.web.cern.ch/luscher/openQCD/

  37. [37]

    Sheikholeslami and R

    B. Sheikholeslami and R. Wohlert,Improved Continuum Limit Lattice Action for QCD with Wilson Fermions,Nucl. Phys.B259(1985) 572

    1985

  38. [38]

    L¨ uscher, S

    M. L¨ uscher, S. Sint, R. Sommer and P. Weisz,Chiral symmetry and O(a) improvement in lattice QCD,Nucl. Phys.B478(1996) 365 [hep-lat/9605038]. [41]JLQCD, CP-PACScollaboration,Non-perturbative O(a)-improvement of Wilson quark action in three-flavor QCD with plaquette gauge action,Phys. Rev. D71(2005) 054505 [hep-lat/0406028]

    Pith/arXiv arXiv 1996

  39. [39]

    Rabinowitz,The exact degree of precision of generalized gauss-kronrod integration rules,Mathematics of Computation35(1980) 1275

    P. Rabinowitz,The exact degree of precision of generalized gauss-kronrod integration rules,Mathematics of Computation35(1980) 1275. 36

    1980

  40. [40]

    Burgio, S

    G. Burgio, S. Caracciolo and A. Pelissetto,Algebraic algorithm for the computation of one loop Feynman diagrams in lattice QCD with Wilson fermions, Nucl. Phys.B478(1996) 687 [hep-lat/9607010]

    Pith/arXiv arXiv 1996

  41. [41]

    Capitani and G

    S. Capitani and G. Rossi,Deep inelastic scattering in improved lattice QCD. 1. The First moment of structure functions,Nucl. Phys.B433(1995) 351 [hep-lat/9401014]

    Pith/arXiv arXiv 1995

  42. [42]

    Capitani,Lattice perturbation theory,Phys

    S. Capitani,Lattice perturbation theory,Phys. Rept.382(2003) 113 [hep-lat/0211036]

    Pith/arXiv arXiv 2003

  43. [43]

    Y.-B. Yang, M. Glatzmaier, K.-F. Liu and Y. Zhao,The 1-loop correction of the QCD energy momentum tensor with the overlap fermion and HYP smeared Iwasaki gluon,1612.02855

    Pith/arXiv arXiv

  44. [44]

    Cheng and N.J

    S.H. Cheng and N.J. Higham,A Modified Cholesky Algorithm Based on a Symmetric Indefinite Factorization,SIAM Journal on Matrix Analysis and Applications19(1998) 1097 [https://doi.org/10.1137/S0895479896302898]. 37

    work page doi:10.1137/s0895479896302898 1998