计算化学公社

标题: vasp寻找过渡态，先用CINEB粗收敛完成，DIMER高精度收敛不了 [打印本页]

作者
Author: kknllsx 时间: 2023-6-17 18:35
标题: vasp寻找过渡态，先用CINEB粗收敛完成，DIMER高精度收敛不了
先用CINEB将力收敛标准为0.2跑完了迁移路径，然后利用nebresults.pl 和neb2dim.pl 进行DIMER方法进行高精度过渡态的寻找，参数按照vaspkit所给的提示设置了，但是一直没有收敛，目前DIMER还在继续计算，我是按照下图的方法进行DIMER搜寻的，希望老师们能指点一下迷津，感谢感谢。

作者
Author: ahxb 时间: 2023-6-18 22:31
本帖最后由 ahxb 于 2023-6-18 22:43 编辑

Dimer方法计算过渡态的过程中，Curvature一列应该绝大多数为负，所以你这个计算肯定不行了，停了吧。
如果你在-NEB中只是用了较低的力收敛判据，但EDIFF和K点等能量计算并未打折扣，那么更合理的做法是不是用neb2dim.pl脚本，而是直接使用CI-NEB计算中能量最高点的结构作为Dimer计算的初始结构，MODECAR通过另外做一次频率计算（用低ENCUT+Gamma点都可以），选择与预想的过渡态对应的振动模式，然后用OUTCAR中的对应的振动模式准备MODECAR。原因有两点：第一是因为CI-NEB计算中过渡态对应的结构多数情况下收敛得比其它结构好，第二是neb2dim.pl大概是用提供的两个结构通过有限差分得到MODECAR，各image之间间距一般比较大，和做一次频率计算再得到MODECAR相比还是太粗糙了，初始结构不错的话，用不算准确的MODECAR太浪费了。不会手动操作的话，也可以用http://bbs.keinsci.com/thread-7493-1-1.html里的脚本。
另外DRotMax不宜过大，否则会在Dimer旋转上浪费过多的计算量。按我个人经验，初始结构还不错的话，DRotMax一般取3~5就可以。NSW取200也少了，特别是在Dimer方法里用IOPT=2选CG方法优化过渡态的时候。在DRotMax取3时，按我个人经验，一般得取到500才能较大概率在EDIFFG=-2E-2的力收敛判据下一次计算完毕。

作者
Author: kknllsx 时间: 2023-6-19 11:58

ahxb 发表于 2023-6-18 22:31
Dimer方法计算过渡态的过程中，Curvature一列应该绝大多数为负，所以你这个计算肯定不行了，停了吧。
如果 ...

好的，谢谢，我试试

作者
Author: qyw 时间: 2024-6-3 19:53

ahxb 发表于 2023-6-18 22:31
Dimer方法计算过渡态的过程中，Curvature一列应该绝大多数为负，所以你这个计算肯定不行了，停了吧。
如果 ...

请问这个MODECAR咋得到的，通过频率

作者
Author: vasp冲冲 时间: 2024-6-3 21:46

ahxb 发表于 2023-6-18 22:31
Dimer方法计算过渡态的过程中，Curvature一列应该绝大多数为负，所以你这个计算肯定不行了，停了吧。
如果 ...

牛

作者
Author: vasp冲冲 时间: 2024-6-3 21:47

ahxb 发表于 2023-6-18 22:31
Dimer方法计算过渡态的过程中，Curvature一列应该绝大多数为负，所以你这个计算肯定不行了，停了吧。
如果 ...

感觉大哥你说的和IDM方法相似了

作者
Author: 后到咸阳斯蒂庞 时间: 2024-9-29 17:06

qyw 发表于 2024-6-3 19:53
请问这个MODECAR咋得到的，通过频率

完成频率计算之后，在OUTCAR中选取对应虚频的后三列，详见白皛白大佬的：http://bbs.keinsci.com/forum.php ... mp;page=1#pid276642

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