计算化学公社

标题: vasp过渡态计算求助，cineb在接近收敛的时候收敛曲线会回弹，然后切应力变的特别大 [打印本页]

作者
Author: chy914681393 时间: 2024-9-16 10:38
标题: vasp过渡态计算求助，cineb在接近收敛的时候收敛曲线会回弹，然后切应力变的特别大
我是刚入门vasp的初学者，最近在算过渡态相关内容，有些问题不懂，因此想请教一下各位大佬。我开始是用nebmake做表面到体相内的线性插值，然后进行cineb计算，结果就是在接近收敛的时候收敛曲线会回弹，然后切应力变的特别大，我考虑可能是因为线性插值结构不合理导致，随后想着使用idpp进行非线性插点，但使用过程中发现循环到1000步的时候依然没有新点生成，这是什么原因呢

作者
Author: kirSTY 时间: 2024-9-16 15:24
看看回弹附近的构象有没有原子距离过近的问题，或者看看初始构象是不是合理吧，最好是手动猜一个过渡态，然后两边分别插点。另外插点可以试下社长的sobNEB那个小程序

作者
Author: 我爱学习d 时间: 2024-9-19 11:36
可以试试改一下算法我之前试IOPT=7的时候也会出现你这种情况，把IOPT换成3就收敛了

作者
Author: chy914681393 时间: 2024-9-19 13:00

我爱学习d 发表于 2024-9-19 11:36
可以试试改一下算法我之前试IOPT=7的时候也会出现你这种情况，把IOPT换成3就收敛了

好的我试试我之前是跳过粗算过程直接用IOPT=2精算，结果不收敛，然后又改成7想着来粗算一下的谢谢大佬

作者
Author: chy914681393 时间: 2024-9-19 13:03

kirSTY 发表于 2024-9-16 15:24
看看回弹附近的构象有没有原子距离过近的问题，或者看看初始构象是不是合理吧，最好是手动猜一个过渡态，然 ...

是的我感觉也是线性插点导致初始构象不合理导致的后来用makeneb的idpp进行非线性插点现在重新跑上了我之前用CASTEP搜索到了过渡态能否直接用CASTEP的过渡态结果来做插值计算呢

作者
Author: 我爱学习d 时间: 2024-9-23 21:26

chy914681393 发表于 2024-9-19 13:00
好的我试试我之前是跳过粗算过程直接用IOPT=2精算，结果不收敛，然后又改成7想着来粗算一下的谢谢大佬

我看官网也是推荐IOPT=7 结果发现非常不好用

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