arxiv:hep-ph/ v3 20 Feb 1998

Transcription

1 ZU-TH 24/96 TUM-HEP-263/96 MPI/PhT/ hep-ph/ December 1996 Weak Radiative B-Meson Decay Beyond Leading Logarithms arxiv:hep-ph/ v3 20 Feb 1998 Konstantin Chetyrkin 1,, Miko laj Misiak 2, and Manfred Münz 3 1 Max-Planck-Institut für Physik, Werner-Heisenberg-Institut, D München, Germany 2 Institut für Theoretische Physik der Universität Zürich, CH-8057 Zürich, Switzerland 3 Physik Department, Technische Universität München, D Garching, Germany Abstract We present our results for three-loop anomalous dimensions necessary in analyzing B X s γ decay at the next-to-leading order in QCD. We combine them with other recently calculated contributions, obtaining a practically complete next-to-leading order prediction for the branching ratio B(B X s γ) = (3.28 ± 0.33) The uncertainty is more than twice smaller than in the previously available leading order theoretical result. The Standard Model prediction remains in agreement with the CLEO measurement at the 2σ level. 1 This work was partially supported by INTAS under Contract INTAS Supported in part by Schweizerischer Nationalfonds and Polish Commitee for Scientific Research (grant 2 P03B , ). 3 Supported in part by Deutsches Bundesministerium für Bildung und Forschung under contract 06 TM 743. Permanent address: Institute of Nuclear Research, Russian Academy of Sciences, Moscow , Russia. Permanent address: Institute of Theoretical Physics, Warsaw University, Hoża 69, 00-6 Warsaw, Poland.

2 1. Weak radiative B-meson decay is known to be a very sensitive probe of new physics [1, 2]. Heavy Quark Effective Theory tells us that inclusive B-meson decay rate into charmless hadrons and the photon is well approximated by the corresponding partonic decay rate Γ(B X s γ) Γ(b X s γ). (1) The accuracy of this approximation is expected to be better than 10% [3]. The inclusive branching ratio B(b X s γ) was extracted two years ago from CLEO measurements of weak radiative B-meson decay. It amounts to [4] B(b X s γ) = (2.32 ± 0.57 ± 0.35) (2) In the forthcoming few years, much more precise measurements of B(B X s γ) are expected from the upgraded CLEO detector, as well as from the B-factories at SLAC and KEK. Acquiring experimental accuracy of below 10% is conceivable. Thus, the goal on the theoretical side is to reach the same accuracy in perturbative calculations of b X s γ decay rate. This requires performing a complete calculation at the next-to-leading order (NLO) in QCD [1, 5]. All the next-to-leading contributions to b X s γ are collected for the first time in the present letter. 2 One of the most complex ingredients of the NLO calculation consists in finding threeloop anomalous dimensions in the effective theory used for resummation of large logarithms ln(mw 2 /m2 b ). This is the main new result presented in this letter. However, we have recalculated all the previously found anomalous dimensions, too. The remaining ingredients of the NLO calculation are already present in the literature [6] [17]. We use the published results for the two-loop matrix elements calculated by Greub, Hurth and Wyler [9], Bremsstrahlung corrections obtained by Ali, Greub [15, 16] and Pott [17], and the matching conditions found by Adel and Yao [13]. The following three sections of the present paper are devoted to presenting the complete NLO formulae for b X s γ. Next, we analyze them numerically. In the end, we discuss predictions for the branching ratio B(B X s γ) including a discussion of the relevant uncertainties. 2 except for the unknown but negligible two-loop matrix elements of the penguin operators which we denote further by P 3,...,P 6. 1

