vanhees71 said:
Well, but that's the problem. Wave functions do not make too much sense for relativistic QT, at least not the same sense as they make in non-relativistic QM, and that's the problem here. All you can say from relativistic QFT are detection probabilities of a particle (i.e., an asymptotic free state, where a particle interpretation is possible) with a detector at a given position. This information is formally encoded in the field-operator autorcorrelation functions (aka ##N##-point Green's functions).
In post 56 I wrote a single particle state
$$|\Psi \rangle = \int d^3 {\bf x}\, \psi (t,{\bf x}) a^\dagger ({\bf x})|0 \rangle, $$
where ##\psi (t,{\bf x}) = \langle {\bf x}|\Psi (t,{\bf x}) \rangle ## is a complex valued wave packet profile, that is a wave function in the x-representation.
It satisfies the relativistic Schroedinger equation (known as the Foldy equation):
$$i {\dot \psi} = \omega_{\bf x} \psi (t,{\bf x}),$$
where ##\omega_{\bf x} = \sqrt{m^2 - \nabla^2}##. Expanding
$$\sqrt{m^2 - \nabla^2} = m -\frac{\nabla^2}{2 m} + ... $$
and taking the ansatz
$$\psi (t,{\bf x}) = {\rm e}^{i m t} {\tilde \psi} ({\bf x})$$
we obtain the nonrelativistic Schroedinger equation as an approximation:
$$i {\dot {\tilde \psi}} = -\frac{\nabla^2}{2 m} {\tilde \psi}.$$
On the other hand, by inserting ##\psi = \psi_R + i \psi_I## into the equation ##i {\dot \psi} = \omega_{\bf x} \psi (t,{\bf x})##, we obtain the following system of first order differential equation for the real and imaginary components of ##\psi:##
$${\dot \psi} _R= \omega_{\bf x} \psi_I ,~~~~~~~~
{\dot \psi} _I= -\omega_{\bf x} \psi_R .$$
From the first equation we have
$$\psi_I = \omega_{\bf x}^{-1} {\dot \psi}_R .$$
Insertinf this into the second equation, we obtain
$$ \omega_{\bf x}^{-1}{\ddot \psi}_R = -\omega_{\bf x} \psi_R .$$
Multiplying the latter equation from the left by ##\omega_{\bf x} ##,
we obtain
$${\ddot \psi_R} + \omega_{\bf x}^2 \psi_R = 0,$$
i.e.,
$${\ddot \psi_R} + (m^2 - \nabla^2 )\psi_R = 0 ,$$
which is the Klein-Gordon equation for the real field ##\phi \equiv \psi_R##. So we have found that a real Kleing-Gordon (c-number ) field takes place within the framework of the relativistic quantum field theory of a Hermitian scalar operator-valued field ##\varphi##.
In other words, we have quantized a real scalar field by promoting it to the operator-valued field, constructed with them the creation operators, and position states, then taken a superposition of such position single particle state, where the superposition "coefficients"
were ##\psi (t,{\bf x})##, satisfying the first order Foldy equation, which can be rewritten as the second order Klein-Gordon equation.
Thus, if we have a real Kleing-Gordon field (c-number, not operator) ##\phi##, the role of the wave function has the superposition
$$\psi = \phi + i \omega_I^{-1} {\dot \phi},$$
of the field ##\phi## and its conjugated momentum ##{\dot \phi}##.