Reconstruction of the vector meson propagator using a generalized eigenvalue problem

Fabian Frech; Finn Stokes; Kalman Szabo

REVIEW 4 major objections 5 minor 24 references

Reconstruction of the vector meson propagator using a generalized eigenvalue problem

T0 review · 4 major / 5 minor · reviewed 2026-08-09 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗

Pith's one-line read By replacing nearly all sequential Dirac inversions with 1500 precomputed low eigenmodes plus stochastic residual sources, this paper removes the $T/a$ factor from the cost of connected two-pion propagators and reconstructs the…

desk verdict A useful low-mode-averaging trick for connected two-pion propagators, with a real but unquantified truncation bias that a referee should ask to be pinned down. read the letter →

arxiv 2501.19186 v1 pith:2SONR4MR submitted 2025-01-31 hep-lat

Fabian Frech , Finn Stokes , Kalman Szabo This is my paper · ORCID

classification hep-lat

keywords latticeQCDstaggeredfermionsrhomesontwo-pionstatesgeneralizedeigenvalueproblemlow-modeaveraginghadronicvacuumpolarizationmuong-2

verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

Long-distance measurements of the rho-meson correlator in lattice QCD lose signal exponentially, and the rho channel is also fed by ten two-pion states. The connected two-pion diagrams that feed this channel are expensive because they require a quark-propagation solve on every time slice, a cost proportional to the time extent $T/a$. This paper claims that splitting the quark propagator into 1500 low eigenmodes plus a stochastic residual removes the $T/a$ factor from the inversion count. It then feeds the cheap two-pion propagators into a generalized eigenvalue problem together with the rho propagators, resolves the contributing energy states, and reconstructs the long-distance rho signal with much smaller noise. The payoff is a cheaper route to the long-distance part of the hadronic vacuum polarization relevant to the muon $g-2$ anomaly.

What carries the argument

The load-bearing object is the low-mode-averaging decomposition of the staggered Dirac propagator in Eqs. (4)-(5): $M^{-1}=\sum_i \frac{1}{\lambda_i}|i\rangle\langle i|+M_r^{-1}$, with $M_r^{-1}$ estimated by stochastic sources $\sigma$. Its job is to replace the time-slice-dependent inversion in the connected two-pion diagram with quantities computed once, so the total inversion count loses its $T/a$ factor. The second mechanism is the generalized eigenvalue problem (GEVP), a variational method that extracts energy eigenstates from a matrix of correlators sharing quantum numbers: solving $C(t_0+dt)V=\lambda C(t_0)V$ on the 13-dimensional matrix of shifted $\rho$ and two-pion correlators separates the vector meson from the ten lower two-pion states. Together they convert a noisy late-time signal into a sum of well-separated exponentials.

What would settle it

Take one fixed two-pion channel (say Goldstone taste, $\vec p^2=2$) and compare the low-mode-averaging estimator of Eq. (2) using 1500 eigenmodes and a set number of residual sources against a full sequential inversion on the same configurations at all $t$ out to $T/2$. If the ratio deviates from unity by more than statistics, or if the deviation grows with $t$ or with $|\vec p|^2$, the truncated eigenmode-plus-stochastic sum is biased for that channel; repeating the GEVP reconstruction with different eigenmode counts (for example 1000 versus 2000) and source numbers would then show whether the extracted $M_i$ and $a_i$ shift, which would settle the central claim.

Watch

Extended reading notes

Core claim

