= \left( \hat{\mathbf{p}} -\frac{e}{c} \mathbf{A} \right)^2 + i \hat{\boldsymbol{\sigma}} \cdot