Communication: multireference equation of motion coupled cluster: a transform and diagonalize approach to electronic structure.

The novel multireference equation-of-motion coupled-cluster (MREOM-CC) approaches provide versatile and accurate access to a large number of electronic states. The methods proceed by a sequence of many-body similarity transformations and a subsequent diagonalization of the transformed Hamiltonian over a compact subspace. The transformed Hamiltonian is a connected entity and preserves spin- and spatial symmetry properties of the original Hamiltonian, but is no longer Hermitean. The final diagonalization spaces are defined in terms of a complete active space (CAS) and limited excitations (1h, 1p, 2h, …) out of the CAS. The methods are invariant to rotations of orbitals within their respective subspaces (inactive, active, external). Applications to first row transition metal atoms (Cr, Mn, and Fe) are presented yielding results for up to 524 electronic states (for Cr) with an rms error compared to experiment of about 0.05 eV. The accuracy of the MREOM family of methods is closely related to its favorable extensivity properties as illustrated by calculations on the O2-O2 dimer. The computational costs of the transformation steps in MREOM are comparable to those of closed-shell Coupled Cluster Singles and Doubles (CCSD) approach.

SF-[2]R12: a spin-adapted explicitly correlated method applicable to arbitrary electronic states.

The [2](R12) method [M. Torheyden and E. F. Valeev, J. Chem. Phys. 131, 171103 (2009)] is an explicitly correlated perturbative correction that can greatly reduce the basis set error of an arbitrary electronic structure method for which the two-electron density matrix is available. Here we present a spin-adapted variant (denoted as SF-[2](R12)) that is formulated completely in terms of spin-free quantities. A spin-free cumulant decomposition and multi-reference generalized Brillouin condition are used to avoid three-particle reduced density matrix completely. The computational complexity of SF-[2](R12) is proportional to the sixth power of the system size and is comparable to the cost of the single-reference MP2-R12 method. The SF-[2](R12) method is shown to decrease greatly the basis set error of multi-configurational wave functions.

A novel interpretation of reduced density matrix and cumulant for electronic structure theories.

We propose a novel interpretation of the reduced density matrix (RDM) and its cumulant that combines linear and exponential parametrizations of the wavefunction. Any n-particle RDM can be written as a weighted average of "configuration interaction" amplitudes. The corresponding n-particle cumulant is represented in terms of two types of contributions: "connected" (statistical averages of substitution amplitudes) and "disconnected" (cross-correlations of substitution amplitudes). A diagonal element of n-RDM represents the average occupation number of the orbital n-tuple. The diagonal elements of 2- and 3-cumulants take particularly elegant forms in the natural spin-orbital basis: they represent the covariances (correlated fluctuations) of the occupation numbers of the orbital pair and triples, respectively. Thus, the diagonal elements of the cumulants quantify the correlation between the orbital occupation numbers. Our interpretation is used to examine the weak to strong correlation transition in the "two electrons in two orbitals" problem.

A state-specific partially internally contracted multireference coupled cluster approach.

A state-specific partially internally contracted multireference coupled cluster approach is presented for general complete active spaces with arbitrary number of active electrons. The dominant dynamical correlation is included via an exponential parametrization of internally contracted cluster operators ( ̂T) which excite electrons from a multideterminantal reference function. The remaining dynamical correlation and relaxation effects are included via a diagonalization of the transformed Hamiltonian ̅H =e(- ̂T)He( ̂T) in the multireference configuration interaction singles space in an uncontracted fashion. A new set of residual equations for determining the internally contracted cluster amplitudes is proposed. The second quantized matrix elements of ̅H , expressed using the extended normal ordering of Kutzelnigg and Mukherjee, are used as the residual equations without projection onto the excited configurations. These residual equations, referred to as the many-body residuals, do not have any near-singularity and thus, should allow one to solve all the amplitudes without discarding any. There are some relatively minor remaining convergence issues that may arise from an attempt to solve all the amplitudes and an initial analysis is provided in this paper. Applications to the bond-stretching potential energy surfaces for N(2), CO, and the low-lying electronic states of C(2) indicate clear improvements of the results using the many-body residuals over the conventional projected residual equations.

Perturbative correction for the basis set incompleteness error of complete-active-space self-consistent field.

To reduce the basis set incompleteness of the complete-active-space self-consistent field (CASSCF) wave function and energy we develop a second-order perturbation correction due to single excitations to complete set of unoccupied states. Other than the one- and two-electron integrals, only one- and two-particle reduced density matrices are required to compute the correction, denoted as [2](S). Benchmark calculations on prototypical ground-state bond-breaking problems show that only the aug-cc-pVXZ basis is needed with the [2](S) correction to match the accuracy of CASSCF energies of the aug-cc-pV(X+1)Z quality.

An algebraic proof of generalized Wick theorem.

The multireference normal order theory, introduced by Kutzelnigg and Mukherjee [J. Chem. Phys. 107, 432 (1997)], is defined explicitly, and an algebraic proof is given for the corresponding contraction rules for a product of any two normal ordered operators. The proof does not require that the contractions be cumulants, so it is less restricted. In addition, it follows from the proof that the normal order theory and corresponding contraction rules hold equally well if the contractions are only defined up to a certain level. These relaxations enable us to extend the original normal order theory. As a particular example, a quasi-normal-order theory is developed, in which only one-body contractions are present. These contractions are based on the one-particle reduced density matrix.

State specific equation of motion coupled cluster method in general active space.

The state specific equation of motion coupled cluster (SS-EOMCC) method is an internally contracted multireference approach, applicable to both ground and excited states. Attractive features of the method are as follows: (1) the SS-EOMCC wave function is qualitatively correct and rigorously spin adapted, (2) both orbitals and dynamical correlation are optimized for the target state, (3) nondynamical correlation and differential orbital relaxation effects are taken care of by a diagonalization of the transformed Hamiltonian in the multireference configuration-interaction singles space, (4) only one- and two-particle density matrices of a complete-active-space self-consistent-field reference state are needed to define equations for the cluster amplitudes, and (5) the method is invariant with respect to orbital rotations in core, active, and virtual subspaces. Prior applications focused on biradical-like systems, in which only one extra orbital is needed to construct the active space, and similarly, single bond breaking processes. In this paper, the applicability of the method is extended to systems of general active spaces. Studies on F(2), H(2)O, CO, and N(2) are carried out to gauge its accuracy. The convergence strategy is discussed in detail.

Application of double ionization state-specific equation of motion coupled cluster method to organic diradicals.

The state-specific equation of motion coupled cluster method is applied to three systems of diradical character: automerization of cyclobutadiene, singlet-triplet gaps of trimethylmethylene, and Bergman reaction. The aim of the paper is to assess the performance of the method and test numerically the importance of orbital optimization, three-body terms in transformed Hamiltonian, and the choice of cluster equations.