+P'_l\left(\frac{\vec{r}'\vec{r}''}{r'r''}\right)\cdot
+\left(\frac{\vec{r}'\times\vec{r}''}{r'r''}\right)\cdot
+\frac{E^{\text{SO,KB}}_l\delta V_l^{\text{SO}}(r'')\frac{u_l(r'')}{r''}}
+ {\int_{r}{\delta V_l^{\text{SO}}}^2(r) u_l^2(r) dr} \,
+\frac{2l+1}{4\pi}\text{ ,}\nonumber\\
+\label{eq:solid:so_fin}
+\end{eqnarray}
+since derivatives of functions that depend on the absolute value of $\vec{r}'$ do not contribute due to the cross product as is illustrated below (equations \eqref{eq:solid:rxp1} and \eqref{eq:solid:rxp2}).
+\begin{eqnarray}
+\left(\vec{r}\times\nabla_{\vec{r}}\right)f(r)&=&
+\left(\begin{array}{l}
+r_y\frac{\partial}{\partial r_z}f(r)-r_z\frac{\partial}{\partial r_y}f(r)\\
+r_z\frac{\partial}{\partial r_x}f(r)-r_x\frac{\partial}{\partial r_z}f(r)\\
+r_x\frac{\partial}{\partial r_y}f(r)-r_y\frac{\partial}{\partial r_x}f(r)
+\end{array}\right)
+\label{eq:solid:rxp1}
+\end{eqnarray}
+\begin{eqnarray}
+r_i\frac{\partial}{\partial r_j}f(r)-r_j\frac{\partial}{\partial r_i}f(r)&=&
+r_if'(r)\frac{\partial}{\partial r_j}(r_x^2+r_y^2+r_z^2)^{1/2}-
+r_jf'(r)\frac{\partial}{\partial r_i}(r_x^2+r_y^2+r_z^2)^{1/2}\nonumber\\
+&=&
+r_if'(r)\frac{1}{2r}2r_j-r_jf'(r)\frac{1}{2r}2r_i=0
+\label{eq:solid:rxp2}