REVIEW 2 major objections 5 minor 27 references
$\mathrm{O}(a)$ improvement of the flavour singlet scalar density in a setup with Wilson fermions
T0 review · 2 major / 5 minor · reviewed 2026-08-07 · deepseek-v4-flash
Pith's one-line read The paper reports the first non-perturbative estimates of the O(a) improvement coefficient $g_S$ of the flavour singlet scalar density, finding values compatible with zero at the smallest couplings.
desk verdict First non-perturbative g_S estimates in the CLS range, clean Ward identity, but the b_g connection is still a plan and the omitted clover term leaves the O(a) ambiguity unquantified. 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 central object is the O(a)-improved flavour singlet scalar density $(S_I)^0$, with improvement pattern containing a power-divergent term proportional to $e_S$ and an O(a) term proportional to $g_S$ times the field-strength insertion $e\mathrm{Tr}[F_{\mu\nu}F_{\mu\nu}]$, discretised in the local, site-symmetric form $\{e\mathrm{Tr}[FF]\}_{S_g}$. The argument is carried by a Ward identity obtained by subtracting two axial Ward identities, one with the operator $S_0 O_{ext}$ and one with $1\cdot O_{ext}$; this subtraction cancels the divergent $e_S$ term and leaves $g_S$ as the only unknown once the renormalisation factors $Z$, $r_m$ and the improvement coefficient $c_A$ are supplied. The route to $b_g$ runs through $b_g = 2g_0^2 g_S$, which requires replacing $\{e\mathrm{Tr}[FF]\}_{S_g}$ by $\{e\mathrm{Tr}[FF]\}_{S_g} + a\,g_0^2\,(\partial c_{sw}/\partial g_0^2)\,(i/2)O_{clover}$ in the Ward identity.
What would settle it
Evaluate the extended Ward identity with the clover-term replacement of eq. (14) on the same five ensembles: if the resulting $g_S$ differs from the values in Fig. 1 by more than the quoted errors, the current estimates are only an intermediate quantity, and if the values remain zero-compatible after the clover term, the conclusion that $g_S$ is small in this coupling range is confirmed.
Extended reading notes
Core claim
On its own terms, the paper claims: non-perturbative estimates for $g_S$ have been determined for the first time. The extraction uses the chiral Ward identity (10), solved on five ensembles at $g_0^2 \in [1.5,1.77]$ with nearly massless, O(a)-improved Wilson fermions and the tree-level Symanzik improved gauge action. The estimates, shown in Fig. 1 for two choices of the time interval, approach zero for small couplings and are compatible with zero within 1-$\sigma$ at the two smallest couplings, confirming that $g_S$ is small in this range. The paper further establishes that the relation $b_g = 2g_0^2 g_S$ holds only after an additional, clover-term-like contribution is included in the Ward identity, and it outlines the implementation of that contribution as the next step toward non-perturbative $b_g$.
Load-bearing premise
The result rests on the assumption that the field-strength term used to improve the singlet scalar density, built only from the gauge action, is the complete O(a) improvement, even though the paper itself notes that a clover-term piece must be added before $g_S$ can be translated into $b_g$; if that missing piece moves the numbers by more than the quoted errors, the values reported here are not the final $b_g$-relevant ones.
Editorial extensions
If this is right
- Sigma-term determinations with Wilson fermions at these couplings can continue to set $g_S=0$ without introducing a bias larger than the present statistical uncertainty.
- The clover-term extension of the Ward identity gives a concrete route to non-perturbative $b_g$ values in the coupling range $g_0^2 \in [1.5,1.77]$, where only one-loop perturbation theory is currently available.
- Because the $g_S$ estimates are compatible with zero at the two finest lattices, the perturbative expectation of a small $g_S$ appears adequate near the continuum limit in this range.
- The correlation functions and disconnected-diagram estimators developed here can be reused for other flavour-singlet Ward identities on the same ensembles.
Reading between the lines
- If the missing clover-term contribution shifts $g_S$ by an amount comparable to or larger than the quoted errors, earlier sigma-term analyses that neglected $g_S$ would need a revised O(a) uncertainty estimate; the paper does not quantify that shift.
- The same Ward identity could be evaluated at additional lattice spacings to test whether the trend toward zero continues monotonically; a non-monotonic behaviour would signal that the remaining O(a) ambiguity is not negligible.
- The stochastic-estimator technology for the disconnected diagrams transfers to other flavour-singlet quantities, such as the topological-charge density or the singlet axial current, where analogous improvement coefficients appear.
- A non-perturbative $b_g$ obtained through $b_g = 2g_0^2 g_S$ would also test the one-loop assumption currently used in the scale and mass settings of large-volume Wilson-fermion simulations.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper reports a first non-perturbative determination of the O(a) improvement coefficient g_S of the flavour singlet scalar density in Nf=3 lattice QCD with Wilson-clover fermions and tree-level Symanzik-improved gauge action. The authors derive a lattice Ward identity (eq. (10)) from two subtracted chiral Ward identities, using Schrödinger functional boundary conditions and previously determined renormalisation/improvement parameters Z, r_m, and c_A. The identity is solved on five CLS ensembles with g0^2 in [1.5,1.77], yielding g_S estimates that approach zero at small coupling and are compatible with zero at the two smallest couplings. The paper also discusses the relation b_g = 2 g0^2 g_S, pointing out that an additional clover-term contribution (eq. (14)) is required for this relation and that the current results cannot yet be used to extract b_g.
Significance. If the results hold, they provide the first direct non-perturbative estimates of g_S in the coupling range used in large-volume CLS simulations, which is relevant for O(a) improvement of meson and baryon sigma terms. The derivation is careful and the numerical analysis uses standard tools (Gamma-method errors, trivial-topology projection, LCP ensembles with nearly massless quarks). The paper is transparent about the missing clover term and clearly states that the b_g extraction is future work. The main value is a proof-of-principle and a set of finite-lattice-spacing estimates that can guide sigma-term calculations, though continuum extrapolation and the clover-term correction remain to be completed.
major comments (2)
- [Sections 4 and 5, eq. (14)] The reported g_S is defined without the clover-term contribution that eq. (14) shows is necessary for the b_g relation. The paper itself calls this "merely an O(a) ambiguity" (Section 4). At the lattice spacings used (a ~ 0.06-0.1 fm), an O(a) ambiguity can be comparable to the quoted statistical errors and to the small extracted values near the chiral point, so the numbers in Fig. 1 may not be the g_S needed for sigma-term improvement if a different discretisation of the field strength is ultimately adopted. A quantitative estimate of the difference between the two definitions (e.g., a perturbative estimate or an evaluation of the clover correlation functions on one ensemble) should be provided, or the abstract should make clear that the presented g_S is an intermediate, definition-dependent quantity and that the b_g extraction is not yet realised.
- [Introduction and Section 5] The claim that this is the first non-perturbative determination of g_S needs qualification with respect to Ref. [7], which determines b_g non-perturbatively in the neighbouring coupling range g0^2 in [0.4,1.5]. Since b_g = 2 g0^2 g_S holds when the clover term of eq. (14) is included, Ref. [7] effectively provides a non-perturbative g_S in a different scheme/coupling range. The paper should explicitly state why the present determination is novel (direct Ward identity, different definition, different coupling range) or soften the "first time" wording to avoid an overstated novelty claim.
minor comments (5)
- [Abstract] The phrase "a relation to b_g ... can also be established, allowing for its non-perturbative extraction as well" is forward-looking; the abstract should make clear that the extraction of b_g is not performed in this paper and that the presented g_S cannot yet be used for that purpose.
- [Section 4, Fig. 1] The figure caption notes that the perturbative prediction is not directly comparable to the data, but the main text would benefit from repeating this caveat when the data are discussed, to prevent the reader from inferring a discrepancy.
- [Footnote 1 and Section 4] The sign convention g_S = -d_S with respect to Ref. [7] is only given in a footnote; it should also appear in the main text where g_S is first introduced, because the comparison with perturbative expressions and with future work depends on this sign.
- [Table 1] The column heading "Nsep [MDU]" is ambiguous; the text explains it, but the table would be clearer with a heading such as "Nsep [MDU]" followed by a footnote, or an explicit statement that Nsep is the separation in molecular dynamics units.
- [Section 4, eqs. (11)-(12)] The definitions of the connected correlation functions are deferred to Ref. [11]; since this is a proceedings paper, a brief restatement or at least a clear mapping of the notation (e.g., "con", "disc") to the definitions in the appendix of [11] would improve readability.
Circularity Check
No circularity: g_S is solved from a Ward identity whose inputs come from independent non-singlet determinations; the omitted clover term is a scheme limitation, not a circular reduction.
full rationale
The central claim — a first non-perturbative estimate of g_S — is obtained by solving the Ward identity (10) for g_S as the only unknown. The inputs Z, r_m and c_A are taken from refs. [11] and [12], which are published non-perturbative determinations based on flavour non-singlet Ward identities and a separate axial-current improvement condition; neither assumes the value of g_S, so these self-citations are independent support rather than circular inputs. The perturbative expression (3) is used only as an order-of-magnitude reference, and the paper explicitly states that the non-perturbative g_S results cannot be directly compared with it because of the missing clover-term contribution in eq. (14). That omission is a scheme/completeness limitation: the reported g_S is defined with the gauge-action-based discretisation {eTr[FF]}_{Sg}, and the paper transparently says that an additional clover-term-like discretisation is needed before b_g can be extracted via b_g = 2 g0^2 g_S. This is not a circular reduction — the Ward identity is a genuine constraint relating measured correlation functions to g_S, with no fitted parameter being relabelled as a prediction. The caveat in Section 5 that the present g_S cannot yet be used for b_g is an honest limitation, not a sign that the derivation reduces to its own inputs. Hence the circularity score is 0.
Assumptions & free parameters
assumptions (6)
- domain assumption Wilson-clover fermions and tree-level Symanzik improved gauge action describe QCD at these couplings.
- domain assumption Schrödinger functional boundary conditions with one-loop boundary improvement coefficients c_t are sufficient for O(a) improvement of the boundary.
- domain assumption Ward identities hold in a fixed topological sector, so projection to the trivial sector is valid.
- domain assumption O(am) terms can be neglected at the simulated quark masses.
- ad hoc to paper The improvement pattern of the singlet scalar density, eq. (8), with the specific discretisation {eTr[FF]}_{Sg}, is complete up to O(a^2).
- standard math The renormalisation parameters Z, r_m from [11] and c_A from [12] are used as external inputs.
Cite this review
Pith. "Pith review of $\mathrm{O}(a)$ improvement of the flavour singlet scalar density in a setup with Wilson fermions." pith.science (2026). https://pith.science/paper/VJYZUZJY
@misc{pith2026250208797,
author = {Pith},
title = {Pith review of: $\mathrmO(a)$ improvement of the flavour singlet scalar density in a setup with Wilson fermions},
year = {2026},
howpublished = {\url{https://pith.science/paper/VJYZUZJY}},
note = {Machine review of arXiv:2502.08797}
}
abstract
We report on our Ward identity determination of the $\mathrm{O}(a)$ improvement coefficient for the flavour singlet scalar density, namely $g_\mathrm{S}$, from three-flavour lattice QCD with Wilson-clover fermions and the tree-level Symanzik improved gauge action. We employ five couplings, $g_0^2 \in [1.5,1.77]$, that cover the range used in large-volume CLS simulations. While $g_\mathrm{S}$ itself is for instance relevant for the $\mathrm{O}(a)$ improvement of meson and baryon sigma terms, a relation to $b_\mathrm{g}$, the $\mathrm{O}(a)$ improvement parameter of the gauge coupling, can also be established, allowing for its non-perturbative extraction as well. With Wilson fermions, $b_\mathrm{g}$ is in principle required for full $\mathrm{O}(a)$ improvement at non-vanishing sea quark masses. We outline our procedure for extracting $b_\mathrm{g}$.
Figures
Reference graph
Works this paper leans on
-
[11]
J. Heitger, F. Joswig, P. L. J. Petrak and A. Vladikas,Ratio of flavour non-singlet and singlet scalar density renormalisation parameters in𝑁f = 3QCD with Wilson quarks, Eur. Phys. J. C 81(2021) 606 [2101.10969], [Erratum: Eur. Phys. J. C 82, 104 (2022)]
arXiv 2021
-
[7]
ALPHAcollaboration, M. Dalla Brida, R. Höllwieser, F. Knechtli, T. Korzec, S. Sint and R. Sommer,Heavy Wilson quarks and O(a) improvement: nonperturbative results for b𝑔, JHEP 01 (2024) 188 [2401.00216]
arXiv 2024
-
[1]
T. Bhattacharya, R. Gupta, W. Lee, S. R. Sharpe and J. M. Wu,Improved bilinears in lattice QCD with non-degenerate quarks, Phys. Rev. D73(2006) 034504 [hep-lat/0511014]
arXiv 2006
-
[2]
Aoki et al.,FLAG Review 2021, Eur
Flavour Lattice A veraging Group (FLAG)collaboration, Y. Aoki et al.,FLAG Review 2021, Eur. Phys. J. C82(2022) 869 [2111.09849]
arXiv 2022
-
[3]
Aoki et al.,FLAG Review 2024, 2411.04268
Flavour Lattice A veraging Group (FLAG)collaboration, Y. Aoki et al.,FLAG Review 2024, 2411.04268
arXiv 2024
-
[4]
A.Agadjanov, D.Djukanovic, G.vonHippel, H.B.Meyer, K.OttnadandH.Wittig, Nucleon Sigma Terms with𝑁f = 2+ 1 Flavors of O(a)-Improved Wilson Fermions, Phys. Rev. Lett. 131 (2023) 261902 [2303.08741]
arXiv 2023
-
[5]
M. Lüscher, S. Sint, R. Sommer and P. Weisz,Chiral symmetry andO(𝑎) improvement in lattice QCD, Nucl. Phys. B478 (1996) 365 [hep-lat/9605038]
arXiv 1996
-
[6]
Further results on O(a) improved lattice QCD to one-loop order of perturbation theory
S. Sint and P. Weisz,Further results on O(a) improved lattice QCD to one loop order of perturbation theory,Nucl. Phys. B502 (1997) 251 [hep-lat/9704001]
work page Pith review arXiv 1997
Show all 27 references
-
[8]
Dalla Brida, R
ALPHAcollaboration, M. Dalla Brida, R. Höllwieser, F. Knechtli, T. Korzec, A. Nada, A. Ramos, S. Sint and R. Sommer,Determination of𝛼𝑠(𝑚𝑍) by the non-perturbative decoupling method,Eur. Phys. J. C82 (2022) 1092 [2209.14204]
2022 arXiv
-
[9]
RQCDcollaboration, G. S. Bali, S. Collins, P. Georg, D. Jenkins, P. Korcyl, A. Schäfer, E. E. Scholz, J. Simeth, W. Söldner and S. Weishäupl,Scale setting and the light baryon spectrum in𝑁f = 2+ 1QCD with Wilson fermions,JHEP 05 (2023) 035 [2211.03744]
2023 arXiv
-
[10]
Sint,Renormalization of the iso-singlet scalar density in lattice QCD with Wilson quarks, Nucl
S. Sint,Renormalization of the iso-singlet scalar density in lattice QCD with Wilson quarks, Nucl. Phys. B Proc. Suppl.73(1999) 282
1999
-
[12]
Bulava, M
J. Bulava, M. Della Morte, J. Heitger and C. Wittemeier,Non-perturbative improvement of the axial current in𝑁f = 3 lattice QCD with Wilson fermions and tree-level improved gauge action,Nucl. Phys. B896(2015) 555 [1502.04999]
2015 arXiv
-
[13]
Fritzsch, J
P. Fritzsch, J. Heitger and S. Kuberski,O(𝑎) improved quark mass renormalization for a non-perturbative matching of HQET to three-flavor QCD, PoS LATTICE2018(2018) 218 [1811.02591]. 9 O(𝑎) improvement of the flavour singlet scalar density Pia Leonie Jones Petrak
2018 arXiv
-
[14]
Conigli, J
A. Conigli, J. Frison, P. Fritzsch, A. Gérardin, J. Heitger, G. Herdoiza, S. Kuberski, C. Pena, H. Simma and R. Sommer,𝑚𝑏 and 𝑓𝐵(∗) in 2+ 1 flavour QCD from a combination of continuum limit static and relativistic results, PoS LATTICE2023(2024) 237 [2312.10017]
2024 arXiv
-
[15]
Bulava and S
J. Bulava and S. Schaefer,Improvement of𝑁f = 3 lattice QCD with Wilson fermions and tree-level improved gauge action,Nucl. Phys. B874 (2013) 188 [1304.7093]
2013 arXiv
-
[16]
Dalla Brida, P
ALPHAcollaboration, M. Dalla Brida, P. Fritzsch, T. Korzec, A. Ramos, S. Sint and R. Sommer,Slow running of the Gradient Flow coupling from 200 MeV to 4 GeV in𝑁f = 3 QCD, Phys. Rev. D95(2017) 014507 [1607.06423]
2017 arXiv
-
[17]
Giusti, T
L. Giusti, T. Harris, A. Nada and S. Schaefer,Frequency-splitting estimators of single-propagator traces, Eur. Phys. J. C79(2019) 586 [1903.10447]
2019 arXiv
-
[18]
Lüscher and S
M. Lüscher and S. Schaefer.http://luscher.web.cern.ch/luscher/openQCD
-
[19]
Fritzsch, A
P. Fritzsch, A. Ramos and F. Stollenwerk,Critical slowing down and the gradient flow coupling in the Schrödinger functional,PoS Lattice2013 (2014) 461 [1311.7304]
2014 arXiv
-
[20]
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]
2023 arXiv
-
[21]
Wolff,Monte Carlo errors with less errors,Comput
U. Wolff,Monte Carlo errors with less errors,Comput. Phys. Commun.156 (2004) 143 [hep-lat/0306017], [Erratum: Comput. Phys. Commun. 176, 383 (2007)]
2004 arXiv
-
[22]
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]
2019 arXiv
-
[23]
Schaefer, R
ALPHAcollaboration, S. Schaefer, R. Sommer and F. Virotta,Critical slowing down and error analysis in lattice QCD simulations,Nucl. Phys. B845 (2011) 93 [1009.5228]
2011 arXiv
-
[24]
Bulava, M
J. Bulava, M. Della Morte, J. Heitger and C. Wittemeier,Nonperturbative renormalization of the axial current in𝑁f = 3 lattice QCD with Wilson fermions and a tree-level improved gauge action,Phys. Rev. D93 (2016) 114513 [1604.05827]
2016 arXiv
-
[25]
Heitger, F
J. Heitger, F. Joswig and A. Vladikas,Ward identity determination of𝑍S/𝑍P for𝑁f = 3 lattice QCD in a Schrödinger functional setup,Eur. Phys. J. C80(2020) 765 [2005.01352]
2020 arXiv
-
[26]
Heitger and F
ALPHAcollaboration, J. Heitger and F. Joswig,The renormalisedO(𝑎) improved vector current in three-flavour lattice QCD with Wilson quarks, Eur. Phys. J. C81 (2021) 254 [2010.09539]
2021 arXiv
-
[27]
Chimirri, P
L. Chimirri, P. Fritzsch, J. Heitger, F. Joswig, M. Panero, C. Pena and D. Preti, Non-perturbative renormalisation and improvement of non-singlet tensor currents in Nf = 3 QCD, JHEP 07(2024) 089 [2309.04314]. 10
2024 arXiv
Reviewed August 7, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.