Corchado Jose C, Bravo Jose L, Espinosa-Garcia Joaquin
Departamento de Química Física, Universidad de Extremadura, 06071 Badajoz, Spain.
J Chem Phys. 2009 May 14;130(18):184314. doi: 10.1063/1.3132223.
A new analytical potential energy surface is presented for the reaction of hydrogen abstraction from methane by a hydrogen atom. It is based on an analytical expression proposed by Jordan and Gilbert [J. Chem. Phys. 102, 5669 (1995)], and its fittable parameters were obtained by a multibeginning optimization procedure to reproduce high-level ab initio electronic structure calculations obtained at the CCSD(T)/cc-pVTZ level. The ab initio information employed in the fit includes properties (equilibrium geometries, relative energies, and vibrational frequencies) of the reactants, products, saddle point, points on the reaction path, and points on the reaction swath. No experimental information is used. By comparison with the reference results we show that the resulting surface reproduces well not only the ab initio data used in the fitting but also other thermochemical and kinetic results computed at the same ab initio level, such as equilibrium constants, rate constants, and kinetic isotope effects, which were not used in the fit. In this way we show that the new potential energy surface is correctly fitted and almost as accurate as the CCSD(T)/cc-pVTZ method in describing the kinetics of the reaction. We analyze the limitations of the functional form and the fitting method employed, and suggest some solutions to their drawbacks. In a forthcoming communication, we test the quality of the new surface by comparing its results with experimental values.
本文给出了一个新的氢原子从甲烷中夺取氢反应的解析势能面。它基于Jordan和Gilbert [J. Chem. Phys. 102, 5669 (1995)] 提出的一个解析表达式,其可拟合参数通过多起点优化程序获得,以重现CCSD(T)/cc-pVTZ水平下的高水平从头算电子结构计算结果。拟合中使用的从头算信息包括反应物、产物、鞍点、反应路径上的点以及反应带内的点的性质(平衡几何结构、相对能量和振动频率)。未使用实验信息。通过与参考结果比较,我们表明所得势能面不仅能很好地重现拟合中使用的从头算数据,还能重现相同从头算水平下计算的其他热化学和动力学结果,如平衡常数、速率常数和动力学同位素效应,这些结果在拟合中未被使用。通过这种方式,我们表明新的势能面拟合正确,在描述反应动力学方面几乎与CCSD(T)/cc-pVTZ方法一样准确。我们分析了所采用的函数形式和拟合方法的局限性,并针对其缺点提出了一些解决方案。在即将发表的一篇通讯中,我们将通过将新势能面的结果与实验值进行比较来检验其质量。