| 查看: 7336 | 回复: 16 | |||
[交流]
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,重新投任务呢？求高人指点迷津~我这个低精度的可行吗？大家在算的时候，怎么设置这个精度滴？ |
回复此楼» 猜你喜欢
295求调剂。一志愿报考郑州大学化学工艺学硕，总分295分
已经有4人回复
-
289求调剂
已经有4人回复
-
江苏省农科院招调剂1名
已经有4人回复
-
290分材料工程085601求调剂 数二英一
已经有3人回复
-
0856调剂
已经有7人回复
-
化工专硕348，一志愿985求调剂
已经有9人回复
-
0856化工专硕求调剂
已经有14人回复
-
282求调剂
已经有5人回复
-
诚招农业博士
已经有3人回复
-
材料学硕318求调剂
已经有9人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 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人回复
» 抢金币啦！回帖就可以得到:
-
Chemosensors 期刊Special Issue 约稿
+2/90
-
青岛科技大学JIEQING团队招收材料科学、化学化工相关领域的硕/博研究生数名
+1/88
-
北京信息科技大学联合清华大学某国家重点实验室，联合培养硕生、博生，并长期招博士后
+1/88
-
大理大学博士研究生招生
+1/61
-
哈工大郑州研究院查正宝研究员招收博士生加盟
+1/43
-
西安工业大学庞教授课题组招收申请考核制 2026级学术博士生（电子/材料/物理/化学等）
+1/38
-
look 26届 || 计算机、电子信息类、电科、控制、通信考研T j信息醇：1o33 7747 o3
+1/20
-
长江大学农学院 2026 硕士招生
+1/16
-
look 26届 || 计算机、电子信息类、电科、控制、通信考研T j信息醇：1o33 7747 o3
+1/16
-
【加拿大卡尔顿大学】招博士生【心理学人工智能】
+1/15
-
吉林大学张越涛教授团队长期招收博士后、2026年入学博士研究生和科研助理
+1/14
-
电子科技大学材料学院SFT创新中心招收2026年博士生及科研助理
+1/13
-
中南林业科技大学生物质绿色转化与功能材料课题组现招收2026年秋季入学博士研究生
+1/9
-
西安电子科学大学杭州研究院刘丽香教授招收智能多模态传感器和微型储能器件方向博士
+1/8
-
环境微生物与污染控制研究所招收2026年硕士研究生
+1/8
-
深圳理工大学刘鑫课题组高薪招聘博士后
+1/7
-
look 26届计算机、电子信息类、电科、控制、通信考研T j信息醇：10+33+77+47+03
+1/7
-
26届计算机、电子信息类、电科、控制、通信考研T j信息pp骏：74+08+76+6+20/
+1/5
-
西交利物浦大学（苏州）/刘雯老师课题组/招博士研究生
+1/5
-
食品考研交流
+1/4
★ ★ ★
小木虫: 金币+0.5, 给个红包，谢谢回帖
贺仪: 金币+2, 谢谢指教！ 2012-05-31 08:34:33
小木虫: 金币+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的其他优化算法来算。 |
10楼2012-05-30 23:14:21
★ ★
小木虫: 金币+0.5, 给个红包，谢谢回帖
贺仪: 金币+1, 谢谢指教！ 2012-05-31 08:34:43
小木虫: 金币+0.5, 给个红包，谢谢回帖
贺仪: 金币+1, 谢谢指教！ 2012-05-31 08:34:43
|
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.” |
11楼2012-05-30 23:21:04
2楼2012-05-22 22:42:08
3楼2012-05-23 09:52:13
4楼2012-05-23 10:33:00
5楼2012-05-23 10:34:55
★ ★ ★
lsh1986914(金币+1): 谢谢参与
franch: 金币+2, 谢谢回帖交流 2012-05-23 19:43:56
lsh1986914(金币+1): 谢谢参与
franch: 金币+2, 谢谢回帖交流 2012-05-23 19:43:56
|
本帖内容被屏蔽 |
6楼2012-05-23 17:36:20
|
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
8楼2012-05-23 20:15:25
9楼2012-05-27 04:51:34
12楼2012-06-01 15:00:31
13楼2012-06-01 15:04:29
|
哦，谢谢了，我改变一下，参数设置。还有一些虫友说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
15楼2012-06-01 22:25:52
16楼2016-07-05 11:05:20
17楼2016-07-06 20:25:16
