Calculate relativistic energies from the Dirac equation
Free particle: E^2 = (pc)^2 + (mc^2)^2
Hydrogen: E = mc^2 / sqrt(1 + (alpha*Z)^2 / (n_r + gamma_j)^2)
gamma_j = sqrt((j + 1/2)^2 - (alpha*Z)^2)
