Pith. sign in

REVIEW 3 major objections 7 minor 1 cited by

Split-even approach to the rare kaon decay $K \to \pi \ell^+ \ell^-$

T0 review · 3 major / 7 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 The split-even estimator cuts the dominant statistical noise in lattice rare-kaon decay by 5 to 25 times at equal cost.

desk verdict Useful exploratory methods note on split-even for rare kaon decays; variance reduction is real but its size is under-quantified — referee it, but don't take the 4–10x at face value. read the letter →

arxiv 2501.18358 v1 pith:FT6L2C6Y submitted 2025-01-30 hep-lat hep-ph

Raoul Hodgson , Vera Gülpers , Ryan Hill , Antonin Portelli This is my paper · ORCID
classification hep-lathep-ph
keywords rarekaondecaylatticeQCDsplit-evenestimatorGIMcancellationvariancereductionformfactorfour-pointcorrelationfunctionfrequencysplitting
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

The paper argues that a stochastic estimator designed to compute differences of quark propagators, the split-even estimator, can cure the dominant statistical error in the lattice QCD calculation of the rare kaon decay $K \to \pi \ell^+ \ell^-$. In this decay the light and charm quark loops almost cancel (a GIM subtraction), and the residual noise swamps the signal; the existing physical-point lattice result $a_+ = -0.87(4.44)$ is about eight times less precise than the experimental central value. The paper shows that applying the split-even estimator to the loop and loop-insertion diagrams reduces the error of the four-point correlation function by 5 to 25 times and of the integrated correlators by 4 to 10 times at the same computational cost. If these gains propagate to the final form factor, the calculation would become sensitive to the current theory-experiment discrepancy. The study is presented as exploratory, based on only ten configurations and 32 noise hits.

What carries the argument

The load-bearing object is the split-even estimator for the difference of two quark propagators, built on the algebraic identity $\Delta_L(x) = D_l^{-1}(x|x)-D_c^{-1}(x|x) = (m_c-m_l)D_l^{-1}(x|x)D_c^{-1}(x|x)$, which holds for Wilson and Domain-Wall fermions. For a Loop-Insertion diagram, the difference of products of propagators is rewritten by adding and subtracting a flavour-changing term, yielding a sum of two terms each containing the product $D_l^{-1}D_c^{-1}$ and hence admitting the same insertion of noise fields. Frequency-splitting then replaces $l-c$ by a chain of intermediate mass differences. The estimator uses the same propagators as the standard subtraction, so the cost is unchanged.

What would settle it

Recompute the standard and split-even variances on a few hundred configurations with more noise hits and extract the resulting form-factor uncertainty; if the split-even advantage drops below roughly 3 times, or if the final $a_+$ uncertainty is no better than the earlier physical-point result, the paper's central claim would be undermined.

Watch

Extended reading notes

Core claim

The central claim is that the split-even estimator is applicable to the Loop and Loop-Insertion diagrams in the lattice computation of the rare kaon decay $K^+ \to \pi^+ \ell^+ \ell^-$, which are the dominant source of statistical error in the existing physical-point calculation. Using the identity $D_l^{-1}-D_c^{-1}=(m_c-m_l)D_l^{-1}D_c^{-1}$, it stochastically estimates the light-charm loop difference directly, rather than as a subtraction of two noisy estimates, yielding approximately 5–25 times error reduction on the four-point correlation function and 4–10 times on the integrated correlators at equal computational cost. Frequency-splitting further shows that the residual variance is dominated by the lightest mass difference, so that additional noise hits should be spent there. The paper concludes that, if these reductions propagate to the final form factor, the calculation would become sensitive to the discrepancy between the theoretical and experimental values of $a_+$.

Load-bearing premise

The variance reductions are measured on only 10 configurations with 32 noise hits, and the load-bearing assumption is that this small sample represents the true statistical behaviour and that the gain survives the amplitude extraction into the final form factor.

Editorial extensions

