Help with the specific improper integral for finding galaxy potential

I'm working on finding the potential of a razor-thin disk for a galaxy. I have encountered a specific improper integral for which I'm looking for an analytical expression. The integral is given below:

$$\Phi(R, z) = -2\pi G \Sigma_0\alpha\int_0^\infty e^{-k|z|} J_0(kR)\left(\alpha^2 + k^2 \right)^{(-3/2)} dk\,,$$

where $$J_0$$ is the Bessel function of the first kind and $$\alpha$$ is some constant.

I'm looking for a solution expressed in terms of any special functions that might appear, no matter how complex they are, but not in terms of power series. I tried using WolframAlpha and couple of other online tools, but none of them find any solution. Could somebody help or point me out to some powerful tool for evaluating integrals like this?

• With binomial expansion of $(a^2+k^2)^{-\frac32}$ and the series form of $J_0$, you may find the Cauchy of two series and then apply Laplace transform each term $x^k$. – Nosrati Dec 7 '18 at 10:23
• This integral does not seem to have a valid solution for any $R$ at $z=0$. Is that what you are expecting from a physical point of view? – James Arathoon Dec 9 '18 at 12:00