QCD radiative corrections to physical quantities are studied using Stevenson's principle of minimal sensitivity (PMS) to define the renormalization. We examine several naive potentials (Cornell group, power law and logarithmic), as well as the more sophisticated Richardson model in order to determine the spectra for the non-relativistic heavy charmonium and bottomonium systems. Predictions are made for the values of hyperfine splittings, leptonic and hadronic decay widths and E1 transition rates for these families of mesons. It is shown that good agreement with experimental data may be achieved by using a constant value of ..lambda../sub QCD/, which is determined by the PMS scheme and the potential model.