24小时热门版块排行榜    

CyRhmU.jpeg
查看: 6982  |  回复: 16
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

lsh1986914

新虫 (小有名气)


[交流] 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,重新投任务呢?求高人指点迷津~我这个低精度的可行吗?大家在算的时候,怎么设置这个精度滴?
回复此楼

» 猜你喜欢

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

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

查看全部散金贴

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

lsh1986914

新虫 (小有名气)


引用回帖:
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 !!!
14楼2012-06-01 15:15:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 17 个回答

gfei2008

木虫 (小有名气)


★ ★ ★
lsh1986914(金币+1): 谢谢参与
lsh1986914: 金币+2 2012-05-23 09:52:41
个人觉得在你第一次计算的精度就可以了。后面增加后最多相差0.1ev左右(可能,不确定)。要是就是用来看migration energy的话就够精度了。
2楼2012-05-22 22:42:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lsh1986914

新虫 (小有名气)


引用回帖:
2楼: Originally posted by gfei2008 at 2012-05-22 22:42:08:
个人觉得在你第一次计算的精度就可以了。后面增加后最多相差0.1ev左右(可能,不确定)。要是就是用来看migration energy的话就够精度了。

嗯,谢谢,我个人也觉得应该可以了,本身这个方法就不是很精确,以前算过,如果精度取得越高,势垒也就越低,但是相差并不是很大~
3楼2012-05-23 09:52:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

qfwuli

金虫 (正式写手)



lsh1986914(金币+1): 谢谢参与
最好能把初末态精度和C-NEB
4楼2012-05-23 10:33:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 高级回复(可上传附件)
信息提示
请填处理意见