Bond functions, covalent potential curves, and the basis set superposition error

In the current practice of quantum chemistry, it is not clear whether corrections for basis set superposition errors should be applied to the calculation of potential energy curves, in order to improve agreement with experimental data. To examine this question, spectroscopic parameters derived from theoretical potential curves are reported for the homonuclear diatomics C2, N2, O2, and F2, using a configuration interaction method. Three different basis sets were used, including double zeta plus polarization, triple zeta plus double polarization, and double zeta polarization augmented by bond functions. The bond function basis sets, which were optimized in the preceding paper to obtain accurate dissociation energies, also gave the most accurate parameters. The potential curves were then corrected for basis set superposition error using the counterpoise correction, and the spectroscopic parameters were computed again. The BSSE‐corrected curves showed worse agreement with experiment for all properties than the original (uncorrected) curves. The reasons for this finding are discussed. In addition to the numerical results, some problems in the application of the BSSE correction to basis sets containing bond functions are shown. In particular, there is an overcounting of the lowering due to the bond functions, regardless of which type of correction is applied. Also, genuine BSSE affects cannot be separated from energy‐lowering effects due to basis set incompleteness, and we postulate that it is the latter which is strongly dominant in the calculation of covalent potential curves. Based on these arguments, two conclusions follow: (1) application of BSSE corrections to potential curves should not be routinely applied in situations where the bonding is strong, and (2) appropriate use of bond functions can lead to systematic improvement in the quality of potential curves.

[1]  James S. Wright,et al.  MRD‐CI potential surfaces using balanced basis sets. VI. Correlation of bond order with bond function composition for first‐row diatomic molecules , 1991 .

[2]  R. Gijbels,et al.  Combined bond‐polarization basis sets for accurate determination of dissociation energies. II. Basis set superposition error as a function of the parent basis set , 1989 .

[3]  R. Gijbels,et al.  Combined bond polarization function basis sets for accurate ab initio calculation of the dissociation energies of AHn molecules (A=LI to F) , 1989 .

[4]  James S. Wright,et al.  MRD CI potential surfaces using balanced basis sets. V. Second-row diatomic hydrides , 1988 .

[5]  James S. Wright.  Accuracy of diatomic potential functions , 1988 .

[6]  James S. Wright,et al.  Approaching the exact energy for H2: Bond functions vs polarization functions , 1987 .

[7]  E. Kruus,et al.  MRD‐CI potential surfaces using balanced basis sets. IV. The H2 molecule and the H3 surface , 1986 .

[8]  David Feller,et al.  Basis Set Selection for Molecular Calculations , 1986 .

[9]  C. E. Dykstra,et al.  Improved counterpoise corrections for the abinitio calculation of hydrogen bonding interactions , 1986 .

[10]  R. Phillips,et al.  Implementation of the table CI method: Matrix elements between configurations with the same number of open-shells , 1985 .

[11]  Donald G. Truhlar,et al.  Systematic study of basis set superposition errors in the calculated interaction energy of two HF molecules , 1985 .

[12]  Richard J. Williams,et al.  MRD‐CI potential surfaces using balanced basis sets. I. First‐row diatomic hydrides , 1983 .

[13]  G. Herzberg,et al.  Constants of diatomic molecules , 1979 .

[14]  A. C. Roach,et al.  Ion-molecule reactions of the rare gases with hydrogen. Part 1.—Diatomics-in-molecules potential energy surface for ArH+2 , 1972 .

[15]  S. F. Boys,et al.  The calculation of small molecular interactions by the differences of separate total energies. Some procedures with reduced errors , 1970 .