Pith. Sign in

REVIEW

Dynamical Twisted Mass Fermions with Light Quarks

PaperPDF from arXiv

data files

ReviewThe objections and verdict, as JSON Paper breakdownThe authors' claims, conditions and measurements, as JSON Signed receiptProof this review has not been altered
X Bluesky LinkedIn Reddit Hacker News
Ph. Boucaud, P. Dimopoulos, F. Farchioni, R. Frezzotti, V. Gimenez, G. Herdoiza, K. Jansen, V. Lubicz
show 14 more authors
G. Martinelli, C. McNeile, C. Michael, I. Montvay, D. Palao, M. Papinutto, J. Pickavance, G.C. Rossi, L. Scorzato, A. Shindler, S. Simula, C. Urbach, U. Wenger, ETM collaboration

High Energy Physics - Lattice · Submitted 15 Jan 2007 · 78 references

Pith's verdict: A genuine milestone dataset for twisted-mass QCD, but the headline LECs and lattice scale come from a chiral extrapolation that runs a factor of five in mPS² below the simulated region, anchored by the single most contaminated point.

T0 · AI review

How this review was made: An AI referee (selfhosted/deepseek-v4.1-flash) read the paper on November 11, 2026 and answered a fixed list of questions about it. T0 is the starting level every paper gets. The levels describe how deeply Pith has checked the paper, never the authors.

  1. T0 AI review An AI read the paper and answered a fixed list of questions.
  2. T1 Audit Automatic checks that the references and files the paper relies on exist and match.
  3. T2 Compute The paper's code or calculations re-run with the same results.
  4. T3 Formal statement The claims written out in a formal math language a computer can read.
  5. T4 Proof Those claims proved, and the proof checked by computer.

How Pith checks papers →

· 5 major objections · 7 minor · Ask about this paper · The compiled record · Read the full review →
Challenge this review Comment (0) Read on arXiv ↗
Cite ▾
arXiv IDhep-lat/0701012 (version 2)
Pith IDLLTZ432N, Pith's permanent record of this paper and its review
CitationA. Shindler, C. McNeile, C. Michael, C. Urbach, D. Palao, ETM Collaboration, F. Farchioni, G.C. Rossi, G. Herdoiza, G. Martinelli, I. Montvay, J. Pickavance, K. Jansen, L. Scorzato, M. Papinutto, P. Dimopoulos, Ph. Boucaud, R. Frezzotti, S. Simula, U. Wenger, V. Gimenez, V. Lubicz. "Dynamical Twisted Mass Fermions with Light Quarks." arXiv:hep-lat/0701012 (2007). https://pith.science/paper/LLTZ432N
Cite this reviewBibTeX for Pith's review of this paper
Claim paper

Why Pith says this

Highlight any text to comment on it.

Referee report

How Pith checks papers → · Full audit of this paper →

The paper

The paper claims

Two-flavour twisted-mass lattice QCD reaches pion masses near 300 MeV and extracts the chiral low-energy constants F, l3 and l4 with small statistical error.

The reading

The paper asks whether lattice QCD with two flavours of Wilson twisted mass quarks, tuned to maximal twist, can be pushed down to pseudo scalar masses near 300 MeV and still be described by next-to-leading-order chiral perturbation theory. If it can, then the chiral extrapolation every lattice calculation needs in order to reach physical quark masses is under control in this discretisation, and the low-energy constants of the chiral Lagrangian can be read off with small statistical error. The authors simulate at a single lattice spacing ($a$ just under 0.1 fm, $a^{-1}$ set by $r_0/a = 5.22(2)$) on a $24^3$ volume at five twisted masses, tune the PCAC quark mass to zero at the lightest one to define maximal twist, and fit the Gasser-Leutwyler formulae for $m_{PS}^2$ and $f_{PS}$ - including a continuum finite-size correction - to the four lightest points. They report $F = 121.3(7)$ MeV, $l_3 = 3.65(12)$ and $l_4 = 4.52(6)$, a lattice spacing of $0.087(1)$ fm fixed by $f_pi$, and a neutral pseudo scalar meson about 20% lighter than the charged one.

What carries the argument

