On the applicability of CCSD(T) for dispersion interactions in large conjugated systems

In light of the recent discrepancies reported between fixed node diffusion Monte Carlo and local natural orbital coupled cluster with single, double, and perturbative triples [CCSD(T)] methodologies for non-covalent interactions in large molecular systems [Al-Hamdani et al., Nat. Commun. 12, 3927 (2021)], the applicability of CCSD(T) is assessed using a model framework. The use of the semi-empirical π-space only Pariser-Parr-Pople (PPP) model for studying large molecules is critically examined and is shown to recover both bandgap closure as system size increases and long range dispersive behavior of r-6 with increasing separation between monomers. Since bandgap closure in systems with long-range Coulomb interactions is problematic for perturbative methods, such as CCSD(T), this model, therefore, serves as a testing ground for such methods, enabling them to be benchmarked with high-order CC methods, which are not possible with ab initio Hamiltonians. Using the PPP model, coupled cluster methodologies, CCSDTQ and CCSDT(Q), are then used to benchmark CCSDT and CCSD(T) methodologies for non-covalent interactions in large one- and two-dimensional molecular systems up to the dibenzocoronene dimer. We show that CCSD(T) demonstrates no signs of overestimating the interaction energy for these systems. Furthermore, by examining the Hartree-Fock HOMO-LUMO gap of these large molecules, the perturbative treatment of the triples contribution in CCSD(T) is not expected to cause problems for accurately capturing the interaction energy for system sizes up to at least circumcoronene.