24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

lsh1986914

新虫 (小有名气)

应助: 1 (幼儿园)
金币: 90.6
帖子: 68
在线: 18.1小时
虫号: 1145561

[交流] NEB计算扩散路径的力收敛精度设置问题

各位大侠们：
大家好！
我用NEB的C-NEB计算扩散路径，开始的时候收敛精度较小,EDIFF=1E-4; 而力收敛精度EDIFFFG=-0.5; 然后把扩散路径算出来，显示势垒高度达到了1.3eV，但是我个人觉得这个精度有点低了，因为我的初末态弛豫的时候都是采用EDIFF=1E-6; EDIFFG=-0.01; 。所以为了提高精度，运行脚本vfin.pl ,把所有的CONTCAR 拷贝为POSCAR，设置EDIFFG=1E-5; EDIFFG=-0.1; 结果跑了几天，离子步已经有了100多步了，可是还是不收敛，然后查看一下扩散路径，发现扩散路径都变了啊，势垒也有60eV了，这是怎么回事啊，是不是有问题呢？要不要kill,重新投任务呢？求高人指点迷津~我这个低精度的可行吗？大家在算的时候，怎么设置这个精度滴？

回复此楼

» 猜你喜欢

不自信的我已经有5人回复
磺酰氟产物，毕不了业了！已经有4人回复
论文终于录用啦！满足毕业条件了已经有16人回复
投稿Elsevier的杂志（返修），总是在选择OA和subscription界面被踢皮球已经有8人回复

» 本主题相关价值贴推荐，对您同样有帮助:

NEB计算过渡态和反应路径已经有9人回复
计算neb出错，请高手帮忙看看已经有27人回复
新手请教用NEB方法计算energy barrier. 已经有15人回复
CI-neb 收敛问题求教已经有15人回复
帮我看看 NEB的INCAR POSCAR 为什么总掉？已经有14人回复
cl-neb结果不收敛已经有19人回复
NEB过渡态搜索遇到的问题已经有33人回复
NEB计算的的收敛问题已经有6人回复
【文献学习】JACS上一篇纯DFT计算的文献： CI-NEB计算，Leuwdin布局分析，态密度分析已经有30人回复
【其他】vasp neb 计算讨论已经有20人回复
【求助成功】怎么样提高NEB计算速度? 已经有19人回复
【求助结束】结构优化的时候如何提高力的收敛精度？已经有10人回复

» 抢金币啦！回帖就可以得到:

查看全部散金贴

1楼 2012-05-22 16:08:54

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

lsh1986914

新虫 (小有名气)

应助: 1 (幼儿园)
金币: 90.6
帖子: 68
在线: 18.1小时
虫号: 1145561

引用回帖:

10楼: Originally posted by huangguiyang at 2012-05-30 23:14:21
力的精度应该要-0.01.至少也要-0.02。EDIFF=1E-5就可以。
-0.5绝对是不行的。势垒会改变很大。

至于扩散路径变了，有很多原因。需要具体分析。可以把POTIM减小，例如POTIM=0.01，或者0.05. 到力较小的时候，再改 ...

哦，谢谢了，我改变一下，参数设置。还有一些虫友说NEBCELL=.F. 只能这样设置，不支持.T.，是这样吗？由于我的初末态在做优化的时候，都放开了晶格长ISIF=4；然后按照我的理解，把NEB计算的时候，应该也要放开晶格，对吧。由于初末态的优化都放开晶格常数，导致了初末态的晶格常数不一致，然后运行nebmake.pl 脚本的时候，出现如下的警告：（是不是真的如有些人说的那样，不支持NEBCELL=.T.）
$ /home/myl/vasp/vtstscripts/nebmake.pl POSCAR1 POSCAR2 4
WARNING: BASIS VECTORS ARE NOT THE SAME
BASIS ELEMENT 0 0 ...
... IS IN FILE 1: 8.89385080815016
... IS IN FILE 2: 8.7371463018205
I HOPE YOU KNOW WHAT YOU ARE DOING

OK, ALL SETUP HERE
FOR LATER ANALYSIS, PUT OUTCARs IN FOLDERS 00 and 05 !!!