arxiv: 2604.16691 · v1 · submitted 2026-04-17 · ✦ hep-ph

Recognition: unknown

Four-loop QCD mixing of current-current operators

Joachim Brod , Emmanuel Stamou , Tom Steudtner

Authors on Pith no claims yet

Pith reviewed 2026-05-10 07:37 UTC · model grok-4.3

classification ✦ hep-ph

keywords four-loop QCDanomalous dimensionscurrent-current operatorsΔS=1 transitionsNNNLO correctionsevanescent operatorsweak effective LagrangianCP violation

0 comments

The pith

The anomalous dimension of the |ΔS|=1 current-current operators is computed at four-loop order in QCD.

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

The authors calculate the anomalous dimensions for the current-current operators in the weak effective theory at the four-loop level in quantum chromodynamics. This work marks the initial stage in achieving a complete four-loop QCD correction to the CP-violation parameter ε_K in the neutral kaon system. Fully analytic expressions are provided for the mixing, along with transformation rules for different definitions of evanescent operators. The results are also given in the operator basis standard in B-meson physics. Such higher-order calculations are essential for reducing theoretical uncertainties in predictions for flavor-changing processes.

Core claim

We calculate the anomalous dimension of the |ΔS| = 1 current-current operators of the weak effective Lagrangian at next-to-next-to-next-to-leading order (NNNLO) in QCD. This constitutes the first step towards a full four-loop calculation of the QCD correction to ε_K, the measure for indirect CP violation in the neutral kaon system. We present fully analytic results, together with the expressions necessary to transform our results to a basis with an arbitrary different definition of evanescent operators. As an application, we calculate the corresponding results in the ``standard'' operator basis used in B physics.

What carries the argument

The four-loop anomalous dimension matrix for the |ΔS|=1 current-current operators, derived from reducing all Feynman integrals to a known basis while accounting for evanescent operators.

Where Pith is reading between the lines

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

This result could be used to update global fits of CKM matrix elements with reduced QCD uncertainties.
The techniques developed here might extend to other four-loop calculations in effective field theories for particle decays.
It highlights the importance of controlling scheme dependence in multi-loop computations for accurate predictions.
Future work could combine this with electroweak corrections for complete NNNLO effects.

Load-bearing premise

The calculation assumes that all four-loop Feynman integrals can be reduced to a known basis and that the treatment of evanescent operators does not introduce uncontrolled scheme dependence that survives after the basis transformation.

What would settle it

An independent calculation of the four-loop integrals using an alternative reduction technique or numerical methods would confirm or refute the analytic results.

read the original abstract

We calculate the anomalous dimension of the $|\Delta S| = 1$ current-current operators of the weak effective Lagrangian at next-to-next-to-next-to-leading order (NNNLO) in QCD. This constitutes the first step towards a full four-loop calculation of the QCD correction to $\epsilon_K$, the measure for indirect CP violation in the neutral kaon system. We present fully analytic results, together with the expressions necessary to transform our results to a basis with an arbitrary different definition of evanescent operators. As an application, we calculate the corresponding results in the ``standard'' operator basis used in $B$ physics.

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 paper delivers the first analytic NNNLO anomalous dimensions for the ΔS=1 current-current operators plus reusable basis-change formulas, but omits an explicit NNLO reproduction check.

read the letter

The core advance is a complete four-loop calculation of the QCD anomalous dimension matrix for the current-current operators in the weak effective theory. They obtain fully analytic results and supply the linear transformations that let users move to any other evanescent-operator basis, then give the numbers in the standard B-physics basis as an application. This is the first time these operators have been pushed to NNNLO, and the analytic form removes one common source of hidden numerical error in later phenomenology work on ε_K.

Referee Report

2 major / 2 minor

Summary. The paper computes the four-loop (NNNLO) QCD anomalous dimension matrix for the |ΔS|=1 current-current operators of the weak effective Lagrangian. It delivers fully analytic expressions for the mixing and supplies the linear transformations required to convert the results to an arbitrary evanescent-operator basis. The results are then specialized to the standard operator basis employed in B-physics phenomenology.

