Interaction of nonlinear Schrodinger solitons with an external potential

Employing a particularly suitable higher-order symplectic integration algorithm, we integrate the one-dimensional nonlinear Schrodinger equation numerically for solitons moving in external potentials. In particular, we study the scattering off an interface separating two regions of constant potential. We find that the soliton can break up into two solitons, eventually accompanied by radiation of non-solitary waves. Reflection coefficients and inelasticities are computed as functions of the height of the potential step and of its steepness.