I am doing a computation in quantum field theory and the following integral occurred to me $$ I(a)=\int_{-\infty}^{+\infty}\frac{e^{-a\sqrt{x^2+1}}dx}{x^2+1} \qquad a\ge 0. $$ I would like to know a closed form for it and the eventual steps to achieve the result.
Thanks.
Let me denote your integral by $I(a)$. Making the change of variables $x=\sinh\varphi$, we get $$I(a)=\int_{-\infty}^{\infty}\frac{e^{-a\cosh\varphi}d\varphi}{\cosh\varphi}.$$ Next, differentiating this expression with respect to $a$, we get $$I'(a)=-\int_{-\infty}^{\infty}e^{-a\cosh\varphi}d\varphi=- 2K_0(a),$$ where $K_0(a)$ denotes the Macdonald function. Integrating back, we find that $$I(a)=I(0)-2\int_0^aK_0(x)dx=\pi-2\int_0^aK_0(x)dx.$$ According to Prudnikov-Brychkov-Marychev (Vol.2, formula 1.12.1.4), the remaining integral is expressed in terms of Macdonald and modified Struve functions, so that $$I(a)=\pi-2a\left\{K_0(a)+\frac{\pi}{2}\left[K_0(a)\mathbf{L}_1(a)+K_1(a)\mathbf{L}_0(a)\right]\right\}.$$ Presumably this answer does not add anything useful to your integral except that the corresponding special functions have a name, but at least one can be sure that it cannot be simplified further.