CP asymmetry in flavour-specific B decays beyond leading logarithms arxiv:hep-ph/ v2 21 Aug 2003

Transcription

1 hep-ph/00744 PITHA 0/05 LMU 16/0 FERMILAB-Pub-0/14-T July 00 CP asymmetry in flavour-specific B decays beyond leading logarithms arxiv:hep-ph/00744v 1 Aug 00 Martin Beneke 1, Gerhard Buchalla, Alexander Lenz and Ulrich Nierste 4 1 Institut für Theoretische Physik E, RWTH Aachen, Sommerfeldstraße 8, D-5074 Aachen, Germany. Ludwig-Maximilians-Universität München, Sektion Physik, Theresienstraße 7, D-80 München, Germany. Fakultät für Physik, Universität Regensburg, D-9040 Regensburg, Germany. 4 Fermi National Accelerator Laboratory, Batavia, IL , USA. Abstract We compute next-to-leading order QCD corrections to the CP asymmetry a fs = Im(Γ 1 /M 1 ) in flavour-specific B d,s decays such as B d Xlν l or B s D s π +. The corrections reduce the uncertainties associated with the choice of the renormalization scheme for the quark masses significantly. In the Standard Model we predict a d fs = (5.0±1.1) As a by-product we also obtain the width difference in the B d system at next-to-leading order in QCD. PACS numbers: 11.0.Er, 1.8.Bx, 1.5.Hw, Nd Keywords: asymmetry, CP; B0 anti-b0, mixing angle; B/s0, semileptonic decay 1 Preliminaries B d and B s mesons mix with their antiparticles. The time evolution of the B q B q system (with q = d or s) is characterized by two hermitian matrices, the mass matrix M q and the decay matrix Γ q. The oscillations between the flavour eigenstates B q and B q involve the three physical quantities M q 1, Γ q 1 and φ q = arg( M q 1/Γ q 1) (see e.g. 1). They are related to the mass and width differences of the B q system as M q = M q 1, Γ q = Γ q L Γ q H = Γ q 1 cosφ q, (1)

2 CP asymmetry in flavour-specific B decays beyond leading logarithms q b q u,c b u,c b u,c q b u,c q Figure 1: Leading order contribution to Γ 1 (left) and a sample NLO diagram (right). The crosses denote effective B = 1 operatorstriggering the b decay. The full set of NLO diagrams can be found in 7. where Γ q L and Γq H denote the widths of the lighter and heavier mass eigenstate, respectively. Here and in the following we neglect tiny corrections of order Γ q 1/M1 q. The CP-violating phase φ q can be measured through the CP asymmetry a q fs in flavourspecific B q f decays, which means that the decays B q f and B q f are forbidden : a q fs = Γ(B q(t) f) Γ(B q (t) f) Γ(B q (t) f)+γ(b q (t) f) = Im Γq 1 M q 1 = Γ q M q tanφ q. () Here B q (t) and B q (t) denote mesons which are tagged as a B q and B q at time t = 0, respectively. An additional requirement in Eq. () is the absence of direct CP violation in B q f, which is equivalent to f B q = f B q. For example, a s fs can be obtained through B s Ds π+. The standard way to access a q fs uses B q Xl ν l decays, which justifies the name semileptonic CP asymmetry for a q fs. The measurement of aq fs does not require tagging (seee.g.). Afurthermethodtoaccess a q fs usesthefullyinclusive, taggedb decayasymmetry discussed in 4. a q fs is small because of two suppression factors: First Γ 1/M 1 = O(m b/mw) suppresses a q fs to the percent level. Second there is a GIM suppression factor m c /m b reducing aq fs by another order of magnitude. This GIM suppression is lifted if new physics contributes to argm 1. Therefore a q fs is very sensitive to new CP phases 1,5. Up to now, the Standard Model (SM) prediction for a q fs was only known in the leading-logarithmic approximation. The unknown next-to-leading order (NLO) QCD corrections were identified as the largest theoretical uncertainty in a q fs 5. While NLO corrections were calculated long ago for Mq 1 6, only certain portions of the QCD corrections to Γ q 1 (relevant to Γ s ) were known so far 7. In Sect. we compute the missing pieces of the latter. Predictions for a q fs and Γ d can be found in Sect.. Γ q 1 at next-to-leading order in QCD In this section we specify the discussion to the case q = d and omit the index q. The generalization of our results to Γ s 1 is straightforward. Γ 1 is an inclusive quantity stemming from decays into final states common to B and B. It can be computed with the help of the heavy quark expansion (HQE) 8 from diagrams like those in Figure 1. The HQE is a simultaneous expansion in Λ QCD /m b and α s (m b ). Corrections of order Λ QCD /m b to Γ 1 have been calculated in 9,10 and applied to a fs in 5.

