计算化学公社

 找回密码 Forget password
 注册 Register

单线态氧跑过渡态不收敛

查看数: 448 | 评论数: 3 | 收藏 Add to favorites 1
关灯 | 提示:支持键盘翻页<-左 右->
    组图打开中,请稍候......
发布时间: 2025-3-24 13:55

正文摘要:

各位老师,在前一个帖子里面听了大神的方法,使用wb97xd跑单线态氧的过渡态,但是结果不收敛,请问一下各位大佬是什么原因呀,结构部分我是看文献中,单线态氧对苯环的进攻方式进行的初猜 %mem=1GB %nprocshared ...

回复 Reply

Uus/pMeC6H4-/キ 发表于 Post on 2025-3-24 22:49:31
本帖最后由 Uus/pMeC6H4-/キ 于 2025-3-25 00:50 编辑

作为练习跑了跑体系,先用小一点且Gaussian针对性优化过的6-311G(d,p)基组在同泛函下预搜索过渡态,再切换到def-TZVP正式优化并检验波函数稳定性,给12核24 GB内存的情况下2个小时足以顺利收敛完成计算(不包含IRC)。后两个任务的相关文件已放在压缩包里。
phenol_1O2_ts.tar.gz (1.46 MB, 下载次数 Times of downloads: 25)
最终单一虚频为412.19i cm-1且振动方向符合单重态氧与苯环发生一步[4+2]环加成的预期,单重态体系波函数也确认稳定。

虽然体系不算大,但既然要用def-TZVP这个档次的基组就得给够并行核数和内存,4核1 GB是肯定不够的,平均分给每核的内存太少了(http://sobereva.com/439)。“不收敛”应当区分SCF不收敛(http://sobereva.com/61)和几何优化不收敛(http://sobereva.com/164)来解决。写IOP(1/8)等价于opt里写maxstep,有关键词能用就不必写抽象的IOP了(http://sobereva.com/331)。

编辑:等下,为什么坐标里有机分子是邻苯二甲酸二甲酯但图片里是苯酚啊(尴尬)且不说Diels-Alder反应的区域选择性对不对,本来和苯环共轭的酯基为什么偏离平面这么多?
wzkchem5 发表于 Post on 2025-3-24 18:08:59
确定这个反应不是无垒反应吗?
确定这个反应不是按电子转移机理反应的吗?
wal 发表于 Post on 2025-3-24 15:35:35
震荡了吗?要是没震荡,初猜还是手摆的,还调小步长搜过渡态,那步数八成不够用

单线态氧用rwb97xd搜完过渡态建议做个波函数稳定性测试

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

GMT+8, 2025-8-15 20:39 , Processed in 0.388720 second(s), 25 queries , Gzip On.

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