The potential energy function of Cr 2 (X 1 AE þ g ) is computed using the MRCI þ Q (multi-reference configuration interaction with Davidson correction) and CASPT2 (complete active space second-order multi-reference perturbation theory) methods and very large basis sets. In addition, a new hybrid method denoted CIPT2, in which excitations solely from the active space are treated by MRCI and the remaining ones by perturbation theory, is proposed and applied. These methods are used to analyse valence and core correlation effects on the potential energy function. MRCI calculations which fully include correlation of the 3p, 3d, and 4s electrons yield excellent agreement with experimental values for the equilibrium distance R e and the harmonic vibrational constant ! e . Correlation of the (3p) orbitals is found to be very important. At the MRCI þ Q level this reduces R e by 0.08 Å and increases ! e by about 200 cm À1 . These effects are not correctly reproduced in CASPT2 and CIPT2. It is shown that in CASPT2 the valence correlation effects are significantly overestimated and the distance dependence of the core correlation is incorrect. The basis set limit of the dissociation energies for MRCI þ Q and CASPT 2 is estimated to be 1.3 eV and 1.9 eV, respectively. The remaining error (% 0:2 À 0:3 eV) of the MRCI þ Q dissociation energy is attributed to unlinked cluster and higher excitation effects, which are only approximately accounted for by the Davidson correction.
IntroductionThe calculation of the dissociation energy and spectroscopic properties of the Cr 2 dimer is a severe test for ab initio electron correlation methods. Its (formally) hextuple 1 AE þ g bond dissociates along the ground-state potential curve to two high-spin Cr atoms with d 5 s 1 ( 7 S) configuration. This leads to very large differential electron correlation effects. In addition, near the equilibrium distance there is a subtle competition between bonding of the 4s and 3d electrons, and repulsion of the 3p-electrons, which makes it necessary to correlate all electrons in the 3p, 3d, and 4s orbitals. To accurately account for these effects has proven to be an intricate computational task. Many different theoretical methods have been used in the past in attempts to reproduce the spectroscopic quantities experimentally measured: generalized valence bond [1, 2], density functional theory [3][4][5][6][7][8][9][10][11][12], coupled cluster methods [3], multi-reference perturbation theory [13][14][15][16][17][18][19][20][21][22], multi-reference configuration