Pith. Sign in

REVIEW

Light hadrons with improved staggered quarks: approaching the continuum limit

High Energy Physics - Lattice · Submitted 24 Feb 2004 · 29 references

Pith's verdict: A solid, unflashy MILC simulation report whose value is in the systematic checks and the new fine ensembles, with the r0/r1 continuum extrapolation the one place a referee should push.

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 →

· 3 major objections · 9 minor · Ask about this paper · The compiled record · Read the full review →
Challenge this review Comment (0) Read on arXiv ↗
Cite ▾
arXiv IDhep-lat/0402030 (version 1)
Pith IDI3QTEYWI, Pith's permanent record of this paper and its review
CitationC. Aubin, C. Bernard, C. DeTar, D. Toussaint, E.B. Gregory, J.E. Hetrick, J. Osborn, R. Sugar, Steven Gottlieb, U.M. Heller. "Light hadrons with improved staggered quarks: approaching the continuum limit." arXiv:hep-lat/0402030 (2004). https://pith.science/paper/I3QTEYWI
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
The paper

The paper claims

Adding a finer lattice lets three-flavor QCD simulations extrapolate to zero spacing, with a stated error budget.

The reading

This paper reports the extension of a long-running program of lattice QCD simulations with three flavors of dynamical quarks, using an improved staggered (Kogut-Susskind) action, to a second and finer lattice spacing of about 0.09 fm, alongside extended runs at the older 0.12 fm spacing and lighter quark masses. The point is not a single number but the ability to do a controlled continuum extrapolation: with two lattice spacings in hand, the authors can ask how much of their answer is an artifact of the discretized grid. They use the static quark potential to define a scale, $r_1$ (fixed by $r_1^2 F(r_1) = 1.00$), and quote $r_1 = 0.317(7)(3)$ fm and a potential shape parameter $r_0/r_1 = 1.474(7)(18)$ at the physical pion mass in the continuum limit. Along the way they check the things that could quietly corrupt such a program: single versus double precision arithmetic, the molecular-dynamics step size, the physical spatial volume, and autocorrelations in simulation time. They also report that at their lightest quark masses the meson propagators in the $0^{++}$, $1^{+}$ and $0^{-}$ channels show the influence of two-meson decay thresholds, and they extract preliminary masses for radially excited pseudoscalar states.

What carries the argument

The load-bearing device is the scale $r_1$, defined implicitly by $r_1^2 F(r_1) = 1.00$ on the static quark potential. It is cheap to compute, insensitive to valence quark mass, and therefore serves as the common ruler that lets runs at different couplings and quark masses be compared and converted to physical units. Around it sit three supporting pieces: the smoothed fit to $r_1/a$ across couplings and masses, which suppresses per-run fitting noise; the interpolating form for the potential that includes a lattice correction term $λ(V_{free}(r) - 1/r)$ used at short distances; and the two-point linear extrapolation in an assumed error variable ($α_s a^2$, or $α_s^2 a^2$, or the measured taste-splitting ratio) to reach zero lattice spacing. The second mechanism carrying physical weight is the two-exponential and three-exponential fit to meson propagators, $C(t) = A_0(e^{-M_0 t} + e^{-M_0(T-t)}) + A_1(e^{-M_1 t} + e^{-M_1(T-t)})$ and its oscillating-state extension, which is what makes the excited-state and decay-channel claims visible at all rather than buried in a single-exponential plateau.

Extended reading notes

Core claim

