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

Abstract

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

accepted to Communication of Nonlinear Science and Numerical Simulations

References in corpus (1)

• Control and Detection of Discrete Spectral Amplitudes in Nonlinear Fourier Spectrum

Cited by in corpus (9)

• Fast Nonlinear Fourier Transform Algorithms Using Higher Order Exponential Integrators
• Exponential Fourth Order Schemes for Direct Zakharov-Shabat problem
• Direct scattering transform: catch soliton if you can
• Introducing phase jump tracking -- a fast method for eigenvalue evaluation of the direct Zakharov-Shabat problem
• Novel Numerical Algorithm with Fourth-Order Accuracy for the Direct Zakharov-Shabat Problem
• Bound-state soliton gas as a limit of adiabatically growing integrable turbulence
• Successive Eigenvalue Removal for Multi-Soliton Spectral Amplitude Estimation
• Fast sixth-order algorithm based on the generalized Cayley transform for the Zakharov-Shabat system in optical applications
• An Efficient Nonlinear Fourier Transform Algorithm for Detection of Eigenvalues from Continuous Spectrum