[论文解读] Split-step Fourier methods for the Gross-Pitaevskii equation

Juha Javanainen, Janne Ruostekoski|arXiv (Cornell University)|Nov 5, 2004

Photonic and Optical Devices参考文献 2被引用 28

一句话总结

该论文证明，对于 Gross-Pitaevskii 方程（GPE）的分裂步傅里叶方法，只要在非线性项中使用最新更新的波函数，其时间步长精度可与线性薛定谔方程的对应方法保持一致。通过符号计算，作者证明了这一简单规则可保持所有测试的分裂步格式（包括高阶方法）的算法精度阶数，并推测该规则在任意空间维数下对所有最小分裂步方法均普遍成立。

ABSTRACT

We perform a systematic study of the accuracy of split-step Fourier transform methods for the time dependent Gross-Pitaevskii equation using symbolic calculation. Provided the most recent approximation for the wave function is always used in the nonlinear atom-atom interaction potential energy, every split-step algorithm we have tried has the same-order time stepping error for the Gross-Pitaevskii equation and the Schroedinger equation.

研究动机与目标

• 系统分析时间依赖 Gross-Pitaevskii 方程（GPE）的分裂步傅里叶方法的精度，GPE 是超冷原子物理中的关键方程。
• 解决 GPE 数值实现中长期存在的模糊性问题，许多研究组采用临时性、未经验证的方法，缺乏明确的收敛保证。
• 确定线性薛定谔方程中分裂步傅里叶方法所具有的良好精度特性，是否可在扩展至非线性 GPE 时得以保留。
• 识别一种普遍适用、计算成本低廉的规则，确保无论使用何种具体算法，分裂步 GPE 求解器均能实现高阶时间精度。

提出的方法

• 采用时间步长 $ h $ 的符号幂级数展开，系统比较 GPE 的精确时间演化算符与分裂步近似。
• 使用算子代数展开非交换动能与势能算符之和的指数，按 $ h $ 的阶数逐阶追踪对换子项。
• 将分裂步序列中系数 $ \alpha_i $ 和 $ \beta_i $ 的多元多项式方程公式化，以强制在期望阶数内消除误差项。
• 通过将指数展开为幂级数并按顺序作用于波函数来实现分裂步算法，使用 Mathematica 进行符号运算。
• 测试了多种分裂步格式，包括 $ \mathcal{O}(h^3) $、$ \mathcal{O}(h^4) $ 和 $ \mathcal{O}(h^5) $ 方法，且分别以位置步长和动量步长作为起始步骤。
• 通过将一维二阶导数替换为拉普拉斯算子，将该方法推广至多维 GPE，确认其在不同维度下的一致性。

实验结果

研究问题

• RQ1 应用于非线性 Gross-Pitaevskii 方程的分裂步傅里叶方法的时间步长误差阶数是多少？
• RQ2 当相同算法应用于非线性 GPE 时，线性薛定谔方程中分裂步方法的精度特性是否能够保留？
• RQ3 是否存在一种处理 GPE 中非线性项 $ g|\psi|^2 $ 的特定规则，可确保其精度阶数与线性情况相同？
• RQ4 精度是否依赖于步骤顺序（例如，以位置或动量步骤开始）或分裂步格式中系数的选择？
• RQ5 所观察到的精度保持特性是否可推广至任意空间维数和更高阶分裂步方法？

主要发现

• 只要在非线性势能项 $ g|\psi|^2 $ 中使用最新可用的波函数，GPE 的分裂步傅里叶方法即可达到与线性薛定谔方程相同的时序误差阶数。
• 对于所有测试的分裂步格式——包括 $ \mathcal{O}(h^3) $、$ \mathcal{O}(h^4) $ 和 $ \mathcal{O}(h^5) $——条件 $ c_0 = 0 $、$ |c_1| = 1 $ 在非线性项中确保了 $ \mathcal{O}(h^3) $ 的精度，且 $ |\psi|^2 $ 中使用的是最新 $ \psi $。
• 使用最新波函数的要求不仅充分，而且在所测试的格式中似乎也是必要条件，以在非线性情况下维持与线性情况相同的误差阶数。
• 作者推测，这一规则——即在 $ |\psi|^2 $ 中使用最新 $ \psi $——对所有最小分裂步方法普遍成立，无论阶数或维度如何。
• 该方法在二维和三维空间中也保持精度，通过显式验证 $ \mathcal{O}(h^3) $ 的三指数格式得到确认。
• 这种鲁棒性的根本原因仍为代数性且尚未被解释，提示在保持范数的非线性微分方程中可能存在更深层次的数学结构。