EFFECTIVE-HAMILTONIANS FOR DELTA-S=1 AND DELTA-B=1 NONLEPTONIC DECAYS BEYOND THE LEADING LOGARITHMIC APPROXIMATION

We calculate the effective hamiltonians for DELTA-B = 1 and DELTA-S = 1 decays including next-to-leading QCD corrections. In particular we present the two-loop 6 x 6 anomalous dimension matrices describing the mixing of current-current and QCD penguin operators calculated in the dimensional regularization scheme with anticommuting gamma-5 and in the 't Hooft-Veltman scheme. The renormalization scheme dependences and their cancellation for physical quantities are discussed in detail. We point out that beyond the leading order, matching between coefficients and hadronic matrix elements of local operators both in the renormalization scale-mu and in the renormalization scheme for the operators has to be made. The next-to-leading corrections enhance substantially the Wilson coefficients of QCD penguin operators and these effects are stronger than the ones found for current-current operators. Consequently, the contributions of QCD penguins to B-decays, the DELTA-I = 1/2 rule and epsilon'/epsilon are enlarged. For mu almost-equal-to 1 GeV and LAMBDA(MSBAR) almost-equal-to 300 MeV, these effects amount to 20-30%. In particular we emphasize that for mu < 1 GeV the next-to-leading corrections to the coefficients of penguin operators are large, reaching 40% for mu almost-equal-to 0.8 GeV. For such low values of mu the renormalization group improved perturbation theory cannot be trusted any longer. Moreover the results are very sensitive functions of LAMBDA(MSBAR).