[论文解读] Understanding Discrepancies of Wavefunction Theories for Large Molecules

Tobias Schäfer, Andreas Irmler+2 more|arXiv (Cornell University)|Jul 1, 2024

Quantum optics and atomic interactions被引用 10

一句话总结

本文指出，由于多体微扰级数截断，CCSD(T) 理论对大型非共价分子复合物存在过度成键问题，并提出 CCSD(cT) 作为校正方法，使该方法在最多 132 个原子的体系中与扩散量子蒙特卡罗 (DMC) 方法达到化学精度（误差在 1 kcal/mol 以内）。通过采用平面波基组和系统化的相关处理，作者解决了 CCSD(T) 与 DMC 之间长期存在的分歧，为大规模材料和药物设计应用提供了可靠的基准预测方法。

ABSTRACT

Quantum mechanical many-electron calculations can predict properties of atoms, molecules and even complex materials. The employed computational methods play a quintessential role in many scientifically and technologically relevant research fields. However, a question of paramount importance is whether approximations aimed at reducing the computational complexity for solving the many-electron Schrödinger equation, are accurate enough. Here, we investigate recently reported discrepancies of noncovalent interaction energies for large molecules predicted by two of the most widely-trusted many-electron theories: diffusion quantum Monte Carlo and coupled-cluster theory. We are able to unequivocally pin down the source of the puzzling discrepancies and present modifications to widely-used coupled-cluster methods needed for more accurate noncovalent interaction energies of large molecules on the hundred-atom scale. This enhances the reliability of predictions from quantum mechanical many-electron theories across a wide range of critical applications, including drug design, catalysis, and the innovation of new functional materials, such as those for renewable energy technologies.

研究动机与目标

• 解决 CCSD(T) 与 DMC 在大型分子复合物中非共价相互作用能长期存在的分歧问题。
• 识别这些分歧的根本原因，特别是在含数百个原子的大体系中的表现。
• 开发并验证一种改进的耦合簇方法 CCSD(cT)，以校正 CCSD(T) 中的过度成键问题。
• 建立一种可靠且可系统性改进的计算框架，采用平面波基组和完整相关处理，实现基准质量的结果。
• 为药物设计、催化和可再生能源材料等大规模应用提供准确、高保真的预测能力。

提出的方法

• 采用平面波基组以消除原子中心高斯函数带来的基组不完备性和偏差。
• 使用标准 CCSD(T) 实现方式，采用完整非局域相关处理，避免先前研究中使用的局域相关近似。
• 提出 CCSD(cT)，一种改进的耦合簇方法，通过紧凑且系统化的校正调整微扰三激发项，从而校正过度成键问题。
• 应用线性回归模型（CCSD(cT)-fit）从 LNO-CCSD(T) 数据中估算相互作用能，并与全 CCSD(T) 和 DMC 结果进行验证。
• 执行大规模 DMC 计算并附带误差条，与 CCSD(T) 和 CCSD(cT) 结果进行比较，确保多种实现方式之间的一致性。
• 采用系统性误差估计和交叉验证评估可靠性，尤其针对 C60@[6]CPPA 等高复杂度体系。

实验结果

研究问题

• RQ1 为何 CCSD(T) 与 DMC 对大型分子复合物的非共价相互作用能结果显著不同，尽管两者均被视为高度精确？
• RQ2 该分歧是源于 CCSD(T) 中的近似（如微扰级数截断），还是源于 DMC 中的近似（如有限尺寸效应或基组限制）？
• RQ3 能否通过改进的 CCSD 方法——CCSD(cT)——解决大型非共价体系中 CCSD(T) 的过度成键问题？
• RQ4 CCSD(cT) 方法在多大程度上能对多样化的大分子体系实现化学精度（与 DMC 误差在 1 kcal/mol 以内）？
• RQ5 当应用于超过 100 个原子的体系（如冠状烯二聚体和 C60@[6]CPPA）时，CCSD(cT) 方法的鲁棒性如何？

主要发现

• CCSD(T) 与 DMC 在大型分子中非共价相互作用能的分歧，主要源于 CCSD(T) 中多体微扰级数的截断，导致过度成键。
• CCSD(cT) 作为 CCSD(T) 的校正版本，显著降低了过度成键误差，并在所有测试体系中（包括 132 原子的 C60@[6]CPPA 复合物）与 DMC 结果一致，误差在 1 kcal/mol 以内。
• 对于平行位移的苯二聚体，平面波基组的 CCSD(T) 结果与高斯基组基准结果高度一致，验证了计算框架的可靠性。
• 基于 LNO-CCSD(T) 数据构建的 CCSD(cT)-fit 模型表现出优异的预测能力，对 GGG、GCGC 和 C2C2PD 体系的误差均在 0.5 kcal/mol 以内。
• 在多个案例中，不同研究对 DMC 的估计值差异超过了其报告的误差条，表明大型体系中 DMC 结果存在不确定性。
• 剩余的微小分歧（如 C60@[6]CPPA 的 1.1 kcal/mol）可能源于 DMC、CCSD(cT)-fit 或 LNO-CCSD(T) 误差估计的不确定性，而非 CCSD(cT) 本身的根本缺陷。