arxiv: 2604.19430 · v1 · submitted 2026-04-21 · ✦ hep-lat

Recognition: unknown

An Update on the Isospin-Breaking Effects in the Pion Decay Constant with Staggered Quarks

Alessandro Cotellucci , Davide Giusti

Authors on Pith no claims yet

Pith reviewed 2026-05-10 01:16 UTC · model grok-4.3

classification ✦ hep-lat

keywords isospin breakingpion decay constantstaggered quarkslattice QCDQED_Laxial correlatorgradient flow scale

0 comments

The pith

Lattice calculations with staggered quarks are isolating isospin-breaking corrections to the pion decay constant.

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

The paper reports progress on an ongoing BMW Collaboration effort to compute electromagnetic isospin-breaking effects in the pion decay constant. It updates the isosymmetric baseline value and the gradient-flow scale w0, then shows early valence-valence pieces of the axial-pseudoscalar correlator on ensembles that vary in volume and lattice spacing. These pieces are computed with Nf=2+1+1 staggered quarks at near-physical pion mass and QED_L. A reader would care because the size of these corrections directly affects the precision with which lattice results can be compared to experiment in tests of the Standard Model.

Core claim

The authors present an update on the computation of isospin-breaking effects in the pion decay constant using Nf=2+1+1 staggered quarks at near-physical pion mass with QED_L. They provide the current isosymmetric value and its use in determining the gradient-flow scale w0, display preliminary valence-valence contributions to the axial-pseudoscalar correlator across several volumes and lattice spacings, and outline the remaining steps needed to complete the calculation.

What carries the argument

The valence-valence contribution to the axial-pseudoscalar correlator evaluated on staggered-quark lattices that include QED_L.

If this is right

A completed calculation would supply a controlled lattice value for the electromagnetic shift in the pion decay constant.
The updated w0 scale would tighten the overall calibration of the lattice ensembles.
Control over volume and spacing dependence in the correlator would reduce systematic errors in the final isospin-breaking correction.
The framework would be ready to incorporate sea-quark isospin breaking in subsequent work.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

The same valence-valence technique could be applied to other leptonic decay constants to obtain consistent electromagnetic corrections.
Once available, the result would allow a more precise extraction of CKM matrix elements from pion and kaon decays.
Comparison with analytic predictions at this stage could already test whether the dominant effects are being captured before full sea-quark inclusion.

Load-bearing premise

The staggered action with QED_L and near-physical masses captures the dominant isospin-breaking contributions to the correlator without large uncontrolled discretization or finite-volume artifacts.

What would settle it

If the continuum-extrapolated valence-valence contribution fails to match the size and sign predicted by chiral perturbation theory for the electromagnetic correction to f_pi, the current approach would require substantial revision.

Figures

Figures reproduced from arXiv: 2604.19430 by Alessandro Cotellucci, Davide Giusti.

**Figure 1.** Figure 1: Status of the determination of the 𝑉𝑢𝑑 and 𝑉𝑢𝑠 elements of the CKM matrix as reported by the Flavour Lattice Averaging Group [1]. In Isosymmetric QCD the pion decay constant 𝑓𝜋 is defined by the hadronic matrix element: ⟨0| 𝑢𝛾¯ 4𝛾5𝑑 |𝜋(0®)⟩ = 𝑀𝜋 𝑓𝜋. (1) To compute the pion decay constant in QCD+QED 𝐹𝜋 one needs to consider the leptonic decay rate 𝜋 + → 𝑙 + 𝜈𝑙(𝛾), in the PDG parametrization it is: 𝐹 2 𝜋 = Γ… view at source ↗

**Figure 2.** Figure 2: Landscape of the ensembles used for the isoQCD determination of 𝑤0 𝑓𝜋 around the physical point from [6]. The analysis follows the strategy presented in [8]. To have a better control over the continuum extrapolation we add the determination of 𝑤0 using the Zeuthen-flow [9] in addition to the Wilsonflow [10], [PITH_FULL_IMAGE:figures/full_fig_p004_2.png] view at source ↗

