[论文解读] On soliton (in-)stability in multi-dimensional cubic-quintic nonlinear Schrödinger equations

Rémi Carles, Christian Klein+1 more|arXiv (Cornell University)|Dec 21, 2020

Advanced Mathematical Physics Problems参考文献 33被引用 5

一句话总结

本文研究了二维和三维立方-五次非线性薛定谔方程中基态孤立波的轨道稳定性，其中立方项为聚焦项，五次项为非聚焦项。通过时间分裂谱方法的数值模拟，作者表明当质量超过临界阈值时，作用量最小化的基态是轨道稳定的，而亚临界质量扰动则导致完全色散解。关键发现是，质量过剩的不稳定基态通过辐射演化为稳定、频率更低的孤立波，证实了不同非线性模型中普遍存在的不稳定性场景。

ABSTRACT

We consider the nonlinear Schrödinger equation with a focusing cubic term and a defocusing quintic nonlinearity in dimensions two and three. The core of this article is the notion of stability of solitary waves. We recall the two standard notions of orbital stability in the context of nonlinear Schrödinger equations, and show that they must be considered as independent from each other. We investigate numerically the notion of orbital stability of ground states in the radially symmetric case, confirming existing conjectures or leading to new ones.

研究动机与目标

• 通过区分能量最小化解与作用量最小化解，解决多维立方-五次NLS中基态稳定性分类长期存在的模糊性。
• 在径向对称设置下，数值研究作用量最小化基态的轨道稳定性，特别关注质量相对于临界阈值的作用。
• 验证或修正关于扰动孤立波长期动力学的现有猜想，特别是其渐近收敛至稳定基态的行为。

提出的方法

• 本研究采用时间分裂谱方法，数值求解二维和三维的立方-五次非线性薛定谔方程。
• 基态作为静止方程 $-\frac{1}{2}\Delta\phi + \omega\phi - |\phi|^2\phi + |\phi|^4\phi = 0$ 的解进行计算，频率 $\omega \in (0, 3/16)$。
• 通过施加小的、局域的 $H^1$-正则扰动于已知的作用量最小化基态 $Q_\omega$，并追踪其随时间演化的 $L^\infty$-范数，评估轨道稳定性。
• 将解的最终状态与预先计算的基态库 $Q_{\underline{\omega}}$ 进行比较，以确定稳定孤立波的渐近频率 $\underline{\omega}$。
• 将初始数据的质量与不稳定 $Q_\omega$ 的质量和最终稳定 $Q_{\underline{\omega}}$ 的质量进行比较，以量化辐射损失。

实验结果

研究问题

• RQ1 在多维立方-五次NLS中，作用量最小化基态在何种条件下是轨道稳定的？
• RQ2 初始数据的质量是否是决定扰动孤立波是否保持局域化或完全色散的关键因素？
• RQ3 质量过剩的不稳定基态是否可演化为频率更低、质量更小的稳定基态？
• RQ4 辐射在该系统中扰动孤立波的长期动力学中起何种作用？
• RQ5 是否存在一个普遍的选择机制，决定渐近稳定孤立波的最终频率 $\underline{\omega}}$？

主要发现

• 当初始数据质量超过不稳定基态 $Q_\omega$ 的质量时，解演化为稳定孤立波 $e^{i\underline{\omega}t}Q_{\underline{\omega}}$，且 $\underline{\omega} < \omega$，表明通过辐射损失了质量。
• 当初始质量低于 $M(Q_\omega)$ 时，解完全色散，证实了纯粹的色散行为。
• 数值模拟显示，解的最终 $L^\infty$-范数与稳定基态 $Q_{\underline{\omega}}$ 的 $L^\infty$-范数高度一致，当 $\omega = 0.01$ 时，$\underline{\omega} \approx 0.047$。
• 最终稳定态的质量始终低于初始质量，当 $\omega = 0.01$ 时，$M(Q_{\underline{\omega}}) \approx 77.05$，而 $M(Q_{\omega=0.01}) \approx 79.44$，证实了显著的辐射损失。
• 该不稳定性场景对不同类型的扰动均具有鲁棒性，包括对称与非对称的高斯扰动，且适用于各种 $\omega < \omega_{\text{crit}}$。
• 最终渐近态始终为频率和质量均低于初始不稳定基态的稳定基态，支持关于普遍动力学弛豫的猜想5.1。