Gluonic and leptonic decays of heavy quarkonia and the determination of αs(mc) and αs(mb)
Kuang-Ta Chao,Han-Wen Huang,Yu-Quan Liu
DOI: https://doi.org/10.1103/PhysRevD.53.221
IF: 5.407
1996-01-01
Physical Review D
Abstract:The determination of the QCD running coupling constants alpha(s)(m(c)) and alpha(s)(m(b)) is studied with heavy quarkonia c (c) over bar and b (b) over bar decays. The decay rates of V --> 3g and V --> e(+)e(-) for V = J/psi and Upsilon are given in terms of the Bethe-Salpeter amplitudes. To the first-order relativistic correction of on-shell quarks, for the leptonic decay, we have Gamma(V --> e(+)e(-)) = (4 pi alpha(2)e(Q)(2)/m(Q)(2))\integral d(3)q[1-(2 (q) over right arrow(2)/3m(Q)(2))]psi(Sch)((q) over right arrow\(2), which agrees with the NRQCD result, while for the gluonic decay we find Gamma(V --> 3g) = [40(pi(2)-9)alpha(s)(3)(m(Q))/81m(Q)(2)]\ integral d(3)q[1-(2.95 (q) over right arrow(2)/m(Q)(2))]psi(Sch)((q) over right arrow)\(2). Here psi(Sch)((q) over right arrow is the Q (Q) over bar bound-state wave function in momentum space, and m(Q) is the heavy quark mass. This result clearly shows that the relativistic correction (due to the (q) over right arrow(2)/m(Q)(2) term in the decay widths) suppresses the gluonic decay more severely than the leptonic decay. We then estimate these decay widths by further including the first-order QCD radiative corrections (given in the <(MS)over bar> scheme and at the heavy quark mass scale) on the basis of the factorization assumption, and using the meson wave functions which are obtained with a QCD-inspired interquark potential. Using the experimental values of the ratio R(g) = Gamma(V --> 3g)/Gamma(V --> e(+)e(-)) approximate to 10, 32 for V = J/psi, Upsilon, respectively, and the calculated widths, we find alpha(s)(m(c)) = 0.26-0.29 and alpha(s)(m(b)) = 0.19-0.21 at m(c) = 1.5 GeV and m(b) = 4.9 GeV. These values for the QCD running coupling constant are substantially enhanced, as compared with the ones obtained without relativistic corrections, and are potentially consistent with the QCD scale parameter Lambda(<(MS)over bar>)((4)) approximate to 200 MeV. We emphasize, however, that our numerical results of the running coupling constant mainly serve as an improved estimate rather than a precise determination for which the existing theoretical uncertainties due to higher order relativistic corrections and the scheme dependence of the radiative corrections should be further clarified, and a first principles estimate of the nonperturbative bound-state effects should be further studied.