We have formulated and implemented the multireference Mukherjee's coupled cluster method with connected singles, doubles, and perturbative triples [MR MkCCSD(T)] in the ACES II program package. Assessment of the new method has been performed on the first three electronic states of the oxygen molecule and on the automerization barrier of cyclobutadiene, where a comparison with other multireference CC treatments and with experimental data where available. The MR MkCCSD(T) method seems to be a promising candidate for an accurate, yet computationally tractable, treatment of systems where the static correlation plays an important role.

Extensions of multireference equation of motion coupled cluster theory (MR-EOMCC) [D. Datta and M. Nooijen, J. Chem. Phys. 137, 204107 (2012)] are presented that include additional correlation effects into the global, internally contracted similarity transformation, induced by the cluster operators. As a result the final uncontracted diagonalization space can be more compact than in the parent MR-EOMCC approach. A wide range of applications, including transition metal atomic excitation spectra, a large set of valence excited states of organic compounds, and potential energy surfaces of ground and excited states of butadiene, is presented to benchmark the applicability of the parent MR-EOMCC methodology and its new variations.

This paper reports the development of a local variant of Mukherjee's state-specific multireference coupled cluster method based on the pair natural orbital approach (LPNO-MkCC). The current implementation is restricted to single and double excitations. The performance of the LPNO-MkCCSD method was tested on calculations of naphthyne isomers, tetramethyleneethane, and β-carotene molecules. The results show that 99.7-99.8% of correlation energy was recovered with respect to the MkCC method based on canonical orbitals. Moreover, the errors of relative energies between different isomers or along a potential energy curve (with respect to the canonical method) are below 0.4 kcal/mol, safely within the chemical accuracy. The computational efficiency of our implementation of LPNO-MkCCSD in the ORCA program allows calculation of the β-carotene molecule (96 atoms and 1984 basis functions) on a single CPU core.

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 multireference equation of motion coupled-cluster theory to transition metal complexes and an orbital selection scheme for the efficient calculation of excitation energies