Abstract
We discuss the scattering theory pertaining to the coupled πNNNN system and approach the problem in two independent ways. The first one starts from a Hamiltonian formalism and coupled Schrödinger equations, whereas the second one employs an off-mass-shell relativistic theory of classifying perturbation diagrams. Both ways lead to connected equations among transition operators in which πNN vertices, as well as nucleon propagators, are completely dressed and renormalized. Furthermore, the physical amplitudes obey two- and three-body unitarity relations. The resultant equations form a sound theoretical basis for subsequent numerical calculations leading to the evaluation of physical observables in the reactions π+dπ+d, π+dN+N, and N+N+N+N.