Determining hybrid content of heavy quarkonia using lattice nonrelativistic QCD

Using lowest-order lattice NRQCD to create heavy meson propagators and applying the spin-dependent interaction, $c_B^{} \frac{-g}{2m_q}\vec\sigma\cdot\vec{B}$, at varying intermediate time slices, we compute the off-diagonal matrix element of the Hamiltonian for the quarkonium-hybrid two-state system. Diagonalizing this two-state Hamiltonian, the admixture of hybrid ($|Q\bar{Q}g>$) in the ground state is found. We present results from a set of quenched lattices with an interpolation in quark mass to match the bottomonium spectrum.