Numerical modeling considerations for an applied nonlinear Schrödinger equation.

A model for nonlinear optical propagation is cast into a split-step numerical framework via a variable stencil-size Crank-Nicolson finite-difference method for the linear step and a choice of two different nonlinear integration schemes for the nonlinear step. The model includes Kerr, Raman scattering, and ionization effects (as well as linear and nonlinear shock, diffraction, and dispersion). We demonstrate the practical importance of numerical effects when interpreting computational studies of high-intensity optical pulse propagation in physical materials.