The bilateral Laplace transform of a Gaussian function could be established as: $$e^{x^2/2}=\frac{1}{\sqrt{2 \pi}} \int_{-\infty}^{\infty}e^{-xy}e^{-y^2/2} dy$$
Then what should be a similar relation for a Gaussian function with complex variable:$e^{z^2/2}$?
Let $z$ be a complex number. The integral is convergent for all $z \in \mathbb{C}$. Completing the square we get: $$ \int_{-\infty}^\infty \mathrm{e}^{-z \cdot y} \mathrm{e}^{-y^2/2} \mathrm{d} y = \mathrm{e}^{z^2/2}\int_{-\infty}^\infty \mathrm{e}^{-(y+z)^2/2} \mathrm{d} y $$ The latter integral multiplying the $\mathrm{e}^{z^2/2}$ is a constant, independent of $z$. Evaluating it at $z=0$ establishes the result.