Fourier transform of the Coulomb potential When trying to find the Fourier transform of the Coulomb potential
$$V(\mathbf{r})=-\frac{e^2}{r}$$
one is faced with the problem that the resulting integral is divergent. Usually, it is then argued to introduce a screening factor $e^{-\mu r}$ and take the limit $\lim_{\mu \to 0}$ at the end of the calculation.
This always seemed somewhat ad hoc to me, and I would like to know what the mathematical justification for this procedure is. I could imagine the answer to be along the lines of: well, the Coulomb potential doesn't have a proper FT, but it has a weak FT defined like this …
Edit: Let me try another formulation:
What would a mathematician (being unaware of any physical meanings) do when asked to find the Fourier transform of the Coulomb potential?
 A: I am personally fond of this short and sweet argument.
If you believe the relatively easy to prove fact that in 3 dimensions
$$
\nabla^2\frac{1}{r} = -4\pi \delta^{(3)}(\bf{r})
$$
then taking the transform of both sides gives
$$
-q^2\widehat{\frac{1}{r}} = - 4\pi
\quad\longrightarrow\quad
\widehat{\frac{1}{r}} = \frac{4\pi}{q^2}
$$
A: If the final result is unique and independent of the regularization method, then it is a mathematically correct answer.
A: I really appreciate the physical explanations made in other answers, but I want to add that Fourier transform of the Coulomb potential makes mathematical sense, too. This answer is meant to clarify on what sense the standard calculation is valid mathematically.
Firstly, and maybe more importantly, I want to emphasize that
The Fourier transform of f is not simply just $\int{f(x)e^{-ikx}dx}$.
For an $L^1$ function (a function which is norm integrable), this is always the case but Coulomb potential is definitely not in $L^1$. So the Fourier transform of it, if it ever exists, is not expected to be the integral above.
So here comes the second question: can Fourier transformation be defined on functions other than $L^1$?
The answer is "yes", and there are many Fourier transformations. Here are two examples.

*

*Fourier transformation on $L^2$ functions (i.e., square integrable functions).

It turns out that the Fourier transform behaves more nicely on $L^2$ than on $L^1$, thanks to the Plancherel's theorem. However, as we mentioned above, if an $L^2$ function is not in $L^1$, then the above integral may not exist and Fourier transform is not given by that integral, either. (However, it has a simple characterization theorem, saying that in this case the Fourier transform is given by the principle-value integration of the above integral.)


