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

24小时热门版块排行榜    

Znn3bq.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » NEB计算扩散路径的力收敛精度设置问题
171/1返回列表
查看: 7432  |  回复: 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返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 335求调剂 +10 想上岸呀！！ 2026-04-12 10/500 2026-04-12 09:04 by 314126402
[考研] 307求调剂 +10 tzq94092 2026-04-10 10/500 2026-04-12 08:18 by wise999
[考研] 一志愿新疆大学085401，314分 +4 咔咔咔咔9 2026-04-05 4/200 2026-04-12 02:21 by 秋豆菜芽
[考研] 070300化学279求调剂 +19 哈哈哈^_^ 2026-04-08 20/1000 2026-04-11 20:43 by stoner78
[考研] 电气工程专硕320求调剂 +4 小麻子111 2026-04-10 4/200 2026-04-11 10:19 by Delta2012
[考研] 282，求调剂 +12 jggshjkkm 2026-04-09 14/700 2026-04-11 09:39 by 猪会飞
[考研] 0854调剂 +4 长弓傲 2026-04-09 4/200 2026-04-11 09:18 by 猪会飞
[材料工程] 材料调剂推荐 +8 蛋糕x2 2026-04-07 8/400 2026-04-10 23:13 by Ftglcn90
[考研] 调剂 +12 卷卷卷心菜_ 2026-04-09 13/650 2026-04-10 22:36 by Ftglcn90
[考研] 0858求调剂 5+5 Gky09300550， 2026-04-10 8/400 2026-04-10 19:13 by chemisry
[考研] 调剂 +19 不逢春 2026-04-05 20/1000 2026-04-10 10:15 by may_新宇
[考研] 085404，285分求调剂 +12 薇薇考研 2026-04-07 14/700 2026-04-09 23:10 by parmtree
[考研] 生物学调剂，一志愿西南大学348，Top期刊一区二作、二区三作，三等奖学金三次 +4 candyyyi 2026-04-09 4/200 2026-04-09 18:39 by l_paradox
[考研] 求调剂材料科学与工程一志愿985初试365分 +5 材化李可 2026-04-08 5/250 2026-04-09 17:00 by Lilly_Li
[考研] 一志愿武理车辆 281 求调剂 +5 上岸研究生. 2026-04-07 5/250 2026-04-09 15:56 by only周
[考研] 368化学求调剂 +13 wwwwabcde 2026-04-07 14/700 2026-04-09 14:47 by heaven_jay
[考研] 336求调剂，一志愿中科大 +9 墨彧 yuyu 2026-04-06 9/450 2026-04-08 11:24 by 想读书的菌菌
[考研] 调剂求助（生物与医药） +6 @6952 2026-04-06 6/300 2026-04-07 23:52 by lys0704
[考研] 331求调剂 +5 张元一 2026-04-07 6/300 2026-04-07 22:13 by hemengdong
[考研] 一志愿河北工业大学材料工程，初试344求专硕调剂 +6 15933906766 2026-04-05 6/300 2026-04-06 13:21 by 无际的草原
信息提示
关闭
请填处理意见
关闭