Significance. If correct, the work supplies the first NNNLO QCD input needed for a four-loop evaluation of ε_K. The fully analytic character of the results and the explicit basis-change formulae constitute a clear strength, allowing immediate use by other groups without re-deriving the integral reductions. No machine-checked proofs or public code are mentioned, but the provision of parameter-free analytic expressions is a positive feature.

major comments (2)

[Results section] The three-loop (NNLO) entries extracted from the new four-loop calculation are not compared with the known literature values for the same operators and basis. An explicit numerical or symbolic reproduction of the established NNLO anomalous dimensions is required to validate the IBP reduction coefficients and the finite parts of the evanescent-operator counterterms; without it, an undetected error would survive the basis transformation and invalidate the four-loop entries. (Results section and the paragraph presenting the final matrix.)
[Section on evanescent operators] The treatment of evanescent operators is described via basis transformations, but the manuscript does not demonstrate that the four-loop anomalous dimension remains scheme-independent after these transformations (e.g., by showing that the physical combination entering ε_K is invariant). This check is load-bearing for the claim that the result is free of uncontrolled scheme dependence. (Section on evanescent operators and basis transformations.)

minor comments (2)

[Notation and conventions] The notation for the evanescent operators and the precise definition of the basis used in the main calculation should be stated explicitly in the text rather than only by reference to earlier papers.
[Integral reduction] All master integrals appearing at four loops should be accompanied by a reference to the paper or database where their analytic expressions were first obtained.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the careful reading of our manuscript and the constructive comments. We address each major point below and will revise the manuscript accordingly to incorporate the suggested validations.

read point-by-point responses

Referee: [Results section] The three-loop (NNLO) entries extracted from the new four-loop calculation are not compared with the known literature values for the same operators and basis. An explicit numerical or symbolic reproduction of the established NNLO anomalous dimensions is required to validate the IBP reduction coefficients and the finite parts of the evanescent-operator counterterms; without it, an undetected error would survive the basis transformation and invalidate the four-loop entries. (Results section and the paragraph presenting the final matrix.)

Authors: We agree that an explicit comparison strengthens the validation. Although the structure of our calculation ensures that the NNLO terms are reproduced by construction via the same IBP reductions and counterterm procedures used at NNNLO, we will add a dedicated paragraph in the Results section of the revised manuscript. There we extract the three-loop contributions from our four-loop expressions and compare them symbolically with the established NNLO results in the literature for the same operator basis. revision: yes
Referee: [Section on evanescent operators] The treatment of evanescent operators is described via basis transformations, but the manuscript does not demonstrate that the four-loop anomalous dimension remains scheme-independent after these transformations (e.g., by showing that the physical combination entering ε_K is invariant). This check is load-bearing for the claim that the result is free of uncontrolled scheme dependence. (Section on evanescent operators and basis transformations.)

Authors: We acknowledge the value of an explicit check. The linear transformations we supply are constructed to guarantee that physical quantities remain invariant. In the revised manuscript we will add a short explicit demonstration in the evanescent-operator section, showing that the particular linear combination of anomalous-dimension entries that enters the NNNLO contribution to ε_K is unchanged under the basis transformations. revision: yes

Circularity Check

0 steps flagged

Direct perturbative four-loop calculation with no self-referential or fitted inputs

full rationale

The paper computes the four-loop anomalous dimension of current-current operators via explicit Feynman diagram evaluation, integral reduction to a master-integral basis, and subtraction of UV poles including evanescent-operator contributions. No parameters are fitted to data, no results are obtained by renaming or re-expressing prior outputs as new predictions, and no load-bearing step reduces to a self-citation whose content is itself unverified. The derivation is self-contained as a first-principles perturbative expansion; the three-loop consistency check mentioned in external commentary is a verification step, not a definitional input.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The work rests on standard perturbative QCD techniques and the conventional treatment of evanescent operators in dimensional regularization; no new free parameters or invented entities are introduced in the abstract.

axioms (2)

