Exploring the MRCI method for calculating interaction energies: application to the HeNe potential curve

Abstract

A multi-reference configuration interaction (MRCI) method is described, which is devised for the calculation of interaction energies of van der Waals complexes and applied to calculating the HeNe potential energy curve. The MRCI calculations make use of a generalized Pople-correction in order to account for the lack of size consistency. The orbital space is partitioned into three subspaces: the first active space (AS1), which contains the strongly occupied orbitals; the second active space (AS2), which contains the main intra-correlating orbitals; and the external space (ES). It is shown that, to keep the error below ± 0.2 K in the excitation scheme and the active orbital space it is sufficient to include only σ-orbitals in AS2 and to use an excitation scheme (labelled Qq-MRCI) that encompasses only up to quadruply excited configurations. The final active orbital space (AS2) turned out to be 2s(He), 2pσ(He), 3s(Ne), 3σ(Ne) and 3dσ(Ne). Other MRCI variants, in which most or all quadruply excited configurations were deleted from the CI expansion (Qt- and Tt-MRCI), were found to be inadequate. Using the Qq-MRCI scheme together with a 197-orbital ‘interaction optimized’ basis set (IO197), the MRCI interaction energy at R = 5.7 a0 was calculated to be -21.12K. The corresponding values at the MP4 and CCSD(T) levels of theory are -20.06 K and -20.99 K, respectively, indicating that the MP4 method is inappropriate for highly accurate calculations on this system. Fitting the calculated data using a generalized Morse function, including an additional C6/R6 term to account for a correct long-range behaviour of the potential, the MRCI well depth was calculated to be -21.16K at Req = 5.73a 0. The MRCI and CCSD(T) potentials have the same quality and are found to be in good agreement with the Hartree-Fock dispersion (HFD-B) potential of Keil, M., Danielson, L. J., and Dunlop, P. J., 1991, J. Chem. Phys., 94, 296. It is concluded that, for basis IO197, the CCSD(T) method is sufficiently accurate for calculating the HeNe interaction. To recover the small, missing contributions (a few tenths of a Kelvin), MRCI should be used.