[论文解读] High order exponential integrators for nonlinear Schrödinger equations with application to rotating Bose-Einstein condensates

Christophe Besse, Guillaume Dujardin+1 more|arXiv (Cornell University)|Jul 2, 2015

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

一句话总结

本文提出高阶指数积分格式——具体为指数 Runge-Kutta (ERK) 与 Lawson 方法——用于求解描述旋转玻色-爱instein凝聚体的非线性薛定谔方程。通过变量变换消除旋转项，作者实现了高阶精度与优异的收敛特性，在长时模拟中，尤其在误差常数较低时，其精度与计算效率均优于经典分裂步方法。

ABSTRACT

This article deals with the numerical integration in time of nonlinear Schrödinger equations. The main application is the numerical simulation of rotating Bose-Einstein condensates. The authors perform a change of unknown so that the rotation term disappears and they obtain as a result a nonautonomous nonlinear Schrödinger equation. They consider exponential integrators such as exponential Runge--Kutta methods and Lawson methods. They provide an analysis of the order of convergence and some preservation properties of these methods in a simplified setting and they supplement their results with numerical experiments with realistic physical parameters. Moreover, they compare these methods with the classical split-step methods applied to the same problem.

研究动机与目标

• 为旋转玻色-爱因斯坦凝聚体 (BEC) 动力学中出现的非自治非线性薛定谔方程开发高阶时间积分方法。
• 在简化的 Sobolev 空间框架下，分析指数 Runge-Kutta 与 Lawson 方法的收敛性与守恒性质。
• 从精度、守恒不变量与计算成本三个方面，比较这些指数方法与经典分裂步方法的性能。
• 在具有旋转与各向异性势阱的二维真实 BEC 模拟中，展示方法的有效性。
• 证明指数积分格式可避免分裂方法在长时间积分中常见的共振现象。

提出的方法

• 通过变量变换将旋转 BEC 方程转化为消除旋转项的形式，得到非自治非线性薛定谔方程。
• 对变换后的方程应用指数 Runge-Kutta (ERK) 与 Lawson 方法，利用线性算子的指数函数精确积分线性部分。
• 采用周期性环面域 $\mathbb{T}_\delta^d$ 以避免边界条件，并通过傅里叶方法简化拉普拉斯算子的谱计算。
• 实现阶数为 6 的高斯型 ERK 与 Lawson 方法，且具有高阶段阶数，以确保低误差常数与高精度。
• 采用伪谱空间半离散化方法，结合傅里叶谱方法，以实现高空间精度，尤其适用于二维模拟。
• 通过误差范数 $\mathcal{E}_{P,h}$、$\mathcal{E}_{M,h}$ 与 $\mathcal{E}_{E,h}$ 分别衡量动量、质量与能量的误差，在一维与二维测试案例中进行方法比较。

实验结果

研究问题

• RQ1 高阶指数积分格式是否能在旋转 BEC 模拟中实现优于经典分裂步方法的精度？
• RQ2 尽管并非精确守恒，指数 Runge-Kutta 与 Lawson 方法在质量及其他不变量的保持上是否能与分裂格式达到相当的精度？
• RQ3 在非自治、具有物理意义的设定下，ERK 与 Lawson 方法的误差常数与收敛率与分裂方法相比如何？
• RQ4 指数积分格式是否能避免分裂方法在长时间积分中出现的共振现象？
• RQ5 对于具有现实参数的二维 BEC 问题，隐式指数方法的计算效率相对于显式分裂方法如何？

主要发现

• 高斯-ERK 方法的误差常数显著低于 Lawson 与分裂方法，在所有测试配置中均表现出更优的有限时间精度。
• 当时间步长为 $10^{-4}$ 量级时，高斯-ERK 方法的相对质量误差与质量守恒的 Lawson 及分裂格式相当，尽管其本身并非精确守恒质量。
• 尽管具有隐式性质，高斯-ERK 方法在 1D 模拟中的计算时间仍低于 Lawson 与分裂方法。
• 在长时间模拟中，高斯-ERK 方法未表现出共振现象，而分裂方法则出现了此类不稳定性。
• 在具有真实参数（$\beta=1000$，$\Omega=0.9$，$\gamma_x=1.05$，$\gamma_y=0.95$）的二维模拟中，高斯-ERK 方法成功再现了先前工作中观察到的涡旋晶格结构。
• 该方法对一般约束势（包括非对称二次势与二次-四次势）均有效，而伪谱方法在此类情形下会失效。