Fan Peng-Dong, Kamiya Muneaki, Hirata So
Quantum Theory Project, Department of Chemistry, University of Florida, Gainesville, Florida 32611-8435.
J Chem Theory Comput. 2007 May;3(3):1036-46. doi: 10.1021/ct600270c.
Several variants of the equation-of-motion coupled-cluster (EOM-CC) method with singles (one-hole or one-particle), doubles (two-hole-one-particle or two-particle-one-hole), and a selected set of triples (three-hole-two-particle or three-particle-two-hole) and/or quadruples (four-hole-three-particle or four-particle-three-hole) have been implemented by computerized symbolic algebra. They are applicable to excitation energies (EE), ionization potentials (IP), and electron affinities (EA), excited-state dipole moments, and transition dipole moments of both closed- and open-shell species and are abbreviated as EE/IP/EA-EOM-CCSDt, EE/IP/EA-EOM-CCSDtq, and EE/IP/EA-EOM-CCSDTq, where the small letters indicate the use of active-space cluster and EE/IP/EA operators. They are also parallel executable and accelerated by the use of spin, spatial, and permutation symmetries. The remarkable effectiveness of the methods in capturing nondynamical correlation effects has been demonstrated by their applications to the vertical excitation energies of C2, the adiabatic excitation energies and dipole moments of the CH radical, the adiabatic excitation energies of the CH2 diradical, the adiabatic excitation energies and dipole moments of formaldehyde, the vertical ionization energies of N2, and the vertical electron affinities of C2. The effectiveness is found to decline when the basis set is extended, causing the active space to become relatively small and also less well-defined. As a remedy, we propose a composite method that combines higher-rank active-space methods with smaller basis sets for nondynamical correlation and lower-rank nonactive-space methods with larger basis sets for dynamical correlation, which is shown to work well for an excited-state potential energy curve of hydrogen fluoride.
通过计算机符号代数实现了运动方程耦合簇(EOM - CC)方法的几种变体,这些变体包括单激发(一孔或一粒子)、双激发(两孔 - 一粒子或两粒子 - 一孔)以及一组选定的三激发(三孔 - 两粒子或三粒子 - 两孔)和/或四激发(四孔 - 三粒子或四粒子 - 三孔)。它们适用于闭壳层和开壳层物种的激发能(EE)、电离势(IP)和电子亲和能(EA)、激发态偶极矩以及跃迁偶极矩,分别简称为EE/IP/EA - EOM - CCSDt、EE/IP/EA - EOM - CCSDtq和EE/IP/EA - EOM - CCSDTq,其中小写字母表示使用活性空间簇和EE/IP/EA算符。这些方法还可以并行执行,并通过利用自旋、空间和置换对称性进行加速。通过将这些方法应用于C2的垂直激发能、CH自由基的绝热激发能和偶极矩、CH2双自由基的绝热激发能、甲醛的绝热激发能和偶极矩、N2的垂直电离能以及C2的垂直电子亲和能,证明了这些方法在捕捉非动力学相关效应方面具有显著效果。当基组扩展时,发现这种效果会下降,导致活性空间变得相对较小且定义也不太明确。作为一种补救措施,我们提出了一种复合方法,该方法将高阶活性空间方法与较小基组用于非动力学相关,将低阶非活性空间方法与较大基组用于动力学相关,结果表明该方法对于氟化氢的激发态势能曲线效果良好。