Zhang Yue, Zhang Shaowen, Li Qian Shu
School of Science, Beijing Institute of Technology, Beijing 100081, People's Republic of China.
J Comput Chem. 2004 Jan 30;25(2):218-26. doi: 10.1002/jcc.10373.
We present a direct ab initio and hybrid density functional theory dynamics study of the thermal rate constants of the unimolecular decomposition reaction of C2H5O-->CH2O + CH3 at a high-pressure limit. MPW1K/6-31+G(d,p), MP2/6-31+G(d,p), and MP2(full)/6-31G(d) methods were employed to optimize the geometries of all stationary points and to calculate the minimum energy path (MEP). The energies of all the stationary points were refined at a series of multicoefficient and multilevel methods. Among all methods, the QCISD(T)/aug-cc-pVTZ energies are in good agreement with the available experimental data. The rate constants were evaluated based on the energetics from the QCISD(T)/aug-cc-pVTZ//MPW1K/6-31+G(d,p) level of theory using both microcanonical variational transition state theory (microVT) and RRKM theory with the Eckart tunneling correction in the temperature range of 300-2500 K. The calculated rate constants at the QCISD(T)/aug-cc-pVTZ/MPW1K/6-31+G(d,p) level of theory are in good consistent with experimental data. The fitted three-parameter Arrhenius expression from the microVT/Eckart rate constants in the temperature range 200-2500 K is k = 2.52 x 10(12)T(0.41)e(-8894.0/T) s(-1). The falloff curves of pressure-dependent rate constants are performed using master-equation method within the temperature range of 391-471 K. The calculated results are in good agreement with the available experimental data.
我们展示了一项关于C2H5O→CH2O + CH3单分子分解反应在高压极限下热速率常数的直接从头算和混合密度泛函理论动力学研究。采用MPW1K/6-31+G(d,p)、MP2/6-31+G(d,p)和MP2(full)/6-31G(d)方法优化所有驻点的几何结构并计算最小能量路径(MEP)。所有驻点的能量在一系列多系数和多水平方法下进行了精修。在所有方法中,QCISD(T)/aug-cc-pVTZ能量与现有实验数据吻合良好。基于QCISD(T)/aug-cc-pVTZ//MPW1K/6-31+G(d,p)理论水平的能量学,使用微正则变分过渡态理论(microVT)和带有埃卡特隧穿校正的RRKM理论,在300 - 2500 K的温度范围内评估速率常数。在QCISD(T)/aug-cc-pVTZ/MPW1K/6-31+G(d,p)理论水平下计算得到的速率常数与实验数据高度一致。在200 - 2500 K温度范围内,由microVT/埃卡特速率常数拟合得到的三参数阿伦尼乌斯表达式为k = 2.52 x 10(12)T(0.41)e(-8894.0/T) s(-1)。在391 - 471 K的温度范围内,使用主方程方法计算压力依赖速率常数的衰减曲线。计算结果与现有实验数据吻合良好。