MRD-CI potential surfaces using balanced basis sets. IV. The H2 molecule and the H3 surface

Abstract

The utility of midbond functions in molecular calculations was tested in two cases where the correct results are known: the H2 potential curve and the collinear H3 potential surface. For H2, a variety of basis sets both with and without bond functions was compared to the exact nonrelativistic potential curve of Kolos and Wolniewicz [J. Chem. Phys. 43, 2429 (1965)]. It was found that optimally balanced basis sets at two levels of quality were the double zeta single polarization plus sp bond function basis (BF1) and the triple zeta double polarization plus two sets of sp bond function basis (BF2). These gave bond dissociation energies De =4.7341 and 4.7368 eV, respectively (expt. 4.7477 eV). Four basis sets were tested for basis set superposition errors, which were found to be small relative to basis set incompleteness and therefore did not affect any conclusions regarding basis set balance. Basis sets BF1 and BF2 were used to construct potential surfaces for collinear H3, along with the corresponding basis sets DZ*P and TZ*PP which contain no bond functions. Barrier heights of 12.52, 10.37, 10.06, and 9.96 kcal/mol were obtained for basis sets DZ*P, TZ*PP, BF1, and BF2, respectively, compared to an estimated limiting value of 9.60 kcal/mol. Difference maps, force constants, and relative rms deviations show that the bond functions improve the surface shape as well as the barrier height.