3 2. The analysis of b X s γ decay begins with introducing an effective hamiltonian H eff = 4G F Vts V tb 2 8 C i (µ)p i (µ) (3) i=1 where V ij are elements of the CKM matrix, P i (µ) are the relevant operators and C i (µ) are the corresponding Wilson coefficients. The complete set of physical operators necessary for b X s γ decay is the following: P 1 = ( s L γ µ T a c L )( c L γ µ T a b L ) P 2 = ( s L γ µ c L )( c L γ µ b L ) P 3 = ( s L γ µ b L ) q( qγ µ q) P 4 = ( s L γ µ T a b L ) q( qγ µ T a q) P 5 = ( s L γ µ1 γ µ2 γ µ3 b L ) q( qγ µ 1 q) P 6 = ( s L γ µ1 γ µ2 γ µ3 T a b L ) q( qγ µ 1 T a q) P 7 = P 8 = e m 16π 2 b ( s L σ µν b R )F µν g m 16π 2 b ( s L σ µν T a b R )G a µν where T a stand for SU(3) color generators. The small CKM matrix element V ub as well as the s-quark mass are neglected in the present paper. Our basis of four-quark operators is somewhat different than the standard one used in refs. [1, 6, 18], although the two bases are physically equivalent. We would like to stress that in our basis, Dirac traces containing γ 5 do not arise in effective theory calculations performed at the leading order in the Fermi coupling G F (but to all orders in QCD). This allows us to consistently use fully anticommuting γ 5 in dimensional regularization, which greatly simplifies multiloop calculations. Our choice of color structures is dictated only by convenience in computer algebra applications. Throughout the paper, we use the MS scheme with fully anticommuting γ 5. Such a scheme is not uniquely defined until one chooses a specific form for the so-called evanescent operators which algebraically vanish in four dimensions [19, 20]. In our three-loop calculation, eight evanescent operators were necessary. We list them explicitly in Appendix A. Resummation of large logarithms ln(mw 2 /m2 b ) is achieved by evolving the coefficients C i (µ) from µ = M W to µ = µ b m b according to the Renormalization Group Equations (RGE). Instead of the original coefficients C i (µ), it is convenient to use certain linear combinations of them, called the effective coefficients [1, 9] C eff i (µ) = C i (µ) for i = 1,..., 6 C 7 (µ) + 6 i=1 y i C i (µ) for i = 7 C 8 (µ) + 6 i=1 z i C i (µ) for i = 8. 2 (4) (5)

4 The numbers y i and z i are defined so that the leading-order b sγ and b s gluon matrix elements of the effective hamiltonian are proportional to the leading-order terms in C eff 7 and C eff 8, respectively [1]. In the MS scheme with fully anticommuting γ 5, we have y = (0, 0, 1, 3 4, 9 20, 3 80) and z = (0, 0, 1, 9 1, 20, ).3 The leading-order contributions to the effective coefficients C eff i (µ) are regularization- and renormalization-scheme independent, which would not be true for the original coefficients C 7 (µ) and C 8 (µ). The effective coefficients evolve according to their RGE µ d dµ Ceff i (µ) = C eff j (µ)γ eff ji (µ) (6) driven by the anomalous dimension matrix ˆγ eff (µ). One expands this matrix perturbatively as follows: ˆγ eff (µ) = α s(µ) 4π ˆγ(0)eff + α2 s (µ) (4π) 2 ˆγ(1)eff +... (7) The matrix ˆγ (0)eff is renormalization-scheme independent, while ˆγ (1)eff is not. In the MS scheme with fully anticommuting γ 5 (and with our choice of evanescent operators) we obtain ˆγ (0)eff = (8) and ˆγ (1)eff = (9) 3 These numbers are different than in section 4 of ref. [1] because we use a different basis of four-quark operators here. 3

5 The matrix ˆγ (1)eff is the main new result of this paper. Its evaluation required performing three-loop renormalization of the effective theory (3). Details of this calculation will be given in a forthcoming publication [21]. We expand the coefficients in powers of α s as follows: C eff i (µ) = C (0)eff i (µ) + α s(µ) 4π C(1)eff i (µ) +... (10) At µ = M W, the values of the coefficients are found by matching the effective theory amplitudes with the full Standard Model ones. The well-known leading-order results read [1, 6, 22] where C (0)eff i (M W ) = C (0) i (M W ) = x = m2 t,ms (M W) M 2 W 0 for i = 1, 3, 4, 5, 6 1 for i = 2 3x 3 2x 2 ln x + 8x3 5x 2 +7x for i = 7 4(x 1) 4 24(x 1) 3 3x 2 4(x 1) 4 ln x + x3 +5x 2 +2x 8(x 1) 3 for i = 8, = m2 t,pole M 2 W ( )24 αs (M W ) ( α s (m t ) α s (m t ) 3 π ) (11). (12) The next-to-leading contributions to the four-quark operator coefficients at µ = M W can be found either by transforming the results of ref. [10] to our operator basis, or by a direct computation. We find where C (1)eff i (M W ) = C (1) i (M W ) = 15 for i = 1 0 for i = 2, 3, 5, 6 E(x) 2 3 for i = 4, (13) E(x) = x(18 11x x2 ) + x2 (15 16x + 4x 2 ) ln x 2 ln x. (14) 12(1 x) 3 6(1 x) 4 3 The next-to-leading contributions to C eff 7 (M W ) and C eff 8 (M W ) can be extracted from ref. [13] according to the following equations: where the quantities C (1)eff 7 (M W ) = 8/3 π 4α C(1) O ( E(x) 2 ), (15) s 9 3 C (1)eff 8 (M W ) = 8 π 4α C(1) O ( E(x) 2 ), (16) s 6 3 C (1) O 51 and C (1) O 52 are given in eqn. (5) of ref. [13] with µ = M W. The terms denoted by 7 and 8 stand for shifts one needs to make when passing from the so-called R renormalization scheme used in ref. [13] to the MS scheme used here i = 4π [ α s C (0) i (M W )(x x 2α s π ] x ln x) C(0) i (M W )(x) 4 for i = 7, 8. (17)

