A perturbative approach to multireference equation-of-motion coupled cluster

We introduce a variant of the multireference equation-of-motion coupled-cluster (MR-EOMCC) method where the amplitudes used for the similarity transformations are estimated from perturbation theory. Consequently, the new variant retains the many-body formalism, a reliance on at most two-body densities, and the state-universal character. As a non-iterative variant, computational costs are reduced, and no convergence difficulties with near-singular amplitudes can arise. Its performance was evaluated on several test sets covering transition metal atoms, small diatomics, and organic molecules against (near-)full CI quality reference data. We further highlight its efficacy on the weakly avoided crossing of LiF and place MR-EOMCC and the new variant into context with linear response theory. The accuracy of the variant was found to be at least on par with expectations for multireference perturbation theories, judging by the NEVPT2 method. The variant can be especially useful in multistate situations where the high accuracy of the iterative MR-EOMCC method is not required.