Multiprecision computation of bright and dark solitons in the discrete nonlinear Schr\"odinger equation
In the authors' words
We study the spectral stability of bright and dark solitons in the discrete nonlinear Schr\"odinger (DNLS) equation using multiprecision arithmetic. The eigenvalues governing stability are exponentially small in the lattice spacing and cannot be resolved with standard double precision. To address this, we develop a computational framework combining multiprecision arithmetic, an exact Jacobian for the stationary problem, and a squared-operator formulation for spectral analysis. This enables accurate resolution of exponentially small eigenvalues and direct comparison with exponential-asymptotic predictions. Our results show that onsite bright solitons are spectrally stable, whereas intersite bright solitons and both onsite and intersite dark solitons are unstable. Bright solitons require only a few eigenvalues and allow efficient large-scale computations, while dark solitons demand higher precision due to their proximity to the continuous spectrum. Simulations up to \(N=65{,}250\) grid points (31.7 GB RAM) highlight the necessity of multiprecision arithmetic for capturing beyond-all-orders spectral effects.
Appeared: Thursday, September 24. arXiv. Preprint with a published version.
DOI: 10.1016/j.physleta.2026.132170
Published version: 10.1016/j.physleta.2026.132170
Authors' comment: Published