CCSD coupled-cluster singles, doubles calculations

There have been many calculations of a for neon itself at various levels of approximation. Some of these are summarized in Ref. [26] and a basis set analysis has been made in Ref. [25]. Of the electron-correlated results[42]-[47], there appears to be a consensus that the best value using coupled-cluster-single-doubles wavefunctions with a perturbational estimate of connected triplet excitations, CCSD(T), is 2.680-2.690 au[25],[44],[46],[47]. 

In Table 6.3, the values of De for RfCU are compared with those obtained within various approximations using relativistic effective core potentials (RECP) Kramers-restricted Hartree-Fock (KRHF) (Han et al 1999), averaged RECP including second-order M0ller-Plesset perturbation theory (AREP-MP2) for the correlation part (Han et al. 1999), RECP coupled-cluster single double (triple) [CCSD(T)] excitations (Han et al. 1999), and a Dirac-Fock-Breit (DFB) method (Malli and Styszynski 1998). The AREP-MP2 calculation of De gives 20.4 eV, while the RECP-CCSD(T) method with correlation leads to 18.8 eV. Our value of De of 19.5 eV is just between these calculated values. 

Atomic calculations. Most atomic calculations for the heaviest elements were performed by using Dirac-Fock (DF) and Dirac-Slater (DS) methods [20-24,58] and later by using multiconfiguration Dirac-Fock (MCDF) [64-72] and Dirac-Coulomb-Breit Coupled Cluster Single Double excitations (DCB CCSD) [73-85] methods, with the latter being presently the most accmate one. 

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. 

As has been demonstrated by numerous studies, the accuracy of properties calculated using the GGA DFT methods is, in most cases, comparable to or better than those from ab initio MP2 (Mpller-Plesset second-order perturbation) or CISD (configuration interaction with single and double excitations) methods. In fact, the accuracy of the DFT results in some instances matches those obtained from the much more costly (but, in principle, more exact) CCSD(T) (coupled cluster singles and doubles with a perturbative inclusion of connected triple excitations) method (29) and the ab initio G1 procedure (30). 

Analytic Energy Gradients for Open-Shell Coupled-Cluster Singles and Doubles (CCSD) Calculations Using Restricted Open-Shell Hartree-Fock (ROHF) Reference Functions. 

---