[论文解读] Efficiency of a Multi-Reference Coupled Cluster method

Emmanuel Giner, Grégoire David+2 more|arXiv (Cornell University)|Sep 10, 2015

Advanced Chemical Physics Studies参考文献 58被引用 7

一句话总结

该论文提出了一种状态特定的、未收缩的多参考耦合簇（MRCC）方法，通过使用 ${\hat{T}}^2$ 算符迭代地对 MRSDCI 矩阵进行修正，以估算三激发和四激发贡献，从而在精度上超越传统 CASSDCI 方法。该方法在键断裂和多组态体系中实现了化学精度（误差约 1–2 mH），并具备完全的并行性以及在参考态和激发态选择上的灵活性。

ABSTRACT

The multi-reference Coupled Cluster method first proposed by Meller et al (J. Chem. Phys. 1996) has been implemented and tested. Guess values of the amplitudes of the single and double excitations (the ${\hat T}$ operator) on the top of the references are extracted from the knowledge of the coefficients of the Multi Reference Singles and Doubles Configuration Interaction (MRSDCI) matrix. The multiple parentage problem is solved by scaling these amplitudes on the interaction between the references and the Singles and Doubles. Then one proceeds to a dressing of the MRSDCI matrix under the effect of the Triples and Quadruples, the coefficients of which are estimated from the action of ${\hat T}^2$. This dressing follows the logics of the intermediate effective Hamiltonian formalism. The dressed MRSDCI matrix is diagonalized and the process is iterated to convergence. The method is tested on a series of benchmark systems from Complete Active Spaces (CAS) involving 2 or 4 active electrons up to bond breakings. The comparison with Full Configuration Interaction (FCI) results shows that the errors are of the order of a few milli-hartree, five times smaller than those of the CASSDCI. The method is totally uncontracted, parallelizable, and extremely flexible since it may be applied to selected MR and/or selected SDCI. Some potential generalizations are briefly discussed.

研究动机与目标

• 开发一种针对强关联体系的、大小一致、大小广延且灵活的多参考耦合簇方法。
• 通过基于参考态与单激发、双激发之间相互作用强度的幅度缩放，解决 MRCC 中的多重母体问题。
• 通过迭代地对 CI 矩阵进行修正，利用 $\hat{T}^2$ 估算的三激发和四激发贡献，提升 CASSDCI 的精度。
• 通过状态特定的形式化方法，实现对激发态和选定参考态空间的应用。
• 提供一种计算效率高、未收缩且可并行化的框架，用于高精度量子化学计算。

提出的方法

• 该方法首先从 MRSDCI 系数中提取单激发和双激发的初始幅度（$\hat{T}$）。
• 通过参考态与单激发、双激发之间相互作用强度的缩放，解决多重母体问题。
• 通过 $\hat{T}^2$ 算符的作用，估算三激发和四激发的贡献，对 MRSDCI 矩阵进行迭代修正。
• 修正过程遵循中间有效哈密顿量形式化，保持大小一致性并改善相关能估计。
• 在每次迭代中对修正后的矩阵进行对角化，当幅度稳定时过程收敛。
• 该方法为完全未收缩形式，支持灵活选择参考态和激发态，且在最耗时的步骤（三激发/四激发计算）中天然具备并行性。

实验结果

研究问题

• RQ1 状态特定的、未收缩的 MRCC 方法是否能在多组态体系中实现比 CASSDCI 更高的精度？
• RQ2 如何通过基于物理动机的幅度缩放，解决 MRCC 中的多重母体问题？
• RQ3 通过 $\hat{T}^2$ 估算的三激发和四激发对 CI 矩阵进行迭代修正，能在多大程度上改善相关能估计？
• RQ4 该形式化方法是否可扩展至激发态，同时保持大小一致性和精度？
• RQ5 在应用选择性参考态和激发态空间截断时，计算权衡和性能提升如何？

主要发现

• 与全配置相互作用（Full CI）相比，该方法的误差约为 1–2 毫哈特ree，比 CASSDCI 的误差小五倍。
• 在一系列基准体系中，包括键断裂过程和具有 2 或 4 个活性电子的体系中，该方法保持了高精度。
• 该方法为完全未收缩且可并行化，最耗时的步骤（三激发/四激发）可进行分布式计算。
• 基于 $\hat{T}^2$ 的迭代修正过程有效捕捉了高阶相关效应，而无需显式引入三激发或四激发幅度。
• 该方法具有灵活性，可应用于选定的参考态空间和/或选定的 SDCI 构型，从而实现计算成本的降低。
• 通过状态特定或状态平均的程序，该形式化方法可扩展至激发态，且计算后可对本征求解矢量进行正交化处理。