24小时热门版块排行榜    

查看: 5219  |  回复: 34
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

小熊尼尼

木虫 (正式写手)

[求助] VASP静态计算中体系不收敛,恳请各位前辈指教 已有1人参与

体系含32个原子,从木虫上搜索过各种处理方法,经修改测试,结果还是不收敛。
INCAR文件如下:
SYSTEM = ********
ISTART = 0
ICHARG = 2
PREC   = Accurate
LREAL  = .F.
ENCUT  = 220
ISPIN  = 2
ISYM   = 0
NELM   = 200
EDIFF  = 1E-04
ISMEAR = 1
SIGMA  = 0.5
IBRION = -1
NSW    = 300
POTIM  = 0.5
EDIFFG = -0.05
ALGO   = Fast
NELMDL = -5
在上面参数中,尝试修改过NELM,ISMEAR,SIGMA,ALGO,NELMDL等参数设置,均未能达到收敛。ENCUT设置为220,是因为现在在做前期的测试。K点为3x3x3。POSCAR由驰豫收敛后的CONTCAR得到,但驰豫是使用的ENCUT为370(小虫想这里应该没有问题)。
在以上述INCAR参数进行计算时,使用grep 'energy without entropy' OUTCAR命令,抓到的结果如下图。恳请前辈指教。如需要其他文件,小虫继续上传。

VASP静态计算中体系不收敛,恳请各位前辈指教
grep energy 截图.png
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

小熊尼尼

木虫 (正式写手)

引用回帖:
15楼: Originally posted by magicmonk at 2016-09-06 19:38:45
我想问一下,你认为静态计算没有收敛的依据是什么……还有以后计算尽量让截断能保持一致

看OUTCAR得情况,应该就不是能量没有收敛了;剩下的就是力收敛问题。今天做过测试,分别设置如下参数:EDIFF=1E-03;EDIFFG=-0.1。结果还是未能收敛。然后我请师兄给我发了一个他做的可以收敛的算例,先正在测试,其中参数主要变化有:NSW采用默认(0),IBRION采用默认(-1),ISIF采用默认(2),NELM=100;力收敛采用的是EDIFFG=-0.05,EDIFF=1E-04。
很不好的是,在编辑给您的这个回帖的时候,结果出来了 ,虽然OUTCAR已给出了‘Total cpu time used (sec): ’,但是并没有提示‘reach the accuracy’。
之所以在后面的计算中采用了ENCUT=650的设置,是因为昨天晚上经教导,说用大些的ENCUT来算,说ENCUT太小的话可能一辈子都不会收敛。(虽然觉得个人ENCUT不会影响是否收敛,但是为了解决问题,还是采用了尽可能大的ENCUT)
前辈,下一步,我想再刚才的设置的基础上再测试一下EDIFF和EDIFFG,如果还不能收敛,我真的找不到怎样再继续测试以达到收敛了。请前辈不吝指教。
16楼2016-09-06 21:00:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 35 个回答

从心而觅

新虫 (小有名气)

啥叫静态计算不收敛啊……
一切有为法,如梦幻泡影!如露亦如电,应作如是观!
2楼2016-09-05 07:50:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

从心而觅

新虫 (小有名气)


蛋蛋小童鞋: 金币+1, 谢谢参与 2016-09-05 11:49:24
一切有为法,如梦幻泡影!如露亦如电,应作如是观!
3楼2016-09-05 08:22:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

小熊尼尼

木虫 (正式写手)

引用回帖:
2楼: Originally posted by 从心而觅 at 2016-09-05 07:50:12
啥叫静态计算不收敛啊……

单点能计算呀

发自小木虫Android客户端
4楼2016-09-05 10:26:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 投票:  有多少人是今天查系统知道结果的? +17 爱看书的可乐 2026-08-26 19/950 2026-09-02 00:19 by xiangy672
[基金申请] 面上函评意见出来了,像什么等级? 20+4 Tsingking1 2026-08-27 17/850 2026-09-01 19:51 by 超级无敌华子
[基金申请] 为什么国自然不能直接公布 +5 bjdxyxy 2026-08-26 5/250 2026-09-01 17:54 by 小伟大博士
[基金申请] 学科评审组评审是指会评吗? +4 瞬息宇宙 2026-08-31 4/200 2026-09-01 14:58 by jiaoxg
[论文投稿] 小白求助 投论文要求的highlights应该如何写 5+3 l1963982152 2026-08-29 4/200 2026-09-01 09:04 by 北京莱茵编辑
[基金申请] 麻烦专家们看看评委们的意见(F口面上) +7 gdd2018 2026-08-28 12/600 2026-09-01 08:32 by 尼古拉斯小虫
[基金申请] 怎么看青基中了没有啊 +6 叶九微 2026-08-26 6/300 2026-08-31 23:54 by yudaoqian88
[基金申请] 国社科又开始会评了,不知道这次命运如何 +7 雨打竹帘 2026-08-30 11/550 2026-08-31 23:16 by hittle2008
[基金申请] 面上意见出来了 +12 黄鸟于飞Chao 2026-08-29 23/1150 2026-08-31 18:57 by 黄鸟于飞Chao
[基金申请] 中青基了要发朋友圈吗? +7 349506619 2026-08-28 7/350 2026-08-31 13:39 by 冼亮淀粉酶
[基金申请] 基金不中,共勉 +12 eulota 2026-08-26 12/600 2026-08-31 08:42 by ZJTJZ
[基金申请] 有没有仍没收到信息的 +7 德尚中行 2026-08-27 8/400 2026-08-30 20:52 by purplejack
[基金申请] 我就是申请一个面上项目而已,这评审意见是按照杰青的条件评的吧? +6 gouxfjh 2026-08-28 11/550 2026-08-30 07:57 by gouxfjh
[基金申请] 国自然面上复盘~欢迎讨论 (金币+15) +15 晴天加油 2026-08-26 16/800 2026-08-29 18:28 by symmetry
[基金申请] 2026年叶企孙基金 +4 bud_bud 2026-08-27 7/350 2026-08-29 07:23 by foolishmani
[基金申请] 怎么查啊 +6 huang1991js 2026-08-26 6/300 2026-08-28 08:42 by winsaint
[基金申请] 看板上这么多中的,有点像50人群里49个人都是骗子的那种感觉…… +5 a089 2026-08-26 6/300 2026-08-27 14:05 by jonewore
[基金申请] 为什么 国际(地区)合作与交流项目 没有放榜? 10+3 majunge000 2026-08-26 11/550 2026-08-27 08:42 by 北京莱茵编辑
[基金申请] 出来了 +9 trojank 2026-08-26 9/450 2026-08-26 14:25 by 宝贝虫子
[基金申请] 国际合作可查了,中了面上 (EPI+1)(金币+50) +18 Ldrop2023 2026-08-26 18/900 2026-08-26 11:15 by cmrandy
信息提示
请填处理意见