Abstract

We obtain analytical conditions for the occurrence of numerical instability (NI) of a split-step method when the simulated solution of the nonlinear Schrödinger equation is close to a plane wave with nonzero carrier frequency. We also numerically study such an instability when the solution is a sequence of pulses rather than a plane wave. The plane-wave-based analysis gives reasonable predictions for the frequencies of the numerically unstable Fourier modes but overestimates the instability growth rate. The latter is found to be strongly influenced by the randomness of the signal’s profile: The more randomly it varies during the propagation, the weaker is the NI. Using an example of a realistic transmission system, we demonstrate that our single-channel results can be used to predict occurrences of NI in multichannel simulations. We also give an estimate for the integration step size for which NI, while present, will not affect simulation results for such systems. Using that estimate may lead to a significant saving of computational time.

© 2013 Optical Society of America

Solitons in regular and random split-step systems

Rodislav Driben, Boris A. Malomed, and Pak L. Chu
J. Opt. Soc. Am. B 20(11) 2338-2348 (2003)

Comparison of numerical methods for modeling laser mode locking with saturable gain

Shaokang Wang, Andrew Docherty, Brian S. Marks, and Curtis R. Menyuk
J. Opt. Soc. Am. B 30(11) 3064-3074 (2013)

Numerical modeling of (D+1)-dimensional solitons in a sign-alternating nonlinear medium with an adaptive fast Hankel split-step method

Georges Nehmetallah and Partha P. Banerjee
J. Opt. Soc. Am. B 22(10) 2200-2207 (2005)

Robust solitary waves in fiber transmission lines with semiconductor optical amplifiers

M. Gölles, I. M. Uzunov, and F. Lederer
J. Opt. Soc. Am. B 16(5) 689-694 (1999)

Modulational instability due to cross-phase modulation versus multiple four-wave mixing: the normal dispersion regime

Andrea Armaroli and Stefano Trillo
J. Opt. Soc. Am. B 31(3) 551-558 (2014)