  Optimized meson operators in the distillation framework are used to study the
charmonium spectrum in two ensembles with two heavy dynamical quarks at half
the physical charm quark mass but different lattice spacings. The use of
optimal meson distillation profiles is shown to increase the overlap with the
ground state significantly, as well as grant access to excited states, for
multiple quantum numbers including hybrid states with very little additional
cost. These same operators are also employed for the calculation of
meson-glueball mixing.
