[論文レビュー] Accurate calculation of eigenvalues and eigenfunctions. I: Symmetric potentials

Francisco M. Fernández|ArXiv.org|Jul 3, 2008

Quantum, superfluid, helium dynamics参考文献 18被引用数 6

ひとこと要約

本稿では、調和振動子や二重井戸型振動子などの対称的ポテンシャルを有する1次元量子系の固有値および固有関数を高精度に計算するためのリッカティ・パデ法（Riccati-Padé Method, RPM）を紹介する。波動関数の対数微分をテイラー級数に展開し、得られる漸化式を解くことで、エネルギー準位および期待値に対する収束が速い上下界が得られ、非調和振動子や二重井戸型ポテンシャルにおいて20桁以上の精度を達成する。

ABSTRACT

We describe a method for the calculation of accurate energy eigenvalues and expectation values of observables of separable quantum-mechanical models. We discuss the application of the approach to one-dimensional anharmonic oscillators with symmetric potential-energy functions.

研究の動機と目的

• 分離可能な量子力学的系において、対称的ポテンシャルを有する系の高精度なエネルギー固有値および期待値を計算する体系的な手法を開発すること。
• 行列要素の計算を要しない、従来の摂動法や行列法の代替手法を提供すること。
• 記号計算を用いて、非調和振動子や二重井戸型振動子に対して、RPMの収束速度の高さと信頼性を示すこと。
• 計算代数系（CAS）を用いた教育的価値を強調すること。

提案手法

• シュレーディンガー方程式から導かれるリッカティ方程式を用い、中心変数として $ f(x) = \frac{s}{x} - \frac{\psi'(x)}{\psi(x)} $ を採用する。
• 変数 $ z = x^2 $ を用いて $ f(x) $ を $ x=0 $ の周りでテイラー級数展開すると、係数 $ f_j(E) $ に対する漸化式が得られる。
• 漸化式の係数 $ f_j(E) $ はエネルギー $ E $ に関する多項式であり、これによりパデ近似を構築し、固有値に対応する根を特定できる。
• 固有値は、$ f_j $ の係数から構成されるヘンケル行列式の根として求められ、上界と下界 $ E_n^{[D,0]} < E_n < E_n^{[D,1]} $ が得られる。
• 収束性は $ \beta $ を変化させながら、$ V(x) = x^4 + \beta x^2 $ のような対称的ポテンシャルに対して検討される。
• 期待値 $ \langle x^2 \rangle $ は系列係数から抽出可能であり、基底状態では $ \langle x^2 \rangle = e^{(1)} $ と表される。

実験結果

リサーチクエスチョン

• RQ1 リッカティ・パデ法は、行列要素の計算を伴わず、対称的量子ポテンシャルの固有値を高精度に提供できるか？
• RQ2 非調和振動子に対して、RPMの収束速度は、モーメント法などの従来手法と比較してどのように異なるか？
• RQ3 障壁の深さを変化させた二重井戸型ポテンシャルに適用した際、この手法の挙動はいかなるものか？
• RQ4 RPMではなぜ偽の根や複数の根が現れるのか、それらはどのように特定・除外できるか？
• RQ5 この手法を用いて、系列展開から直接期待値を計算できる範囲はどの程度か？

主な発見

• 非調和振動子 $ V(x) = x^4 + \beta x^2 $ に対して、$ \beta = 0 $ の場合、RPMは固有値を20桁の精度で達成し、$ \langle x^2 \rangle = 0.3620226487886768452 $ を得る。
• $ \beta = -15 $ の場合、基底状態エネルギーは $ -50.8413872843819547 < E_0 < -50.8413872843819543 $ と境界づけられ、励起状態に対しても同様の精度を示す。
• RPMはエネルギー準位に対して厳密な上界と下界を提供し、$ E_n^{[D,0]} < E_n < E_n^{[D,1]} $ を満たし、次数 $ D $ の増加に伴い収束性が向上する。
• RPMの収束速度は、特に深い二重井戸型ポテンシャルに対して、モーメント法よりも顕著に速い。
• 根の特定プロセスにおいて偽の根や複数の根が現れるが、物理的に意味のある根として正しい固有値が一貫して特定可能である。
• 期待値 $ \langle x^2 \rangle $ は系列係数から直接得られ、基底状態では $ \langle x^2 \rangle = e^{(1)} $ と表される。