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
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
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
- 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.
The abstract, sentence by sentence
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
| # | sentence | role | the authors' own work | states the main conclusion | strength |
|---|---|---|---|---|---|
| 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 was examined, and what was not
- arXiv listing abstract v1
- 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
Figures
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
3 later papers by other authors are listed below: 3 used its data or values. 3 more papers by the same authors are listed: 3 used its data or values.
Uses its data or values 6 of 24 papers with this label · 3 by the same authors
-
Used r_0=0.465±0.012 to relate to this scale determination through the Sommer scale for matching purposes.
Aspects of SU(3) baryon extrapolation (2009) · 0.99 · paraphrased: the licence does not allow quoting
-
Referred to this paper for the reported lattice scales and quark masses of all sets.
Running mass, effective energy and confinement: the lattice quark propagator in Coulomb gauge (2012) · 0.98 · paraphrased: the licence does not allow quoting
-
Performed its numerical calculations using 400 coarse MILC gauge configurations with V=20^3×64 and a^-1≈1.6 GeV.
Form factors for rare B decays: strategy, methodology, and numerical study (2009) · 0.89 · paraphrased: the licence does not allow quoting
By the same authors
-
Used gauge configurations provided by the MILC Collaboration that incorporate the effects of three light quarks: up, down, and strange.
Quarkonium mass splittings in three-flavor lattice QCD (2009) (same authors) · 0.89 · paraphrased: the licence does not allow quoting
-
Calculated the lattice spacing by measuring the ratio r_1/a from the static quark potential and using the value r_1=0.318 fm.
HISQ action in dynamical simulations (2009) (same authors) · 0.94 · paraphrased: the licence does not allow quoting
-
Used the values r_0 = 0.469(7) fm and r_1 = 0.318(7) fm, known to about 2 percent.
Toward a precise determination of Tc with 2+1 flavors of quarks (2007) (same authors) · 0.98 · paraphrased: the licence does not allow quoting
How the literature uses it
53 arXiv papers cite this in 85 sentences
Papers citing it · 42 sentences without a confident label
Uses its method or tool 3 papers
-
In practice, $ππ$ four-point functions were fit by varying the minimum fitting distances $D_min$, and with the maximum distance $D_max$ either at $T/2$ or where the fractional statistical errors surpassed about $20%$ for two sequential time slices [ref] .
Lattice QCD study of the s-wave $\pi\pi $ scattering lengths in the I=0 and 2 channels · $I=2$ channel · 0.98
-
Also in Table [ref] are the tadpole factors $u_0$ [ref] , determined from the mean plaquette and used to improve the gauge-configuration actions [ref] .
Tuning Fermilab Heavy Quarks in 2+1 Flavor Lattice QCD with Application to Hyperfine Splittings · Gauge Configurations and Related Parameters · 0.96
-
Our code is based on the publicly available MILC package[] for improved staggered quarks [ref] on a Symanzik gauge background.
More about QCD on compact spaces · Introduction · 1.00
Uses its data or values 24 papers
-
Our calculations are carried out on several ensembles of lattice gauge fields with a realistic sea of 2+1 flavors, made available by the MILC Collaboration [ref] .
B-Bbar Mixing and Matching with Fermilab Heavy Quarks · Long-Distance Matrix Elements · 0.98
-
For all hadrons and resonances, we use the pion mass- and lattice spacing-dependence given in Refs. [ref] and parametrized in Ref. [ref] .
Is there still any Tc mystery in lattice QCD? Results with physical masses in the continuum limit III · Comparison between HRG model and lattice results · 0.89
-
We have fitted the quark (pion) mass and lattice spacing dependence of the $ρ, K^*, ϕ, N$ and $Ω$ masses obtained in Refs. [ref] by a simple Ansatz [equation] The values of the fit parameters $m_0$, $a_1, a_2, b_1$ and $b_2$ are given in Table [ref].
QCD Equation of State and Hadron Resonance Gas · Hadron masses on the lattice · 0.96
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
-
The Fermilab Lattice and MILC Collaborations have reported several calculations [ref] based on ensembles of lattice gauge fields with 2+1 flavors of sea quarks, generated by the MILC Collaboration [ref] .
Tuning Fermilab Heavy Quarks in 2+1 Flavor Lattice QCD with Application to Hyperfine Splittings · Introduction · 0.88
-
The dynamical lattice QCD calculations of hadron interactions to date, have either been performed with Wilson fermions or a mixed lattice action [ref] of domain-wall valence fermions [ref] and the Asqtad improved [ref] rooted staggered MILC configurations [ref] .
pi-pi Scattering in Twisted Mass Chiral Perturbation Theory · Introduction · 0.95
-
Previous work on excited and high spin mesons in the light quark sector can be found in Refs. [ref] for quenched calculations and in Refs. [ref] for the dynamical case.
Low and high spin mesons from $N_f=2$ Clover-Wilson lattices · 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.