I find myself confused as to whether relativistic propagators can be understood as probability amplitudes or not. This Wikipedia article for example, under 'Relativistic propagators', states that:
In relativistic quantum mechanics and quantum field theory the propagators are Lorentz-invariant. They give the amplitude for a particle to travel between two spacetime events.
Peskin and Schroeder similarly write on p.27 above eq. (2.50) that
$D(x-y)=\langle0|\phi(x)\phi(y)|0\rangle$ is the amplitude for a particle to propagate from $y$ to $x$.
On the other hand the answer to this question claims essentially that the propagator in relativistic QFT does not produce a probability because it is not conserved.
The problem is that in order to make sense as a probability, we should have something like (where $x=(t,\vec x)$) $$\int d^3 \vec x |D(x-y)|^2 = constant,$$
i.e. the total probability should not change over time (otherwise the particle might vanish), but this does not seem to be satisfied by the propagator (it is not even clear if this expression is Lorentz-invariant).
For example looking at the explicit expression for the Feynman propagator for timelike intervals ($\tau^2_{xy} = (x-y)^2 > 0$) $$ D(x-y) = -\frac{1}{4 \pi} \delta(\tau_{xy}^2) + \frac{m}{8 \pi \tau_{xy}} H_1^{(1)}(m \tau_{xy}) $$
trying to square this will yield terms having $\delta$-function squared, which I can't imagine can be properly integrated over.
So which one is correct? is the argument about probability conservation wrong and if so why, or is it actually satisfied but just hidden somewhere in the math, or is the statement in Wikipedia wrong?