Numerical inverse scattering transform for the derivative nonlinear Schrodinger equation

Shikun Cui,Zhen Wang

2023-08-21

Abstract:In this paper, we develop the numerical inverse scattering transform (NIST) for solving the derivative nonlinear Schrodinger (DNLS) equation. The key technique involves formulating a Riemann-Hilbert problem (RHP) that is associated with the initial value problem and solving it numerically. Before solving the RHP, two essential operations need to be carried out. Firstly, high-precision numerical calculations are performed on the scattering data. Secondly, the RHP is deformed using the Deift-Zhou nonlinear steepest descent method. The DNLS equation has a continuous spectrum consisting of the real and imaginary axes and features three saddle points, which introduces complexity not encountered in previous NIST approaches. In our numerical inverse scattering method, we divide the $(x,t)$-plane into three regions and propose specific deformations for each region. These strategies not only help reduce computational costs but also minimize errors in the calculations. Unlike traditional numerical methods, the NIST does not rely on time-stepping to compute the solution. Instead, it directly solves the associated Riemann-Hilbert problem. This unique characteristic of the NIST eliminates convergence issues typically encountered in other numerical approaches and proves to be more effective, especially for long-time simulations.

Numerical Analysis,Mathematical Physics

What problem does this paper attempt to address?

The problem that this paper attempts to solve is to develop the Numerical Inverse Scattering Transform (NIST) method to solve the Derivative Nonlinear Schrödinger (DNLS) equation. Specifically, the goals of the paper are:

1. **Construct the Riemann - Hilbert problem**: Transform the initial value problem into an associated Riemann - Hilbert problem (RHP) and solve it numerically. This is a crucial step in solving the DNLS equation.

2. **High - precision numerical calculation**: Before solving the Riemann - Hilbert problem, high - precision numerical calculations of the scattering data are required. This includes calculating discrete eigenvalues, reflection coefficients, and normalization constants.

3. **Deform the Riemann - Hilbert problem**: Since the DNLS equation has a continuous spectrum (the real and imaginary axes) and three saddle points, this brings complexity to the traditional NIST method. Therefore, the paper proposes a deformation strategy based on the Deift - Zhou nonlinear steepest - descent method, divides the \((x, t)\) plane into three regions, and proposes specific deformation methods for each region to reduce the computational cost and minimize the computational error.

4. **Numerical inverse scattering**: Utilize the deformed Riemann - Hilbert problem and solve this problem by the Chebyshev collocation method to recover the solution of the DNLS equation.

5. **Verify the effectiveness of the method**: Verify the effectiveness and accuracy of the NIST method in long - time simulations through numerical results and comparisons with other methods (such as the Fourier spectral method).

In summary, this paper aims to develop an efficient and accurate numerical method for solving the Derivative Nonlinear Schrödinger equation, especially for solving long - time evolution problems. This method not only avoids the time - stepping problem in traditional numerical methods but also solves the problems of convergence and error accumulation.