Quantum wave packet and quasiclassical trajectory studies of OH+CO: influence of the reactant channel well on thermal rate constants.

We study the OH+CO-->H+CO2 reaction with both six-dimensional quantum wave packets (QM) and quasiclassical trajectories (QCT), determining reaction probabilities and thermal rate constants (or coefficients), and studying the influence of the reactant channel hydrogen-bonded complex well on the reaction dynamics. The calculations use the recently developed Lakin-Troya-Schatz-Harding (LTSH) ground electronic state potential energy surface, along with a modified surface developed for this study (mod-LTSH), in which the reactant channel well is removed. Our results show that there can be significant differences between the QM and QCT descriptions of the reaction for ground-state reactants and for energies important to the thermal rate constants. Zero-point energy violation plays an important role in the QCT results, and as a result, the QCT reaction probability (for ground-state reactants and zero impact parameter) is much higher than its QM counterpart at moderate to low reagent translational energies. The influence of the reactant channel well in the QCT results is to enhance reactivity at moderate energies and to suppress reactivity at the very lowest collision energies. The QM results also show the enhancement at moderate energies but, while the very lowest translational energies cannot be adequately converged, they do not indicate any tendency toward suppression as energy is reduced. QCT calculations for excited rotational states of the reactants show that the suppression of reactivity associated with the reactant channel well is less important when the reactants are rotating, and as a result, the influence of the reactant channel well on the thermal rate coefficients is relatively small, being important below 200 K. Our results indicate that there still remain important discrepancies between experiment and theory in this low temperature regime and that further improvements of the potential are needed.

[1]  Trevor J. Sears,et al.  A theoretical study of the potential energy surface for the reaction OH+CO→H+CO2 , 2001 .

[2]  Michael Frenklach,et al.  OH(OD) + CO: Measurements and an Optimized RRKM Fit , 1998 .

[3]  David C. Clary,et al.  Potential optimized discrete variable representation , 1992 .

[4]  S. Gray,et al.  Quantum reactive scattering of O( 1 D) + H 2 and O( 1 D) + HD , 1998 .

[5]  G. Schatz,et al.  Quasiclassical Trajectory Study of Energy and Angular Distributions for the H + CO2 → OH + CO Reaction† , 2002 .

[6]  Robert J. Kee,et al.  Chemical Kinetics and Combustion Modeling , 1990 .

[7]  Bian,et al.  van der waals interactions in the Cl + HD reaction , 1999, Science.

[8]  G. Kroes,et al.  Erratum to: ‘Full-dimensional quantal initial state-selected reaction probabilities (J=0) for the reaction OH(v=0,j=0)+CO(v=0,j=0)→CO2+H’ [Chem. Phys. Lett. 352 (2002) 281–287] , 2003 .

[9]  Stephen K. Gray,et al.  Dispersion fitted finite difference method with applications to molecular quantum mechanics. , 2001 .

[10]  R. Levine,et al.  Molecular Reaction Dynamics and Chemical Reactivity , 1987 .

[11]  G. Schatz,et al.  A reduced dimension quantum reactive scattering study of OH + CO⇌H + CO2 , 1992 .

[12]  G. Schatz,et al.  Theoretical studies of polyatomic bimolecular reaction dynamics. , 1995, Annual review of physical chemistry.

[13]  David T. Anderson,et al.  Exploring the OH+CO reaction coordinate via infrared spectroscopy of the OH–CO reactant complex , 2000 .

[14]  G. Kroes,et al.  Converged five-dimensional quantum calculations for OH+CO→H+CO2 , 2002 .

[15]  Geert-Jan Kroes,et al.  Full-dimensional quantal initial state-selected reaction probabilities (J=0) for the reaction OH(v=0,j=0)+CO(v=0,j=0)→CO2+H , 2002 .

[16]  G. Schatz,et al.  A quasiclassical trajectory study of the OH+CO reaction , 1991 .

[17]  Ian W. M. Smith,et al.  Reaction between hydroxyl (deuteroxyl) radicals and carbon monoxide at temperatures down to 80 K: experiment and theory , 1993 .

[18]  G. Schatz,et al.  A quasiclassical trajectory study of OH rotational excitation in OH+CO collisions using ab initio potential surfaces , 1992 .

[19]  G. Schatz,et al.  A quasiclassical trajectory study of H+CO2: Angular and translational distributions, and OH angular momentum alignment , 1997 .

[20]  H. Hamann,et al.  High pressure range of addition reactions of HO. II. Temperature and pressure dependence of the reaction HO+CO⇔HOCO→H+CO2 , 1996 .

[21]  J. Cullum,et al.  Computing eigenvalues of very large symmetric matrices—An implementation of a Lanczos algorithm with no reorthogonalization , 1981 .

[22]  Richard G. Compton,et al.  Research in Chemical Kinetics , 1993 .

[23]  Ian W. M. Smith,et al.  Role of Hydrogen-Bonded Intermediates in the Bimolecular Reactions of the Hydroxyl Radical , 2002 .

[24]  R. Macdonald,et al.  A crossed‐beam study of the state‐resolved integral cross sections for the inelastic scattering of OH(X 2Π) with CO and N2 , 1991 .

[25]  A. Wagner,et al.  Mapping the OH + CO-->HOCO reaction pathway through IR spectroscopy of the OH-CO reactant complex. , 2001, Faraday discussions.

[26]  George C. Schatz,et al.  Quantum and quasiclassical calculations on the OH+CO→CO2+H reaction , 1993 .

[27]  John Z. H. Zhang,et al.  Quantum calculations of reaction probabilities for HO + CO→ H + CO2 and bound states of HOCO , 1995 .

[28]  Stephen K. Gray,et al.  A quantum dynamics study of H2+OH→H2O+H employing the Wu–Schatz–Lendvay–Fang–Harding potential function and a four-atom implementation of the real wave packet method , 2002 .

[29]  George C. Schatz,et al.  Quantum dynamics of a planar model for the complex forming OH+CO→H+CO2 reaction , 1995 .

[30]  Lawrence B. Harding,et al.  A quasiclassical trajectory study of the reaction OH+CO→H+CO2 , 2003 .

[31]  Rosendo Valero,et al.  Theoretical rate constants for the OH+CO→H+CO2 reaction using variational transition state theory on analytical potential energy surfaces , 2002 .

[32]  S. Gray,et al.  Quantum dynamics with real wave packets, including application to three-dimensional (J=0) D + H2 → HD + H reactive scattering , 1998 .

[33]  A. Ravishankara,et al.  Kinetic study of the reaction of OH with CO from 250 to 1040 K , 1983 .

[34]  Tucker Carrington,et al.  The discrete variable representation of a triatomic Hamiltonian in bond length–bond angle coordinates , 1992 .