The central claim is that low-mode averaging makes the connected $\pi\pi\to\pi\pi$ part of the vector correlator cheap enough for a state-by-state reconstruction. Writing $M^{-1}=\sum_i \lambda_i^{-1}|i\rangle\langle i|+M_r^{-1}$ and approximating $M_r^{-1}$ by stochastic sources removes the time-dependent inversion in Eq. (2), so the number of Dirac inversions no longer grows with the lattice time extent $T/a$. With those two-pion propagators in hand, the paper assembles the correlation matrix of Eq. (7) -- the $\rho\rho$, $\rho\pi\pi$, $\pi\pi\rho$ and $\pi\pi\pi\pi$ blocks, extended by pencil-of-functions shifts to a 13-dimensional problem -- and solves the GEVP $C(t_0+dt)V=\lambda C(t_0)V$ at $t_0=5a$, $dt=a$. The resulting eigenvalues $\tilde\lambda_i(t)$ give effective energies and the eigenvectors give coefficients $a_i$; the reconstructed $\rho_{\rm rec}(t)=\sum_i a_i e^{-M_i t}$ matches the directly simulated rho propagator in the shown range while carrying less noise at large $t$. The authors frame this as the first step toward more precise staggered-fermion calculations of the hadronic vacuum polarization relevant to the muon $g-2$ discrepancy.

Load-bearing premise

The load-bearing premise is that 1500 low eigenmodes plus the stochastic residual sources reproduce the exact quark propagator accurately for every pion operator, momentum, and taste that enters the correlator, with no bias that grows with time or momentum; the paper shows the comparison for only one momentum/taste channel and does not include a convergence study, so if the approximation drifts, the GEVP inputs and the reconstructed rho propagator inherit a systematic error.

Editorial extensions

If this is right

The cost of connected two-pion diagrams becomes independent of the lattice time extent, so longer boxes and longer Euclidean time separations are accessible without multiplying the Dirac inversion count.
The rho propagator can be decomposed into its physical constituents: the vector meson and the two-pion states below it, each with a fitted energy and coefficient.
The reconstructed long-distance rho propagator carries much smaller statistical noise than direct simulation, which is precisely the noisy part of the hadronic vacuum polarization that matters for the muon $g-2$.
Because the 13-dimensional GEVP filters out two noisy high-excitation states with the pencil-of-functions shifts, early-time data alone are sufficient to determine the late-time behavior.
The same pipeline with smeared vector operators -- the authors' own planned next step -- should improve the overlap of the $\rho$ and $\pi\pi$ interpolators.

Reading between the lines

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

A direct corollary the paper does not pursue is that the same low-mode-averaging trick should apply to other four-point functions whose sequential inversion runs over time slices, not just the $\rho$ channel.
The observed momentum and taste dependence of the noise penalty suggests an optimization the paper does not spell out: use more eigenmodes for high-momentum or nonlocal-taste operators and more stochastic sources for low-momentum ones, or choose the split per channel.
If the reconstruction remains stable when the eigenmode count is varied, the method could become a practical reconstruction-first strategy for the long-distance window of $a_\mu^{\rm HVP}$, reducing statistical effort on the expensive late-time region.
The GEVP decomposition is a lattice analogue of separating $\rho$ and nonresonant $\pi\pi$ contributions in experiment; with enough states it could produce the timelike pion form factor directly from lattice data, though the paper does not attempt that.

Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Desk Editor's Note

read the letter

Colleague—

The one thing to know: this is a short proceedings paper whose real contribution is a low-mode averaging (LMA) trick for the connected two-pion propagator that removes the T/a factor from the inversion count. That is useful if you work on staggered HVP. The GEVP reconstruction of the rho is carried over from their earlier work, and the reconstruction itself is a demonstration rather than a new method.

What is solid: the scaling argument in Fig. 2 is concrete and the asymptotic advantage is real; the idea of using precomputed eigenvectors from the Dirac preconditioner plus stochastic residual sources is sensible and cheap to test. Fig. 3 shows the method reproduces the full-inversion propagator in one channel, and the text is honest that the comparison is representative, with larger momenta and non-local tastes showing more difference.

The soft spots are the ones the stress-test flags, and they are genuine. The single validation channel is a Goldstone taste at |p|^2=2. The GEVP input matrix includes momenta up to |p|^2=4 and multiple tastes (Table 1), and the text concedes those are exactly the channels where LMA and full inversion differ more. There is no convergence study over the 1500 eigenvectors or the T/2 stochastic sources, and the GEVP plots have no error bars. The concern is not that the method is wrong; it is that the size of the bias in the reconstructed rho propagator is unknown. Using early-time slices in the GEVP does not fix a bias already present in C(t0) and C(t0+dt). This is an addressable gap, not a fatal flaw. Also, this is a proceedings paper: no ensemble details, no analysis code, no jackknife or bootstrap description. For a proceedings that is acceptable, but it limits how much independent checking a referee can do.

