应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » NEB计算扩散路径的力收敛精度设置问题
171/1返回列表
查看: 7338  |  回复: 16
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

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

» 猜你喜欢

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

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

查看全部散金贴
1楼 2012-05-22 16:08:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

huangguiyang

新虫 (小有名气)
★ ★ ★
小木虫: 金币+0.5, 给个红包，谢谢回帖
贺仪: 金币+2, 谢谢指教！ 2012-05-31 08:34:33
力的精度应该要-0.01.至少也要-0.02。EDIFF=1E-5就可以。
-0.5绝对是不行的。势垒会改变很大。

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

肯定要kill了。

提交之前，应用可视化软件检查一下路径是否合理。而且路径不能太过复杂。如果太过复杂的话，或许还要你自己去调整一下中间镜像的结构，让初始路径合理一点。

复杂的路径，采用默认的优化算法会比较稳定。

只要最简单的路径，用CI-NEB的其他优化算法才会速度更快。例如空位扩散，这个就可以用CI-NEB的其他优化算法来算。
赞一下(2人) 回复此楼
10楼2012-05-30 23:14:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

huangguiyang

新虫 (小有名气)
★ ★
小木虫: 金币+0.5, 给个红包，谢谢回帖
贺仪: 金币+1, 谢谢指教！ 2012-05-31 08:34:43
引用回帖:
7楼: Originally posted by lsh1986914 at 2012-05-23 20:12:47
NSW=600
IBRION=3 ; POTIM=0 ; IOPT=3;
ISIF=4
EDIFFG=-0.1  ; LMAXMIX=4; IALGO=38

IMAGES=4; SPRING=-5; ICHAIN=0; LCLIMB=.T.; LTANGENTOLD=.F.; LDNEB=.F.
NEBCELL=.T.
这个是我的INCAR中关于NEB计算 ...

EDIFF  = 1E-5
   EDIFFG = -0.01
  LCLIMB=.TRUE.
   IBRION = 3
   POTIM  = 0.1    time-step for ion-motion
   ISIF=2
   IMAGES=4
不要改变cell的形状，改变cell形状的neb，理论上就有问题。正式版的CI-NEB中没有这部分的。
除非路径简单，否则默认的优化算法最稳定和快速。
避免使用IOPT=3;
删去“ SPRING=-5; ICHAIN=0; LCLIMB=.T.; LTANGENTOLD=.F.; LDNEB=.F .NEBCELL=.T.”
赞一下(1人) 回复此楼
11楼2012-05-30 23:21:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

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的回帖

qfwuli

金虫 (正式写手)
★ ★
franch: 金币+2, 鼓励交流，， 呵呵 2012-05-23 13:23:30
最好能把初末态精度和NEB中的精度设为一致，这样最终处理数据时加上初末态数值的图形会比较一致。
（刚才没有写完就误发送了，嘻嘻，无需给金币了）
赞一下 回复此楼
5楼2012-05-23 10:34:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

happyweb

禁虫 (初入文坛)
★ ★ ★
lsh1986914(金币+1): 谢谢参与
franch: 金币+2, 谢谢回帖交流 2012-05-23 19:43:56
本帖内容被屏蔽
6楼2012-05-23 17:36:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lsh1986914

新虫 (小有名气)
引用回帖:
6楼: Originally posted by happyweb at 2012-05-23 17:36:20:
60eV跑远了吧,改一下优化算法吧,IOPT=3感觉比较正常,另外精度0.1就够了,反正还要算频率来检验

NSW=600
IBRION=3 ; POTIM=0 ; IOPT=3;
ISIF=4
EDIFFG=-0.1  ; LMAXMIX=4; IALGO=38

IMAGES=4; SPRING=-5; ICHAIN=0; LCLIMB=.T.; LTANGENTOLD=.F.; LDNEB=.F.
NEBCELL=.T.
这个是我的INCAR中关于NEB计算的相关参数，这个NEB网页上推荐的比较容易收敛的方法，请帮我看看，我的参数设置对不？尤其是那几个没有打开的参数，现在都还不知道什么意思。另外我做的是表面物理吸附，杂质在不同位置吸附时vasp弛豫都稳定，只是想看看这两个态之间的势垒高度是多少？
赞一下 回复此楼
7楼2012-05-23 20:12:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lsh1986914

新虫 (小有名气)
引用回帖:
5楼: Originally posted by qfwuli at 2012-05-23 10:34:55:
最好能把初末态精度和NEB中的精度设为一致，这样最终处理数据时加上初末态数值的图形会比较一致。
（刚才没有写完就误发送了，嘻嘻，无需给金币了）

那就是不收敛怎么办？
赞一下 回复此楼
8楼2012-05-23 20:15:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

suntag

金虫 (初入文坛)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
你为什么要改变cellshape（ISIF=4）？这样收敛很慢的的，收敛标准的EDIFFG=-0.5或者-0.1实在有点宽。。。
赞一下 回复此楼
9楼2012-05-27 04:51:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lsh1986914

