|
|
本帖最后由 forpaper 于 2022-11-29 10:19 编辑
各位老师好,学生在学习多参考计算时,找了一篇文献(J. Am. Chem. Soc. 2018, 140, 2906−2918)重复他的多参考计算工作。反应为OH自由基提取甲醇氧上的氢原子,作者采用CASPT2/MG3S//M08-HX/MG3S方法计算的这个反应的能垒。
因为MG3S基组对于CHO原子的定义等价于6-311+G(2df,2p)。我采取CASPT2/6-311+G(2df,2p)方法计算的单点能。
过渡态的结构直接从支撑材料中提取的
反应物是把甲醇和OH自由基分隔在10A以上的距离,结构如图所示
文献中测试了(3,3) (5,5) (7,7) (9,9) (11,11)这五种活性空间,我选择重复(3,3)这个活性空间的结果。
下图是文献中的(3,3)活性空间包含的轨道
下面是我用MOKIT计算得到的轨道(手动调整了轨道顺序)
文献中采取不同活性空间计算得到的能垒都在2-4 kcal/mol。而我计算得到的能垒为8 kcal/mol。所以想请教下各位老师,是不是我的活性空间选取不合理。
附件是我的反应物和过渡态的输入输出文件。
|
-
-
Reactants.zip
911.29 KB, 下载次数 Times of downloads: 7
反应物输入文件
-
-
TS.zip
1.5 MB, 下载次数 Times of downloads: 8
过渡态输入文件
|