Post-CCSD(T) corrections in the S66 noncovalent interactions benchmark
Emmanouil Semidalas, A. Daniel Boese, Jan M. L. Martin
Abstract
For noncovalent interactions, it is generally assumed that CCSD(T) approaches the exact solution within the basis set. For most of the S66 benchmark, we present CCSDT and CCSDT(Q) corrections with a DZP basis set. For hydrogen bonds, pure London, and mixed-influence complexes, CCSD(T) benefits from error cancellation between (repulsive) higher-order triples, T 3 − ( T ) , and (attractive) connected quadruples, (Q). For π -stacking complexes, this cancellation starts breaking down and CCSD(T) overbinds; CCSD(T) Λ corrects the problem at the expense of London complexes. Simple two- or three-parameter models predict CCSDT(Q)–CCSD(T) differences to 0.01 kcal mol − 1 RMS, requiring no calculations with steeper scaling than O ( N 7 ) . • Post-CCSD(T) corrections for most of the S66 noncovalent interactions benchmark • For most systems, higher-order triples and connected quadruples nicely cancel • For π -stacking systems, however, higher-order triples are much more important • Hence, the cancellation breaks down, and CCSD(T) tends to overbind • This corroborates recent claims based on diffusion Monte Carlo • Simple estimation formulas for post-CCSD(T) contributions are proposed.