To follow up on the unexpectedly-good performance of several coupled-cluster models with approximate inclusion of 3-body clusters [J. Chem. Phys. 151, 064102 (2019)] we performed a more complete assessment of the 3CC method [J. Chem. Phys. 125, 204105 (2006)] for accurate computational thermochemistry in the standard HEAT framework. New spin-integrated implementation of the 3CC method applicable to closed- and open-shell systems utilizes a new automated toolchain for derivation, optimization, and evaluation of operator algebra in many-body electronic structure. We found that with a double-zeta basis set the 3CC correlation energies and their atomization energy contributions are almost always more accurate (with respect to the CCSDTQ reference) than the CCSDT model as well as the standard CCSD(T) model. The mean absolute errors in cc-pVDZ {3CC, CCSDT, and CCSD(T)} electronic (per valence electron) and atomization energies relative to the CCSDTQ reference for the HEAT dataset [J. Chem. Phys. 121, 11599 (2004)], were {24, 70, 122} $μE_h/e$ and {0.46, 2.00, 2.58} kJ/mol, respectively. The mean absolute errors in the complete-basis-set limit {3CC, CCSDT, and CCSD(T)} atomization energies relative to the HEAT model reference, were {0.52, 2.00, and 1.07} kJ/mol, The significant and systematic reduction of the error by the 3CC method and its lower cost than CCSDT suggests it as a viable candidate for post-CCSD(T) thermochemistry applications, as well as the preferred alternative to CCSDT in general.
- 论文ID: 2407.08859
- 标题: "Best" iterative coupled-cluster triples model: More evidence for 3CC
- 作者: Nakul K. Teke, Ajay Melekamburath, Bimal Gaudel, Edward F. Valeev
- 机构: Department of Chemistry, Virginia Tech, Blacksburg, VA 24061
- 分类: physics.chem-ph (化学物理)
- 发表时间: 2024年7月11日
- 论文链接: https://arxiv.org/abs/2407.08859
本文对3CC (3-body coupled cluster)方法进行了更全面的评估,以评估其在标准HEAT框架下精确计算热化学的性能。研究实现了适用于闭壳层和开壳层系统的新的自旋积分3CC方法,利用了用于多体电子结构中算符代数推导、优化和评估的新自动化工具链。研究发现,在双zeta基组下,3CC相关能和原子化能贡献几乎总是比CCSDT模型和标准CCSD(T)模型更准确。在HEAT数据集上,相对于CCSDTQ参考值,{3CC, CCSDT, CCSD(T)}的平均绝对误差分别为{24, 70, 122} μEh/e和{0.46, 2.00, 2.58} kJ/mol。
耦合簇(Coupled-Cluster, CC)方法是分子电子结构理论的金标准,在材料科学和核物理中应用日益广泛。它能够紧凑且系统地改进单行列式主导态的电子相关描述。
- CCSD(T)的不足: 尽管CCSD(T)被称为量子化学的"金标准",但在需要高精度或单行列式不再是良好参考的情况下仍不够用
- CCSDT的成本问题: 完整的CCSDT方法虽然更严格,但计算成本高达O(N⁸),且通常精度不如CCSD(T)
- 后CCSD(T)方法的需求: 对于预测性化学能量学等高精度应用,需要超越CCSD(T)的方法
作者在之前的研究中发现3CC方法对小分子集合表现出意外的优异性能,产生的相关能显著接近CCSDTQ参考值。本研究旨在对更大的闭壳层和开壳层系统集合进行更全面的3CC方法评估。
- 实现了新的自旋积分3CC方法: 适用于闭壳层和开壳层系统,基于新的自动化工具链
- 全面的基准测试: 在标准HEAT热化学基准集上对3CC方法进行了系统评估
- 性能优势验证: 证明3CC在精度和计算成本之间提供了比CCSDT更好的平衡
- 技术创新: 开发了SeQuant符号张量代数引擎和相关的自动化实现框架
3CC方法属于nCC方法族,是内部修正的耦合簇方法之一。其设计目标是对n电子系统精确,通过从CC振幅方程中移除具有非空穴共轭(NHCJ)结构的二体及高体簇的二次贡献。
耦合簇波函数定义为:
∣ΨCC⟩≡exp(T^)∣0⟩
其中簇算符包含k体相关子:
T^=∑k=1KT^k
3CC方法通过从CCSDT的三重激发振幅方程中消除特定项来实现:
对于双激发方程,移除的NHCJ贡献为:
A^i1i2a1a2(21gˉi3i4a3a4ti1i3a1a3ti2i4a2a4−41gˉi3i4a3a4ti3i4a1a3ti1i2a2a4)
- 符号算符代数: 高效处理标量和算符场上的张量代数表达式
- 自旋积分: 符号化处理自旋,提高计算效率
- 张量网络优化: 自动确定最优收缩顺序,保证正确的渐近标度
- 在线符号变换: 在执行期间而非代码生成期间进行符号操作
- 张量代数解释器: 直接解释张量表达式,支持分布式内存平台
- MO积分DSL: 简化分子轨道积分的实现
使用HEAT (High accuracy Extrapolated Ab initio Thermochemistry)基准数据集,包含HEAT数据集中除H和H₂外的所有分子种类。
HEAT模型的0K焓定义为:
EHEAT≡EHFCBS+δECCSD(T)CBS+δECCSDTCBS+δECCSDTQ+δEZPVE+δEREL+δESO+δEDBOC
- 绝对相关能误差: 相对于CCSDTQ和CCSDTQP参考值
- 原子化能误差: 化学键形成/断裂能量的预测精度
- 平均绝对误差(MAE): 系统性能的主要指标
- 每价电子误差: 标准化的能量误差指标
- CCSD(T): 标准的单双激发加微扰三重激发
- CCSDT: 完整的单双三重激发耦合簇
- CCSDTQ: 包含非迭代四重激发修正
- CCSDT(Q): 微扰四重激发修正
相对于CCSDTQ参考值的平均绝对误差(mEh):
- CCSD(T): 1.216 (每价电子: 0.122)
- CCSDT: 0.750 (每价电子: 0.070)
- 3CC: 0.261 (每价电子: 0.024)
- CCSDTQ: 0.186 (每价电子: 0.017)
- CCSDT(Q): 0.114 (每价电子: 0.011)
相对于CCSDTQ参考值的平均绝对误差(kJ/mol):
- CCSD(T): 2.58
- CCSDT: 2.00
- 3CC: 0.46
- CCSDTQ: 0.53
- CCSDT(Q): 0.37
相对于HEAT电子能量参考值的原子化能误差(kJ/mol):
- CCSD(T): 1.07 (MAE)
- CCSDT: 2.00 (MAE)
- 3CC: 0.52 (MAE)
- 系统性改进: 3CC在绝对相关能和原子化能方面都显著优于CCSD(T)和CCSDT
- 开壳层优势: 对于开壳层系统,3CC相比CCSD(T)和CCSDT的优势更加明显
- 成本效益: 3CC与CCSDT具有相同的O(N⁸)复杂度,但精度显著更高
- 接近高阶方法: 3CC的原子化能精度接近更昂贵的CCSDTQ方法
对于一些具体分子的表现:
- CN分子: 3CC相对于CCSDTQ的误差仅为0.024 mEh,而CCSDT为1.405 mEh
- CO₂分子: 3CC原子化能误差为2.57 kJ/mol,显著优于CCSDT的4.12 kJ/mol
- 氟化物系统: 3CC在处理强相关的氟化物系统时表现尤为出色
3CC属于Bartlett和Musiał提出的nCC方法族,这些方法通过移除特定的非空穴共轭贡献来改进传统CC方法。
- ACP-D45/ACCD方法: 早期的类似思想,移除CCD双激发方程中的特定项
- pCCSD方法: Huntington和Nooijen的参数化CCSD方法
- 区分簇近似: Kats和Manby的DCSD方法及其扩展
- HEAT协议: 高精度外推从头算热化学
- W3/W4理论: Martin等人的加权理论
- ccCA方法: 相关一致复合方法
- 3CC是最优的迭代三重激发方法: 在相同的O(N⁸)成本下,3CC比CCSDT提供显著更高的精度
- 推荐的后CCSD(T)方法: 对于需要超越CCSD(T)精度的应用,3CC是理想选择
- 系统性改进: 3CC对闭壳层和开壳层系统都提供一致的改进
- 高端热化学: 3CC可以替代CCSDT作为高精度热化学协议的基础
- 四重激发的起点: 3CC应作为包含四重激发修正(如CCSDT(Q))的起始方法
- 成本效益: 相比完整的CCSDTQ,3CC提供了良好的精度-成本平衡
- 基组收敛性: CCSDT和CCSD(T)差异的基组误差不可忽略,需要更大基组的评估
- 实现效率: 虽然形式复杂度相同,但3CC的实际实现可以进一步优化
- 更大系统: 需要开发降阶方法使3CC适用于更大的分子系统
- 算法优化: 利用张量分解和局域化技术降低实际计算成本
- 显式相关: 结合F12技术加速基组收敛
- 自适应方法: 开发自适应CC方法控制精度和成本
- 更大基准: 扩展到更大分子系统的基准测试
- 严格的理论基础: 3CC方法有清晰的理论动机和数学推导
- 全面的基准测试: 使用标准的HEAT数据集进行系统评估
- 技术创新: SeQuant工具链代表了量子化学软件开发的重要进展
- 实用价值: 为高精度量子化学计算提供了新的实用选择
- 系统规模限制: 目前的实现仅适用于小到中等分子(最多10个原子)
- 基组依赖性: 需要更系统地研究基组收敛性质
- 开壳层分析: 对开壳层系统的特殊行为需要更深入的理论分析
- 学术贡献: 为耦合簇理论的发展提供了重要证据
- 软件影响: SeQuant工具链可能影响未来量子化学软件的开发
- 应用前景: 为需要高精度的化学应用提供了新工具
- 需要化学精度(1 kcal/mol)或更高精度的热化学计算
- 过渡金属化合物等需要精确处理电子相关的系统
- 作为更高阶方法(如包含四重激发)的基础方法
- 基准计算和方法开发研究
本文引用了93篇相关文献,主要包括:
- 耦合簇理论的基础文献(Coester, Čížek, Paldus等)
- HEAT热化学协议的相关工作
- nCC方法族的原始论文(Bartlett & Musiał)
- 现代量子化学软件和算法发展
总体评价: 这是一篇高质量的计算化学论文,通过严格的基准测试证明了3CC方法的优越性,为高精度量子化学计算提供了重要的方法学贡献。技术实现上的创新和全面的数值验证使其成为该领域的重要进展。