[论文解读] Stationary Schr\\"odinger equation in the semi-classical limit: numerical coupling of oscillatory and evanescent regions

Anton Arnold, Claudia Negulescu|arXiv (Cornell University)|Jun 16, 2016

Electromagnetic Simulation and Numerical Methods参考文献 21被引用 6

一句话总结

本文提出了一种新颖的非重叠域分解方法，用于求解半经典极限下的1D定态薛定谔方程，结合了在振荡区域使用的WKB推进格式与在倏逝区域使用的基于WKB的有限元法（FEM）。该方法首次提供了对WKB-FEM的完整误差分析，并证明了混合方法具有与刚性无关的渐近保持性，从而可在不解析快速振荡或边界层的情况下，于粗网格上获得精确解。

ABSTRACT

This paper is concerned with a 1D Schr\\"odinger scattering problem involving both oscillatory and evanescent regimes, separated by jump discontinuities in the potential function, to avoid "turning points". We derive a non-overlapping domain decomposition method to split the original problem into sub-problems on these regions, both for the continuous and afterwards for the discrete problem. Further, a hybrid WKB-based numerical method is designed for its efficient and accurate solution in the semi-classical limit: a WKB-marching method for the oscillatory regions and a FEM with WKB-basis functions in the evanescent regions. We provide a complete error analysis of this hybrid method and illustrate our convergence results by numerical tests.

研究动机与目标

• 为解决当势函数变号时，1D定态薛定谔方程在半经典极限下所面临的数值挑战，此时会同时出现振荡区域与倏逝区域。
• 开发一种非重叠域分解方法，将问题划分为振荡区域（初值问题）与倏逝区域（两点边值问题）子域。
• 设计一种混合数值格式，在振荡区域使用WKB推进法，在倏逝区域使用WKB-FEM，确保在粗网格上保持精度。
• 为倏逝区域中的WKB-FEM方法提供完整的误差分析，该分析在以往文献中尚属空白。
• 证明混合方法具有渐近保持性，且收敛性不依赖于小参数 $\varepsilon$。

提出的方法

• 采用非重叠域分解方法，根据势函数 $a(x)$ 的符号将问题划分为振荡区域与倏逝区域子域，避免在拐点处处理。
• 在振荡区域（$a(x) > 0$），采用WKB推进格式以解析方式消除主导振荡，得到一个更平滑的残差问题，可在粗网格上求解。
• 在倏逝区域（$a(x) < 0$），采用使用WKB基函数（指数型）的有限元法，以捕捉边界层而无需精细网格划分。
• 精心构造界面条件，以确保解及其导数在子域边界处连续，从而支持单次扫描求解策略。
• 通过最后的缩放步骤，在连续层面恢复精确的解析解，确保各区域间的一致性。
• 通过与刚性无关的误差估计建立了理论收敛性，并由数值实验加以验证。

实验结果

研究问题

• RQ1 非重叠域分解方法能否在不解析小尺度结构的前提下，准确耦合半经典薛定谔方程中振荡与倏逝区域？
• RQ2 在倏逝区域中，基于WKB的有限元法是否能实现与 $\varepsilon$ 无关的最优收敛性？
• RQ3 混合WKB推进法与WKB-FEM方法能否实现渐近保持性，且误差界在 $\varepsilon$ 上一致有界？
• RQ4 应采用何种界面条件，才能确保在单次扫描中精确恢复连续解？
• RQ5 如何为倏逝区域中的WKB-FEM建立完整的误差分析，而该分析在以往文献中尚属缺失？

主要发现

• 所提出的非重叠域分解方法在应用于连续问题时，通过适当的界面条件与缩放处理，可在单次扫描中精确恢复解析解。
• 混合WKB推进法与WKB-FEM方法实现了渐近保持性，其收敛误差不依赖于半经典参数 $\varepsilon$。
• 首次为倏逝区域中的WKB-FEM提供了完整的误差分析，结果表明在粗网格上可实现最优收敛率。
• 数值实验验证了理论误差界，并表明该方法在 $\varepsilon$ 极小时仍具有鲁棒性与高精度。
• 通过利用WKB渐近性质，该方法避免了对精细网格（$h = \mathcal{O}(\varepsilon)$）的需求，显著降低了计算成本。
• 结果表明解的振幅与导数在两个区域中均保持一致有界，支持了混合格式的稳定性和收敛性。