[发明专利]计算机计算芳烃加氢反应活化能及反应速率常数的方法

权利要求书

1.一种利用计算机计算不同温度下的芳烃加氢反应活化能及反应速率常数的方法，该方法包括以下步骤：

1)利用MS软件的Visualizer模块画出芳烃加氢反应的反应物分子式和生成物分子式，得到与所述反应物分子式和生成物分子式对应的反应物结构数据文件和生成物结构数据文件；

2)利用MS软件的Forcite模块中的Geometry Optimization功能对所述反应物结构数据文件和生成物结构数据文件进行分子力学能量最小化处理，得到能量最小反应物结构数据文件和能量最小生成物结构数据文件；

3)利用MS软件的DMol³模块中的Geometry Optimization功能对所述能量最小反应物结构数据文件和能量最小生成物结构数据文件进行量子力学结构优化处理，得到能量最小反应物构象文件和能量最小生成物构象文件，并且利用MS软件的DMol³模块中的Energy功能对所述能量最小反应物构象文件进行计算，得到不同温度下的反应物能量；

4)利用MS软件的Reaction Preview功能生成反应的轨迹文件，并利用MS软件的DMol³模块中的TS Search功能搜索反应过渡态构型，得到过渡态构型文件；

5)利用MS软件的DMol³模块中的Energy功能对所述过渡态构型文件进行计算，得到不同温度下的过渡态能量；

6)计算不同温度下的过渡态能量与对应温度下的反应物能量的差值，得到所述不同温度下的芳烃加氢反应活化能；以及

7)利用所述不同温度下的芳烃加氢反应活化能计算得到对应温度下的芳烃加氢反应速率常数。

2.根据权利要求1所述的方法，其中，步骤2)中对所述反应物结构数据文件和生成物结构数据文件进行分子力学能量最小化处理时，选用力场Universial或力场COMPASS或力场cvff。

3.根据权利要求1所述的方法，其中，下列任一种情况下，选用函数LDA-VWN或函数GGA-BLYP或函数RPBE，选用基组DNP或基组DND，将函数Charge取值为0-0.3之间的任意值：

步骤3)中对所述能量最小化的反应物结构数据文件和能量最小化的生成物结构数据文件进行量子力学结构优化处理时；

步骤3)中计算不同温度下的反应物能量时；

步骤4)中进行反应过渡态构型搜索时；或

步骤5)中计算不同温度下的过渡态能量时。