Bound states in relativistic quantum field theory
Just as in QM, a bound state in QFT is a pole of the nonperturbative T-matrix (, i.e., the S-matrix minus 1 stripped from a delta function), and its possible wave functionals are the functionals in the range of the linear operator obtained as the residual of the T-matrix at the pole. (A proof is given in my answer to this question at Physics Overflow.
Perturbative QED, even in the fully rigorous treatment in
I haven't seen a single rigorous treatment of such an issue in quantum field theory.
Weinberg states in his QFT book (Vol. I) repeatedly that bound state problems (and this includes the Lamb shift) are still very poorly understood (though the Lamb shift is one of the most accurately predicted physical quantity). On p.564 he says,
There are, of course, methods for approximating bound state problems, based on Bethe-Salpeter equations, Schwinger-Dyson equations, and some other approaches. See, e.g., the review
Approximate Bethe-Salpeter equations for the bound state dynamics can be obtained by resumming infinite families of Feynman diagrams. See, e.g., Chapter 14 of Weinberg's QFT I.
But all of this is done in completely uncontolled approximations, and to get numerically consistent results is currently more an art than a science.
This leaves plenty of scope for interesting (but hard) new work on bound states on both the physical and mathematical side.
Arnold Neumaier (Arnold.Neumaier@univie.ac.at) A theoretical physics FAQ