Is It Possible To Obtain Coupled Cluster Quality Energies at near Density Functional Theory Cost? Domain-Based Local Pair Natural Orbital Coupled Cluster vs Modern Density Functional Theory

Abstract

The recently developed domain-based local pair natural orbital coupled cluster theory with single, double, and perturbative triple excitations (DLPNO-CCSD(T)) delivers results that are closely approaching those of the parent canonical coupled cluster method at a small fraction of the computational cost. A recent extended benchmark study established that, depending on the three main truncation thresholds, it is possible to approach the canonical CCSD(T) results within 1 kJ (default setting, TightPNO), 1 kcal/mol (default setting, NormalPNO), and 2-3 kcal (default setting, LoosePNO). Although thresholds for calculations with TightPNO are 2-4 times slower than those based on NormalPNO thresholds, they are still many orders of magnitude faster than canonical CCSD(T) calculations, even for small and medium sized molecules where there is little locality. The computational effort for the coupled cluster step scales nearly linearly with system size. Since, in many instances, the coupled cluster step in DLPNO-CCSD(T) is cheaper or at least not much more expensive than the preceding Hartree-Fock calculation, it is useful to compare the method against modern density functional theory (DFT), which requires an effort comparable to that of Hartree-Fock theory (at least if Hartree-Fock exchange is part of the functional definition). Double hybrid density functionals (DHDF's) even require a MP2-like step. The purpose of this article is to evaluate the cost vs accuracy ratio of DLPNO-CCSD(T) against modern DFT (including the PBE, B3LYP, M06-2X, B2PLYP, and B2GP-PLYP functionals and, where applicable, their van der Waals corrected counterparts). To eliminate any possible bias in favor of DLPNO-CCSD(T), we have chosen established benchmark sets that were specifically proposed for evaluating DFT functionals. It is demonstrated that DLPNO-CCSD(T) with any of the three default thresholds is more accurate than any of the DFT functionals. Furthermore, using the aug-cc-pVTZ basis set and the LoosePNO default settings, DLPNO-CCSD(T) is only about 1.2 times slower than B3LYP. With NormalPNO thresholds, DLPNO-CCSD(T) is about a factor of 2 slower than B3LYP and shows a mean absolute deviation of less than 1 kcal/mol to the reference values for the four different data sets used. Our conclusion is that coupled cluster energies can indeed be obtained at near DFT cost.