화학공학소재연구정보센터
Journal of Physical Chemistry, Vol.98, No.38, 9488-9497, 1994
Ab-Initio Potential-Energy Surface for H+ocs Reactions - Extended Basis-Sets and Correlation Treatment
Ab initio calculations are presented for the potential energy surface (PES) of H + OCS. The two triple-zeta AO basis sets used are the 6-311+G(2df,2p) and aug-cc-pVTZ. The highest levels of correlation employed are MP4 and QCISD(T) starting with UHF zeroth-order wave functions. There are two major reaction channels on the PES : reaction I is H(S-2) + OCS((1) Sigma) --> OH((II)-I-2) + CS((1) Sigma), and reaction II is H(S-2) + OCS((1) Sigma) --> SH((II)-I-2) + CO((1) Sigma). Multiple pathways leading from reactants to (I) and (II) were determined. Results of this study substantiate findings from an earlier quantum chemical study using a lower level of theory, including (1) the existence of 12 transition states and six stable four-body intermediates; (2) the qualitative description of the PES; i.e., geometries, relative barriers, and well depths are similar to those in the earlier study; and (3) the entrance channel transition states leading to (II) are tight, as suggested by experiment. The results presented here also support the explanations of observed product energy distributions for (I) and (II) based on the earlier ab initio study. An additional transition state connecting the cis-HOCS and cis-HSCO minima was located, confirming a previous suggestion that reaction II could result from hydrogen migration after HOCS formation. The size of the barrier to hydrogen migration on this PES, however, indicates that it is a high-energy reaction. The low-energy pathway leading to (II) is through direct formation of HSCO. The current results show a substantial improvement in the quantitative agreement with experiment over the previously calculated values. Our best calculations predict reaction enthalpies for (I) and (II) to be 57.3 and -11.5 kcal/mol, respectively, well within the uncertainty of the experimental data. The measured activation energy for (II), which includes all dynamic and quantum effects such as tunneling, is 3.9 kcal/mol, Our best theoretical value for the lowest energy barrier leading to (II) is 5.7 kcal/mol (including zero-point corrections). We emphasize that corrections due to tunneling are not included in this value.