3 Γ q 1 at next-to-leading order in QCD We decompose Γ 1 as Γ 1 = λ c Γcc 1 + λ cλ u Γ uc 1 + λ u Γuu 1 with the CKM factors λ i = Vid V ib for i = u,c,t. The coefficients Γ ab 1, a,b = u,c in Eq. (), which are computed from diagrams like those in Figure 1, are positive. We present the new NLO expressions for the coefficient Γ uc 1 in the appendix. Γcc 1 has already been given at NLO in 7, and Γ uu 1 can be inferred by taking the limit z 0 in Γcc 1. It is convenient to write Γ 1 = λ t Γ cc 1 + (Γ uc 1 Γ cc M 1 M 1) λ u + (Γ uc 1 Γ cc 1 Γ uu 1 λ 1) λ u t λ t ( ( ) λ = 10 4 u c 1 +c B S B +c m + a 1 +a B S B +a m ) λu λ t + b 1 +b B S B +b m λ t (). (4) Here B = B(µ = m b ) and B S = B S (µ = m b ) parameterize the hadronic matrix elements of the local B= operators Q and Q S : Q = qγ µ (1 γ 5 )bqγ µ (1 γ 5 )b, B d Q(µ ) B d = 8 f B d M B d B(µ ), Q S = q(1+γ 5 )bq(1+γ 5 )b, B d Q S (µ ) B d = 5 f B d M B d B S (µ ) = 5 f B d M B d M B d m b (µ )+m d (µ ) B S(µ ). (5) The mass M Bd and decay constant f Bd cancel from Eq. (4). B and B S depend on the scale µ and the renormalization scheme used in the computation of the matrix elements in Eq. (5). When combining values for B S /B with our results for c 1,, a 1, and b 1, below, one must verify that they correspond to the same scheme. Details on the renormalization scheme used by us can be found in 7. Often the parameter B S rather than B S is chosen to parameterize Q S. As shown in Eq. (5), they differ by a factor involving MS masses. m b (m b ) is smaller than the pole mass m b by roughly 0.4GeV. For the evaluation of Eq. (4) we also need the SM prediction for M 1 : M 1 = λ t with the QCD factors η B = and b B (µ) = α s (µ) 6/ 1+ α s(µ) 4π G F 1π M B d η B B(µ )b B (µ )f B d M W S ( m t M W 5165, b B (m b ) = 1.5± Note that results from lattice gauge theory are often quoted for the scale and scheme invariant parameter B = b B (µ )B(µ ) rather than B(m b ) entering Eq. (4). We use the following input for the physical parameters (where m i m i (m i )): m b = (4.5±0.08)GeV, α s (M Z ) = 0.118±0.00, B S/B = 1.4±0., ) m c = (1.0±0.05)GeV, (6) m t = (167±5)GeV, m pow b = (4.8±0.)GeV. (7)

4 4 CP asymmetry in flavour-specific B decays beyond leading logarithms The top mass mainly enters the result through S(m t/mw) in Eq. (6), which evaluates to S(m t /M W ) =.40±0.11. In the power corrections a m, b m, c m the renormalization scheme is not fixed, because corrections of order α s /m b are unknown. The expansion parameter of the HQE is the pole mass and we use m pow b = 4.8±0.GeV (and m d = 7MeV) in a m, b m and c m. For the determination of a = a 1 +a B S B +a m (8) and the analogously defined quantities b and c we take B S /B = 1.4 ± 0., which covers the range of recent lattice computations 1. We estimate the accuracy of our calculation by computing the coefficients in two schemes for the quark masses (pole and MS), as explained in the appendix. Further we vary the renormalization scale µ 1 between one half and twice the b quark mass in the corresponding scheme. The result is shown in Figure for the coefficient a, which is most relevant to a fs : While the dependence on µ 1 is small in both LO and NLO, the scheme dependence is huge in LO and reduced by roughly a factor of 4 in NLO. We quote our coefficients for the two schemes and add the errors from Eq. (7), and the uncertainty from the µ 1 -dependence in quadrature: LO, MS LO, pole NLO, MS NLO, pole a a b b c c a m b m c m (9) In the case of a m,...,c m the difference between the LO and NLO columns stems solely from the QCD factor η B. The reduction of the scheme dependence of a 1,...,c is evident from the comparison of the last two columns with the first two ones. Our final values for a, b, and c are at NLO (LO results in parentheses): a = 1.0±.4 (14.7±6.7) b = 0.±0.1 (0.6±0.5) c = 40.1±15.8 ( 6.±15.6) (10) They have been obtained by averaging the results in the pole scheme and the MS scheme for central values of the input parameters. The error from scheme dependence was taken to

