[論文レビュー] Using tensor hypercontraction density fitting to achieve an O(L^4) CISD algorithm

Neil Shenvi, Helen van Aggelen+1 more|arXiv (Cornell University)|Sep 13, 2012

Tensor decomposition and applications被引用数 3

ひとこと要約

本稿では、テンソルハイパコントラクション密度適合（THC-DF）を用いた、単一および二重励起を伴う配置相互作用（CISD）のO(L⁴)アルゴリズムを提示する。これにより、標準的なCISDにおけるO(L⁶)の計算スケーリングが低減される。O(L)の補助関数を用いた低ランクテンソル因子分解により電子反発積分を分解することで、計算コストを著しく削減しつつ、近似的にCISDに近い精度を達成し、正確なCISD極限への形式的接近を保ったままとなる。

ABSTRACT

Recently, Hohenstein et al[1] introduced tensor hypercontraction density fitting to decompose the rank-4 electron repulsion integral tensor as the product of five rank-2 tensors. In this paper, we use this methodology to construct an algorithm which calculates the approximate ground state energy in O(L^4) operations. We test our method using several small molecules and show that we quickly approach the CISD limit with a small number of auxiliary functions.

研究の動機と目的

• 従来のCISD手法におけるO(L⁶)スケーリングを克服し、小分子にしか適用できない制限を解消すること。
• テンソルハイパコンタクション密度適合（THC-DF）を活用することで、形式的にO(L⁴)スケーリングを有する波動関数ベースの手法を開発し、高い精度を維持すること。
• 補助関数の数を増やすに従って、この手法が正確なCISDエネルギーに収束することを示し、CISDの既知の限界にもかかわらず物理的挙動を保つこと。
• 従来の相関手法が計算的に不可能とされる系において、効率的かつ高精度な電子構造計算を可能にすること。

提案手法

• 本手法は、O(L)の補助関数を用いて、ランク4の電子反発積分（ERI）テンソルを5つのランク2テンソルの積に分解するテンソルハイパコンタクション密度適合（THC-DF）を採用する。
• ERIテンソルはε^{ik}_{jl} = Σ_{a,b} x_{ia} x_{ka} Z_{ab} x_{jb} x_{lb}として因子分解され、基底間の変換や効率的な計算が可能になる。
• CISDエネルギーは、励起振幅を用いた波動関数のパrametrizationに基づき計算され、ERIテンソルはTHC-DFを介して表現されることで計算コストが低減される。
• THCテンソルZ_{ab}に関するエネルギー勾配を解析的に導出し、補助基底パラメータの反復最適化が可能になる。
• 本手法は、ERIが既に正規直交分子軌道基底に変換されていると仮定しており、THC-DFを用いたERIテンソルのO(L³)変換が可能になる。
• 繰り返しのテンソル畳み込み結果をキャッシュすることで性能を向上させ、メモリを消費して計算速度を向上させる。

実験結果

リサーチクエスチョン

• RQ1 テンソルハイパコンタクション密度適合（THC-DF）を用いることで、CISDのスケーリングをO(L⁶)からO(L⁴)に低減しつつ、精度を保持できるか？
• RQ2 補助関数の数を増やすに従って、THC-DFに基づくCISD手法が正確なCISDエネルギーにどの程度速やかに収束するか？
• RQ3 標準的なCISDや他の相関手法と比較して、このO(L⁴) CISD手法の計算コストと精度のトレードオフはいかなるものか？
• RQ4 ERIテンソルのTHC-DF因子分解が、顕著な誤差を導入せずに、配置相互作用手法に効果的に適用可能か？

主な発見

• 提案されたアルゴリズムは、CISDエネルギー評価においてO(L⁴)スケーリングを達成し、標準的なO(L⁶)スケーリングに比べて顕著な低減を実現した。
• 補助関数の数を増やすに従って、この手法は正確なCISDエネルギーに収束することが確認され、完全なCISD極限に近づく能力が裏付けられた。
• わずかにO(L)程度の補助関数を用いるだけで、近似的にCISDに近い精度が達成され、急速な収束が示された。
• ERI因子分解に正確なTHC-DF分解（ERIテンソルの固有値分解を介して）を用いることで、精度評価時に因子分解による近似誤差が生じないことが保証された。
• O(L⁶)スケーリングのため従来は計算不可能とされていた系、特に大きな分子に対しても、CISDエネルギーの効率的計算が可能になった。