I am interested in the value of
$$\int_0^\infty e^{-\alpha t}\frac{e^{-\frac{|x-y|^2}{2t}}}{\sqrt{2\pi t}}\, dt $$
this is the laplace transform of the Heat kernel (changing the time variable)
This integral came up in this exercise:
Exercise 3.7. Show that the resolvent for Brownian motion is given by $$U(\alpha)f(x)=\frac{1}{\sqrt{2\alpha}}\int_{-\infty}^\infty f(y)e^{-\sqrt{2\alpha}|x-y|}dy.$$
where
$$U(\alpha)f=\int_0^\infty e^{-\alpha t}T(t)f\,dt,\quad\alpha>0,\tag{3.6}$$
First setps:
\begin{align} U(\alpha)f(x) &= \int_0^\infty e^{-\alpha t} T(t) f (x)\, dt \\ &= \int_0^\infty e^{-\alpha t} \int_{-\infty}^\infty \frac{e^{-\frac{|x-y|^2}{2t}}}{\sqrt{2\pi t}} f (y)\, dy\, dt \\ &= \int_{-\infty}^\infty f(y)\int_0^\infty e^{-\alpha t} \frac{e^{-\frac{|x-y|^2}{2t}}}{\sqrt{2\pi t}} \, dt \, dy \\ \end{align}
So one can guess that
$$ \int_0^\infty e^{-\alpha t} \frac{e^{-\frac{|x-y|^2}{2t}}}{\sqrt{2\pi t}} \, dt = \frac{\exp{- \sqrt{2\alpha}|x-y|}}{\sqrt{2 \alpha}}$$
But How do we get such result?
Here is my attempt at proving this:
$$\int_0^\infty e^{-\alpha t} \frac{e^{-\frac{|x-y|^2}{2t}}}{\sqrt{2\pi t}} \, dt = \int_0^\infty \exp\bigg\{-\alpha t - \frac{|x-y|^2}{2t}\bigg\} \frac{1}{\sqrt{2\pi t}} \, dt \\ = \int_0^\infty \exp\bigg\{-\bigg(\sqrt{\alpha t} - \frac{|x-y|}{\sqrt{2t}}\bigg)^2 + \sqrt{2\alpha}|x-y|\bigg\} \frac{1}{\sqrt{2\pi t}} \, dt\\ = \exp\bigg\{- \sqrt{2\alpha}|x-y|\bigg \}\int_0^\infty \exp\bigg\{-\bigg(\sqrt{\alpha t} - \frac{|x-y|}{\sqrt{2t}}\bigg)^2 \bigg\} \frac{1}{\sqrt{2\pi t}} \, dt$$
So now I am at the point that I need to prove
$$\int_0^\infty \exp\bigg\{-\bigg(\sqrt{\alpha t} - \frac{|x-y|}{\sqrt{2t}}\bigg)^2 \bigg\} \frac{1}{\sqrt{2\pi t}} \, dt = \frac{1}{\sqrt{2 \alpha}} $$
But I don't see how to do it
If $A,B>0$ $$\int_{0}^{+\infty}\exp\left(-A^2 x^2-\frac{B^2}{x^2}\right)\,dx =\sqrt{\frac{B}{A}}\int_{0}^{+\infty}\exp\left(-ABx^2-\frac{AB}{x^2}\right)\,dx\tag{1}$$ hence it is enough to compute: $$ I(C)=\int_{0}^{+\infty}\exp\left(-C^2 x^2-\frac{C^2}{x^2}\right)\,dx. \tag{2}$$ By splitting the integration range as $(0,1)\cup(1,+\infty)$ and setting $t+\frac{1}{t}=u$ we have: $$ I(C) = \int_{2}^{+\infty}\frac{t}{\sqrt{t^2-4}} e^{-C^2(t^2-2)}\,dt = \int_{0}^{+\infty}\frac{1}{\sqrt{u}}e^{-C^2(4u+2)}\,du\tag{3}$$ then setting $u=v^2$ we get: $$ I(C) = 2e^{-2C^2}\int_{0}^{+\infty}e^{-4C^2 v^2}\,dv = \color{red}{\frac{\sqrt{\pi}}{2C}\cdot e^{-2C^2}}.\tag{4}$$