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

A. Vasylchenkova, J.e. Prilepsky, D. Shepelsky, A. Chattopadhyay

Research output: Contribution to journalArticlepeer-review


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.

Original languageEnglish
Pages (from-to)347-371
Number of pages25
JournalCommunications in Nonlinear Science and Numerical Simulation
Early online date6 Sept 2018
Publication statusPublished - 1 Mar 2019

Bibliographical note

© 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.


  • Inverse scattering method
  • Nonlinear Schrödinger equation
  • Numerical algorithms
  • Signal processing


Dive into the research topics of 'Direct nonlinear Fourier transform algorithms for the computation of solitonic spectra in focusing nonlinear Schrödinger equation'. Together they form a unique fingerprint.

Cite this