Solution (source code)

= Solution

$L\cdot S=\tfrac12(J^2-L^2-S^2)$. Hence the coupled $|n\ell sjj_z\rangle$ <basis> diagonalizes the spin-orbit term.