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

ひとこと要約

本稿では、1次元シュレーディンガー方程式の半古典的極限における2次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$ 依存の数値位相積分誤差を克服する。
• 位相および変換ステップにおける数値誤差を含めた誤差解析を精緻化することで、2次WKB近似の漸近的精度を向上させる。
• 低次の数値積分（例：シンプソン則）を代替して、指数的収束と条件数の低減を実現するスレクトル法を位相計算に導入する。
• スレクトル法による位相計算が、機械精度の限界に近い $\varepsilon$ 値を含む広範な $\varepsilon$ 値範囲で一様に高精度な解を可能にすることを示す。

提案手法

• 元の高振動的シュレーディンガー方程式を、解析的位相および振幅分解を用いて滑らかさの高い系に変換する2次WKB近似を適用する。
• WKBアンサツト $\varphi(x) \sim \exp\left(\frac{1}{\varepsilon}\sum_{p=0}^{\infty} \varepsilon^p \phi_p(x)\right)$ を用いて、一次近似の位相および振幅補正を導出する。
• スレクトルコロケーション（バーサインチ・カーティス法にバーサインチ補間を組み合わせ）を用いて、位相関数 $\tilde{\phi}(x) = \int_0^x \left(\sqrt{a(\tau)} - \varepsilon^2 \beta(\tau)\right) d\tau$ の高精度積分を実装する。
• WKB近似誤差と位相における数値積分誤差の両方を考慮した精緻な誤差解析を導入し、スレクトル法を用いることで、後者に $\mathcal{O}(1/\varepsilon)$ の誤差項が生じないことを示す。
• 位相除去によって滑らかさが向上した変数 $Z(x)$ について、粗いメッシュ上で標準的な有限差分法を用いて変換系を解く。
• スレクトル位相計算と低次の数値積分（シンプソン則）を比較し、特に $\varepsilon$ が小さい場合に、より良好な条件数と低減された数値誤差を示す。

実験結果

リサーチクエスチョン

• RQ1 1次元シュレーディンガー方程式のWKBベーススキームにおいて、位相計算の $1/\varepsilon$ 依存誤差を排除または著しく低減できるか？
• RQ2 数値積分法の選択（スレクトル法対数値積分）が、半古典的極限におけるWKB法の収束速度および条件数にどのように影響するか？
• RQ3 位相および変換ステップにおける数値誤差を含めた精緻な誤差解析は、より正確で信頼性の高い収束見積もりをもたらすか？
• RQ4 スレクトル法による位相計算により、機械精度の限界に近い $\varepsilon$ 値を含む広範な $\varepsilon$ 値範囲で一様に高精度な解が得られるか？
• RQ5 振動的ではあるが振幅が小さい性質を持つ $Z(x)$ の変換系は、さらなるスレクトル処理に適しているか？

主な発見

• スレクトル位相計算により、従来の低次の数値積分を用いたWKBスキームで問題となっていた $\mathcal{O}(1/\varepsilon)$ 誤差項が排除され、一様に高精度な解が得られるようになった。
• $\varepsilon = 10^{-5}$ の場合、スレクトル法における誤差は $h \in [10^{-5}, 10^{-2}]$ の範囲ですでに機械精度以下に低下し、倍精度計算で最適な精度が達成されたことを示している。
• WKB誤差の飽和レベルを超える $h$ に対して、誤差プロットの傾き解析により、$h$ について2次収束が確認された。
• 図5の右図（スレクトル位相）では、非常に小さな $h$ に対しても人工的な誤差増加が見られず、シンプソン則のケースとは異なり、より良好な条件数と低減された丸め誤差に起因する。
• 変換変数 $Z$ の誤差は、問題となる $E/\varepsilon$ 項を含まない境界 (3.51) を満たしており、スレクトル法のロバスト性が裏付けられている。
• 1次WKBスキームに比べ、本手法は精度および収束速度の両面で優れており、2次WKB誤差は $h \leq 10^{-4} \text{--} 10^{-3}$ の範囲でのみ飽和に達するが、スレクトル法ではより小さな $h$ に対しても高精度を維持している。