Yao Yuxuan, Lakshmanan Sandhiya, Pratihar Subha, Hase William L
Department of Chemistry and Biochemistry, Texas Tech University, Lubbock, Texas 79409, United States.
CSIR - National Institute of Science, Technology and Development Studies, New Delhi 110060, India.
J Phys Chem A. 2020 Mar 5;124(9):1821-1828. doi: 10.1021/acs.jpca.9b11513. Epub 2020 Feb 20.
The CH + O reaction has a quite complex ground state singlet potential energy surface (PES). There are multiple minima and transition states before forming the 10 possible reaction products. A previous direct chemical dynamics simulation at the UM06/6-311++G(d,p) level of theory ( 2019, 123, 4360-4369) found that reaction on this PES is predominantly direct without trapping in the potential minima. The first minima CH + O encounters is that for the CHO Criegee intermediate and statistical theory assumes the reactive system is trapped in this intermediate with a lifetime given by Rice-Ramsperger-Kassel-Marcus (RRKM) theory. In the work presented here, a direct dynamics simulation is performed with the above UM06 theory, with the trajectories initialized in the CHO intermediate with a random distribution of vibrational energy as assumed by RRKM theory. There are substantial differences between the dynamics for CHO dissociation and CH + O reaction. For the former there are four product channels, while for the latter there are seven in agreement with experiment. Product energy partitioning for the two simulations are in overall good agreement for the CO + H and CO + HO product channels, but in significant disagreement for the HCO + OH product channel. Though CHO is excited randomly in accord with RRKM theory, its dissociation probability is biexponential and not exponential as assumed by RRKM. In addition, the CHO dissociation dynamics follow non-intrinsic reaction coordinate (non-IRC) pathways. An important finding is that the nonstatistical dynamics for the CH + O reaction give results in agreement with experiment.
CH与O的反应具有相当复杂的基态单重态势能面(PES)。在形成10种可能的反应产物之前存在多个极小值和过渡态。先前在UM06/6 - 311++G(d,p)理论水平下进行的直接化学动力学模拟(J. Phys. Chem. A 2019, 123, 4360 - 4369)发现,该势能面上的反应主要是直接反应,不会被困在势能极小值处。CH与O遇到的第一个极小值是CHO克里奇中间体的极小值,统计理论假设反应体系被困在这个中间体中,其寿命由赖斯 - 拉姆齐 - 卡斯尔 - 马库斯(RRKM)理论给出。在本文所展示的工作中,使用上述UM06理论进行了直接动力学模拟,轨迹按照RRKM理论假设在CHO中间体中以振动能量的随机分布进行初始化。CHO解离动力学和CH与O反应动力学之间存在显著差异。对于前者有四个产物通道,而对于后者有七个产物通道,与实验结果一致。对于CO + H和CO + HO产物通道,两次模拟的产物能量分配总体上吻合良好,但对于HCO + OH产物通道则存在显著差异。尽管CHO按照RRKM理论被随机激发,但其解离概率是双指数的,而不是RRKM理论假设的指数形式。此外,CHO解离动力学遵循非本征反应坐标(non - IRC)路径。一个重要的发现是,CH与O反应的非统计动力学给出的结果与实验相符。