REVIEW 3 major objections 3 minor 52 references
Toward inclusive observables with staggered quarks: the smeared $R$~ratio
T0 review · 3 major / 3 minor · reviewed 2026-08-12 · deepseek-v4-flash
Pith's one-line read Smeared R ratio can be computed from staggered-quark lattices by splitting even and odd time slices, with low-energy results matching the standard parameterization.
desk verdict A clearly written pilot showing the first staggered-quark HLT smeared R ratio, with an honest list of limitations and one real, unquantified systematic in the parity-separation trick. 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 load-bearing object is the HLT algorithm, a convex-optimization spectral reconstruction that turns a Euclidean correlator $C(\tau)=\int_0^\infty d\omega\, e^{-\omega\tau}\rho_L(\omega)$ into a smeared spectral function $\tilde{\rho}_L(\omega^*)=\sum_\tau g_\tau(\omega^*)C(\tau)$. The coefficients $g_\tau$ are chosen to make the effective smearing kernel approximate a Gaussian while balancing statistical error through the covariance matrix and a tuned hyperparameter $\lambda$. The new piece for staggered quarks is the even/odd time-slice split: separate reconstructions are run on even and odd $\tau$ data, and Eq. (11) combines them as $\rho_\pm=\frac12(\rho_{\rm even}\pm\rho_{\rm odd})$ to eliminate the oscillating $1^+$ contributions.
What would settle it
Take a staggered ensemble with high statistics and a known R ratio, or a synthetic correlator built from a known two-parity spectrum, and reconstruct the smeared spectral function using only even, only odd, and combined even/odd time slices. If the linear combination $\rho_\pm = \frac{1}{2}(\rho_{\rm even}\pm\rho_{\rm odd})$ does not reproduce the known definite-parity smeared spectrum within errors, or if splitting the data in different ways changes the result by more than the quoted uncertainties, the even/odd separation is biased.
Extended reading notes
Core claim
On the paper's own terms, the central discovery is that the opposite-parity contamination in staggered vector-vector correlators can be handled by applying the HLT spectral reconstruction to even and odd time slices separately and then forming $\rho_\pm = \frac{1}{2}(\rho_{\rm even}\pm \rho_{\rm odd})$. Because a staggered correlator has spectral decomposition $C(\tau)=\sum_n (-1)^{n(\tau+1)} |\langle \Omega|O|n\rangle|^2 e^{-E_n\tau}/(2E_n)$, even and odd time slices carry the combinations $\rho_-+\rho_+$ and $\rho_--\rho_+$; the linear combination cancels the unwanted parity. The paper reports that the smeared $R$ ratio reconstructed this way on two physical-mass staggered ensembles agrees well with the comparison parameterization at low energies and departs above about 1 GeV, and it presents this as progress, not a final continuum result. The domain-wall reconstruction on one ensemble is included for qualitative comparison.
Load-bearing premise
The paper's method assumes that separately reconstructing the even and odd time-slice data and then combining the results cleanly removes the opposite-parity oscillating states, even though each reconstruction uses only half the time points and the paper itself warns that this may be statistically prohibitive.
Editorial extensions
If this is right
- If the even/odd separation works, every existing staggered-quark ensemble becomes usable for smeared spectral quantities, not just the two shown here.
- The same pipeline can be aimed at the hadronic tensor for inclusive pion scattering, the paper's stated next target.
- A reliable smeared $R$ ratio from staggered quarks would provide a lattice cross-check of muon $g-2$ hadronic vacuum polarization through the dispersion relation connecting $R(s)$ to the vector correlator.
- The observed low-energy agreement suggests that smearing with width roughly 300 to 500 MeV suppresses the inverse problem enough that discretization and finite-volume effects, though present, do not dominate below 1 GeV.
- Above 1.25 GeV the deviations define a targeted problem, quantifying lattice artifacts and finite-volume effects, that the authors plan to address in follow-up work.
Reading between the lines
- A testable extension is to vary the smearing width $\sigma$: if the parity-separation logic is correct, the definite-parity spectral densities should be $\sigma$-independent after smearing within errors, while any residual opposite-parity leakage would show up as a systematic drift.
- The even/odd trick likely transfers to nucleon inclusive structure functions relevant to neutrino experiments, since the same parity-mixing structure appears in staggered nucleon correlators; the paper does not make this claim.
- If halved time points prove prohibitive, the paper's interpolation and oscillating-state subtraction ideas could recover full statistics; comparing those three methods on the same ensemble would isolate the systematic error of the even/odd split.
- The comparison to a single domain-wall ensemble at different lattice spacing cannot separate quark-mass effects from discretization effects; a dedicated same-spacing comparison would sharpen the conclusion.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. This LATTICE2024 proceedings reports a pilot study of the smeared R ratio using staggered-quark correlation functions and the Hansen-Lupo-Tantalo (HLT) spectral reconstruction. The authors compute vector-current correlators on two MILC HISQ ensembles and one RBC/UKQCD domain-wall ensemble, and compare the reconstructed smeared spectral functions to the Bernecker-Meyer parameterization. To handle the opposite-parity (oscillating) contributions in staggered correlators, they separate even and odd time slices, reconstruct a smeared spectral function for each subset, and then combine them via Eq. (11) to isolate the J^P=1^- contribution. The paper reports 'promising agreement' with the parameterization below about 1 GeV and deviations at higher energies, and it discusses several alternative strategies (oscillating-state subtraction, correlator interpolation, and modified basis functions) for future work.
Significance. If the method is validated, this would be a valuable step toward computing inclusive hadronic observables such as the hadronic tensor from lattice QCD using staggered fermions, which are computationally cheaper and widely used. The paper is honest about its preliminary nature: it presents a single figure, explicitly lists uncontrolled systematics, and frames the result as 'progress' rather than a final determination. The explicit discussion of opposite-parity contamination strategies is useful. However, the central methodological claim — that the even/odd time-slice procedure yields the smeared 1^- spectral function — is not quantitatively established because the separate HLT reconstructions do not share a common smearing kernel. This is a load-bearing issue that must be addressed before the staggered curves in Fig. 1 can be interpreted as the desired smeared R ratio.
major comments (3)
- [Section 4.1, Eq. (11)]
- [Section 5, Fig. 1]
- [Section 5, Fig. 1]
minor comments (3)
- [Throughout]
- [Section 3, Table 1]
- [Section 4, Eq. (10)]
Circularity Check
No significant circularity: the smeared R-ratio analysis is benchmarked against an external parameterization, and the parity decomposition is an algebraic identity rather than a fitted input.
full rationale
The paper's central pipeline applies the independent HLT spectral-reconstruction method to staggered and domain-wall correlators and compares the resulting smeared R ratio to the Bernecker-Meyer parameterization [51]. No parameter of the reconstruction is fitted to that benchmark: the Gaussian smearing width sigma is chosen a priori, the HLT coefficients solve the convex optimization of Eqs. (7)-(8), and the quoted errors combine jackknife and HLT systematic errors following Ref. [34]. The even/odd parity separation in Eq. (11) is an exact algebraic consequence of the staggered spectral decomposition in Eq. (10), not a fitted quantity, and it is applied before comparison with the external parameterization. The paper explicitly labels its own result preliminary and flags unquantified artifacts, including opposite-parity effects, lattice artifacts, and finite-volume effects (Section 6). The citations to prior work by the same authors or collaborators (Refs. [29], [40], [49]) concern analysis perspective, the fat-link-only valence action, and interpolation in a related context; none supplies the paper's central claim or forbids alternatives. The skeptic concern about different even/odd HLT kernels is a possible systematic bias, not a circularity: it does not make the prediction equivalent to an input by construction. Overall the derivation is self-contained against external benchmarks.
Assumptions & free parameters
free parameters (3)
- HLT regularization parameter lambda =
not reported
- Gaussian smearing widths sigma =
530 and 280 MeV
- HLT basis parameters (tau_max, omega_0, alpha) =
not reported
assumptions (4)
- standard math The Euclidean correlator is the Laplace transform of a spectral density, and the smeared finite-volume spectral function converges to the infinite-volume one in the ordered limit sigma to 0 after L to infinity.
- domain assumption Staggered vector-current correlators admit the spectral decomposition of Eq. (10), with opposite-parity states appearing with factor (-1)^{tau+1}.
- domain assumption The HLT algorithm of Ref. [7], with systematic errors defined as in Ref. [34], produces a reliable smeared spectral reconstruction for the present data.
- domain assumption The Bernecker-Meyer parameterization (Ref. [51]) is an adequate benchmark for comparison at the quoted quark content and smearing widths.
Cite this review
Pith. "Pith review of Toward inclusive observables with staggered quarks: the smeared $R$~ratio." pith.science (2026). https://pith.science/paper/CA2TYXXV
@misc{pith2026241114300,
author = {Pith},
title = {Pith review of: Toward inclusive observables with staggered quarks: the smeared $R$~ratio},
year = {2026},
howpublished = {\url{https://pith.science/paper/CA2TYXXV}},
note = {Machine review of arXiv:2411.14300}
}
abstract
Inclusive hadronic observables are ubiquitous in particle and nuclear physics. Computation of these observables using lattice QCD is challenging due the presence of a difficult inverse problem. As a stepping stone to more complicated observables, we report on progress to compute the smeared $R$~ratio with staggered quarks using the spectral reconstruction algorithm of Hansen, Lupo, and Tantalo. We compare staggered-quark results on two ensembles to domain-wall results on a single ensemble and to the Bernecker-Meyer parameterization. This work utilizes two ensembles generated by the MILC collaboration using highly improved staggered quarks and one ensemble generated by the RBC/UKQCD collaboration using domain-wall quarks. Possible strategies for controlling opposite-parity effects associated with staggered quarks are discussed.
Figures
Reference graph
Works this paper leans on
-
[1]
S. Hashimoto,Inclusive semi-leptonic b meson decay structure functions from lattice qcd, Progress of Theoretical and Experimental Physics2017(2017) 053B03 [1703.01881]
arXiv 2017
-
[2]
P. Gambino and S. Hashimoto,Inclusive semileptonic decays from lattice qcd,Phys. Rev. Lett.125 (2020) 032001 [2005.13730]
arXiv 2020
-
[3]
P. Gambino, S. Hashimoto, S. Mächler, M. Panero, F. Sanfilippo, S. Simula et al.,Lattice qcd study of inclusive semileptonic decays of heavy mesons,JHEP07(2022) 083 [2203.11762]
arXiv 2022
- [4]
-
[5]
A.V. Manohar,An Introduction to spin dependent deep inelastic scattering, inLake Louise Winter Institute: Symmetry and Spin in the Standard Model, 3, 1992 [hep-ph/9204208]
arXiv 1992
-
[6]
Pijpers and M
F. Pijpers and M. Thompson,Faster formulations of the optimally localized averages method for helioseismic inversions, Astronomy and Astrophysics262 (1992) L33
1992
- [7]
-
[8]
H. Shinaoka, J. Otsuki, M. Ohzeki and K. Yoshimi,Compressing Green’s function using intermediate representation between imaginary-time and real-frequency domains,Phys. Rev. B 96(2017) 035147
work page 2017
Show all 52 references
-
[9]
Otsuki, M
J. Otsuki, M. Ohzeki, H. Shinaoka and K. Yoshimi,Sparse modeling approach to analytical continuation of imaginary-time quantum monte carlo data,Phys. Rev. E95(2017) 061302
2017
-
[10]
Itou and Y
E. Itou and Y. Nagai,Sparse modeling approach to obtaining the shear viscosity from smeared correlation functions, JHEP07(2020) 007 [2004.02426]
2020 arXiv
-
[11]
Otsuki, M
J. Otsuki, M. Ohzeki, H. Shinaoka and K. Yoshimi,Sparse modeling in quantum many-body problems, J. Phys. Soc. Jap.89(2020) 012001 [1911.04116]. 8 Toward inclusive observables with staggered quarks Douglas B. A. Stewart
2020 arXiv
-
[12]
J. Li, M. Wallerberger, N. Chikano, C.-N. Yeh, E. Gull and H. Shinaoka,Sparse sampling approach to efficient ab initio calculations at finite temperature, Phys. Rev. B101 (2020) 035144 [1908.07575]
2020 arXiv
-
[13]
Chikano, J
N. Chikano, J. Otsuki and H. Shinaoka,Performance analysis of a physically constructed orthogonal representation of imaginary-time Green’s function, Phys. Rev. B98(2018) 035104 [1803.07257]
2018 arXiv
-
[14]
S. Shi, L. Wang and K. Zhou,Rethinking the ill-posedness of the spectral function reconstruction: Why is it fundamentally hard and how artificial neural networks can help, Comput. Phys. Commun.282 (2023) 108547 [2201.02564]
2023 arXiv
-
[15]
Chen, H.T
S.Y. Chen, H.T. Ding, F.Y. Liu, G. Papp and C.B. Yang,Machine learning spectral functions in lattice QCD,PoSLATTICE2021(2021) 148 [2110.13521]
2021 arXiv
-
[16]
Kades, J.M
L. Kades, J.M. Pawlowski, A. Rothkopf, M. Scherzer, J.M. Urban, S.J. Wetzel et al.,Spectral reconstruction with deep neural networks,Phys. Rev. D102 (2020) 096001 [1905.04305]
2020 arXiv
-
[17]
K. Zhou, L. Wang, L.-G. Pang and S. Shi,Exploring QCD matter in extreme conditions with machine learning,Progress in Particle and Nuclear Physics135 (2023) 104084 [2303.15136]
2023 arXiv
-
[18]
Lechien and D
T. Lechien and D. Dudal,Neural network approach to reconstructing spectral functions and complex poles of confined particles, SciPost Phys.13(2022) 097 [2203.03293]
2022 arXiv
-
[19]
L. Wang, S. Shi and K. Zhou,Reconstructing spectral functions via automatic differentiation, Phys. Rev. D106 (2022) L051502 [2111.14760]
2022 arXiv
-
[20]
Valentine and M
A.P. Valentine and M. Sambridge,Gaussian process models I: A framework for probabilistic continuous inverse theory,Geophys. J. Int.220 (2020) 1632
2020
-
[21]
Del Debbio, T
L. Del Debbio, T. Giani and M. Wilson,Bayesian approach to inverse problems: an application to NNPDF closure testing,Eur. Phys. J. C82(2022) 330 [2111.05787]
2022 arXiv
-
[22]
Candido, L
A. Candido, L. Del Debbio, T. Giani and G. Petrillo,Inverse Problems in PDF determinations,PoSLATTICE2022(2023) 098 [2302.14731]
2023 arXiv
-
[23]
Horak, J.M
J. Horak, J.M. Pawlowski, J. Rodríguez-Quintero, J. Turnwald, J.M. Urban, N. Wink et al., Reconstructing QCD spectral functions with Gaussian processes,Phys. Rev. D105(2022) 036014 [2107.13464]
2022 arXiv
-
[24]
Pawlowski, C.S
J.M. Pawlowski, C.S. Schneider, J. Turnwald, J.M. Urban and N. Wink,Yang-Mills glueball masses from spectral reconstruction, Phys. Rev. D108 (2023) 076018 [2212.01113]
2023 arXiv
-
[25]
Horak, J.M
J. Horak, J.M. Pawlowski, J. Turnwald, J.M. Urban, N. Wink and S. Zafeiropoulos, Nonperturbative strong coupling at timelike momenta,Phys. Rev. D107(2023) 076019 [2301.07785]. 9 Toward inclusive observables with staggered quarks Douglas B. A. Stewart
2023 arXiv
-
[26]
Rothkopf,Inverse problems, real-time dynamics and lattice simulations, EPJ Web Conf
A. Rothkopf,Inverse problems, real-time dynamics and lattice simulations, EPJ Web Conf. 274 (2022) 01004 [2211.10680]
2022 arXiv
-
[27]
Rothkopf,Bayesian inference of real-time dynamics from lattice QCD,Front
A. Rothkopf,Bayesian inference of real-time dynamics from lattice QCD,Front. Phys.10 (2022) 1028995 [2208.13590]
2022 arXiv
-
[28]
Bulava,The spectral reconstruction of inclusive rates, PoSLATTICE2022(2023) 231 [2301.04072]
J. Bulava,The spectral reconstruction of inclusive rates, PoSLATTICE2022(2023) 231 [2301.04072]
2023 arXiv
-
[29]
Bergamaschi, W.I
T. Bergamaschi, W.I. Jay and P.R. Oare,Hadronic structure, conformal maps, and analytic continuation, Phys. Rev. D108 (2023) 074516 [2305.16190]
2023 arXiv
-
[30]
Bailas, S
G. Bailas, S. Hashimoto and T. Ishikawa,Reconstruction of smeared spectral functions from Euclidean correlation functions, Progress of Theoretical and Experimental Physics2020 (2020) 043B07
2020
-
[31]
Lehner and A.S
C. Lehner and A.S. Meyer,Consistency of hadronic vacuum polarization between lattice QCD and the R ratio,Phys. Rev. D101 (2020) 074515 [2003.04177]
2020 arXiv
-
[32]
Aoyama et al.,The anomalous magnetic moment of the muon in the Standard Model,Phys
T. Aoyama et al.,The anomalous magnetic moment of the muon in the Standard Model,Phys. Rept.887 (2020) 1 [2006.04822]
2020 arXiv
-
[33]
De Santis et al.,Lattice calculation of the𝑅 ratio smeared with Gaussian kernels,PoS LATTICE2022(2023) 307 [2212.12493]
A. De Santis et al.,Lattice calculation of the𝑅 ratio smeared with Gaussian kernels,PoS LATTICE2022(2023) 307 [2212.12493]
2023 arXiv
-
[34]
Extended Twisted Mass (ETM) collaboration, Probing the energy-smeared𝑅 ratio using lattice QCD, Phys. Rev. Lett.130 (2023) 241901 [2212.08467]
2023 arXiv
-
[35]
DeGrand and C.E
T. DeGrand and C.E. Detar,Lattice Methods for Quantum Chromodynamics, World Scientific, Singapore (2006)
2006
-
[36]
Huang, E
Z. Huang, E. Gull and L. Lin,Robust analytic continuation of Green’s functions via projection, pole estimation, and semidefinite relaxation, Phys. Rev. B107 (2023) 075151 [2210.04187]
2023 arXiv
-
[37]
Fei, C.-N
J. Fei, C.-N. Yeh and E. Gull,Nevanlinna analytical continuation, Phys. Rev. Lett.126 (2021) 056402 [2010.04572]
2021 arXiv
-
[38]
Hansen, H.B
M.T. Hansen, H.B. Meyer and D. Robaina,From deep inelastic scattering to heavy-flavor semileptonic decays: Total rates into multihadron final states from lattice QCD,Phys. Rev. D 96 (2017) 094513 [1704.08993]
2017 arXiv
-
[39]
MILCcollaboration, Lattice QCD ensembles with four flavors of highly improved staggered quarks,Phys. Rev. D87(2013) 054505 [1212.4768]
2013 arXiv
-
[40]
Aubin, T
C. Aubin, T. Blum, M. Golterman and S. Peris,Muon anomalous magnetic moment with staggered fermions: Is the lattice spacing small enough?, Phys. Rev. D106 (2022) 054503 [2204.12256]. 10 Toward inclusive observables with staggered quarks Douglas B. A. Stewart
2022 arXiv
-
[41]
RBC, UKQCD collaboration, Domain wall QCD with physical quark masses, Phys. Rev. D 93 (2016) 074505 [1411.7017]
2016 arXiv
-
[42]
RBC, UKQCD collaboration, Update of Euclidean windows of the hadronic vacuum polarization,Phys. Rev. D108 (2023) 054507 [2301.08696]
2023 arXiv
-
[43]
Foley, K
J. Foley, K. Jimmy Juge, A. O’Cais, M. Peardon, S.M. Ryan and J.-I. Skullerud,Practical all-to-all propagators for lattice QCD,Comput. Phys. Commun.172 (2005) 145 [hep-lat/0505023]
2005 arXiv
-
[44]
Boyle, A
P.A. Boyle, A. Yamaguchi, A. Portelli and G. Cossu,Grid: A next generation data parallel C++ QCD library, PoSLATTICE2015(2016) 023 [1512.03487]
2016 arXiv
-
[45]
Portelli, N
A. Portelli, N. Lachini, F. Erben, F. Joswig, R.R. Hodgson, F. Ó hÓgáin et al., aportelli/hadrons: Hadrons v1.4, June, 2023. 10.5281/zenodo.8023716
2023 doi
-
[46]
Lepage,gplepage/corrfitter: corrfitter version 8.1.1, Nov., 2020
G.P. Lepage,gplepage/corrfitter: corrfitter version 8.1.1, Nov., 2020. 10.5281/zenodo.4281296
2020 doi
-
[47]
Bennett et al.,Meson spectroscopy from spectral densities in lattice gauge theories, Phys
E. Bennett et al.,Meson spectroscopy from spectral densities in lattice gauge theories, Phys. Rev. D110 (2024) 074509 [2405.01388]
2024 arXiv
-
[48]
Lahert,The light-quark connected contribution to the muon’s anomalous magnetic moment, Ph.D
S. Lahert,The light-quark connected contribution to the muon’s anomalous magnetic moment, Ph.D. thesis, University of Illinois Urbana-Champaign, 2023
2023
-
[49]
Lahert, C
S. Lahert, C. DeTar, A.X. El-Khadra, S. Gottlieb, A.S. Kronfeld and R.S. Van de Water,The two-pion contribution to the hadronic vacuum polarization with staggered quarks, 2409.00756
-
[50]
Lepage, C
G.P. Lepage, C. Gohlke and D. Hackett,gplepage/gvar: gvar version 13.1.1, July, 2024. 10.5281/zenodo.12675777
2024 doi
-
[51]
Bernecker and H
D. Bernecker and H. Meyer,Vector correlators in lattice QCD: methods and applications, Eur. Phys. J. A47(2011) 148 [1107.4388]
2011 arXiv
-
[52]
Fermilab Lattice, HPQCD, MILC collaboration, Light-quark connected intermediate-window contributions to the muon𝑔− 2 hadronic vacuum polarization from lattice QCD, Phys. Rev. D107 (2023) 114514 [2301.08274]. 11
2023 arXiv
Reviewed August 12, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.