\hat{H}=\sum_{q=x,y,z}(\frac{\hat{p}^{2}_{q}}{2m}+\frac{1}{2}m\omega^{2}_{q}\hat{q}^2)