  When the cancellation of the leading renormalon contributions is
incorporated, the total energy of a b bbar system, E_{tot,bbbar}(r) = 2
m_{pole,b} + V_QCD(r), agrees well with the potentials used in phenomenological
models for heavy quarkonia in the range 0.5 GeV^-1 < r < 3 GeV^-1. We provide a
connection between the conventional potential-model approaches to the
quarkonium spectroscopy and the recent computation based on perturbative QCD.
