New results for the OH (nu = 0,j = 0) + CO (nu = 0,j = 0) --> H + CO2 reaction: Five- and full-dimensional quantum dynamical study on several potential energy surfaces.

Abstract

Full- [six-dimensional (6-D)] and reduced-dimensional [five-dimensional (5-D)] quantum wave packet calculations have been performed for the title reaction to obtain reaction probabilities deriving from the ground rovibrational states of OH and CO with total angular momentum J = 0. Three potential energy surfaces (PES) are studied, namely, those of Bradley and Schatz (BS), Yu, Muckerman, and Sears (YMS), and Lakin, Troya, Schatz, and Harding (LTSH). 6-D calculations are performed only for the BS PES, while 5-D results are reported for all three PES'. The 6-D results obtained in the present work improve on those previously reported, since a larger vibrational basis and a better representation of the OH and CO bonds has been introduced. In particular, we now employ a generalized Lanczos-Morse discrete variable representation for both the OH and CO vibrations. In a further improvement, the generalized discrete variable representation of the CO vibration is based on different CO intramolecular potentials for the asymptotic and product grids employed in our projection formalism. This new treatment of the vibrational bases allows for a large reduction in computation time with respect to our previous implementation of the wave packet method, for a given level of accuracy. As a result, we have been able to extend the range of collision energies for which we can obtain converged 6-D results to a higher energy (0.8 eV) than was possible before (0.5 eV). The comparison of the new 6-D and previous 5-D results for the BS PES shows good agreement of the general trend in the reaction probabilities over all collision energies considered (0.1-0.8 eV), while our previous 6-D calculation showed reaction probabilities that differed from the 5-D results by up to 10% between 0.5 and 0.8 eV. The 5-D reaction probabilities reveal interesting trends for the different PES'. In particular, at low energies (< 0.2 eV) the LTSH PES gives rise to much larger reactivity than the other PES', while at high energies (> 0.3 eV) its reaction probability decreases with respect to the BS and YMS PES', being more than a factor of 2 smaller at 0.8 eV. A 5-D calculation on a modified version of the LTSH surface shows that the van der Waals interaction in the entrance channel, which is not correctly described in the other PES' is largely responsible for its larger reactivity at low energies. The large difference between the 5-D reaction probabilities for the YMS and LTSH PES' serves to emphasize the importance of the van der Waals interaction for the reactivity at low energies, because most of the stationary point energies on the YMS and LTSH PES are rather similar, being in line with high-level ab initio information.