REVIEW 2 major objections 4 minor 1 cited by
$\Delta I =1/2$ process of $K\to\pi\pi$ decay on multiple ensembles with periodic boundary conditions
T0 review · 2 major / 4 minor · reviewed 2026-08-10 · deepseek-v4-flash
Pith's one-line read The paper reports the first periodic-boundary-condition lattice-QCD continuum extrapolation of ε'/ε, finding Re(ε'/ε)=17.5(6.8)(4.9)(5.0)×10^-4, in agreement with experiment.
desk verdict First PBC continuum extrapolation of ε' is a real step, but the new 1.4 GeV point has a postponed on-shell interpolation systematic that makes the central value conditional. 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 objects are the $\Delta I = 1/2$ amplitude $A_0$ and the ratio $\varepsilon'/\varepsilon$, built from renormalized matrix elements $M_i^{\overline{\mathrm{MS}}}(\mu)$ and three-flavor Wilson coefficients $z_i(\mu)$ and $y_i(\mu)$. The two-pion final state is isolated with a variational GEVP over momentum and $\sigma$ operators, the Lellouch-Lüscher factor converts finite-volume states to infinite-volume matrix elements, linear interpolation in $E_{\pi\pi}/m_K$ reaches the on-shell point, and RI/SMOM step scaling to $\mu = 4.0$ GeV supplies the nonperturbative renormalization. The continuum limit is taken by a two-point linear fit of $A_0$ in $a^2$, with the scaling-violation error estimated from the shift to an $O(a^4)$ fit.
What would settle it
Compute the same $\Delta I = 1/2$ amplitude and $\mathrm{Re}(\varepsilon'/\varepsilon)$ on an ensemble with $a^{-1}\approx 2.7$ GeV using the identical periodic-boundary-condition analysis; if the new point lies off the $O(a^2)$ line by more than the quoted scaling-violation error, the continuum extrapolation and its uncertainty would be wrong.
Extended reading notes
Core claim
On its own terms, this paper establishes that the $\Delta I = 1/2$ $K\to\pi\pi$ amplitude and $\varepsilon'$ can be computed with periodic boundary conditions at two lattice spacings and extrapolated to the continuum. The central number is $\mathrm{Re}(\varepsilon'/\varepsilon) = 17.5(6.8)(4.9)(5.0)\times 10^{-4}$, which agrees with both the experimental value $16.6(2.3)\times 10^{-4}$ and the earlier G-parity lattice result $21.7(2.6)(6.2)(5.0)\times 10^{-4}$. The paper also reports that the bare matrix elements of the electroweak penguin operators $Q_7$ and $Q_8$ show some tension between the periodic and G-parity calculations at the shared inverse spacing of about 1.4 GeV, while the other operators are consistent. All of these results are labeled preliminary.
Load-bearing premise
The continuum numbers rest on the assumption that discretization errors behave as the square of the lattice spacing across the two spacings used here, so adding a third, finer lattice could move the central value.
Editorial extensions
If this is right
- The value Re(ε'/ε)=17.5(6.8)(4.9)(5.0)×10^-4 is compatible with the experimental world average, so the Standard Model remains consistent with direct CP violation in kaon decays at current precision.
- The consistency of most bare matrix elements between the periodic and G-parity calculations at the shared 1.4 GeV spacing supports periodic boundary conditions as a viable route for this process.
- The continuum extrapolation, though based on two spacings, gives the first estimate of ε' in which a discretization error is assigned rather than quoted as a single-spacing uncertainty.
- The paper's stated outlook is that adding finer periodic-boundary-condition ensembles with inverse spacings up to 2.7 GeV should reduce the finite-lattice-spacing error, and that the dominant remaining within-isospin error is the ~12% perturbative truncation of the Wilson coefficients.
Reading between the lines
- If the linear $a^2$ extrapolation holds on finer lattices, the periodic-boundary-condition route could reach the experimental precision for $\varepsilon'$ using already-generated domain-wall ensembles, avoiding the overhead of G-parity boundary conditions.
- The reported tension in the $Q_7$ and $Q_8$ bare matrix elements between the two boundary-condition treatments suggests that finite-volume or boundary-condition effects on the left-right penguin operators deserve a dedicated study before the final error budget is closed.
- A direct test of the claimed order-independence of interpolation and renormalization would be to repeat the earlier procedure, renormalizing before the energy interpolation, on the new 1.4 GeV ensemble and compare the on-shell matrix elements.
- Once a third spacing is included, the current $O(a^2)$ scaling-violation estimate, based on only two relatively coarse points, is likely to be revised; the fine ensembles already generated provide the data for that test.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. This proceedings paper reports preliminary lattice QCD results for the ΔI = 1/2 K→ππ matrix elements and the direct CP-violation parameter ε′/ε, computed with periodic boundary conditions (PBC) on two ensembles with inverse lattice spacings a⁻¹ ≈ 1.0 GeV and 1.4 GeV. The coarser ensemble updates an earlier PBC calculation with roughly doubled statistics, while the finer ensemble is new. The authors perform a GEVP analysis to isolate the first-excited two-pion state, apply Lellouch–Lüscher factors, interpolate bare matrix elements to the on-shell energy, renormalize in RI/SMOM with step scaling, and combine the renormalized amplitudes with external Wilson coefficients. They present the first continuum extrapolation of the ΔI = 1/2 amplitude and of ε′/ε from PBC data, obtaining Re(ε′/ε) = 17.5(6.8)(4.9)(5.0)×10⁻⁴, consistent with experiment and with the earlier GPBC result, and estimate an O(a²) scaling-violation systematic by comparing with an O(a⁴) fit. The paper is explicitly preliminary, with several systematic effects deferred to a later full publication.
Significance. If correct, this work provides the first continuum-extrapolated value of ε′/ε from lattice QCD with periodic boundary conditions, a milestone that could eventually reduce the dominant finite-spacing systematic in this quantity. The analysis follows well-established procedures from the RBC/UKQCD program: the GEVP state isolation, A2A/AMA methods, RI/SMOM renormalization with step scaling, and Lellouch–Lüscher factor treatment are all standard and internally consistent. The paper also honestly acknowledges the limitations of a two-point continuum extrapolation on coarse lattices and the preliminary nature of the result. The main risk is that the central claim — the continuum-extrapolated ε′/ε — depends on a new finer ensemble whose on-shell interpolation systematic is unquantified and explicitly postponed to the full paper. That omission directly affects the headline number and the conclusions drawn from it.
major comments (2)
- [Section 3, Figure 3 and text on the 32³ lattice] The on-shell interpolation systematic for the new finer ensemble is left unquantified. The paper states: 'On the 32³ lattice, the n = 1 two-pion energy is about 11% larger than the kaon mass and this could enhance the systematic error, but we postpone the estimation of this effect until the full paper.' This ensemble is the finer point in the continuum extrapolation (Figure 4), so the central result Re(ε′/ε) = 17.5(6.8)(4.9)(5.0)×10⁻⁴ does not include this systematic, despite the interpolation gap being roughly twice as large as on the 1.0 GeV lattice (6%). The omission is load-bearing: if the interpolation bias at 1.4 GeV is not negligible, the continuum extrapolation and thus the quoted ε′/ε would shift by an amount not covered by the quoted errors. This concern is reinforced by the 'some tension' observed for Q₇ and Q₈, which are important for Im A₀, at this same lattice spacing. The authors should either provide an estimate of this systematic (e.g., from a comparison of linear and quadratic interpolation, or from the GPBC data) or explicitly exclude it from the quoted central value and error budget, with a clear statement that the presented continuum result is conditional on this systematic being small.
- [Section 3, Figure 4 and continuum limit discussion] The continuum extrapolation uses only two lattice spacings with an assumed O(a²) form, and the scaling-violation systematic is estimated by comparing with an O(a⁴) fit. With only two data points, an O(a⁴) alternative is not a test of the extrapolation; it simply shifts the curve by a fixed amount. The paper acknowledges this ('somewhat uncertain since we use only two lattices with relatively large spacings'), but the quoted 11% scaling-violation error for Im A₀ should be understood as a qualitative measure rather than a robust estimate. The authors should make this limitation more prominent in the abstract and summary, and ideally present the O(a⁴) central value explicitly so that readers can see the sensitivity. This issue is separate from the interpolation systematic of the previous comment but compounds its impact on the final ε′/ε.
minor comments (4)
- [Section 2] The phrase 'the244 lattice' appears in the discussion of the coarser ensemble; this is likely a typo for 'the 24³×64 lattice' or 'the 24³ lattice'. Please clarify.
- [Section 2, Eq. (6)] In the displayed equation for C³ᵖᵗₙ,ᵢ, the angle brackets are rendered as 'D' and 'E' (apparently a LaTeX/encoding issue). The expectation value notation should be fixed in the final version.
- [Section 4] The word 'Möbuis' in the sentence about existing ensembles is a typo for 'Möbius'.
- [References] Reference [9] is a PoS proceedings contribution by C. Kelly; it would be helpful to include the full author list and arXiv identifier if available, to improve traceability of the ensemble parameters.
Circularity Check
No significant circularity: the central ε′ result is a direct lattice computation with external Wilson coefficients and benchmark comparisons.
full rationale
The derivation chain is self-contained: the new ingredient that dominates the ε′ prediction, Im A0, is obtained from lattice three-point functions, GEVP-optimized two-pion operators, Lellouch-Lüscher factors, and RI/SMOM renormalization with step scaling, none of which is defined in terms of the final ε′ value. Wilson coefficients are taken from external references [21,22] and are not fitted here. The ε′ formula (Eq. 8) uses experimental Re A0 and Re A2 and the earlier ΔI=3/2 result for Im A2 [4]; this is a standard input decomposition, not a definition of the target in terms of itself, and the dominant Im A0 term is computed in this paper. The PBC-vs-GPBC comparison is used as a consistency check; the GPBC 2020 result is not an input to the continuum fit. The continuum extrapolation (Figure 4) is a two-point linear fit in a² with a scaling-violation estimate from an O(a⁴) alternative; it is statistically weak, as the paper acknowledges, but not circular. The postponed estimate of the on-shell interpolation systematic on the 32³ lattice (Section 3) and the observed Q7/Q8 tension with GPBC are accuracy and systematic-error concerns, not instances where a prediction reduces by construction to an input. No self-citation chain is load-bearing: prior RBC/UKQCD work supplies methodology, ensembles, and the small Im A2 input, but the central continuum-extrapolated result is independently computed here.
Assumptions & free parameters
free parameters (3)
- O(a^2) slope of A0 continuum extrapolation
- Linear interpolation coefficient to the kaon mass for finite-volume matrix elements
- O(a^4) alternative extrapolation value used for the scaling-violation systematic
assumptions (5)
- domain assumption The Lellouch-Lüscher factor converts finite-volume two-pion matrix elements to infinite volume for the two lowest states in the elastic region.
- domain assumption O(a^2) scaling is the leading discretization error for these domain wall fermion ensembles.
- domain assumption Linear interpolation of matrix elements to the on-shell point introduces a controlled systematic; the 6% offset on the 243 lattice is small and the 11% offset on the 323 lattice is postponed.
- domain assumption RI/SMOM step-scaling to 4 GeV with perturbative conversion to MS is valid.
- domain assumption The GEVP with optimized two-pion operators isolates the first-excited two-pion state.
Cite this review
Pith. "Pith review of $\Delta I =1/2$ process of $K\to\pi\pi$ decay on multiple ensembles with periodic boundary conditions." pith.science (2026). https://pith.science/paper/2QRCBBXY
@misc{pith2026250118077,
author = {Pith},
title = {Pith review of: $\Delta I =1/2$ process of $K\to\pi\pi$ decay on multiple ensembles with periodic boundary conditions},
year = {2026},
howpublished = {\url{https://pith.science/paper/2QRCBBXY}},
note = {Machine review of arXiv:2501.18077}
}
abstract
We present our preliminary results for the $\Delta I = 1/2$ matrix elements of $K\to\pi\pi$ decay and $\varepsilon'$, the measure of direct $CP$ violation in $K\to\pi\pi$, computed on multiple ensembles with periodic boundary conditions (PBC) at the inverse lattice spacings of $a^{-1} \approx 1.0$~GeV and 1.4~GeV. The finer lattice ensemble is newly introduced as an extension to the first PBC calculation [1], while the calculation on the coarser ensemble is updated with the approximately doubled statistics. Our first attempt to take the continuum limit is also discussed with acknowledging potential significance of the $O(a^2)$ scaling violation on these coarse lattices.
Figures
Figures from the paper (2 more)
Forward citations
Cited by 1 Pith paper
-
From scattering towards multi-hadron weak decays
A review of current lattice QCD scattering calculations shows that finite-volume formalisms now enable multi-hadron weak decay studies with direct relevance to flavour physics.
Reference graph
Works this paper leans on
-
[1]
T. Blum et al. [RBC and UKQCD], Phys. Rev. D 108, no.9, 094517 (2023) doi:10.1103/PhysRevD.108.094517 [arXiv:2306.06781 [hep-lat]]
arXiv 2023
-
[2]
T. Blum, P. A. Boyle, N. H. Christ, N. Garron, E. Goode, T. Izubuchi, C. Jung, C. Kelly, C. Lehner and M. Lightman, et al. Phys. Rev. Lett. 108, 141601 (2012) doi:10.1103/PhysRevLett.108.141601 [arXiv:1111.1699 [hep-lat]]
work page Pith review arXiv 2012
-
[3]
T. Blum, P. A. Boyle, N. H. Christ, N. Garron, E. Goode, T. Izubuchi, C. Jung, C. Kelly, C. Lehner and M. Lightman, et al. Phys. Rev. D 86, 074513 (2012) doi:10.1103/PhysRevD.86.074513 [arXiv:1206.5142 [hep-lat]]
work page Pith review arXiv 2012
-
[4]
T. Blum, P. A. Boyle, N. H. Christ, J. Frison, N. Garron, T. Janowski, C. Jung, C. Kelly, C. Lehner and A. Lytle, et al. Phys. Rev. D 91, no.7, 074502 (2015) doi:10.1103/PhysRevD.91.074502 [arXiv:1502.00263 [hep-lat]]
arXiv 2015
- [5]
-
[6]
R. Abbott et al. [RBC and UKQCD], Phys. Rev. D 102, no.5, 054509 (2020) doi:10.1103/PhysRevD.102.054509 [arXiv:2004.09440 [hep-lat]]
arXiv 2020
-
[7]
QED-corrected Lellouch-L\"uscher formula for $K \rightarrow \pi\pi$ decay
Y. Cai and Z. Davoudi, PoS LATTICE2018, 280 (2018) doi:10.22323/1.334.0280 [arXiv:1812.11015 [hep-lat]]
work page Pith review arXiv 2018
-
[8]
$\pi-\pi$ scattering, QED and finite-volume quantization
N. Christ, X. Feng, J. Karpie and T. Nguyen, Phys. Rev. D106, no.1, 014508 (2022) doi:10.1103/PhysRevD.106.014508 [arXiv:2111.04668 [hep-lat]]
work page Pith review arXiv 2022
Show all 23 references
-
[9]
Kelly, PoSLATTICE2024, 343 (2024)
C. Kelly, PoSLATTICE2024, 343 (2024)
2024
-
[10]
Foley, K
J. Foley, K. Jimmy Juge, A. O’Cais, M. Peardon, S. M. Ryan and J. I. Skullerud, Comput. Phys. Commun.172, 145-162 (2005) doi:10.1016/j.cpc.2005.06.008 [arXiv:hep-lat/0505023 [hep-lat]]
2005 arXiv
-
[11]
Blum et al
T. Blum et al. [RBC and UKQCD], Phys. Rev. D 104, no.11, 114506 (2021) doi:10.1103/PhysRevD.104.114506 [arXiv:2103.15131 [hep-lat]]
2021 arXiv
-
[12]
T. Blum, T. Izubuchi and E. Shintani, Phys. Rev. D 88, no.9, 094503 (2013) doi:10.1103/PhysRevD.88.094503 [arXiv:1208.4349 [hep-lat]]
2013 arXiv
-
[13]
Shintani, R
E. Shintani, R. Arthur, T. Blum, T. Izubuchi, C. Jung and C. Lehner, Phys. Rev. D91, no.11, 114511 (2015) doi:10.1103/PhysRevD.91.114511 [arXiv:1402.0244 [hep-lat]]
2015 arXiv
-
[14]
Blumet al
T. Blumet al. [RBC and UKQCD], Phys. Rev. D107, no.9, 094512 (2023) [erratum: Phys. Rev. D 108, no.3, 039902 (2023)] doi:10.1103/PhysRevD.107.094512 [arXiv:2301.09286 [hep-lat]]
2023 arXiv
-
[15]
Mcglynn, PoSLATTICE2015, 019 (2016) doi:10.22323/1.251.0019 9 Δ𝐼 = 1/2 process of𝐾→𝜋𝜋 decay on multiple ensembles with periodic boundary conditions
G. Mcglynn, PoSLATTICE2015, 019 (2016) doi:10.22323/1.251.0019 9 Δ𝐼 = 1/2 process of𝐾→𝜋𝜋 decay on multiple ensembles with periodic boundary conditions
2016 doi
-
[16]
Yin and R
H. Yin and R. D. Mawhinney, PoSLATTICE2011, 051 (2011) doi:10.22323/1.139.0051 [arXiv:1111.5059 [hep-lat]]
2011 arXiv
-
[17]
Luscher and U
M. Luscher and U. Wolff, Nucl. Phys. B 339, 222-252 (1990) doi:10.1016/0550- 3213(90)90540-T
1990 doi
-
[18]
B.Blossier,M.DellaMorte,G.vonHippel,T.MendesandR.Sommer,JHEP 04,094(2009) doi:10.1088/1126-6708/2009/04/094 [arXiv:0902.1265 [hep-lat]]
2009 arXiv
-
[19]
Lellouch and M
L. Lellouch and M. Luscher, Commun. Math. Phys. 219, 31-44 (2001) doi:10.1007/s002200100410 [arXiv:hep-lat/0003023 [hep-lat]]
2001 arXiv
-
[20]
Arthur et al
R. Arthur et al. [RBC and UKQCD], Phys. Rev. D 83, 114511 (2011) doi:10.1103/PhysRevD.83.114511 [arXiv:1006.0422 [hep-lat]]
2011 arXiv
-
[21]
A. J. Buras, M. Jamin and M. E. Lautenbacher, Nucl. Phys. B 408, 209-285 (1993) doi:10.1016/0550-3213(93)90535-W [arXiv:hep-ph/9303284 [hep-ph]]
1993 arXiv
-
[22]
Buchalla, A
G. Buchalla, A. J. Buras and M. E. Lautenbacher, Rev. Mod. Phys.68, 1125-1144 (1996) doi:10.1103/RevModPhys.68.1125 [arXiv:hep-ph/9512380 [hep-ph]]
1996 arXiv
-
[23]
Tomii, PoSLATTICE2019, 174 (2020) doi:10.22323/1.363.0174 10
M. Tomii, PoSLATTICE2019, 174 (2020) doi:10.22323/1.363.0174 10
2020 doi
Reviewed August 10, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.