Integrating factor techniques applied to the Schr{ö}dinger-like equations. Comparison with Split-Step methods

The nonlinear Schr{\"o}dinger and the Schr{\"o}dinger-Newton equations model many phenomena in various fields. Here, we perform an extensive numerical comparison between splitting methods (often employed to numerically solve these equations) and the integrating factor technique, also called Lawson method. Indeed, the latter is known to perform very well for the nonlinear Schr{\"o}dinger equation, but has not been thoroughly investigated for the Schr{\"o}dinger-Newton equation. Comparisons are made in one and two spatial dimensions, exploring different boundary conditions and parameters values. We show that for the short range potential of the nonlinear Schr{\"o}dinger equation, the integrating factor technique performs better than splitting algorithms, while, for the long range potential of the Schr{\"o}dinger-Newton equation, it depends on the particular system considered.