Long-time asymptotic for the derivative nonlinear Schr"odinger equation with decaying initial value

We present a new Riemann-Hilbert problem formalism for the initial value problem for the derivative nonlinear Schr"odinger (DNLS) equation on the line. We show that the solution of this initial value problem can be obtained from the solution of some associated Riemann-Hilbert problem. This new Riemann-Hilbert problem for the DNLS equation will lead us to use nonlinear steepest-descent/stationary phase method or Deift-Zhou method to derive the long-time asymptotic for the DNLS equation on the line.

💡 Research Summary

The paper develops a novel Riemann–Hilbert (RH) problem framework for the Cauchy problem of the derivative nonlinear Schrödinger (DNLS) equation on the real line. Starting from the Lax pair associated with DNLS, the authors construct the direct scattering problem, define Jost solutions, and introduce the scattering matrix whose entries consist of the reflection coefficient (r(k)), the transmission coefficient (a(k)), and a finite set of discrete eigenvalues ({k_j}) with associated norming constants. These data are then encoded into a matrix‑valued RH problem whose jump matrix on the real axis contains the oscillatory exponential factor (\exp(\pm 2i t\theta(k))), where the phase (\theta(k)=k^{2}+\alpha/k^{2}) reflects the derivative nonlinearity.

The central technical contribution is the application of the Deift–Zhou nonlinear steepest‑descent method to this RH problem. By locating the stationary points of (\theta(k)) (solutions of (\theta’(k)=0)) the authors design a contour deformation that opens lenses around the real axis, thereby converting the original jump matrix into a product of a near‑identity matrix and a model jump supported only near the stationary points. In the neighborhoods of these points, the RH problem is reduced to a solvable model expressed in terms of Airy functions (or, equivalently, parabolic cylinder functions), yielding explicit local parametrices. For each discrete eigenvalue a small circular contour is introduced, and a separate “soliton” RH problem is solved, providing the contribution of isolated solitary waves. Matching the global solution with the local parametrices produces an asymptotic expansion whose error is bounded by (O(t^{-1/2})).

The final long‑time asymptotics consist of two distinct parts. The continuous spectrum generates a radiative component that decays like (t^{-1/2}) and oscillates with phase (\theta(k_{0})t), where (k_{0}) is the stationary point. The discrete spectrum yields soliton terms that retain a fixed amplitude and acquire a phase shift determined by the norming constants; these solitons travel without dispersion. Interaction terms between solitons and radiation appear only at higher order and are suppressed by the steepest‑descent analysis.

Overall, the work provides a rigorous, systematic derivation of the long‑time behavior of DNLS solutions for rapidly decaying initial data, improving upon earlier inverse‑scattering approaches by delivering precise error estimates and a clear separation of soliton and radiation dynamics. The methodology is readily adaptable to other integrable equations with derivative nonlinearities and has potential implications for physical contexts such as nonlinear optics, plasma waves, and Bose–Einstein condensates where DNLS serves as a model.