6 Only the leading (zeroth-order) term in α s needs to be retained in the above equation. The resulting explicit expressions for C eff 7 (M W ) and C eff 8 (M W ) in the MS scheme read C (1)eff 7 (M W ) = 16x4 122x x 2 ( 8x Li 9(x 1) ) + 6x4 + 46x 3 28x 2 ln 2 x x 3(x 1) x5 588x x x x ln x (x 1) x x x x 436 (18) 486(x 1) 4 C (1)eff 8 (M W ) = 4x4 + 40x x 2 ( + x Li 6(x 1) ) + 17x3 31x 2 ln 2 x x 2(x 1) x x x x x ln x 216(x 1) x x x x (x 1) 4 (19) Having specified the initial conditions at µ = M W, we are ready to write the solution to the Renormalization Group Equations (6) for C eff 7 (µ b ) C (0)eff 7 (µ b ) = η 16 (0) 23 C 7 (M W ) + 8 ( ) (0) η 23 η 23 C 8 (M W ) h i η a i, (20) i=1 C (1)eff 7 (µ b ) = η C (1)eff + ( ) (1)eff η 23 η 23 C 8 (M W ) 7 (M W ) ( η η ( ) (0) η 23 η 23 C 7 (M W ) + where η = α s (M W )/α s (µ b ), and the numbers a i h i are as follows ) η η C (0) 8 (M W ) 8 (e i ηe(x) + f i + g i η)η a i, (21) 14 a i = (, 16, 6, , , , , ) e i = (, 0, 0, , , , ), f i = ( , , , , , , , ) g i = ( , , , , , , , ) h i = (, , , 1, , , , ) As far as the other Wilson coefficients are concerned, only the leading-order contributions i=1 5

