REVIEW
Improved bilinears in lattice QCD with non-degenerate quarks
High Energy Physics - Lattice · Submitted 4 Nov 2005 · 28 references
Pith's verdict: A careful, honest bookkeeping paper that maps which Ward identities fix which O(a) improvement coefficients for non-degenerate Wilson fermions — the algebra holds up, and the self-correction in Appendix C is the most valuable part.
T0 · AI review
Why Pith says this
Highlight any text to comment on it.
Referee report
The paper
The paper claims
The paper extends Symanzik O(a) improvement of Wilson-fermion bilinears from degenerate to non-degenerate dynamical quarks, and argues that for three or four flavors the resulting proliferation of improvement coefficients is almost entirely
The reading
What carries the argument
The engine is a spurion enumeration combined with Ward identity bookkeeping. The mass matrix $M$ is treated as transforming in the adjoint of the flavor group, which restricts the allowed $O(a)$ terms to the anticommutator structure $tr({{lambda, M} O})$ and forbids the commutator by CP; this yields the master formulas eqs. (15), (23), (24) and (25) for non-singlet bilinears, singlet bilinears, non-singlet masses and the singlet mass. The second ingredient is the relation $partial m_H / partial hat{m}_j = V langle H | hat{S}^{(jj)} | H rangle / langle H | H rangle$, enforced not only in the continuum limit but at finite lattice spacing; because the fermion action depends on the bare masses only through $a m_j S^{(jj)}$, the derivative of the hadron mass with respect to a bare mass is a known scalar matrix element, and the gluon term supplies $b_g$. Chaining Jacobians from improved to bare masses and matching coefficients produces Table I. The third ingredient is the classification of contact terms: when the axial variation of the action meets another operator at the same point, extra operators of the appropriate dimension and quantum numbers appear, and the paper's rule is to enumerate them by symmetry rather than by the equations-of-motion mnemonic used earlier, which the authors state cannot be trusted because contact terms do not factorize.
Extended reading notes
Core claim
The central claim is a counting and determination result. Away from the chiral limit, each non-singlet bilinear needs three improvement coefficients beyond its chiral-limit normalization, written $b_O$, $b_O$-bar and $f_O$ in the general ansatz of eq. (15); each singlet bilinear needs two, $d_O$ and $d_O$-bar, plus a ratio $r_O$ of singlet to non-singlet normalization; and the quark masses need their own set, $r_m$, $b_m$, $b_m$-bar, $d_m$, $d_m$-bar. The authors show that for $N_f = 3$ the vector charge normalization fixes all vector coefficients except $c_V$ and $c_V$-bar, that matching the mass derivative of hadron masses to the scalar matrix element fixes the relations collected in Table I (including $b_S = -2b_m$, $b_S$-bar $= -b_m$-bar, and $b_g = 2g_0^2 g_S$), and that two-point and three-point axial Ward identities then determine essentially everything else. Only two scale-independent combinations survive: $d_A$-bar and $b_T$-bar $+ d_T$-bar; the scale-dependent $Z_T$, $Z_S Z_P$ and $r_A$ are also out of reach, as they must be, since Ward identities carry no renormalization scale. For $N_f = 2$ the same machinery leaves eight combinations undetermined, including $b_A$-bar, so that even the normalization of the non-singlet axial current away from the chiral limit is not fixed. The paper is explicitly a correction of an earlier version of the same analysis: enforcing vector Ward identities was previously claimed to determine more than it does, and Appendix C traces the error to a contact term proportional to $tr(O)$ that had been missed, which is exactly what allows the $f_O$ term to exist without being forced to vanish.
Load-bearing premise
The whole determination rests on being able to list every $O(a)$ contact term that can appear when the axial variation of the action meets another operator, and the authors show their own earlier list was incomplete, so the corrected list, obtained by counting symmetries rather than by deriving the terms, is the least controlled step.
Editorial extensions
If this is right
- For a 2+1 flavor simulation, the coefficients needed for the flavor off-diagonal axial current ($Z_A$, $b_A$, $b_A$-bar) and for the electromagnetic current ($Z_V$, $b_V$, $b_V$-bar, $f_V$) are all obtainable from Ward identities, with the vector ones requiring quark-disconnected contractions.
- Because $b_g = 2g_0^2 g_S$, determining $g_S$ from the axial transformation of the singlet scalar gives an independent, self-contained determination of $b_g$, so the improvement programme does not have to assume $b_g$ is known in advance.
- The improved quark condensate can be obtained by evaluating the axial variation of the singlet pseudoscalar at $tilde{m}_{jk} = 0$, which removes the power-divergent identity-operator contribution by explicit subtraction.
- Two scale-independent combinations, $d_A$-bar and $b_T$-bar $+ d_T$-bar, plus the scale-dependent $Z_T$, $Z_S Z_P$ and $r_A$, must come from non-perturbative renormalization or from matching short-distance correlators to perturbation theory.
- For two flavors, eight combinations of scale-independent coefficients remain unknown, so the normalization of the non-singlet axial, pseudoscalar and tensor bilinears away from the chiral limit is not determined by Ward identities alone.
Reading between the lines
- The determination scheme's real vulnerability is not the algebra but the completeness of the contact-term list: a numerical 2+1 flavor check that the tabulated coefficient set makes three-point axial Ward identity ratios independent of source and separation to $O(a^2)$ would validate the whole chain, and a residual $O(a)$ drift would locate exactly where the enumeration is short.
- The contrast between three and two flavors suggests the operative resource is the number of independently variable mass handles rather than the number of flavors, which makes a graded or partially quenched extension the natural next test, though the paper itself cautions that a non-unitary sea weakens the Symanzik basis.
- If $d_A$-bar and $b_T$-bar $+ d_T$-bar are genuinely invisible to Ward identities, a two-loop perturbative estimate of both would settle whether the gap matters numerically at $a^{-1}$ between 2 and 4 GeV, in the same spirit as the paper's own tree-level estimate for $b_A$.
- The paper's remark that Wilson averaging and maximally twisted mass make these coefficients unnecessary at the cost of requiring an even number of flavors frames the result as a benchmark: the improvement coefficients computed here are precisely the size of the effect those approaches are designed to cancel.
Original abstract
We describe the extension of the improvement program for bilinear operators composed of Wilson fermions to non-degenerate dynamical quarks. We consider two, three and four flavors, and both flavor non-singlet and singlet operators. We find that there are many more improvement coefficients than with degenerate quarks, but that, for three or four flavors, nearly all can be determined by enforcing vector and axial Ward identities. The situation is worse for two flavors, where many more coefficients remain undetermined.
The compiled record
The paper's own sentences, arranged, kept beside the review and separate from it. Each label is drawn from a closed vocabulary; the number beside it is the probability the classifier gave that label. Numbers in the tables are the authors', checked by code. No sentence in this record was written by a model.
The abstract, sentence by sentence
We describe the extension of the improvement program for bilinear operators composed of Wilson fermions to non-degenerate dynamical quarks. We consider two, three and four flavors, and both flavor non-singlet and singlet operators. We find that there are many more improvement coefficients than with degenerate quarks, but that, for three or four flavors, nearly all can be determined by enforcing vector and axial Ward identities. The situation is worse for two flavors, where many more coefficients remain undetermined.
background contribution method result interpretation main claim
| # | sentence | role | the authors' own work | states the main conclusion | strength |
|---|---|---|---|---|---|
| s0 | We describe the extension of the improvement program for bilinear operators composed of Wilson fermions to non-degenerate dynamical quarks. | contribution 1.00 | 0.89 | 0.09 | asserted 0.94 |
| s1 | We consider two, three and four flavors, and both flavor non-singlet and singlet operators. | contribution 0.67 | 0.49 | 0.04 | asserted 0.87 |
| s2 | We find that there are many more improvement coefficients than with degenerate quarks, but that, for three or four flavors, nearly all can be determined by enforcing vector and axial Ward identities. | result 1.00 | 0.95 | 0.67 | asserted 1.00 |
| s3 | The situation is worse for two flavors, where many more coefficients remain undetermined. | result 0.63 | 0.79 | 0.10 | asserted 1.00 |
What was examined, and what was not
- arXiv listing abstract v3
- full text
- figures and tables
- references and citing sentences
- relations to other papers
- P abstract
- S sections
- T body sentences
- C citations
- R cross-paper relations
References
How later papers used it
60 later papers cite it
3 later papers by other authors are listed below: 1 built on it and 2 used its method or tool. 1 more paper by the same authors is listed: 1 used its method or tool.
Builds on it 1 paper
-
Combined the argument for (a) improvement for TM with an extension of this paper to Wilson-like theories, extending work originally formulated only for Wilson fermions.
Heavy-quark physics with a tmQCD valence action (2018) · 0.98 · paraphrased: the licence does not allow quoting
Uses its method or tool 3 of 15 papers with this label · 1 by the same authors
-
(23) of [ref] , we obtain [equation] Again, we replaced the lattice decay constants within the $𝒪(a)$ terms by partially improved ones.
Masses and decay constants of the $\eta$ and $\eta^\prime$ mesons from lattice QCD (2021) · Renormalization and $𝒪(a)$ improvement · 0.99
-
The renormalised $O(a)$ improved current reads [ref] [equation] where $m_qc$ and $M$ denote the bare vector Ward identity quark mass combinations [equation] The hopping parameter for quark flavour $q$ is denoted by $κ_q$ and $κ_crit$ labels its critical value.
$\mathrm{D}$ and $\mathrm{D_s}$ decay constants in $N_{\rm f}=2+1$ QCD with Wilson fermions (2024) · Observables · 0.84
By the same authors
-
We will do this using the notation and results given in Ref. [ref] for an $N_f$ flavor theory with two light degenerate flavors, $m_u = m_d ≡m_l$, and $N_f-2$ heavier flavors denoted generically by $m_s$.
The pion-nucleon sigma term from lattice QCD (2021) (same authors) · Renormalization · 0.90
How the literature uses it
58 arXiv papers cite this in 118 sentences
Papers citing it · 89 sentences without a confident label
Contests it 3 papers
-
We would therefore like to point out a subtlety that was overlooked in [ref] .
Heavy Wilson Quarks and O($a$) Improvement: Nonperturbative Results for $b_{\rm g}$ · Chiral Ward identities and $$ · contest 0.80
-
For a detailed derivation, which also corrects an oversight in [ref] , we refer to our paper [ref] .
A non-perturbative determination of $b_{\rm g}$ · Improvement condition for $$ · contest 0.79
-
This introduces an additional term $∝(g_0^2)$ in the derivation which was possibly overlooked in the original publication [ref] .
$\mathrm{O}(a)$ improvement of the flavour singlet scalar density in a setup with Wilson fermions · Conclusion and Outlook · contest 0.76
Uses its method or tool 13 papers
-
For the currents relevant for this study, we perform the substitutions [ref] ($a=1,…,8$): [equation] where we re-expressed the last equation using the singlet AWI in the massless case.[Note that $c_P^s=-16π^2 g_P/√(6)$ to leading order in $g^2$.
Masses and decay constants of the $\eta$ and $\eta^\prime$ mesons from lattice QCD · Gauge ensembles · 0.93
-
This ratio can be determined non-perturbatively from the slope of the axial Ward identity quark mass ($m̃_q$) as a function of the vector Ward identity mass ($m_q$), see, e.g, Ref. [ref] : [equation] where the quark masses are defined as [equation] $A_4^I$ and $P$ denote the $𝒪(a)$ improved axial-vector current and the pseudoscalar operator, respectively, with corresponding renormalization factors $Z_A^ns$ and $Z_P^ns$. $κ_c,sea$ denotes the critical mass parameter along the isosymmetric line for which the quark mass is zero.
Direct determinations of the nucleon and pion $\sigma$ terms at nearly physical quark masses · Renormalization · 0.93
-
The renormalised $O(a)$ improved current reads [ref] [equation] where $m_qc$ and $M$ denote the bare vector Ward identity quark mass combinations [equation] The hopping parameter for quark flavour $q$ is denoted by $κ_q$ and $κ_crit$ labels its critical value.
$\mathrm{D}$ and $\mathrm{D_s}$ decay constants in $N_{\rm f}=2+1$ QCD with Wilson fermions · Observables · 0.84
Compares with it 1 paper
-
Finally we compare these results with those obtained in [ref] , to see whether we can match the $4$ improvement terms found in eq.
Flavour blindness and patterns of flavour symmetry breaking in lattice simulations of up, down and strange quarks · $𝐎(𝐚)$ improvement of quark masses · 1.00
Cites it as background 6 papers
-
We follow the standard notation in the literature [ref] where $c_*$ terms cancel $(a)$ cutoff effects and $b_*$, $b̅_*$ the leading mass-dependent effects of order $am$, induced from the valence or sea quark sector, respectively.
Mass-improvement of the vector current in three-flavor QCD · Prelude · 0.85
-
The $𝒪(a)$ improvement of quark bilinears has been worked out for Wilson fermions in [ref] .
Masses and decay constants of the $\eta$ and $\eta^\prime$ mesons from lattice QCD · Renormalization and $𝒪(a)$ improvement · 0.80
-
For those cases, the mass-dimension 4 operator basis relevant for $(a)$ corrections has been discussed earlier [ref] with explicit improvement in mind and not its impact on scaling due to the powers $Γ̂_i^$.
Lattice artifacts of local fermion bilinears up to $\mathrm{O}(a^2)$ · Introduction · 0.83
Discussion
Disagree with the review, or know something it missed? Leave a comment, highlight any passage above to comment on it, or challenge this review formally.