分数阶非线性薛定谔方程数值逼近的高阶方法

《Mathematics and Computers in Simulation》：A high-order method for the numerical approximation of fractional nonlinear Schr?dinger equations

　　

在本论文中，分数阶非线性薛定谔（fNLS）方程的周期初值问题在空间上通过傅里叶谱伽辽金方法进行离散，在时间上通过基于隐式中点法则（IMR）组合的对角隐式高阶龙格-库塔格式进行离散。研究人员证明了空间半离散化和全离散化的一些性质和误差估计。通过多个数值实验展示了该格式的收敛结果和整体性能。

分数阶非线性薛定谔（fNLS）方程最初由Laskin在分数阶量子力学中引入，并应用于非线性光学和水波模型。该方程的周期初值问题在数值逼近方面，现有方法如Crank-Nicolson（CN）有限差分、傅里叶谱方法等，虽能保持部分守恒量，但在时间精度上多为二阶，且缺乏高精度且保持几何结构（如辛结构）的数值方案。此外，对于分数阶情形（0
研究人员开展的研究包括：对fNLS方程的周期初值问题，在空间上采用傅里叶谱伽辽金（FSG）方法进行半离散化，在时间上采用Yoshida提出的基于IMR组合的高阶对角隐式RK方法（系列包括q=1阶的IMR本身和q=3阶的四阶方法）。通过理论分析，证明了半离散系统在L²和H^s（s∈(1/2,1]）范数下的误差估计分别为O(N^-μ)和O(N^s-μ)，其中N为三角多项式逼近次数，μ为解的正则性指数。对于全离散格式，在足够小的时间步长k下，证明了L²和H^s范数下的误差估计分别为O(k^α+N^-μ)和O(k^α+N^s-μ)，其中α由时间方法的阶数决定（α=2或4）。数值实验验证了这些收敛阶，并展示了长期模拟中误差的线性增长（与孤立波参数误差相关）以及三个守恒量（质量、动量、哈密顿量）的几乎保持。

主要关键技术方法包括：傅里叶谱伽辽金（FSG）空间离散方法，将方程投影到三角多项式空间S_N，利用FFT实现；基于隐式中点法则（IMR）的组合高阶龙格-库塔（RK）时间离散方法（Yoshida方法，q=1或3阶），每个阶段通过迭代求解固定点系统（如(3.43)所示）来获得内部解；迭代求解采用对角化处理，依赖于傅里叶表示。论文中未涉及生物样本或队列来源。

**2. 谱半离散化（空间）**
**2.1 半离散化的呈现与性质**
通过傅里叶谱伽辽金方法，将周期初值问题离散为有限维系统（2.7）。研究人员证明了半离散解保持三个守恒量：质量I₁、动量I₂和哈密顿量H（命题2.1）。在s>1/2条件下，通过先验估计得到半离散解在H^s范数下的有界性（引理2.1）。

**2.2 半离散化的收敛性**
假设精确解u∈H^μ（μ>s），利用投影性质和Gronwall引理，研究人员证明了半离散解在L²范数下的误差为O(N^-μ)（定理2.2），在H^s范数下的误差为O(N^s-μ)（定理2.1），体现了谱精度。进而得到半离散解时间导数的有界性（命题2.2），为全离散分析提供基础。

**3. 全离散化**
**3.1 全离散格式与初步性质**
采用Yoshida提出的基于IMR组合的高阶RK方法（如(3.2)所示），并证明了全离散解保持质量I₁和动量I₂（引理3.1），但哈密顿量H的保持不是理论保证的，期望误差呈指数小。通过Brouwer不动点定理证明了内部阶段解的存在性（命题3.1），并在条件（3.7）下证明了唯一性及H^s范数的有界性（引理3.2）。

**3.2 全离散格式的收敛性**
通过引入局部时间误差Θ_n，并假设其满足O(k^α+1)（α=2或4，由方法阶数决定），研究人员在L²和H^s范数下建立了全离散解与半离散解之间的误差传递，结合半离散误差，得到全离散解的总体误差估计：在L²范数下为O(k^α+N^-μ)（定理3.1），在H^s范数下为O(k^α+N^s-μ)（定理3.2）。这些结果依赖于解的正则性足够高以及时间步长k足够小以满足条件（3.7）。

**4. 数值实验**
以非分数阶情形（s=1）的孤立波解（4.5）为参考，验证了时间精度：IMR（q=1）达到二阶，三阶段方法（q=3）达到四阶（表1、表2）。长期模拟（T=100）中，L²误差随时间线性增长，符合理论预期（由于参数误差的相位漂移）。对于分数阶情形（s=0.75），采用数值生成的近似孤立波作为初始条件，展示了守恒量误差随时间几乎保持（图3、图4），且振幅和速度误差几乎恒定（图5、图6），表明该方法在长期模拟中具有高精度和几何保持性质。

讨论部分总结：数值实验证实了理论收敛阶，并揭示了长期模拟中误差的线性增长源于孤立波参数（振幅、速度、相位）的误差，而质量、动量守恒意味着振幅和速度几乎保持，只有相位误差随时间线性增长。对于分数阶情形，类似性质成立，且哈密顿量误差也几乎保持。这些结果证明该全离散格式对于fNLS方程的长期数值模拟是可靠且高效的。

研究结论部分翻译：本文提出的全离散格式在L²和H^s（s∈(1/2,1]）范数下分别达到O(k^α+N^-μ)和O(k^α+N^s-μ)的收敛阶，其中α由时间方法阶数决定（α=2或4）。数值实验验证了二阶和四阶时间精度，以及长期模拟中误差的线性增长和守恒量（质量、动量、哈密顿量）的几乎保持，为分数阶非线性薛定谔方程（尤其是孤立波动力学）的长期高精度模拟提供了可靠数值工具。