计算化学公社

标题: OH提取甲醇氢原子，活性空间选取问题 [打印本页]

作者
Author: forpaper 时间: 2022-11-29 10:19
标题: OH提取甲醇氢原子，活性空间选取问题
本帖最后由 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以上的距离，结构如图所示[attach]59664[/attach]
文献中测试了(3,3) (5,5) (7,7) (9,9) (11,11)这五种活性空间，我选择重复(3,3)这个活性空间的结果。
下图是文献中的(3,3)活性空间包含的轨道
[attach]59668[/attach][attach]59667[/attach][attach]59669[/attach]下面是我用MOKIT计算得到的轨道(手动调整了轨道顺序）
[attach]59670[/attach][attach]59670[/attach][attach]59671[/attach]
文献中采取不同活性空间计算得到的能垒都在2-4 kcal/mol。而我计算得到的能垒为8 kcal/mol。所以想请教下各位老师，是不是我的活性空间选取不合理。

附件是我的反应物和过渡态的输入输出文件。

欢迎光临计算化学公社 (http://bbs.keinsci.com/)