Theoretische Chemie, Fakultät für Chemie, Universität Bielefeld, Postfach 100131, D-33501 Bielefeld, Germany.
J Chem Phys. 2010 Feb 28;132(8):084103. doi: 10.1063/1.3304920.
Full-dimensional quantum dynamics calculations for the H+CH(4)-->H(2)+CH(3) reaction using curvilinear coordinates are presented. A curvilinear coordinate system to describe reactions of the type X+YCH(3)-->XY+CH(3) is developed which facilitates efficient calculations using the multiconfigurational time-dependent Hartree (MCTDH) approach. To describe the bending motion of the X and Y atoms relative to the axis defined by the CH(3) fragment, coordinates based on stereographic projection are introduced. These coordinates yield a kinetic energy operator free of singularities within the dynamically relevant region. Employing this curvilinear coordinate system, full-dimensional and reduced dimensional MCTDH calculations study the cumulative reaction probability (for J=0) and the thermal rate constant for the H+CH(4) reaction on the Jordan-Gilbert potential energy surface [J. Chem. Phys. 102, 5669 (1995)]. The full-dimensional results agree very well with previous full-dimensional MCTDH results which used transition state based normal coordinates. The results of our eight-dimensional (8D) calculations are in reasonable agreement with the full-dimensional ones. They deviate significantly from older 8D results of Zhang et al. [J. Chem. Phys. 127, 234213 (2007)] but agree well with more recent results from the same group.
本文使用曲线坐标系统对 H+CH(4)-->H(2)+CH(3) 反应进行了全维量子动力学计算。开发了一种曲线坐标系来描述 X+YCH(3)-->XY+CH(3) 类型的反应,这有助于使用多组态含时哈特里(MCTDH)方法进行高效计算。为了描述 X 和 Y 原子相对于由 CH(3)片段定义的轴的弯曲运动,引入了基于立体投影的坐标。这些坐标在动态相关区域内产生无奇点的动能算子。本文采用该曲线坐标系统,通过全维及降维 MCTDH 计算,研究了 Jordan-Gilbert 势能面上 H+CH(4)反应的累积反应概率(J=0)和热速率常数[J. Chem. Phys. 102, 5669 (1995)]。全维结果与先前使用过渡态基正则坐标的全维 MCTDH 结果非常吻合。八维(8D)计算结果与全维结果基本一致。它们与 Zhang 等人的旧 8D 结果有显著差异[J. Chem. Phys. 127, 234213 (2007)],但与同一组的最新结果吻合较好。