The circularity point is mild. The GEVP states are extracted from the same correlator matrix used for the reconstruction, but the masses and amplitudes are fit to early-time plateaus rather than tuned to match the late-time propagator, so it is not a fitted tautology. Still, a fully independent check would need a second ensemble or a different method.

Who gets value: lattice people working on the two-pion contribution to HVP, especially in the BMW-style staggered program. It is not a paper that changes the broader g-2 picture on its own. I would send it to a referee rather than desk reject, expecting that referee to ask for a convergence study and error bars on the GEVP quantities. The core efficiency claim is plausible and worth engaging.

Referee Report

4 major / 5 minor

Summary. The paper proposes an efficient way to compute the connected two-pion contribution to the staggered vector-meson correlator by splitting the Dirac inverse into 1500 low eigenmodes plus a stochastic residual estimate (Eqs. (4)-(5)), thereby reducing the inversion count relative to the sequential-inversion formula in Eq. (3). The resulting two-pion and rho correlators are assembled into a 13-dimensional correlation matrix, and a generalized eigenvalue problem (GEVP, Eqs. (6)-(8)) is used to extract energy states and reconstruct the rho propagator. Results are presented for a single (6^3 x 8.5) fm^4 ensemble, including cost scaling, a one-channel comparison of low-mode averaging with full inversion, and a reconstructed-versus-simulated rho correlator.

Significance. If the low-mode averaging estimator is unbiased in all channels entering the GEVP, the method would be a practically useful reduction of the cost of connected two-pion diagrams in staggered hadronic vacuum polarization calculations, and the GEVP reconstruction provides a physically instructive decomposition of the rho propagator. The paper is transparent about the construction, gives explicit Clebsch-Gordan tables and cost formulas, and states the approximation in Eq. (5) clearly. The main significance, however, depends on validation that is currently missing: the truncation at 1500 eigenmodes is checked only in one channel, and no statistical errors are shown for the GEVP results that underpin the central reconstruction claim.

major comments (4)

[§3, Eq. (5), Fig. 3] The estimator in Eq. (5) is the basis for the efficiency gain, but its accuracy is established only by the single Goldstone-taste, |p|^2=2 channel shown in Fig. 3. The ratio panel has no error bars and visibly departs from unity at large t, and the text itself states that larger momenta and non-local tastes show larger differences. Those channels are exactly the ones entering the 13x13 GEVP (Table 1 includes momenta up to |p|^2=4 and all taste structures). A bias in C(t0) or C(t0+dt) from the eigenmode truncation is inherited by the GEVP eigenvalues and amplitudes in Eqs. (6)-(8), so the central reconstruction claim requires a convergence study in the eigenvector count and the number of residual sources, plus a channel-by-channel LMA-versus-full-inversion comparison at least for the states used in the GEVP.
[Abstract and §3, Eq. (3), Fig. 2] The abstract and §3 state that the method removes the T/a factor from the inversion count in Eq. (3). The implementation described in the caption of Fig. 2, however, uses T/2 stochastic residual sources sigma for the residual part of Eq. (5), which introduces an O(T) cost unless it is explicitly shown that these sources are reused across all time separations and operators without further inversions. As written, the asymptotic statement should be that the method replaces N_O x T/a with an O(T) term (with a better prefactor), not that the T/a factor is removed. The cost claim in its current form is therefore overstated.
[§4, Figs. 4-5] The central numerical claim, that the reconstructed rho propagator has much smaller noise at large Euclidean times, is not supported by a statistical analysis. Figures 4 and 5 show no error bars on the GEVP eigenvalues, effective energies, effective coefficients, or the reconstructed-versus-simulated comparison, and no bootstrap, jackknife, or covariance treatment is described. Without errors, the claim of smaller noise and the reliability of the reconstruction cannot be assessed. The paper should provide error bands or error bars on M_eff, a_eff, and the late-time ratio of reconstructed to simulated correlators.
[§4, Eqs. (6)-(8), Fig. 5] The reconstruction is an in-sample check: the eigenvectors V_i are obtained from the same correlation matrix C(t) whose rho-column is later compared with rho_rec(t), and the masses and amplitudes are fitted to early-time plateaus of that same matrix. The agreement in Fig. 5 therefore does not by itself validate the decomposition into physical energy states. To support the physical claim, the authors should test stability under changes of t0, dt, the number of states kept, and the state-selection cutoff, and ideally compare with an independent full-inversion rho correlator or a different source setup.

