TY - JOUR

T1 - Fast sixth-order algorithm based on the generalized Cayley transform for the Zakharov-Shabat system associated with nonlinear Schrodinger equation

AU - Medvedev, Sergey

AU - Chekhovskoy, Igor

AU - Vaseva, Irina

AU - Fedoruk, Mikhail

N1 - Funding Information:
S. Medvedev and I. Vaseva were supported by the Russian Science Foundation grant number 17-72-30006 . I. Chekhovskoy and M. Fedoruk were supported by the Russian Science Foundation grant number 20-11-20040 .
Publisher Copyright:
© 2021 Elsevier Inc.

PY - 2022/1/1

Y1 - 2022/1/1

N2 - The nonlinear Schrödinger equation (NLSE) is widely used in telecommunication applications, since it allows one to describe the propagation of pulses in an optical fiber. Recently some new approaches based on the nonlinear Fourier transform (NFT) have been actively explored to compensate for fiber nonlinearity and to exceed the limitations of nonlinearity-imposed limits of linear transmission methods. Despite the fact that the numerical solution of NLSE is a general problem, nevertheless, the optical community has been focusing on this issue. Improving the accuracy of the NFT algorithms remains an urgent problem in optics. In particular, it is important to increase the approximation order of the methods, especially in problems where it is necessary to analyze the structure of complex waveforms. To correctly describe them and their spectral parameters, more accurate and fast numerical methods are needed. We propose a novel general approach for constructing sixth-order (with respect to an integration step) finite-difference schemes for first-order linear differential systems. These schemes are based on the generalized Cayley transform and include exponential integrators as a special case. If the system has a time-dependent skew-hermitian matrix then the schemes conserve the quadratic first integral automatically. Then we apply our method to solve the direct spectral problem for the Zakharov-Shabat system. New schemes with fractional rational transition matrix allow the use of fast algorithms to solve the initial problem for a large number of values of the spectral parameter.

AB - The nonlinear Schrödinger equation (NLSE) is widely used in telecommunication applications, since it allows one to describe the propagation of pulses in an optical fiber. Recently some new approaches based on the nonlinear Fourier transform (NFT) have been actively explored to compensate for fiber nonlinearity and to exceed the limitations of nonlinearity-imposed limits of linear transmission methods. Despite the fact that the numerical solution of NLSE is a general problem, nevertheless, the optical community has been focusing on this issue. Improving the accuracy of the NFT algorithms remains an urgent problem in optics. In particular, it is important to increase the approximation order of the methods, especially in problems where it is necessary to analyze the structure of complex waveforms. To correctly describe them and their spectral parameters, more accurate and fast numerical methods are needed. We propose a novel general approach for constructing sixth-order (with respect to an integration step) finite-difference schemes for first-order linear differential systems. These schemes are based on the generalized Cayley transform and include exponential integrators as a special case. If the system has a time-dependent skew-hermitian matrix then the schemes conserve the quadratic first integral automatically. Then we apply our method to solve the direct spectral problem for the Zakharov-Shabat system. New schemes with fractional rational transition matrix allow the use of fast algorithms to solve the initial problem for a large number of values of the spectral parameter.

KW - Direct scattering transform

KW - Fast numerical methods

KW - Nonlinear Fourier transform

KW - Nonlinear Schrödinger equation

KW - Zakharov-Shabat problem

UR - http://www.scopus.com/inward/record.url?scp=85117019406&partnerID=8YFLogxK

U2 - 10.1016/j.jcp.2021.110764

DO - 10.1016/j.jcp.2021.110764

M3 - Article

AN - SCOPUS:85117019406

VL - 448

JO - Journal of Computational Physics

JF - Journal of Computational Physics

SN - 0021-9991

M1 - 110764

ER -