The load-bearing device is maximal twist defined operationally: the bare untwisted mass is tuned so that the PCAC quark mass $m_{PCAC} = sum_x <d_0 A_0^a(x) P^a(0)> / (2 sum_x <P^a(x) P^a(0)>)$ vanishes at the smallest twisted mass, $a mu = 0.004$. Because the operator in the numerator becomes parity-odd in the physical basis, this single tuning condition delivers automatic $O(a)$ improvement without any operator-specific improvement coefficients, and the residual $mu$-dependence of $m_{PCAC}$ seen in figure 1(a) is argued to be an $O(a)$ effect that only reshuffles $O(a^2)$ artifacts. A second device is the exact lattice Ward identity $f_{PS} = 2 mu / m_{PS}^2 |<0|P^1(0)|pi>|$, which needs no renormalisation constant because $Z_P = 1/Z_mu$; this is what makes $f_{PS}$ a clean input to the chiral fit. The third is the pair of finite-volume chiral formulae, with the Gasser-Leutwyler correction $g_1(lambda)$ at $lambda = sqrt(2 B_0 mu) L^2$, fitted jointly to $m_{PS}$ and $f_{PS}$ so that their statistical correlation is carried through a bootstrap.

Extended reading notes

Core claim

The central claim is that at maximal twist the twisted-mass formulation delivers a set of pseudo scalar masses and decay constants, spanning roughly 300 to 500 MeV, that next-to-leading-order $N_f = 2$ chiral perturbation theory describes well - and that this description yields the low-energy constants. Concretely, fitting Eqs. (5) and (6) simultaneously to the four lightest twisted masses gives $2aB_0 = -4.99(6)$, $aF = -0.0534(6)$, $log(a^2 Lambda_3^2) = -1.93(10)$ and $log(a^2 Lambda_4^2) = -1.06(4)$, from which the authors extract $a mu_pi = 0.00078(2)$, $l_3 = 3.65(12)$, $l_4 = 4.52(6)$ and $F = 121.3(7)$ MeV. The same fit, anchored to the physical $f_pi$, gives $a = 0.087(1)$ fm, about 10% below the value obtained by setting $r_0 = 0.5$ fm. Adding the fifth and heaviest mass still describes $m_{PS}^2$ but misses $f_{PS}$, which the authors read as the expected edge of the chiral regime rather than a failure of the tuning.

Load-bearing premise

The load-bearing premise is that continuum next-to-leading-order chiral perturbation theory, together with its continuum finite-size correction, adequately describes data taken at a single lattice spacing and a single volume up to pseudo scalar masses around 500 MeV - an assumption the authors themselves flag in section 3.3(iii), where they note that the 20% neutral-charged pion splitting may make the continuum treatment of finite-size effects inadequate.

Editorial extensions

If this is right

  • The chiral extrapolation from $m_{PS} approx 300$ MeV to the physical point is under control in this discretisation, so the same tuning and fitting machinery can be carried to other lattice spacings and volumes.
  • Automatic $O(a)$ improvement at maximal twist is confirmed in the dynamical two-flavour case, extending what quenched scaling tests had shown.
  • The low-energy constants $F$, $l_3$ and $l_4$ become lattice-determined quantities with competitive statistical errors, providing an independent cross-check of values extracted from experimental data.
  • Isospin breaking from the twisted-mass term is a real but modest $O(a^2)$ effect, about 20% in the neutral-charged pion splitting at $a mu = 0.004$ and smaller than the quenched result with opposite sign, so it can be treated as a systematic rather than a blocker.
  • The tree-level Symanzik improved gauge action, chosen because it weakens the first-order phase transition that forbids light pions at coarser spacings, is a workable compromise for reaching small quark masses.

Reading between the lines

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

  • Because the whole analysis rests on one lattice spacing and one volume, the reported constants are best read as an existence proof of the method rather than as final numbers; the announced $32^3$ run and the second lattice spacing are what will actually decide whether the quoted errors are honest.
  • If the neutral pion really is 20% lighter than the charged one at these masses, then the finite-volume correction applied to $f_{PS}$ is being computed with the wrong pion mass in the loop; a mixed-action or partially quenched chiral treatment that allows unequal charged and neutral masses would be the natural repair.
  • The 10% discrepancy between the $f_pi$-based scale and the $r_0 = 0.5$ fm scale is a useful handle: it says the two common calibration routes disagree at roughly the size of the residual cutoff effects the paper hopes are small, and a second lattice spacing would separate them.
  • The same tuning-by-parity-odd-operator recipe should transfer unchanged to $N_f = 2+1+1$ twisted-mass simulations, where the lightest pseudo scalar mass is the binding constraint on cost.
Original abstract

We present results of dynamical simulations with 2 flavours of degenerate Wilson twisted mass quarks at maximal twist in the range of pseudo scalar masses from 300 to 550 MeV. The simulations are performed at one value of the lattice spacing a \lesssim 0.1 fm. In order to have O(a) improvement and aiming at small residual cutoff effects, the theory is tuned to maximal twist by requiring the vanishing of the untwisted quark mass. Precise results for the pseudo scalar decay constant and the pseudo scalar mass are confronted with chiral perturbation theory predictions and the low energy constants F, \bar{l}_3 and \bar{l}_4 are evaluated with small statistical errors.

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.

