Well, I guess here we have a typical "clash of civilizations" between mathematicians and physicists :-)).
Strictly speaking you are right: So far there is no example for a real-world relativistic QFT of interacting particles that could be strictly formulated in terms of mathematical rigor. Nevertheless relativistic QFT, as used in its perturbative setup to calculate S-matrix elements (cross sections) is a pretty successful theory.
Also I don't see, what divergent series have to do with the original question. As was already pointed out in this thread the definition of the operator is best established in momentum representation
\tilde{\psi}(\vec{p}) \mapsto \sqrt{1+\vec{p}^2} \tilde{\psi}(\vec{p}).
There is no need to expand this in a power series of \vec{p}^2.
The issue in relativistic field theory is that causality has a more restrictive meaning than in non-relativistic field theory: In the former case there must not be "faster-than-light propagation", i.e., a signal must always propagate with a speed less than the speed of light, i.e., the time evolution of a wave packet with compact support in space at an initial time must always give a wave packet with compact spatial support at any later time.
At the same time, the energy should be bounded from below. From that idea, one comes to the idea to investigate the non-local wave equation
\mathrm{i} \partial_t \phi(t,\vec{x}) = \sqrt{1-\vec{\nabla}^2} \phi(t,\vec{x}),
interpreting the operator on the right-hand side as the Hamilton operator for a non-interacting spin-less particle of mass m=1.
As emphasized above the operator is better defined in momentum space, and in momentum space the equation reads
\mathrm{i} \partial_t \tilde{\phi}(t,\vec{p})=\sqrt{1+\vec{p}^2} \tilde{\phi}(t,\vec{p}).
This implies
\tilde{\phi}(t,\vec{p})=\exp(-\mathrm{i} \omega_p t) \tilde{\phi}_0(\vec{p}).
In spatial representation this means
\phi(t,\vec{x})=\int_{\mathbb{R}^3} \mathrm{d}^3 \vec{x}' U(t,\vec{x}-\vec{x}') \phi_0(\vec{x}')
with
U(t,\vec{y})=\int_{\mathbb{R}^3} \frac{\mathrm{d}^3 \vec{p}}{(2 \pi)^3} \exp(-\mathrm{i} \sqrt{1+\vec{p}^2}) \exp(\mathrm{i} \vec{p} \cdot \vec{y}).
It can be easily shown that this does not vanish for space-like arguments. For \vec{y}^2 \gg t^2 we can estimate the integral with the stationary-phase method giving
U(t,\vec{y}) \propto \exp[-\sqrt{\vec{y}^2-t^2}],
which is small but nowhere 0.
This means that parts of the wave propagate faster than the speed of light, and this violates causality in a relativistic framework. The only known solution within local relativistic QFT is to introduce creation operators in addition to annihilation in the mode decomposition of the Klein-Gordon field, i.e., to introduce anti-particles into the game to fulfill the boundedness of energy and causality at the same time. See Weinberg, The Quantum Theory of Fields, Vol. I for further details.