If this is right

  • The dominant stochastic error in the physical-point lattice calculation of $K^+ \to \pi^+ \ell^+ \ell^-$ can be reduced by an order of magnitude at equal computational cost, removing the main obstacle to a precise first-principles prediction.
  • A future full-statistics calculation using the split-even estimator could bring the lattice value of $a_+$ close to experimental precision, testing whether the current discrepancy stems from missing contributions or from new physics.
  • Because Loop-Insertion diagrams contribute little to signal and noise, resources can be redirected from those expensive diagrams to Loop diagrams and to the light-quark-mass part of the spectrum.
  • The frequency-splitting results indicate where additional noise hits are most effective: the lightest mass difference $l-c_1$ dominates the remaining variance.
  • The same estimator is applicable to the neutral kaon mode $K_S \to \pi^0 \ell^+ \ell^-$ and to disconnected diagrams, so the method generalises beyond the charged channel studied here.

Reading between the lines

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

  • If the 5–25 times variance reduction persists at full statistics, the bottleneck shifts from statistical noise to systematics such as the zMöbius-to-Möbius bias correction and disconnected diagrams, which the paper lists as remaining limitations.
  • The same algebraic trick could be applied to other flavour-changing neutral-current processes on the lattice where GIM subtraction creates a loop-noise cancellation, such as rare $B$-meson decays or hadronic vacuum polarization contributions.
  • A decisive test would be to measure the variance ratio on a larger ensemble before committing to a production run; a ten-configuration sample cannot reliably distinguish a 10-fold from a 3-fold improvement.
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.

Referee Report

3 major / 7 minor

Summary. This LATTICE2024 proceedings contribution investigates a variance-reduction strategy for the dominant statistical error in the lattice computation of the rare kaon decay K+→π+ℓ+ℓ−. The target is the GIM-type difference of light and charm loop propagators, D_l^{-1}−D_c^{-1}, which the split-even estimator rewrites as (m_c−m_l)D_l^{-1}D_c^{-1} so that stochastic noise propagates between the two propagators (Eqs. 12–13); the same construction is extended to the Loop-Insertion diagrams by adding and subtracting a flavor-changing current term (Eqs. 14–17). On the RBC-UKQCD physical-point ensemble, using 10 configurations, 6 time translations and 32 noise hits per mass splitting, the authors report an error reduction of 5–25× for the 4-point correlation function and 4–10× for the integrated correlators, find that the Loop-Insertion diagrams contribute subdominantly to the variance, and observe that the lightest mass splitting l−c1 dominates the remaining noise. The stated outlook is that, assuming propagation to the final form factor, the method reaches the level needed to become sensitive to the present theory–experiment discrepancy in a_+.

Significance. If the reported variance reductions are robust, this is a valuable step for rare-kaon physics: the existing physical-point result a_+ = −0.87(4.44) (Eq. 8) is roughly an order of magnitude too imprecise to test the theory–experiment discrepancy, and the method attacks precisely the identified bottleneck, the stochastic light-charm GIM subtraction. The derivation is a genuine strength: Eq. (12) is a standard identity, the LI rewriting in Eqs. (15)–(17) is a correct and transparent extension verifiable by direct algebra, and the derivation contains no fitted parameters or ad hoc assumptions. The paper is also commendably explicit about its limitations (10 configurations, uncorrected zMöbius bias, neglected disconnected diagrams, untested propagation to the form factor), and the observed dominance of the lightest splitting is a concrete, falsifiable prediction that can guide future resource allocation. The main reservation is that the headline 4–10× factors are measured on a sample far too small to pin down their magnitude, so the central claim that the method is 'at the level required' is plausible but not yet established.

