An efficient closed-shell singles and doubles coupled-cluster method

A reformulated set of equations for the closed-shell singles and doubles coupled-cluster (CCSD) method is presented. A computational cost of n_v⁴n₀²+7n_v³n₀³+1n_v²n₀⁴ for the n⁶ steps is obtained, where n_v is the number of virtual molecular orbitals included in the CCSD procedure, n₀ is the number of doubly occupied molecular orbitals and n=n₀+n_v. Test calculations for the cis and trans isomers of FNNF and planar and pyramidal CH₃⁻ are presented. Equilibrium structures determined with large Gaussian basis sets at the second-order Møller-Plesset (MP2) perturbation level of theory are reported and used for the other electron correlation methods. With the largest one-particle basis set (144 contracted Gaussian functions), the equilibrium geometries of cis- and trans-FNNF agree with experiment. Based on analyses of planar and pyramidal CH₃⁻ wavefunctions and the calculated inversion barrier, it is suggested that the molecular anion may not exist in a planar configuration but that autodetachment of an electron occurs before the transition state is reached. Comparisons of our new CCSD procedure demonstrate that coupled-cluster methods are not significantly more expensive than similar electron correlation techniques.