The t decay rate into hadrons of invariant mass smaller than p s 0 ¿ L QCD can be calculated in QCD assuming global quark-hadron duality. It is shown that this assumption holds for s 0 . 0.7 GeV 2 . From measurements of the hadronic mass distribution, the running coupling constant a s ͑s 0 ͒ is extracted in the range 0.7 GeV 2 , s 0 , m 2 t . At s 0 m 2 t , the result is a s ͑m 2 t ͒ 0.329 6 0.030. The running of a s is in good agreement with the QCD prediction.PACS numbers: 12.38. Aw, 11.10.Hi, 12.38.Qk, 13.35.Dx The scale dependence of coupling constants is one of the key features of renormalizable quantum field theories. In QCD, the effective coupling constant a s ͑Q 2 ͒ is predicted to decrease with the momentum transfer Q 2 , a property referred to as asymptotic freedom [1]. This prediction has been tested by comparing data obtained from experiments operating at different energies [2]; it has also been studied in single high-energy experiments at ep and pp colliders, where a large range in Q 2 can be probed simultaneously [3]. Here we propose a test of the scale dependence of a s ͑Q 2 ͒ in the low-energy regionOur method is based on integrals of the invariant mass distribution in hadronic t decays. It provides a unique opportunity to test one of the most important predictions of QCD in a single experiment and at low energies, where the effect of the running of a s is most pronounced.We shall consider the t decay rate into hadrons of invariant mass squared smaller than s 0 , normalized to the leptonic decay rate,where dR t ͞ds is the inclusive hadronic spectrum. As long as s 0 ¿ L 2 QCD , the quantity R t ͑s 0 ͒ can be calculated in QCD using the operator product expansion (OPE) [4,5]. Applying the OPE in the physical region assumes global quark-hadron duality, i.e., that decay rates admit a QCD description after a "smearing" over a sufficiently wide energy interval has been performed [6], which in the present case is provided by the integration over the range 0 , s , s 0 . The question of how accurate this assumption is and for what values of s 0 it applies is a phenomenological one; it cannot be answered yet from theoretical grounds. Below, we shall investigate this question, comparing data with theoretical predictions based on the duality assumption. A similar test of duality has been performed in Ref. [7], using data on the e 1 e 2 ! hadrons cross section.The t decay rate into hadrons can be written in terms of moments M ͑J͒ k of the absorptive part of current-current correlation functions of angular momentum J [8,9]. The quantity R t ͑s 0 ͒ is given bywhere S EW Ӎ 1.0194 accounts for electroweak radiative corrections [10]. The moments can be written as contour integrals along a circle of radius s 0 in the complex plane. Since the only large mass scale in these integrals is s 0 , the OPE provides an expansion in powers of 1͞s 0 ,The leading term is given by perturbation theory alone. Terms suppressed by powers of 1͞s 0 consist of perturbative coefficients c ͑J͒ n multiplying dimensionful parameters ͗O...