REVIEW 3 major objections 5 minor 2 cited by
Long-distance reconstruction of QED corrections to the hadronic vacuum polarization for the muon g-2
T0 review · 3 major / 5 minor · reviewed 2026-08-05 · deepseek-v4-flash
Pith's one-line read The long-distance QED correction to the hadronic vacuum polarization can be reconstructed from the rho and pion-photon finite-volume states, cutting statistical noise on the hardest diagram by more than a factor of five.
desk verdict The isospin relations are new and the method is promising, but the individual-diagram reconstruction leans on an untested D1-dominance assumption that the noise-reduction claim depends on. 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 mechanism is a two-state reconstruction using a correlation matrix built from the vector current operator and a pion-photon operator. The pion-photon operator is formed by combining a charged pion with a transversely polarized photon in a definite octahedral-group channel, so that it can be used in a GEVP-like fit. The load-bearing identities are the diagram relations of Eq. (40), derived by comparing correlators with different isospin and quark-charge assignments; they relate the pion-photon contributions of diagrams (F), (V), (S), and (D1) and imply a strong cancellation in the physical combination. These relations let the authors reconstruct individual noisy diagrams from fitt
What would settle it
Compute the pion-photon contribution of diagram (D1) on the same ensemble; if it is comparable in size to the combination (V) + 2(S), the reconstruction is incomplete. Alternatively, perform a direct, high-statistics inclusive calculation of the long-distance G^(2)(t) integrand at physical pion mass and compare it with the reconstructed rho-plus-pion-photon sum.
Extended reading notes
Core claim
The paper claims that, at the pion mass studied, the long-distance part of the second-order QED contribution to the vector-vector correlator is dominated by two finite-volume states: the rho and the lowest pion-photon state. The authors construct a two-operator correlation matrix containing the physical vector current and a pion-photon operator, fit the energies and overlaps of both states, and use isospin and charge-assignment relations to connect the pion-photon contributions of diagrams (V), (S), (F), and (D1). These relations, summarized in Eq. (40), allow the pion-photon part of the noisy diagram (S) to be reconstructed from the cleaner pion-photon operator data. The numerical demonstra
Load-bearing premise
The numerical reconstruction assumes that diagram (D1) is negligible; if it is not, the relations of Eq. (40) do not by themselves determine the individual pion-photon pieces, and the agreement shown in the figures is conditional on that assumption.
Editorial extensions
If this is right
- The long-distance tail of e^2 G^(2)(t) can be computed without relying on model-based replacements or truncations, improving the first-principles status of the QED correction to the hadronic vacuum polarization.
- The pion-photon contribution can be separated from the rest of the QED corrections, making it possible to subtract the finite-volume QED_L pion-photon part and add back an infinite-volume QED calculation.
- The noise reduction, demonstrated as greater than a factor of five on diagram (S), directly translates into a reduction in the computational cost or uncertainty of the final g-2 contribution.
- At physical pion mass, the method extends by adding two-pion operators and states, so the approach is not limited to the heavy-pion demonstration presented here.
- The diagram relations could also be used to construct correlator combinations in which the pion-photon contribution is absent, an alternative to reconstruction that may be useful in other QCD+QED observables.
Reading between the lines
- The same reconstruction strategy, with appropriate operators, is likely applicable to the QED corrections of other long-distance-dominated observables, such as hadronic light-by-light pieces, where pion-photon states play a similar role.
- A testable extension is to apply the method at physical pion mass with multi-pion operators; if the reconstructed and direct inclusive long-distance integrands agree there, the method's practical validity would be established beyond the single heavy-pion ensemble.
- The strong suppression of the pion-photon contribution in the physical combination suggests that model-based estimates of this long-distance piece, currently used in some calculations, may be less numerically important than previously thought—though the size of the neglected sub-leading diagrams still needs to be checked.
- If diagram (D1) is found to be non-negligible, the relations of Eq. (40) would need to be extended with additional independent measurements; the current demonstration would then serve as a proof of the reconstruction machinery rather than a complete numerical result.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper proposes a method to reconstruct the long-distance part of the O(e^2) QED corrections to the hadronic vacuum polarization (HVP) from exclusive finite-volume states. Working in QCD+QED with a two-state truncation (rho and pi-gamma), the authors derive isospin/charge relations, Eq. (40), linking the pi-gamma contributions of diagrams (V), (S), (D1), and (F). They extract the pi-gamma energy and overlap factors from off-diagonal correlators and, after subtracting the pi-gamma contribution from the vector-vector correlator, fit the QED corrections to the rho parameters. Numerical results on one RBC/UKQCD ensemble at a^{-1}≈1.73 GeV and m_pi≈275 MeV are shown for the integrand t^4 G^(2)(t), including a claimed statistical noise reduction by more than a factor of five for diagram (S).
Significance. If the approach is viable, it addresses a known bottleneck in lattice HVP calculations: the long-distance QED corrections are currently modeled or truncated rather than computed from first principles. The isospin decomposition leading to Eq. (40) is explicit and checkable, and the numerical demonstration is internally consistent, with fitted effective masses, overlaps, and subtracted fits displayed. The main idea of reconstructing the pi-gamma and rho contributions to individual diagrams is attractive and potentially impactful for the g-2 precision program.
major comments (3)
- [Sec. III, Eq. (40)] The per-diagram pi-gamma reconstruction sets (D1)=0 through the assumption stated as 'assuming dominance over diagram (D1) for now.' Because Eq. (40) gives 2(S)-(D1)=(V), this forces (S)=V/2. The coefficient of (D1) in Eq. (3) is 25/162, comparable to the coefficients of (V) and (S), so the assumption is not a priori safe. Consequently the 'pi-gamma reconstruction' curves in Figs. 4 and 6 for (S), and hence the factor-of-five noise-reduction claim for diagram (S) in Fig. 7, are conditional on an untested assumption. The paper would need a numerical estimate or bound for the (D1) pi-gamma amplitude, or an independent means of fixing the split, before the individual-diagram reconstruction can be regarded as validated.
- [Sec. II A, Eq. (14)] The overlap |V_{1,pi-gamma}^{(1)}| is extracted from the same two-state effective model used to form the reconstructed pi-gamma contribution to G^(2). Thus the agreement in Figs. 4 and 6 is a self-consistency check of the model together with Eq. (40) under the D1 assumption, not an independent prediction of the vector-vector correlator. The paper should state this clearly and provide a direct test of the two-state ansatz, for example by comparing the reconstructed inclusive G^(2) against the raw correlator in the mid-time window or by adding a third operator to the GEVP.
- [Sec. III, Figs. 4-7] The numerical demonstration is at a single lattice spacing, a single unphysical pion mass (m_pi≈275 MeV), and m_pi L=3.8, with no systematic assessment of excited-state contamination or finite-volume QED effects. The fits include excited states, but no stability analysis or alternate fit ranges are shown, and the claimed noise reduction is not quantified with statistical errors before and after reconstruction. These limitations are acknowledged in the text as future work, but they should be reflected in the abstract and conclusion, where the strength of 'high-precision' and 'more than a factor of five' currently exceeds what one ensemble can establish.
minor comments (5)
- [Sec. I] Typo: 'correctons' should be 'corrections'; also 'Mobius' should be 'Möbius' if a unified spelling is desired.
- [Figs. 4 and 6] The legend labels are confusing: some panels are labelled '(V)' and others '(L)', and the meaning of the open symbols versus lines is not fully defined. Please use consistent labels and clarify which curves correspond to the direct lattice data and which to the reconstruction.
- [Sec. III, Fig. 7] The 'more than a factor of five' noise reduction is stated without a precise statistical comparison. A small table with the errors of the partial sums with and without reconstruction at representative t values would make the claim quantitative and reproducible.
- [Sec. II C] The third relation in Eq. (40) involving (T), (Td), (D1d), (D3), (D2d), and (D2) is derived but not used in the reconstruction. A short comment on why it is not exploited and whether it could provide a cross-check would be helpful.
- [Sec. III] The statement that the comparison in Fig. 4 is incomplete without vector-current renormalization is made but not elaborated. Since the pi-gamma contribution is unaffected by Z_V^QED, it would be useful to state explicitly which reconstructed quantities do and do not depend on the missing renormalization.
Circularity Check
No circularity: the diagram relations are derived from first principles and the numerical reconstruction is a transparent fit-based method, not a prediction from its own inputs.
full rationale
The central derivation (Sec. II C, Eq. (40)) obtains the relations (F)=(V) and 2(S)-(D1)=(V) by comparing isospin/charge decompositions of vector-vector and pion-photon-vector correlators. This is a first-principles result; it does not assume the target reconstruction. The numerical reconstruction uses |V^(1)_{1,pi-gamma}| fitted from the off-diagonal correlator C21 (Sec. II A/III), not from the individual diagram data it is compared against in Figs. 4 and 6. The comparison is therefore a genuine cross-check of the relations and the D1≈0 assumption, not a fit to the same data relabeled as a prediction. The paper explicitly labels the D1 neglect as an assumption ('assuming dominance over diagram (D1) for now') and lists it as future work; this is a limitation, not a circular step. The leading-order rho parameters are taken from a separate published lattice study [8]; this is external support, not a self-citation chain. No fitted parameter is renamed as a prediction, and no uniqueness theorem is imported from the authors' prior work. The method is self-contained and the numerical demonstration is honestly presented as a reconstruction.
Assumptions & free parameters
free parameters (6)
- E^(0)_(pi-gamma) = E_pi + E_gamma =
fitted from C22 effective masses (Fig. 2)
- |V^(1)_(2,pi-gamma)| =
fitted from C22 (Fig. 2)
- |V^(1)_(1,pi-gamma)| =
fitted from C21 (Fig. 3)
- c^(2)_rho = 2Re(V^(0)_(1,rho) V^(2)*_(1,rho)) =
fit to C11 after pion-photon subtraction (Fig. 5)
- c^(0)_rho E^(2)_rho = |V^(0)_(1,rho)|^2 E^(2)_rho =
fit to C11 after pion-photon subtraction (Fig. 5)
- V^(0)_(1,rho) and E^(0)_rho =
from leading-order GEVP of Ref. [8]
assumptions (6)
- domain assumption At order e^2, the finite-volume spectral decomposition approximates the long-distance part of G by the rho and pion-photon states only.
- domain assumption The pion-photon operator counts as order e in the power counting, and cross terms between rho and pion-photon vanish at leading order.
- ad hoc to paper Diagram (D1) is negligible relative to (V), (S), and (F) for the pion-photon reconstruction.
- domain assumption Vector-current renormalization and quark-mass shift counterterms are not needed for the long-distance reconstruction.
- domain assumption The QED_L finite-volume regulator has a well-defined transfer matrix so that the spectral form Eq. (4) holds.
- standard math The diagram relations Eq. (40), derived for specific quark charge assignments, apply to the physical charge assignment through isospin decomposition.
Cite this review
Pith. "Pith review of Long-distance reconstruction of QED corrections to the hadronic vacuum polarization for the muon g-2." pith.science (2026). https://pith.science/paper/VTG3AZ73
@misc{pith2026250821685,
author = {Pith},
title = {Pith review of: Long-distance reconstruction of QED corrections to the hadronic vacuum polarization for the muon g-2},
year = {2026},
howpublished = {\url{https://pith.science/paper/VTG3AZ73}},
note = {Machine review of arXiv:2508.21685}
}
abstract
The long-distance contribution of QED corrections to the hadronic vacuum polarization is particularly challenging to compute in lattice QCD+QED. Currently, it is one of the limiting factors towards matching the precision of the recent result by the Fermilab E989 experiment for the muon g-2. In this work, we present a method for obtaining high-precision results for this contribution by reconstructing exclusive finite-volume state contributions. We find relations between the pion-photon contributions of individual diagrams and demonstrate the reconstruction method with lattice QCD+QED data at a single lattice spacing of $a^{-1} \approx 1.73$ GeV and $m_\pi \approx 275$ MeV.
Figures
Figures from the paper (4 more)
Forward citations
Cited by 2 Pith papers
-
Isospin-breaking effects in inclusive hadronic $\tau$ data for the muon $(g-2)$ from first principles
A lattice QCD+QED strategy is outlined for calculating isospin-breaking effects in inclusive tau decays to support high-precision HVP contributions to muon g-2.
-
Field-theoretic versus data-driven evaluations of electromagnetic corrections to hadronic vacuum polarization in $(g-2)_\mu$
Virtual electromagnetic corrections largely cancel radiative-channel contributions in data-driven HVP evaluations for muon g-2, reconciling timelike and spacelike methods via a VMD model.
Reference graph
Works this paper leans on
-
[1]
R. Aliberti et al. , The anomalous magnetic moment of the muon in the Standard Model: an update, (2025), arXiv:2505.21476 [hep-ph]
arXiv 2025
-
[2]
T. Blum, P. A. Boyle, V. G¨ ulpers, T. Izubuchi, L. Jin, C. Jung, A. J¨ uttner, C. Lehner, A. Portelli, and J. T. Tsang (RBC, UKQCD), Calculation of the hadronic vacuum polarization contribution to the muon anomalous magnetic moment, Phys. Rev. Lett. 121, 022003 (2018), arXiv:1801.07224 [hep-lat]
arXiv 2018
- [3]
-
[4]
S. Bors´ anyiet al., Leading hadronic contribution to the muon magnetic moment from lattice QCD, Nature 593, 51 (2021), arXiv:2002.12347 [hep-lat]
arXiv 2021
-
[5]
A. Boccaletti et al. , High precision calculation of the hadronic vacuum polarisation contribution to the muon anomaly, (2024), arXiv:2407.10913 [hep-lat]
arXiv 2024
-
[6]
D. Djukanovic, G. von Hippel, S. Kuberski, H. B. Meyer, N. Miller, K. Ottnad, J. Parrino, A. Risch, and H. Wittig, The hadronic vacuum polarization contribution to the muon g − 2 at long distances, JHEP 04, 098, arXiv:2411.07969 [hep-lat]
-
[7]
J. Parrino, V. Biloshytskyi, E.-H. Chao, H. B. Meyer, and V. Pascalutsa, Computing the UV-finite electromagnetic corrections to the hadronic vacuum polarization in the muon (g − 2) from lattice QCD, JHEP 07, 201, arXiv:2501.03192 [hep-lat]
-
[8]
T. Blum et al. (RBC, UKQCD), The long-distance window of the hadronic vacuum polarization for the muon g-2, Phys. Rev. Lett. 134, 201901 (2025), arXiv:2410.20590 [hep-lat]
arXiv 2025
Show all 18 references
-
[9]
Bruno, T
M. Bruno, T. Izubuchi, C. Lehner, and A. S. Meyer, Exclusive Channel Study of the Muon HVP, PoS LA TTICE2019, 239 (2019), arXiv:1910.11745 [hep-lat]
2019 arXiv
-
[10]
Blum et al
T. Blum et al. (RBC, UKQCD), Update of Euclidean windows of the hadronic vacuum polarization, Phys. Rev. D 108, 054507 (2023), arXiv:2301.08696 [hep-lat]
2023 arXiv
-
[11]
Hayakawa and S
M. Hayakawa and S. Uno, QED in finite volume and finite size scaling effect on electromagnetic properties of hadrons, Prog. Theor. Phys. 120, 413 (2008), arXiv:0804.2044 [hep-ph]
2008 arXiv
-
[12]
Di Carlo, M
M. Di Carlo, M. T. Hansen, N. Hermansson-Truedsson, and A. Portelli, QED r: a finite-volume QED action with redis- tributed spatial zero-momentum modes, (2025), arXiv:2501.07936 [hep-lat]
2025
-
[13]
T. Blum, N. Christ, M. Hayakawa, T. Izubuchi, L. Jin, C. Jung, C. Lehner, and C. Tu (RBC, UKQCD), Hadronic light- by-light contribution to the muon anomaly from lattice QCD with infinite volume QED at physical pion mass, Phys. Rev. D 111, 014501 (2025), arXiv:2304.04423 [hep-lat]
2025 arXiv
-
[14]
Bernecker and H
D. Bernecker and H. B. Meyer, Vector Correlators in Lattice QCD: Methods and applications, Eur. Phys. J. A 47, 148 (2011), arXiv:1107.4388 [hep-lat]
2011 arXiv
-
[15]
Lehner et al., Grid Python Toolkit (GPT)
C. Lehner et al., Grid Python Toolkit (GPT)
-
[16]
Boyle et al., Grid
P.A. Boyle et al., Grid. 11
-
[17]
P. A. Boyle, G. Cossu, A. Yamaguchi, and A. Portelli, Grid: A next generation data parallel C++ QCD library, PoS LA TTICE2015, 023 (2016)
2016
-
[18]
Yamaguchi, P
A. Yamaguchi, P. Boyle, G. Cossu, G. Filaci, C. Lehner, and A. Portelli, Grid: OneCode and FourAPIs, PoS LA T- TICE2021, 035 (2022), arXiv:2203.06777 [hep-lat]
2022 arXiv
Reviewed August 5, 2026 · model on record in the stance chip above.
Discussion (0). Sign in to comment.