minor comments (5)

[§3, Eq. (5)] The placement of M_r^{-1} relative to the stochastic sources in Eq. (5) is ambiguous; the text should state explicitly that the stochastic sources approximate the identity in the residual subspace and how the operator product acts on a source vector.
[§2, Table 1 and Table 2] There are several typos: 'Multiplicities of the the diagram', 'construcing', and 'approxiamtion' should be corrected.
[Fig. 2 caption] The caption says 'On the right-hand side' for both panels; one of these should refer to the left-hand panel.
[Fig. 3 caption] The caption says the full propagators are on the right-hand side and the ratio on the left-hand side, but the figure arrangement appears to be the opposite; the caption should match the panels.
[§4, Eq. (9)] The fit ranges, the choice of Delta=2a, and the criterion for filtering two noisy states with the pencil-of-functions method are not stated; these details are needed to reproduce the effective-energy and coefficient plateaus.

Circularity Check

0 steps flagged · score 2.0 of 10

No significant circularity; minor self-citation is non-load-bearing and independently confirmed.

full rationale

The central derivation is self-contained. The low-mode-averaging estimator in Eqs. (4)-(5) is an identity plus a stochastic approximation: M^-1 is split into an exact low-eigenmode sum and a residual estimated with random sources, giving an unbiased estimator. The cost-scaling claim (removal of the T/a factor from Eq. (3)) follows directly from the form of Eq. (2) and the estimator, not from any fitted parameter. The GEVP reconstruction in Eqs. (6)-(8) is a conventional spectral decomposition of the correlation matrix C(t). The reconstructed rho propagator is obtained by fitting effective energies and coefficients to early-time plateaus, then compared with the simulated C_rho_rho at all t in Fig. 5; because the fit uses only early-time data, the late-time comparison is a consistency check, not a quantity forced to match by construction. The only self-citation is to the authors' prior work [13] for staggered Clebsch-Gordan coefficients and two-pion operator construction, and the paper explicitly cites independent confirmations [7,22]; the LMA efficiency claim does not depend on that construction. Incomplete validation of LMA in all GEVP channels is a correctness risk, not circularity.

Assumptions & free parameters 5 free parameters · 4 assumptions · 0 invented entities

The central claim rests on the standard lattice QCD framework plus the low-mode averaging approximation. All nontrivial symmetry input is imported from [13]. The key unquantified choice is the 1500-eigenmode truncation and stochastic residual, which is an ad hoc approximation for this calculation.

free parameters (5)

Number of low eigenmodes = 1500
Chosen for the low-mode averaging split in Eq. (5); no convergence scan is shown, so the truncation error is uncontrolled.
GEVP reference times t0 and dt = t0=5a, dt=a
Fixed by hand in Eq. (6); results for effective energies and amplitudes depend on this choice and no stability check is shown.
Effective-mass shift Delta = 2a
Used in Eq. (9); chosen without justification or sensitivity test.
Stochastic source counts = 18 wall sources, 32 residual sources
The number of random sources affects the residual variance in Eq. (5); no noise-bias tradeoff study is shown.
State-selection cutoff = 2(p_i^2 + m_pi^2) <= m_rho^2
Defines which 10 two-pion states enter the generalized eigenvalue problem; omitted states are assumed negligible but no systematic check of this cutoff is reported.

assumptions (4)

