H=\frac{1}{2m}\left(\vec{p}-q\vec{A}\right)^2+e\phi