State-to-State Time-Dependent Quantum Dynamics Studies of the Si(³P) + OH(X ²Π) → OSi(X ¹Σ_g⁺) + H(²S) Reaction Based...

Quantum and quasi-classical dynamics calculations were conducted for the reaction of Si with OH on the latest potential energy surface (PES), which is obtained by fitting tens of thousands of ab initio energy points by using the many-body expansion formula. To obtain an accurate PES, all energy points calculated with aug-cc-pVQZ and aug-cc-pV5Z basis sets were extrapolated to the complete basis set limit. The accuracy of our new PES was verified by comparing the topographic characteristics and contour maps of potential energy with other works. In addition, the anharmonic vibrational frequencies of HOSi and HSiO based on the present ab initio and PES by means of quantum dynamics methods were calculated. Dynamics information such as reaction probability, integral cross sections (ICS), product distribution, and rate constants was obtained on the new HOSi(X²A') PES. The dynamic information calculated using the quasi-classical trajectory method and time-dependent wave packet method is generally in good agreement, except for the vibrational state-resolved ICSs of product. The calculated differential cross section and capture time reveal that the reaction is primarily governed by the complex formation mechanism.