NRQCD and Static Systems -- A General Variational Approach
  We present initial results from Monte Carlo simulations of NRQCD-light, static-light, and NRQCD-NRQCD mesons, using a variational technique (MOST), as part of our ongoing calculation of the $f_{B}$ decay constant. The basis states for the variational calculation are quark-antiquark operators separated by all possible relative distances not equivalent under the cubic group (for example, for a $20^{3}$ lattice there are 286 operators). The efficacy of the method is demonstrated by the good plateaus obtained for the ground state and the clean extraction of the wave functions of the ground and first radially excited state. 
Terrence Draper, Craig McNeile and Constantine Nenkov