*Fourier transform of distributions (generalized functions)
It is in this sense that the Forier transform of Coulomb potential holds. The Coulomb potential, although not an $L^1$ or $L^2$ function, is a distribution. So we need to use the definition of the Fourier transform to distributions in this case. Indeed, one can check the definition and directly calculate the Fourier transform of it. However, the physicists' calculation illustrates another point.
Fourier transformation on distributions (however it is defined) is continuous (under a certain topology on the distribution space, but let's not be too specific about it).
Remember that if f is continuous, then $x_\epsilon\rightarrow x$ implies that $f(x_\epsilon)\rightarrow f(x)$.
Now $\frac{1}{r}e^{-\mu r}\rightarrow\frac{1}{r}$ when $\mu\rightarrow 0$ (again, under the "certain topology" mentioned above), and therefore continuity implies
$$\operatorname{Fourier}\left\{\frac{1}{r}e^{-\mu r}\right\}\rightarrow \operatorname{Fourier}\left\{\frac{1}{r}\right\}.$$
However, $\frac{1}{r}e^{-\mu r}$ is in $L^1$ and therefore its Fourier transformation can be computed using the integral $\int{f(x)e^{-ikx}dx}$.
Therefore those physicists' computations make perfect mathematical sense, but it's on Fourier transform of distributions, which is much more general than that on $L^1$ functions.
Wish this answer can build people's confidence that the Coulomb potential Fourier transformation problem is not only physically reasonable but also mathematically justifiable.
A: This is not really a complete answer, but more of an interesting observation which might help to understand the regularized result better. First observe that
\begin{equation}
\int_V d^3 r \frac{e^{i\mathbf q \cdot \mathbf r}}{r} = 2\pi \int_0^R dr \int_{-1}^1 du \, r e^{iqru} = \frac{4\pi}{q} \int_0^R dr \, \sin qr = \frac{4 \pi}{q^2} \left( 1 - \cos qR \right).
\end{equation}
In the limit $R \rightarrow \infty$ the integral diverges since it fluctuates between $0$ and $8\pi/q^2$ indefinitely. The regularized Fourier transform is equal to the average value of these extremes. This is similar to the divergent series $1-1+1-1+\cdots$ which can be regularized to $1/2$.
A: I think you need this integral
$$\int {\frac{1}{{\left( {2\pi } \right)^3 }}\frac{{e^{i\vec k \cdot \vec r} }}{{k^2  + a^2 }}d^3 \vec k}  = \frac{{e^{ - ar} }}{{4\pi r}}$$
We suppoes initially that $a \ne 0$. If we use polar coordinates  with $\vec k$ in the z-axis we have:
$$\int {\frac{1}{{\left( {2\pi } \right)^3 }}\frac{{e^{i\vec k \cdot \vec r} }}{{k^2  + a^2 }}d^3 \vec k}  = \frac{1}{{\left( {2\pi } \right)^3 }}\int\limits_0^\infty  {\frac{{k^2 }}{{k^2  + a^2 }}dk} \int\limits_0^{2\pi } {d\varphi } \int\limits_{ - \pi /2}^{\pi /2} {\cos \theta e^{ikr\sin \theta } } d\theta  = $$
$$ = \frac{1}{{\left( {2\pi } \right)^2 }}\int\limits_0^\infty  {\frac{{k^2 }}{{k^2  + a^2 }}dk} \int\limits_{ - 1}^1 {e^{ikrx} dx}  = \frac{2}{{\left( {2\pi } \right)^2 }}\int\limits_0^\infty  {\frac{{k^2 }}{{k^2  + a^2 }}\frac{{i\sin \left( {kr} \right)}}{{ikr}}dk} $$
$$ = \frac{1}{{\left( {2\pi } \right)^2 }}\frac{1}{r}\int\limits_{ - \infty }^\infty  {\frac{{k\sin \left( {kr} \right)}}{{k^2  + a^2 }}} dk$$
It is the immaginary part of:
$$  \frac{1}{{\left( {2\pi } \right)^2 }}\frac{1}{r}\int\limits_{ - \infty }^\infty  {\frac{{ke^{ikr} }}{{k^2  + a^2 }}} dk = \frac{1}{{\left( {2\pi } \right)^2 }}\frac{1}{r}2\pi i\frac{{iae^{ - ar} }}{{2ia}} = \frac{1}{{4\pi r}}e^{ - ar} i$$
so:
$$\int {\frac{1}{{\left( {2\pi } \right)^3 }}\frac{{e^{i\vec k \cdot \vec r} }}{{k^2  + a^2 }}d^3 \vec k}  = \frac{1}{{4\pi r}}e^{ - ar} $$
The unique singular point inside the circumference is $k=ia$. 
We now determine the case $a=0$. We have:
$$\int {\frac{1}{{\left( {2\pi } \right)^3 }}\frac{{e^{i\vec k \cdot \vec r} }}{{k^2 }}d^3 \vec k}  = \frac{1}{{\left( {2\pi } \right)^2 }}\frac{1}{r}\int\limits_{ - \infty }^\infty  {\frac{{\sin \left( {kr} \right)}}{k}} dk = $$
$$\frac{1}{{\left( {2\pi } \right)^2 }}\frac{1}{r}\int\limits_{ - \infty }^\infty  {\frac{{\sin t}}{t}dt = \frac{1}{{4\pi r}}} $$
A: I'm not sure that a justification is required. 
We don't feel the need to justify solutions obtained by making a change of variables, nor those obtained by differentiating under the integral sign, nor those obtained by analytic continuation and application of contour integrals and so on for a host of other tricks of symbol manipulations.
But if you insist...
You can interpret this as the potential due to a massive scaler boson (the Yukawa potential), and then take the limit of vanishing mass.
A: Dear Emerson, the "weak Fourier transform" of the potential, which is simply $4\pi / q^2$, is the "proper physical" value of the Fourier transform.
In physics, it simply makes no sense to say that a Fourier transform "doesn't exist": the Coulomb potential clearly does exist, and when we're calculating something that depends on its Fourier transform, the answer also has to exist because the world has to behave in some way. It's just generally impossible in physics to say that a valid theory produces an "ill-defined answer". We must always struggle to get an actual, well-defined answer.
And quantum electrodynamics, even though it has a superficially divergent Fourier transform of the Coulomb potential, is an excellent theory. Still, we need the Fourier transform to answer many questions. So the only question can be how to find out the right answer - not whether an answer exists. Of course that it does. The physically correct answer is most easily obtained by regulating the potential - as a limit of the Yukawa potential - and by exchanging the limit with the integration in the middle of the calculation, to guarantee that the result is sensible.
This is like giving a mass to the photon - we call it an infrared regulator. Such things are omnipresent in quantum field theory. If an expression doesn't naively make sense, one tries to express it as a limit of a more general expression, and calculate a more general calculation, which is more well-defined, and try to take the required limit at the very end. This is the right physicist's approach to seemingly ill-defined expressions and it always leads us out of the trouble and it always tells us something more meaningful.
A failure - a refusal to calculate the answers - is just not acceptable in physics.
