This is Problem (6.1) from Schwartz's QFT and the Standard Model. I am trying to directly calculate, by performing the integral over momenta, the explicit position-space expression of the Feynman propagator. That is, performing the integral
\begin{equation} D_{F} = \int \, \frac{d^4 k}{(2 \pi)^4} \, \frac{i}{k^2 - m^2 + i \varepsilon} e^{i k (x_{1} - x_{2})}. \end{equation} I know that one can quote the solution that this equals some function including a modified Bessel function of the second kind as
\begin{equation} D_{F} = \frac{m}{4 \pi^2 r + i \varepsilon} K_{1} (mr), \end{equation} for spacelike separations, where $r = \vec{x}_{1} - \vec{x}_{2}$.
You are asked, after working this expression out, to take the massless limit, which evaluates to
\begin{equation} \lim_{m \rightarrow 0} D_{F} = \frac{1}{4 \pi^2 r^2 + i \varepsilon}. \end{equation} I understand the logic in these posts (The poles of Feynman propagator in position space, How to obtain the explicit form of Green's function of the Klein-Gordon equation?) in how to manipulate the integral to get it into an integral representation of $K_{1}$, but Question 1 would be how is it that this limit arrived at? I have looked at some of the different representation of $K_{1}$ and am still unsure how this limit is calculated. My shaky first thoughts is that it approaches $\frac{1}{mr}$, and that you just say $mr(i \varepsilon) = i \varepsilon$ as this $\varepsilon$ is only used for keeping track of contours, etc., and will be taken to zero. This is what seems to be occurring in the top answer in this post: Propagator of a scalar in position space. $K_{1}$ does seem to visually appear like a $1/r$ type function when plotted.
My other question (I will be happy to move this to separate post if inappropriate) involves an apparently different derivation of the above, some of the steps of which I would appreciate being explained, or hinted at. For completeness, I want to derive the dependence on $K_{1}$, not just state it and take the limit:
Question 2: Specifically, I don't quite understand how the integral over $d k^{0}$ is being performed from the second the third line (have we transformed to a frame in which $\vec{k} = 0$?). The question mentions use of Schwinger parameters, which I think is perhaps treated in the linked posts above, but I'm not sure as they don't seem to be in the form stated in Schwartz.
I'm a bit unsure on the last couple of steps but I assume that I'll find some integral trig representation of $K_{1}$ somewhere that'll match up.
Question 3: Why can we just invoke Lorentz symmetry here to reinstate the rest of the coordinate dependence? Is this just because we transformed frames to ease calculation?