Quantum dynamics studies of the <inline-formula><tex-math id="M2">\begin{document}$\rm D+SiD^+ \to D_2+Si^ +$\end{document}</tex-math><alternatives><graphic xmlns:xlink="http://www.w3.org/1999/xlink" xlink:href="22-20221155_M2.jpg"/><graphic xmlns:xlink="http://www.w3.org/1999/xlink" xlink:href="22-20221155_M2.png"/></alternatives></inline-formula> reaction
Authors: Zhao Wen-Li , Sun Feng-Wei , Zhang Hong , Wang Yong-Gang , Gao Feng , Meng Qing-Tian ,   School of Information Science and Engineering, Shandong Agricultural University, Taian 271018, China,   School of Physics and Electronics, Shandong Normal University, Jinan 250358, China
Published: August 2022

The quantum dynamics calculations are carried out for the title reaction D +SiD<sup>+</sup>→D<sub>2</sub>+Si<sup>+</sup> to obtain the initial (<inline-formula><tex-math id="M8">\begin{document}$ \nu = 0{\text{ }},j = 0 $\end{document}</tex-math><alternatives><graphic xmlns:xlink="http://www.w3.org/1999/xlink" xlink:href="22-20221155_M8.jpg"/><graphic xmlns:xlink="http://www.w3.org/1999/xlink" xlink:href="22-20221155_M8.png"/></alternatives></inline-formula>)reaction probability, integral cross section (ICS) and rate constant on the potential energy surface (PES) of Gao, Meng and Song. A total of 110 partial waves are calculated by using the Chebyshev wave packet method with full Coriolis coupling (CC) and centrifugal sudden (CS) approximation in a collision energy range from 1.0 × 10<sup>–3 </sup>to 1.0 eV. The calculated probability decreases with the collision energy increasing except for <i>J≤</i>40. The calculation results indicate that the CS approximation will overestimate or underestimate the reaction probability . The ICS decreases with the collision energy increasing and shows an oscillatory structure due to the<inline-formula><tex-math id="Z-20221114130407">\begin{document}$\rm{SiH_2^+} $\end{document}</tex-math><alternatives><graphic xmlns:xlink="http://www.w3.org/1999/xlink" xlink:href="22-20221155_Z-20221114130407.jpg"/><graphic xmlns:xlink="http://www.w3.org/1999/xlink" xlink:href="22-20221155_Z-20221114130407.png"/></alternatives></inline-formula>well on the reaction path. The results show that the neglect of the Coriolis coupling leads to the overestimation of the cross section and the rate constant. Besides, the discrepancy between the integral cross sections from the CC and CS calculations decreases clearly with collision energy increasing. Comparison with the corresponding results of H+CH<sup>+</sup> reaction indicates that isotope substitution reaction makes the cross section and the rate constant underestimated. The resulting integral reaction cross section displays less oscillatory structure, especially in the exact quantum calculation with the full Coriolis coupling effect taken into consideration. The kinetic isotope effect <inline-formula><tex-math id="Z-20221117061024-1">\begin{document}$(\kappa_{\rm H+SiH^+}(T)/\kappa_{\rm D+SiD^+}(T))$\end{document}</tex-math><alternatives><graphic xmlns:xlink="http://www.w3.org/1999/xlink" xlink:href="22-20221155_Z-20221117061024-1.jpg"/><graphic xmlns:xlink="http://www.w3.org/1999/xlink" xlink:href="22-20221155_Z-20221117061024-1.png"/></alternatives></inline-formula>is found to decrease with temperature increasing. It can be seen that the reduced mass of reactant can exert a certain effect on dynamic behavior.