计算化学公社

标题: CI-NEB方法搜索过渡态的问题 [打印本页]

作者
Author: 美美pig 时间: 2018-11-16 20:24
标题: CI-NEB方法搜索过渡态的问题
大家好，    刚接触了一段时间VASP，目前在用CI-NEB算过渡态的时候遇到了一些问题，还望得到大家的知道和帮助！
   问题1：在搜索过渡态的时候出现了如下图1这种现象，00表示初始结构，01和02这两个image的能量比初始结构还低，是因为初始结构不是最稳定结构的原因吗？这种情况我是不是应该考虑去将01和02的CONTCAR变为POSCAR进行优化呢？
   问题2：在优化结构的时候我们一般用IBRION=2（共轭梯度算法）来优化，但是算过渡态的时候是用IBRION=3（最速下降算法）来进行优化，所以直接用初始的能量作为标准是不是会有一点误差呢？毕竟优化的方法不一样，能量也会有一点变化吧？
   问题3：在计算过渡态的时候一般很难达到收敛标准，是我的设置的问题吗？还是有更好的方法？希望得到大家的指导！

我的INCAR，如图2。
希望能得到大家的帮助！谢谢大家了！

作者
Author: 美美pig 时间: 2018-11-16 20:28
补充一下，图一中第三列为能量

作者
Author: 啦啦黑还黑 时间: 2018-11-17 00:42
问题1：有两种可能，(1)你的初态不是一个极小值，(2)在初态和02之间还存在一个过渡态。建议，把00和01这两个结构用严格的收敛方法重新优化一下。

问题2：你的题目描述里有好几个错，
1. IBRION=3不是最速下降算法，没记错应该是阻尼md方法。
2. IBRION=3和POTIM=0就关闭了VASP自带的优化算法，实际上用的是IOPT=1，vtst的LBFGS优化算法。
3. 现在这情况，肯定不能用初态的能量做标准，按照我问题1中说的先重做优化。

问题3：INCAR中的问题：EDIFF至少要在1E-6以下。用1E-4算出力非常粗糙，所以过渡态不收敛
过渡态计算收敛技巧很多。一个比较普适的方法就是先用CINEB把力收敛到0.2左右，然后用neb2dim.pl，把能量最高的点单独拿出来跑dimer。这样找过渡态成功率很高。

作者
Author: 美美pig 时间: 2018-11-17 17:43

啦啦黑还黑发表于 2018-11-17 00:42
问题1：有两种可能，(1)你的初态不是一个极小值，(2)在初态和02之间还存在一个过渡态。建议，把00和01这两 ...

首先，非常感谢您的指导！

您说的"用严格的收敛方法重新优化一下”指的是优化结构的方法吗，就是用IBRION=2?

作者
Author: steve 时间: 2020-1-3 12:32
本帖最后由 steve 于 2020-10-15 21:45 编辑

啦啦黑还黑发表于 2018-11-17 00:42
问题1：有两种可能，(1)你的初态不是一个极小值，(2)在初态和02之间还存在一个过渡态。建议，把00和01这两 ...

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