Coupled Cluster Doubles theory, linear

The ab initio HF calculations reported below have been performed with the GAUSSIAN 76 [26] program package. The atomic basis sets applied are a minimal (STO-3G [26]) one, a split valence (6-31G [26]) one, a split-valence one plus a set of five d-functions on carbon (6-31G [26]), and one with an additional set of p-functions on hydrogen (6-31G [26]). The correlation energy has been computed using Mpller-Plesset many body perturbation theory of second order (MP2) [27], the linear approximation of Coupled Cluster Doubles theory (L-CCD)... 

Another possibility is to increase the order of PT while restricting the excitation level of the intervening intermediate states. This is most easily done when the excitations are limited to doubles, yielding the rath order MBPT with doubles, or DMBPT(n). In fact, in this case, the summation can be carried out to infinite order, yielding the DMBPT(oo) method [27]. However, this result can be easily seen to be equivalent to the linear version of coupled-cluster (CC) theory that is restricted to two-body amplitudes, namely, to the L-CPMET (linear coupled-pair many-electron theory) or, more succinctly, to L-CCD (linear CC with doubles) in current terminology. 

In the following, the theory of Kutzelnigg s linear R12 functions shall be presented and analyzed in the framework of the coupled-cluster doubles (CCD) method, To illustrate the ideas and approximations employed in the linear R12 methods, it is sufficient to consider the CCD model, as the corresponding CCD-R12 theory exhibits all properties of the R12 theories. It is a relatively simple matter to include singles (CCSD) or even triples (for example in the CCSD(T) method), and CI-R12-type wave functions or MPn-R12 energies require essentially the same computational procedures as the CCD-R12 approach. 

However, until today no systematic comparison of methods based on MpUer-Plesset perturbation (MP) and Coupled Cluster theory, the SOPPA or multiconfigurational linear response theory has been presented. The present study is a first attempt to remedy this situation. Calculations of the rotational g factor of HF, H2O, NH3 and CH4 were carried out at the level of Hartree-Fock (SCF) and multiconfigurational Hartree-Fock (MCSCF) linear response theory, the SOPPA and SOPPA(CCSD) [40], MpUer-Plesset perturbation theory to second (MP2), third (MP3) and fourth order without the triples contributions (MP4SDQ) and finally coupled cluster singles and doubles theory. The same basis sets and geometries were employed in all calculations for a given molecule. The results obtained with the different methods are therefore for the first time direct comparable and consistent conclusions about the performance of the different methods can be made. 

Pedersen, T. B., Sanchez de Meras, A. M.and Koch, H. (2004). Polarizability and optical rotation calculated from the approximate coupled cluster singles and doubles CC2 linear response theory using Cholesky decompositions, /. Chem. Phys. 120, pp. 8887-8897, doiilO.1063/1. 1705575. 

It is important to note that, at each level of coupled-cluster theory, we include through the exponential parameterization of Eq. (28) all possible determinants that can be generated within a given orbital basis, that is, all determinants that enter the FCI wave function in the same orbital basis. Thus, the improvement in the sequence CCSD, CCSDT, and so on does not occur because more determinants are included in the description but from an improved representation of their expansion coefficients. For example, in CCS theory, the doubly-excited determinants are represented by ]HF), whereas the same determinants are represented by (T2 + Tf) HF) in CCSD theory. Thus, in CCSD theory, the weight of each doubly-excited determinant is obtained as the sum of a connected doubles contribution from T2 and a disconnected singles contribution from Tf/2. This parameterization of the wave function is not only more compact than the linear parameterization of configuration-interaction (Cl) theory, but it also ensures size-extensivity of the calculated electronic energy. 

The perturbation operator is reminiscent of the doubles cluster operator (13.2.7) of coupled-cluster theory. It is a linear combination of excitation operators, each multiplied by a first-order amplitude of the form... 

These equations should be interpreted to indicate that the double amplitudes are at least linear in the fluctuation potential and that the singles are at least quadratic in the fluctuation potential. The order equations (14.6.14) and (14.6.15) should not be interpreted as giving any direct indication of the accuracy of the singles and doubles amplitudes. To establish their accuracy, we go to the next level of coupled-cluster theory (CCSDT) and determine the lowest orders of the corrections introduced. Referring to (14.6.11)-(14.6.13), we find... 

---