Edward Witten's 1977 paper 'Short Distance Analysis of Weak Interactions' combined the operator product expansion with renormalization-group improvement to compute how weak amplitudes - nonleptonic and rare kaon decays - depend on the W-boson mass MWM_W and on heavy-quark masses, and it stated precisely what the framework could not yet deliver: an all-orders rigorous foundation for products of several weak currents, the nonperturbative hadronic matrix elements and finite matching constants needed for absolute predictions, control of renormalization-scheme dependence, and the behavior of mass-matching functions in confining theories. Most of these items have since been closed: two- and three-loop anomalous dimensions and one-loop matching are complete for ∣ΔS∣=1|\Delta S|=1 and ∣ΔS∣=2|\Delta S|=2, general effective-field-theory running covers all chiral structures including right-handed currents, and the mass-scheme worry was resolved - in the negative - by renormalon theory. Two precise cores remain open: (1) an all-orders-in-perturbation-theory operator-product-expansion theorem for multi-current products, including bilocal time-ordered operators, in asymptotically free non-Abelian gauge theories (proved only for φ44\varphi_4^4-type scalar theories); and (2) first-principles nonperturbative determinations of the long-distance-dominated rare-kaon amplitudes, such as the two-photon contribution to KL→μ+μ−K_L\to\mu^+\mu^- and the K→πγ∗γ∗K\to\pi\gamma^*\gamma^* form factors, which today rest on dispersive estimates, an exploratory lattice bound, and preliminary coarse-lattice computations of the two-photon amplitude (2024).

Background

Nonleptonic and rare decays of kaons probe the weak interaction at two widely separated mass scales: the heavy mediators (the WW boson of mass MW≈80M_W\approx 80 GeV, and any heavy quarks) and the light hadronic scale of a few hundred MeV at which the observed mesons live. Amplitudes therefore mix perturbatively computable short-distance physics with nonperturbative long-distance physics. The standard separation tool is the effective Hamiltonian: heavy fields are integrated out, leaving a sum Heff=∑iCi(μ) QiH_{\rm eff}=\sum_i C_i(\mu)\,Q_i over local four-quark operators QiQ_i (which change the total strangeness ΔS\Delta S by one unit in nonleptonic kaon decays, or by two units in kaon mixing) with short-distance coefficient functions Ci(μ)C_i(\mu) (Wilson coefficients), so that an amplitude factorizes as ∑iCi(μ)⟨Qi(μ)⟩\sum_i C_i(\mu)\langle Q_i(\mu)\rangle between hadronic states. The Wilson coefficients are computed in perturbation theory using the operator product expansion (OPE) - the statement that a product of local operators at short relative distance expands into a series of local operators with coefficient functions depending on the separation - and the large logarithms αsln⁡(MW/μ)\alpha_s\ln(M_W/\mu) generated by the scale hierarchy are summed by the renormalization group (RG), the framework relating a theory at one momentum scale to the same theory at another. When several operators mix under scaling, their scale dependence is governed by an anomalous-dimension matrix; the finite constants relating operator normalizations or mass parameters between theories with and without a heavy field are called matching constants, with heavy-particle decoupling - the statement that very heavy fields affect low-energy physics only through such finite corrections - dating to Appelquist and Carazzone (1975).

Witten (1977) gave the first systematic, model-independent formulation of this short-distance analysis for weak amplitudes. He treated products of weak currents, including the bilocal time-ordered operator ∫d4y T{P(0)P(y)}\int d^4y\,T\{P(0)P(y)\} that represents two-WW intermediate states, derived RG-improved asymptotic formulae for the dependence of amplitudes on MWM_W and on a heavy-quark mass mp′m_{p'}, and applied them to channels such as KL→μ+μ−K_L\to\mu^+\mu^- and K+→π+e+e−K^+\to\pi^+e^+e^-, where the GIM mechanism (the Glashow-Iliopoulos-Maiani cancellation: the vanishing of leading flavor-changing neutral-current contributions once the charm quark and its weak-doublet partner are included) removes the leading terms and exposes the residual heavy-mass dependence. The characteristic results are power laws with fractional exponents generated by one-loop running of the strong coupling, for example a two-WW contribution to KL→μ+μ−K_L\to\mu^+\mu^- proportional to 25 x24/25(1−x1/25)25\,x^{24/25}(1-x^{1/25}) with x=αˉ(MW)/αˉ(mp′)x=\bar\alpha(M_W)/\bar\alpha(m_{p'}), where αˉ\bar\alpha is the running strong coupling. Witten emphasized that powers such as (gˉ2)1/25(\bar g^2)^{1/25} are not small for real-world couplings, so the asymptotic formulae are qualitative rather than quantitative in nature.

