[論文レビュー] Numerov numerical method applied to the Schrödinger equation

Francisco Caruso, V. Oguri|arXiv (Cornell University)|Mar 17, 2014

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

ひとこと要約

本稿は、1次元調和振動子および水素原子に対する量子力学における固有値問題を解くために、Numerov法の教育的応用を提示している。標準的なNumerov法を、シュレーディンガー方程式に含まれる1階微分項を扱えるように拡張し、数値的結果を解析的解と比較して検証した。2次精度の有限差分スキームを用いて、エネルギー準位および径方向波動関数の計算において高い精度を達成した。

ABSTRACT

In this paper it is shown how to solve numerically eigenvalue problems associated to second order linear ordinary differential equations, containing also terms which depend on the variable. A didactic presentation of the Numerov Method is given and, in the sequel, it is applied to two quantum non-relativistic problems with well known analytical solutions: the simple harmonic oscillator and the hydrogen atom. The numerical results are compared to those obtained analytically.

研究の動機と目的

• 量子力学における2階線形常微分方程式を解くNumerov法の教育的導入を提供すること。
• 波動関数の1階微分項を含む場合に対応できるように、標準的なNumerov法を拡張すること。
• 2つのベンチマーク量子系（調和振動子および水素原子）の時間に依存しないシュレーディンガー方程式を数値的に解くこと。
• 数値的固有値および波動関数を解析的解と比較し、手法の精度を検証すること。
• 現代の量子力学的計算における、この手法の頑健さおよび一般応用可能性を示すこと。

提案手法

• シュレーディンガー方程式を、1階微分の明示的使用を避ける2次精度の有限差分スキームとして離散化するためにNumerov法を適用した。
• 波動関数の1階微分に依存する項を含めるために、修正されたNumerov公式を導出した。これにより、角運動量を有する水素原子への適用が可能になった。
• 3点ステンシルを用いて初期条件から波動関数を前向きおよび後向きに発展させ、安定性と精度を確保した。
• 空間領域の両端（例えば、原点および無限遠）に境界条件を適用し、固有値問題を解いた。
• コードはROOTフレームワークを用いたC++で実装され、公開リンクから入手可能である。
• 数値的エネルギー準位および径方向波動関数を正確な解析的結果と比較することで、手法の妥当性を検証した。

実験結果

リサーチクエスチョン

• RQ1 1階微分項を含むように拡張されたNumerov法は、1次元調和振動子のエネルギー固有値をどの程度正確に計算できるか？
• RQ2 修正されたNumerov法は、3次元の水素原子の径方向波動関数およびエネルギー準位を正しく再現できるか？
• RQ3 調和振動子および水素原子の両者について、数値的固有値と解析的解との間の数値誤差はどの程度か？
• RQ4 角運動量がゼロでない場合のシュレーディンガー方程式を含む特異境界値問題を解く際、この手法はどの程度の性能を示すか？
• RQ5 Numerov法は、標準的な教科書的ケースを超えて、現代の量子力学的問題へどの程度一般化され、応用可能か？

主な発見

• 1次元調和振動子および水素原子（l=1）の数値的エネルギー準位は、解析的解と非常に近く一致しており、偏差は2％未満であった。
• 水素原子（l=1）について、計算されたエネルギー準位はE2 = -3.47 eV（解析的：-3.40 eV）、E3 = -1.54 eV（解析的：-1.51 eV）、E4 = -0.83 eV（解析的：-0.85 eV）であった。
• 数値的に得られた径方向波動関数は、形状および節の構造において解析的解と一致しており、手法の実装の正しさを裏付けた。
• 有効ポテンシャルおよび径方向微分を組み込むことで、修正されたNumerov法は水素原子のシュレーディンガー方程式における1階微分項を効果的に処理できた。
• この手法は、原点における特異性を含む広大な空間領域においても2次収束性および安定性を示した。
• ROOTフレームワークを用いたC++による実装は、教育的および研究的応用に適した効率的で信頼性の高いものであった。