Green Function expressed in terms of Hankel function (of the second kind) I am reading Davies' and Birrell's book on QFT in Curved Spacetimes. In Chapter 2 and specifically in the subchapter 2.7, the authors derive an expression for the Feynman propagator
$$G_F(x,x')=\frac{-i\pi}{(4\pi i)^{n/2}}\Big(\frac{2m^2}{-\sigma+i\epsilon}\Big)^{(n-2)/4}
H^{(2)}_{\frac{1}{2}n-1}\Big\{[2m^2(\sigma-i\epsilon)]^{1/2}\Big\}\tag{2.77}$$
by performing the contour integration and by choosing one of the poles to lie on the upper half of the $\text{Im}(k^0)$ - $\text{Re}(k^0)$ plane and the other one to lie on the lower half of the abovementioned plane (just like any similar derivation in QFT textbooks). This is the first time I see the Feynman propagator in this form and I do not know how to derive it, so any tip would be appreciated...
My attempt revolves around performing the $k^0$ integral first, from the expression
$$G_F(x,x')=\int\frac{d^nk}{(2\pi)^n}\frac{e^{-ik\cdot x}}{k^2-m^2}$$
(where $k^2=k\cdot k$ with the $\cdot$ denotes a four-vector multiplication) and then trying to bring the resulting expression to the form of the integral representation of the Hankel function, namely
$$H^{(2)}_{\nu}(x)=\frac{1}{\pi i}\int_0^{\infty e^{i\pi}}dt 
\frac{e^{x/2(t-1/t)}}{t^{\nu+1}}$$
So, I wanted to ask whether or not

*

*I am headed to the right direction and if yes, what would my $t$ variable be? The three momentum vector magnitude perhaps? (I know this sounds like a "do my derivation for me" question, but I wouldn't want that. Some tips would be nice... Specifically, there is a gap between the integral representation expression and the result I get from performing the $k^0$ contour integral and it would be nice if somehow someone would breach that gap)


*is there another approach for deriving the propagator without the integral representation of the Hankel function... For instance, I thought of writting the Kleing-Gordon equation for the propagator in spherical coordinates and then solve it. But then I would need boundary conditions and hence another dead end...
Any help/comments on my attempts would be appreciated.
 A: There is a useful class of dummy integration variable trick, which applies also later when calculating Feynman diagrams. You can make the following observation:
$$
\frac{1}{x} = \int_0^{+\infty}dt e^{-xt}
$$
Thanks to the definition of the Feynman propagator, which makes you avoid the poles in a particular way, using contour integration, you can modify the contour to integrate $k^0$ over the imaginary axis without crossing any poles. You get:
$$
G_F = \int\frac{d^nk}{(2\pi)^n}\frac{e^{-ikx}}{k^2-m^2+i\epsilon} \\
 = -i\int\frac{d^nk_E}{(2\pi)^n}\frac{e^{ik_Ex_E}}{k_E^2+m^2}
$$
with $k^0 = ik_E^0$ and $x^0 = ix_E^0$ and the spatial components unchanged, so that the new dot product in the exponential is the usual euclidean dot product. Now you can apply the trick:
$$
G_E=-i\int_0^{+\infty}dt\int\frac{d^nk_E}{(2\pi)^n}\exp(ik_Ex_E-t(k_E^2+m^2)) \\
=-i\int_0^{+\infty}dt\frac{1}{(2\pi)^n}\sqrt{\frac{\pi}{t}}^ne^{-m^2t-x_E^2/4t}
$$
where I completed the square and performed the gaussian integral (and swapped the integrals). By rescaling, you'll obtain the integral representation of a Bessel function. If you're keeping the same $x$, ie $x_E$ is imaginary, you'll get Hankel's function, but you'll need to be careful in the definition of the above integrals and modify the contour. The easiest way it to Wick rotate, and assume $x_E$ is real. You will get rigorously the  modified Bessel functions of the second kind. Extending $G_E$ thus calculated on the imaginary axis of $x^0$ will give you Hankel's function.
If you try to directly solve the KG equation, you need to implement the pole avoidance of the Feynman propagator somehow, which is not easy to do. You could try to do it in a similar fashion by calculating $G_E$ on the imaginary axis and recognising it as the Green's function of $-\Delta_E+m^2$. This is where you can use polar coordinates, but the calculation is tricky due to the complex exponential. It's doable in small dimensions, especially easy in 3D euclidean (ie 1+2D spacetime) which gives you the Yukawa potential, but in general dimension $n$, the calculation is way harder, only to get back to the same integral representation.
Hope this helps.
