Renormalization group evolution for theΔS=1effective Hamiltonian withNf=2+1

Abstract

We discuss the renormalization group (RG) evolution for the $\ensuremath{\Delta}S=1$ operators in unquenched QCD with ${N}_{f}=3$ (${m}_{u}={m}_{d}={m}_{s}$) or, more generally, ${N}_{f}=2+1$ (${m}_{u}={m}_{d}\ensuremath{\ne}{m}_{s}$) flavors. In particular, we focus on the specific problem of how to treat the singularities which show up only for ${N}_{f}=3$ or ${N}_{f}=2+1$ in the original solution of Buras et al. for the RG evolution matrix at next-to-leading order. On top of the original treatment of Buras et al., we use a new method of analytic continuation to obtain the correct solution in this case. It is free of singularities and can therefore be used in numerical analysis of data sets calculated in lattice QCD.