24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

lsh1986914

新虫 (小有名气)

应助: 1 (幼儿园)
金币: 90.6
帖子: 68
在线: 18.1小时
虫号: 1145561

[交流] 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楼 2012-05-22 16:08:54

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

huangguiyang

新虫 (小有名气)

1ST强帖: 1
应助: 4 (幼儿园)
金币: 176.1
帖子: 84
在线: 12.9小时
虫号: 483011

★ ★
小木虫: 金币+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.”