Compiled by classifier-1.13-20260917 via OpenRouter, question bank abstract-p1-20260923, 2026-09-23: 1 request, 3,690 input tokens, $0.0002 reported.

What the abstract itself states, as probabilities: claims first 0.03 · states limitation 0.47 · human participants 0.01 · controlled comparison 0.07 · randomized assignment 0.01 · states quantitative result 0.83 · compares with prior results 0.49 · mentions code or data availability 0.02.

The abstract, sentence by sentence

arXiv listing abstract v2. Results and the main claim are underlined; .

We present results of dynamical simulations with 2 flavours of degenerate Wilson twisted mass quarks at maximal twist in the range of pseudo scalar masses from 300 to 550 MeV. The simulations are performed at one value of the lattice spacing a \lesssim 0.1 fm. In order to have O(a) improvement and aiming at small residual cutoff effects, the theory is tuned to maximal twist by requiring the vanishing of the untwisted quark mass. Precise results for the pseudo scalar decay constant and the pseudo scalar mass are confronted with chiral perturbation theory predictions and the low energy constants F, \bar{l}_3 and \bar{l}_4 are evaluated with small statistical errors.

background contribution method result interpretation main claim

What the classifier reports for each sentence. The number is the probability it gave the label.
#sentencerolethe authors' own workstates the main conclusionstrength
s0 We present results of dynamical simulations with 2 flavours of degenerate Wilson twisted mass quarks at maximal twist in the range of pseudo scalar masses from 300 to 550 MeV. contribution 0.93 0.96 0.09 asserted 1.00
s1 The simulations are performed at one value of the lattice spacing a \lesssim 0.1 fm. method 0.99 0.94 0.03 asserted 0.98
s2 In order to have O(a) improvement and aiming at small residual cutoff effects, the theory is tuned to maximal twist by requiring the vanishing of the untwisted quark mass. method 1.00 0.73 0.05 asserted 0.98
s3 Precise results for the pseudo scalar decay constant and the pseudo scalar mass are confronted with chiral perturbation theory predictions and the low energy constants F, \bar{l}_3 and \bar{l}_4 are evaluated with small statistical errors. result 1.00 0.95 0.61 asserted 1.00

What kind of paper the classifier reports this is: simulation study 1.00 · evidence from simulation or computation 1.00.

What was examined, and what was not

sources read
  • arXiv listing abstract v2
not examined
  • full text
  • figures and tables
  • references and citing sentences
  • relations to other papers
levels run
  • P abstract
not run
  • S sections
  • T body sentences
  • C citations
  • R cross-paper relations

Every quoted sentence on this record is the authors'. Every label is from Pith's closed vocabulary, with the probability shown beside it. Every number in the measurements table is the authors', checked by code. Record as JSON

Figures
Figure 1. (a) PCAC quark mass amPCAC as function of aµ and (b) Sommer parameter (r0/a) as functions of (aµ) 2 . The solid line in subfigure (b) represents a linear fit in (aµ) 2 to the data. called Γ-method as described in Ref. [69] and cross-checked with a bootstrap analysis and a jack-knife analysis of blocked data. For a detailed description of our error analysis we refer to a forthcoming paper of our collaboration [34]. 3.2 Force parameter r0 In simulations of the quenched approximation of lattice QCD, the Sommer pa￾rameter r0 [70] with a value of 0.5 fm, was widely used to set the lattic
Figure 2. In (a) we show (amPS) 2/(aµ) as a function of aµ. In addition we plot the χPT fit with Eq. (5) to the data from the lowest four µ-values. In (b) we show (amPS) 2 as a function of aµ. Here we present two χPT fits with Eq. (5), one taking all data points and one leaving out the point at the largest value aµ = 0.015. In both figures (a) and (b) we show finite size corrected (L → ∞) data points. (mPS) 2 − (m0 PS) 2 = O(a 2Λ 4 QCD) = O(p 4 ), from which it follows that to the order we have been working the effects of the pion mass splitting do not affect, in par￾ticular, the finite
Figure 3. We show afPS as a function of aµ together with fits to χPT formula Eq. (6). We present two fits, one taking all data and one leaving out the point at the largest value aµ = 0.015. We show finite size corrected (L → ∞) data points. into account, shows that the uncharged pseudo scalar meson is about 20% lighter than the charged one. We obtain am± PS = 0.1359(7), am0 PS = 0.111(11), or, expressed differently, r 2 0 ((m0 PS) 2 − (m± PS) 2 ) = c(a/r0) 2 with c = −4.5(1.8). This coefficient is a factor of 2 smaller than the value found in quenched investi￾gations [28]. Note that t

