REVIEW 2 major objections 5 minor 2 cited by
Computing $1/m_Q$ and $1/m_Q^2$ corrections to the static potential with lattice gauge theory using gradient flow
T0 review · 2 major / 5 minor · reviewed 2026-08-12 · deepseek-v4-flash
Pith's one-line read Tree-level improvement removes the dominant lattice and flow-time errors from the Wilson loops behind the 1/m_Q^2 corrections to the static potential.
desk verdict A solid, clearly preliminary proceedings paper showing tree-level improvement works for the static potential and plausibly extends to field-strength correlators; the main open question is whether the improvement coefficient transfers. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The central object is the generalized Wilson loop with two field-strength insertions, Eq. (3), whose ground-state limit yields the correlator of clover-defined chromoelectric or chromomagnetic fields in the flux tube of a static quark-antiquark pair. The argument is carried by tree-level improvement. For the static potential the ansatz (7) adds the correction term $\tilde c^{(0)}(4\pi G(r,t_f)-1/r)$, with $G(r,t_f)$ the lattice tree-level gluon propagator at flow time $t_f$; subtracting this term and the fitted $r$-independent shift $V_c(t_f,a)$ removes most of the discretization and flow dependence. For two-field-strength correlators the same idea enters through Eq. (12), where $\tilde c (4\pi/C_F g^2)(C^{\mathrm{lattice}}_{\mathrm{tree\ level}} - C^{\mathrm{continuum}}_{\mathrm{tree\ level}})$ is subtracted from the numerical correlator, with the coefficient $\tilde c$ taken from the static-potential fit. Gradient flow is the enabling regulator: flowed correlators need no extra renormalization, and the flow time $t_f$ regulates the logarithmic divergence of chromomagnetic insertions.
What would settle it
Compute the tree-level-improved correlator (12) on a fourth lattice ensemble with roughly half the lattice spacing of the coarsest ensemble used here at the same physical $r$ and $t_f$; if the small-$t$ data still shift by more than the statistical errors when $a$ and $t_f$ are reduced, the tree-level subtraction is not capturing the dominant artifact.
Extended reading notes
Core claim
The central claim is that the dominant finite-lattice-spacing and finite-flow-time artifacts in generalized Wilson loops with two field-strength insertions are tree-level effects, and that subtracting the difference between the lattice and continuum tree-level results removes them. For the static potential the paper fits the ansatz (7), $V^{(0)}(r,t_f) = -c^{(0)}/r + \sigma r + V_c(t_f,a) + \tilde c^{(0)}(4\pi G(r,t_f) - 1/r)$, where $G(r,t_f)$ is the flowed lattice tree-level propagator; removing the shift $V_c$ and the correction term collapses data from three ensembles and four flow times onto one curve down to $r \approx 0.05$ fm. The same strategy is applied to the two-chromoelectric-field correlator $\langle \Sigma^+_{g,r}|E_z(t,0)E_z(0,0)|\Sigma^+_{g,r}\rangle_c$ via Eq. (12), producing improved correlators whose small-$t$ behavior is essentially independent of flow time. The paper also reports that the Gromes and first BBMP relations are satisfied within statistical errors at finite flow time, in contrast to an earlier multilevel computation that used an approximate multiplicative renormalization; this supports the interpretation that the earlier violations came mostly from the renormalization prescription rather than from lattice artifacts.
Load-bearing premise
The method assumes the dominant finite-lattice-spacing and finite-flow-time artifacts are exactly the tree-level lattice-versus-continuum difference, with the remaining deviation absorbed by one r-independent shift and a single fitted coefficient that is the same for the static potential and for the two-field-strength correlators.
Editorial extensions
If this is right
- The static potential can be extrapolated to $a \to 0$ and $t_f \to 0$ while including data at $r \approx 0.05$ fm, where unimproved data are dominated by flow effects.
- The tree-level improved $\langle E_z E_z \rangle_c$ correlators agree across flow times at small $t$, so the integrals defining the potential corrections can be evaluated reliably at short separations.
- The Gromes and first BBMP relations, which need no matching coefficients, are satisfied within statistical errors at finite flow time, indicating that earlier violations came mostly from the approximate renormalization prescription.
- Spin-independent $1/m_Q^2$ corrections can be obtained directly from a combined continuum and zero-flow-time extrapolation, while spin-dependent corrections still require NLO matching coefficients to convert from the gradient-flow scheme to the $\overline{\mathrm{MS}}$ scheme.
Reading between the lines
- If the fitted coefficient $\tilde c^{(0)}$ is genuinely transferable from the static potential to every two-field-strength correlator, the same subtraction can be applied to the chromomagnetic correlators that enter the spin-dependent potentials, extending the method beyond the $\langle E_z E_z \rangle_c$ example shown.
- The agreement of improved data at $r \approx 0.05$ fm suggests the residual nonperturbative lattice artifacts are subleading at surprisingly short distances; a test would be to check whether $\tilde c^{(0)}$ stays constant when computed on a finer lattice.
- Because gradient-flow regularization avoids the leftover $O(g^4)$ and $O(g^6)$ terms of the earlier multiplicative renormalization, the Gromes and BBMP relations could be promoted from consistency checks to constraints used in the fits that determine the spin-dependent potentials.
Signed reviews
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper reports preliminary lattice QCD results for the spin-dependent and spin-independent O(1/m_Q^2) corrections to the static potential, obtained from generalized Wilson loops with two field-strength insertions renormalized via gradient flow. The main methodological content is a tree-level improvement prescription: for the static potential, the lattice data are fitted with the Cornell ansatz plus a correction term proportional to the difference between lattice and continuum tree-level expressions, Eq. (7); for field-strength correlators, an analogous subtraction is proposed in Eq. (12), with the coefficient \tilde c either fitted globally or identified with the static-potential parameter \tilde c^{(0)}. The static-potential improvement is demonstrated on three ensembles and four flow times. For the 1/m_Q^2 potentials, the paper shows fits for four selected contributions on one ensemble (and, for one of them, all ensembles and flow times) and checks the Gromes and BBMP relations at one flow time. The authors explicitly state that the results are preliminary and that continuum and zero-flow-time extrapolations of the potential corrections are still ongoing.
Significance. If the tree-level improvement program works as claimed, it would provide a practical route to controlled small-r and small-t data for the 1/m_Q and 1/m_Q^2 static-potential corrections, addressing a long-standing difficulty with renormalization and signal-to-noise in field-strength correlators. The static-potential demonstration in Sec. 4.2, where flow-time dependence is visibly reduced after the subtraction, is encouraging and is supported by data on three ensembles. The use of gradient flow also avoids the approximate Huntley-Michael renormalization that was a concern in earlier work. However, the central new step for the correlators, Eq. (12), rests on an assumption about the transferability of the coefficient \tilde c from the static potential to the field-strength correlators, and the manuscript does not yet provide the quantitative evidence needed to validate that assumption. As it stands, this is a promising proceedings contribution rather than a completed validation of the method.
major comments (2)
- [Sec. 4.4, Eq. (12)] The identification of the correlator-improvement coefficient \tilde c with the static-potential fit parameter \tilde c^{(0)} from Eq. (7) is not demonstrated. The text states that using \tilde c = \tilde c^{(0)} leads to improved correlators of 'similar quality' to a global fit of \tilde c, but it reports neither the fitted value of \tilde c, its uncertainty, nor a quantitative comparison with \tilde c^{(0)}. Since Eq. (12) is an exact relation only at tree level, and the flowed field-strength insertions carry their own flow-time and renormalization dependence, the equality \tilde c = \tilde c^{(0)} is a load-bearing assumption for the claim that the improved small-t correlator data are trustworthy. Please report the global correlator fit for \tilde c and compare it with \tilde c^{(0)} explicitly, or recast the improved-correlator results as conditional on this transferability assumption.
- [Sec. 4.4, first paragraph] The statement that NRQCD matching coefficients 'differ from 1 at O(α2)' appears inconsistent with the cited NLO matching calculation of Ref. [19], which computes one-loop, i.e. O(α_s), corrections to the relevant matching coefficients. If the sentence is intended to justify ignoring matching coefficients in a tree-level improvement, it should be phrased as a tree-level statement; as written, it is a statement about the matching coefficients themselves and is inaccurate.
minor comments (5)
- [Sec. 4.3, Fig. 3] The paper claims that the long-range parameter g_Λ' is 'determined quite accurately' and that the long-range term is 'crucial' for V_LS^(1,1), but it does not report the fitted parameter values, uncertainties, or χ²/DOF for any of the ansaetze in Fig. 3. Please provide these numbers in a table or in the text so that the claims can be checked.
- [Sec. 4.3, Fig. 4 and surrounding text] The conclusion that violations of the Gromes relation in Ref. [14] 'were mostly caused by approximations in the Huntley-Michael renormalization prescription' goes beyond what can be concluded from a single ensemble and a single flow time without continuum extrapolation. The authors correctly note the lack of extrapolation, but the interpretive statement should be softened or accompanied by a controlled comparison.
- [Sec. 4.4, Fig. 5] The tree-level improved correlator data are shown for only one spatial separation, r ≈ 0.24 fm, and the claim about reliable data at 'significantly smaller r and t' is explicitly deferred to future work. It would be helpful to state clearly in Sec. 4.4 that the small-r part of the claim is not yet demonstrated in this paper.
- [Sec. 4.2 and Fig. 2] The two-step fitting procedure involving first a fit with fixed c~(0) = 0 (Fig. 2, left) and then the full ansatz with c~(0) free (Fig. 2, right) is described only briefly. A more explicit description of how c~(0) is determined and whether the small-r points are included in that determination would help the reader assess the significance of the collapse shown in the right panel.
- [Throughout] There are minor typographical errors, e.g. 'apprroximately' in Sec. 4.3, and the caption of Fig. 2 refers to a vertical grey line without defining its meaning in the caption; these should be corrected.
Circularity Check
No circularity found: the tree-level improvement is fit-based, but no predicted quantity is defined by its own fit input, and the c~ transfer to correlators is an untested assumption rather than a tautology.
full rationale
Walking the derivation chain: Eq. (7) is a fit ansatz for the static potential in which c~(0) is a nuisance parameter fitted together with c(0), sigma, and V_c. The resulting collapse of the improved static-potential data is therefore partly a consequence of the fit, but the paper does not present that collapse as an independent prediction, and the central claim is checked against external relations (Gromes and BBMP relations, and the tree-level result c(0)/4r^3 at intermediate r). Eq. (12) defines the tree-level improved field-strength correlator by subtracting the lattice-minus-continuum tree-level difference with a coefficient c~; using c~=c~(0) from the static-potential fit is a transferability assumption. It is not circular by construction because the correlator data are not used to fix c~(0), and the potential corrections extracted from these correlators are separate targets. The only self-citation, Ref. [15] for the fitting strategy of the integrals in Eq. (2), is not load-bearing for the new tree-level improvement claim and does not close a definitional loop. No equation defines its target in terms of the quantity it is supposed to predict, and no fitted parameter is relabeled as a prediction. The unquantified transferability of c~ in Sec. 4.4 is a potential systematic-error concern, but it is not a circularity. Score 0.
Assumptions & free parameters
free parameters (6)
- c^(0) =
not quoted
- sigma =
not quoted
- V_c(t_f,a) =
one value per ensemble and flow time
- c~(0) =
not quoted
- g_Lambda' =
not quoted
- additional f(r) fit coefficients for r V_LS^(2,0), r V_LS^(1,1), V_p2^(1,1), V_S12^(1,1) =
unknown
assumptions (5)
- domain assumption The heavy-quark potential admits the 1/m_Q expansion of Eq. (1), and V^(1), V^(2)_SD and V^(2)_SI can be obtained from Wilson-loop correlators via Eqs. (2)-(3).
- domain assumption Gradient-flowed gauge fields produce correlators that need no additional renormalization for t_f > 0 and have a well-defined small-flow-time continuum limit.
- domain assumption The ground state |Sigma+_g,r> dominates the large-Delta_t limit in Eq. (3), so the generalized Wilson-loop ratio isolates the desired correlator.
- ad hoc to paper Tree-level lattice perturbation theory expressions, Eq. (8) and Ref. [32], capture the dominant finite-a and finite-t_f artifacts, so the correction terms in Eqs. (7) and (12) remove most systematic errors.
- domain assumption Physical units are set by r_0 = 0.5 fm.
Cite this review
Pith. "Pith review of Computing $1/m_Q$ and $1/m_Q^2$ corrections to the static potential with lattice gauge theory using gradient flow." pith.science (2026). https://pith.science/paper/NSNCLG34
@misc{pith2026241111640,
author = {Pith},
title = {Pith review of: Computing $1/m_Q$ and $1/m_Q^2$ corrections to the static potential with lattice gauge theory using gradient flow},
year = {2026},
howpublished = {\url{https://pith.science/paper/NSNCLG34}},
note = {Machine review of arXiv:2411.11640}
}
abstract
We present selected preliminary lattice gauge theory results for $O(1/m_Q)$ and $O(1/m_Q^2)$ corrections to the static potential. These results are based on Wilson loops with two field strength insertions, which we renormalize using gradient flow. We explore tree level improvement to reduce systematic errors in the Wilson loops due to the finite lattice spacing and flow time, in particular at small temporal and spatial separations.
Figures
Figures from the paper (2 more)
Forward citations
Cited by 2 Pith papers
-
Hybrid spin-dependent and hybrid-quarkonium mixing potentials at order $(1 /m_Q)^1$ from SU(3) lattice gauge theory
First SU(3) lattice computation of the four order-(1/m_Q)^1 hybrid spin-dependent and hybrid-quarkonium mixing potentials at a single lattice spacing.
-
Spectra of light and heavy mesons with $J \le 5$ in a relativistic Bethe-Salpeter approach
First J=4,5 meson Bethe-Salpeter tensor bases plus exploratory rainbow-ladder spectra and Regge trajectories from light through bottom quarkonia.
Reference graph
Works this paper leans on
-
[19]
N. Brambilla and X. P. Wang, JHEP06, 210 (2024) [arXiv:2312.05032 [hep-ph]]
arXiv 2024
- [1]
-
[2]
A.Barchielli,E.MontaldiandG.M.Prosperi,Nucl.Phys.B 296,625(1988)[erratum: Nucl. Phys. B303, 752 (1988)]
work page 1988
-
[3]
N. Brambilla, A. Pineda, J. Soto and A. Vairo, Nucl. Phys. B566, 275 (2000) [arXiv:hep- ph/9907240 [hep-ph]]
-
[4]
N. Brambilla, A. Pineda, J. Soto and A. Vairo, Phys. Rev. D63, 014023 (2001) [arXiv:hep- ph/0002250 [hep-ph]]
-
[5]
A. Pineda and A. Vairo, Phys. Rev. D63, 054007 (2001) [erratum: Phys. Rev. D64, 039902 (2001)] [arXiv:hep-ph/0009145 [hep-ph]]
arXiv 2001
-
[6]
N. Brambilla, A. Pineda, J. Soto and A. Vairo, Rev. Mod. Phys.77, 1423 (2005) [arXiv:hep- ph/0410047 [hep-ph]]
arXiv 2005
- [7]
Show all 32 references
-
[8]
Anzai, D
C. Anzai, D. Moreno and A. Pineda, Phys. Rev. D98, 114034 (2018) [arXiv:1810.11031 [hep-ph]]
2018 arXiv
-
[9]
Peset, A
C. Peset, A. Pineda and J. Segovia, Phys. Rev. D98, 094003 (2018) [arXiv:1809.09124 [hep-ph]]
2018 arXiv
-
[10]
de Forcrand and J
P. de Forcrand and J. D. Stack, Phys. Rev. Lett.55, 1254 (1985)
1985
-
[11]
Campostrini, K
M. Campostrini, K. Moriarty and C. Rebbi, Phys. Rev. Lett.57, 44 (1986)
1986
-
[12]
G. S. Bali, K. Schilling and A. Wachter, Phys. Rev. D56, 2566-2589 (1997) [arXiv:hep- lat/9703019 [hep-lat]]
1997
-
[13]
G. S. Bali and P. Boyle, Phys. Rev. D59, 114504 (1999) [arXiv:hep-lat/9809180 [hep-lat]]
1999 arXiv
-
[14]
Koma and M
Y. Koma and M. Koma, Nucl. Phys. B769, 79-107 (2007) [arXiv:hep-lat/0609078 [hep-lat]]
2007 arXiv
-
[15]
Eichberg and M
M. Eichberg and M. Wagner, PoSLATTICE2023, 068 (2024) [arXiv:2311.06560 [hep-lat]]
2024 arXiv
-
[16]
Lüscher, JHEP08, 071 (2010) [erratum: JHEP03, 092 (2014)] [arXiv:1006.4518 [hep- lat]]
M. Lüscher, JHEP08, 071 (2010) [erratum: JHEP03, 092 (2014)] [arXiv:1006.4518 [hep- lat]]
2010 arXiv
- [17]
-
[18]
Brambilla, V
N. Brambilla, V. Leino, J. Mayer-Steudte and A. Vairo, Phys. Rev. D109, 114517 (2024) [arXiv:2312.17231 [hep-lat]]
2024 arXiv
-
[20]
A. M. Eller and G. D. Moore, Phys. Rev. D97, 114507 (2018) [arXiv:1802.04562 [hep-lat]]. 9 Computing 1/𝑚𝑄 and 1/𝑚2 𝑄 corrections to the static potential Michael Eichberg
2018 arXiv
-
[21]
Rev.D105, 054503(2022)[arXiv:2111.00741[hep-lat]]
C.Schlosserand M.Wagner, Phys. Rev.D105, 054503(2022)[arXiv:2111.00741[hep-lat]]
2022 arXiv
-
[22]
Sciarra, C
A. Sciarra, C. Pinke, M. Bach, F. Cuteri, L. Zeidlewicz, C. Schäfer, T. Breitenfelder, C. Cza- ban, S. Lottini, P.F. Depta, CL2QCD (v1.1), Zenodo (2021),https://doi.org/10.5281/ zenodo.5121917
2021
-
[23]
Joswig, S
F. Joswig, S. Kuberski, J. T. Kuhlmann and J. Neuendorf, Comput. Phys. Commun.288, 108750 (2023) [arXiv:2209.14371 [hep-lat]]
2023 arXiv
-
[24]
Wolff [ALPHA], Comput
U. Wolff [ALPHA], Comput. Phys. Commun.156, 143-153 (2004) [erratum: Comput. Phys. Commun.176, 383 (2007)] [arXiv:hep-lat/0306017 [hep-lat]]
2004 arXiv
-
[25]
Ramos, Comput
A. Ramos, Comput. Phys. Commun.238, 19-35 (2019) [arXiv:1809.01289 [hep-lat]]
2019 arXiv
-
[26]
Barchielli, N
A. Barchielli, N. Brambilla and G. M. Prosperi, Nuovo Cim. A103, 59 (1990)
1990
-
[27]
Brambilla, M
N. Brambilla, M. Groher, H. E. Martinez and A. Vairo, Phys. Rev. D90, 114032 (2014) [arXiv:1407.7761 [hep-ph]]
2014 arXiv
-
[28]
Oncala and J
R. Oncala and J. Soto, Phys. Rev. D96, 014004 (2017) [arXiv:1702.03900 [hep-ph]]
2017 arXiv
-
[29]
Brambilla, D
N. Brambilla, D. Gromes and A. Vairo, Phys. Lett. B576, 314-327 (2003) [arXiv:hep- ph/0306107 [hep-ph]]
2003
-
[30]
Gromes, Z
D. Gromes, Z. Phys. C26, 401 (1984)
1984
-
[31]
Huntley and C
A. Huntley and C. Michael, Nucl. Phys. B286, 211-230 (1987)
1987
-
[32]
Fritzsch and A
P. Fritzsch and A. Ramos, JHEP10, 008 (2013) [arXiv:1301.4388 [hep-lat]]. 10
2013 arXiv
Reviewed August 12, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.