Active-space Coupled-cluster Methods Research Articles

Novel strategy to implement active-space coupled-cluster methods.

A new approach is presented for the efficient implementation of coupled-cluster (CC) methods including higher excitations based on a molecular orbital space partitioned into active and inactive orbitals. In the new framework, the string representation of amplitudes and intermediates is used as long as it is beneficial, but the contractions are evaluated as matrix products. Using a new diagrammatic technique, the CC equations are represented in a compact form due to the string notations we introduced. As an application of these ideas, a new automated implementation of the single-reference-based multi-reference CC equations is presented for arbitrary excitation levels. The new program can be considered as an improvement over the previous implementations in many respects; e.g., diagram contributions are evaluated by efficient vectorized subroutines. Timings for test calculations for various complete active-space problems are presented. As an application of the new code, the weak interactions in the Be dimer were studied.

Hybrid coupled cluster methods: Combining active space coupled cluster methods with coupled cluster singles, doubles, and perturbative triples

Based on the coupled-cluster singles, doubles, and a hybrid treatment of triples (CCSD(T)-h) method developed by us [J. Shen, E. Xu, Z. Kou, and S. Li, J. Chem. Phys. 132, 114115 (2010); and ibid. 133, 234106 (2010); and ibid. 134, 044134 (2011)], we developed and implemented a new hybrid coupled cluster (CC) method, named CCSD(T)q-h, by combining CC singles and doubles, and active triples and quadruples (CCSDtq) with CCSD(T) to deal with the electronic structures of molecules with significant multireference character. These two hybrid CC methods can be solved with non-canonical and canonical MOs. With canonical MOs, the CCSD(T)-like equations in these two methods can be solved directly without iteration so that the storage of all triple excitation amplitudes can be avoided. A practical procedure to divide canonical MOs into active and inactive subsets is proposed. Numerical calculations demonstrated that CCSD(T)-h with canonical MOs can well reproduce the corresponding results obtained with non-canonical MOs. For three atom exchange reactions, we found that CCSD(T)-h can offer a significant improvement over the popular CCSD(T) method in describing the reaction barriers. For the bond-breaking processes in F(2) and H(2)O, our calculations demonstrated that CCSD(T)q-h is a good approximation to CCSDTQ over the entire bond dissociation processes.

Biorthogonal moment expansions in coupled-cluster theory: Review of key concepts and merging the renormalized and active-space coupled-cluster methods

After reviewing recent progress in the area of the development of coupled-cluster (CC) methods for quasi-degenerate electronic states that are characterized by stronger non-dynamical correlation effects, including new generations of single- and multi-reference approaches that can handle bond breaking and excited states dominated by many-electron transitions, and after discussing the key elements of the left-eigenstate completely renormalized (CR) CC and equation-of-motion (EOM) CC methods, and the underlying biorthogonal method of moments of CC (MMCC) equations [P. Piecuch, M. Włoch, J. Chem. Phys. 123 (2005) 224105; P. Piecuch, M. Włoch, J.R. Gour, A. Kinal, Chem. Phys. Lett. 418 (2006) 467; M. Włoch, M.D. Lodriguito, P. Piecuch, J.R. Gour, Mol. Phys. 104 (2006) 2149], it is argued that it is beneficial to merge the CR-CC/EOMCC and active-space CC/EOMCC [P. Piecuch, Mol. Phys. 108 (2010) 2987, and references therein] theories into a single formalism. In order to accomplish this goal, the biorthogonal MMCC theory, which provides compact many-body expansions for the differences between the full configuration interaction and CC or, in the case of excited states, EOMCC energies, obtained using conventional truncation schemes in the cluster operator T and excitation operator Rμ, is generalized, so that one can correct the CC/EOMCC energies obtained with arbitrary truncations in T and Rμ for the selected many-electron correlation effects of interest. The resulting moment expansions, defining the new, Flexible MMCC (Flex-MMCC) formalism, and the ensuing CC(P;Q) hierarchy, proposed in the present work, enable one to correct energies obtained in the active-space CC and EOMCC calculations, in which one selects higher many-body components of T and Rμ via active orbitals and which recover much of the relevant non-dynamical and some dynamical electron correlation effects in applications involving potential energy surfaces (PESs) along bond breaking coordinates, for the effects of higher-order, primarily dynamical, correlations missing in the active-space CC/EOMCC considerations. The Flex-MMCC corrections to the active-space CC/EOMCC energies are mathematically similar to the non-iterative energy corrections defining the existing left-eigenstate CR-CC and CR-EOMCC methods, such as CR-CC(2,3) and CR-EOMCC(2,3). The potential advantages of the Flex-MMCC and CC(P;Q) formalisms are illustrated by describing the initial implementation and numerical tests of the novel CC hybrid scheme, abbreviated as CC(t;3), in which one corrects the results of the CC calculations with singles, doubles, and active-space triples, termed CCSDt, for the remaining effects due to connected triple excitations that are missing in the CCSDt considerations, but are present in the MMCC-based CR-CC(2,3) approach. By examining bond breaking in the HF, F2, and F2+ molecules, it is demonstrated that the CC(t;3) method improves the CCSDt and CR-CC(2,3) results, providing PESs that agree with those obtained with the full CC theory with singles, doubles, and triples (CCSDT) to within small fractions of a millihartree, at the fraction of the computer costs of the CCSDT calculations. Different strategies for defining active-space triples within the CC(t;3) scheme and the underlying CCSDt method are discussed. When limited to the ground-state problem, the CC(t;3) approach can be regarded as an improved and rigorously derived extension of the recently proposed CCSD(T)-h method [J. Shen, E. Xu, Z. Kou, S. Li, J. Chem. Phys. 132 (2010) 114115], in which triples corrections of the CCSD(T) type are replaced by their more robust CR-CC(2,3)-style analogs.