Figures reproduced from arXiv:hep-lat/0701012 by the authors.

References
How later papers used it

61 later papers cite it

Evidence: substantial. 32 of the citing papers share no author with it, and 7 of those carry a confident label from Pith's citation classifier.

Pith read 20 of the 104 citing sentences, choosing only those whose citation markers point at this paper alone, and wrote the paraphrases from them. The summary, the counts, the evidence level, the links and the quotes come from code; the summary counts the papers listed below. A later paper is listed under a label when Pith's citation classifier gives one of its sentences that label at probability 0.8 or more (0.95 for support), Pith's reading finds that the citing paper itself did it, and a second check confirms the paraphrase against the sentence. This section lists uses only. It does not list criticism or disagreement, so their absence here says nothing about how the paper has been received. Papers that share an author with this one are listed apart, under "By the same authors". A quote is the citing authors' own sentence, shown when their paper's licence is CC BY, CC BY-SA or CC0. JSON

Uses its method or tool 1 of 6 papers with this label · 1 by the same authors

By the same authors

Uses its data or values 2 of 10 papers with this label · 1 by the same authors

By the same authors

  • Used the physical light quark mass previously determined in this paper to extrapolate its pseudoscalar meson mass down to that point.

    Twisted mass QCD in the charm sector (2007) (same authors) · 0.91 · paraphrased: the licence does not allow quoting

Compares with it 2 of 5 papers with this label · 2 by the same authors

By the same authors

How the literature uses it

61 arXiv papers cite this in 104 sentences

Each citing sentence is taken from the citing paper's LaTeX source and typed by Pith's classifier. A label appears at probability 0.8 or more (0.95 for support), a contest at 0.7 or more; the quoted words are the citing authors'. How each label was checked · JSON

Papers citing it · 61 sentences without a confident label

Uses its method or tool 6 papers
Uses its data or values 10 papers
Compares with it 5 papers
  • In this approach $f_π$ turns out to be about 10-15 % larger than in [ref] , at the matching point, but also the error-bars are of the same order of magnitude, and therefore still compatible.

    Results from overlap valence quarks on a twisted mass sea · Results in the pion sector · 1.00

  • The main differences are: i) the use of a single lattice spacing ($a ≃0.09$ fm) both in the present work and in Refs. [ref] , while the results obtained from two lattice spacings ($a ≃0.07$ and $0.09$ fm) are taken into account in Ref. [ref] ; ii) a better statistical accuracy of the data for $M_π$ and $f_π$ in Refs. [ref] due to the use of all the (correlated) trajectories produced by the ETM collaboration with respect to the present work in which only a subset of 240 (uncorrelated) trajectories are employed; iii) the presence of the results of the run $R_1$ at $M_π≃260$ in the present work and in Ref. [ref] at variance with Refs. [ref] .

    Electromagnetic form factor of the pion from twisted-mass lattice QCD at Nf=2 · Chiral fits · 1.00

  • By looking at table [ref], we see that the value of the pion decay constant in lattice units as obtained from the C-PQChPT fit is in agreement, at the level of 1.4 $σ$, with the result of our previous study, $af_π=0.0576(7)$ [ref] .

    Light quark masses and pseudoscalar decay constants from Nf=2 Lattice QCD with twisted mass fermions · Physical results · 1.00

Cites it as support 1 paper
Cites it as background 10 papers
  • Twisted mass lattice QCD [ref] has recently emerged as a viable fermion discretization method for computing gauge configurations with two flavors of light quarks, ($up$ and $down$) [ref] and hopeful prospects of $2+1+1$ ($up$, $down$, $strange$ and $charm$) flavors of dynamical sea fermions [ref] in the chiral regime.

    pi-pi Scattering in Twisted Mass Chiral Perturbation Theory · Introduction · 0.94

  • A comprehensive simulation programme has been undertaken by the European Twisted Mass (ETM) collaboration, including simulations with dynamical up and down quarks ($N_f=2$) [ref] and with dynamical up, down, strange and charm quarks ($N_f=2+1+1$) [ref] .

    Overlap valence quarks on a twisted mass sea: a case study for mixed action Lattice QCD · Wilson Twisted Mass Fermions · 0.81

  • The twisted mass formulation of Lattice QCD (see for example references [ref] ) is being studied extensively with $N_f=2$ dynamical flavours by the European Twisted Mass (ETM) collaboration [ref] .

    Status of ETMC simulations with Nf=2+1+1 twisted mass fermions · Introduction · 0.99

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.