[论文解读] High order explicit symplectic integrators for the Discrete Non Linear Schrödinger equation

Jehan Boreux, Timotéo Carletti+1 more|arXiv (Cornell University)|Dec 15, 2010

Numerical methods for differential equations参考文献 11被引用 3

一句话总结

本文提出了一类用于离散非线性薛定谔方程（DNLS）的高阶显式辛积分器，基于 Yoshida 的组合方法，实现了长期数值模拟中极佳的能量与质量守恒。该方法可达到最高 10 阶精度，同时保持显式与计算高效，相较于非辛及隐式格式在长时间模拟中展现出更高的精度与稳定性。

ABSTRACT

We propose a family of reliable symplectic integrators adapted to the Discrete Non-Linear Schrödinger equation; based on an idea of Yoshida (H. Yoshida, Construction of higher order symplectic integrators, Physics Letters A, 150, 5,6,7, (1990), pp. 262.) we can construct high order numerical schemes, that result to be explicit methods and thus very fast. The performances of the integrators are discussed, studied as functions of the integration time step and compared with some non symplectic methods.

研究动机与目标

• 为离散非线性薛定谔方程（DNLS）开发高阶显式辛积分器，该方程缺乏可分离的哈密顿结构，因此难以应用标准辛方法。
• 克服隐式辛方法（如 Gauss-Radau、基于生成函数的方案）的局限性，这些方法需要昂贵的牛顿迭代，不适用于大规模 DNLS 模拟。
• 确保长期保持关键物理不变量——能量、质量以及共轭关系 $\overline{p} = q$——这对研究孤子态和负温度区域等渐近动力学至关重要。
• 通过优化时间步长与积分器阶数，在计算效率与精度之间取得平衡，最小化 CPU 成本的同时保持机器精度级别的能量守恒。
• 证明尽管系统不具备准可积性，基于 Yoshida 组合法构造的显式高阶辛方案在 DNLS 中依然可行且表现更优。

提出的方法

• 将 Yoshida 的组合方法改进以构造 DNLS 的高阶辛积分器，通过递归组合二阶辛映射实现 2、4、6、8、10 阶精度。
• 将 DNLS 系统的哈密顿量定义为 $ H = -i \sum_{l=1}^{N} \left[ \frac{(p_{l+1}-p_l)(q_{l+1}-q_l)}{h^2} - p_l^2 q_l^2 \right] $，该表达式无法分离为动能与势能项，因此标准分裂方法不可行。
• 构造积分器映射 $ \exp(L_A) $ 与 $ \exp(L_B) $，其中 $ L_A $ 和 $ L_B $ 分别代表哈密顿流的线性部分与非线性部分，利用向量-矩阵乘积计算 $ \exp(L_B) $。
• 通过使用优化系数 $ (c_j, d_j) $ 的 $ \exp(L_A) $、$ \exp(L_B) $ 及其时间反演形式的组合实现积分器，确保辛结构与高阶精度。
• 通过采用 Yoshida 方法的改进版本（$ Y^{\text{opt}}_{2m} $）优化方案，减少组合步数，从而降低 CPU 成本，同时不损失阶数。
• 利用 $ \exp(L_B) $ 的类傅里叶结构，通过基于快速傅里叶变换（FFT）的方法加速计算，将复杂度从 $ \mathcal{O}(N^2) $ 降低至 $ \mathcal{O}(N \log N) $。

实验结果

研究问题

• RQ1 能否为缺乏可分离哈密顿结构、因而难以应用标准辛分裂方法的 DNLS 系统构造高阶显式辛积分器？
• RQ2 所提出的辛积分器在长时间积分下，其能量与质量守恒性能相较于非辛与隐式格式如何？
• RQ3 在实现期望精度与能量守恒的前提下，积分器阶数、时间步长与计算成本（CPU 时间）之间的权衡关系如何？
• RQ4 若初始时刻满足共轭关系 $ \overline{p(t)} = q(t) $，数值方案在多大程度上能保持该关系？
• RQ5 所提方法能否在使用大时间步长时实现机器精度级别的能量守恒，从而高效模拟 DNLS 中的长时间渐近动力学？

主要发现

• 所提出的高阶显式辛积分器（$ Y_{2m} $）在时间步长 $ \tau \sim 10^{-3} $ 下使用 $ Y_4 $，或在 $ \tau \sim 10^{-2} $ 下使用 $ Y_6 $ 时，总能量的相对误差低于 $ 10^{-15} $，展现出极佳的长期能量守恒性能。
• 质量积分 $ I = \sum_l p_l q_l $ 以接近机器精度的方式被保持，证实了该方法在维持关键物理不变量方面的有效性。
• 若初始时刻满足 $ \overline{p(t)} = q(t) $，则该共轭关系在数值上被保持，这对模拟玻色-爱因斯坦凝聚体等物理系统至关重要。
• CPU 时间随时间步长 $ \tau $ 的倒数近似线性增长（固定阶数），随阶数 $ 2m $ 呈 $ 3^m $ 的指数增长，但优化后的 $ Y^{\text{opt}}_{2m} $ 版本相比标准 $ Y_{2m} $ 将成本分别降低 1.5 倍（$ Y_6 $）与 2.2 倍（$ Y_8 $）。
• 在等效精度下，$ Y_6(\tau = 10^{-2}) $ 所需 CPU 时间少于 $ Y_4(\tau = 10^{-3}) $，且 $ Y_8(\tau = 2 \times 10^{-2}) $ 比 $ Y_6(\tau = 10^{-2}) $ 更高效，表明高阶格式在较大时间步长下更具效率。
• 修正的四阶 Runge-Kutta 方法比 $ Y_4 $ 更快，但比 $ Y^{\text{opt}}_6 $ 更慢；而标准 RK4 显著慢于 $ Y_4 $，凸显了辛方法在效率与精度方面的优越性。