7.9 Coupled-Cluster Excited-State and Open-Shell Methods

7.9.2 EOM-XX-CCSD and CI Suite of Methods

Q-Chem features the most complete set of EOM-CCSD models,Krylov:2008 enabling accurate, robust, and efficient calculations of electronically excited states (EOM-EE-CCSD or EOM-EE-OD); Sekino:1984, Koch:1990b, Stanton:1993, Krylov:2000, Levchenko:2004; ground and excited states of diradicals and triradicals (EOM-SF-CCSD and EOM-SF-OD);Krylov:2001a, Levchenko:2004 ionization potentials and electron attachment energies, as well as problematic doublet radicals and cation or anion radicals (EOM-IP/EA-CCSD).Sinha:1989, Stanton:1994, Nooijen:1995 The EOM-DIP-CCSD, EOM-2SF-CCSD, and EOM-DEA-CCSD methods are available as well. Conceptually, EOM is very similar to configuration interaction (CI): target EOM states are found by diagonalizing the similarity transformed Hamiltonian H¯=e-THeT,

H¯R=ER, (7.45)

where T and R are general excitation operators with respect to the reference determinant |Φ0. In the EOM-CCSD models, T and R are truncated at single and double excitations, and the amplitudes T satisfy the CC equations for the reference state |Φ0:

Φia|H¯|Φ0 = 0 (7.46)
Φijab|H¯|Φ0 = 0 (7.47)

The computational scaling of EOM-CCSD and CISD methods is identical, i.e., 𝒪(N6), however EOM-CCSD is numerically superior to CISD because correlation effects are “folded in” in the transformed Hamiltonian, and because EOM-CCSD is rigorously size-intensive.

By combining different types of excitation operators and references |Φ0, different groups of target states can be accessed as explained in Fig. 7.1. For example, electronically excited states can be described when the reference |Φ0 corresponds to the ground state wave function, and operators R conserve the number of electrons and a total spin.Stanton:1993 In the ionized/electron attached EOM models,Stanton:1994, Nooijen:1995 operators R are not electron conserving (i.e., include different number of creation and annihilation operators)—these models can accurately treat ground and excited states of doublet radicals and some other open-shell systems. For example, singly ionized EOM methods, i.e., EOM-IP-CCSD and EOM-EA-CCSD, have proven very useful for doublet radicals whose theoretical treatment is often plagued by symmetry breaking. Finally, the EOM-SF methodKrylov:2001a, Levchenko:2004 in which the excitation operators include spin-flip allows one to access diradicals, triradicals, and bond-breaking.Krylov:2017

Q-Chem features EOM-EE/SF/IP/EA/DIP/DSF-CCSD methods for both closed and open-shell references (RHF/UHF/ROHF), including frozen core/virtual options. For EE, SF, IP, and EA, a more economical flavor of EOM-CCSD is available (EOM-MP2 family of methods). All EOM models take full advantage of molecular point group symmetry. Analytic gradients are available for RHF and UHF references, for the full orbital space, and with frozen core/virtual orbitals.Levchenko:2005 Properties calculations (permanent and transition dipole moments and angular momentum projections, S2, R2, etc.) are also available. The current implementation of the EOM-XX-CCSD methods enables calculations of medium-size molecules, e.g., up to 15–20 heavy atoms. Using RI approximation 6.8.5 or Cholesky decomposition 6.8.6 helps to reduce integral transformation time and disk usage enabling calculations on much larger systems. EOM-MP2 and EOM-MP2t variants are also less computationally demanding. The computational cost of EOM-IP calculations can be considerably reduced (with negligible decline in accuracy) by truncating virtual orbital space using FNO scheme (see Section 7.9.9). Legacy features available in CCMAN

The CCMAN module of Q-Chem includes two implementations of EOM-IP-CCSD. The proper implementationPieniazek:2008 is used by default is more efficient and robust. The EOM_FAKE_IPEA keyword invokes is a pilot implementation in which EOM-IP-CCSD calculation is set up by adding a very diffuse orbital to a requested basis set, and by solving EOM-EE-CCSD equations for the target states that include excitations of an electron to this diffuse orbital. The implementation of EOM-EA-CCSD in CCMAN also uses this trick. Fake IP/EA calculations are only recommended for Dyson orbital calculations and debug purposes. (CCMAN2 features proper implementations of EOM-IP and EOM-EA (including Dyson orbitals)).

A more economical CI variant of EOM-IP-CCSD, IP-CISD is also available in CCMAN. This is an N5 approximation of IP-CCSD, and can be used for geometry optimizations of problematic doublet states.Golubeva:2009