[論文レビュー] Very-High-Precision Normalized Eigenfunctions For A Class Of Schrödinger Type Equations

Amna Noreen, Kåre Olaussen|arXiv (Cornell University)|Apr 22, 2011

Electromagnetic Simulation and Numerical Methods参考文献 3被引用数 5

ひとこと要約

本稿では、非常に高精度な正規化定数（最大100万桁以上）を備えたシュレーディンガー型方程式の正規化固有関数を計算する線形スケーリングアルゴリズムを提示する。台形則とWKB近似、およびポアソン再summationを組み合わせることで、必要な固有関数評価回数が所望の精度 $P$ に対して線形に増加することを保証し、量子行列要素に不可欠な正規化積分の高精度計算を効率的に行える。

ABSTRACT

We demonstrate that it is possible to compute wave function normalization constants for a class of Schr¨odinger type equations by an algorithm which scales linearly (in the number of eigenfunction evaluations) with the desired precision P in decimals.

研究の動機と目的

• 波動関数正規化定数を極めて高精度（最大100万桁以上）で計算する手法を開発すること。
• 高精度な量子力学的計算に不可欠な正規化積分を正確に評価する計算課題に対処すること。
• 必要な固有関数評価回数が所望の精度 $P$ に対して線形に増加することを示し、高精度計算を効率的に行えるようにすること。
• 台形則とWKB近似を活用し、誤差を推定し、最小誤差を達成するための最適なステップサイズ $h$ と積分範囲 $x_{\text{max}}$ を最適化すること。
• 量子系における行列要素を極めて高精度に評価するために不可欠な正規化積分を信頼性高く計算可能にする。

提案手法

• 数値積分に台形則を用い、滑らかで急速に減衰する被積分関数に対して、オイラー＝マクローリンおよびポアソン再summation解析により、これが最適であると正当化する。
• WKB近似を用いて $\psi(x)^2$ のフーリエ変換および $\psi(x)^2$ の漸近的減衰を推定し、有限な $h$ および $x_{\text{max}}$ に対する誤差推定を可能にする。
• 切捨て誤差と尾部誤差をバランスさせることで最適なステップサイズ $h$ を導出し、$x^4$ ポテンシャルに対しては $h \approx 1.58 (M+1)^{-2/3}$、二重井戸ポテンシャルに対しては $h \approx 1.73\, s^{1/3} P^{-2/3}$ が得られる。
• 積分範囲 $x_{\text{max}}$ は $\psi(x)^2 \lesssim 10^{-P}$ となるように選ぶ。二重井戸の場合、$x_{\text{max}} \approx 1.51\, s^{1/3} P^{1/3}$ が得られる。
• 必要な評価回数 $M+1$ は、$x^4$ ポテンシャルおよび二重井戸ポテンシャルの両方で、漸近的に $0.87\,P$ に比例する。これにより $P$ に対して線形スケーリングが保証される。
• 高精度算術を全般に使用し、乗算コストは $P^{1.6}$ から $P \log P \log \log P$ の間で変動し、合計時間は $P^3$ 稍速くスケーリングする。

実験結果

リサーチクエスチョン

• RQ1 シュレーディンガー固有関数の正規化積分を、固有関数評価回数に比例して精度を上げられるか？
• RQ2 高精度正規化において誤差を最小化するための最適なステップサイズ $h$ および積分範囲 $x_{\text{max}}$ をどのように決定できるか？
• RQ3 WKB近似が、誤差推定のための $\psi(x)^2$ のフーリエ変換および漸近的減衰をどれほど正確に記述できるか？
• RQ4 異なるポテンシャル（$x^4$ および $(x^2-1)^2$）に対して、所望の精度 $P$ に対する必要な評価回数はどのようにスケーリングするか？
• RQ5 本手法は、$x^4$ および二重井戸系の例で示されたように、高励起状態および degenerate 固有状態に対しても高精度で拡張可能か？

主な発見

• 所望の精度 $P$ に達するための固有関数評価回数は $P$ に対して線形に増加し、$x^4$ ポテンシャルおよび二重井戸ポテンシャルの両方で $M+1 \approx 0.87\,P$ に比例する。
• $x^4$ ポテンシャルにおいて、基底状態の固有エネルギーを100万桁以上、50,000番目の励起状態を50,000桁以上で計算した。
• $s = 1/50,000$ の二重井戸ポテンシャルにおいて、2つの最低固有エネルギーを30,000桁以上で計算し、それらの差は10次WKB予測と48桁一致した。
• WKBおよび鞍点解析から導かれた誤差推定値 $\varepsilon(M) \approx \text{e}^{-2.64(M+1)}$ は、数値結果と整合的である。
• 計算された正規化積分の精度は、推定値から1〜2小数桁以内に収束しており、誤差モデルの妥当性が裏付けられた。
• 高精度算術の影響により、アルゴリズムの合計実行時間は $P^3$ 稍速くスケーリングするが、独立した固有関数評価の並列化が可能である。