On the efficiency of different numerical methods for the calculation of intrapulse Raman scattering of optical solitons

Abstract

In this paper we compare the performance of different numerical methods for the calculation of the asymptotic evolution of soliton self-frequency shift in the presence of intrapulse Raman scattering (IRS) in optical fibers. First we have calculated the order of global accuracy for the fundamental soliton and the second-order bound state of the unperturbed nonlinear Schrodinger equation for the following numerical methods: the simple split step (SS) method, the full SS method, the reduced SS method, the reduced SS method with fourth order Runge–Kutta (RK4), the Blow–Wood method, the fourth order Runge–Kutta in the interaction picture (RK4IP) method and the Agrawal SS method with one and two iterations. We have shown that the asymptotic evolution of soliton self-frequency shift in the presence of IRS in optical fiber can be best described by the Agrawal SS method (compared to the Blow–Wood method and the RK4IP method). The obtained numerical results for the soliton position and frequency are in agreement with the predictions for these parameters according to the perturbation theory. We have shown that in the presence of IRS the fundamental soliton quickly develops an oscillating tail on its left. The generated tail hardly influences the soliton propagation at large distances. At large distances the form of the soliton gradually becomes asymmetric. The increase of the frequency resolution leads to a notable increase of the maximum propagation distance of the numerical soliton under the influence of IRS.