English

Factorized Quadruples and a Predictor of Higher-Level Correlation in Thermochemistry

Chemical Physics 2025-03-26 v1

Abstract

Coupled cluster theory has had a momentous impact on the ab initio prediction of molecular properties, and remains a staple ingratiate in high-accuracy thermochemical model chemistries. However, these methods require inclusion of at least some connected quadruple excitations, which generally scale at best as O(N9)\mathcal{O}(N^9) with the number of basis functions. It very difficult to predict, a priori, the effect correlation past CCSD(T) has on a give reaction energies. The purpose of this work is to examine cost-effective quadruple corrections based on the factorization theorem of many-body perturbation theory that may address these challenges. We show that the O(N7)\mathcal{O}(N^7), factorized CCSD(TQf{}_\text{f}) method introduces minimal error to predicted correlation and reaction energies as compared to the O(N9)\mathcal{O}(N^9) CCSD(TQ). Further, we examine the performance of Goodson's continued fraction method in the estimation of CCSDT(Q)Λ{}_\Lambda contributions to reaction energies, as well as a "new" method related to %TAE[(T)] that we refer to as a scaled perturbation estimator. We find that the scaled perturbation estimator based upon CCSD(TQf{}_\text{f})/cc-pVDZ is capable of predicting CCSDT(Q)Λ{}_\Lambda/cc-pVDZ contributions to reaction energies with an average error of 0.07 kcal mol1{}^{-1} and a RMST of 0.52 kcal mol1{}^{-1} when applied to a test-suite of nearly 3000 reactions. This offers a means by which to reliably ballpark how important post-CCSD(T) contributions are to reaction energies while incurring no more than CCSD(T) formal cost and a little mental math.

Keywords

Cite

@article{arxiv.2407.03519,
  title  = {Factorized Quadruples and a Predictor of Higher-Level Correlation in Thermochemistry},
  author = {James H. Thorpe and Zachary W. Windom and Rodney J. Bartlett and Devin A. Matthews},
  journal= {arXiv preprint arXiv:2407.03519},
  year   = {2025}
}
R2 v1 2026-06-28T17:28:34.911Z