Coupled cluster

About: Coupled cluster is a research topic. Over the lifetime, 6280 publications have been published within this topic receiving 301055 citations.

Extensive generalization of renormalized coupled-cluster methods.

Abstract: The recently developed completely renormalized (CR) coupled-cluster (CC) methods with singles, doubles, and noniterative triples or triples and quadruples [CR-CCSD(T) or CR-CCSD(TQ), respectively], which are based on the method of moments of CC equations (MMCC) [K. Kowalski and P. Piecuch, J. Chem. Phys. 113, 18 (2000)], eliminate the failures of the standard CCSD(T) and CCSD(TQ) methods at larger internuclear separations, but they are not rigorously size extensive. Although the departure from strict size extensivity of the CR-CCSD(T) and CR-CCSD(TQ) methods is small, it is important to examine the possibility of formulating the improved CR-CC methods, which are as effective in breaking chemical bonds as the existing CR-CCSD(T) and CR-CCSD(TQ) approaches, which are as easy to use as the CR-CCSD(T) and CR-CCSD(TQ) methods, and which can be made rigorously size extensive. This may be particularly useful for the applications of CR-CC methods and other MMCC approaches in calculations of potential energy surfaces of large many-electron systems and van der Waals molecules, where the additive separability of energies in the noninteracting limit is very important. In this paper, we propose different types of CR-CC approximations, termed the locally renormalized (LR) CCSD(T) and CCSD(TQ) methods, which become rigorously size extensive if the orbitals are localized on nointeracting fragments. The LR-CCSD(T) and LR-CCSD(TQ) methods rely on the form of the energy expression in terms of the generalized moments of CC equations, derived in this work, termed the numerator-denominator-connected MMCC expansion. The size extensivity and excellent performance of the LR-CCSD(T) and LR-CCSD(TQ) methods are illustrated numerically by showing the results for the dimers of stretched HF and LiH molecules and bond breaking in HF and H2O.

Aspects of separability in the coupled cluster based direct methods for energy differences

Abstract: In this paper we have discussed in detail the aspects of separability of the energy differences obtained from coupled cluster based “direct” methods such as the open-shell Coupled Cluster (CC) theory and the Coupled Cluster based Linear Response Theory (CC-LRT). It has been emphasized that, unlike the state energiesper se, the energy differences have a semi-local character in that, in the asymptotic limit of non-interacting subsystemsA, B, C, etc., they are separable as ΔE A , ΔE B , ΔE A + ΔE B , etc. depending on the subsystems excited. We classify the direct many-body methods into two categories: core-extensive and core-valence extensive. In the former, we only implicitly subtract the ground state energy computed in a size-extensive manner; the energy differences are not chosen to be valence-extensive (separable) in the semi-local sense. The core-valence extensive theories, on the other hand, are fully extensive — i.e., with respect to both core and valence interactions, and hence display the semi-local separability. Generic structures of the wave-operators for core-extensive and core-valence extensive theories are discussed. CC-LRT is shown to be core-extensive after a transcription to an equivalent wave-operator based form. The emergence of valence disconnected diagrams for two and higher valence problems are indicated. The open-shell CC theory is shown to be core-valence extensive and hence fully connected. For one valence problems, the CC theory and the CC-LRT are shown to be equivalent. The equations for the cluster amplitudes in the Bloch equation are quadratic, admitting of multiple solutions. It is shown that the cluster amplitudes for the main peaks, in principle obtainable as a series inV from the zeroth order roots of the model space, are connected, and hence the energy differences are fully extensive. It is remarkable that the satellite energies obtained from the alternative solutions of the CC equations are not valence-extensive, indicating the necessity of a formal power series structure inV of the cluster amplitudes for the valence-extensivity. The alternative solutions are not obtainable as a power series inV. The CC-LRT is shown to have an effective hamiltonian structure respecting “downward reducibility”. A unitary version of CC-LRT (UCC-LRT) is proposed, which satisfy both upward and downward reducibility. UCC-LRT is shown to lead to the recent propagator theory known as the Algebraic Diagrammatic Construction. It is shown that both the main and the satellite peaks from UCC-LRT for the one valence problems are core-valence extensive owing to the hermitized nature of the underlying operator to be diagonalized.

