计算化学公社

 找回密码 Forget password
 注册 Register
Views: 1700|回复 Reply: 0
打印 Print 上一主题 Last thread 下一主题 Next thread

[新手求助] OH提取甲醇氢原子,活性空间选取问题

[复制链接 Copy URL]

13

帖子

0

威望

1372

eV
积分
1385

Level 4 (黑子)

跳转到指定楼层 Go to specific reply
楼主
本帖最后由 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: 5

反应物输入文件

TS.zip

1.5 MB, 下载次数 Times of downloads: 6

过渡态输入文件

本版积分规则 Credits rule

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2024-11-25 13:41 , Processed in 0.242695 second(s), 23 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list