Solution of Definite integral:$\int_{-\infty}^{\infty} \int_{-\infty}^{\infty} \frac{1}{\sqrt{x^{2}+y^{2}+z^{2}}}e^{i(k_{1}x+k_{2}y)}dxdy$ I'm trying to evaluate the following two dimensional integral:
$\int_{-\infty}^{\infty} \int_{-\infty}^{\infty} \frac{1}{\sqrt{x^{2}+y^{2}+z^{2}}}e^{i(k_{1}x+k_{2}y)}dxdy$
The paper that i'm following reports the following solution:
$\int_{-\infty}^{\infty} \int_{-\infty}^{\infty} \frac{1}{\sqrt{x^{2}+y^{2}+z^{2}}}e^{i(k_{1}x+k_{2}y)}dxdy = \frac{1}{\sqrt{k_{1}^{2}+k_{2}^{2}}}e^{-\sqrt{k_{1}^{2}+k_{2}^{2}}z}$
it might be useful using the following definitions:
$r = \sqrt{x^{2}+y^{2}+z^{2}}$
and
$k = \sqrt{k_{1}^{2}+k_{2}^{2}}$
I think this should follow from some kind of substitution;
anyway i tried to solve it with mathematica but i'm not getting any solution...
Any idea on how to obtain the result shown by the authors?
Thanks in advance
 A: Ok let's see:
Denote the Integral by $I(k_x,k_y,z)$
Going over to Polar coordinates we have:
$x=r\sin(\phi),r\cos(\phi)$, $|det(J)|=r$
\begin{align}
I(k_x,k_y,z)=\int_0^{\infty} dr \frac{r}{\sqrt{r^2+z^2}}\int_0^{2\pi} d\phi e^{-ir[k_x sin(\phi)+k_y cos(\phi)]}
\end{align}
It can now be shown that the inner integral equals $2\pi J_0[r\sqrt{k_x^2+k_y^2}]$
So we are left with (1)
\begin{align}
I(k_x,k_y,z)=2\pi\int_0^{\infty} dr \frac{r}{\sqrt{r^2+z^2}}J_0[r\sqrt{k_x^2+k_y^2}]
\end{align} 
Using now Formula 6.554.1 from Gradstheyn/Ryzhik which states (2)
that $\int_0^{\infty}dq\frac{q}{\sqrt{s^2+q^2}}J_0[qp]=\frac{e^{-sp}}{p}$
We get 
\begin{align}
I(k_x,k_y,z)=\frac{2\pi}{\sqrt{k_x^2+k_y^2}}e^{-z\sqrt{k_x^2+k_y^2}}
\end{align}
which is the stated result apart from a factor of $2\pi$ which supposedly stems from some on of the definitions of Fouriertransform  
Unluckily i'm very busy so i can't fill in the proof which leds us to (1) or try to proof (2) so i would appreciate if someone could fill this gaps
A: Let me fill in the details of how to perform the angular integration in tired's answer. First, we can make our lives easier by notating that the integrand $e^{i \mathbf{k}\cdot\mathbf{x}}$ with $\mathbf{k}=k_1 \hat{x}+k_2\hat{y}=k\hat{k}$ is invariant under rotations about the $z$-axis. So we're free to orient the coordinates such that $\hat{k}$ as the positive $x$-axis i.e. $\mathbf{k}\mapsto k\hat{x}$. Then, defining $s=\sqrt{x^2+y^2}$, the integral becomes
$$I(k,z)=\int_{-\infty}^{\infty} \int_{-\infty}^{\infty} \frac{1}{\sqrt{x^{2}+y^{2}+z^{2}}}e^{i k x}\,dxdy=\int_{0}^{\infty}\int_{-\pi}^{\pi}\frac{e^{ik s\cos\theta }}{\sqrt{s^2+z^2}}s\,d\theta \,ds.$$
While this integral hardly seems better, the $\theta$-integral is recognizable to Mathematica in terms of a Bessel function. Indeed, if uses the integral representation $J_n(x)=\int\limits_{-\pi}^{\pi} e^{i (n \theta -x\sin \theta)}\,d\theta$ (see eq. 71 on Mathworld's Bessel functions page, or derive by hand using an Jacobi-Anger identity) then the $n=0$ case allows us to write the integration as
$$I(k,z)=2\pi\int_{0}^{\infty}\frac{s\,J_0(ks)}{\sqrt{s^2+z^2}}\,ds=\frac{2\pi}{k}\int_{0}^{\infty}\frac{u\,J_0(u)}{\sqrt{u^2+k^2 z^2}}\,du\quad\text{with }u=ks.$$
What remains is to justify the integral identity cited by tired which allows us to conclude $I(k,z)=2\pi k^{-1} e^{-k |z|}.$ I'll see if I can recall the details.
A: In my first answer, I proceeded by trying to directly compute the momentum integrals. In doing so, I'd forgotten that there is a rather less direct approach which yet has the advantage of not requiring Bessel function identities.
First, note that our task is to take the $2$D Fourier transform of the function $G(\mathbf{r})=(x^2+y^2+z^2)^{-1/2}=r^{-1}$. My choice of $G$ as symbol for this function is not done idly, for this is in fact the free-space Green's function of the $3$D  Laplacian i.e. $\nabla^2 G(\mathbb{r})=-4\pi \delta^{3}(\mathbf{r})$. This is a standard Dirac-delta identity which can be verified via the divergence theorem and direct differentiation.
The $3$D Fourier transform of this is readily computed to be $(k_x^2+k_y^2+k_z^2)G(\mathbf{k})=4\pi.$ To get the planar Fourier transform, we can then take the inverse transform with respect to $k_z$. That iss,
\begin{align}
G(k_x,k_y,z)=\mathcal{F}^{-1}[G(\mathbf{k})](z)
&=\frac{1}{2\pi}\int_{-\infty}^\infty e^{-i k_z z}G(\mathbf{k})\,dk_z\\
&=\int_{-\infty}^\infty \frac{2e^{-i k_z z}}{k_x^2+k_y^2+k_z^2}\,dk_z\\
&=\frac{2}{k}\int_{-\infty}^\infty \frac{e^{-i k z u}}{1+u^2}\,du
\end{align}
where in the last line I have introduced $k=\sqrt{k_x^2+k_y^2}$ and $k_z=k u$. This last integral can be evaluated via methods of contour integration (see Wikipedia of such) to obtain the desired result
$$\boxed{G(k_x,k_y,z)=2\pi\dfrac{e^{-k |z|}}{k}=2\pi\left(\dfrac{e^{-|z|\sqrt{k_x^2+k_y^2}}}{\sqrt{k_x^2+k_y^2}}\right)}$$
