Los armónicos sólidos fueron explicados por Misha, así que déjame completar el resto de los detalles.
El operador laplaciano viene dado por $\Delta = \partial_x^2 + \partial_y^2 + \partial_z^2$ . En primer lugar, supongamos que sólo nos interesa la parte radial. Utilizando la regla de la cadena (y dejando que $D \equiv \partial_r$ como en sus referencias), podemos escribir $$\partial_x = r_{,x} D, \quad \partial_x^2 = r_{,xx} D + (r_{,x} D)^2.$$ Tenemos que calcular $$ r_{,x} = {\partial \sqrt{x^2 + y^2 + z^2} \over \partial x } = {x \over r} $$ y $$ r_{,xx} = \partial_x {x \over r} = {1 \over r} - {x^2 \over r^3}.$$ Las expresiones para $\partial_y$ y $\partial_z$ son, por supuesto, similares, por lo que $$\Delta = ({3 \over r} - {r^2 \over r^3})D + D^2 = {2 \over r}D + D^2.$$
Ahora, la parte esférica del operador laplaciano viene dada por $-{L^2 \over r^2}$ donde $\mathbf L$ es el operador de momento angular. Si utilizamos armónicos esféricos de nivel $l$ (que son estados propios de L^2 correspondientes al valor propio $l(l+1)$ ) y hacer una sustitución $\phi \to {\phi \over r}$ obtenemos el resultado de Griffith.