H=\frac{1}{2m}(\vec{p}-q\vec{A})^2+\mu_B\vec{\sigma}\cdot\vec{B}+q\phi