[论文解读] Stability of Solitary Waves for a Generalized Derivative Nonlinear Schrödinger Equation

Xiao Liu, Gideon Simpson+1 more|arXiv (Cornell University)|Jun 15, 2012

Advanced Mathematical Physics Problems参考文献 31被引用 4

一句话总结

本文研究了广义导数非线性薛定谔方程（gDNLS）中具有幂次非线性项 $\sigma > 0$ 的孤立波解的轨道稳定性。通过谱分析与数值模拟，结果表明稳定性关键取决于 $\sigma$ 与孤立子速度 $c$：当 $\sigma \geq 2$ 时孤立波不稳定，当 $\sigma < 1$ 时稳定，而在 $1 < \sigma < 2$ 范围内则表现出与速度相关的稳定性阈值，即慢速孤立子（$c < 2z_0\sqrt{\omega}$）稳定，快速孤立子不稳定。

ABSTRACT

We consider a derivative nonlinear Schrödinger equation with a general nonlinearity. This equation has a two parameter family of solitary wave solutions. We prove orbital stability/instability results that depend on the strength of the nonlinearity and, in some instances, their velocity. We illustrate these results with numerical simulations.

研究动机与目标

• 分析广义导数非线性薛定谔方程（gDNLS）中幂次非线性项 $\sigma > 0$ 的孤立波解的轨道稳定性。
• 确定非线性强度（$\sigma$）与孤立子参数（$\omega$，$c$）如何影响稳定性。
• 将先前关于标准 DNLS 方程（$\sigma = 1$）的研究结果扩展至更广泛的非线性类别。
• 为 $1 < \sigma < 2$ 范围内速度依赖的稳定性阈值提供数值证据。
• 基于谱分析与变分分析，建立孤立波轨道稳定或不稳定的条件。

提出的方法

• 研究考虑广义 DNLS 方程 $i\partial_t\psi + \partial_x^2\psi + i|\psi|^{2\sigma}\partial_x\psi = 0$，其存在双参数孤立波解。
• 孤立波解以形式 $\psi_{\omega,c}(x,t) = \phi_{\omega,c}(x - ct)e^{i\omega t}$ 构造，其中 $\phi_{\omega,c}$ 满足非线性常微分方程。
• 通过哈密顿量及守恒量（质量 $Q$，动量 $P$）的变分结构分析轨道稳定性，利用其对参数 $\omega$ 与 $c$ 的导数。
• 稳定性判据依赖于作用泛函二阶变分的符号，该符号由孤立波周围线性化算符的谱性质导出。
• 通过渐近展开与双曲函数相关的积分恒等式，推导出 $\partial_c Q$，$\partial_\omega Q$，$\partial_c P$ 与 $\partial_\omega P$ 的解析表达式。
• 数值模拟用于验证 $1 < \sigma < 2$ 范围内稳定性阈值的存在，识别出分隔稳定与不稳定区域的临界速度 $c = 2z_0\sqrt{\omega}$。

实验结果

研究问题

• RQ1 广义 DNLS 方程的孤立波解在哪些 $\sigma$ 值下是轨道稳定的或不稳定的？
• RQ2 当 $\sigma \in (1,2)$ 时，孤立子速度 $c$ 如何影响孤立波的稳定性？
• RQ3 在 $1 < \sigma < 2$ 范围内，分隔稳定与不稳定孤立波的临界速度阈值 $c_0(\sigma)$ 是什么？
• RQ4 为何 $\sigma < 1$ 的孤立波即使在不同速度下也始终轨道稳定？
• RQ5 gDNLS 方程的稳定性行为与经典 DNLS 方程（$\sigma = 1$）有何不同？

主要发现

• 当 $\sigma \geq 2$ 时，所有孤立波解 $\psi_{\omega,c}$ 均为轨道不稳定，这是由于强非线性导致波包结构失稳。
• 当 $\sigma < 1$ 时，所有孤立波解均为轨道稳定，该结论由变分分析支持，且作用泛函的二阶变分保持正定。
• 在 $1 < \sigma < 2$ 范围内，存在与速度相关的稳定性阈值：当 $c < 2z_0\sqrt{\omega}$ 时孤立波稳定，而当 $c > 2z_0\sqrt{\omega}$ 时不稳定，其中 $z_0 \in (-1,1)$ 为依赖于 $\sigma$ 的常数。
• 数值模拟证实了在 $1 < \sigma < 2$ 情况下存在临界速度 $c_0(\sigma)$，该值将稳定与不稳定区域分隔开，且随着 $\sigma$ 趋近于 2，$z_0(\sigma)$ 递减。
• 稳定性分析依赖于质量 $Q$ 与动量 $P$ 对 $\omega$ 与 $c$ 的导数，这些导数以涉及双曲函数的积分形式表达，并显示与 $\tilde{\kappa}$ 成正比，$\tilde{\kappa}$ 为依赖于 $\sigma$、$\omega$ 与 $c$ 的缩放因子。
• 本文给出了 $\partial_c Q$，$\partial_\omega Q$，$\partial_c P$ 与 $\partial_\omega P$ 的显式解析表达式，其形式涉及 $\alpha_n = \int_0^\infty h^{-1/\sigma - n} dx$，从而可推导出稳定性阈值。