计算化学公社

标题: CINEB寻找过渡态问题求助 [打印本页]

作者
Author: wrtgcn 时间: 2019-7-12 17:35
标题: CINEB寻找过渡态问题求助
初学VASP，现计算一个气体分子在催化剂表面吸附后断键过程的过渡态。过程中遇到几个问题一直没能解决，真心向各位老师或同学求教，在此先谢谢大家了。

1. 插了4个images，结构由nebmake生成，粗算时跑出来的能量曲线趋势如图： (, 下载次数 Times of downloads: 34) 。怀疑是之前末态找的不准所致，于是用图中能量最低点提高精度后重新跑了优化，也确实得到了能量更低的新末态。重新插点跑CINEB后发现新生成的能量曲线还是这个趋势，于是再次重算末态。然而目前已经算了3次末态，每次确实能量都更低了，但是跑CINEB的能量趋势始终是还是如图。请问大家，为什么会出现这样的情况？应该如何解决？

2. 用VASP自带的优化器IBRION=1，POTIM=0.1跑CINEB粗算，实际上跑了不到100步四个images的力都已经很解决0.01的收敛限了，但是直到跑了1500+步后还是无法收敛。查看结果后发现各images的磁矩相差的很大。请问如果想通过固定磁矩帮助收敛，应该如何设置输入文件，又如何判断选择哪个磁矩才是合理的呢？

再次感谢各位了！

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