5 Γ q 1 at next-to-leading order in QCD ÄÇ ÔÓÐ 15 ÈË Ö Ö ÔÐ Ñ ÒØ 10 ÆÄÇ ÔÓÐ ÆÄÇ ÅË ÄÇ ÅË ½ Figure : Dependence of a on the scale µ 1. The solid (dashed) lines show the NLO (LO) results. be half the difference between the results in the two schemes. The errors quoted in Eq. (10) were obtained by combining in quadrature the latter error with the uncertainties in the MS scheme from scale dependence (µ 1 ), m c, m b, m t, α s (M Z ), B S /B, and the b-mass in the power corrections. In order to understand the size of the coefficients a, b, c at leading and next-to-leading order and the impact of various uncertainties, it is instructive to expand in the small parameter z = m c/m b 0.1. The leading terms in this expansion behave as follows: a 1 a b 1,b c 1,c a m b m c m LO z z z 1 z z 1 (11) NLO α s z,α s zlnz α s z α s z α s Here we have displayed the coefficients a i, b i and c i separately, indicating the leading order terms and the NLO corrections. In the SM the CP asymmetry a fs does not depend on c i, but only on a i and b i, on which we shall focus for the moment. Both a and b exhibit an interesting pattern of GIM suppression, which leads to a pronounced hierarchy among the different contributions. All of the coefficients ofa fs have tovanishasz 0. Thedominant termisa 1, while a issuppressed byone, b 1, even by two additional powers of z at LO. This strong hierarchy is alleviated at NLO, where the z and z terms receive corrections of order α s z and α s z. Hence they are still parametrically smaller than a 1, which remains the most important coefficient. As a consequence of this pattern, the coefficients b 1, get larger relative corrections at NLO, but remain strongly suppressed in comparison to a 1. This suppression is also not changed by the power corrections b m. Thus b has only a minor impact on a fs. An additional welcome feature is the suppression of a, which considerably reduces the dependence on the hadronic matrix elements B S/B.

6 6 CP asymmetry in flavour-specific B decays beyond leading logarithms We emphasize that the dominant term a 1 is free of hadronic uncertainties since the matrix element B in Γ 1 cancels against the identical quantity in M 1. It can be seen from Eq. (11) that power corrections to a are suppressed by an additional factor of z. As a result of all these properties, a fs is quite accurately known in the SM, once the NLO QCD effects are taken into account. Note that the latter are important to eliminate the sizable scheme ambiguity of the leading order calculation. We remark that the α s zlnz term in a 1 is peculiar to the choice of pole masses z = m c,pole /m b,pole, which at one-loop order is equivalent to z = m c (m c)/m b (m b). Expressing the results in terms of z = m c (m b)/m b (m b), the zlnz term is eliminated. As discussed in 11 the absence of these terms holds to all orders in α s. Finally, at NLO the overall uncertainty in a and b comes predominantly from m c and from the residual scheme dependence. The situation is different for c, which is enhanced relative to a, b. Here sizable uncertainties are still present at NLO from the dependence on B S /B, power corrections and, to a lesser extent, also from residual scale and scheme dependence. The parameter c enters the width difference Γ d and, in general, the expression for a fs in the presence of new physics. In these cases one has larger theoretical uncertainties than in the SM analysis of a fs. Phenomenology In the SM the CP asymmetry for the B d system reads a d fs = Im Γ 1 M 1 = aim λ u + bim λ u 10 4, (1) λ t λ t where a and b are given in Eq. (10). In terms of Wolfenstein parameters ρ and η the CKM quantities in Eq. (1) are λ u λ t = Im λ u λ t = sinβ R t, 1 ρ i η cosβ isinβ 1 = 1 (1) (1 ρ) + η R t Im ( λu λ t ) = sinβ R t sinβ R t where β = arg( λ t /λ c ) and R t (1 ρ) + η are one angle and one side of the usual unitarity triangle. A future measurement of a d fs will allow us to constrain ρ and η within the SM using the theoretical values for a and b. This is illustrated in Figure. Using Eq. (10) and 1 (14) R t = 0.91±0.05, β = (.4±1.4) (15) we predict for a d fs in the SM a d fs = (5.0±1.1) 10 4 (16)

