结合Cholesky分解与点群对称性实现CCSD解析梯度的高效新算法及其在分子几何优化中的应用

《The Journal of Physical Chemistry A》：A Novel Implementation of CCSD Analytic Gradients Using Cholesky Decomposition of the Two-Electron Integrals and Abelian Point-Group Symmetry

本文介绍了一种基于电子排斥积分（ERI）Cholesky分解（CD）与阿贝尔点群对称性相结合的耦合簇单双（CCSD）解析梯度创新实现方案。该方案针对中等及大型对称分子体系，通过对称适配的两步算法进行双电子积分CD，并在求解CCSD密度矩阵时对对称适配的核位移导数进行Cholesky向量微分。文档通过评估包含数百个基函数的对称体系几何优化效率，定量展示了利用点群对称性带来的计算增益，为精确预测分子性质提供了高效工具。

　　

2. 理论

文章的理论部分详述了Cholesky分解（CD）在量子化学中的应用，特别是针对双电子积分及其导数的处理，并概要推导了CCSD解析梯度相关的表达式。

2.1. 双电子积分的Cholesky分解

由于电子排斥积分（ERI）矩阵（采用Mulliken表示法）是对称且半正定的，因此可以通过Cholesky分解来表示。该分解利用了ERI矩阵的秩亏缺特性，用Cholesky向量（CV）表示，从而有效压缩了全张量中存储的信息。由于ERI张量本身不是满秩的，其非零特征值的数量仅随基组大小线性增长，因此实际需要计算的Cholesky向量数量N_ch本身也按O(N)标度，远小于所有唯一基对的数量，这降低了内存需求。

对于ERI导数的处理，由于导数矩阵不再是半正定的，无法直接进行CD。文中介绍了两种策略：一是对未微分的ERI的CD表达式直接求导，但这需要处理依赖于所有扰动的微分Cholesky向量，在几何扰动情况下较为繁琐；二是利用CD与分辨率恒等式/密度拟合（RI/DF）近似之间的形式等价性，通过引入辅助基组并利用度量矩阵的Cholesky分解，可以避免构造微分Cholesky向量。后一种方法通过定义变换后的Cholesky向量，得到了ERI张量一阶导数的CD型表达式，为后续梯度计算提供了便利。

2.2. CCSD的拉格朗日量

由于耦合簇（CC）方法的非变分性质，直接对电子能量求导会增加计算的复杂性。为了消除求解扰动振幅方程的需要，文章采用了拉格朗日方法。通过定义一个CC拉格朗日量，其中引入Λ算符（包含λ振幅）作为拉格朗日乘子来强制t振幅满足通常的CC振幅方程，并引入其他乘子来强制执行Brillouin条件和分子轨道的正交归一性。对于CCSD，Λ算符展开为Λ?₁+ Λ?₂。该方法简化了梯度表达式的推导。

2.3. 关于簇振幅的平稳性：Λ方程

要求CC拉格朗日量对t振幅和λ振幅都是平稳的。对λ振幅的平稳性条件等价于通常的t振幅方程。而对t振幅的平稳性条件则产生一组关于λ振幅的线性方程，即Λ方程。文中给出了针对闭壳层（RHF参考波函数）情况、完全自旋适应的单重和双重λ振幅的具体方程，其实现基于Gauss, Stanton和Bartlett报道的方程。这些方程涉及Fock矩阵元、双电子积分变换后的中间体（如W?, W）以及由CCSD密度矩阵和Cholesky向量构成的中间体G等，需要迭代求解以获得λ振幅。

总结

总的来说，本文提出了一种新颖的CCSD解析梯度实现方案，其核心是结合了双电子积分的Cholesky分解近似与阿贝尔点群对称性的利用。理论部分系统阐述了如何通过CD高效表示ERI及其导数，并采用拉格朗日框架推导CCSD梯度，避免了直接求解扰动方程。通过求解Λ方程获得λ振幅，最终能有效计算能量对核坐标的导数（即几何梯度）。该方案特别适用于中等大小的对称分子体系，能在保持CCSD方法高精度的同时，显著降低计算成本和内存需求，为精确的分子几何优化和性质预测提供了有力的计算工具。