Consider a neutral spin one-half Dirac particle with mass and with an intrinsic magnetic moment, and assume the Hamiltonian
$$
H=c \boldsymbol{\alpha} \cdot \frac{\hbar}{i} \nabla+\beta m c^2+\lambda B \beta \Sigma_z
$$
in the presence of a uniform constant magnetic field along the $z$ axis. Determine the important constants of the motion, and derive the energy eigenvalues. Show that orbital and spin motions are coupled in the relativistic theory but decoupled in a nonrelativistic limit. The coefficient $\lambda$ is a constant, proportional to the gyromagnetic ratio.