[論文レビュー] Algorithm for reduction of boundary-value problems in multistep adiabatic approximation

А. А. Гусев, О. Чулуунбаатар+4 more|arXiv (Cornell University)|May 12, 2010

Differential Equations and Boundary Problems参考文献 17被引用数 5

ひとこと要約

本稿では、マルチステップ一般化カントロヴィチ法に基づく多段断熱近似を用いて、多次元境界値問題を低減するための記号的数値的アルゴリズム、MSGKM を提示する。この手法により、ヘリウム原子を例に挙げて、高速変数を段階的に除去することで、少数粒子量子系における基底状態および励起状態の効率的計算が可能となり、パラメトリック固有値問題とMPI並列化された連立常微分方程式の解法によって収束性が向上し、計算コストが低減される。

ABSTRACT

The adiabatic approximation is well-known method for effective study of few-body systems in molecular, atomic and nuclear physics, using the idea of separation of "fast" and "slow" variables. The generalization of the standard adiabatic ansatz for the case of multi-channel wave function when all variables treated dynamically is presented. For this reason we are introducing the step-by-step averaging methods in order to eliminate consequently from faster to slower variables. We present a symbolic-numerical algorithm for reduction of multistep adiabatic equations, corresponding to the MultiStep Generalization of Kantorovich Method, for solving multidimensional boundary-value problems by finite element method. An application of the algorithm to calculation of the ground and first exited states of a Helium atom is given.

研究の動機と目的

• すべての変数を動的に扱う多チャンネル波動関数への標準断熱アンザッツの一般化を目的とする。
• 少数粒子量子系に生じる多次元固有値問題を解くための体系的な記号的数値的アルゴリズムの開発を目的とする。
• パラメトリック基底関数を用いた高速変数の段階的除去によって、高次元問題における収束性の向上と計算コストの低減を目的とする。
• ヘリウム原子、量子ドット、複雑なポテンシャルを有する量子井戸などの原子およびナノスケール系への応用を可能にするための目的とする。
• 将来のMPIおよびGRID環境への実装を支援し、大規模パラメトリック固有値問題のスケーラブルな解法を実現するための目的とする。

提案手法

• 本手法は、全波動関数を高速な変数から始める階層的なパラメトリック解に分解する多段断熱アンザッツを採用する。
• 各ステップに対して、逐次的に遅い変数に依存する補助的パラメトリック固有値問題を定式化する。
• アルゴリズムは、パラメトリック依存性を解析的に制御できるように、Mapleを用いた記号計算により有効ポテンシャルおよび行列要素を生成する。
• 有効ポテンシャルの行列要素は、より速い変数に対するネストされた積分を再帰的に計算することで得られ、直前のステップで得られた改良された基底関数が用いられる。
• 境界条件および正規直交制約が各段階で強制され、解の整合性と数値的安定性が維持される。
• MPIを用いた並列化が設計されており、結果として得られる2N−1個の結合された常微分方程式のスケーラブルな解法が可能となる。

実験結果

リサーチクエスチョン

• RQ1 標準断熱近似は、すべての変数を動的に扱う多チャンネル波動関数にどのように一般化できるか？
• RQ2 多次元量子系において収束速度を最大化する最適な変数除去の順序は何か？
• RQ3 パラメトリック固有値問題を体系的に導出し、有効ポテンシャルの正確かつ安定した計算を保証するにはどうすればよいか？
• RQ4 MPIまたはGRID技術を用いて、大規模量子系に対して本手法を効率的に並列化できるか？
• RQ5 ヘリウム原子に対して、本多段アプローチは、従来の断熱法や基底関数展開法と比較して、精度および収束性において優れているか？

主な発見

• MSGKMアルゴリズムは、元の多次元境界値問題を、2N−1個のパラメトリック固有値問題の列に成功して低減した。各問題は独立に解ける。
• 特異点（三重衝突点 R→0）に近い領域でも、有効ポテンシャルおよびエネルギー曲線が滑らかに振る舞い、信頼性の高い収束性を保証する。
• 展開に二パラメトリック基底関数を用いることで、ルジャンドル多項式展開と比較して、行列要素計算の収束速度と精度が向上した。
• 計算されたエネルギー固有値は上界にあり、特に R→0 および R→∞ の漸近的極限において、既知の物理的期待と良好に一致している。
• 三体閾値以下の散乱問題に必要な漸近状態と整合性があることが、フォック展開および先行研究との比較により確認された。
• 本アルゴリズムは、変数係数および外部場を有する複雑な量子系を解くためのフル機能プログラム（KANTBP 3.0）への拡張が進行中である。