[论文解读] Resolution of the identity approximation applied to PNOF correlation calculations

Juan Felipe Huan Lew-Yee, Mario Piris+1 more|arXiv (Cornell University)|Dec 31, 2020

Advanced Chemical Physics Studies参考文献 78被引用 20

一句话总结

本文提出 DoNOF-RI，一种应用于 PNOF7 电子相关方法的解析-因子（RI）近似，将算术复杂度从五次方降至四次方，内存复杂度从四次方降至三次方，从而实现显著的计算加速。RI 收敛后，通过使用四中心电子排斥积分重启，可恢复精确结果，该方法在环烷烃及多种分子上均验证了高精度与高效率。

ABSTRACT

In this work, the required algebra to employ the resolution of the identity approximation within Piris Natural Orbital Functional (PNOF) is developed, leading to an implementation named DoNOF-RI. The arithmetic scaling is reduced from fifth-order to fourth-order, and the memory scaling is reduced from fourth-order to third-order, allowing significant computational time savings. After the DoNOF-RI calculation has fully converged, a restart with four-center electron repulsion integrals can be performed to remove the effect of the auxiliary basis set incompleteness, quickly converging to the exact result. The proposed approach has been tested on cycloalkanes and other molecules of general interest to study the numerical results as well as the speed-ups achieved by PNOF7-RI when compared with PNOF7.

研究动机与目标

• 为降低 PNOF7 电子相关计算的高成本，其复杂度因四中心电子排斥积分变换而达到 O(N^5)。
• 在 PNOF 框架内实现解析-因子（RI）近似，以降低算术与内存复杂度。
• 通过 RI 实现更快收敛，同时通过使用精确四中心积分的后处理重启保留精度。
• 在环烷烃及其他分子上验证该方法，评估其效率提升与能量精度。

提出的方法

• 应用 RI 近似，以辅助基集展开替代四中心电子排斥积分（ERI），降低积分变换成本。
• 采用密度拟合方法，将 AO-ERI 表示为辅助基函数的线性组合，从而高效计算库仑与交换矩阵。
• 在轨道与占据数优化步骤中均实现 RI，将 AO 到 MO 的 ERI 变换复杂度从 O(N^5) 降低至 O(N^4)。
• 在 PNOF7 中引入电子对约束，进一步缩小活跃轨道空间，提升效率。
• 在 RI 收敛后，使用完整的四中心 ERI 进行后处理重启，以消除辅助基集不完整的影响，恢复精确结果。
• 在优化过程中采用有界 Bunch-Kaufman 对角主元分解方法进行对称矩阵分解。

实验结果

研究问题

• RQ1 RI 近似能否成功应用于 PNOF7 电子相关计算，以降低计算复杂度？
• RQ2 在 PNOF7 中应用 RI 后，算术与内存复杂度的降低程度如何？
• RQ3 RI 近似是否引入显著的关联能误差，这些误差是否可被修正？
• RQ4 在环烷烃及其他测试体系中，实际计算速度提升了多少？
• RQ5 使用完整四中心积分的后处理 RI 重启能否高效恢复精确的 PNOF7 结果？

主要发现

• DoNOF-RI 方法将算术复杂度从 O(N^5) 降低至 O(N^4)，内存复杂度从 O(N^4) 降低至 O(N^3)，实现了显著的计算节省。
• 该方法在 PNOF7 计算中实现了显著加速，尤其在较大体系中更为明显，归因于积分变换复杂度的降低。
• 在 RI 完全收敛后，通过使用四中心 ERI 重启，可快速恢复精确的 PNOF7 结果，证实了 RI 近似的准确性。
• 环烷烃（最多九个碳原子）及其他分子（如苯、恶唑、硼嗪）的能量结果与参考 PNOF7 值高度一致。
• 该方法在大幅提高计算效率的同时保持了化学精度，尤其适用于强关联体系。
• 基于电子对的 PNOF7 实现结合 RI，可在显著降低计算成本的同时实现对强关联的精确处理。