The same paper listed what the framework lacked. First, the power-counting arguments rest on Weinberg's theorem (Weinberg 1960), which controls the large-momentum asymptotics of convergent Euclidean Feynman integrals, but no all-orders proof for multi-current operator products is given. Second, the absolute amplitudes require nonperturbative matrix elements of the local operators, not computable in perturbation theory. Third, the finite matching functions (such as B(1,gˉ(MW))B(1,\bar g(M_W)), EE, GG, and S(1,gˉ(mp′))S(1,\bar g(m_{p'})) in the paper) and higher-loop anomalous dimensions were undetermined and scheme-dependent, the renormalization scheme being the prescription for splitting loop integrals into divergent and finite parts. Fourth, right-handed-current extensions were only sketched. Fifth, the assumed benign behavior of the mass-matching function S(1,g)S(1,g) as the coupling g→0g\to 0 in a confining theory was unproven.

The subsequent history is as follows: the two-loop program is complete (Buras-Jamin-Lautenbacher-Weisz 1992; Herrlich-Nierste 1994, 1996; reviewed in Buchalla-Buras-Lautenbacher 1996), next-to-next-to-leading-order (NNLO) corrections followed (Gorbahn-Haisch 2005; Brod-Gorbahn 2012; Brod-Gorbahn-Stamou 2020), general effective-field-theory anomalous dimensions now cover all chiral structures including right-handed currents (Jenkins-Manohar-Stoffer 2018), Caswell-Kennedy (1983) rigorized the underlying Feynman-integral asymptotics, the Hollands-Kopper school proved OPE convergence for perturbative scalar φ44\varphi_4^4 theory (the four-dimensional scalar field theory with quartic self-interaction; Hollands-Kopper 2012; Holland-Hollands 2013; Holland-Hollands-Kopper, arXiv:1411.1785) but not for gauge theories, renormalon theory - the factorial divergence of perturbative series caused by integration over momenta far below the scale of interest, which induces intrinsic ambiguities of order the nonperturbative scale - settled the mass-scheme question in the negative (Beneke-Braun 1995; Beneke 2021), and lattice QCD - quantum chromodynamics, the strong-interaction gauge theory, regularized on a discrete space-time grid so that the path integral can be simulated numerically - delivered direct K→ππK\to\pi\pi amplitudes (Blum et al. 2011; Abbott et al. 2020) while the long-distance-dominated rare modes remain without first-principles determinations (Boyle et al. 2023 obtained only a bound, and the first lattice computations of the two-photon amplitude, in 2024, are still preliminary and coarse, Chao 2024; Chao and Hu 2024; the best KL→ℓ+ℓ−K_L\to\ell^+\ell^- prediction is dispersive, Hoferichter-Hoid-Ruiz de Elvira 2024; the K→πℓ+ℓ−K\to\pi\ell^+\ell^- form-factor computation remains limited by stochastic noise, Hodgson, Gülpers, Hill, and Portelli 2025; and the 2026 lattice calculation of structure-dependent rare-kaon form factors, Di Palma et al., treats the different K→ℓν ℓ+ℓ−K\to\ell\nu\,\ell^+\ell^- channels).

Problem Statement

Witten's short-distance analysis of weak interactions (1977) left open a set of unresolved issues. After the resolutions summarized above, the surviving problem is to complete the two remaining ingredients of that program, which together would make the short-distance analysis of weak amplitudes both rigorously founded and quantitatively reliable.

(1) Rigorous all-orders foundation. Establish - or refute - the multi-current operator product expansion used in the short-distance analysis of weak amplitudes, to all orders in perturbation theory, in an asymptotically free non-Abelian gauge theory (a theory such as QCD - quantum chromodynamics, the strong-interaction gauge theory, whose coupling weakens logarithmically at short distances - or a controlled version retaining its gauge and GIM structures): the expansion of products of two or more weak currents, including bilocal time-ordered operators such as ∫d4y T{P(0)P(y)}\int d^4y\,T\{P(0)P(y)\}, where PP is the weak-current product of the second-order weak interaction whose integration over the relative coordinate yy represents two-WW intermediate states and TT denotes time ordering, into local operators with renormalization-group-improved coefficient functions, with explicit remainder bounds of the kind proved for perturbative scalar φ44\varphi_4^4 theory (the four-dimensional scalar field theory with quartic self-interaction) by Hollands and Kopper. This entails generalizing Weinberg's large-momentum theorem (which fixes the asymptotic behavior of Feynman integrals when some integration momenta grow large) and the existing scalar-theory convergence proofs to gauge theories, controlling gauge fixing, the Ward/Slavnov-Taylor identities (the identities encoding gauge invariance in the gauge-fixed quantization), and exceptional-momentum configurations (regions where subsets of external momenta become simultaneously soft).

