| 查看: 7077 | 回复: 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,重新投任务呢？求高人指点迷津~我这个低精度的可行吗？大家在算的时候，怎么设置这个精度滴？ |
回复此楼» 猜你喜欢
投稿Elsevier的杂志（返修），总是在选择OA和subscription界面被踢皮球
已经有6人回复
-
中科院杭州医学所招收博士生一名（生物分析化学、药物递送）
已经有4人回复
-
求个博导看看
已经有18人回复
-
自荐读博
已经有6人回复
-
青基代表作，AAAI之类的A会的special track在国内认可度高吗？还是归为workshop之流？
已经有3人回复
-
上海工程技术大学【激光智能制造】课题组招收硕士
已经有6人回复
-
上海工程技术大学张培磊教授团队招收博士生
已经有4人回复
-
临港实验室与上科大联培博士招生1名
已经有9人回复
-
写了一篇“相变储能技术在冷库中应用”的论文，论文内容以实验为主，投什么期刊合适？
已经有6人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 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/165
-
经济学博士（金融方向）招生，211重点大学，2026年9月入学，申请-考核制。
+1/67
-
香港理工大学-应用生物与化学科技学系 招收2025年博士研究生
+2/62
-
所有人都在说要坚强
+4/44
-
征女友 @长安
+1/43
-
深圳理工大学梁国进课题组招聘研究助理教授、博后多名（电化学储能方向）
+1/41
-
西南交通大学前沿院碳中和与物质循环利用课题组招收博士生
+1/29
-
浙江师范大学申利国教授招聘博士后研究人员
+1/28
-
暨南大学理工学院 光子技术研究院段宣明团队申请制读博招生
+1/28
-
【AI、水文方向】香港科技大学（广州）研究助理招聘
+1/26
-
深圳大学李天任博士课题组研究生招生信息
+1/24
-
南京工业大学招收2026年全日制学术博士（供热、供燃气通风与空调）
+1/13
-
吉林大学材料物理本科生求问调剂信息
+1/13
-
青岛大学 丁欣 课题组 招收2026秋化学博士1名
+1/8
-
哈尔滨工业大学招收硕士研究生（欢迎环境、市政、生物、化学、农业等专业，长期有效）
+1/7
-
上海理工大学 生物医学工程专业 招收2026年全日制博士生一名
+1/5
-
武汉大学郭宇铮教授课题组招收博士后等研究人员【先进封装/芯片/人工智能等方向】
+1/4
-
福州大学梁宇航副教授招收2026年申核制博士研究生/硕士研究生（理论计算方向）
+1/3
-
湖南大学机械与运载工程学院赵岩副教授课题组招生2026级普通博士生1名
+1/2
-
有多余纯化系统，20-200mm高压制备分离系统，配套齐全可对外代工、委托加工、项目合作
+1/1
★ ★ ★
小木虫: 金币+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