I. INTRODUCTION can be written in a way similar to that of the paraxial Helmholtz equation that describes the evolution of the optical beams in a resonator cavity. The time-dependent solution of the 1-D Schrödinger equation11. W. Louisell, Quantum Statistical Properties of Radiation (Chichister-Wiley, 1973). We will pursue the method proposed in Ref. 22. O. Morandi, “Quantum motion with trajectories: beyond the Gaussian beam approximation,” J. Phys. A: Math. Theor.