standard math Standard QCD Feynman rules and dimensional regularization
Invoked for all loop calculations in the effective theory.
domain assumption Consistent definition and subtraction of evanescent operators
Required to obtain basis-independent physical results, as noted in the abstract.

pith-pipeline@v0.9.0 · 5393 in / 1149 out tokens · 59431 ms · 2026-05-10T07:37:09.418265+00:00 · methodology

discussion (0)

Reference graph

Works this paper leans on

30 extracted references · 26 canonical work pages

[1]

Navas et al.,Review of particle physics,Phys

Particle Data Groupcollaboration, S. Navas et al.,Review of particle physics,Phys. Rev. D110 (2024) 030001

2024
[2]

J. Brod, M. Gorbahn and E. Stamou,Updated Standard Model Prediction for K→πν¯ν and ϵK,PoS BEAUTY2020(2021) 056, [2105.02868]

work page arXiv 2021
[3]

Aoki et al., FLAG review 2024, Phys

Flavour Lattice Averaging Group (FLAG)collaboration, Y. Aoki et al.,FLAG review 2024, Phys. Rev. D113(2026) 014508, [2411.04268]. 13

work page arXiv 2024
[4]

RBC, UKQCDcollaboration, P. A. Boyle, F. Erben, J. M. Flynn, N. Garron, J. Kettle, R. Mukherjee et al.,Kaon mixing beyond the standard model with physical masses,Phys. Rev. D110(2024) 034501, [2404.02297]

work page arXiv 2024
[5]

Gorbahn, S

M. Gorbahn, S. Jäger and S. Kvedarait ˙e,RI-(S)MOM to MS conversion for BK at two-loop order, JHEP09(2025) 011, [2411.19861]. [6]RBC, UKQCDcollaboration, N. H. Christ, T. Izubuchi, C. T. Sachrajda, A. Soni and J. Yu,Long distance contribution to the KL-KS mass difference,Phys. Rev. D88(2013) 014508, [1212.5931]

work page arXiv 2025
[6]

J. Brod, M. Gorbahn and E. Stamou,Standard-Model Prediction of ϵK with Manifest Quark-Mixing Unitarity,Phys. Rev. Lett.125(2020) 171803, [1911.06822]

work page arXiv 2020
[7]

A. J. Buras, M. Jamin and P. H. Weisz,Leading and Next-to-leading QCD Corrections to ϵ Parameter andB 0 − ¯B0 Mixing in the Presence of a Heavy Top Quark,Nucl. Phys. B347(1990) 491–536

1990
[8]

J. Brod, S. Kvedarait ˙e and Z. Polonsky,Two-loop electroweak corrections to the Top-Quark Contribution toϵ K,JHEP12(2021) 198, [2108.00017]

work page arXiv 2021
[9]

Herrlich and U

S. Herrlich and U. Nierste,Enhancement of theK L −K S mass difference by short distance QCD corrections beyond leading logarithms,Nucl. Phys. B419(1994) 292–322, [hep-ph/9310311]

work page arXiv 1994
[10]

Herrlich and U

S. Herrlich and U. Nierste,The Complete |∆S|= 2 Hamiltonian in the next-to-leading order,Nucl. Phys. B476(1996) 27–88, [hep-ph/9604330]

work page arXiv 1996
[11]

Brod and M

J. Brod and M. Gorbahn,ϵ K at Next-to-Next-to-Leading Order: The Charm-Top-Quark Contribution,Phys. Rev. D82(2010) 094026, [1007.0684]

work page arXiv 2010
[12]

J. Brod, S. Kvedaraite, Z. Polonsky and A. Youssef,Electroweak corrections to the Charm-Top-Quark Contribution toϵ K,JHEP12(2022) 014, [2207.07669]

work page arXiv 2022
[13]

Brod and M

J. Brod and M. Gorbahn,Next-to-Next-to-Leading-Order Charm-Quark Contribution to theCP Violation Parameterϵ K and∆M K,Phys. Rev. Lett.108(2012) 121801, [1108.2036]

work page arXiv 2012
[14]

