[论文解读] Numerical computation of soliton dynamics for NLS equations in a driving potential

Marco Caliari, Marco Squassina|ArXiv.org|Aug 25, 2009

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

一句话总结

本文对具有时间不变外部势的非线性薛定谔方程中的孤立子动力学进行了数值研究，表明在半经典极限（$\varepsilon \to 0$）下，解的质心遵循由牛顿方程 $\ddot{x} = -\nabla V(x)$ 所控制的经典轨迹。通过谱方法和指数积分器，作者验证了孤立子动力学行为，并确认了在二维谐振势下，$\mathbb{H}_\varepsilon$ 范数中的 $\mathcal{O}(\varepsilon)$ 误差界。

ABSTRACT

We provide some numerical computations for the soliton dynamics of the nonlinear Schrödinger equation with an external potential. After computing the ground state solution $r$ of a related elliptic equation we show that, in the semi-classical regime, the center of mass of the solution with initial datum modelled on $r$ is driven by the solution of a Newtonian type law. Finally, we provide some examples and analyze the numerical errors in the two dimensional case when $V$ is an harmonic potential.

研究动机与目标

• 在半经典区域研究带外部势的非线性薛定谔方程中的孤立子动力学。
• 通过变分最小化计算相关椭圆方程的基态解。
• 验证解析预测：解的质心遵循经典轨迹 $\ddot{x} = -\nabla V(x)$。
• 分析二维谐振势下模拟中的数值误差。

提出的方法

• 通过在固定 $L^2$-范数约束下对能量泛函进行变分最小化，数值计算基态 $r$。
• 在空间中使用谱配点法，结合傅里叶谱方法以实现高精度。
• 基于 $\varphi_1$ 和 $\varphi_2$ 函数应用指数积分器进行时间积分，确保二阶精度。
• 通过 $\hat{R}_{n+1} - \hat{R}_{n2}$ 的误差估计实现自适应时间步长，以提高效率。
• 通过变量变换对问题进行缩放，以在 $\varepsilon \to 0$ 极限下实现稳定的数值计算。
• 通过 $\mathbb{H}_\varepsilon$ 范数进行误差分析，其定义为 $\|\phi\|_{\mathbb{H}_\varepsilon}^2 = \varepsilon^{2-N}\|\nabla\phi\|_{L^2}^2 + \varepsilon^{-N}\|\phi\|_{L^2}^2$，以量化收敛性。

实验结果

研究问题

• RQ1 在半经典极限下，带势的非线性薛定谔方程解的质心是否遵循经典轨迹 $\ddot{x} = -\nabla V(x)$？
• RQ2 相关椭圆方程的基态解在数值上能多准确地计算？
• RQ3 在 $\mathbb{H}_\varepsilon$ 范数下，用孤立子轮廓 $r((x - x(t))/\varepsilon)$ 近似解的数值误差量级是多少？
• RQ4 势的选择（例如，频率比为有理数或无理数的谐振势）如何影响孤立子质心的长期动力学？
• RQ5 基于残差误差估计的自适应时间步长是否能提高长时间模拟中的计算效率？

主要发现

• 初始数据以基态 $r$ 建模时，解的质心在半经典区域遵循经典轨迹 $\ddot{x} = -\nabla V(x)$。
• 数值模拟证实，真实解与孤立子轮廓 $r((x - x(t))/\varepsilon)$ 之间的 $\mathbb{H}_\varepsilon$-范数误差被限制在 $\mathcal{O}(\varepsilon)$ 以内，与解析理论预测一致。
• 对于频率比为有理数的二维谐振势，孤立子轨迹是周期性的；对于无理数频率比，轨迹是遍历的。
• 双峰初始条件导致孤立子碰撞，数值格式准确捕捉了动力学行为。
• 采用自适应时间步长的指数积分器在空间上实现谱精度，在时间上实现二阶精度，误差估计使计算更加高效。