(2) First-principles long-distance input. Determine from first principles - by lattice QCD or another fully nonperturbative method with no phenomenological tuning to the target observable - the long-distance-dominated contributions to rare kaon amplitudes on which the short-distance framework is silent: the two-photon intermediate-state contribution to KL→μ+μ−K_L\to\mu^+\mu^- and the γ∗γ∗\gamma^*\gamma^* form factors governing K+→π+ℓ+ℓ−K^+\to\pi^+\ell^+\ell^- (where γ∗\gamma^* denotes a virtual photon and ℓ\ell a charged lepton), with controlled uncertainties sufficient for quantitative Standard-Model predictions.

References

  1. Steven Weinberg, 'High-energy behavior in quantum field theory', Physical Review 118 (1960) 838, DOI 10.1103/PhysRev.118.838, https://doi.org/10.1103/PhysRev.118.838
  2. Thomas Appelquist and J. Carazzone, 'Infrared singularities and massive fields', Physical Review D 11 (1975) 2856, DOI 10.1103/PhysRevD.11.2856, https://doi.org/10.1103/PhysRevD.11.2856
  3. Edward Witten, 'Short distance analysis of weak interactions', Nuclear Physics B 122 (1977) 109, DOI 10.1016/0550-3213(77)90428-x, https://doi.org/10.1016/0550-3213(77)90428-x
  4. William E. Caswell and A. D. Kennedy, 'Asymptotic behavior of Feynman integrals: convergent integrals', Physical Review D 28 (1983) 3073, DOI 10.1103/PhysRevD.28.3073, https://doi.org/10.1103/PhysRevD.28.3073
  5. Andrzej J. Buras, Matthias Jamin, Markus E. Lautenbacher, and Peter H. Weisz, 'Effective Hamiltonians for ΔS=1\Delta S = 1 and ΔB=1\Delta B = 1 non-leptonic decays beyond the leading logarithmic approximation', Nuclear Physics B 370 (1992) 69, DOI 10.1016/0550-3213(92)90345-C, https://doi.org/10.1016/0550-3213(92)90345-C
  6. Stefan Herrlich and Ulrich Nierste, 'Enhancement of the KLK_L-KSK_S mass difference by short-distance QCD corrections beyond leading logarithms', Nuclear Physics B 419 (1994) 292, DOI 10.1016/0550-3213(94)90044-2, arXiv:hep-ph/9310311, https://arxiv.org/abs/hep-ph/9310311
  7. M. Beneke and V. M. Braun, 'Naive nonabelianization and resummation of fermion bubble chains', Physics Letters B 348 (1995) 513, DOI 10.1016/0370-2693(95)00184-M, arXiv:hep-ph/9411229, https://arxiv.org/abs/hep-ph/9411229
  8. Stefan Herrlich and Ulrich Nierste, 'The complete ∣ΔS∣=2|\Delta S| = 2 Hamiltonian in the next-to-leading order', Nuclear Physics B 476 (1996) 27, DOI 10.1016/0550-3213(96)00324-0, arXiv:hep-ph/9604330, https://arxiv.org/abs/hep-ph/9604330
  9. Gerhard Buchalla, Andrzej J. Buras, and Markus E. Lautenbacher, 'Weak decays beyond leading logarithms', Reviews of Modern Physics 68 (1996) 1125, DOI 10.1103/RevModPhys.68.1125, arXiv:hep-ph/9512380, https://arxiv.org/abs/hep-ph/9512380
  10. Martin Gorbahn and Ulrich Haisch, 'Effective Hamiltonian for non-leptonic ∣ΔF∣=1|\Delta F| = 1 decays at NNLO in QCD', Nuclear Physics B 713 (2005) 291, DOI 10.1016/j.nuclphysb.2005.01.047, arXiv:hep-ph/0411071, https://arxiv.org/abs/hep-ph/0411071
  11. T. Blum, P. A. Boyle, N. H. Christ, N. Garron, E. Goode, T. Izubuchi, C. Lehner, Q. Liu, R. D. Mawhinney, C. T. Sachrajda, A. Soni, C. Sturm, H. Yin, and R. Zhou (RBC/UKQCD), 'K to ππ\pi\pi decay amplitudes from lattice QCD', Physical Review D 84 (2011) 114503, DOI 10.1103/PhysRevD.84.114503, arXiv:1106.2714, https://arxiv.org/abs/1106.2714
  12. Joachim Brod and Martin Gorbahn, 'Next-to-next-to-leading-order charm-quark contribution to the CP violation parameter ϵK\epsilon_K and ΔMK\Delta M_K', Physical Review Letters 108 (2012) 121801, DOI 10.1103/PhysRevLett.108.121801, arXiv:1108.2036, https://arxiv.org/abs/1108.2036
  13. Stefan Hollands and Christoph Kopper, 'The operator product expansion converges in perturbative field theory', Communications in Mathematical Physics 313 (2012) 257, DOI 10.1007/s00220-012-1457-4, arXiv:1105.3375, https://arxiv.org/abs/1105.3375
  14. Jan Holland and Stefan Hollands, 'Operator product expansion algebra', Journal of Mathematical Physics 54 (2013) 072302, DOI 10.1063/1.4810780, arXiv:1205.4904, https://arxiv.org/abs/1205.4904
  15. Jan Holland, Stefan Hollands, and Christoph Kopper, 'The operator product expansion converges in massless φ44\varphi_4^4-theory', Communications in Mathematical Physics 342 (2016) 385, DOI 10.1007/s00220-015-2486-6, arXiv:1411.1785, https://arxiv.org/abs/1411.1785
  16. Elizabeth E. Jenkins, Aneesh V. Manohar, and Peter Stoffer, 'Low-energy effective field theory below the electroweak scale: anomalous dimensions', Journal of High Energy Physics 01 (2018) 084, DOI 10.1007/JHEP01(2018)084, arXiv:1711.05270, https://arxiv.org/abs/1711.05270
  17. R. Abbott, T. Blum, P. A. Boyle, M. Bruno, N. H. Christ, D. Hoying, C. Jung, C. Kelly, C. Lehner, R. D. Mawhinney, D. J. Murphy, C. T. Sachrajda, A. Soni, M. Tomii, and T. Wang (RBC/UKQCD), 'Direct CP violation and the ΔI=1/2\Delta I = 1/2 rule in K→ππK\to\pi\pi decay from the Standard Model', Physical Review D 102 (2020) 054509, DOI 10.1103/PhysRevD.102.054509, arXiv:2004.09440, https://arxiv.org/abs/2004.09440
  18. Joachim Brod, Martin Gorbahn, and Emmanuel Stamou, 'Standard-Model prediction of ϵK\epsilon_K with manifest quark-mixing unitarity', Physical Review Letters 125 (2020) 171803, DOI 10.1103/PhysRevLett.125.171803, arXiv:1911.06822, https://arxiv.org/abs/1911.06822
  19. Martin Beneke, 'Pole mass renormalon and its ramifications', European Physical Journal Special Topics 230 (2021) 2565, DOI 10.1140/epjs/s11734-021-00268-w, arXiv:2108.04861, https://arxiv.org/abs/2108.04861
  20. P. A. Boyle, F. Erben, J. M. Flynn, V. Gülpers, R. C. Hill, R. Hodgson, A. Jüttner, F. Ó hÓgáin, A. Portelli, and C. T. Sachrajda (RBC/UKQCD), 'Simulating rare kaon decays K+→π+ℓ+ℓ−K^{+}\to\pi^{+}\ell^{+}\ell^{-} using domain wall lattice QCD with physical light quark masses', Physical Review D 107 (2023) L011503, DOI 10.1103/PhysRevD.107.L011503, arXiv:2202.08795, https://arxiv.org/abs/2202.08795
  21. Martin Hoferichter, Bai-Long Hoid, and Jacobo Ruiz de Elvira, 'Improved Standard-Model prediction for KL→ℓ+ℓ−K_L\to\ell^+\ell^-', Journal of High Energy Physics 04 (2024) 071, DOI 10.1007/JHEP04(2024)071, arXiv:2310.17689, https://arxiv.org/abs/2310.17689
  22. En-Hung Chao, 'Progress on KL→μ+μ−K_L\to\mu^+\mu^- from lattice QCD', arXiv:2403.18885 (2024), https://arxiv.org/abs/2403.18885
  23. En-Hung Chao and Ceran Hu, 'Two-photon contribution to the KL→μ+μ−K_L\to\mu^+\mu^- decay amplitude from a 1/a≈11/a\approx 1 GeV lattice', arXiv:2411.00710 (2024), https://arxiv.org/abs/2411.00710
  24. Raoul Hodgson, Vera Gülpers, Ryan Hill, and Antonin Portelli, 'Split-even approach to the rare kaon decay K→πℓ+ℓ−K\to\pi\ell^+\ell^-', arXiv:2501.18358 (2025), https://arxiv.org/abs/2501.18358
  25. R. Di Palma, R. Frezzotti, G. Gagliardi, V. Lubicz, G. Martinelli, C. T. Sachrajda, F. Sanfilippo, S. Simula, and N. Tantalo, 'Rare kaon decays K−→ℓ−νˉℓ ℓ′+ℓ′−K^-\to\ell^-\bar\nu_\ell\,\ell'^+\ell'^-: Standard Model predictions from lattice QCD', arXiv:2605.22727 (2026), https://arxiv.org/abs/2605.22727