  We compute the energy levels of some of the lower-lying heavy quarkonium
states perturbatively up to O(alpha_s^5*m) and O(alpha_s^5*m*log[alpha_s]).
Stability of the predictions depends crucially on the unknown 4-loop pole-MSbar
mass relation. We discuss the current status of the predictions with respect to
the observed bottomonium spectrum.