**Figure 3.** Figure 3: Continuum limit extrapolation of the iso QCD value of 𝑤0 𝑓𝜋. The upper plot represents the continuum extrapolation both in the lattice spacing and in the taste-symmetry breaking terms Δ𝐾𝑆 (computed by comparing the masses of different meson tastes on our ensembles); in the lower plot the probability distribution function (PDF) including both statistical and systematic variations is shown in red. Median 0.1… view at source ↗

**Figure 4.** Figure 4: Renormalized sea quark electromagnetic derivative of 𝑤0𝐹𝜋 divided by its QCD value. The left panel shows the continuum, the right panel the infinite volume extrapolation [PITH_FULL_IMAGE:figures/full_fig_p006_4.png] view at source ↗

**Figure 5.** Figure 5: shows the preliminary result for the effective factorizable valence derivative for the four different volumes considered. −0.2 −0.1 0.0 ∂ 2 w 0 Feff π /∂e 2 v L =3.2 fm Preliminary −0.02 −0.01 0.00 0.01 0.02 L =4.2 fm Preliminary 1 2 3 4 x4/fm −0.01 0.00 0.01 0.02 ∂ 2 w 0 Feff π /∂e 2 v L =6.3 fm Preliminary 1 2 3 4 x4/fm 0.005 0.010 0.015 0.020 L =8.4 fm Preliminary E.M. derivative a = 0.1315 fm [PITH_FU… view at source ↗

**Figure 6.** Figure 6: Factorizable valence electromagnetic derivative of 𝑤0𝐹𝜋 for 𝑎 = 0.0952 fm. where the precision is dominated by the continuum extrapolation of the isosymmetric part. The determination of isospin-breaking effects was done by combining sea quark effects from BMW ensembles and the valence-quark effects from the RM123 work. We are working on improving the continuum extrapolation (by adding a finer lattice spaci… view at source ↗

read the original abstract

We present an update on the ongoing computation of the isospin-breaking effects in the Pion Decay Constant from the BMW Collaboration. The calculation is carried out with N$_f$=2+1+1 staggered quarks with a near-physical pion mass and QED$_{\text{L}}$. We give an update on the isosymmetric value and the current determination used to compute the gradient-flow scale $w_{0}$, then we present some preliminary results of the valence-valence contribution to the axial-pseudoscalar correlator for different volumes and lattice spacings. We also discuss the next steps and plans.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Desk Editor's Note private letter to a colleague

This is a brief status update on BMW's ongoing isospin-breaking pion decay constant work, with an isosymmetric w0 refresh and some preliminary valence-valence correlator data but no final numbers or extrapolations.

read the letter

This is an interim update from the BMW group on their calculation of isospin-breaking effects in the pion decay constant. They report an update to the isosymmetric value and w0 scale, plus some preliminary valence-valence axial-pseudoscalar correlator results at varying volumes and spacings, but no final extrapolated correction or error analysis. The paper does a decent job of showing the current status of their staggered quark setup with QED_L at near-physical pion mass. Continuing from prior work, the preliminary data suggest they are making progress on the correlators before tackling the full isospin-breaking contributions. Where it falls short is in the lack of any concrete numbers or assessment of systematics. Without fit details or error bars, it's impossible to evaluate whether the chosen lattice spacings and volumes are adequate to keep artifacts small, which is the key assumption here. The work inherits parameters from earlier isosymmetric calculations, but the abstract doesn't indicate any new circularity issues. This paper is primarily for specialists in lattice flavor physics who follow the BMW collaboration's isospin-breaking program. It won't offer much to someone outside that niche looking for a new result on f_pi or CKM-related quantities. I wouldn't cite this in my work, as it contains no citable outcome. It probably doesn't merit sending to peer review yet, since it's more of a status report than a finished piece of research. The authors should complete the analysis first.

Referee Report

0 major / 3 minor

Summary. The manuscript is an update from the BMW Collaboration on their ongoing lattice QCD computation of isospin-breaking effects in the pion decay constant. It employs N_f=2+1+1 staggered quarks at near-physical pion mass with QED_L, reports an update to the isosymmetric value and the gradient-flow scale w_0, presents preliminary valence-valence contributions to the axial-pseudoscalar correlator across a few volumes and lattice spacings, and outlines next steps.

Significance. If completed with controlled extrapolations and uncertainties, the calculation would provide a valuable lattice determination of isospin-breaking corrections to f_π, relevant for precision tests of the Standard Model and CKM phenomenology. The staggered+QED_L setup is a standard choice for such studies, and the preliminary correlator data at multiple spacings and volumes represent incremental progress toward controlling discretization and finite-volume effects.

minor comments (3)

Abstract: the preliminary results for the valence-valence axial-pseudoscalar correlator are described without any error bars, fit details, or numerical values, which limits the reader's ability to assess the current precision or trends with volume and spacing.
The update on the isosymmetric w_0 determination (mentioned after the abstract) should explicitly state the fitting procedure or reference the prior isosymmetric work from which parameters are inherited, to clarify any potential circularity in the scale setting.
The discussion of next steps would benefit from a brief outline of the planned chiral/continuum extrapolations and how the full isospin-breaking correction (including sea-quark contributions) will be assembled from the valence-valence piece shown.

Simulated Author's Rebuttal

0 responses · 0 unresolved

We thank the referee for the positive evaluation of our manuscript and the recommendation for minor revision. The work is indeed a preliminary update on our ongoing N_f=2+1+1 staggered calculation with QED_L of isospin-breaking corrections to f_π, including updates to the isosymmetric value, w_0 scale setting, and sample valence-valence correlator data.

Circularity Check

0 steps flagged

No significant circularity detected

full rationale

The paper is a status update on an ongoing lattice QCD computation using staggered quarks with QED_L. It reports an update to the isosymmetric w0 scale (from prior BMW work) and preliminary valence-valence axial-pseudoscalar correlator results at a few volumes and spacings, without any final extrapolated value, derivation, or prediction of isospin-breaking corrections. No equation or claim reduces by construction to a fitted input, self-definition, or self-citation chain; the work remains computational and preliminary with no load-bearing circular steps.

Axiom & Free-Parameter Ledger

0 free parameters · 0 axioms · 0 invented entities

No explicit free parameters, axioms, or invented entities are identifiable from the abstract. The calculation relies on standard staggered QCD and QED_L techniques from prior literature.

pith-pipeline@v0.9.0 · 5399 in / 1212 out tokens · 25750 ms · 2026-05-10T01:16:23.127929+00:00 · methodology

discussion (0)

Reference graph

Works this paper leans on

15 extracted references · 15 canonical work pages · 3 internal anchors

[1]

FLAG Review 2024

Aoki, Y. and others. FLAG Review 2024. 2024. arXiv:2411.04268

work page internal anchor Pith review arXiv 2024
[2]

and Giusti, D

Di Carlo, M. and Giusti, D. and Lubicz, V. and Martinelli, G. and Sachrajda, C. T. and Sanfilippo, F. and Simula, S. and Tantalo, N. Light-meson leptonic decay rates in lattice QCD+QED. Phys. Rev. D. 2019. doi:10.1103/PhysRevD.100.034514. arXiv:1904.08731

work page doi:10.1103/physrevd.100.034514 2019
[3]

and Lubicz, V

Giusti, D. and Lubicz, V. and Martinelli, G. and Sachrajda, C. T. and Sanfilippo, F. and Simula, S. and Tantalo, N. and Tarantino, C. First lattice calculation of the QED corrections to leptonic decay rates. Phys. Rev. Lett. 2018. doi:10.1103/PhysRevLett.120.072001. arXiv:1711.06537

work page doi:10.1103/physrevlett.120.072001 2018
[4]

and Hansen, M

Di Carlo, M. and Hansen, M. T. and Portelli, A. and Hermansson-Truedsson, N. Relativistic, model-independent determination of electromagnetic finite-size effects beyond the pointlike approximation. Phys. Rev. D. 2022. doi:10.1103/PhysRevD.105.074509. arXiv:2109.05002

work page doi:10.1103/physrevd.105.074509 2022
[5]

Lattice Calculation of Light Meson Radiative Leptonic Decays

Christ, Norman H. and Feng, Xu and Izubuchi, Taku and Jin, Luchang and Sachrajda, Christopher T. and Tuo, Xin-Yu. Lattice Calculation of Light Meson Radiative Leptonic Decays. 2025. arXiv:2510.26993

work page internal anchor Pith review Pith/arXiv arXiv 2025
[6]

Isospin-breaking corrections to light-meson leptonic decays from lattice simulations at physical quark masses

Boyle, Peter and others. Isospin-breaking corrections to light-meson leptonic decays from lattice simulations at physical quark masses. JHEP. 2023. doi:10.1007/JHEP02(2023)242. arXiv:2211.12865

work page doi:10.1007/jhep02(2023)242 2023
[7]

and Martinelli, G

Lubicz, V. and Martinelli, G. and Sachrajda, C. T. and Sanfilippo, F. and Simula, S. and Tantalo, N. Finite-Volume QED Corrections to Decay Amplitudes in Lattice QCD. Phys. Rev. D. 2017. doi:10.1103/PhysRevD.95.034504. arXiv:1611.08497

work page doi:10.1103/physrevd.95.034504 2017
[8]

Hybrid calculation of hadronic vacuum polarization in muon g-2 to 0.48\%

Boccaletti, A. and others. High precision calculation of the hadronic vacuum polarisation contribution to the muon anomaly. 2024. arXiv:2407.10913

work page internal anchor Pith review Pith/arXiv arXiv 2024
[9]

Properties and uses of the Wilson flow in lattice QCD

L. Properties and uses of the Wilson flow in lattice QCD. JHEP. 2010. doi:10.1007/JHEP08(2010)071. arXiv:1006.4518

work page doi:10.1007/jhep08(2010)071 2010
[10]

and Sint, S

Ramos, A. and Sint, S. Symanzik improvement of the gradient flow in lattice gauge theories. Eur. Phys. J. C. 2016. doi:10.1140/epjc/s10052-015-3831-9. arXiv:1508.05552

work page doi:10.1140/epjc/s10052-015-3831-9 2016
[11]

and others

Aoki, Y. and others. FLAG Review 2021. Eur. Phys. J. C. 2022. doi:10.1140/epjc/s10052-022-10536-1. arXiv:2111.09849

work page doi:10.1140/epjc/s10052-022-10536-1 2021
[12]

and others

Borsanyi, Sz. and others. Leading hadronic contribution to the muon magnetic moment from lattice QCD. Nature. 2021. doi:10.1038/s41586-021-03418-1. arXiv:2002.12347

work page doi:10.1038/s41586-021-03418-1 2021
[13]

and others

Stokes, Finn M. and others. BMW/DMZ calculation of the hadronic vacuum polarisation for the muon magnetic moment. 2026. arXiv:2603.03835

work page arXiv 2026
[14]

Workman, R. L. and others. Review of Particle Physics. PTEP. 2022. doi:10.1093/ptep/ptac097

work page doi:10.1093/ptep/ptac097 2022
[15]

Leading isospin-breaking corrections to pion, kaon, and charmed-meson masses with twisted- mass fermions

Giusti, D. and Lubicz, V. and Tarantino, C. and Martinelli, G. and Sanfilippo, F. and Simula, S. and Tantalo, N. Leading isospin-breaking corrections to pion, kaon and charmed-meson masses with Twisted-Mass fermions. Phys. Rev. D. 2017. doi:10.1103/PhysRevD.95.114504. arXiv:1704.06561

work page doi:10.1103/physrevd.95.114504 2017