[論文レビュー] Quantum three body problems using harmonic oscillator bases with different sizes

B. Silvestre-Brac, R. Bonnaz+2 more|arXiv (Cornell University)|Mar 24, 2020

Stellar, planetary, and galactic studies参考文献 6被引用数 6

ひとこと要約

本稿では、ジャコビ座標における異なるサイズパラメータを有する調和振動子基底関数を用いて波動関数を展開することで、量子三体問題を解くための新規手法を提案する。この手法により収束が速くなり、計算コストが低減される。行列要素はタルミの積分と最適化された数値技術を用いて正確に計算可能であり、精度を損なわずに大幅に計算効率が向上する。

ABSTRACT

We propose a new treatment for the quantum three-body problem. It is based on an expansion of the wave function on harmonic oscillator functions with different sizes in the Jacobi coordinates. The matrix elements of the Hamiltonian can be calculated without any approximation and the precision is restricted only by the dimension of the basis. This method can be applied whatever the system under consideration. In some cases, the convergence property is greatly improved in this new scheme as compared to the old traditional method. Some numerical tricks to reduce computer time are also presented.

研究の動機と目的

• 質量が不等な三体量子系において、従来の調和振動子基底法の収束が遅い問題を解決する。
• ジャコビ座標における異なる長さスケールを許容することで、計算効率を向上させ、物理的な非対称性をよりよく反映する。
• 解析的タルミの積分と数値最適化を用いて、さまざまなオシレーターサイズにおいても正確な行列要素計算を可能にする。
• メモリ使用量と計算時間を削減しながらも、特に粒子質量の差が大きい系においても高い精度を維持する。
• 非相対論的および半相対論的運動論に普遍的に適用可能な体系的で包括的なフレームワークを提供する。

提案手法

• 三体波動関数を、2つのジャコビ座標における異なるサイズパラメータ（b_x, b_y）を有する調和振動子関数で展開する。
• 運動エネルギーおよびポテンシャルエネルギーの行列要素に対して、ポテンシャルVおよびオシレーターサイズbに依存するタルミの積分I_p(V,b)を用いた正確な解析的表現を用いる。
• 係数B(n,n',l,p)を事前に計算し、連続的なメモリ配置に格納することで、効率的な順次アクセスを実現し、行列要素の計算を高速化する。
• 試行関数（例：y = αx⁻² + βx¹ + γx⁰）を用いた修正ブレント法を適用し、非線形パラメータb_xおよびb_yに関するエネルギー最小化を実行する。
• 運動量空間表現を用いることで、半相対論的運動論を自然に扱える。なぜなら、HO関数のフーリエ変換もHO関数のまま保たれるからである。
• 再帰関係および漸近展開を用いて、タルミの積分に含まれる特殊関数（例：クンマーのU）の正確な評価を実現する。

実験結果

リサーチクエスチョン

• RQ1 質量が不等な三体量子系は、ジャコビ座標における異なるサイズパラメータを有する調和振動子基底関数を用いることで、より効率的に解けるか？
• RQ2 本手法と従来の固定サイズの調和振動子基底法との間で、エネルギー固有値の収束速度にはどのような差があるか？
• RQ3 異なるオシレーターサイズを用いることで、与えられた精度における計算コストおよびメモリ要件にどのような影響があるか？
• RQ4 非局所的または半相対論的ポテンシャルの行列要素は、サイズが変化する基底関数においても正確かつ効率的に計算可能か？
• RQ5 修正ブレント法などの数値最適化技術は、サイズパラメータに関するエネルギー最小化をどの程度高速化できるか？

主な発見

• 本手法は、質量差が大きい系において特に顕著な収束速度の向上を達成する。
• ハミルトニアンの行列要素は、タルミの積分を用いて近似なしに正確に計算可能であり、精度は基底の次元にのみ制限される。
• 事前に計算した係数B(n,n',l,p)と連続的なメモリアクセスにより、和分演算の最適化が図られ、計算時間が著しく短縮される。
• 試行関数（例：x⁻² + x¹ + 定数）を用いた修正ブレント法により、サイズパラメータに関するエネルギー最小化がより高速かつ安定して実行可能となる。
• 本手法は非相対論的および半相対論的運動論の両方に対して普遍的に適用可能であり、運動量空間における行列要素は解析的に計算可能である。
• 与えられた精度において、必要な基底状態の数は劇的に削減される（例：N=1000からN=10に）。これにより、観測量の評価における計算コストが100倍低減される。