[论文解读] 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⁵)。该方法在绝对能量误差上达到亚毫电子伏特级别，相对能量误差低于0.1 kcal/mol，与CCSD(T)相比，且可通过投影算子秩或特征值容差实现系统性误差控制。

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⁷)计算复杂度的瓶颈，该瓶颈限制了其在大分子体系中的应用。
• 通过投影算子秩或特征值容差实现系统性误差控制，同时保持体系尺寸广延性。
• 通过试点实现展示在绝对和相对能量上均达到化学精度。

提出的方法

• 该方法采用张量超收缩（THC）技术，通过CANCENCOMP/PARAFAC（CP）分解的正交投影算子对三体振幅 $ t_{ijk}^{abc} $ 进行因式分解。
• 三体振幅表示为 $ 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 该方法在多种分子体系中对绝对和相对能量的重现程度如何？
• RQ4 随着体系尺寸增大，误差如何变化？该方法是否保持亚线性至线性的误差增长？

主要发现

• THC-CCSD(T)方法通过利用三体振幅的低秩张量因式分解，将（T）能的计算复杂度从O(N⁷)降低至O(N⁵)。
• 当使用足够高的投影算子秩时，该方法在绝对能量误差上达到亚毫电子伏特（≤1 mEh）级别，相对能量误差低于0.1 kcal/mol，与CCSD(T)相比。
• 随着投影算子特征值容差减小或投影算子秩增加，能量误差平滑收敛至精确的CCSD(T)值。
• 该方法在体系尺寸增大时表现出亚线性至线性的误差增长，表明其在更大体系中具有鲁棒性和可扩展性。
• 当投影算子秩足够高时，THC-CCSD(T)生成的势能面连续且与CCSD(T)高度一致。
• 在中等投影算子秩（122–156）区间观察到误差收敛的拐点，归因于CP分解算法的数值伪影，而非方法本身。