Pilot applications of internally contracted multireference coupled cluster theory, and how to choose the cluster operator properly

The internally contracted multireference coupled cluster (icMRCC) method allows a highly accurate description of both static and dynamic correlation with a computational scaling similar to single reference coupled cluster theory. The authors show that the method can lose its orbital invariance and size consistency when no special care is taken in the elimination of redundant excitations. Using the BeH(2) model system, four schemes are compared which differ in their treatment of linear dependencies between excitations of different rank (such as between singles and doubles). While the energy curves agree within tens of μE(h) when truncating the cluster operator at double excitations (icMRCCSD), inclusion of triple excitations (icMRCCSDT) leads to significant differences of more than 1 mE(h). One scheme clearly yields the best results, while the others even turn out to be not size consistent. The former procedure uses genuine single and double excitations and discards those linear combinations of (spectator) double and triple excitations which have the same effect on the reference function. With this approach, the equilibrium structure and harmonic vibrational frequencies of ozone obtained with icMRCCSDT are in excellent agreement with CCSDTQ. The authors further apply icMRCC methods to potential energy surfaces of HF, LiF, N(2), and to the singlet-triplet splitting of benzynes. In particular, the latter calculations have been made possible by implementing the method with the proper formal scaling using automated techniques.