7 Phenomenology ÈË Ö Ö ÔÐ Ñ ÒØ Figure : Constraints in the ( ρ, η) plane implied by given values of the CP asymmetry a d fs. The area between the solid pair of curves on the right represents the theoretical uncertainty at NLO, assuming a d fs = 10. Similarly, the curves on the left indicate the uncertainty for a d fs = both at NLO (solid) and at LO (dashed). The currently favoured solution for the unitarity triangle is also shown. This result is entirely dominated by the a-term in Eq. (1) since the small contribution from b is further suppressed by its CKM coefficient, which is small for standard CKM parameters. Our results can also be applied to the case of B s mesons, where Eq. (1) holds with obvious replacements. Here the term proportional to b is strongly CKM suppressed and can be neglected. SU() breaking in a is negligible as well and the result in Eq. (10) may be used. We then find (V us = 0.) a s fs = a V us R t sinβ 10 4 = (0.1±0.04) 10 4 (17) The width difference in the B d system is given by Γ d / M d = Re(Γ 1 /M 1 ). The real part of Γ 1 /M 1 can be found using Eqs. (4), (10), (1) and (15). It turns out that for the parameters in Eq. (15) the c-term yields the full result to within about %. In view of the large uncertainty of c, the contributions from a and b can be safely neglected. We then obtain the SM prediction Γ d M d = (4.0±1.6) 10, Γ d Γ d = (.0±1.) 10 (18) where the second expression follows with the experimental value M d /Γ d = This result for Γ d /Γ d is in agreement with 1,10. To the extent that SU() breaking in the ratio of bag factors B S /B can be neglected, the number for Γ/ M in Eq. (18) applies to the B s system as well. The effects of new physics in M 1 on a d fs have been discussed in 5. If magnitude and phase of M 1 are parameterized as M 1 = r de iθ d M SM 1 (19)

8 8 CP asymmetry in flavour-specific B decays beyond leading logarithms one obtains 5 ( ) ( ) a d fs = Re Γ1 sinθ d Γ1 cosθ d +Im M 1 SM rd M 1 SM rd Since the real part of Γ 1 /M 1 in the SM is much larger than the imaginary part, a fs is particularly sensitive to new physics. In this more general context our results can also be used. However, it has to be kept in mind that the SM analysis leading to Eq. (15) may no longer be true in the presence of new physics and the determination of CKM quantities then needs to be modified. To summarize, we have computed the CP violating observables a q fs at next-to-leading order in QCD. We include the effect of penguin operators in the weak Hamiltonian and the power corrections of relative order Λ QCD /m b. Our SM predictions are given in Eqs. (16) and (17). We emphasize that within the heavy-quark expansion the a q fs can be reliably computed in the SM as functions of CKM parameters. A crucial element is the small sensitivity to hadronic parameters, which enter only as the ratio B S /B and only with a suppression factor of z = (m c /m b ). After including the NLO corrections, the theoretical error ona q fs isreduced to about 0%. This is largely due to a reduction of the scheme ambiguity in the definition of quark masses by a factor of 4 in comparison with the LO result. The remaining uncertainty is larger for Γ q. The result at next-to-leading order in QCD is given in Eq. (18). The measurement of a q fs is possible using suitable flavour-specific decay modes of neutral B mesons. If it can be performed with sufficient accuracy, it will provide a significant test of the Standard Model. The large sensitivity of a q fs to new physics is reinforced by the improved theoretical analysis presented here. Note added The topic of this paper has also been addressed by Ciuchini et al. 14, who pointed out an error in an earlier preprint version of this paper. Our analytical results in Eq. (5) now agree with those in Eqs. (4-45) of 14. We thank the authors of 14 for clarifying communication. Acknowledgement The work of M.B. is supported in part by the Bundesministerium für Bildung und Forschung, Project 05 HT1PAB/, and by the DFG Sonderforschungsbereich/Transregio 9 Computergestützte Theoretische Teilchenphysik. (0) A NLO coefficients Here we collect more detailed results for the coefficients in Eq. (). The HQE expresses Γ ab 1 for the B d system as (F Γ ab 1 = G F m b 4π f B d M ab Bd (z)+p ab (z) ) 8 B ( FS ab (z)+ps ab (z) ) 5 B S +Γ ab 1,1/m b.(1)

