  Using the Bodwin--Braaten--Lepage factorization theorem in heavy quarkonium
decay and production processes, we calculated matrix elements associated with
S- and P-wave bottomonium decays via lattice QCD simulation methods. In this
work, we report preliminary results on the operator matching between the
lattice expression and the continuum expression at one loop level.
Phenomenological implications are discussed using these preliminary
\(\overline{MS}\) matrix elements.
