\hat{H}=\frac{\vec{\hat{p}}^2}{2m}+\frac{1}{2}m\omega^2\vec{\hat{r}^2}