7 to them are necessary in the complete NLO analysis of b X s γ. We obtain 4 C (0) 1 (µ b ) = η η 6 23, C (0) 2 (µ b ) = 1 3 η η 23, 6 3 C (0) 3 (µ b ) = 1 27 η η η η η η , 63 C (0) 4 (µ b ) = 1 9 η η η η η η , 21 C (0) 5 (µ b ) = η η η η η η , 126 C (0) 6 (µ b ) = 1 C (0)eff 8 (µ b ) = 36 η η η η η η , 84 ) η ( C (0) 8 (M W ) η η η η (22) In our numerical analysis, the values of α s (µ) in all the above formulae are calculated with use of the NLO expression for the strong coupling constant where β 0 = 23 3 and β 1 = α s (µ) = α s(m Z ) v(µ) [ 1 β 1 α s (M Z ) β 0 4π v(µ) = 1 β 0 α s (M Z ) 2π ln ( ) MZ µ ] ln v(µ), (23) v(µ), (24) As easily can be seen from eqns. (10) (24), the Wilson coefficients at µ = µ b depend on only five parameters: M Z, M W, α s (M Z ), m t,pole and µ b. In our numerical analysis, the first two of them are fixed to M Z = GeV and M W = GeV [23]. For the remaining three, we take α s (M Z ) = ± [24], (25) m t,pole = 175 ± 6 GeV [25] (26) and µ b [2.5 GeV, 10 GeV]. (27) For µ b = 5 GeV and central values of the other parameters, we find α s (µ b ) = 0.212, C (0)eff (µ b ) = ( 0.480, 1.023, , , , , 0.312, 0.148), C (1)eff 7 (µ b ) = = (28) In the intermediate step of the latter equation, we have shown relative importance of two contributions to C (1)eff 7 (µ b ). The first of them originates from terms proportional to C (1)eff 7 (M W ) 4 Analogous expressions in ref. [1] are somewhat different, because we use a different basis of four-quark operators here. 6

8 and C (1)eff 8 (M W ) in eqn. (21), i.e. it is due to the NLO matching. The second of them comes from the remaining terms, i.e. it is driven by the next-to-leading anomalous dimensions. Quite accidentally, the two terms tend to cancel each other. In effect, the total influence of C (1)eff 7 (µ b ) on the b X s γ decay rate does not exceed 6%. (Explicit expressions for the decay rate are given below). However, if the relative sign between the two contributions was positive, then the decay rate would be affected by almost 30%. This observation illustrates the importance of explicit calculation of both the NLO matching conditions and the NLO anomalous dimensions. 3. At this point, we are ready to use the calculated coefficients at µ = µ b m b as an input in the NLO expressions for b X s γ decay rate given by Greub, Hurth and Wyler [9]. However, we make an explicit lower cut on the photon energy in the Bremsstrahlung correction We write the decay rate as follows: E γ > (1 δ)e max γ (1 δ) m b 2. (29) Γ[b X s γ] Eγ>(1 δ)emax γ = Γ[b sγ] + Γ[b sγ gluon] Eγ>(1 δ)emax γ = where and D = C (0)eff 7 (µ b ) + α s(µ b ) 4π = G2 F α em 32π 4 V tsv tb 2 m 3 b,polem 2 b,ms (m b) ( D 2 + A ), (30) { C (1)eff 7 (µ b ) + 8 i=1 A = ( e αs(µ b)lnδ(7+2 lnδ)/3π 1 ) C (0)eff 7 (µ b ) 2 + α s(µ b ) π [ C (0)eff i (µ b ) r i + γ (0)eff i7 ln m ]} b µ b 8 i,j=1 i j (31) C (0)eff i (µ b )C (0)eff j (µ b )f ij (δ). (32) The contribution denoted by D 2 is independent of the cutoff parameter δ. By differentiating D with respect to µ b and using the RGE (6), one can easily find out that the leading µ b -dependence cancels out in D. When calculating D 2 in our numerical analysis, we consistently set the O(α 2 s) term to zero. The contribution denoted by A vanishes in the limit α s 0. The first term in A contains the exponentiated (infrared) logarithms of δ which remain after the IR divergences cancel between the virtual and Bremsstrahlung corrections to b X s γ. The terms denoted by r i in D depend on what convention is chosen for additive constants in the functions f i7 (δ) and f 78 (δ). We fix this convention by requiring that all f ij (δ) vanish in the formal limit δ 0. 7

9 In order to complete the full presentation of b X s γ results, it remains to give explicitly the constants r i and the functions f ij (δ). The constants r 8 and r 2 are exactly as given in ref. [9] r 8 = 4 27 ( π2 6iπ), (33) r 2 = 2 { π 2 z 3/ [ π ζ(3) + ( π 2 )L + 108L L 3] z + 16πi + [ π 2 + ( π 2 )L + 36L 3] z 2 + [ 54 84π L 756L 2] z } 3 { [ π 2 + 9L + 9L 2] z + [ 3π 2 + 9L 2] z 2 + [28 12L] z 3 } + O(z 4 L 4 ), (34) where z = m 2 c /m2 b and L = ln z. In our operator basis (4), the constant r 1 does not vanish, but is proportional to r 2 r 1 = 1 6 r 2. (35) The constant r 7 turns out to be 5 r 7 = π2. (36) The quantities r 3,..., r 6 originate from two-loop matrix elements of the penguin operators P 3,..., P 6. They remain the only still unknown elements of the formally complete NLO analysis of b X s γ. However, the Wilson coefficients C i are very small for i = 3,..., 6. In consequence, setting r 3,..., r 6 to zero (which we do in our numerical analysis) has a negligible effect, i.e. it can potentially affect the decay rate by only about 1%. Let us now turn to the functions f ij (δ). One can easily obtain f 77 (δ) from eqns. (23) and (24) of ref. [16]. (These equations were also the source for our r 7 and the exponent in eqn. (32) above.) f 77 (δ) = 10 3 δ δ2 2 9 δ3 + 1 δ(δ 4) lnδ. (37) 3 The function f 88 (δ) can be found by integrating eqn. (18) of the same paper [16] f 88 (δ) = 1 27 { 2 ln m b m s [ δ 2 + 2δ + 4 ln(1 δ) ] 5 It differs from the one given in ref. [9] where a formal limit δ 1 was taken and the contribution due to f 77 (1) was absorbed into r 7. 8

10 +4Li 2 (1 δ) 2π2 3 δ(2 + δ) lnδ + 8 ln(1 δ) 2 } 3 δ3 + 3δ 2 + 7δ. (38) The logarithm of m s in the above equation is the only point in our analysis, at which the s-quark mass cannot be neglected. Collinear divergences which arise in the massless s-quark limit have been discussed in ref. [26]. The properly resumed photon spectrum was there found to be finite in the m s 0 limit, and suppressed with respect to the case when ln(m b /m s ) in f 88 (δ) is used. Even with m s as small as 0.1 GeV, the term proportional to ln(m b /m s ) is negligible, i.e. it affects the decay rate by less than 1%. Therefore, for simplicity, we will just use the naive expression (38) with m b /m s = 50 in our numerical analysis. The functions f 22 (δ), f 27 (δ), f 28 (δ) and f 78 (δ) can be found from appendix B of ref. [9]. After performing some of the phase-space integrations, one finds f 22 (δ) = 16z 27 f 27 (δ) = 8z2 9 δ [ (1 δ)/z δ 0 (1 δ)/z 0 G(t) dt (1 zt) + 1 t 2 ( dt Re G(t) + t ) /z 1/z (1 δ)/z dt (1 zt) 2 G(t) + 1 2, t 2 ( dt (1 zt)re G(t) + 2) t ], (1 δ)/z f 28 (δ) = 1 3 f 27(δ), f 78 (δ) = 8 [ Li 2 (1 δ) π2 9 6 δ lnδ δ 1 4 δ2 + 1 ] 12 δ3, (39) where, as before, z = m 2 c /m2 b and 2 arctan 2 t/(4 t), for t < 4 G(t) = π 2 /2 + 2 ln 2 [( t + t 4)/2] 2iπ ln[( t + t 4)/2], for t 4. (40) In our operator basis (4), the functions f 1j (δ) do not vanish, but are proportional to the functions f 2j (δ) f 11 (δ) = 1 36 f 22(δ), f 12 (δ) = 1 3 f 22(δ), f 17 (δ) = 1 6 f 27(δ), f 18 (δ) = 1 6 f 28(δ). The only functions f ij (δ) which have not yet been given explicitly are the ones with at least one index corresponding to the penguin operators P 3,..., P 6. As checked in ref. [17], these functions cannot affect the decay rate by more than 2%. Therefore, we neglect them in our numerical analysis. However, for completeness, a formula from which they can be obtained is given in our Appendix B. 4. The decay rate given in eqn. (30) suffers from large uncertainties due to m 5 b,pole and the CKM angles. One can get rid of them by normalizing b X s γ decay rate to the semileptonic 9

11 decay rate of the b-quark Γ[b X c e ν e ] = G2 F m5 b,pole κ(z) V 192π 3 cb 2 (41) g(z) where is the phase-space factor, and g(z) = 1 8z + 8z 3 z 4 12z 2 lnz κ(z) = 1 2α s(m b ) h(z) 3π g(z) ( ) z = m2 c,pole m 2 b,pole (42) (43) is a sizable next-to-leading QCD correction to the semileptonic decay [27]. The function h(z) has been given analytically in ref. [28] ( 25 h(z) = (1 z 2 ) z + 25 ) 4 z2 ( 17 +(1 z 2 ) z + 17 ) 3 z2 ln(1 z) 4(1 + 30z 2 + z 4 ) ln z ln(1 z) [ (1 + 16z 2 + z 4 )[6Li 2 (z) π 2 ] 32z 3/2 (1 + z) Thus, the final perturbative quantity we consider is the ratio ( + z ln z z 4 3 z ) 3 z3 + z 2 ln 2 z (36 + z 2 ) π 2 4Li 2 ( z) + 4Li 2 ( z) 2 ln z ln ( )] 1 z 1 +. z R quark (δ) = Γ[b X sγ] Eγ>(1 δ)emax γ Γ[b X c e ν e ] = V ts V tb 2 V cb 2 6α em πg(z) F ( D 2 + A ). (44) where D and A are given in eqns. (31) and (32), respectively, and F = 1 κ(z) ( mb (µ = m b ) m b,pole ) 2 = 1 ( 1 8 κ(z) 3 ) α s (m b ). (45) π 5. Let us now turn to the numerical results. Besides the five parameters listed above eqn. (25), the quantity R quark (δ) in eqn. (44) depends on a few more SM parameters. They are the following: (i) The ratio m c,pole /m b,pole. Similarly to ref. [9], we will use m c,pole m b,pole = 0.29 ± 0.02, (46) which is obtained from m b,pole = 4.8 ± 0.15 GeV and m b,pole m c,pole = 3.40 GeV. (ii) The b-quark mass itself. Apart from the ratio m c /m b, the b-quark mass enters our formulae only in two places: in the explicit logarithm in D (31) and in the factor F (45), as the argument of α s. In both cases, the m b -dependent terms are next-to-leading, and their 10

12 m b -dependence is logarithmic. Changing m b from 4.6 GeV to 5 GeV in these places affects b X s γ decay rate by less than 2%. Thus, it does not really matter whether one uses the pole mass or the MS mass there. We choose to use m b,pole = 4.8 ± 0.15 GeV. (iii) The electromagnetic coupling constant α em. It is not a priori known whether this constant should be renormalized at µ m b or µ M W, because the QED corrections have not been included. We allow α em to vary between α em (m b ) and α em (M W ), i.e. we take αem 1 = ± 2.3. (iv) The ratio of CKM factors. We will use V ts V tb = ± (47) V cb given in ref. [29]. It corresponds to gaussian error analysis, which is in accordance with interpreting the quoted error as a 1σ uncertainty. The dependence of R quark (δ) on δ is shown in Fig. 1. The middle curve is the central value. The uncertainty of the prediction is described by the shaded region. It has been found by adding in squares all the parametric uncertainties mentioned above in points (i) (iv) and in eqns. (25) (27). R quark (δ) [10 3 ] δ Figure 1: R quark (δ) as a function of δ. The ratio R quark (δ) is divergent in the limit δ 1, which is due to the function f 88 (δ). In 11

13 this limit, the physical quantity to consider would be only the sum of b X s γ and b X s decay rates, in which this divergence would cancel out. The divergence at δ 1 is very slow. It manifests itself only by a small kink in Fig. 1 where values of δ up to were included. For δ = 0.99, the contribution from f 88 (δ) to R quark (δ) is well below 1%. In order to allow easy comparison with previous experimental and theoretical publications, we choose δ max = 0.99 as the first particular value of δ to consider. This value of δ will be assumed in discussing the total B X s γ decay rate below. A value of δ which is closer to what is actually measured is δ = z (m c /m b ) 2. It corresponds to counting only the photons with energies above the charm production threshold in the b-quark rest frame. This is the second value of δ we will consider below. We find R quark (δ max ) = (3.12 ± 0.29) 10 3, (48) R quark (δ = z) = (2.86 ± 0.24) (49) In both cases, the uncertainty is below 10%. The dominant sources of it are m c /m b and µ b. The relative importance of various uncertainties is shown in Table. 1. Source α s (M Z ) m t,pole µ b m c,pole /m b,pole m b,pole α em CKM angles δ = δ max 2.5% 1.7% 6.6% 5.2% 0.5% 1.9% 2.1% δ = z 2.2% 1.7% 3.3% 6.6% 0.6% 1.9% 2.1% Table I. Uncertainties in R quark (δ) due to various sources. Let us note that setting δ = z introduces an additional spurious inaccuracy due to m c /m b in R quark. However, this additional uncertainty is significantly smaller than the original inaccuracy due to m c /m b for δ uncorrelated with m c /m b. It is quite surprising that the µ b -dependence is weaker for δ = z than for δ = δ max. Naively, one would expect larger µ b -dependence for small values of δ, when the first term in eqn. (32) becomes more and more important. Clearly, δ = z is not small enough for this effect to show up. The weaker µ b -dependence for δ = z seems to be quite accidental. 6. In the end, we need to pass from the calculated b-quark decay rates to the B-meson decay rates. Relying on the Heavy Quark Effective Theory (HQET) calculations we write B(B X s γ) = B(B X c e ν e ) R quark (δ max ) 12 ( 1 δnp sl m 2 b ) + δnp rad, (50) m 2 b

14 where δ NP sl and δ NP rad parametrize nonperturbative corrections to the semileptonic and radiative B-meson decay rates, respectively. Following refs. [3, 30], we express δ NP sl and δ NP rad in terms of the HQET parameters λ 1 and λ 2 δsl NP = 1 ( ) 3 2 λ 6(1 z)4 1 + λ 2, (51) 2 g(z) δ NP rad = 1 2 λ λ 2, (52) where g(z) is the same as in eqn. (42). The poorly known parameter λ 1 cancels out in the r.h.s. of eqn. (50) which depends only on the difference δ NP sl from B B mass splitting δ NP rad. The value of λ 2 is known λ 2 = 1 4 (m2 B m2 B ) 0.12 GeV2. (53) The two nonperturbative corrections in eqn. (50) are both around 4% in magnitude, and they tend to cancel each other. In effect, they sum up to only around 1%. Such a small number has to be taken with caution, because the four-quark operators P 1,..., P 6 have not been included in the calculation of δrad NP. Contributions from these operators could potentially give one- or two-percent effects. Nevertheless, it seems reasonable to conclude that the total nonperturbative correction to eqn. (50) is well below 10%, i.e. it is smaller than the inaccuracy of the perturbative calculation of R quark (δ max ). When the above caution is ignored and B(B X c e ν e ) = (10.4 ± 0.4)% [23] is used, one obtains the following numerical prediction for B X s γ branching ratio B(B X s γ) = (3.28 ± 0.33) (54) The central value of this prediction is outside the 1σ experimental error bar in eqn. (2). 6 However, the experimental and theoretical error bars practically touch each other. Therefore, we conclude that the present B X s γ measurement remains in agreement with the Standard Model. This conclusion holds in spite of that the theoretical uncertainty is now more than twice smaller than in the previously available leading order prediction [1]. In the measurement of B X s γ branching fraction, one needs to choose a certain lower bound on photon energies. Instead of talking about the total decay rate, it is convenient to 6 We identify the b-quark decay rate given by CLEO with the B-meson decay rate. Large statistical errors in the CLEO result make this identification acceptable. 13

15 count only the photons with energies above the charm production threshold in the B-meson rest frame E γ > m2 B m2 D 2m B 2.31 GeV. (55) Most of the photons in B X s γ survive this energy cut, and a huge background from charm production is removed. B(B X s γ) above. Let us denote the corresponding branching fraction by If the b-quark was infinitely heavy, it would not move inside the B-meson. The restriction (55) would then be equivalent to setting δ = z at the quark level. Since the Fermi motion of the b-quark inside the B-meson is a O( Λ/m b ) effect, we can write similarly to eqn. (50) B(B X s γ) above = B(B X c e ν) R quark (δ = z) (1 + O( Λ/m b )). (56) The O( Λ/m b ) part of the latter equation can be calculated using specific models for the Fermi motion of the b-quark inside the B-meson, as it has been done e.g. in ref. [16] (see also ref. [31]). Performing a similar calculation with use of the NLO values of the Wilson coefficients is beyond the scope of the present paper. Such an analysis will be necessary when more statistics allows to reduce experimental errors in measurements of B(B X s γ) above. 7. To conclude, we have presented practically complete NLO formulae for b X s γ decay in the Standard Model. They include previously published contributions as well as our new results for three-loop anomalous dimensions. Our prediction for B X s γ branching fraction in the Standard Model is (3.28 ± 0.33) 10 4 which remains in agreement with the CLEO measurement at the 2σ level. Clearly, an interesting test of the SM will be provided when more precise experimental data are available. More importantly, the B X s γ mode will have more exclusion power for extensions of the SM. 14

16 Appendix A. Here, we give the eight evanescent operators we have used in our anomalous dimension computation. Giving them explicitly is necessary in order to fully specify our renormalization scheme, and thus give meaning to the anomalous dimension matrix presented in eqn. (9). P 11 = ( s L γ µ1 γ µ2 γ µ3 T a c L )( c L γ µ 1 T a b L ) 16P 1 P 12 = ( s L γ µ1 γ µ2 γ µ3 c L )( c L γ µ 1 b L ) 16P 2 P 15 = ( s L γ µ1 γ µ2 γ µ3 γ µ4 γ µ5 b L ) q ( qγ µ 1 γ µ 4 γ µ 5 q) 20P P 3 P 16 = ( s L γ µ1 γ µ2 γ µ3 γ µ4 γ µ5 T a b L ) q ( qγ µ 1 γ µ 4 γ µ 5 T a q) 20P P 4 P 21 = ( s L γ µ1 γ µ2 γ µ3 γ µ4 γ µ5 T a c L )( c L γ µ 1 γ µ 4 γ µ 5 T a b L ) 20P P 1 (57) P 22 = ( s L γ µ1 γ µ2 γ µ3 γ µ4 γ µ5 Tc L )( c L γ µ 1 γ µ 4 γ µ 5 b L ) 20P P 2 P 25 = ( s L γ µ1 γ µ2 γ µ3 γ µ4 γ µ5 γ µ6 γ µ7 b L ) q ( qγ µ 1 γ µ 4 γ µ 5 γ µ 6 γ µ 7 q) 336P P 3 P 26 = ( s L γ µ1 γ µ2 γ µ3 γ µ4 γ µ5 γ µ6 γ µ7 T a b L ) q ( qγ µ 1 γ µ 4 γ µ 5 γ µ 6 γ µ 7 T a q) 336P P 4. Appendix B. Here, we express the functions f ij (δ) present in eqn. (32) in terms of the quantities M kl (t, u) given in eqn. (27) of ref. [17]. 7 For each (ij) (77) and i j, f ij (δ) is given by f ij (δ) = 2 3 (2 δ ij) 8 [ 1 W ik k,l=1 k l 1 δ 1 ] du dt M kl (t, u) W jl, (58) 1 u where Ŵ = (59) In four space-time dimensions, the matrix Ŵ transforms our operator basis (4) into the one used in ref. [17] P i D=4 === 8 W ik Ô k. (60) k=1 7 There is a global factor of 2 missing in the expression for M 78 given in eqn. (27) of ref. [17]. 15

17 References [1] A. J. Buras, M. Misiak, M. Münz, and S. Pokorski, Nucl. Phys. B424 (1994) 374. [2] J. Hewett, Top ten models constrained by b sγ, published in proceedings Spin structure in high energy processes, p. 463, Stanford, [3] A.F. Falk, M. Luke and M. Savage, Phys. Rev. D49 (1994) [4] M.S. Alam et al., Phys. Rev. Lett. 74 (1995) [5] A. Ali and C. Greub Phys. Lett. B293 (1992) 226. [6] B. Grinstein, R. Springer and M.B. Wise, Nucl. Phys. B339 (1990) 269. [7] R. Grigjanis, P.J. O Donnell, M. Sutherland and H. Navelet, Phys. Lett. B213 (1999) 355, Phys. Lett. B286 (1992) 413(E). [8] M. Ciuchini, E. Franco, G. Martinelli, L. Reina and L. Silvestrini, Phys. Lett. B316 (1993) 127; M. Ciuchini, E. Franco, L. Reina and L. Silvestrini, Nucl. Phys. B421 (1994) 41. [9] C. Greub, T. Hurth and D. Wyler, Phys. Lett. B380 (1996) 385, Phys. Rev. D54 (1996) [10] A. J. Buras, M. Jamin, M. E. Lautenbacher and P.H. Weisz, Nucl. Phys. B370 (1992) 69, Nucl. Phys. B375 (1992) 501 (addendum). [11] M. Ciuchini, E. Franco, G. Martinelli and L. Reina, Nucl. Phys. B415 (1994) 403. [12] M. Misiak and M. Münz, Phys. Lett. B344 (1995) 308. [13] K. Adel and Y.P. Yao, Phys. Rev. D49 (1994) [14] P. Cho, B. Grinstein, Nucl. Phys. B365 (1991) 279. [15] A. Ali and C. Greub, Z. Phys. C49 (1991) 431, Phys. Lett. B259 (1991) 182. [16] A. Ali and C. Greub, Phys. Lett. B361 (1995) 146 (equation numbers refer to the preprint version hep-ph/ ). [17] N. Pott, Phys. Rev. D54 (1996) 938. [18] F.J. Gilman and M.B. Wise Phys. Rev. D20 (1979) [19] M. J. Dugan and B. Grinstein, Phys. Lett. B256 (1991) 239. [20] S. Herrlich and U. Nierste, Nucl. Phys. B455 (1995) 39. [21] K. Chetyrkin, M. Misiak and M. Münz, in preparation. [22] T. Inami and C. S. Lim, Prog. Theor. Phys. 65 (19) 297, 65 (19) 1772E. [23] Particle Data Group, Phys. Rev. D54 (1996) 1. [24] M. Schmelling, plenary talk given at the ICHEP96 conference, Warsaw, July

18 [25] P. Tipton, plenary talk given at the ICHEP96 conference, Warsaw, July [26] A. Kapustin, Z. Ligeti and H. D. Politzer, Phys. Lett. B357 (1995) 653. [27] N. Cabibbo and L. Maiani, Phys. Lett. B79 (1978) 109. [28] Y. Nir, Phys. Lett. B221 (1989) 184. [29] A.J. Buras, plenary talk given at the ICHEP96 conference, Warsaw, July 1996 (hep-ph/ ). [30] M. Neubert, Proc. of the 17th Symposium on Lepton-Photon Interactions, Beijing, 1995, ed. by Z. Zhi-Peng and C. He-Sheng, World Scientific 1996 (and references therein). [31] M. Neubert, Phys. Rev. D49 (1994)