NLO QCD Renormalization Group Evolution for Non-Leptonic Delta F=2 Transitions in the SMEFT

We present for the first time NLO QCD Renormalization Group (RG) evolution matrices for non-leptonic $\Delta F=2$ transitions in the Standard Model Effective Field Theory (SMEFT). To this end we transform first the known two-loop QCD anomalous dimension matrices (ADMs) of the BSM operators in the so-called BMU basis into the ones in the common Weak Effective Theory (WET) basis (the so-called JMS basis) for which tree-level and one-loop matching to the SMEFT are already known. This allows us subsequently to find the two-loop QCD ADMs for the SMEFT non-leptonic $\Delta F=2$ operators in the Warsaw basis. Having all these ingredients we investigate the impact of these NLO QCD effects on the QCD RG evolution of SMEFT Wilson coefficients for non-leptonic $\Delta F=2$ transitions from the new physics scale $\Lambda$ down to the electroweak scale $\mu_\text{ew}$. The main benefit of these new contributions is that they allow to remove renormalization scheme dependences present both in the one-loop matchings between the WET and SMEFT and also between SMEFT and a chosen UV completion. But the NLO QCD effects, calculated here in the NDR scheme, turn out to be small, in the ballpark of a few percent but larger than one-loop Yukawa top effects when only the $\Delta F=2$ operators are considered. The technology developed in our paper allows to obtain the ADMs in the SMEFT from the ones of the BMU basis also for non-leptonic $\Delta F=1$ decays and the results of this more involved analysis will be presented soon in another publication.