It is well known that the Feynman propagator can be expressed as a Green function of the Klein-Gordon operator:
$$G_F(x'-x)=i\int\frac{d^4p}{(2\pi)^4}\frac{e^{-ip\cdot(x'-x)}}{p^2-m^2+i\epsilon}$$
On the other hand, we can define the retarded and advanced propagators such as the former exists only if $$(t'-t)>0$$ and the second exists only if $$(t'-t)<0$$. For example, the retarded propagator has the expression:
$$G_R(x'-x)=i\int\frac{d^4p}{(2\pi)^4}\frac{e^{-ip\cdot(x'-x)}}{(p^0-E_p+i\delta)(p^0+E_p+i\delta)}$$
Or performing the temporal integral:
$$G_R(x'-x)=\theta(x'^{0}-x^0)\int\frac{d^3p}{(2\pi)^3}\frac{e^{-i\vec{p}\cdot(\vec{x'}-\vec{x})}}{2E_p}(e^{iE_p(x'^{0}-x^0)}-e^{-iE_p(x'^{0}-x^0)})$$
It is well known that the Feynman propagator is a Green function of the Klein-Gordon operator. However, how can we verify that the retarded propagator is a Green function of the Klein-Gordon operator?
From what I've read, the Feynman and the retarded propagator differ by a solution for the Klein-Gordon equation, but all the books and articles I've seen skip this demonstration.