Tractability gains in symmetry-adapted perturbation theory including coupled double excitations: CCD+ST(CCD) dispersion with natural orbital truncations

Abstract

This work focuses on efficient and accurate treatment of the intermolecular dispersion interaction using the CCD+ST(CCD) dispersion approach formulated by Williams et al. [J. Chem. Phys. 103, 4586 (1995)]. We apply natural orbital truncation techniques to the solution of the monomer coupled-cluster double (CCD) equations, yielding substantial accelerations in this computationally demanding portion of the SAPT2+(CCD), SAPT2+(3)(CCD), and SAPT2+3(CCD) analyses. It is shown that the wholly rate-limiting dimer-basis particle-particle ladder term can be computed in a reduced natural virtual space which is essentially the same size as the monomer-basis virtual space, with an error on the order of a few thousandths of 1 kcal mol(-1). Coupled with our existing natural orbital techniques for the perturbative triple excitation contributions [E. G. Hohenstein and C. D. Sherrill, J. Chem. Phys. 133, 104107 (2010)], this technique provides speedups of greater than an order of magnitude for the evaluation of the complete SAPT2+3(CCD) decomposition, with a total error of a few hundredths of 1 kcal mol(-1). The combined approach yields tractability gains of almost 2× in the system size, allowing for SAPT2+3(CCD)/aug-cc-pVTZ analysis to be performed for systems such as adenine-thymine for the first time. Natural orbital based SAPT2+3(CCD)/aug-cc-pVTZ results are presented for stacked and hydrogen-bonded configurations of uracil dimer and the adenine-thymine dimer.