[論文レビュー] Observations on Gaussian bases for Schrodinger's equation

Ilan Degani|ArXiv.org|Jul 31, 2007

Spectroscopy and Quantum Chemical Studies参考文献 19被引用数 3

ひとこと要約

本稿では、急勾配のポテンシャル領域における過剰に広がった希薄なガウス関数を避けることで、多次元シュレーディンガー固有値問題を解くために修正されたガウス基底法を提案する。特に高いポテンシャル障壁付近で精度が向上する。この手法はコロケーション法およびガレルキン法の両方を強化し、数値結果では、直交性のずれよりもハミルトニアンのトレース最小化基準を用いた幅スケーリングを行うことで、著しい誤差低減が得られている。

ABSTRACT

One of the few methods for generating efficient function spaces for multi-D Schrodinger eigenproblems is given by Garashchuk and Light in J.Chem.Phys. 114 (2001) 3929. Their Gaussian basis functions are wider and sparser in high potential regions, and narrower and denser in low ones. We suggest a modification of their approach based on the following observation: In very steep potential regions, wide, sparse, Gaussians should be avoided even if their centers have high potential values. Our numerical results illustrate that a dramatic improvement in accuracy may be obtained in this way. We also compare the errors of collocation to those of a Galerkin approach, test a criterion for scaling Gaussian widths based on deviation from orthogonality of collocation eigenfunctions, and suggest a criterion for scaling Gaussian widths based on Hamiltonian trace minimization.

研究の動機と目的

• 急勾配のポテンシャル領域における既存のガウス基底法の非効率性（広がった希薄なガウス関数による精度の低下）を是正すること。
• 多変数シュレーディンガー固有値問題におけるコロケーション法およびガレルキン法の両方の精度を向上させること。
• ガウス関数の幅スケーリングのための新しい基準を開発・検証すること、特にハミルトニアンのトレース最小化に基づくもの。
• 各手法に最適化された基底関数を用いた場合の、コロケーション法とガレルキン法の相対的性能を評価すること。
• 固有関数の直交性のずれが、コロケーション法における幅選択の指針として有効であるかを評価すること。

提案手法

• ガラシチュクとライトのガウス基底構築法を修正し、高ポテンシャルで勾配が急な領域における広がった希薄なガウス関数の使用を制限する。
• 適応的スケーリングされたガウス基底を用いて、1次元のモースポテンシャル問題にガレルキン法およびコロケーション法を適用する。
• 基底関数の広がりを制御するグローバルな幅パラメータ c を導入し、平均固有関数誤差を最小化するように最適化する。
• ハミルトニアン行列のトレースを最小化する基準に基づく、新たな幅スケーリング基準を導入し、安定性と精度を向上させる。
• コロケーション法で得られた固有関数の直交性のずれをヒューリスティックとして用い、幅選択の指針とする。
• 重なり行列 S と基底関数行列 Φ の条件数を比較し、数値的安定性を評価する。

実験結果

リサーチクエスチョン

• RQ1 急勾配のポテンシャル領域における広がった希薄なガウス関数の回避が、多変数シュレーディンガー問題におけるガウス基底法の精度を顕著に向上させ得るか？
• RQ2 各手法に別個に最適化された基底関数を用いた場合、コロケーション法とガレルキン法の性能はどのように異なるか？
• RQ3 コロケーション法で得られた固有関数の直交性のずれが、最適なガウス幅選択の信頼できる指標として機能するか？
• RQ4 ハミルトニアンのトレース最小化に基づく幅スケーリング基準が、従来の方法よりも優れた精度をもたらすか？
• RQ5 基底関数の分布および幅が、コロケーション法およびガレルキン法におけるシステム行列の条件数に与える影響は何か？

主な発見

• 急勾配領域における広がったガウス関数の回避を実施した修正ガウス基底法は、元のガラシチュク＝ライト法と比較して、著しい精度向上を達成している。
• 各手法に別個に最適化された幅を用いた場合、ガレルキン法は固有値および固有関数誤差の両面でコロケーション法を常に上回り、ガレルキン／基底集合3の組み合わせが最小の誤差を達成した。
• ハミルトニアンのトレース最小化基準に基づく幅スケーリングは、直交性のずれヒューリスティックよりも優れた精度とより良好な条件数を持つ行列をもたらした。
• モースAポテンシャルにおいて、ガレルキン法と基底集合3を用いた場合、最適幅で平均固有関数誤差は10^-10未満に低下し、SおよびΦの条件数は10^10未満にまで低下した。
• コロケーション法の性能は、この手法に最適化された幅を用いることで顕著に向上したが、依然としてガレルキン法に劣り、特に高精度領域では顕著であった。
• トレース最小化基準を用いることで、重なり行列Sの条件数は桁単位で低下し、数値的安定性の向上が示された。