X

Search Results

Searching....

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

7.11.3 EOM-XX-CCSD and CI Suite of Methods

(April 13, 2024)

Q-Chem features the most complete set of EOM-CCSD models, 660 Krylov A. I.
Annu. Rev. Phys. Chem.
(2008), 59, pp. 433.
Link
enabling accurate, robust, and efficient calculations of electronically excited states (EOM-EE-CCSD or EOM-EE-OD); 1106 Sekino H., Bartlett R. J.
Int. J. Quantum Chem. Symp.
(1984), 18, pp. 255.
Link
, 631 Koch H. et al.
J. Chem. Phys.
(1990), 93, pp. 3345.
Link
, 1159 Stanton J. F., Bartlett R. J.
J. Chem. Phys.
(1993), 98, pp. 7029.
Link
, 656 Krylov A. I., Sherrill C. D., Head-Gordon M.
J. Chem. Phys.
(2000), 113, pp. 6509.
Link
, 727 Levchenko S. V., Krylov A. I.
J. Chem. Phys.
(2004), 120, pp. 175.
Link
; ground and excited states of diradicals and triradicals (EOM-SF-CCSD and EOM-SF-OD); 657 Krylov A. I.
Chem. Phys. Lett.
(2001), 338, pp. 375.
Link
, 727 Levchenko S. V., Krylov A. I.
J. Chem. Phys.
(2004), 120, pp. 175.
Link
ionization potentials and electron attachment energies, as well as problematic doublet radicals and cation or anion radicals (EOM-IP/EA-CCSD). 1133 Sinha D. et al.
Chem. Phys. Lett.
(1989), 154, pp. 544.
Link
, 1161 Stanton J. F., Gauss J.
J. Chem. Phys.
(1994), 101, pp. 8938.
Link
, 894 Nooijen M., Bartlett R. J.
J. Chem. Phys.
(1995), 102, pp. 3629.
Link
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.54)

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.55)
Φijab|H¯|Φ0 = 0 (7.56)

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. 1159 Stanton J. F., Bartlett R. J.
J. Chem. Phys.
(1993), 98, pp. 7029.
Link
In the ionized/electron attached EOM models, 1161 Stanton J. F., Gauss J.
J. Chem. Phys.
(1994), 101, pp. 8938.
Link
, 894 Nooijen M., Bartlett R. J.
J. Chem. Phys.
(1995), 102, pp. 3629.
Link
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 method 657 Krylov A. I.
Chem. Phys. Lett.
(2001), 338, pp. 375.
Link
, 727 Levchenko S. V., Krylov A. I.
J. Chem. Phys.
(2004), 120, pp. 175.
Link
in which the excitation operators include spin-flip allows one to access diradicals, triradicals, and bond-breaking.

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. 728 Levchenko S. V., Wang T., Krylov A. I.
J. Chem. Phys.
(2005), 122, pp. 224106.
Link
Properties calculations (permanent and transition dipole moments and angular momentum projections, S^2, R^2, 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 (Section 6.10.8) or Cholesky decomposition (Section 6.10.9) 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.11.13).

7.11.3.1 EOM-CC and projection-based embedding

Due to the high computational cost of this method, the application of EOM-CCSD to large systems is difficult in terms of both computational time and resources. One strategy to overcome this issue is to combine EOM-CC with projection-based embedding 791 Manby F. R. et al.
J. Chem. Theory Comput.
(2012), 8, pp. 2564.
Link
, 933 Parravicini V., Jagau T. C.
Mol. Phys.
(2021), 119, pp. e1943029.
Link
. This method performs remarkably well for the calculation of ionized, core-ionized and valence excitation energies. Rydberg states are also described with sufficient accuracy, while it is not advised to employ projection-based embedding combined with EOM-EA-CCSD. Theory and job commands for projection-based embedding are reported in Section 11.6.

7.11.3.2 Legacy features available in CCMAN

The CCMAN module of Q-Chem includes two implementations of EOM-IP-CCSD. The proper implementation 980 Pieniazek P. A., Bradforth S. E., Krylov A. I.
J. Chem. Phys.
(2008), 129, pp. 074104.
Link
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. 418 Golubeva A. A., Pieniazek P. A., Krylov A. I.
J. Chem. Phys.
(2009), 130, pp. 124113.
Link