domain assumption Staggered fermion symmetry group and Clebsch-Gordan coefficients from [13]
Used to construct two-pion states with vector quantum numbers in Section 2 and Tables 1-2; correctness of these coefficients is cited rather than derived.
domain assumption Euclidean lattice QCD and time-reversibility of QCD
Eq. (1) interchanges rho and two-pion states using time reversal to avoid separate inversions; this is a standard QCD property.
standard math GEVP spectral decomposition converges for the 13x13 matrix
Eq. (6) assumes the eigenvalues of the generalized problem isolate the physical states with exponentially suppressed contamination; standard variational method assumption.
ad hoc to paper Low-mode averaging identity approximation
Eq. (5) approximates the identity with stochastic sources sigma; the paper gives no proof of convergence or bias bound for this approximation.

how reviews work

0 comments

Cite this review

Pith. "Pith review of Reconstruction of the vector meson propagator using a generalized eigenvalue problem." pith.science (2026). https://pith.science/paper/2SONR4MR

@misc{pith2026250119186,
  author       = {Pith},
  title        = {Pith review of: Reconstruction of the vector meson propagator using a generalized eigenvalue problem},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/2SONR4MR}},
  note         = {Machine review of arXiv:2501.19186}
}

abstract

For long distances in the euclidean time the vector-vector correlator ($\rho$) has an exponentially decreasing signal-to-noise ratio. However, the vector correlator not only consists of the vector meson but also receives contributions from a two-pion system with the same quantum numbers. We measure all two-pion propagators with an energy lower than the mass of the resting vector meson and employ a generalized eigenvalue problem (GEVP) to resolve the different contributing energy states. Using those we can reconstruct the propagator with a much smaller noise at large euclidean time distances. In this work we present an efficient way to measure two-pion propagators and our results on reconstruction of the vector meson propagator with staggered fermions in a $(6^3\times 8.5)\,\textrm{fm}^4$ box.

Figures

Figures reproduced from arXiv: 2501.19186 by the authors.

**Figure 1.** The diagrams that have to be implemented for GEVP. Diagram c) vanishes identically in the 𝐼 = 1 case. 3. Computation of the connected 𝜋𝜋 → 𝜋𝜋 diagram The general strategy to compute the diagrams that are shown in [PITH_FULL_IMAGE:figures/full_fig_p003_1.png] view at source ↗

