[論文レビュー] 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

ひとこと要約

本稿では、半古典的極限における1次元定常シュレーディンガー方程式を解くための、非重複領域分割法を提案する。この手法は、振動領域では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.

研究の動機と目的

• ポテンシャル関数が符号を変えることで振動領域と減衰領域を有する1次元定常シュレーディンガー方程式を半古典的極限で解く数値的課題に対処すること。
• 問題を振動領域（初期値問題）と減衰領域（境界値問題）に分ける非重複領域分割法を構築すること。
• 振動領域では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)$）を必要とせず、計算コストを顕著に削減できる。
• 両領域にわたる解の振幅と微分が一様に有界であることが示され、ハイブリッドスキームの安定性と収束性を裏付ける。