Computation of NRQCD Parameters in Potential Models

This paper gives details of computing the numerical values of the radial wave function at the origin of heavy quark-antiquark system, jRn`(0)j 2, and its derivatives in potential quark models. These quan-tities provide numerical estimates of the non-perturbative parameters in Non-Relativistic QCD (NRQCD), an eective eld theory describ-ing decay and production of heavy quarkonium states. In particular, we computed wave functions at the origin for three dierent spheri-cally symmetric potentials describing cc, bb and cb bound states. The methods used for calculating the radial wave function at the origin are covered, along with the code used to produce the values given. Comparisons are made with values obtained by Eichten and Quigg [1] and Quigg and Rosner [2].