REVIEW 3 major objections 5 minor 39 references
The paper reports the first lattice QCD+QED calculation of the structure-dependent radiative correction to the electron/muon leptonic decay ratio for the pion and kaon at physical quark masses and in the continuum limit, yielding the most p
Reviewed by Pith at T0; open to challenge. T0 means a machine referee read the full paper against a public rubric. the ladder, T0–T4 →
T0 review · deepseek-v4-flash
2026-08-01 05:00 UTC pith:2PCXIZKL
load-bearing objection First lattice calculation of the SD correction to R_e/mu is real work and probably basically right, but the 0.01% error bar is not yet earned: the continuum limit is a two-point linear fit with an unquantified O(a^4) term, and disconnected diagrams are dropped without a numerical bound. the 3 major comments →
First Lattice QCD Determination of Lepton-Flavor-Universality Ratios in Light-Meson Leptonic Decays
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
The central result is a first-principles value for the O(α) structure-dependent radiative correction to Rₑ/μ = Γ(P→eνγ)/Γ(P→μνγ) for P = π, K. The paper computes this correction on physical-mass domain-wall fermion ensembles, using two lattice spacings for a linear-in-a² continuum extrapolation, with a partially quenched correction for the small kaon-mass mismatch between ensembles. The lattice pion real-photon SD term lies noticeably below the previous ChPT value, a shift driven by the photon-momentum dependence of the vector and axial form factors, while the virtual corrections agree within errors. In Coulomb gauge with an error-cancellation ratio RfP, the final predictions are Rₑ/μ = 1.23
What carries the argument
The calculation rests on the infinite-volume reconstruction (IVR) method, extended here to Coulomb-gauge photons. IVR removes the power-law finite-volume errors caused by the massless photon, replacing them with exponential suppression, by splitting the loop integral into a short-distance part evaluated from lattice data and a long-distance part reconstructed from a single-meson intermediate state. The Coulomb gauge is the key new element: because the point-like part of the hadronic tensor satisfies the same Ward identity as the full tensor, the difference is gauge invariant, and Coulomb-gauge photons have no infrared divergence in the diagram of interest, so the point-like contribution is r
Load-bearing premise
The continuum limit is obtained by assuming discretization errors are strictly linear in a² using only two lattice spacings, and an O(a²) short-distance contribution is set to zero, so a third, finer lattice spacing is needed to confirm that residual lattice effects do not shift the central values beyond the quoted uncertainties.
What would settle it
Compute the same SD corrections on a third, finer lattice spacing (around a⁻¹ ≈ 3 GeV) with the same physical pion mass and similar volume; if the new point deviates from the two-point a² line by more than the quoted statistical error, the linear continuum extrapolation is falsified.
If this is right
- The structure-dependent radiative correction is now known from first principles at a precision that no longer dominates the error budget of the Standard-Model Rₑ/μ prediction.
- The lattice value of the pion real-photon SD term implies that the leading ChPT prediction misses a photon-momentum dependence of the form factors; a matched higher-order ChPT calculation should reproduce the shift if the lattice form factors are right.
- Because the Coulomb-gauge IVR scheme suppresses the point-like background, the same machinery can be applied to other QED corrections in meson decays where the signal is a small remainder after large cancellations.
- The quoted O(0.003%) higher-order QED uncertainty will become the limiting theory error once lattice statistics improve; a direct two-loop calculation of the non-logarithmic pieces would then be needed.
Where Pith is reading between the lines
- If the linear-in-a² continuum extrapolation holds, the residual discretization uncertainty is several times smaller than the current statistical error; a third lattice spacing is the cleanest check that the quoted central values do not drift.
- The same Coulomb-gauge IVR framework could be applied to the isospin-breaking correction in the Kμ2/πμ2 ratio, where the point-like subtraction is an order of magnitude larger; success there would extend first-principles control to CKM unitarity tests.
- A direct comparison of the lattice form factors FV(xγ) and FA(xγ) with new high-statistics measurements of radiative pion and kaon decay spectra would test the momentum dependence that drives the pion shift.
- If a future experiment reaches 0.01% precision and agrees with the lattice value but not with ChPT, that would be evidence for the momentum-dependent form-factor effect and would strengthen the case for using lattice input as the benchmark for lepton-universality tests.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. This paper reports the first lattice QCD+QED calculation of the structure-dependent (SD) O(alpha) radiative correction to the lepton-flavor-universality ratio R_{e/mu} for pi and K, using the infinite-volume reconstruction method with Coulomb-gauge photons on RBC/UKQCD physical-pion ensembles. The authors compute the virtual and real-photon SD corrections separately, extrapolate the 48I and 64I results linearly in a^2 to the continuum, and combine them with the known point-like and higher-order QED corrections to obtain R_{e/mu} = 1.23501(10)x10^-4 (pi) and 2.47653(34)x10^-5 (K). The paper claims a 0.01%-level hadronic uncertainty, competitive with PIONEER's projected precision, and presents the results as first-principles Standard-Model benchmarks.
Significance. If the 0.01% precision claim is substantiated, this is an important milestone: it replaces the ChPT-based estimate of the dominant hadronic uncertainty in R_{e/mu} with a first-principles lattice determination at the physical point and in the continuum limit. The paper is carefully structured: the supplemental material contains detailed derivations of the IVR weight functions, analytic point-like subtraction in both gauges, per-ensemble results for all schemes, a channel-by-channel finite-volume decomposition, and the kaon-mass-correction procedure. The Coulomb/Feynman gauge agreement and the R_fP error-cancellation technique are genuine internal cross-checks. However, the central precision claim rests on a two-point linear-in-a^2 continuum extrapolation and on an unquantified omission of quark-disconnected diagrams, so the significance is conditional until those systematics are addressed.
major comments (3)
- [Fig. 3 and Sec. S2.1] The continuum limit is obtained from a linear-in-a^2 extrapolation using only the 48I and 64I ensembles. The paper explicitly defers a third-spacing check to future work, and the t=0, x=0 contribution to the weight functions is set to zero with only the argument that it is O(a^2), with no numerical verification. With two points, any O(a^4) curvature is unconstrained by the fit. A rough estimate using the measured a^2 slope (B ≈ 0.25 GeV^-2 for the pion δ_vir^SD) gives an O(a^4) bias in the intercept of order B a_48I^2 a_64I^2 ≈ 0.015–0.03 × 10^-3, comparable to the statistical error of the extrapolated δ_vir^SD (≈0.05 × 10^-3). The final R_{e/mu} uncertainties (0.008% for pi, 0.014% for K) therefore likely omit a discretization systematic of at least comparable size. The Coulomb/Feynman gauge agreement does not test the linear ansatz, because both gauges use the same two-point extrapolat
- [Numerical analysis (ensemble paragraph)] The calculation omits quark-disconnected contractions without assigning a numerical uncertainty. The statement that these vanish in the SU(3)-flavor limit and are 'expected to be suppressed' is not a quantitative bound. At the claimed 0.01% precision, this omission is a systematic effect that must be at least estimated or bounded; otherwise the error budget is incomplete. Since a direct calculation is left to future work, the paper should include a conservative estimate, for example from a partially quenched or exploratory disconnected-diagram calculation, or an explicit phenomenological bound, before the central precision claim can be taken as final.
- [Eqs. (12)–(13) and Sec. S5] The final R_{e/mu} error is quoted as a single 'stat' error, but δ_real^SD is computed from form factors determined on the same 48I and 64I ensembles as the hadronic functions used for δ_vir^SD (Ref. [28]). The paper does not state whether the statistical errors of δ_vir^SD and δ_real^SD are combined with their correlations taken into account. Since both quantities are derived from the same gauge configurations, treating them as independent could either over- or under-estimate the combined uncertainty. Please clarify the covariance treatment, or justify independence.
minor comments (5)
- [Abstract vs. Eqs. (12)–(13)] The abstract quotes R_{e/mu} = 1.23501(10)x10^-4, while Eq. (12) gives 1.23501(9)_stat(4)_{αn≥2}; the combined one-sigma error is indeed 10 in the last digit. Please add a sentence indicating how the abstract error is obtained from the component errors.
- [Table S3] For the kaon δ_real^SD, the ChPT row is listed as '—' because it is not used in R_{e/mu}. It would be helpful to state explicitly in the table caption that the kaon δ_real^SD is shown only for completeness and is not part of the final R_K result.
- [Sec. S2.1] The statement 'we set the weight function at this point to zero' should specify that this is only the single point t=0, x=0, and should note in the main text that this is an O(a^2) effect removed by the continuum extrapolation, with residual uncertainty covered by the discretization systematic requested above.
- [Sec. S6] The kaon-mass correction uses a linear slope in m_K^2 and neglects the lattice-spacing dependence of the slope; the paper states the shift is small compared to the statistical error. A one-sentence quantitative justification (e.g., the slope is measured on 64I and the shift is below 0.05σ) would make this more transparent.
- [Fig. S3] The figure would benefit from labels showing the statistical error on the corrected 64I point, so the reader can visually compare the size of the kaon-mass correction with the statistical uncertainty.
Circularity Check
No significant circularity: the lattice inputs are not fitted to the target R_e/mu; the self-citation to Ref. [28] is independent support.
full rationale
The paper's derivation chain is self-contained with respect to its target observable. The final R_e/mu is assembled from Eq. (2) as R^(0) times [1 + Delta_alpha,pt + Delta_alpha,SD + Delta_alpha,n>=2]. Only Delta_alpha,SD is new lattice input; the point-like correction and the leading-log higher-order QED terms are taken from external references [33], not from fits. The lattice SD correction is obtained directly from the hadronic tensor H_mu_rho via the IVR weight functions, with the point-like piece subtracted analytically (Sec. S3). No parameter is fitted to the quoted R_e/mu values. The real-photon correction uses the form factors FV and FA from the authors' prior lattice calculation [28]; that is a self-citation, but it is a separately published lattice determination of radiative leptonic decay form factors on the same ensembles, not a quantity constructed to reproduce R_e/mu. The R_fP ratio used for error cancellation is defined from independent decay-constant components (Eq. (10)) and is not a fit to the target ratio. The two-spacing linear-in-a^2 continuum extrapolation is a systematic limitation, but it is a correctness risk rather than circularity: nothing in that extrapolation is equivalent by construction to the result it claims to predict. No step reduces to its own inputs by definition, so the circularity score is 0.
Axiom & Free-Parameter Ledger
free parameters (4)
- f_P (meson decay constant normalization) =
not quoted in text
- R_fP ratio =
≈1 (continuum/infinite-volume limit)
- Meson charge radius (IVR single-particle FV correction) =
not quoted
- Kaon-mass slope ∂δ_vir/∂m_K² =
not quoted; from 64I vs 64I-pq2 (31 configs)
axioms (7)
- standard math O(α) factorization of R_{e/μ} into point-like, structure-dependent, and higher-order QED parts (Eq. 2), with point-like formula from Marciano-Sirlin
- standard math Point-like hadronic tensor H_pt satisfies the same Ward identity as the full QCD hadronic tensor, making the SD subtraction gauge invariant
- domain assumption IVR long-distance contribution is saturated by the ground-state P meson at t=-t_s; residual heavier J^P=1^- states are suppressed
- domain assumption Quark-disconnected contractions vanish in the SU(3)-flavor limit and are negligible at the quoted precision; omitted without a numerical error
- domain assumption Continuum extrapolation is linear in a² using exactly two lattice spacings (48I, 64I); t=0,x=0 weight set to zero as O(a²)
- domain assumption Partially quenched kaon-mass correction is linear in m_K² and neglects valence-sea mismatch
- domain assumption Higher-order QED uncertainty is estimated as 0.003% from two-loop RG running plus a [1-5]×α²/π² estimate
read the original abstract
The ratio of electronic to muonic leptonic decay widths, $R_{e/\mu}$, for the light mesons $\pi$ and $K$, provides a clean test of lepton flavor universality (LFU) and a sensitive probe of physics beyond the Standard Model. Its Standard-Model prediction is exceptionally precise, with the leading uncertainty associated with the structure-dependent (SD) radiative correction of $O(0.1\%)$. As experiments such as PIONEER and NA62 aim for unprecedented precision, this SD correction has become an essential ingredient in precision experiment--theory comparisons. We present the first lattice QCD$+$QED calculation of this SD correction at the physical pion mass and in the continuum limit. We employ the infinite-volume reconstruction (IVR) method with Coulomb-gauge photons, significantly reducing both statistical errors and finite-volume effects. We obtain the Standard-Model predictions, $R_{e/\mu}=1.23501(10)\times10^{-4}$ for $\pi$ and $R_{e/\mu}=2.47653(34)\times10^{-5}$ for $K$. Our results reduce the hadronic uncertainty in $R_{e/\mu}$, provide the most precise Standard-Model predictions to date, and establish first-principles benchmarks for future high-precision tests of LFU.
Figures
Reference graph
Works this paper leans on
-
[1]
D. Bryman, V. Cirigliano, A. Crivellin, and G. Inguglia, Testing Lepton Flavor Universality with Pion, Kaon, Tau, and Beta Decays, Ann. Rev. Nucl. Part. Sci.72, 69 (2022), arXiv:2111.05338 [hep-ph]
Pith/arXiv arXiv 2022
-
[2]
V. Cirigliano and I. Rosell, Two-loop effective theory analysis ofπ(K)→e¯ν e[γ]branching ratios, Phys. Rev. Lett.99, 231801 (2007), arXiv:0707.3439 [hep-ph]
Pith/arXiv arXiv 2007
-
[3]
Impeduglia, R
G. Impeduglia, R. Plano, A. Prodell, N. Samios, M. Schwartz, and J. Steinberger,βDecay of the Pion, Phys. Rev. Lett.1, 249 (1958)
1958
-
[4]
Fazzini, G
T. Fazzini, G. Fidecaro, A. W. Merrison, H. Paul, and A. V. Tollestrup, Electron Decay of the Pion, Phys. Rev. Lett.1, 247 (1958)
1958
-
[5]
R. P. Feynman and M. Gell-Mann, Theory of Fermi in- teraction, Phys. Rev.109, 193 (1958)
1958
-
[6]
E. C. G. Sudarshan and R. e. Marshak, Chirality invari- anceandtheuniversalFermiinteraction,Phys.Rev.109, 1860 (1958)
1958
-
[7]
Fidecaro, The Discoveries of Rare Pion Decays at the CERN Synchrocyclotron, Adv
G. Fidecaro, The Discoveries of Rare Pion Decays at the CERN Synchrocyclotron, Adv. Ser. Direct. High Energy Phys.23, 397 (2015)
2015
-
[8]
O. U. Shanker,π ℓ2,K ℓ3 andK 0− ¯K 0 Constraints on Leptoquarks and Supersymmetric Particles, Nucl. Phys. B204, 375 (1982)
1982
-
[9]
B. A. Campbell and D. W. Maybury, Constraints on scalar couplings fromπ ± →l ±νl, Nucl. Phys. B709, 419 (2005), arXiv:hep-ph/0303046
Pith/arXiv arXiv 2005
-
[10]
B. A. Campbell and A. Ismail, Leptonic Pion Decay And Physics Beyond The Electroweak Standard Model, (2008), arXiv:0810.4918 [hep-ph]
Pith/arXiv arXiv 2008
-
[11]
V. Cirigliano, M. Gonzalez-Alonso, and M. L. Graesser, Non-standardChargedCurrentInteractions: betadecays versus the LHC, JHEP02, 046, arXiv:1210.4553 [hep- ph]
-
[12]
V. Cirigliano and I. Rosell,π/K→e¯νe branching ratios toO(e 2p4)inChiralPerturbationTheory,JHEP10,005, arXiv:0707.4464 [hep-ph]
-
[13]
Aguilar-Arevaloet al.(PiENu), Improved Measure- ment of theπ→eνBranching Ratio, Phys
A. Aguilar-Arevaloet al.(PiENu), Improved Measure- ment of theπ→eνBranching Ratio, Phys. Rev. Lett. 115, 071801 (2015), arXiv:1506.05845 [hep-ex]
Pith/arXiv arXiv 2015
-
[14]
C. Lazzeroniet al.(NA62), Precision Measurement of the Ratio of the Charged Kaon Leptonic Decay Rates, Phys. Lett. B719, 326 (2013), arXiv:1212.4012 [hep-ex]
Pith/arXiv arXiv 2013
-
[15]
W. Altmannshoferet al.(PIONEER), PIONEER: Stud- ies of Rare Pion Decays, (2022), arXiv:2203.01981 [hep- ex]
Pith/arXiv arXiv 2022
-
[16]
W. Altmannshoferet al.(PIONEER), Testing Lepton Flavor Universality and CKM Unitarity with Rare Pion Decays in the PIONEER experiment, inSnowmass 2021 (2022) arXiv:2203.05505 [hep-ex]
Pith/arXiv arXiv 2021
-
[17]
A. Adelmannet al.(PIONEER), European Strategy for Particle Physics Update – PIONEER: a next generation rare pion decay experiment, (2025), arXiv:2504.06375 [hep-ex]
Pith/arXiv arXiv 2025
-
[18]
K. Massri, Unleashing NA62’s full physics potential: Sta- tus and prospects, Talk presented on behalf of the NA62 Collaboration at the XIII International Conference on Kaon Physics (KAON2025), Mainz, Germany (2025), 11 September 2025
2025
-
[19]
N. Carrasco, V. Lubicz, G. Martinelli, C. T. Sachrajda, N. Tantalo, C. Tarantino, and M. Testa, QED Correc- tions to Hadronic Processes in Lattice QCD, Phys. Rev. D91, 074506 (2015), arXiv:1502.00257 [hep-lat]
Pith/arXiv arXiv 2015
-
[20]
D. Giusti, V. Lubicz, C. Tarantino, G. Martinelli, C. T. Sachrajda, F. Sanfilippo, S. Simula, and N. Tantalo, First Lattice Calculation of the QED Corrections to Lep- tonic Decay Rates, Phys. Rev. Lett.120, 072001 (2018), 7 arXiv:1711.06537 [hep-lat]
Pith/arXiv arXiv 2018
-
[21]
M. Di Carlo, G. Martinelli, D. Giusti, V. Lubicz, C. T. Sachrajda, F. Sanfilippo, S. Simula, and N. Tantalo, Light-meson leptonic decay rates in lattice QCD+QED, Phys. Rev. D100, 034514 (2019), arXiv:1904.08731 [hep- lat]
Pith/arXiv arXiv 2019
-
[22]
P. Boyleet al., Isospin-breaking corrections to light- meson leptonic decays from lattice simulations at physi- cal quark masses, JHEP02, 242, arXiv:2211.12865 [hep- lat]
-
[23]
A. Desiderioet al., First lattice calculation of radiative leptonic decay rates of pseudoscalar mesons, Phys. Rev. D103, 014502 (2021), arXiv:2006.05358 [hep-lat]
Pith/arXiv arXiv 2021
-
[24]
R. Frezzotti, M. Garofalo, V. Lubicz, G. Martinelli, C. T. Sachrajda, F. Sanfilippo, S. Simula, and N. Tantalo, Comparison of lattice QCD+QED predictions for radia- tive leptonic decays of light mesons with experimental data, Phys. Rev. D103, 053005 (2021), arXiv:2012.02120 [hep-ph]
Pith/arXiv arXiv 2021
-
[25]
G. Gagliardi, F. Sanfilippo, S. Simula, V. Lubicz, F. Mazzetti, G. Martinelli, C. T. Sachrajda, and N. Tan- talo, Virtual photon emission in leptonic decays of charged pseudoscalar mesons, Phys. Rev. D105, 114507 (2022), arXiv:2202.03833 [hep-lat]
Pith/arXiv arXiv 2022
-
[26]
R. Di Palma, R. Frezzotti, G. Gagliardi, V. Lubicz, G. Martinelli, C. T. Sachrajda, F. Sanfilippo, S. Sim- ula, and N. Tantalo, Kaon radiative leptonic decay rates from lattice QCD simulations at the physical point, Phys. Rev. D111, 114523 (2025), arXiv:2504.08680 [hep-lat]
Pith/arXiv arXiv 2025
-
[27]
X.-Y.Tuo, X.Feng, L.-C.Jin,andT.Wang,LatticeQCD calculation ofK→ℓνℓ ′+ℓ′− decay width, Phys. Rev. D 105, 054518 (2022), arXiv:2103.11331 [hep-lat]
Pith/arXiv arXiv 2022
-
[28]
P. Boyle, N. H. Christ, X. Feng, T. Izubuchi, L. Jin, C. T. Sachrajda, and X.-Y. Tuo, Lattice calculation of light meson radiative leptonic decays, Phys. Rev. D114, 014504 (2026), arXiv:2510.26993 [hep-lat]
Pith/arXiv arXiv 2026
-
[29]
R. Di Palma, R. Frezzotti, G. Gagliardi, V. Lubicz, G. Martinelli, C. T. Sachrajda, F. Sanfilippo, S. Simula, and N. Tantalo, Rare kaon decaysK− →ℓ −¯νℓℓ′+ℓ′−: Standard Model predictions from lattice QCD, (2026), arXiv:2605.22727 [hep-ph]
Pith/arXiv arXiv 2026
-
[30]
R. Di Palma, R. Frezzotti, G. Gagliardi, V. Lubicz, G. Martinelli, C. T. Sachrajda, F. Sanfilippo, S. Sim- ula, and N. Tantalo, Complete lattice QCD calcu- lation ofK − →ℓ −¯νℓℓ ′+ℓ ′− form factors, (2026), arXiv:2605.22742 [hep-lat]
Pith/arXiv arXiv 2026
-
[31]
X. Feng and L. Jin, QED self energies from lattice QCD without power-law finite-volume errors, Phys. Rev. D 100, 094509 (2019), arXiv:1812.09817 [hep-lat]
Pith/arXiv arXiv 2019
-
[32]
N. H. Christ, X. Feng, L.-C. Jin, C. T. Sachrajda, and T. Wang, Radiative corrections to leptonic decays using infinite-volume reconstruction, Phys. Rev. D108, 014501 (2023), arXiv:2304.08026 [hep-lat]
Pith/arXiv arXiv 2023
-
[33]
W. J. Marciano and A. Sirlin, Radiative corrections to πl2 decays, Phys. Rev. Lett.71, 3629 (1993)
1993
-
[34]
Blumet al.(RBC and UKQCD), Domain wall QCD with physical quark masses, Phys
T. Blumet al.(RBC and UKQCD), Domain wall QCD with physical quark masses, Phys. Rev. D93, 074505 (2016), arXiv:1411.7017 [hep-lat]
Pith/arXiv arXiv 2016
-
[35]
Navaset al.(Particle Data Group), Review of particle physics, Phys
S. Navaset al.(Particle Data Group), Review of particle physics, Phys. Rev. D110, 030001 (2024)
2024
-
[36]
A. Grozin, Lectures on QED and QCD, in3rd Dubna International Advanced School of Theoretical Physics (2005) arXiv:hep-ph/0508242. 8 Supplemental Material S1. DERIVATION OF IVR FORMULAS This section derives the scalar-function representation used to evaluate the tensor convolutions in Eq. (7). Sec. S1.1 gives the Coulomb-gauge derivation, while Sec. S1.2 p...
Pith/arXiv arXiv 2005
-
[37]
We first correct the 64I result fromm64I K tom 48I K , δvir SD ⏐⏐ m64I K →m48I K =δ vir SD(m64I K ,mℓ) + ∂δ vir SD ∂m2 K [( m48I K )2 − ( m64I K )2] . (S35)
-
[38]
We then perform thea2-linear continuum extrap- olation using the corrected 64I result together with the 48I result, obtaining the continuum- extrapolated resultδ vir,cont SD (m48I K ,mℓ)atm 48I K = 499.2MeV
-
[39]
Finally, we shift the continuum-extrapolated re- sult to the physical charged-kaon mass,m K± = 493.677MeV, using the same slope, δvir,cont SD (mK±,mℓ) =δ vir,cont SD (m48I K ,mℓ) + ∂δ vir SD ∂m2 K [ m2 K±− ( m48I K )2] . (S36) This step neglects the lattice-spacing dependence of the slope; since the shift itself is smaller than thestatisticaluncertainty, ...
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.