Pith. sign in

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 →

arxiv 2502.08797 v1 pith:VJYZUZJY submitted 2025-02-12 hep-lat

classification hep-lat PACS 11.15.Ha12.38.Gc
keywords O(a)improvementflavoursingletscalardensityg_SWardidentityWilsonfermionslatticeQCDb_gsigmaterms
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 aims to determine, for the first time, the non-perturbative value of $g_S$, the O(a) improvement coefficient of the flavour singlet scalar density in lattice QCD with Wilson-clover fermions. It derives a lattice Ward identity in which $g_S$ is the only unknown and evaluates it on five ensembles with $N_f=3$ mass-degenerate quarks at couplings $g_0^2 \in [1.5,1.77]$, the range used in large-volume simulations. The resulting estimates decrease toward zero at small couplings and are compatible with zero within one standard deviation at the two smallest couplings. This supports the common practice of neglecting $g_S$ in meson and baryon $\sigma$-term calculations, and the paper also shows how the same identity must be extended with a clover-term insertion before it can deliver the gauge-coupling improvement coefficient $b_g$.

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.

Watch

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

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

  • 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.
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

2 major / 5 minor

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)
  1. [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.
  2. [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)
  1. [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.
  2. [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.
  3. [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.
  4. [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.
  5. [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

0 steps flagged · score 0.0 of 10

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 0 free parameters · 6 assumptions · 0 invented entities

No new particles or entities are introduced. The only quantity determined is the standard improvement coefficient g_S. The key input that carries risk is the definition of the field-strength discretisation in the improvement pattern, not an invented entity.

assumptions (6)
  • domain assumption Wilson-clover fermions and tree-level Symanzik improved gauge action describe QCD at these couplings.
    This is the discretisation used throughout; it is the physical model of the theory being studied.
  • domain assumption Schrödinger functional boundary conditions with one-loop boundary improvement coefficients c_t are sufficient for O(a) improvement of the boundary.
    Stated in Section 3; relies on the standard ALPHA framework for mass-independent renormalisation and improvement.
  • domain assumption Ward identities hold in a fixed topological sector, so projection to the trivial sector is valid.
    Invoked in Section 3 with reference [19]; this is a standard practice to avoid critical slowing down of topological charge.
  • domain assumption O(am) terms can be neglected at the simulated quark masses.
    PCAC masses in Table 1 are small (e.g. 0.000011(14) at the finest lattice), so the O(am) term in eq. (9f) is dropped and the Ward identity is valid up to O(a^2).
  • 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).
    This is the load-bearing assumption identified in the review; the paper itself shows that an additional clover term is needed for the b_g relation, eq. (14).
  • standard math The renormalisation parameters Z, r_m from [11] and c_A from [12] are used as external inputs.
    These are prior non-perturbative determinations from the same collaboration, not derived in this paper.

how reviews work

0 comments
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

Figures reproduced from arXiv: 2502.08797 by the authors.

Figure 1
Figure 1. 𝑔S results obtained via eq. (10) for 𝑁f = 3 using our 𝑍 and 𝑟m interpolation formulas from [11] to construct the 𝑍𝑟m values at the given couplings needed. The perturbative prediction from eq. (3) (dashed red line) cannot be directly compared to the 𝑔S data points, see text below. The time interval [𝑡1, 𝑡2] (in eq. (10)) is set to [𝑇/4, 3𝑇/4] (labelled by T/4) or [𝑇/3, 2𝑇/3] (labelled by T/3) as shown by different co… view at source ↗

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

27 extracted references · 5 canonical work pages

  1. [11]

    Heitger, F

    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)]

  2. [7]

    Dalla Brida, R

    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]

  3. [1]

    Bhattacharya, R

    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]

  4. [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]

  5. [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

  6. [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]

  7. [5]

    Lüscher, S

    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]

  8. [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]

Show all 27 references
  1. [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]

  2. [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]

  3. [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

  4. [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]

  5. [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

  6. [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]

  7. [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]

  8. [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]

  9. [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]

  10. [18]

    Lüscher and S

    M. Lüscher and S. Schaefer.http://luscher.web.cern.ch/luscher/openQCD

  11. [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]

  12. [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]

  13. [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)]

  14. [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]

  15. [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]

  16. [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]

  17. [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]

  18. [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]

  19. [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

Pith tools

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