Coupled-cluster with active space selected higher amplitudes: Performance of seminatural orbitals for ground and excited state calculations

The active space approach for coupled-cluster models is generalized using the general active space concept and implemented in a string-based general coupled-cluster code. Particular attention is devoted to the choice of orbitals on which the subspace division is based. Seminatural orbitals are proposed for that purpose. These orbitals are obtained by diagonalizing only the hole-hole and particle-particle block of the one-electron density of a lower-order method. The seminatural orbitals are shown to be a good replacement for complete active space self-consistent field orbitals and avoid the ambiguities with respect to the reference determinant introduced by the latter orbitals. The seminatural orbitals also perform well in excited state calculations, including excited states with strong double excitation contributions, which usually are difficult to describe with standard coupled-cluster methods. A set of vertical excitation energies is obtained and benchmarked against full configuration interaction calculations, and alternative hierarchies of active space coupled-cluster models are proposed. As a simple application the spectroscopic constants of the C(2) B (1)Delta(g) and B(') (1)Sigma(g) (+) states are calculated using active space coupled-cluster methods and basis sets up to quadruple-zeta quality in connection with extrapolation and additivity schemes.

On the Nature of Unrestricted Orbitals in Variational Active Space Wave Functions

Active space coupled cluster methods exhibit unusual, nonsmooth spin symmetry-breaking behavior where the unrestricted minimum lies higher in energy at short bond distances and crosses below the restricted solution at longer distances. The restricted solution is also observed to be a stable minimum slightly beyond the symmetry-breaking point. This behavior arises due to differences in the optimal active spaces defining the restricted and unrestricted wave functions and results in unrestricted wave functions that are not strictly size consistent. We suggest a new, size-consistent model that allows the orbitals to break spin symmetry only within the active space.

Active-space coupled-cluster methods through connected quadruple excitations

Coupled-cluster methods that include just a subset of all connected triple, quadruple, or both excitation amplitudes, according to the ansatz of and Adamowicz co-workers [Int. Rev. Phys. Chem. 12, 339 (1993); J. Chem. Phys. 99, 1875 (1993); 100, 5792 (1994)] and Piecuch et al. [J. Chem. Phys. 110, 6103 (1999)], have been implemented into parallel execution programs. They are applicable to closed- and open-shell species and they take advantage of real Abelian point-group symmetry. A symbol manipulation program has been invoked to automate the implementation. These methods have been applied to the singlet-triplet separations of five triatomic hydrides (CH2, NH2+, SiH2, PH2+, and AsH2+) with consideration of scalar relativistic effects. They have been shown to be remarkably effective with errors arising from the use of a very small subset of higher-order excitations being no more than a few tenths of 1 kcal/mol.

A comparison of the renormalized and active-space coupled-cluster methods: Potential energy curves of BH and F 2

The renormalized noniterative CCSD(T) and CCSD(TQ) methods, which are the examples of the method of moments of coupled-cluster (CC) equations, and the active-space CC approach with internal and semi-internal triexcited clusters (CCSDt) are applied to the potential energy curves of BH and F 2. It is shown that the renormalized CCSD(T) and CCSD(TQ) methods are practically as effective in removing the failing of the standard noniterative CC approaches at larger internuclear separations as the CCSDt method. The results of the CCSDt and renormalized CCSD(T) and CCSD(TQ) calculations are also compared with those obtained with other perturbative CC approaches that focus on bond breaking.