[论文解读] WKB-method for the 1D Schrödinger equation in the semi-classical limit: enhanced phase treatment

Anton Arnold, C. Klein. B. Ujvari|arXiv (Cornell University)|Aug 6, 2018

Quantum Mechanics and Non-Hermitian Physics参考文献 25被引用 3

一句话总结

本文通过改进误差分析并用谱方法替代数值积分来计算相位，增强了用于求解一维薛定谔方程在半经典极限下的二阶WKB方法。该方法通过消除相位积分中依赖于$1/\varepsilon$的数值误差，在粗网格上实现了均匀准确的高阶收敛，显著提升了小$\varepsilon$时的鲁棒性和精度。

ABSTRACT

This paper is concerned with the efficient numerical computation of solutions to the 1D stationary Schrödinger equation in the semiclassical limit in the highly oscillatory regime. A previous approach to this problem based on explicitly incorporating the leading terms of the WKB approximation is enhanced in two ways: first a refined error analysis for the method is presented for a not explicitly known WKB phase, and secondly the phase and its derivatives will be computed with spectral methods. The efficiency of the approach is illustrated for several examples.

研究动机与目标

• 解决标准数值格式在$\varepsilon \ll 1$时难以解析一维薛定谔方程高度振荡解的低效问题。
• 克服先前基于WKB的方法中数值相位积分误差存在的$1/\varepsilon$依赖性，从而限制精度的问题。
• 通过细化相位和变换步骤中的数值误差分析，提高二阶WKB近似方法的渐近精度。
• 用谱方法替代低阶求积方法（如Simpson法则）来计算相位，以实现指数收敛并减少条件数问题。
• 证明谱方法计算相位可实现宽范围$\varepsilon$值下的均匀准确解，甚至达到机器精度极限。

提出的方法

• 应用二阶WKB近似，通过解析相位和振幅分解，将原始高度振荡的薛定谔方程转化为更平滑的系统。
• 使用WKB假设$\varphi(x) \sim \exp\left(\frac{1}{\varepsilon}\sum_{p=0}^{\infty} \varepsilon^p \phi_p(x)\right)$推导主导阶相位和振幅校正项。
• 采用谱配点法（Clenshaw-Curtis配合重心插值）对相位函数$\tilde{\phi}(x) = \int_0^x \left(\sqrt{a(\tau)} - \varepsilon^2 \beta(\tau)\right) d\tau$实现高精度积分。
• 引入细化的误差分析，同时考虑WKB近似误差和相位中数值积分误差，表明当使用谱方法时，后者不再具有$\mathcal{O}(1/\varepsilon)$的量级。
• 在粗网格上使用标准有限差分法求解经变换后的变量$Z(x)$，利用相位去除后带来的光滑性优势。
• 将谱相位计算与低阶求积方法（如Simpson法则）进行对比，证明其在条件数和数值误差方面表现更优，尤其在$\varepsilon$较小时。

实验结果

研究问题

• RQ1 在基于WKB的一维薛定谔方程求解中，是否能够消除或显著降低相位计算中依赖于$1/\varepsilon$的数值误差？
• RQ2 在半经典极限下，数值积分方法的选择（谱方法 vs. 求积方法）如何影响WKB方法的收敛速率和条件数？
• RQ3 是否通过包含相位和变换步骤中数值误差的细化误差分析，能够获得更准确且可靠的收敛性估计？
• RQ4 谱方法计算相位是否能够实现宽范围$\varepsilon$值下的均匀准确解，包括接近机器精度的值？
• RQ5 考虑到$Z(x)$具有振荡但振幅较小的特性，其简化系统是否适合进一步采用谱方法处理？

主要发现

• 谱方法计算相位消除了早期基于低阶求积的WKB方法中困扰的$\mathcal{O}(1/\varepsilon)$误差项，从而实现均匀准确的解。
• 当$\varepsilon = 10^{-5}$时，谱方法在$h \in [10^{-5}, 10^{-2}]$范围内的误差低于机器精度，表明该方法在双精度下实现了最优精度。
• 在WKB误差饱和水平以上，误差图中的斜率分析证实方法在$h$上实现了二阶收敛。
• 图5右侧（谱相位）显示，在极小$h$时未出现人工误差增长，而Simpson法则情形则存在此问题，归因于更好的条件数和更少的舍入误差。
• 变换变量$Z$的误差满足不包含问题项$E/\varepsilon$的界(3.51)，证实了谱方法的鲁棒性。
• 该方法在精度和收敛速率上均优于一阶WKB格式；二阶WKB误差仅在$h \leq 10^{-4} \text{--} 10^{-3}$时达到饱和，而谱方法在更小的$h$下仍保持高精度。