计算化学公社

标题: Vasp做CINEB计算寻找过渡态 [打印本页]

作者
Author: 小书童 时间: 2017-4-15 11:25
标题: Vasp做CINEB计算寻找过渡态
本帖最后由小书童于 2017-4-15 15:21 编辑

Vasp做CINEB计算，实际做的是一条反应路径，能量最高点位过渡态。对于计算CINEB有如下问题请教：

1.如何查看CINEB计算的走势。（1）nebef.pl可以查看各个Images的能量和力，但是并不能看到计算的走势如何。（2）进入到各个文件夹内看能量和力的走势图，但也只是这个Image的走势，并不是整体的收敛趋势，而且各个点的收敛趋势似乎并不一致。所以请教各位老师，如何查看CINEB计算的收敛趋势好坏？

2.假如我们可以查看收敛走势了，那么我们发现趋势不好，就应该去调整结构重新计算。vasp的CINEB根据始末结构进行插点，然后进行计算。如果我想调整路径初猜，应该是调整路径上的各个点的POSCAR，一个路径上怎么也3-5个Images,如果手动调整的话，实在是不太现实啊。请教各位老师，对于调整路径初猜有什么好的方法吗？

3.对于计算CINEB计算，各位老师觉得跑多少步收敛能接受？（当然这个同体系、初猜结构等等都有关）

4.vasp做结构优化的时候会指定POTIM这个优化的步长，在做CINEB的时候关闭vasp自带的优化求解方法，POTIM=0，设定IOPT选择vtstcode中的优化方法，但是似乎并未对优化的步长进行设定，请教各位老师，这个CINEB计算的时候的优化步长怎么设定呢？以及POTIM这个参数设定的是优化的步长，还是步长上限？

作者
Author: 卡开发发 时间: 2017-4-16 11:52
2、可以考虑试试nebmake.pl分段来插值，比如其中某个结构做调整，然后你可以对反应物和调整的结构做一次nebmake.pl，然后调整的结构和产物再做一次nebmake.pl，多余的结构手工删除以及文件的编号要手工调整。
4、POTIM会决定步长，但不是步长也不是步长上限，“ For IBRION=1,2 or 3, POTIM serves as a scaling constant for the forces.”vtstcode当中的优化算法应该有个MAXMOVE的参数，这个是用来限制步长上限的。

作者
Author: 小书童 时间: 2017-4-16 16:36

卡开发发发表于 2017-4-16 11:52
2、可以考虑试试nebmake.pl分段来插值，比如其中某个结构做调整，然后你可以对反应物和调整的结构做一次neb ...

多谢老师指点，第二个问题的答案很不错，减轻很多手动操作。您通常怎么看CINEB计算的收敛趋势呢？

作者
Author: 卡开发发 时间: 2017-4-16 17:09

小书童发表于 2017-4-16 16:36
多谢老师指点，第二个问题的答案很不错，减轻很多手动操作。您通常怎么看CINEB计算的收敛趋势呢？

不知道，我一般过渡态做得少，CI-NEB+VASP用的也就很少了，具体模拟方面没啥经验。

作者
Author: 小书童 时间: 2017-4-16 17:34

卡开发发发表于 2017-4-16 17:09
不知道，我一般过渡态做得少，CI-NEB+VASP用的也就很少了，具体模拟方面没啥经验。

嗯嗯，谢谢老师耐心的解答，赞

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