计算化学公社

标题: 求助：Cl-NEB计算出的过渡态没有虚频 [打印本页]

作者
Author: 1234567 时间: 2022-5-23 00:18
标题: 求助：Cl-NEB计算出的过渡态没有虚频
本帖最后由 1234567 于 2022-5-23 00:37 编辑

各位前辈，小弟最近在学习过渡态相关计算，在学习的过程中发现下述问题：
Cl-NEB计算正常收敛，原子扩散路线和能量趋势变化也比较正常；
但是选择能量最高点（暂认为是过渡态）进行频率计算时，发现没有虚频出现；根据网上搜集信息，进行了增大Cl-NEB力收敛精度(EDIFFG增加到-0.01)的计算(此时EDIFF=1E-6)，然后再选择能量最高image计算频率(EDIFF=1E-6)，发现还是没有虚频；【注：此时检查neb.dat中切向力列，是接近0的，如附图】
(, 下载次数 Times of downloads: 28)
过程中考虑过结构选择的不合理性，比如选择的过渡态不是真正的一阶鞍点等，但是根据文献报道和个人吸附能计算结果（代表性的位点），选择的初始点和过渡态点基本跟文献一致，也与个人计算的吸附能结果比较符合；
然后就想着试试公社里前辈们分享的经验，选择dimer对Cl-NEB(EDIFFG=-0.01/EDIFF=1E-6)计算得到的能量最高点进行再次优化：但是，很遗憾DIMCAR显示如下，我就终止了计算;【INCAR中EDIFF=1e-7/EDIFFG=-0.001，设置精度比较高是因为与Cl-NEB一样的话第一个离子步就收敛了，还是没有虚频】

因为对dimer，Cl-NEB计算都是开始学习的阶段，实在有点拿不准最主要的问题是啥，现在就是真的怀疑是不是找的过渡态有问题，需要做势能面计算来确定能量到底怎么变化么？
请各位前辈多多指教，小弟这里先谢谢大家了。

作者
Author: Riszzz 时间: 2022-5-23 16:34
你这个DIMCAR，曲率一直都是正的，显然不太对
这种情况我遇到过很多次，有些时候甚至dimer和ci-neb都能正常收敛，但是过渡态没有虚频，我也尝试提高收敛精度，但仍然没有。最后都是重新选择了初末态结构，才得到了正确的结果。

作者
Author: 1234567 时间: 2022-5-23 17:39

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

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

前辈方便再问一下这种情况下再找合适的初末态结构有比较高效的办法么？我看网上有大神帖子说算吸附势能面来确定能量最低点和次低点的，但是我还没试过，目前找初末态还是借助前面的吸附模型能量对比得来的。
再次感谢前辈：）

作者
Author: Riszzz 时间: 2022-5-23 19:16
本帖最后由 Riszzz 于 2022-5-23 19:20 编辑

1234567 发表于 2022-5-23 17:39
好的谢谢前辈指点~
那我再试试补充一些之前没算过的吸附位的结构优化，看看有没有能量更低的点出现。
...

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

作者
Author: 1234567 时间: 2022-5-23 19:56

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

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

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