Vansteenkiste P, Van Neck D, Van Speybroeck V, Waroquier M
Laboratory of Theoretical Physics, Center for Molecular Modeling, Ghent University, Proeftuinstraat 86, B-9000 Gent, Belgium.
J Chem Phys. 2006 Jan 28;124(4):044314. doi: 10.1063/1.2161218.
Large-amplitude motions, particularly internal rotations, are known to affect substantially thermodynamic functions and rate constants of reactions in which flexible molecules are involved. Up to now all methods for computing the partition functions of these motions rely on the Pitzer approximation of more than 50 years ago, in which the large-amplitude motion is treated in complete independence of the other (vibrational) degrees of freedom. In this paper an extended hindered-rotor model (EHR) is developed in which the vibrational modes, treated harmonically, are correctly separated from the large-amplitude motion and in which relaxation effects (the changes in the kinetic-energy matrix and potential curvature) are taken into account as one moves along the large-amplitude path. The model also relies on a specific coordinate system in which the Coriolis terms vanish at all times in the Hamiltonian. In this way an increased level of consistency between the various internal modes is achieved, as compared with the more usual hindered-rotor (HR) description. The method is illustrated by calculating the entropies and heat capacities on 1,3-butadiene and 1-butene (with, respectively, one and two internal rotors) and the rate constant for the addition reaction of a vinyl radical to ethene. We also discuss various variants of the one-dimensional hindered-rotor scheme existing in the literature and its relation with the EHR model. It is argued why in most cases the HR approach is already quite successful.
已知大幅度运动,尤其是内旋转,会显著影响涉及柔性分子的反应的热力学函数和速率常数。到目前为止,所有计算这些运动配分函数的方法都依赖于50多年前的皮兹近似,其中大幅度运动被完全独立于其他(振动)自由度来处理。本文开发了一种扩展的受阻转子模型(EHR),其中将按简谐处理的振动模式与大幅度运动正确分离,并且在沿着大幅度路径移动时考虑了弛豫效应(动能矩阵和势曲率的变化)。该模型还依赖于一个特定的坐标系,在该坐标系中哈密顿量中的科里奥利项始终为零。通过这种方式,与更常见的受阻转子(HR)描述相比,各种内模式之间的一致性水平得到了提高。通过计算1,3 - 丁二烯和1 - 丁烯(分别具有一个和两个内转子)的熵和热容以及乙烯基自由基与乙烯加成反应的速率常数来说明该方法。我们还讨论了文献中存在的一维受阻转子方案的各种变体及其与EHR模型的关系。阐述了为什么在大多数情况下HR方法已经相当成功。