Prediction of the bottomonium D-wave spectrum from full lattice QCD
  We calculate the full spectrum of D-wave states in the Upsilon system in lattice QCD for the first time, using an improved version of NonRelativistic QCD on coarse and fine "second generation" gluon field configurations from the MILC collaboration that include the effect of up, down, strange and charm quarks in the sea. Taking the 2S-1S splitting to set the lattice spacing, we determine the ^3D_2 - 1S splitting to 2.3%, and find agreement with experiment. Our prediction of the fine structure relative to the ^3D_2 gives the ^3D_3 at 10.181(5) GeV and the ^3D_1 at 10.147(6) GeV. We also discuss the overlap of ^3D_1 operators with ^3S_1 states. 
J. O. Daldrop, C. T. H. Davies, R. J. Dowdall
hep-lat hep-ph