\hat{H}=-\frac{\hbar^2}{2m}\frac{\partial^2}{\partial t^2}+q\phi(\vec{r})