The central claim is that a second, finer lattice spacing turns the collaboration's three-flavor simulation set into a genuine two-point continuum extrapolation with a quantified systematic budget. Using the static potential fitted to $V(r) = C + $σ$ r - $α$/r + $λ$(V_{free}(r) - 1/r)$, the authors obtain $r_1/a$ for each run and then fit a smoothed form $��log(r_1/a) = C_{00} + C_{10}(10/g^2 - 7.0) + C_{01}(2m_{u,d} + m_s) + C_{20}(10/g^2 - 7.0)^2$ so that run-to-run statistical fluctuations do not propagate into the scale conversion. Setting the scale with the Υ 2S–1S and 1P–1S splittings gives $r_1 = 0.324(4)$ fm on coarse lattices and $0.320(4)$ fm on fine lattices; extrapolating away the discretization error gives $r_1 = 0.317(7)(3)$ fm. For the potential shape they fit coarse and fine results simultaneously to a constant plus a linear quark-mass term plus a term in $a^2$α_s, finding $r_0/r_1 = 1.476(7) - 0.049(10)(M_π/M_ρ)^2 - 0.12(4)(a/r_1)^2$α_s(a)/α_s(0.12 fm), and hence $r_0/r_1 = 1.474(7)(18)$ at the physical mass in the continuum. A separate and physically striking claim is that the light pseudoscalar propagators support two-state fits, yielding an excited $0^{-+}$ mass that extrapolates to $1362(41)(247)$ MeV, consistent with the $π(1300)$; the excited kaon extrapolates to $1527(46)(68)$ MeV, consistent with the $K(1460)$; and the excited $s̄s$ state to $1645(40)(145)$ MeV, consistent with the $η(1440)$.

Load-bearing premise

The continuum extrapolation rests on a straight line through two lattice spacings in an error variable the authors choose by hand, taking an average ratio of 0.4 among three plausible scalings, so the quoted continuum value inherits whatever that choice does not cancel.

Editorial extensions

If this is right

  • Downstream analyses that need a lattice scale — quark masses, $f_π$, $f_K$, $α_s$, heavy-light decay constants — can use a scale whose discretization error has been estimated by an actual second lattice spacing rather than by assumption alone.
  • The quoted $r_0/r_1 = 1.474(7)(18)$ supplies a continuum-anchored shape parameter for the static potential, which is what heavy-quark and potential-model calculations compare against.
  • The demonstration that single-precision propagator computations agree with double precision at $am_{u,d} = 0.005$ means the collaboration can push to lighter quark masses without paying a large precision penalty.
  • The explicit null results for finite-spatial-size effects between $L=20$ and $L=28$ (2.4 fm and 3.4 fm, with $M_π L$ of 4.48 and 6.27) justify treating the smaller volume as adequate at these masses.
  • The identification of two-meson threshold effects in the $0^{++}$, $1^{+}$ and $0^{-}$ channels, and the consistency of the excited pseudoscalar masses with $π(1300)$, $K(1460)$ and $η(1440)$, argues that the fourth-root staggered formulation is not introducing an unphysical length scale above roughly two lattice spacings in these channels.

Reading between the lines

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

  • A third lattice spacing, say around 0.06 fm, is the obvious next test: with only two points the linearity of the continuum extrapolation in the chosen error variable is assumed, not checked, so a third point would either confirm the slope or expose curvature that the current error bar does not cover.
  • Because the coarse and fine runs used different strange quark masses ($am_s = 0.050$ versus $0.031$), the extrapolation in lattice spacing is partly doing the work of a strange-mass correction; a reader should expect the two effects to be entangled until a run matches the physical strange mass at both spacings.
  • The excited-state masses come from wall sources chosen to suppress excited-state contamination, so a different source shape (e.g. point or derivative sources) should reproduce the same $π(1300)$-like and $K(1460)$-like masses if the identification is robust rather than a fitting artifact.
  • The reported autocorrelation times, largest for the pion mass and larger on fine lattices than coarse, suggest that future fine-lattice runs at lighter quark masses will need substantially more stored configurations per unit of statistical error than the coarse runs did.
Original abstract

We have extended our program of QCD simulations with an improved Kogut-Susskind quark action to a smaller lattice spacing, approximately 0.09 fm. Also, the simulations with a approximately 0.12 fm have been extended to smaller quark masses. In this paper we describe the new simulations and computations of the static quark potential and light hadron spectrum. These results give information about the remaining dependences on the lattice spacing. We examine the dependence of computed quantities on the spatial size of the lattice, on the numerical precision in the computations, and on the step size used in the numerical integrations. We examine the effects of autocorrelations in "simulation time" on the potential and spectrum. We see effects of decays, or coupling to two-meson states, in the 0++, 1+, and 0- meson propagators, and we make a preliminary mass computation for a radially excited 0- meson.

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, 5,556 input tokens, $0.0002 reported.

