I’m interested in the form of the photon propagator in position space when expressed in a general $R_{\xi}$ gauge. The integral representation of this propagator is usually written as the sum of two Fourier transforms:
$$G_F^{\mu{}\nu{}}\left(x-y,\xi{}\right)=-i{\ g}^{\mu{}\nu{}}\int\frac{d^4k}{{\left(2\pi{}\right)}^4}\frac{e^{-ik\bullet{}\left(x-y\right)}}{k^2+i\epsilon{}}+i\ \xi{}\int\frac{d^4k}{{\left(2\pi{}\right)}^4}\frac{e^{-ik\bullet{}\left(x-y\right)}k^{\mu{}}k^{\nu{}}}{{(k^2+i\epsilon{})}^2}.\tag{1}$$
The value of the gauge parameter $\xi$ is arbitrary in general, but often is selected to impose a specific gauge condition on the propagator (e.g., $\xi$=0 for Feynman or $\xi$=1 for Landau). However, I want to leave $\xi$ unspecified.
Now the first integral in (1) just gives the massless scalar Feynman propagator, which has a well-known form in position space, so the term in (1) proportional to the metric is explicitly:
$$i\ g^{\mu{}\nu{}}G_F(x-y)=g^{\mu{}\nu{}}\left[\frac{1}{{4\pi{}}^2}PV\frac{1}{{(x-y)}^2}+\frac{i}{4\pi{}}\delta{}\left({(x-y)}^2\right)\right].\tag{2}$$
But what is the explicit form of the second photon propagator term proportional to $\xi$? The second integral involves higher-order poles and a more complicated tensor structure. Can anyone point me to a reference which demonstrates how to evaluate this second integral and then displays the result in position space?