[論文レビュー] Tensor Hypercontraction Form of the Perturbative Triples Energy in Coupled-Cluster Theory

Andy Jiang, Justin M. Turney+1 more|arXiv (Cornell University)|Oct 13, 2022

Physics of Superconductivity and Magnetism参考文献 74被引用数 7

ひとこと要約

本稿は、カップルドクラスター理論における摂動的三体項(T)エネルギーのテンソルハイパーコントラクション(THC)形式を導入し、三体振幅を直交射影子を介して因子化し、CANCENCOMP/PARAFAC(CP)分解を用いることで、計算スケーリングをO(N⁷)からO(N⁵)に低減する。この手法は、CCSD(T)と比較して絶対エネルギー誤差が1ミリハートリー未満、相対エネルギー誤差が0.1 kcal/mol未塔を達成し、射影子のランクまたは固有値許容誤差を用いた系統的収束を実現する。

ABSTRACT

We present the working equations for a reduced-scaling method of evaluating the perturbative triples (T) energy in coupled-cluster theory, through the tensor hypercontraction (THC) of the triples amplitudes ($t_{ijk}^{abc}$). Through our method we can reduce the scaling of the (T) energy from the traditional O($N^{7}$) to a more modest O($N^{5}$). We also discuss implementation details to aid future research, development, and software realization of this method. Additionally, we show that this method yields sub-millihartree (mEh) differences from CCSD(T) when evaluating absolute energies, and sub-0.1 kcal/mol energy differences when evaluating relative energies. Finally, we demonstrate that this method converges to the true CCSD(T) energy through the systematic increasing of the rank or eigenvalue tolerance of the orthogonal projector, as well as exhibiting sub-linear to linear error growth with respect to system size.

研究の動機と目的

• カップルドクラスター理論における摂動的三体項(T)エネルギーを計算する低スケーリング手法を開発すること。
• 大規模分子系への応用を制限する従来のCCSD(T)のO(N⁷)スケーリングのボトルネックを解消すること。
• 射影子のランクまたは固有値許容誤差を用いた系統的誤差制御を可能にするとともに、サイズ拡張性を維持すること。
• パイロット実装を用いて、絶対エネルギーおよび相対エネルギーの両方において化学的精度を示すことを目的とする。

提案手法

• 本手法は、CANCENCOMP/PARAFAC(CP)分解を介した直交射影子を用いて、三体振幅$ t_{ijk}^{abc} $をテンソルハイパーコントラクション(THC)で因子化する。
• 三体振幅は$ t_{ijk}^{abc} = \tilde{T}^{XYZ} y_i^X y_j^Y y_k^Z y_a^X y_b^Y y_c^Z $として表現され、ここで$ y $-テンソルと$ \tilde{T} $は因子化された形を表す。
• 直交射影子$ U^V_{ia} $は$ U^V_{ia} = \tau_{VX} y_i^X y_a^X $として分解され、振幅の低ランク表現が可能になる。
• 本手法は、射影子ランク$ n_{\text{proj}} $が系のサイズに線形にスケーリングするランク低減戦略を用い、O(N⁵)スケーリングを実現する。
• 誤差制御は、直交射影子の固有値許容誤差を調整するか、射影子基底関数の数を増やすことで達成される。
• 本手法はサイズ拡張性を維持し、十分な射影子ランクが使用されれば、連続的なポテンシャルエネルギー面がCCSD(T)と一致する。

実験結果

リサーチクエスチョン

• RQ1 CCSD(T)における摂動的三体項(T)エネルギーのO(N⁷)スケーリングを、化学的精度を保持したままO(N⁵)に低減できるか？
• RQ2 THCに基づく三体振幅の因子化は、射影子ランクまたは固有値許容誤差による系統的誤差制御を可能にするか？
• RQ3 本手法は、多様な分子系において、絶対エネルギーおよび相対エネルギーの両方でCCSD(T)エネルギーをどの程度正確に再現できるか？
• RQ4 系のサイズが増加するに従い、誤差はどのように変化するか？また、本手法は誤差の非線形的から線形的成長を維持するか？

主な発見

• THC-CCSD(T)手法は、三体振幅の低ランクテンソル因子化を活用することで、(T)エネルギーの計算スケーリングをO(N⁷)からO(N⁵)に低減する。
• 十分な射影子ランクが使用された場合、絶対エネルギー誤差が1ミリハートリー未塔（≤1 mEh）であり、相対エネルギー誤差が0.1 kcal/mol未塔を達成する。
• 射影子の固有値許容誤差を小さくするか、射影子ランクを増加させることで、誤差は正確なCCSD(T)値へ滑らかに収束する。
• 本手法は、系のサイズに対して非線形的から線形的の誤差増加を示し、大規模系へのロバストでスケーラブルな性能を示す。
• 十分に高い射影子ランクが使用された場合、THC-CCSD(T)が生成するポテンシャルエネルギー面は連続的であり、CCSD(T)と密接に一致する。
• 中間の射影子ランク（122–156）で誤差収束にひずみが生じる現象が観察され、これはCP分解アルゴリズムのアーティファクトに起因するものであり、本手法自体の問題ではないとされる。