  The vNRQCD Lagrangian for colored heavy scalar fields in the fundamental
representation of QCD and the renormalization group analysis of the
corresponding operators are presented. The results are an important ingredient
for renormalization group improved computations of scalar-antiscalar bound
state energies and production rates at next-to-next-to-leading-logarithmic
(NNLL) order.