major comments (3)
  1. [§3, Figs. 2–4] The headline quantitative claims — approximately 5–25× error reduction for the 4-point correlation function and 4–10× for the integrated correlators — are quoted without any statistical uncertainty on the variance ratios themselves. The underlying variances come from a single sample of 10 gauge configurations with 6 time translations that are correlated by construction; with an effective number of independent samples between 10 and 60, the relative uncertainty on one sample variance is between roughly 20% and 50%, and the uncertainty on a ratio of two such variances is considerably larger, so a nominal 4× reduction is statistically consistent with a true reduction of only about 2×. Because the stated motivation (being sensitive to the theory–experiment discrepancy) depends directly on the magnitude of the reduction, I request that the authors (i) attach jackknife or bootstrap errors over configurations to the variance-ratio estimates, (ii) report the effective number of independent samples, and (iii) clearly separate the spread of the reduction factors across t_H bins or integration ranges from the statistical error of each factor. The quoted ranges appear to be the spread across bins or masses rather than confidence intervals, and the text should say so explicitly.
  2. [§4 Conclusions] The conclusion that the observed improvements 'are at the level required in order to be sensitive to the discrepancy between existing theory and experimental results' is stronger than the evidence supports. Three gaps are acknowledged in the manuscript itself: the zMöbius-to-Möbius bias is stated to be 'expected to be small' without any numerical bound (§3); disconnected diagrams are omitted (§1.2); and the penultimate paragraph concedes that 'it is not certain if all of the improvement will be realised' in the final form-factor analysis. The variance comparison itself is largely unaffected by the zMöbius bias, since both estimators use the same action, but the physics claim is not. I recommend reformulating the conclusion as a conditional statement — the improvement would be at the required level if it propagates through the full analysis — and explicitly listing the missing steps (bias correction with an estimate of its size, disconnected diagrams, extrapolation to the physical charm mass, and propagation to a_+) that separate this variance study from a physics result.
  3. [§2–§3, Figs. 2 and 4] The claim that the comparison in Fig. 2 is made 'at equal computational cost' is not verifiable as written. Eq. (13) establishes the equal-cost statement for a single mass difference, but the numerical comparison uses frequency splitting into five mass differences, each with 32 noise hits, requiring inversions for the additional intermediate masses (a cost the text only labels as a 'trade-off'). Please specify the total number of Dirac inversions used in the standard run versus the split-even run, state whether the same noise fields were used for both estimators, and if possible give relative iteration counts for the different masses, so that the cost normalization of the variance comparison can be checked.
minor comments (7)
  1. [§1.2, §3, §4] There are duplicated-word typos in 'a reasonable signal could could be achieved' (§1.2), 'this study study' (§4), and 'figure fig. 2' (§3); these should be corrected.
  2. [§2, Eq. (13)] The left-side notation D_c^{-1}(η_i^†|x) is used without definition; the text defines D^{-1}(x|η) only, so please add a sentence defining the action of D_c^{-1} on η^†.
  3. [Figs. 3 and 4] The variance panels do not state what the plotted variance is averaged over; please specify whether it is the sample variance over the 10 configurations (or 60 samples) and whether it refers to a single time slice or the integrated quantity.
  4. [Table 1] The Theory row leaves a_S and b_S empty; add a sentence or reference explaining why the theoretical values for the neutral kaon mode are not quoted.
  5. [§3] For the statement that the zMöbius bias is 'expected to be small', give at least a qualitative reference or a bound from a comparison at a single mass, so that the reader can gauge the size of the neglected effect.
  6. [§1.2] The statement that 'the LI diagrams are the most expensive to compute' would be substantiated by a one-sentence accounting of the Dirac inversions required per noise hit for the L and LI estimators.
  7. [Abstract and §3] The abstract and the conclusions would benefit from stating explicitly that the numerical results are based on 10 configurations with 32 noise hits, since that context is currently given only in the body of §3.

Circularity Check

0 steps flagged · score 0.0 of 10

No significant circularity: the variance-reduction claims are direct measurements, and the only self-cited identity is a standard algebraic relation that is not load-bearing.

full rationale

The paper does not derive its central quantitative claim from fitted inputs; it measures the variance of the standard and split-even estimators on the same ensemble and compares them. The split-even estimator in Eq. (13) is explicitly constructed from the same propagators as the standard estimator at the same computational cost, and the reported 5-25x and 4-10x reductions are empirical comparisons of measured variances, not predictions forced by construction. The identity in Eq. (12), D_l^{-1} - D_c^{-1} = (m_c - m_l) D_l^{-1} D_c^{-1}, is cited to [11], whose authors overlap with the present paper, but this identity is a direct algebraic consequence of the form of the Dirac operator D_q = D_0 + m_q and is independently verifiable; it is not a contested uniqueness theorem or ansatz smuggled in by citation. The paper applies an existing estimator from [10] rather than renaming a known result, and it does not fit parameters and then call them predictions. The acknowledged limitations, such as the low statistics of 10 configurations and the uncertainty in whether the variance reduction propagates to the final form factor, are statistical robustness concerns rather than circularity. The derivation chain is therefore self-contained with respect to its claims.

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