9 A NLO coefficients 9 The short-distance coefficients F ab (z) contain the contributions from the B = 1 currentcurrent operators Q 1 and Q. The NLO results for F cc (z) and FS cc (z) have been derived in 7, where these coefficients are called F(z) and F S (z), respectively. Further F uu = F cc (0) and FS uu = FS cc (0). The coefficients P(z) and P S (z) contain the contributions from penguin operators. They come with small coefficients, which simplifies the NLO calculation 7. Our new calculation concerns F uc, FS uc, P uc and PS uc. We decompose F uc and FS uc as in 7,11: F uc (z) = C 1F uc 11(z)+C 1 C F uc 1(z)+C F uc (z), F uc ij (z) = F uc,(0) ij (z)+ α s(µ 1 ) 4π Fuc,(1) ij (z,x µ1,x µ )+O(α s) () with x µ = µ/m b and an analogous notation for FS,ij uc. The B =n operators, n = 1,, are defined at the scale µ n = O(m b ). The dependence of Fij uc on µ 1 diminishes order-by-order in α s. Throughout this paper we use the same operator definitions and renormalization schemes as in 7, with one important addition: In a fs the renormalization scheme of the quark masses is an important issue and we choose two different schemes for the computation of the a i, b i, c i in Eq. (4). For both schemes we take the MS masses m c (m c ) and m b (m b ) as the basic input. In the first scheme (pole scheme) we express the observables in terms of m b = m b,pole = m b (1+4α s (m b )/π), using the one-loop relation between pole- and MS-quark mass. In this schemewedefinethevariablez asz = (m c (m c )/m b (m b )), whichtoone-looporderisequivalent to the ratio of pole masses squared. In the second scheme (MS scheme) we take m b = m b (m b ) and replace z by z = (m c (m b )/m b (m b )), where both running masses are defined at the scale m b. The results below for the functions F uc,(1) ij (z) are valid in the pole scheme. The corresponding functions F ab,(1) ij (z) in the MS scheme are obtained via the relation F ab,(1) ij The coefficients read: (z) = F ab,(1) ij (z)+ Fab,(0) ij F uc,(0) 11 (z) = (1 z) (1+ z ) F uc,(0) 1 (z) = (1 z) (1+ z ) F uc,(0) (z) = 1 (1 z) F uc,(0) S,11 (z) = (1 z) (1+z) F uc,(0) S,1 (z) = (1 z) (1+z) (z) 8zlnz Fab,(0) ij (z). () z F uc,(0) S, (z) = (1 z) (1+z) (4) F uc,(1) 11 (z,x µ1,x µ ) = 16(1 z) (+z) Li (z)+ ln(1 z) ln(z) +

10 10 CP asymmetry in flavour-specific B decays beyond leading logarithms 4(1 z) (5+7z) ln(1 z) + ( z 10+14z 15z ) ln(z) + (1 z) (5+z) ln(x µ ) + (1 z) ( z 104z ) 6 F uc,(1) 1 (z,x µ1,x µ ) = (1 z) (+z) Li (z)+ ln(1 z) ln(z) ( (1 z) ( +z +94z )) ( ( z 80+69z 16z )) ln(1 z) + 6z 6 (1 z) (17+4z) (1 z) ( z +z ) 18 4(1 z) (5+z) ln(x µ1 ) + ln(x µ ) + + ln(z) + F uc,(1) (5 8z) (1 z) (1+z) ln(1 z) ln(z) (z,x µ1,x µ ) = Li (z)+ + (1 z) ( 7+z +z ) ( ( z 6 9z 0z +z )) ln(1 z) + 6z 6 (1 z) (5+4z) (1 z) ( 1+4z) (1 z) (4 z) ln(x µ1 ) + ln(x µ ) + π + (1 z) ( 16 95z +44z ) 18 ln(z) + F uc,(1) S,11 (z,x µ 1,x µ ) = (1 z) (1+z) Li (z)+ 8(1 z) ( 4+14z z ) ln(1 z) + (1 z) (1+z) ln(1 z) ln(z) + ( 8z +z 1z +z ) ln(z) + ln(x µ ) + 4 (1 z) ( 10 z +1z ) F uc,(1) S,1 (z,x µ 1,x µ ) = 64(1 z) (1+z) Li (z)+ 4(1 z) ( 1+15z +47z 1z ) ln(1 z) + z ln(1 z) ln(z) +

