Coupled clustered theory

Bartlett R J 1995 Coupled-cluster theory an overview of recent developments Modem Eiectronic Structure Theory vo 2, ed D R Yarkony (Singapore World Scientific) pp 1047-131... 

A number of types of calculations begin with a HF calculation and then correct for correlation. Some of these methods are Moller-Plesset perturbation theory (MPn, where n is the order of correction), the generalized valence bond (GVB) method, multi-conhgurational self-consistent held (MCSCF), conhgu-ration interaction (Cl), and coupled cluster theory (CC). As a group, these methods are referred to as correlated calculations. 

The QCISD method is also very closely related to coupled cluster theory, with singles and doubles (CCSD). In contrast to QCISD. 

Analogously to MP methods, coupled cluster theory may also be based on a UFIF reference wave function. The resulting UCC methods again suffer from spin contamination of the underlying UHF, but the infinite nature of coupled cluster methods is substantially better at reducing spin contamination relative to UMP. Projection methods analogous to those of the PUMP case have been considered but are not commonly used. ROHF based coupled cluster methods have also been proposed, but appear to give results very similar to UCC, especially at the CCSD(T) level. 

This raises a dilemma in treating second- and higher-order properties in coupled-cluster theory. In the EOM-CC approach, which is basically a Cl calculation for a non-Hermitian Hamiltonian H= that incorporates... 

We do not pretend to give here an exhaustive account of all the possible applications ofNSS s into Quantum Chemistry. Some areas, which for sure can be studied from the nested summation point of view, like the Coupled Cluster Theory [14], are not included here. 

T. Daniel Crawford and Henry F. Schaefer III, An Introduction to Coupled Cluster Theory for Computational Chemists. 

In Table 1.2, we have listed the valence cc-pVDZ electronic energies and AEs of N2 and HF at different levels of coupled-cluster theory. The energies are given as deviations from the FCI values. Comparing the different levels of theory, we note that the error is reduced by one order of magnitude at each level. In particular, at the CCSDT level, there is a residual error of the order of a few kJ/mol in the calculated energies and AEs, suggesting that the CCSDTQ model is usually needed to reproduce experimental measurements to within the quoted errors bars (often less than 1 kJ/mol). 

To improve on the CCSD description, we go to the next level of coupled-cluster theory, including corrections from triple excitations -see the third row of Table 1.4, where we have listed the triples corrections to the energies as obtained at the CCSD(T) level. The triples corrections to the molecular and atomic energies are almost two orders of magnitude smaller than the singles and doubles corrections. However, for the triples, there is less cancellation between the corrections to the molecule and its atoms than for the doubles. The total triples correction to the AE is therefore only one order of magnitude smaller than the singles and doubles corrections. 

There is also a hierarchy of electron correlation procedures. The Hartree-Fock (HF) approximation neglects correlation of electrons with antiparallel spins. Increasing levels of accuracy of electron correlation treatment are achieved by Mpller-Plesset perturbation theory truncated at the second (MP2), third (MP3), or fourth (MP4) order. Further inclusion of electron correlation is achieved by methods such as quadratic configuration interaction with single, double, and (perturbatively calculated) triple excitations [QCISD(T)], and by the analogous coupled cluster theory [CCSD(T)] [8],... 

The field of quantum chemistry has seen tremendous development over the last thirty years. Thanks to high-accuracy models such as coupled-cluster theory and standardized, widely available program packages such as Gaussian 98, what was once merely an esoteric tool of a few specialists has evolved into an indispensable source of knowledge for both the prediction and the interpretation of chemical phenomena. With the development of reduced scaling algorithms for coupled cluster... 

Very accurate values of the dipole and quadrupole polarizability for the equilibrium internuclear distance of HF can be found in a review article by Maroulis [71], calculated with finite-field Mpller-Plesset perturbation theory at various orders and coupled cluster theory using a carefully selected basis set. 

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. 

Improving the T and Ti Components by the Extended Coupled-Cluster Theory... 

Static charge-density susceptibilities have been computed ab initio by Li et al (38). The frequency-dependent susceptibility x(r, r cd) can be calculated within density functional theory, using methods developed by Ando (39 Zang-will and Soven (40 Gross and Kohn (4I and van Gisbergen, Snijders, and Baerends (42). In ab initio work, x(r, r co) can be determined by use of time-dependent perturbation techniques, pseudo-state methods (43-49), quantum Monte Carlo calculations (50-52), or by explicit construction of the linear response function in coupled cluster theory (53). Then the imaginary-frequency susceptibility can be obtained by analytic continuation from the susceptibility at real frequencies, or by a direct replacement co ico, where possible (for example, in pseudo-state expressions). 

Use of Equation (1) in numerical work requires a means of generating x(r, r i(o) as well as the average charge density. Direct variational methods are not applicable to the expression for E itself, due to use of the virial theorem. However, both pc(r) and x(r, r ico) (39-42, 109-112) are computable with density-functional methods, thus permitting individual computations of E from Eq. (1) and investigations of the effects of various approximations for x(r, r ico). Within coupled-cluster theory, x(r, r ico) can be generated directly (53) from the definition in Eq. (3) then Eq. (1) yields the coupled-cluster energy in a new form, as an expectation value. 

---