No fitted physical parameters; the only hand-chosen inputs are the intermediate masses. The split-even identity is an established mathematical result. The main assumptions are standard lattice QCD approximations plus the expectation that the variance reduction survives the full amplitude analysis.

free parameters (1)
  • Intermediate charm masses c1-c5 = a m_c = 0.0362, 0.15, 0.25, 0.30, 0.35; M_eta_c = 693(2), 1492(1), 2007(1), 2230(1), 2432(1) MeV
    Chosen intermediate masses for frequency splitting and charm extrapolation; the reported variance reduction for each mass splitting depends on this choice, but the values are not fitted to a final physical result.
assumptions (4)
  • standard math D_l^{-1} - D_c^{-1} = (m_c - m_l) D_l^{-1} D_c^{-1} for Wilson and Domain-Wall fermions
    Used in Eq. (12) to rewrite loop differences; proof referenced to [11].
  • domain assumption zMöbius action tuned to approximate Möbius action with smaller L_s, so the resulting bias is negligible
    Section 3: correction expected small and neglected in this study.
  • domain assumption Disconnected diagrams are color and SU(3) flavor suppressed
    Section 1.2: neglected in this study.
  • domain assumption Variance reduction on integrated correlators propagates to the final form factor
    Conclusions: expected but not certain; the paper explicitly flags this.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Split-even approach to the rare kaon decay $K \to \pi \ell^+ \ell^-$." pith.science (2026). https://pith.science/paper/FT6L2C6Y

@misc{pith2026250118358,
  author       = {Pith},
  title        = {Pith review of: Split-even approach to the rare kaon decay $K \to \pi \ell^+ \ell^-$},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/FT6L2C6Y}},
  note         = {Machine review of arXiv:2501.18358}
}
read the original abstract

In recent years the rare kaon decay has been computed directly at the physical point. However, this calculation is currently limited by stochastic noise stemming from a light and charm quark loop GIM subtraction. The split-even approach is an alternative estimator for such loop differences, and has shown a large variance reduction in certain quantities. We present an investigation into the use of the split-even estimator in the calculation of the rare kaon decay.

Figures

Figures reproduced from arXiv: 2501.18358 by the authors.

Figure 1
Figure 1. Representatives of the different classes of diagrams in the 4-point correlation function. These are the Non-Loop (top left), Loop (top right), Loop-Insertion (bottom) respectively. The Non-Loop and Loop topologies also have additional diagrams with the current inserted on different legs. There are additional disconnected diagrams that are not considered in this study. 1.2 Existing lattice result The framework for co… view at source ↗
Figure 2
Figure 2. Comparison between the 4-point correlation function in eq. (6) for different charm masses using the standard (left) and split-even (right) estimators. The non-loop diagram is shown for reference. and strange quark masses and an inverse lattice spacing of 𝑎 −1 = 1730 MeV. All sea quarks and the strange valance quark are simulated using the Möbius Domain-Wall action, while the light and charm valance quarks are simula… view at source ↗
Figure 3
Figure 3. The 4-point correlation function (left) and the corresponding variance (right), broken down into the NL, L, and LI diagrams. to the variance at a level below that of the L diagrams. Since the LI diagrams are the most expensive to compute, we may therefore benefit greatly by reducing the computational effort spent on these diagrams, for example by performing fewer time translations and/or reducing the number of noise… view at source ↗
Figures from the paper (1 more)
Figure 4
Figure 4. Figure 4: Variance of the two integrated correlators (top and bottom) using the standard (left) and split-even (right) estimators. nature of this study study, measurements will need to have to be made on additional configurations, along with additional noise hits. While it is ex…

Discussion (0). Continue with ORCID to comment.

Forward citations

Cited by 1 Pith paper

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score. Full citation record

  1. From scattering towards multi-hadron weak decays

    hep-lat 2025-01 unverdicted novelty 1.0 of 10

    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

