Benchmark calculations with correlated molecular wave functions. IV. The classical barrier height of the H+H2→H2+H reaction

Using systematic sequences of correlation consistent Gaussian basis sets from double to sextuple zeta quality, the classical barrier height of the H+H2 exchange reaction has been calculated by multireference configuration interaction (MRCI) methods. The MRCI calculations for collinear H3 have also been calibrated against large basis set full CI (FCI) results, which demonstrate that the MRCI treatment leads to energies less than 1 μhartree (≤0.001 kcal/mol) above the FCI energies. The dependence of both the H2 and H3 total energies on the basis set is found to be very regular, and this behavior has been used to extrapolate to the complete basis set (CBS) limits. The resulting estimate of the H–H–H CBS limit yields a classical barrier height, relative to exact H+H2, of 9.60±0.02 kcal/mol; the best directly calculated value for the barrier is equal to 9.62 kcal/mol. These results are in excellent agreement with recent quantum Monte Carlo calculations.

[1]  Peter J. Knowles,et al.  A new determinant-based full configuration interaction method , 1984 .

[2]  D L Diedrich,et al.  An Accurate Quantum Monte Carlo Calculation of the Barrier Height for the Reaction H + H2 → H2 + H , 1992, Science.

[3]  Rick A. Kendall,et al.  BENCHMARK CALCULATIONS WITH CORRELATED MOLECULAR WAVE FUNCTIONS. III: CONFIGURATION INTERACTION CALCULATIONS ON FIRST ROW HOMONUCLEAR DIATOMICS , 1993 .

[4]  E. Clementi,et al.  Selected Topics in ab initio Computational Chemistry in Both Very Small and Very Large Chemical Systems , 1991 .

[5]  H. Monkhorst,et al.  New Born–Oppenheimer potential energy curve and vibrational energies for the electronic ground state of the hydrogen molecule , 1986 .

[6]  R. Raffenetti,et al.  General contraction of Gaussian atomic orbitals: Core, valence, polarization, and diffuse basis sets; Molecular integral evaluation , 1973 .

[7]  M. Blomberg,et al.  The H3 potential surface revisited , 1985 .

[8]  E. Levin,et al.  Ab initio potential energy surface for H–H2 , 1993 .

[9]  J. L. Dunham.  The Energy Levels of a Rotating Vibrator , 1932 .

[10]  Rick A. Kendall,et al.  Benchmark calculations with correlated molecular wave functions. II. Configuration interaction calculations on first row diatomic hydrides , 1993 .

[11]  Thom H. Dunning,et al.  Benchmark calculations with correlated molecular wave functions. I: Multireference configuration interaction calculations for the second row diatomic hydrides , 1993 .

[12]  P. Knowles,et al.  An efficient internally contracted multiconfiguration–reference configuration interaction method , 1988 .

[13]  C. Bauschlicher,et al.  A reevaluation of the H3 potential , 1990 .

[14]  F. B. Brown,et al.  A new semi-empirical method of correcting large-scale configuration interaction calculations for incomplete dynamic correlation of electrons , 1985 .

[15]  Bin Liu,et al.  An accurate three‐dimensional potential energy surface for H3 , 1978 .

[16]  David M. Ceperley,et al.  Quantum Monte Carlo for molecules: Green’s function and nodal release , 1984 .

[17]  C. Horowitz,et al.  Functional representation of Liu and Siegbahn’s accurate ab initio potential energy calculations for H+H2 , 1978 .

[18]  Bowen Liu,et al.  Ab initio potential energy surface for linear H3 , 1973 .

[19]  J. Rychlewski,et al.  Many‐electron explicitly correlated Gaussian functions. I. General theory and test results , 1993 .

[20]  P. Knowles,et al.  An efficient method for the evaluation of coupling coefficients in configuration interaction calculations , 1988 .

[21]  Bowen Liu.  Classical barrier height for H+H2→H2+H , 1984 .

[22]  M. R. Peterson,et al.  An improved H3 potential energy surface , 1991 .

[23]  P. Knowles,et al.  An efficient second-order MC SCF method for long configuration expansions , 1985 .

[24]  P. Tipler,et al.  Modern Physics , 1976 .

[25]  T. H. Dunning.  Gaussian basis sets for use in correlated molecular calculations. I. The atoms boron through neon and hydrogen , 1989 .

[26]  Donald G. Truhlar,et al.  A double many‐body expansion of the two lowest‐energy potential surfaces and nonadiabatic coupling for H3 , 1987 .

[27]  R. Bartlett,et al.  Basis set quantum chemistry and quantum Monte Carlo: Selected atomic and molecular results , 1992 .

[28]  Peter J. Knowles,et al.  A determinant based full configuration interaction program , 1989 .

[29]  T. Dunning,et al.  Electron affinities of the first‐row atoms revisited. Systematic basis sets and wave functions , 1992 .

[30]  P. Knowles,et al.  A second order multiconfiguration SCF procedure with optimum convergence , 1985 .

[31]  E. Clementi,et al.  The Hylleraas‐CI method in molecular calculations. III. Implementation and numerical verification of a three‐electron many‐center theory , 1991 .

[32]  P. Reynolds,et al.  H + H2 reaction barrier: A fixed‐node quantum Monte Carlo study , 1985 .