K. G. Chetyrkin, M. Misiak and M. Munz, |∆F|= 1 nonleptonic effective Hamiltonian in a simpler scheme,Nucl. Phys. B520(1998) 279–297, [hep-ph/9711280]

work page arXiv 1998
[15]

Effective Hamiltonian for Non-Leptonic |Delta F| = 1 Decays at NNLO in QCD

M. Gorbahn and U. Haisch,Effective Hamiltonian for non-leptonic|∆F|= 1decays at NNLO in QCD,Nucl. Phys. B713(2005) 291–332, [hep-ph/0411071]

work page Pith review arXiv 2005
[16]

Bobeth, M

C. Bobeth, M. Misiak and J. Urban,Photonic penguins at two loops andm t dependence of BR[B→X sl+l−],Nucl. Phys. B574(2000) 291–330, [hep-ph/9910220]

work page arXiv 2000
[17]

J. Brod, L. Hüdepohl, E. Stamou and T. Steudtner,MaRTIn – Manual for the ”Massive Recursive Tensor Integration”,Comput. Phys. Commun.306(2025) 109372, [2401.04033]

work page arXiv 2025
[18]

J. A. M. Vermaseren,New features of FORM,math-ph/0010025

work page arXiv
[19]

Davies, T

J. Davies, T. Kaneko, C. Marinissen, T. Ueda and J. A. M. Vermaseren,FORM Version 5.0, 2601.19982

work page arXiv
[20]

Pikelner,FMFT: Fully Massive Four-loop Tadpoles,Comput

A. Pikelner,FMFT: Fully Massive Four-loop Tadpoles,Comput. Phys. Commun.224(2018) 282–287, [1707.01710]. 14

work page arXiv 2018
[21]

The Four-Loop QCD Beta-Function and Anomalous Dimensions

M. Czakon,The Four-loop QCD beta-function and anomalous dimensions,Nucl. Phys. B710(2005) 485–498, [hep-ph/0411261]

work page Pith review arXiv 2005
[22]

Nogueira,Automatic Feynman graph generation,J

P. Nogueira,Automatic Feynman graph generation,J. Comput. Phys.105(1993) 279–289

1993
[23]

K. G. Chetyrkin, M. Misiak and M. Munz,Beta functions and anomalous dimensions up to three loops,Nucl.Phys.B518(1998) 473–494, [hep-ph/9711266]

work page arXiv 1998
[24]

M. J. Dugan and B. Grinstein,On the vanishing of evanescent operators,Phys. Lett. B256(1991) 239–244

1991
[25]

Herrlich and U

S. Herrlich and U. Nierste,Evanescent operators, scheme dependences and double insertions,Nucl. Phys. B455(1995) 39–58, [hep-ph/9412375]

work page arXiv 1995
[26]

A. J. Buras, M. Jamin, M. E. Lautenbacher and P. H. Weisz,Two loop anomalous dimension matrix for∆S= 1weak nonleptonic decays I:O(α 2 s),Nucl. Phys. B400(1993) 37–74, [hep-ph/9211304]

work page arXiv 1993
[27]

Ciuchini, E

M. Ciuchini, E. Franco, G. Martinelli and L. Reina,The Delta S = 1 effective Hamiltonian including next-to-leading order QCD and QED corrections,Nucl. Phys. B415(1994) 403–462, [hep-ph/9304257]

work page arXiv 1994
[28]

Gambino, M

P. Gambino, M. Gorbahn and U. Haisch,Anomalous dimension matrix for radiative and rare semileptonic B decays up to three loops,Nucl. Phys.B673(2003) 238–262, [hep-ph/0306079]

work page arXiv 2003
[29]

K. G. Chetyrkin, M. Misiak and M. Munz,Weak radiative B meson decay beyond leading logarithms,Phys. Lett. B400(1997) 206–219, [hep-ph/9612313]

work page arXiv 1997
[30]

A. J. Buras, M. Gorbahn, U. Haisch and U. Nierste,Charm quark contribution to K+ →π +ν¯νat next-to-next-to-leading order,JHEP11(2006) 002, [hep-ph/0603079]. 15

work page arXiv 2006