+The $E_l^{\text{SO,KB}}$ are given by
+\begin{equation}
+E_l^{\text{SO,KB}}=
+\frac{\braket{\delta V_lu_l}{u_l\delta V_l}}
+ {\bra{u_l}\delta V_l\ket{u_l}}=
+\frac{\int_{r}\delta V^2_l(r)u^2_l(r)}r^2dr
+ {\int_{r'}\int_{r''}\braket{u_l}{r'}\bra{r'}\delta V_l
+\ket{r''}\braket{r''}{u_l}}=
+\end{equation}