No abstract

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 ̅Ĥ =e(- ̂T)Ĥe( ̂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 ̅Ĥ , 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.

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.

Application of multireference equation of motion coupled-cluster theory to transition metal complexes and an orbital selection scheme for the efficient calculation of excitation energies

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.