A full configuration interaction calculation (FCI) ultimately defines the innate molecular orbital description of a molecule. Its density matrix and the natural orbitals obtained from it quantify the difference between having N-dominantly occupied orbitals in a reference determinant for a wavefunction to describe N-correlated electrons and how many of those N-electrons are left to the remaining virtual orbitals. The latter provides a measure of the multi-determinantal character (MDC) required to be in a wavefunction. MDC is further split into a weak correlation part and a part that indicates stronger correlation often called multi-reference character (MRC). If several virtual orbitals have high occupation numbers, then one might argue that these additional orbitals should be allowed to have a larger role in the calculation, as in MR methods, such as MCSCF, MR-CI, or MR-coupled-cluster (MR-CC), to provide adequate approximations toward the FCI. However, there are problems with any of these MR methods that complicate the calculations compared to the uniformity and ease of application of single-reference CC calculations (SR-CC) and their operationally single-reference equation-of-motion (EOM-CC) extensions. As SR-CC theory is used in most of today’s “predictive” calculations, an assessment of the accuracy of SR-CC at some truncation of the cluster operator would help to quantify how large an issue MRC actually is in a calculation, and how it might be alleviated while retaining the convenient SR computational character of CC/EOM-CC. This paper defines indices that identify MRC situations and help assess how reliable a given calculation is.

Solution of dark, doubly excited states usually requires EOM-CCSDT or EOM-CCSDTQ for appropriate treatment. A new route to obtain these doubly excited states using EOM-CCSD is demonstrated. Traditionally, EOM-CC is performed on a closed shell reference state for it has a well-described single reference CC wavefunction. In this note, we attempt to use low spin open-shell states such as MS=0 triplet and open-shell singlet as a reference state. Using this intermediate excited state as a reference state provides us with the benefit of obtaining a doubly excited state, as a single excitation with the cost of EOM-CCSD.

To treat doubly excited states, the treatment of triple excitations is considered necessary in the framework of equation-of-motion coupled cluster (EOM-CC) methods. We investigate models without explicit triples and seek quantitative measure for the performance of EOM based on CC with singles and doubles (CCSD) or modified CCSD (Distinguishable Cluster Approximation) approaches for states with predominant double excitation character. We also test the efficacy of including triples in perturbative manner through EOM-CCSD(T) and in an iterative way through EOM-CCSDT-3 method. Extended similarity transformed EOM-CCSD(EXT-STEOM-CCSD) method is also tested and provides superior quality results at comparatively low cost. We use the QUEST2 benchmark set of double excitations proposed by Loos et al. [J. Chem. Theory Comput.2019151939] to investigate the performance of methods such as EOM-CCSD, EOM-DCSD, EXT-STEOM-CCSD, ΔCCSD, and ΔDCSD. We also test a tailored CC approach, ΔpairCCD-TCCSD.