计算化学公社

 找回密码 Forget password
 注册 Register

求助:Cl-NEB计算出的过渡态没有虚频

查看数: 4951 | 评论数: 4 | 收藏 Add to favorites 1
关灯 | 提示:支持键盘翻页<-左 右->
    组图打开中,请稍候......
发布时间: 2022-5-23 00:18

正文摘要:

本帖最后由 1234567 于 2022-5-23 00:37 编辑 各位前辈,小弟最近在学习过渡态相关计算,在学习的过程中发现下述问题: Cl-NEB计算正常收敛,原子扩散路线和能量趋势变化也比较正常; 但是选择能量最高点(暂认 ...

回复 Reply

1234567 发表于 Post on 2022-5-23 19:56:12
Riszzz 发表于 2022-5-23 19:16
我是个小垃圾,就分享下之前的经历,你可以等等其他老师的回答。
你跑dimer的时候给MODECAR了吗,你这个 ...

哈哈,前辈太谦虚了,经验很宝贵的~特别感谢分享以及给的建议,我会一一尝试的。
其实我在算的时候给MODECAR了,是根据过渡态相邻两个images生成的,现在这个模型算的过渡态路线在x-y面大概就是一条直线的样子;
我之前算Cl-NEB一直认为设计扩散路径的初末态是全局能量最低点,过渡态是能量次低点来着,不知道这个理解是不是有问题,我再去找书看看,真的特别感谢,有前辈指导交流真好~祝前辈有一个开心的一天~

Riszzz 发表于 Post on 2022-5-23 19:16:16
本帖最后由 Riszzz 于 2022-5-23 19:20 编辑
1234567 发表于 2022-5-23 17:39
好的 谢谢前辈指点~
那我再试试补充一些之前没算过的吸附位的结构优化,看看有没有能量更低的点出现。
...

我是个小垃圾,就分享下之前的经历,你可以等等其他老师的回答。
你跑dimer的时候给MODECAR了吗,你这个像是没给的样子。不要有误区呀,不一定是能量最低的才是正确的初末态,这个就跟极小值和最小值的区别似的,极小值不一定是最小值。
他们算的多的人,有的时候都是直接摆结构跑dimer的,因为dimer比neb帧数少,快一些,但是相对不太稳定,可以通过先找到过渡态,再根据振动方向找初末态。

评分 Rate

参与人数
Participants 1
eV +2 收起 理由
Reason
1234567 + 2 谢谢分享

查看全部评分 View all ratings

1234567 发表于 Post on 2022-5-23 17:39:03
Riszzz 发表于 2022-5-23 16:34
你这个DIMCAR,曲率一直都是正的,显然不太对
这种情况我遇到过很多次,有些时候甚至dimer和ci-neb都能正 ...

好的 谢谢前辈指点~
那我再试试补充一些之前没算过的吸附位的结构优化,看看有没有能量更低的点出现。

前辈 方便再问一下这种情况下再找合适的初末态结构有比较高效的办法么?我看网上有大神帖子说算吸附势能面来确定能量最低点和次低点的,但是我还没试过,目前找初末态还是借助前面的吸附模型能量对比得来的。
再次感谢前辈:)
Riszzz 发表于 Post on 2022-5-23 16:34:45
你这个DIMCAR,曲率一直都是正的,显然不太对
这种情况我遇到过很多次,有些时候甚至dimer和ci-neb都能正常收敛,但是过渡态没有虚频,我也尝试提高收敛精度,但仍然没有。最后都是重新选择了初末态结构,才得到了正确的结果。

评分 Rate

参与人数
Participants 1
eV +3 收起 理由
Reason
1234567 + 3 谢谢

查看全部评分 View all ratings

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

GMT+8, 2026-2-20 10:21 , Processed in 0.189427 second(s), 32 queries , Gzip On.

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