Reinisch Guillaume, Miki Kenji, Vignoles Gérard L, Wong Bryan M, Simmons Chris S
J Chem Theory Comput. 2012 Aug 14;8(8):2713-2724. doi: 10.1021/ct300278x. Epub 2012 Jun 12.
We propose a general approach to describe large amplitude motions (LAM) with multiple degrees of freedom (DOF) in molecules or reaction intermediates, which is useful for the computation of thermochemical or kinetic data. The kinetic part of the LAM Lagrangian is derived using a Z-matrix internal coordinate representation within a new numerical procedure. This derivation is exact for a classical system, and the uncertainties on the prediction of observable quantities largely arise from uncertainties on the LAM potential energy surface (PES) itself. In order to rigorously account for these uncertainties, we present an approach based on Bayesian theory to infer a parametrized physical model of the PES using ab initio calculations. This framework allows for quantification of uncertainties associated with a PES model as well as the forward propagation of these uncertainties to the quantity of interest. A selection and generalization of some treatments accounting for the coupling of the LAM with other internal or external DOF are also presented. Finally, we discuss and validate the approach with two applications: the calculation of the partition function of 1,3-butadiene and the calculation of the high-pressure reaction rate of the CH(3) + H → CH(4) recombination.
我们提出了一种通用方法来描述分子或反应中间体中具有多个自由度(DOF)的大幅度运动(LAM),这对于热化学或动力学数据的计算很有用。LAM拉格朗日量的动力学部分是在一种新的数值程序中使用Z矩阵内坐标表示法推导出来的。这种推导对于经典系统是精确的,并且可观测量预测的不确定性很大程度上源于LAM势能面(PES)本身的不确定性。为了严格考虑这些不确定性,我们提出了一种基于贝叶斯理论的方法,使用从头算计算来推断PES的参数化物理模型。该框架允许量化与PES模型相关的不确定性,以及将这些不确定性正向传播到感兴趣的量。还介绍了一些考虑LAM与其他内部或外部自由度耦合的处理方法的选择和推广。最后,我们通过两个应用来讨论和验证该方法:1,3 - 丁二烯配分函数的计算以及CH(3) + H → CH(4) 重组的高压反应速率的计算。