RESUMEN
We present here a massively parallel implementation of the recently developed CPS(D-3) excitation energy model that is based on cluster perturbation theory. The new algorithm extends the one developed in Baudin et al. [J. Chem. Phys., 150, 134110 (2019)] to leverage multiple nodes and utilize graphical processing units for the acceleration of heavy tensor contractions. Furthermore, we show that the extended algorithm scales efficiently with increasing amounts of computational resources and that the developed code enables CPS(D-3) excitation energy calculations on large molecular systems with a low time-to-solution. More specifically, calculations on systems with over 100 atoms and 1000 basis functions are possible in a few hours of wall clock time. This establishes CPS(D-3) excitation energies as a computationally efficient alternative to those obtained from the coupled-cluster singles and doubles model.
RESUMEN
The cluster perturbation series, CPS(D), for coupled cluster singles and doubles excitation energies is considered. It is demonstrated that the second-order model CPS(D-2) is identical to the configuration interaction singles with perturbative doubles, CIS(D) model. The third-order model, CPS(D-3), provides excitation energies of coupled cluster singles and doubles (CCSD) quality in the sense that the difference between CPS(D-3) and CCSD excitation energies is of the same size or smaller than the effect of adding triples corrections to CCSD excitation energies. We further show that the third-order corrections can be efficiently implemented, in particular, when the resolution of the identity approximation is used for integrals. We also show that the CPS(D-3) excitation energies can be determined for system sizes that are far beyond what can be considered in conventional CCSD excitation energy calculations.