**Figure 2.** On the right-hand side, one can see the volume scaling of the LMA method compared to the usual sequential inversion. On the right-hand side, one can see the scaling with the number of operators. Both plots use 1500 eigenvectors and 18 random wall sources 𝜁. To estimate the residual part, we include 𝑇/2 random sources 𝜎. 0 20 40 60 t[a] 103 104 105 106 107 108 109 full inversion 1500 EV + random sources 0 20 40 60 t[… view at source ↗

**Figure 3.** Two Goldstone pions with momenta 𝑝®1 = (0, 1, 1) and 𝑝®2 = (0, −1, −1). The full propagators are on the right-hand side, while their ratio is shown on the left-hand side. We use 18 random wall sources for both methods. For the LMA method, we use 1500 eigenvectors and 32 sources to estimate the residual part. 5 [PITH_FULL_IMAGE:figures/full_fig_p005_3.png] view at source ↗

Figures from the paper (2 more)

**Figure 4.** Figure 4: Right hand side: Eigenvalues 𝜆˜(𝑡) of the GEVP including all the 𝜋𝜋-propagators with 2(𝑝 2 𝑖 + 𝑚 2 𝜋𝜉 ) ≤ 𝑚 2 𝜌 as well as the (shifted) 𝜌-propagators. RIght hand side: Effective energies 𝑀 eff 𝑖 (𝑡) of the states shown on the left hand side. 0 5 10 15 20 25 30 t/a 10−…

**Figure 5.** Figure 5: Left hand side: Coefficents 𝑎 eff. 𝑖 (𝑡) of the considered states. Right hand: Reconstruction of the vector meson propagator 𝜌rec.(𝑡)compared to the simulated propagator 𝜌(𝑡). 5. Conclusion In this work, we presented a refined approach to computing the connected 𝜋𝜋 → 𝜋…

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

24 extracted references · 13 canonical work pages

[1]

Aguillard, D. P. et al. (Muon g-2),Measurement of the Positive Muon Anomalous Magnetic Moment to 0.20 ppm, Phys. Rev. Lett. (2023), 161802. arXiv:2308.06230 [hep-ex]

arXiv 2023
[2]

Aoyama,T.etal., TheanomalousmagneticmomentofthemuonintheStandardModel ,Phys. Rept. (2020), 1. arXiv:2006.04822 [hep-ph]

arXiv 2020
[3]

arXiv:2002.12347 [hep-lat]

Borsanyi,Sz.etal., Leadinghadroniccontributiontothemuonmagneticmomentfromlattice QCD, Nature (2021), 51. arXiv:2002.12347 [hep-lat]

arXiv 2021
[4]

et al.,High precision calculation of the hadronic vacuum polarisation contri- bution to the muon anomaly, arXiv:2407.10913 [hep-lat], (2024)

Boccaletti, A. et al.,High precision calculation of the hadronic vacuum polarisation contri- bution to the muon anomaly, arXiv:2407.10913 [hep-lat], (2024)

arXiv 2024
[5]

Djukanovic, Dalibor et al.,The hadronic vacuum polarization contribution to the muon𝑔− 2 at long distances, arXiv:2411.07969 [hep-lat], CERN-TH-2024-196, MITP-24-080, (2024)

arXiv 2024
[6]

et al., Checks on QED and strong-isospin breaking corrections to𝑎HVP 𝜇 , PoS LATTICE2024 (2024), 271, (2025)

Risch, A. et al., Checks on QED and strong-isospin breaking corrections to𝑎HVP 𝜇 , PoS LATTICE2024 (2024), 271, (2025)

work page 2024
[7]

Lahert, Shaun et al.,The two-pion contribution to the hadronic vacuum polarization with staggered quarks, arXiv:2409.00756 [hep-lat], (2024)

arXiv 2024
[8]

Bazavov, Alexei et al., Hadronic vacuum polarization for the muon𝑔− 2 from lattice QCD: Long-distance and full light-quark connected contribution, arXiv:2412.18491 [hep- lat], (2024)

arXiv 2024

Show all 24 references

[9]

Erben, Felix et al.,Rho resonance, timelike pion form factor and implications for lattice studiesofthehadronicvacuumpolarization ,Phys.Rev.D(2020),054504.arXiv:1910.01083 [hep-lat]

2020 arXiv
[10]

Bruno, Mattia et al.,Exclusive Channel Study of the Muon HVP, PoS LATTICE2019 (2019),

2019
[11]

et al.,The long-distance window of the hadronic vacuum polarization for the muon g-2, arXiv:2410.20590 [hep-lat], (2024)

Blum, T. et al.,The long-distance window of the hadronic vacuum polarization for the muon g-2, arXiv:2410.20590 [hep-lat], (2024)

2024 arXiv
[12]

arXiv:2112.11647 [hep-lat]

Lahert, Shaun et al., Hadronic vacuum polarization of the muon on 2+1+1-flavor HISQ ensembles: an update, PoS LATTICE2021 (2022), 526. arXiv:2112.11647 [hep-lat]

2022 arXiv
[13]

Frech, Fabian J. et al. (Budapest-Marseille-Wuppertal),The mixing of two-pion and vector- meson states using staggered fermions, PoS LATTICE2023 (2024), 251. arXiv:2401.00514 [hep-lat]

2024 arXiv
[14]

and Hoelbling, C.,Locality and symmetry properties of staggered fermions with taste splitting mass term, PoS LATTICE2024 (2024), 370, (2024)

Chreim, N. and Hoelbling, C.,Locality and symmetry properties of staggered fermions with taste splitting mass term, PoS LATTICE2024 (2024), 370, (2024). 8 Reconstruction of the vector meson propagator using a generalized eigenvalue problemFabian J. Frech

2024
[15]

Golterman, Maarten,Staggered fermions, arXiv:2406.02906 [hep-lat], (2024)

2024 arXiv
[16]

Kilcup, G. W. et al.,A Tool Kit for Staggered Fermions, Nucl. Phys. B (1987), 493

1987
[17]

Golterman, Maarten F. L. et al.,Lattice Baryons With Staggered Fermions, Nucl. Phys. B (1985), 328

1985
[18]

Baake,M.etal., StructureandRepresentationsoftheSymmetryGroupoftheFour-dimensional Cube, J. Math. Phys. (1982), 944. [Erratum: J. Math. Phys. 23, 2595 (1982)]

1982
[19]

et al.,Representations of the Rotation Reflection Symmetry Group of the Four-dimensional Cubic Lattice, Nucl

Mandula, Jeffrey E. et al.,Representations of the Rotation Reflection Symmetry Group of the Four-dimensional Cubic Lattice, Nucl. Phys. B (1983), 91

1983
[20]

Golterman, Maarten F. L. et al.,Irreducible Representations of the Staggered Fermion Sym- metry Group, Nucl. Phys. B (1986), 417

1986
[21]

Sakata, I.,A general method for obtaining clebsch-gordan coefficients of finite groups. 1. its application to point and space groups, J. Math. Phys. (1974), 1702

1974
[22]

Lahert, Shaun, The light-quark connected contribution to the muon’s anomalous magnetic moment, PhD thesis, Illinois U., Urbana (main), Illinois U., Urbana, 2023

2023
[23]

DeTar, Carleton and Lee, Song-Haeng,Variational method with staggered fermions, Phys. Rev. D (2015), 034 9

2015
[239]

arXiv:1910.11745 [hep-lat]

1910 arXiv

Pith tools

Reviewed August 9, 2026 · model on record in the stance chip above.

[1] [1]

Aguillard, D. P. et al. (Muon g-2),Measurement of the Positive Muon Anomalous Magnetic Moment to 0.20 ppm, Phys. Rev. Lett. (2023), 161802. arXiv:2308.06230 [hep-ex]

arXiv 2023

[2] [2]

Aoyama,T.etal., TheanomalousmagneticmomentofthemuonintheStandardModel ,Phys. Rept. (2020), 1. arXiv:2006.04822 [hep-ph]

arXiv 2020

[3] [3]

arXiv:2002.12347 [hep-lat]

Borsanyi,Sz.etal., Leadinghadroniccontributiontothemuonmagneticmomentfromlattice QCD, Nature (2021), 51. arXiv:2002.12347 [hep-lat]

arXiv 2021

[4] [4]

et al.,High precision calculation of the hadronic vacuum polarisation contri- bution to the muon anomaly, arXiv:2407.10913 [hep-lat], (2024)

Boccaletti, A. et al.,High precision calculation of the hadronic vacuum polarisation contri- bution to the muon anomaly, arXiv:2407.10913 [hep-lat], (2024)

arXiv 2024

[5] [5]

Djukanovic, Dalibor et al.,The hadronic vacuum polarization contribution to the muon𝑔− 2 at long distances, arXiv:2411.07969 [hep-lat], CERN-TH-2024-196, MITP-24-080, (2024)

arXiv 2024

[6] [6]

et al., Checks on QED and strong-isospin breaking corrections to𝑎HVP 𝜇 , PoS LATTICE2024 (2024), 271, (2025)

Risch, A. et al., Checks on QED and strong-isospin breaking corrections to𝑎HVP 𝜇 , PoS LATTICE2024 (2024), 271, (2025)

work page 2024

[7] [7]

Lahert, Shaun et al.,The two-pion contribution to the hadronic vacuum polarization with staggered quarks, arXiv:2409.00756 [hep-lat], (2024)

arXiv 2024

[8] [8]

Bazavov, Alexei et al., Hadronic vacuum polarization for the muon𝑔− 2 from lattice QCD: Long-distance and full light-quark connected contribution, arXiv:2412.18491 [hep- lat], (2024)

arXiv 2024