11 REFERENCES 11 ( 4z 8+9z 87z +1z ) ln(z) + 16(1 z) (1+z) ln(x µ1 ) + (1 z) ( 10 7z +107z ) 9 64(1 z) (1+z) ln(x µ ) + F uc,(1) 16 (1 4z) (1 z) (1+z) S, (z,x ln(1 z) ln(z) µ 1,x µ ) = Li (z)+ + 4(1 z) (1+z) ( 1+1z +z ) ( 4z z +18z z ) ln(1 z) + ln(z) + z 16(1 z) (1+z) 8 (1 z) (1+z) (1 z) (1+z) ln(x µ1 ) + ln(x µ ) + π + 8 (1 z) ( 5 8z +19z ) 9 In terms of the function P(z) used in 7 the penguin coefficients in Eq. (1) read P cc (z) = P(z), P uu = P(0) and with P uc (z) = P(z)+P(0) P uc = α s(µ 1 ) 4π + P uc, P uc S (z) = P S(z)+P S (0) C (µ 1) 1 (1+z) 1 4z 18 (5) 8 P uc (6) lnz (1+z) 1 4zlnσ 4z and σ = (1 1 4z)/(1+ 1 4z). P uc is of order z and numerically negligible. The power corrections Γ ab 1,1/m b were first obtained for ab = cc,uu in 9 and for ab = uc in 10. We have re-computed the case ab = uc here, confirming the results of 10. In the notation of 9 we find (... B... B ) Γ uc 1,1/m b References = G F m b 4πM B (1 z) (1+z)K R 0 (1+z)(K 1 R 1 +K R 1 ) 1+z +z (K 1 R +K R ) 1z 1 z 1 z (K 1 R +K R ) (7). (8) 1 K. Anikeev et al., B physics at the Tevatron: Run II and beyond, hep-ph/001071, Chapters 1. and 8..

12 1 CP asymmetry in flavour-specific B decays beyond leading logarithms J. S. Hagelin and M. B. Wise, Nucl. Phys. B 189 (1981) 87; J. S. Hagelin, Nucl. Phys. B 19 (1981) 1; A. J. Buras, W. Slominski and H. Steger, Nucl. Phys. B 45 (1984) 69. I. Dunietz, R. Fleischer and U. Nierste, Phys. Rev. D 6 (001) M. Beneke, G. Buchalla and I. Dunietz, Phys. Lett. B 9 (1997) 1. 5 S. Laplace, Z. Ligeti, Y. Nir and G. Perez, Phys. Rev. D 65 (00) A. J. Buras, M. Jamin and P. H. Weisz, Nucl. Phys. B 47 (1990) M. Beneke, G. Buchalla, C. Greub, A. Lenz and U. Nierste, Phys. Lett. B 459 (1999) M.A. ShifmanandM. B.Voloshin, in: Heavy Quarks ed. V. A.KhozeandM.A. Shifman, Sov. Phys. Usp. 6 (198) 87; M. A. Shifman and M. B. Voloshin, Sov. J. Nucl. Phys. 41 (1985) 10 Yad. Fiz. 41 (1985) 187; M. A. Shifman and M. B. Voloshin, Sov. Phys. JETP 64 (1986) 698 Zh. Eksp. Teor. Fiz. 91 (1986) 1180; I. I. Bigi, N. G. Uraltsev and A. I. Vainshtein, Phys. Lett. B 9 (199) 40 Erratum-ibid. B 97 (199) M. Beneke, G. Buchalla and I. Dunietz, Phys. Rev. D 54 (1996) A. S. Dighe, T. Hurth, C. S. Kim and T. Yoshikawa, Nucl. Phys. B 64 (00) M. Beneke, G. Buchalla, C. Greub, A. Lenz and U. Nierste, Nucl. Phys. B 69 (00) S. Aoki et al. JLQCD Collaboration, Phys. Rev. D 67 (00) , hep-lat/00808; D. Becirevic, V. Gimenez, G. Martinelli, M. Papinutto and J. Reyes, Nucl. Phys. Proc. Suppl. 106 (00) 85, hep-lat/ M. Battaglia et al., The CKM matrix and the unitarity triangle, hep-ph/ M. Ciuchini, E. Franco, V. Lubicz, F. Mescia and C. Tarantino, hep-ph/00809.