What the abstract itself states, as probabilities: claims first 0.04 · states limitation 0.74 · human participants 0.01 · controlled comparison 0.06 · randomized assignment 0.01 · states quantitative result 0.45 · compares with prior results 0.48 · mentions code or data availability 0.02.

The abstract, sentence by sentence

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

We have extended our program of QCD simulations with an improved Kogut-Susskind quark action to a smaller lattice spacing, approximately 0.09 fm. Also, the simulations with a approximately 0.12 fm have been extended to smaller quark masses. In this paper we describe the new simulations and computations of the static quark potential and light hadron spectrum. These results give information about the remaining dependences on the lattice spacing. We examine the dependence of computed quantities on the spatial size of the lattice, on the numerical precision in the computations, and on the step size used in the numerical integrations. We examine the effects of autocorrelations in "simulation time" on the potential and spectrum. We see effects of decays, or coupling to two-meson states, in the 0++, 1+, and 0- meson propagators, and we make a preliminary mass computation for a radially excited 0- meson.

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 have extended our program of QCD simulations with an improved Kogut-Susskind quark action to a smaller lattice spacing, approximately 0.09 fm. contribution 0.76 0.95 0.05 asserted 1.00
s1 Also, the simulations with a approximately 0.12 fm have been extended to smaller quark masses. contribution 0.56 0.91 0.03 asserted 1.00
s2 In this paper we describe the new simulations and computations of the static quark potential and light hadron spectrum. contribution 1.00 0.90 0.07 asserted 0.90
s3 These results give information about the remaining dependences on the lattice spacing. interpretation 0.91 0.71 0.16 asserted 0.93
s4 We examine the dependence of computed quantities on the spatial size of the lattice, on the numerical precision in the computations, and on the step size used in the numerical integrations. method 0.74 0.93 0.03 asserted 0.86
s5 We examine the effects of autocorrelations in "simulation time" on the potential and spectrum. method 0.57 0.84 0.03 asserted 0.96
s6 We see effects of decays, or coupling to two-meson states, in the 0++, 1+, and 0- meson propagators, and we make a preliminary mass computation for a radially excited 0- meson. result 1.00 0.97 0.12 asserted 0.99

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 v1
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. FIG. 1: A “shape parameter” for the static potential,
Figure 2. FIG. 2: Pseudoscalar masses as a function of minimum distanc
Figure 3. FIG. 3: Vector meson masses as a function of minimum distance
Figure 4. FIG. 4: Nucleon masses. The blue diamonds and octagons are qu
Figure 5. FIG. 5: Single elimination jackknife masses for the pion, fr
Figure 6. FIG. 6: Normalized autocorrelations for the
Figure 9. FIG. 9: Fit plot showing ground state and excited state masse
Figure 10. FIG. 10: Ground state and excited state 0
Figure 11. FIG. 11: Two-state fits to three-flavor kaon pseudoscalar pro
Figure 12. FIG. 12: Three-state fits to three-flavor kaon pseudoscalar p
Figure 13. FIG. 13: Quenched kaon mass fit plot showing ground state and fi
Figure 14. FIG. 14: Summary of fits of kaon propagators. Ground state and
Figure 15. FIG. 15: The “J” parameter. The red squares and crosses are th
Figure 16. FIG. 16: The “big picture”. Crude continuum and chiral extra

Showing 14 of the paper's 16 figures. See all figures in the PDF ↗ Figures reproduced from arXiv:hep-lat/0402030 by the authors.

References
How later papers used it

58 later papers cite it

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

Pith read 9 of the 92 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 data or values 6 of 24 papers with this label · 3 by the same authors

By the same authors

How the literature uses it

53 arXiv papers cite this in 85 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 · 42 sentences without a confident label

Uses its method or tool 3 papers
Uses its data or values 24 papers
Compares with it 1 paper
  • In Figure [ref] we compare our estimates of the $Δ$, $Λ$ and $Ξ$ baryon masses shown as lines with available lattice data from the MILC collaboration [ref] .

    QCD Equation of State and Hadron Resonance Gas · Hadron masses on the lattice · 0.93

Cites it as background 3 papers

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.