新虫 (小有名气)
但是，我的初始态和末态的晶格常数都变了，NEBCELL=.T.不打开，不是就是错的吗？ 我的计算就还是要看这种模型的结构，我做的是物理吸附，路径不是很复杂，这样不可以吗？
赞一下 回复此楼
12楼2012-06-01 15:00:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lsh1986914

新虫 (小有名气)
引用回帖:
9楼: Originally posted by suntag at 2012-05-27 04:51:34
你为什么要改变cellshape（ISIF=4）？这样收敛很慢的的，收敛标准的EDIFFG=-0.5或者-0.1实在有点宽。。。

我还要计算这种模型的结构，其中就是要看他稳定后的晶格常数，所以我的初始和末态优化的时候，都改变了它的晶格，我想这里肯定也要设置晶格可变的，不是这样吗？
赞一下 回复此楼
13楼2012-06-01 15:04:29
已阅   回复此楼   关注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的回帖

huangguiyang

新虫 (小有名气)
引用回帖:
14楼: Originally posted by lsh1986914 at 2012-06-01 15:15:55
哦，谢谢了，我改变一下，参数设置。还有一些虫友说NEBCELL=.F. 只能这样设置，不支持.T.，是这样吗？由于我的初末态在做优化的时候，都放开了晶格长ISIF=4； 然后按照我的理解，把NEB计算的时候，应该也要放开晶格 ...

不要放开，保持一致才行。
赞一下 回复此楼
15楼2012-06-01 22:25:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

浪心99

新虫 (小有名气)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
你好，我在运行脚本vfin.pl的时候，会出现Error: No [Output Directory] given.
请问是怎么回事啊？
赞一下 回复此楼
16楼2016-07-05 11:05:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

抖抖豆豆

新虫 (小有名气)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
引用回帖:
16楼: Originally posted by 浪心99 at 2016-07-05 11:05:20
你好，我在运行脚本vfin.pl的时候，会出现Error: No  given.
请问是怎么回事啊？

你使用这个命令,在后面要给它一个输出路径的
赞一下 回复此楼
17楼2016-07-06 20:25:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 lsh1986914 的主题更新
171/1返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 0856材料与化工，270求调剂 +11 YXCT 2026-03-01 13/650 2026-03-02 21:38 by sunny81
[基金申请] 成果系统访问量大，请15分钟后再尝试。由此给您造成的不便，敬请谅解。 +6 xhuama 2026-03-02 7/350 2026-03-02 21:23 by kingkocxr
[考研] 材料类考研调剂 +3 gemmgemm 2026-03-01 4/200 2026-03-02 20:02 by hypershenger
[考研] 课题组接收材料类调剂研究生 +4 gaoxiaoniuma 2026-02-28 5/250 2026-03-02 19:33 by Xinyuu
[考研] 求调剂 +10 yunziaaaaa 2026-03-01 11/550 2026-03-02 19:17 by caszguilin
[考研] 289求调剂 +8 yang婷 2026-03-02 9/450 2026-03-02 19:08 by zhukairuo
[考研] 化工京区271求调剂 +6 11ing 2026-03-02 6/300 2026-03-02 18:52 by caszguilin
[考研] 一志愿华中科技大学，化学专业344分，求调剂 +3 邢xing1 2026-03-02 3/150 2026-03-02 17:32 by houyaoxu
[考研] 欢迎采矿、地质、岩土、计算机、人工智能等专业的同学报考 +6 pin8023 2026-02-28 8/400 2026-03-02 17:13 by 0854蹲调剂
[考研] 一志愿华南理工大学材料与化工326分，求调剂 +3 wujinrui1 2026-02-28 3/150 2026-03-02 16:36 by chuocheng
[考研] 303求调剂 +5 今夏不夏 2026-03-01 5/250 2026-03-02 15:01 by 向上的胖东
[考研] 化工专硕348，一志愿985求调剂 +6 弗格个 2026-02-28 9/450 2026-03-02 14:09 by liyongv
[考研] 292求调剂 +7 yhk_819 2026-02-28 7/350 2026-03-02 12:43 by 无际的草原
[考研] 281求调剂 +5 2026计算机_诚心 2026-03-01 8/400 2026-03-02 11:05 by 汪！？！
[基金申请] 成果系统访问量大，请一小时后再尝试。---NSFC啥时候好哦，已经两天这样了 +4 NSFC2026我来了 2026-02-28 4/200 2026-03-01 22:37 by 铁门栓
[考博] 26申博 +4 想申博！ 2026-02-26 6/300 2026-03-01 17:32 by 想申博！
[考研] 328求调剂 +3 aaadim 2026-03-01 5/250 2026-03-01 17:29 by njzyff
[考研] 313求调剂 +3 水流年lc 2026-02-28 3/150 2026-03-01 16:01 by 新能源达人
[考研] 调剂 +3 简木ChuFront 2026-02-28 3/150 2026-03-01 11:46 by 王伟要上岸啊
[硕博家园] 【博士招生】太原理工大学2026化工博士 +4 N1ce_try 2026-02-24 8/400 2026-02-26 08:40 by N1ce_try
信息提示
关闭
请填处理意见
关闭