We saw that the Dirac equation predicted that there is interaction between the electron spin and the magnetic field, SB, with strength μB=2meℏc. When the electron has angular momentum L, such as in an atomic orbital, there is also a BL interaction and a spin-orbit coupling SL. The Dirac equation (along with symmetry arguments) predicts the strength of all three interactions, as well as other corrections. To see the effect of these terms on the hydrogen atom, we have to take the non-relativistic limit.
(a) For the Schr\"odinger equation, we need the Hamiltonian, not the Lagrangian. Find the Dirac Hamiltonian by writing the Dirac equation as i∂tψ=HDψ. Write the Hamiltonian in terms of momenta pi rather than derivatives ∂i.
(b) Calculate (HD+eA0)2 in the Weyl representation for ψ=(ψLψR). Leave in terms of σi, pi and Ai. Put back in the factors of c and ℏ, keeping the charge e dimensionless.
(c) Now take the square root of this result and expand in c1, subtracting off the zero-point energy mc2, i.e. compute H=HD−mc2 to order c0. Looking at the σi term, how big are the electron's electric and magnetic dipole moments?
(d) The size of the terms in this Hamiltonian are only meaningful because the spin and angular momentum operators have the same normalization. Check the normalization of the angular momentum operators Li=εijkxjpk and the spin operators Si=21σi by showing that they both satisfy the rotation algebra: [Ji,Jj]=iεijkJk.
(e) The gyromagnetic ratio, ge (sometimes called the g-factor), is the relative size of the SB and LB interactions. Choose a constant magnetic field in the z direction, then isolate the BzLz coupling in H. Extract the electron gyromagnetic ratio ge by writing the entire coupling to the magnetic field in the Hamiltonian as μBBz(Lz+geSz)=Bzμz, with μ≡μB(L+geS). How could you experimentally measure ge (e.g. with spectroscopy of the hydrogen atom)?
(f) In spherical coordinates, the Schr\"odinger equation has a L2 term. With spin, you might expect that this becomes Lμ=μB(L2+geLS), making the LS term proportional to the g-factor. This is wrong. It misses an important relativistic effect, Thomas precession. It is very hard to calculate directly, but easy to calculate using symmetries. With no magnetic field, the atom, with spin included, is still rotationally invariant. Which of J=L+S or μ=L+geS is conserved (i.e. commutes with H)? Using this result, how does the spin-orbit coupling depend on ge?
(g) There are additional relativistic effects coming from the Dirac equation. Expand the Dirac equation to next order in c21, producing a term that scales as p4.
(h) Now let us do some dimensional analysis -- there is only one scale me. Show that the electron's Compton wavelength, the classical electron radius, re, the Bohr radius, a0, and the inverse-Rydberg constant, Ry−1, are all me times powers of αe. Are the splittings due to the p4 term fine structure (ΔE∼αe2E), hyperfine structure (ΔE∼αe4E) or something else? [Hint: write out a formula for the energy shift using time-dependent perturbation theory, then see which of the above length scales appears.]