15 extracted references · 4 canonical work pages · cited by 1 Pith paper

  1. [1]

    Buchalla, A.J

    G. Buchalla, A.J. Buras and M.E. Lautenbacher,Weak decays beyond leading logarithms, Rev. Mod. Phys.68(1996) 1125 [hep-ph/9512380]

    arXiv 1996

  2. [2]

    Rare Kaon Decays on the Lattice

    G. Isidori, G. Martinelli and P. Turchetti,Rare kaon decays on the lattice, Physics Letters B 633 (2006) 75 [hep-lat/0506026]

    work page Pith review arXiv 2006

  3. [3]

    D’Ambrosio, G

    G. D’Ambrosio, G. Ecker, G. Isidori and J. Portoles,The Decays𝐾→𝜋𝑙+𝑙− beyond leading order in the chiral expansion, JHEP 08(1998) 004 [hep-ph/9808289]

    arXiv 1998

  4. [4]

    D’Ambrosio, D

    G. D’Ambrosio, D. Greynat and M. Knecht,Matching long and short distances at order O(𝛼𝑠) in the form factors for𝐾→𝜋ℓ+ℓ−,Phys. Lett. B797 (2019) 134891 [1906.03046]

    arXiv 2019

  5. [5]

    NA62collaboration, A measurement of the𝐾+→𝜋+𝜇+𝜇− decay,JHEP 11(2022) 011 [2209.05076]

    arXiv 2022

  6. [6]

    NA48/1collaboration, Observation of the rare decay𝐾𝑆→𝜋0𝜇+𝜇−,Phys. Lett. B599 (2004) 197 [hep-ex/0409011]

    work page Pith review arXiv 2004

  7. [7]

    RBC, UKQCD collaboration, Prospects for a lattice computation of rare kaon decay amplitudes: 𝐾→𝜋ℓ+ℓ− decays, Phys. Rev. D92(2015) 094512 [1507.03094]

    arXiv 2015

  8. [8]

    N.H.Christ,X.Feng, A.Juttner,A.Lawson, A.PortelliandC.T.Sachrajda, Firstexploratory calculation of the long-distance contributions to the rare kaon decays𝐾→𝜋ℓ+ℓ−, Phys. Rev. D94(2016) 114516 [1608.07585]

    arXiv 2016

Show all 15 references
  1. [9]

    RBC, UKQCD collaboration, Simulating rare kaon decays𝐾+→𝜋+ℓ+ℓ− using domain wall lattice QCD with physical light quark masses, Phys. Rev. D107 (2023) L011503 [2202.08795]

    2023 arXiv

  2. [10]

    Giusti, T

    L. Giusti, T. Harris, A. Nada and S. Schaefer,Frequency-splitting estimators of single-propagator traces,Eur. Phys. J. C79(2019) 586 [1903.10447]

    2019 arXiv

  3. [11]

    Harris, V

    T. Harris, V. Gülpers, A. Portelli and J. Richings,Efficiently unquenching QCD+QED at O(𝛼),PoS LATTICE2022(2023) 013 [2301.03995]

    2023 arXiv

  4. [12]

    RBC, UKQCD collaboration, Domain wall QCD with physical quark masses, Phys. Rev. D 93 (2016) 074505 [1411.7017]

    2016 arXiv

  5. [13]

    Mcglynn,Algorithmic improvements for weak coupling simulations of domain wall fermions, PoS LATTICE2015(2016) 019

    G. Mcglynn,Algorithmic improvements for weak coupling simulations of domain wall fermions, PoS LATTICE2015(2016) 019

    2016

  6. [14]

    T. Blum, T. Izubuchi and E. Shintani,New class of variance-reduction techniques using lattice symmetries, Phys. Rev. D88 (2013) 094503 [1208.4349]

    2013 arXiv

  7. [15]

    Boyle, A

    P. Boyle, A. Juttner, M.K. Marinkovic, F. Sanfilippo, M. Spraggs and J.T. Tsang,An exploratory study of heavy domain wall fermions on the lattice, JHEP 04(2016) 037 [1602.04118]. 9

    2016 arXiv

Pith tools

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