计算化学公社

 找回密码 Forget password
 注册 Register
Views: 23320|回复 Reply: 4
打印 Print 上一主题 Last thread 下一主题 Next thread

[VASP] Vasp做CINEB计算寻找过渡态

[复制链接 Copy URL]

407

帖子

0

威望

856

eV
积分
1263

Level 4 (黑子)

跳转到指定楼层 Go to specific reply
楼主
本帖最后由 小书童 于 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这个参数设定的是优化的步长,还是步长上限?





3753

帖子

3

威望

1万

eV
积分
19670

Level 6 (一方通行)

围观吃瓜群众

2#
发表于 Post on 2017-4-16 11:52:33 | 只看该作者 Only view this author
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的参数,这个是用来限制步长上限的。

评分 Rate

参与人数
Participants 2
eV +7 收起 理由
Reason
小书童 + 5 谢谢,学习了。
sobereva + 2

查看全部评分 View all ratings

日常打哑谜&&探寻更多可能。
原理问题不公开讨论,非商业性质讨论欢迎私聊。不做培*,不接代*,不接*发谢谢。

407

帖子

0

威望

856

eV
积分
1263

Level 4 (黑子)

3#
 楼主 Author| 发表于 Post on 2017-4-16 16:36:02 | 只看该作者 Only view this author
卡开发发 发表于 2017-4-16 11:52
2、可以考虑试试nebmake.pl分段来插值,比如其中某个结构做调整,然后你可以对反应物和调整的结构做一次neb ...

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

3753

帖子

3

威望

1万

eV
积分
19670

Level 6 (一方通行)

围观吃瓜群众

4#
发表于 Post on 2017-4-16 17:09:44 | 只看该作者 Only view this author
小书童 发表于 2017-4-16 16:36
多谢老师指点,第二个问题的答案很不错,减轻很多手动操作。您通常怎么看CINEB计算的收敛趋势呢?

不知道,我一般过渡态做得少,CI-NEB+VASP用的也就很少了,具体模拟方面没啥经验。
日常打哑谜&&探寻更多可能。
原理问题不公开讨论,非商业性质讨论欢迎私聊。不做培*,不接代*,不接*发谢谢。

407

帖子

0

威望

856

eV
积分
1263

Level 4 (黑子)

5#
 楼主 Author| 发表于 Post on 2017-4-16 17:34:34 | 只看该作者 Only view this author
卡开发发 发表于 2017-4-16 17:09
不知道,我一般过渡态做得少,CI-NEB+VASP用的也就很少了,具体模拟方面没啥经验。

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

本版积分规则 Credits rule

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2025-8-15 10:44 , Processed in 0.294908 second(s), 21 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list