The Hamiltonian I am considering is $$H=\vec{\alpha} \cdot \vec{p} + \beta m,$$ where ##\alpha^i = \gamma^0\gamma^i##, ##\beta=\gamma^0##, ##p_k = i\partial_k##, ##\vec{p} = -i\vec{\partial}##, so that ##i\partial_t \psi=H\psi## gives just the Dirac equation ##(i\gamma^\mu\partial_\mu -m...