I'm trying to understand the use of the Euclidean correlation functions in QFT. I chased down the problems I was having to how they manifest in the simplest example I could think of: the two-point propagator for the Klein-Gordon equation. V. P. Nair (pdf pages 57-58) starts with the Feynman propagator for the Klein Gordon equation,
$$ G(x,y) = \lim_{\epsilon\to0^+}\int_{-\infty}^\infty dk_0 \int_{\mathbb{R}^3}d^3\textbf{k}\; \frac{i}{k_0^2-\textbf{k}^2-m^2+i\epsilon}e^{-ik_0(x_0-y_0)+i\textbf{k}\cdot(\textbf{x}-\textbf{y})}.\tag{4.13} $$
He then argues that you can deform the contour such that the $k_0$ integral goes up the imaginary axis, to get
$$ G(x,y) = \int_{-i\infty}^{i\infty} dk_0 \int_{\mathbb{R}^3}d^3\textbf{k}\; \frac{i}{k_0^2-\textbf{k}^2-m^2}e^{-ik_0(x_0-y_0)+i\textbf{k}\cdot(\textbf{x}-\textbf{y})},\tag{4.17} $$
at which point you're a change of variables away from getting the relationship we want between the Minkowski and Euclidean propagators. Nair says that "there is no crossing of poles of the integrand in this deformation", and I can see that: you're deforming the contour through the upper right and lower left quadrants of the complex plane, so avoid the poles. My issue is what about the quarter-circular contours at infinity? You have to leave the endpoints fixed when you deform the contour, so to get the $k_0$ integral to go along the imaginary line we must have a contour that joins the ends of the imaginary to the real line which vanishes. But surely this can't be the case in both the upper right and lower left contours, as the integrand has a factor $\propto \exp\left(\text{Im}\{k_0\} x_0\right)$, which depending on the sign of $x_0$ will diverge at either large positive imaginary $k_0$ or large negative imaginary $k_0$?
There is a slightly different way of driving at the same problem. Nair arrives at the relation
$$ G(x,y) = G_E(x,y)|_{x^4=-ix^0,y^4=-iy^0},\tag{4.18} $$
where the Euclidean propagator is defined
$$ G_E(x,y) = \int_{\mathbb{R}^4} d^4k\; \frac{1}{\sum_{j=1}^4(k_j)^2+m^2}e^{i\sum_{j=1}^4k_j(x_j-y_j)}.\tag{4.19} $$
The issue here is that if you put imaginary values of $x_4-y_4$ into the defining integral then you get an exponential divergence in the $k_4$ integral, so the result is poorly defined.
So what's going on here? Am I missing something obvious or is Nair doing some egregious handwaving? And, if the latter, could you possibly point me in the direction of a treatment of the relationship between the Euclidean and Minkowski correlation functions which isn't quite as mathematically technical as the Osterwalder and Schrader paper? (Which is all I've managed to find referenced elsewhere!) When I've tried to find the relation in more complicated and general cases - for instance by looking at the partition function expressed as a path integral - I think I've stumbled on more or less the same problem, of this divergence of the exponential factor, so I do feel that if I get this derivation of the KG propagator sorted then the rest ought to fall into place.