arxiv: v2 [hep-ph] 9 Jun 2009

Transcription

1 Semileptonic to P-Wave Charmonia (X c X c h c Transitions within QCD Sum Rules K. Azizi H. Sundu M. Bayar Department of Physics Middle East Technical University 653 Ankara Turkey Physics Department Kocaeli University 438 Izmit Turkey e4634@metu.edu.tr hayriye.sundu@kocaeli.edu.tr melahat.bayar@kocaeli.edu.tr arxiv:9.467v hep-ph 9 Jun 9 The form factors of the semileptonic S(AVlν (l = τµe transitions where S and AV denote the scalar X c and axial vector (X c h c mesons are calculated within the framework of the three-point QCD sum rules. The heavy quark effective theory limit of the form factors are also obtained and compared with the values of the original transition form factors. The results of form factors are used to estimate the total decay widths and branching ratios of these transitions. A comparison of our results on branching ratios with the predictions of other approaches is also presented. PACS numbers:.55.hx 3..He I. INTRODUCTION The meson is the only known meson composed of two heavy quarks of different flavor and charge; a charm quark and a bottom antiquark. It were discovered by the collider detector at Fermilab (CDF Collaboration in pp collision via the decay mode J/ψl ± ν at s = 8 TeV. Discovery of the meson has demonstrated the possibility of the experimental study of the charm-beauty system and has created considerable interest in its spectroscopy When the large hadron collider (LHC runs a plenty number of events which are expected to be about 8 per year with the the luminosity values of L = 34 cm s and s = 4 TeV will be produced 7 8. Therefore not only experimental but also theoretical study on mesons will be of great interests in many respects. Among the B mesons the carries a distinctive signature and has reached great interest recently for the following reasons: Firstly the meson decay channels are expected to be very rich in comparison with other B mesons so investigation of such type of decays can be used in the calculation of the Cabibbo-Kabayashi-Maskawa(CKM matrix elements leptonic decay constant as well as the origin of the CP and T violation. Secondly the meson because of containing the heavy quarks provides more accuracy and confidence in the understanding of the QCD dynamics. The meson can decay via the b udsc and also the c uds transitions. Among those transitions at quark level the tree level b c transition governs the to P-wave charmonia plays a significant role because this is the most dominant transition. In the literature there are several studies on the mesons in different models. Some possible meson decays such as lνγ ρ + γ K + γ and Bu l+ l Bu γ B c D lν P(DD s l + l /ν ν Dsd l+ l Xν ν and Dsγ have been studied in the frame of light-cone QCD and three-point QCD sum rules The weak productions of new charmonium in semileptonic decays of were also studied in the framework of light cone QCD sum rules in 8. In 9 a larger set of exclusive nonleptonic and semileptonic decays of the meson were investigated in the relativistic constituent quark model. Weak decays of the meson to charmonium and D mesons in the relativistic quark model have been discussed in. Moreover the (D Ds νν transitions were also studied within the relativistic constituent quark model in. Present work is devoted to the study of the S(AVlν. The long distance dynamics of such transitions can be parameterized in terms of some form factors which play fundamental role in analyzing such transitions. For evaluation of the form factors the QCD sum rules as nonperturbative approach based on the fundamental QCD Lagrangian is used. The obtained results for the form factors are used to estimate the total decay rate and branching fractions for the related transitions. The heavy quark effective theory (HQET limit of the form factors are also calculated and compared with their values. In these transitions the main contribution comes from the perturbative part since the heavy quark condensates are suppressed by inverse of the heavy quark masses and can be safely omitted and two-gluon condensate contributions are very small and we will ignore them. Note that the to P-wave charmonia transitions have also been investigated in the framework of covariant light-front quark model (CLQM the renormalization group method (RGM relativistic constituent quark model (RCQM and nonrelativistic constituent quark model (NRCQM in For more about those transitions see also

2 The outline of the paper is as follows: In Section II the some rules for the transition form factors relevant to the S(AVlν decays are obtained. Section III encompasses the calculation of the HQET limit of the form factors and section IV is devoted to the numerical analysis of the form factors and their HQET limits decay rates branching ratios conclusion and comparison of our results with the other approaches. II. SUM RULES FOR THE S(AVlν TRANSITION FORM FACTORS The X c (X c h c lν decays proceed via the b c transition at the quark level. The effective Hamiltonian responsible for these transitions can be written as: H eff = G F V cb ν γ µ ( γ 5 l c γ µ ( γ 5 b. ( We need to sandwich Eq. ( between initial and final meson states in order to obtain the matrix elements of S(AVlν. Hence the amplitude of this decay is written as follows: M = G F V cb ν γ µ ( γ 5 l < S(AV(p c γ µ ( γ 5 b (p >. ( It is necessary to calculate the matrix elements < S(AV(p cγ µ ( γ 5 b (p > appearing in Eq. (. In the S case in final state the only axial-vector part of the transition current c γ µ ( γ 5 b contribute to the matrix element stated above. However in the AV case both vector and axial-vector parts have contributions. Considering the parity and Lorentz invariances the aforementioned matrix element can be parameterized in terms of the form factors in the following way: < S(p cγ µ γ 5 b (p >= f (q P µ +f (q q µ (3 < AV(p f V (q ε cγ µ γ 5 b (p >= i (m Bc +m Xc ε µναβε ν p α p β (4 < AV(p ε cγ µ b (p > = i f (q (m Bc +m Xc ε µ f + (q (m Bc +m Xc (ε pp µ f (q (m Bc +m Xc (ε pq µ (5 where f (q f (q f V (q f (q f + (q and f (q are transition form factors and P µ = (p+p µ q µ = (p p µ. From the general philosophy of the QCD sum rules we see a hadron from two different windows. First we see it from the outside so we have a hadron with hadronic parameters such as its mass and leptonic decay constant. Second we see the internal structure of the hadron namely quarks and gluons and their interactions in QCD vacuum. In technique language we start with the main object in QCD sum rules so called the correlation function. The correlation function is calculated in two different ways: From one side it is saturated by a tower of hadrons called the phenomenological or physical side. On the other hand the QCD or theoretical side it is calculated in terms of quark andgluonsinteractingin QCDvacuum bythe help ofthe operatorproductexpansion(opewherethe shortand long distance effects are separated. The farmer is calculated using the perturbation theory (perturbative contribution however the latter is parameterized in terms of vacuum condensates with different mass dimensions. In the present work there is no light quarks which the non-perturbative contributions mainly come from their vacuum condensates and the heavy quark condensate contributions are suppressed by inverse of the heavy quark mass and can be safely removed. The two-gluon contributions are also very small and here we will ignore those contributions. Hence the only contribution comes from the perturbative part. Equating two representations of the correlation function and applying double Borel transformation with respect to the momentum of the initial and final states to suppress the contribution of the higher states and continuum sum rules for the physical quantities form factors are obtained. To proceed we consider the following correlation functions: Π µ (p p = i d 4 xd 4 ye ipx e ipy < TJ S (yjµ V ;A (J Bc (x > (6 Π µν (p p = i d 4 xd 4 ye ipx e ipy < TJ νav (yj V ;A µ (J Bc (x > (7

3 where J S (y = cuc J νav (y = cγ ν γ 5 c J Bc (x = bγ 5 c are the interpolating currents of the S AV and mesons respectively and J V µ ( = cγ µb J A µ = cγ µγ 5 b are the the vector and axial-vector parts of the transition current. In order to calculate the phenomenological or physical part of the correlator given in Eq. (6 two complete sets of intermediate states with the same quantum numbers as the interpolating currents J S(AV and J Bc are inserted As a result the following representations of the above-mentioned correlators are obtained: 3 Π µ (p p = < J S( S(p >< S(p J V ;A µ ( (p >< (p J Bc ( > (p m S (p m + (8 Π µν (p p = < Jν AV ( AV(p ε >< AV(p ε J V;A µ ( (p >< (p J Bc ( > (p m AV (p m + (9 where represents the contributions coming from higher states and continuum. The vacuum to the hadronic state matrix elements in Eq. (8 can be parameterized in terms of the leptonic decay constants as: < J S S(p > = if S < (p J Bc >= i f m m b +m c < J ν AV AV(p ε >= f AV m AV ε ν. ( Using Eqs. (3- the final expressions of the phenomenological side of the correlation functions are obtained as: Π µ (p p f S f Bc m B = c f (p m S (p m (q P µ +f (q q µ + excited states m b +m c Π µν (p p = f m f AV m AV (m b +m c (p m AV (p m f (q g µν (m Bc +m AV f +(q P µ p ν (m Bc +m AV f (q q µ p ν (m Bc +m AV + ε αβµν p α p β f V (q (m Bc +m AV + excited states ( where we will choose the structures P µ q µ ε µναβ p α p β g µν and (p µp ν ±p µp ν to evaluate the form factors f f f V f and f ± respectively. On the QCD side the aforementioned correlation functions can be calculated by the help of the OPE in the deep space-like region where p (m b + m c and p (m c. As we mentioned before the main contributions to the theoretical part of the correlation functions come from bare-loop (perturbative diagrams. To calculate those contributions the correlation functions are written in terms of the selected structures as follows: Π µ = Π per P µ +Π per q µ Π µν = Π per V ε µναβp α p β +Π per g µν + Πper + (p µp ν +p µ p ν+ Πper (p µp ν p µ p ν ( where each Π per i function is written in terms of the double dispersion representation in the following way: Π per i = (π ds ds ρ i (ss q (s p (s p + subtraction terms (3 where the functions ρ i (ss q are called the spectral densities. Using the usual Feynman integral for the bare loop diagram the spectral densities can be calculated with the help of Cutkosky rules i.e. by replacing the quark propagators with Dirac delta functions: p m πδ(p m which implies that all quarks are real. After some straightforward calculations the spectral densities are obtained as follows: ρ (ss q = N c I (ss q ρ (ss q = N c I (ss q (m b 3m c m c +A{(m b m c m c s+b{(m b m c m c s (m b +m c m c +A{(m b m c m c +s B{(m b m c m c +s

4 ρ V (ss q = 4N c I (ss q ρ (ss q = N c I (ss q ρ + (ss q = N c I (ss q ρ (ss q = N c I (ss q where I (ss q = (m c m b A+m c B +m c 4(m c C(m c m b +m c u+{m c (4s+u m b ua+ m b s +m c (s +ub m c +(m b 3m c A m c B +(m b m c D +(m b m c E m c (m c +m b A+m c B +(m b m c D +(m c m b E 4λ / (ss q λ(abc = a +b +c ac bc ab = m b m c s = s u = s+s q A = λ(ss q ( u s B = λ(ss q ( u s C = λ(ss q s+ s u+m c( 4ss +u D = λ(ss q 6 s u+ (ss +u +s (3 s +m c( 4ss +u E = λ(ss q 3 s u+ (ss +u u(3 s+m c ( 4ss +u and N c = 3 is the number of colors. The integration region for the perturbative contribution in the Eq. (3 is determined requiring that the arguments of the three δ functions vanish simultaneously. Therefore the physical region in the s and s plane is described by the following non-equality: f(ss = ss +(s+s q (m b s m c λ / (m b sm c λ/ (ss q 4 (4 (5 +. (6 Equating the coefficient of the selected structures from the phenomenological and the OPE expressions and applying double Borel transformations with respect to the variables p and p (p M p M in order to suppress the contributions of the higher states and continuum the QCD sum rules for the form factors f (q and f (q for the X c lν decay can be acquired: f (q = (m { b +m c f Bc m e m Bc /M e mx c /M f Xc (π s (m b +m c ds s (m c ds ρ (ss q θ f (ss e s/m The form factors f V f f + and f for AVlν decays are also obtained as: f i (q = κ (m b +m c η f Bc m e m Bc /M +m AV /M f AV m AV s s (π ds dsρ i (ss q θ f (ss e s/m (m b +m c 4m c e s /M. (7 s /M (8

5 where i = V± and η = m Bc +m AV for i = V± and η = m Bc +m AV for i = are considered. Here κ = + for i = ± and κ = for i = and V. In the above equations the s and s are continuum thresholds in s and s channels respectively. In order to subtract the contributions of the higher states and continuum the quark-hadron duality assumption is used i.e. it is assumed that ρ higherstates (ss = ρ OPE (ss θ(s s θ(s s. (9 Note that the double Borel transformation used in calculations is written as: ˆB (p m m (p m n ( m+n Γ(m Γ(n e m /M e m /M (M m (M 5 n. ( Now we would like to explain our reason for ignoring the contributions of the gluon condensates to the QCD side of the correlation function. These contributions for the related form factors are obtained as the following orders: f G α s f G V+ α s f G α s π G mn π G mn π G mn3 b mm c M k M l b mm c M k M l b mm3 c M k3 M l3 n +m = k +l n +m = k +l n 3 +m 3 = k 3 +l 3 + ( where α s is the strong coupling constant and M and M are Borel mass parameters. Recalling the magnitude of the αs π G =. GeV 4 36 and considering the working region of the Borel parameters (see numerical analysis section the gluon condensate contributions become very small and here we ignore those small contributions ( maximum contribution is obtained for f G which is not more than few percent. At the end of this section we would like to present the differential decay rates of the S(AVlν in terms of the transition form factors. The differential decay width for Slν is obtained as follows : ( dγ dq = q 9π 3 m 3 G F V cb λ / (m m m { S q l q (q +m l f (q (m +m S q + (m m SRef (q f (q + f (q q + (q +m l q f (q (m m S + (m m Sq Ref (q f (q + f (q q 4 ( and also the differential decay width corresponding to AVlν decays are acquired as: ( { dγ dq = 6π 4 m p G F V cb 4 (A +A q f V (4m p + f { + (A +A q f V (4m p +m p (m m AV q m AV + f f + m p (m +m AV q f m p m AV m p m (Re(f f + +f f +(m m AV f +f AV m B A c p m f +(m m AV f + +q 4 f AV m AV + (m m AVRe(f f + +q f f +q (m m AVRe(f + f q (3

6 6 where p = λ/ (m m AV q m Bc A = q (q m l I A = 6q 4(q m l(q +m li I = π ( m l q. (4 III. HEAVY QUARK EFFECTIVE THEORY LIMIT OF THE FORM FACTORS In this section we calculate the heavy quark effective theory (HQET limits of the transition form factors for S(AVlν. For this aim following references we use the parameterization y = νν = m +m S(AV q m Bc m S(AV (5 where ν and ν are the four-velocities of the initial and final meson states respectively. Next we try to find the y dependent expressions of the form factors by taking m b m c = m b z where z is given by z = y + y. In this limit the new Borel parameters T = M /(m b +m c and T = M /4m c are defined. The new continuum thresholds ν and ν are also parameterized as: and the new integration variables take the following form: The leptonic decay constants are rescaled: ν = s (m b +m c m b +m c ν = s 4m c m c (6 ν = s (m b +m c m b +m c ν = s 4m c m c. (7 ˆf Bc = m b +m c f Bc ˆfS(AV = m c f S(AV. (8 To evaluate the form factors in HQET we also need to redefine the form factors in the following form: f f = (m Bc +m S f V+ = f V+. (9 m Bc +m AV After standard calculations we obtain the y-dependent expressions of the form factors for Slν transition as follows: f = 3( +y 3+z+y( 3 z +z 8 3 e ( Λ T + Λ T π ˆf S ˆfBc z 3/4 (+y + ( +y (+ z z z θ lim mb f (vv { dν dν e ν T ν T (3

7 3( +y +y(+ z +4 z +z f = 8 π ˆf S ˆfBc z 3/4 (+y + z θ lim mb f (vv ( +y (+ z z 3 e ( Λ T + Λ T { dν dν e ν T ν T 7 (3 and for AVlν decay the y-dependent expressions of the form factors are acquired as: 3(3+ z +y + z +y z { f V = 8 ( π ˆf AV ˆfBc z 5/4 (+y(+ e Λ T + Λ T z + z θ lim mb f (vv ( +y (+ z z dν dν e ν T ν T (3 f = 3( +y +y(+ z+3 z 8 π ˆf AV ˆfBc z 5/4 (+y(3+ z + θ lim mb f (vv z ( +y (+ z z ( e Λ T + Λ T { dν dν e ν T ν T (33 3( +y (3+ z +y (+ z +5y( +z z f + = 3 π ˆf AV ˆfBc z 3/4 (+y + ( +y (+ z z z θ lim mb f (vv 3 e ( Λ T + Λ T { dν dν e ν T ν T 3( +y (3+ z +y (+ z +y(3+8 z +5z+ z f = 3 3 π ˆf AV ˆfBc z 3/4 (+y + ( +y (+ z z z θ lim mb f (vv where Λ = m Bc (m b +m c and Λ = m S(AV m c. e ( Λ T + Λ T { (34 dν dν e ν T ν T (35

8 8 IV. NUMERICAL ANALYSIS This section is devoted to the numerical analysis of the form factors their HQET limit and branching ratios. The sum rules expressions for the form factors depict that they mainly depend on the leptonic decay constants continuum thresholds s and s and Borel parameters M and M. In calculations the quark masses are taken to be m c (µ = m c =.75±.5 GeV m b = (4.7±. GeV 4 and the meson masses are chosen as: m Bc = 6.86 GeV m hc = GeV m Xc = GeV m Xc = GeV 4. For the values of the leptonic decay constants we use f Bc = (4±4 MeV and f Xc = f Xc = f hc = (34 +9 MeV3. The two-point QCD sum rules are used to determine the continuum thresholds s and s. These thresholds are not completely arbitrary and they are related to the energy of the exited states. The result of the physical quantities form factors should be stable with respect to the small variation of these parameters. Generally the s are obtained to be (m hadron Here we use s = (45±5 GeV and s = (6± GeV. Since the Borel parameters M and M are not physical quantities the form factors should not depend on them. The reliable regions for the Borel parameters M and M can be determined by requiring that not only the contributions of the higher states and continuum are effectively suppressed but the contribution of the operator with the highest dimension be small. As a result of the above-mentioned requirements the working regions are determined to be 5 GeV M 35 GeV and GeV M GeV. The numerical values of the form factors at q = for X c lν and AVlν transitions are given in the Tables I and II respectively. f ( f ( X clν.673± ±.437 TABLE I: The values of the form factors for the X clν decay at M = 5 GeV M = 5 GeV and q =. f ( f V( f +( f ( X clν.84±.5.949±.6.±.6.586±.79 h clν.84±.5.954±.8.±.6.588±.8 TABLE II: The values of the form factors for the AVlν decays at M = 5 GeV M = 5 GeV and q =. a b m fit f ( X clν f ( X clν TABLE III: Parameters appearing in the form factors of the X clν decay at M = 5 GeV and M = 5 GeV. In order to estimate the decay width of the S(AVlν transitions we need to know the q dependent form factors in the whole physical region m l q (m Bc m S(AV. Our form factors are truncated at about q = 4 GeV. To extend our results to the full physical region we search for parameterization of the form factors in such a way that in the region q 4 GeV this parameterization coincides with the sum rules predictions. The following fit parameterization is chosen for the form factors with respect to q : f i (q = a ( q + b (36 ( q m m fit fit where the values of the parameters a b and m fit for the X c lν and (X c h c lν are given in the Tables III and IV respectively.

9 9 a b m fit f ( X clν f V( X clν f +( X clν f ( X clν f ( h clν f V( h clν f +( h clν f ( h clν TABLE IV: Parameters appearing in the form factors of the X clν and h clν decays at M = 5 GeV and M = 5 GeV. q (GeV y f f (HQET f f (HQET TABLE V: Values of the form factors and their HQET limits for the X clν at M = 5 GeV M = 5 GeV T =.9 GeV and T =.94 GeV. To calculate the numerical values of the form factors at HQET limit the values of Λ =.3GeV and Λ =.86GeV(.96GeV are used for Slν ( AVlν transitions respectively (see In Tables V VI and VII we compare the values of the form factors and their HQET limits for considered transitions in the interval q 7 and corresponding values of the y. Comparing the form factors and their HQET values in those Tables we see that all form factors and their HQET limits have the same behavior with respect to the q i.e. they both growth or fail by increasing the values of q. The HQET limit of the form factors are comparable with their original values and in large q those form factors and their HQET values become very close to each other. The results presented at Tables VI and VII also indicate that the form factors and their HQET limits for X c lν and h c lν have values very close to each other since the X c and h c mesons are both axial vectors i.e. J P = + and have nearly the same mass. At the end of this section we would like to calculate the values of the branching ratios for these decays. Taking into account the q dependency of the form factors and performing integration over q from the differential decay rates in Eqs. ( 3 in the interval m l q (m Bc m S(AV and also using the total life-time of the meson τ Bc =.46±.7 s 4 we obtain the branching ratios of the related transitions as presented in Table VIII. This Table is also contain the predictions of the other approaches such as covariant light-front quark model (CLQM renormalization group method (RGM relativistic constituent quark model (RCQM and nonrelativistic constituent quark model (NRCQM These results can be tested in the future experiments. In conclusion using the QCD sum rules approach we investigated the semileptonic S(AVlν decays. The q dependencies of the transition form factors were calculated. The HQET limits of the form factors were also evaluated and compared with original form factors. The obtained results were used to estimate the total decay widths and branching ratios of these transitions. A comparison of the results for branching fractions was also presented. V. ACKNOWLEDGMENT The authors would like to thank T. M. Aliev and A. Ozpineci for their useful discussions. K. A. would like to thank TUBITAK Turkish Scientific and Research Council for their partial financial support.

10 q (GeV y f f (HQET f V f V(HQET f f (HQET f f +(HQET TABLE VI: Values of the form factors and their HQET limits for the X clν at M = 5 GeV M = 5 GeV T =.9 GeV and T =.94 GeV. q (GeV y f f (HQET f V f V(HQET f f (HQET f f +(HQET TABLE VII: Values of the form factors and their HQET limits for the h clν at M = 5 GeV M = 5 GeV T =.9 GeV and T =.94 GeV. X clν X clν h clν Present work.8±.5.46±.4.4±.4 CLQM RGM RCQM RCQM NRCQM X cτν X cτν h cτν Present work.49±.6.47±.44.37±.38 CLQM RGM RCQM RCQM NRCQM TABLE VIII: Branching ratios of the semileptonic (X cx ch clν (l = eµτ transitions in different approaches.

11 F. Abe et al. CDF Collabration Phys. Rev. D 58 4 (998. A. Adb El-Hady M. A. K. Lodhi and J. P. Vary Phys. Rev. D ( L P Fulcher Phys. Rev. D ( D. Ebert R. N. Faustov and V. O. Galkin Phys. Rev. D (3. 5 S. Godfrey Phys. Rev. D (4. 6 A. K. Rai P. C. Vinodkumar Pramana (6 hep-ph/6694v. 7 S. Stone to appear in proceedings of Heavy Flavor Physics: A Probe of Nature s Grand Design Varenna Italy July 997 hep-ph/ D.S. Du Z. Wang Phys. Rev. D (989; C.H. Chang Y.Q. Chen ibid. 48 ( ; K. Cheung Phys. Rev. Lett (993; E. Braaten K. Cheung T. Yuan Phys. Rev. D 48 R549 (993 9 T. M. Aliev M. Savci Phys. Lett. B 434 ( T. M. Aliev M. Savci J. Phys. G 4 ( T. M. Aliev M. Savci Eur. Phys. J. C 47 (6 43. T. M. Aliev M. Savci Phys. Lett. B 48 ( N. Ghahramany R. Khosravi K. Azizi Phys. Rev. D (8. 4 K. Azizi R. Khosravi Phys. Rev. D (8. 5 K. Azizi F. Falahati V. Bashiry S. M. Zebarjad Phys. Rev. D (8. 6 K. Azizi R. Khosravi V. Bashiry Eur. Phys. J. C (8. 7 K. Azizi V. Bashiry Phys. Rev. D (7. 8 Yu-Ming Wang Cai-Dian Lu Phys. Rev. D (8. 9 M. A. Ivanov J. G. Korner and P. Santorelli Phys. Rev. D (6. D. Ebert R.N. Faustov V.O. Galkin Phys. Rev. D (3. D. Ebert R.N. Faustov V.O. Galkin Eur. Phys. J. C 3 9 (3. A. Faessler Th. Gutsche M. A. Ivanov J. G. Korner V. E. Lyubovitskij Eur. Phys. J.direct C 4 8(. 3 Xiao-Xia Wang Wei Wang Cai-Dian L hep-ph/9.93v. 4 C. H. Chang Y. Q. Chen G. L. Wang and H. S. Zong Phys. Rev. D (. 5 M. A. Ivanov J. G. Korner and P. Santorelli Phys. Rev. D (6. 6 M. A. Ivanov J. G. Korner and P. Santorelli Phys. Rev. D (5 Erratum-ibid. D (7. 7 E. Hernandez J. Nieves and J. M. Verde-Velasco Phys. Rev. D (6 arxiv:hep-ph/ T. Becher H. Boos E. Lunghi JHEP7:6 (7. 9 U. Aglietti L. D. Giustino G. Ferrera A. Renzaglia G. Ricciardi L. Trentadue Phys. Lett. B (7. 3 P. Urquijo Belle Collaboration Proceedings of the XXXIII International Conference of the High Energy Physics (ICHEP 6 Moscow (6 arxiv: hep-ex/ K. Abe Belle Collaboration Proceedings of the XXXIII International Conference of the High Energy Physics (ICHEP 6 Moscow (6 arxiv: hep-ex/ C. Schwanda Belle Collaboration Phys. Rev. D (7. 33 Chao-Hsi Chang Yu-Qi Chen Guo-Li Wang Hong-Shi Zong Phys. Rev. D (. 34 P. Colangelo F. De Fazio T.N. Pham Phys. Rev. D (4. 35 Chao-Hsi Cahng XXXVIIth Rencontres de Moriond on QCD and High Energy Hadronic Interactions arxiv:hep-ph/5. 36 M. A. Shifman A. I. Vainshtein V. I. Zakharov Nucl. Phys. B 47 ( Ming Qiu Huang Phys. Rev. D 69 ( M. Neubert Phys. Rep. 45 ( T. M. Aliev K. Azizi A. Ozpineci Eur. Phys. J. C 5 ( K. Azizi M. Bayar Phys. Rev. D (8. 4 B. L. Ioffe Prog. Part. Nucl. Phys. 56 ( C. Amsler et al. Particle Data Group Phys. Lett. B 667 (8. 43 T. Huang C. W. Luo Phys. Rev. D 5 ( Y. B. Dai C. S. Huang C. Liu and S. L. Zhu Phys. Rev. D 68 4 (3.