Partially linearized, fully size-extensive, and reduced multireference coupled-cluster methods. II. Applications and performance

Abstract

The partially linearized (pl), fully size-extensive multireference (MR) coupled-cluster (CC) method, fully accounting for singles (S) and doubles (D) and approximately for a subset of primary higher than doubles, referred to as plMR CCSD, as well as its plMR CCSD(T) version corrected for secondary triples, as described in Part I of this paper [X. Li and J. Paldus, J. Chem. Phys. 128, 144118 (2008)], are applied to the problem of bond breaking in the HF, F2, H2O, and N2 molecules, as well as to the H4 model, using basis sets of a DZ or a cc-pVDZ quality that enable a comparison with the full configuration interaction (FCI) exact energies for a given ab initio model. A comparison of the performance of the plMR CCSD/CCSD(T) approaches with those of the reduced MR (RMR) CCSD/CCSD(T) methods, as well as with the standard single reference (SR) CCSD and CCSD(T) methods, is made in each case. For the H4 model and N2 we also compare our results with the completely renormalized (CR) CC(2,3) method [P. Piecuch and M. Włoch, J. Chem. Phys. 123, 224105 (2005)]. An important role of a proper choice of the model space for the MR-type methods is also addressed. The advantages and shortcomings of all these methods are pointed out and discussed, as well as their size-extensivity characteristics, in which case we distinguish supersystems involving noninteracting SR and MR subsystems from those involving only MR-type subsystems. Although the plMR-type approaches render fully size-extensive results, while the RMR CCSD may slightly violate this property, the latter method yields invariably superior results to the plMR CCSD ones and is more easy to apply in highly demanding cases, such as the triple-bond breaking in the nitrogen molecule.