Q: Integration function spherical coordinates, convolution How can I calculate the following integral explicitly:
$$\int_{R^3}\frac{f(x)}{|x-y|}dx$$
where $f$ is a function with spherical symmetry that is $f(x)=f(|x|)$?
I tried to use polar coordinates at $x=0$ but it didn't help. Any idea on how to do this? Do you think it is doable somehow?

A: Let $$g(y)= \int_{\mathbb{R}^3}\frac{f(x)}{|x−y|}dx$$
Using Fourier theory or the argument here: Can convolution of two radially symmetric function be radially symmetric? $g$ will be spherically symmetric.
So, it is only necessary to find the value of $g$ for $y$ along the north pole: $y = (0,0,r)$; other $y$ with same magnitude $r$ will have the same value by spherical symmetry.
Setting $x=(s\sin(θ)\cos(ϕ),s\sin(θ)\sin(ϕ),s\cos(θ))$, observe that $|x-y| =\sqrt{r^2+s^2-2rs\cos(\theta)}$ for $y$ along the north pole.
So,
$$g(r)=\int_0^\infty s^2\int_0^{\pi}\sin(\theta)\int_0^{2\pi}\frac{f(s)}{\sqrt{r^2+s^2-2rs\cos(\theta)}}d\phi d\theta ds$$
$$\implies g(r)=2\pi\int_0^\infty f(s)s^2\int_0^{\pi}\sin(\theta)\frac{1}{\sqrt{r^2+s^2-2rs\cos(\theta)}} d\theta ds$$
$$\implies g(r)=2\pi\int_0^\infty f(s)\frac{r}{s}(r+s-|r-s|)ds$$.
So, we have reduced the 3-D convolution to a 1-D integral operator with kernel $K(r,s)=2\pi\frac{r}{s}(r+s-|r-s|)$. Obviously, there's no further simplification without knowing the form of $f$.

A: This is a singular integral and so you can expect some weird behavior but if $f$ has spherical symmetry, then I would change to spherical coordinates. Then you'll have $f(x) = f(r)$. $dx$ will become $r^2\sin(\theta)drd\theta d\phi$. The tricky part is then what becomes of $|x-y|$. Recall that $|x-y| = \sqrt{(x-y)\cdot(x-y)} = \sqrt{|x|^2-2x\cdot y+|y|^2}$. In our case, $|x|^2 = r^2$ and $x = (r\sin(\theta)\cos(\phi), r\sin(\theta)\sin(\phi), r\cos(\theta))$. From here, I'm not sure how much simplification there can be. What more are you looking for?
