Benchmarking for perturbative triple-excitations in EE-EOM-CC methods

Perturbative triples corrections ((T) and (T̃)) to the equation of motion coupled cluster singles and doubles (EOM-CCSD) are rederived and implemented in a pilot parallel code. The vertical excitation energies of molecules in the test set of Sauer et al. [J. Chem. Theor. Comput. 2009, 5, 555] are reported and compared to the iterative EOM-CCSDT-3 method. The average absolute deviations of EOM-CCSD(T) and EOM-CCSD(T̃) from EOM-CCSDT-3 over this wide test set are 0.06 and 0.18 eV, respectively. The poor performance of the latter suggests misbalanced handling of the (T̃) terms. Scaling curves for EOM-CCSD(T) are also presented to demonstrate the performance across multiple compute nodes, thus enabling the routine and accurate study of excited states for ever larger molecular systems.