  We investigate the effect of a non-perturbative gluon condensate on the
spectrum of bound states of extremely heavy quarks using Non-Relativistic-QCD
(NRQCD). Simulations are done with an improved gluon action on small lattices.
We compare our numerical results with perturbative calculations to test the
consistency.
