No matter what I can't seem to arrive at this answer. I've tried $\partial_{i}\left (\partial_{i}\left ( x_{i}x_{j}r^{n} \right ) \right )$ and $\delta _{ij}\partial_{i}\left (\partial_{j}\left ( x_{i}x_{j}r^{n} \right ) \right )$ and I never manage to get the final answer.
Can someone help?
Thanks!
First recall $\partial_i r=x_ir^{-1}$, so $\partial_ir^m=mx_ir^{m-2}$. This will come in handy in our calculation. I'll also assume you are working with $\mathbb{R}^3$, otherwise $\delta_{kk}$ will not be $3$.
We compute $$\require{cancel} \begin{align*} \nabla^2(x_ix_jr^n) &=\partial_k\partial_kx_ix_jr^n\\ &=\partial_k((\partial_kx_i)x_jr^n+x_i(\partial_kx_j)r^n+x_ix_j(\partial_kr^n))\\ &=\partial_k(\delta_{ik}x_jr^n+\delta_{jk}x_ir^n+nx_ix_jx_kr^{n-2})\\ &=\partial_ix_jr^n+\partial_jx_ir^n+n\partial_kx_ix_jx_kr^{n-2}\\ &=\delta_{ij}r^n+x_j\partial_ir^n+\delta_{ji}r^n+x_i\partial_jr^n\\ &\quad+n\left[\delta_{ik}x_jx_kr^{n-2}+x_i\delta_{jk}x_kr^{n-2}+x_ix_j\cancelto{3}{\color{red}{\delta_{kk}}}r^{n-2}+x_ix_jx_k\partial_kr^{n-2}\right]\\ &=\delta_{ij}r^n+nx_ix_jr^{n-2}+\delta_{ij}r^n+nx_ix_jr^{n-2}\\ &\quad+n\left[x_ix_jr^{n-2}+x_ix_jr^{n-2}+3x_ix_jr^{n-2}+(n-2)x_ix_j\cancelto{r^2}{\color{red}{x_kx_k}}r^{n-4}\right]\\ &=2\delta_{ij}r^n+n(n+5)x_ix_jr^{n-2}\\ \end{align*} $$