Accuracy of geometries: Influence of basis set, exchange-correlation potential, inclusion of core electrons, and relativistic corrections

Abstract: The geometries of a set of small molecules were optimized using eight different exchange–correlation (xc) potentials in a few different basis sets of Slater-type orbitals, ranging from a minimal basis (I) to a triple-zeta valence basis plus double polarization functions (VII). This enables a comparison of the accuracy of the xc potentials in a certain basis set, which can be related to the accuracies of wavefunction-based methods such as Hartree–Fock and coupled cluster. Four different checks are done on the accuracy by looking at the mean error, standard deviation, mean absolute error and maximum error. It is shown that the mean absolute error decreases with increasing basis set size, and reaches a basis set limit at basis VI. With this basis set, the mean absolute errors of the xc potentials are of the order of 0.7–1.3 pm. This is comparable to the accuracy obtained with CCSD and MP2/MP3 methods, but is still larger than the accuracy of the CCSD(T) method (0.2 pm). The best performing xc potentials are found to be Becke–Perdew, PBE and PW91, which perform as well as the hybrid B3LYP potential. In the second part of this paper, we report the optimization of the geometries of five metallocenes with the same potentials and basis sets, either in a nonrelativistic or a scalar relativistic calculation using the zeroth-order regular approximation approach. For the first-row transition-metal complexes, the relativistic corrections have a negligible effect on the optimized structures, but for ruthenocene they improve the optimized Ru–ring distance by some 1.4–2.2 pm. In the largest basis set used, the absolute mean error is again of the order of 1.0 pm. As the wavefunction-based methods either give a poor performance for metallocenes (Hartree–Fock, MP2), or the size of the system makes a treatment with accurate methods such as CCSD(T) in a reasonable basis set cumbersome, the good performance of density functional theory calculations for these molecules is very promising; even more so as density functional theory is an efficient method that can be used without problems on systems of this size, or larger.

Selection of active spaces for multiconfigurational wavefunctions

Abstract: The efficient and accurate description of the electronic structure of strongly correlated systems is still a largely unsolved problem. The usual procedures start with a multiconfigurational (usually a Complete Active Space, CAS) wavefunction which accounts for static correlation and add dynamical correlation by perturbation theory, configuration interaction, or coupled cluster expansion. This procedure requires the correct selection of the active space. Intuitive methods are unreliable for complex systems. The inexpensive black-box unrestricted natural orbital (UNO) criterion postulates that the Unrestricted Hartree-Fock (UHF) charge natural orbitals with fractional occupancy (e.g., between 0.02 and 1.98) constitute the active space. UNOs generally approximate the CAS orbitals so well that the orbital optimization in CAS Self-Consistent Field (CASSCF) may be omitted, resulting in the inexpensive UNO-CAS method. A rigorous testing of the UNO criterion requires comparison with approximate full configuration interaction wavefunctions. This became feasible with the advent of Density Matrix Renormalization Group (DMRG) methods which can approximate highly correlated wavefunctions at affordable cost. We have compared active orbital occupancies in UNO-CAS and CASSCF calculations with DMRG in a number of strongly correlated molecules: compounds of electronegative atoms (F2, ozone, and NO2), polyenes, aromatic molecules (naphthalene, azulene, anthracene, and nitrobenzene), radicals (phenoxy and benzyl), diradicals (o-, m-, and p-benzyne), and transition metal compounds (nickel-acetylene and Cr2). The UNO criterion works well in these cases. Other symmetry breaking solutions, with the possible exception of spatial symmetry, do not appear to be essential to generate the correct active space. In the case of multiple UHF solutions, the natural orbitals of the average UHF density should be used. The problems of the UNO criterion and their potential solutions are discussed: finding the UHF solutions, discontinuities on potential energy surfaces, and inclusion of dynamical electron correlation and generalization to excited states.

Kramers-restricted closed-shell CCSD theory

Abstract: A Kramers' restricted version of the closed shell coupled cluster singles doubles theory is presented. The theory may be used in conjunction with 2 or 4-component relativistic reference wavefunctions. The intrinsic treatment of the spin-orbit coupling doubles the number of independent quantities (amplitudes and integrals) relative to a spin-independent formalism. The number of operations required to evaluate the equations is four times larger than in the optimal spin-independent closed shell formalism.