Direct nonlinear Fourier transform algorithms for the computation of solitonic spectra in focusing nonlinear Schrödinger equation

Vasylchenkova, A., Prilepsky, J.e., Shepelsky, D. and Chattopadhyay, A. (2019). Direct nonlinear Fourier transform algorithms for the computation of solitonic spectra in focusing nonlinear Schrödinger equation. Communications in Nonlinear Science and Numerical Simulation, 68 , pp. 347-371.


Starting from a comparison of some established numerical algorithms for the computation of the eigenvalues (discrete or solitonic spectrum) of the non-Hermitian version of the Zakharov–Shabat spectral problem, this article delivers new algorithms that combine the best features of the existing ones and thereby allays their relative weaknesses. Our algorithm is modelled within the remit of the so-called direct nonlinear Fourier transform (NFT) associated with the focusing nonlinear Schrödinger equation. First, we present the data for the calibration of existing methods comparing the relative errors associated with the computation of the continuous NF spectrum. Then each method is paired with different numerical algorithms for finding zeros of a complex-valued function to obtain the eigenvalues. Next we describe a new class of methods based on the contour integrals evaluation for the efficient search of eigenvalues. After that we introduce a new hybrid method, one of our main results: the method combines the advances of contour integral approach and makes use of the iterative algorithms at its second stage for the refined eigenvalues search. The veracity of our new hybrid algorithm is established by estimating the convergence speed and accuracy across three independent test profiles. Along with the development of a new approach for the computation of the eigenvalues, our study also addresses the problem of computation of the so-called norming constants associated with the eigenvalues. We show that our formalism effectively amounts to accurate and fast enough computation of residues of the reflection coefficient in the upper complex half-plane of the spectral parameter.

Publication DOI:
Dataset DOI:
Divisions: Engineering & Applied Sciences > Electrical, electronic & power engineering
Engineering & Applied Sciences > Systems analytics research institute (SARI)
Engineering & Applied Sciences > Institute of Photonics
Engineering & Applied Sciences > Mathematics
Additional Information: © 2018, Elsevier. Licensed under the Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 International Funding: JEP acknowledges the support from the UK EPSRC Programme Grant UNLOC EP/J017582/1. JEP and DS are thankful to the Erasmus+ ERC mobility programme between the Aston University and Kharkiv National University that helped us to launch the collaborative activity. AC acknowledges the RISE-FRAMED grant.
Uncontrolled Keywords: Inverse scattering method,Nonlinear Schrödinger equation,Numerical algorithms,Signal processing,Numerical Analysis,Modelling and Simulation,Applied Mathematics
Full Text Link: https://arxiv.o ... /1708.01144.pdf
Related URLs: https://linking ... 007570418302855 (Publisher URL)
http://www.scop ... tnerID=8YFLogxK (Scopus URL)
Published Online Date: 2019-03
Published Date: 2019-03-01
Authors: Vasylchenkova, A.
Prilepsky, J.e.
Shepelsky, D.
Chattopadhyay, A. ( 0000-0001-5499-6008)

Export / Share Citation


Additional statistics for this record