Summing Lax series for nonparaxial beam propagation

It is shown that the approach proposed by Lax et al. [Phys. Rev. A 11, 1365 (1975)] for studying the propagation of an electromagnetic beam beyond the paraxial approximation can be efficiently employed evenwhen the beam under consideration presents a very nonparaxial character. The method that we present consists of applying a nonlinear resummation scheme, the so-called delta transformation, to the divergent perturbative series arising from the Lax scheme. Numerical results pertinent to the evaluation of transverse and longitudinal components of the electric field are presented for the particular case of vectorial Gaussian beams, showing the effectiveness of the method.