Suppose the matrix elements of the perturbation Hamiltonian H’ = -eE . r are
since \(\langle l'|H'|l\rangle\) = 0 for I', I states of the same parity. Then for energy level P3/2, we have
For energy levels P1/2 and S1/2, we diagonalize the Hamiltonian in the corresponding subspace, i.e, solve
The roots are \(\lambda\) = \(\pm|a|\), which give the new wave functions