 Solution

The Heisenberg picture evolution of the annihilation operator follows from $[H,a_{\mathbf{p}}]=-E_{\mathbf{p}}{a}_{\mathbf{p}}$. With $p^0=E_{\mathbf{p}}$, this gives the mode expansion of a free fieldDefine ${\Delta }^{+}(z)={\int }_{\mathbf{p}}{e}^{-ip\cdot z}/(2E_{\mathbf{p}})$. The time-ordered product then gives the Feynman propagatorThe last expression follows by closing the $p^0$ contour integral below for $z^0>0$ and above for $z^0<0$; the scalar Feynman propagator pole prescription supplies the required positive- and negative-frequency terms.

For $z^0\ne 0$ every term obeys the homogeneous Klein-Gordon equation. At $z^0=0$, the equal-time commutator of the fields vanishes, so the propagator is continuous and there is no derivative of a Dirac delta function. Its first time derivative has the jumpThe second time derivative therefore supplies precisely the source:Thus the defined Feynman propagator is a Green's function with source normalization $-i$; $i{\Delta }_F$ has source normalization $+1$. Keeping this factor avoids confusing the vacuum expectation value convention with conventions which define the propagator after dividing by $i$.^[ref]