24小时热门版块排行榜    

查看: 3486  |  回复: 20

losenq

新虫 (小有名气)

[求助] ENCUT收敛测试时,能量差为什么总不能达到1meV呢?

ENCUT收敛测试结果:

ENCUT    Energy
400       -308.42528
450      -308423939
500      -308.452559
550      -308.51164
600      -308.59528
650     -308.676606
700     -308.744955
750     -308.784962
800    -308.811974
850   -308.824027

计算中用的K点经过测试了。
INCAR文件也没啥特别的,就是静态计算那四个。。。。
不知道怎么回事,能量随着ENCUT的增加有一直减小的趋势。。。差值怎么也不能到1meV啊。
各位大神,我该怎么办呢?增加K点也试过了,结果差不多!!
肿么办,肿么办,好捉急啊。。。。然后POTCAR默认的ENCUT=400。。。。
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

identation

金虫 (著名写手)

★ ★
franch: 金币+2, 谢谢回帖交流 2013-11-11 22:56:17
???????:
11?: Originally posted by jpchou at 2013-11-04 16:52:13
??
????????@??
?????????
??????????
??? identation ????
??F??????????
??]????@????y...

vasp?????????????????????????+3????????Sm3+;?????????????gga??Sm_3?????PAW????f???????о??
12楼2013-11-04 19:40:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

liusj0228

银虫 (小有名气)

【答案】应助回帖

★ ★ ★
losenq: 金币+2, 有帮助 2013-11-04 22:16:07
identation: 金币+1, 谢谢交流 2013-11-11 16:34:30
应该是总能量除以原子数,也就是每个原子的能量差收敛到1meV
请各位高手不吝赐教
15楼2013-11-04 21:38:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

jpchou

专家顾问 (著名写手)

【答案】应助回帖

★ ★ ★
感谢参与,应助指数 +1
sunyang1988: 金币+1, 谢谢交流 2013-11-01 22:56:01
losenq: 金币+2, 有帮助, 我用的EIDFF=1E-6,EDIFFG=-0.01收敛条件已经相当高了 2013-11-02 20:21:15
我猜 ...... EDIFF 加深
預設是 1E-04
算一輩子能量差都不會小於1meV
用力的把它加到 1E-6 吧
但切記
不要低於 1E-8 (or -7 ?)
那是無意義的行為
2楼2013-11-01 04:36:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jpchou

专家顾问 (著名写手)

★ ★
identation: 金币+2, 谢谢指教 2013-11-04 10:28:38
???????:
2?: Originally posted by jpchou at 2013-11-01 04:36:38
??? ...... EDIFF ????
?A?O?? 1E-04
????????????????С?1meV
???????????? 1E-6 ??
?????
?????? 1E-8 (or -7 ?)
????o???x???О?

??
???м????????????data
??????????????c?F??
??????????????????????
??????????
??????y????N?????y??
3楼2013-11-03 00:02:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

losenq

新虫 (小有名气)

???????:
3?: Originally posted by jpchou at 2013-11-03 00:02:21
??
???м????????????data
??????????????c?F??
??????????????????????
??????????
??????y????N?????y??...

????????????????

??в???ENCUT?????????????????????????????????1meV.

????KPOINT???????????????????meV.???????????????
??????????ENCUT????????????????????????????????ENCUT=600????????????ο?????????MS?е?CASETEP??????

CASETEP?????????VASP?????
????????????????????????塣
4楼2013-11-03 21:42:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

losenq

新虫 (小有名气)

???????:
3?: Originally posted by jpchou at 2013-11-03 00:02:21
??
???м????????????data
??????????????c?F??
??????????????????????
??????????
??????y????N?????y??...

???????????????????????????????????????????л????????????????Щ?????????С???????
5楼2013-11-03 22:09:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jpchou

专家顾问 (著名写手)

★ ★ ★
franch: 金币+3, 谢谢回帖交流 2013-11-04 17:26:58
???????:
5?: Originally posted by losenq at 2013-11-03 22:09:10
???????????????????????????????????????????л????????????????Щ?????????С???????...

???
??????}????????
????f?????????????y????N???
INCAR ????N???
K?c?x?
?c VASP ?汾

??????
ENCUT ????y??c????????????????P?S
?????F?????? ENCUT ????y?
?Л]?п??]?????????]???κ??P?S??
6楼2013-11-04 01:44:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

losenq

新虫 (小有名气)

★ ★
franch: 金币+2, 鼓励交流。 2013-11-04 17:27:09
???????:
6?: Originally posted by jpchou at 2013-11-04 01:44:10
???
??????}????????
????f?????????????y????N???
INCAR ????N???
K?c?x?
?c VASP ?汾

??????
ENCUT ????y??c????????????????P?S
?????F?????? ENCUT ????y?
?Л]?п??]?????? ...

INCAR???????????????????????????????
ISTATE=0
ICHRANG=2
ISMEAR=-5
EDIFF=1E-6
PREC=ACCURATE
???????????????????????????ɡ?
K??????????4*4*4

?????K????????????ENCUT=550?????????????????ENMAX=400??????400*1.3=520?????????550????????K???4?????????4*4*4??5*5*5??6*6*6??7*7*7??8*8*8?????????????С??????????meV??????????4*4*4???

????п?????????????????????PAW_PBE

????????????????ο?????????CASETEP???????600?????????1.3*ENMAX?????????VASP???????600????????????????
7楼2013-11-04 10:16:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

losenq

新虫 (小有名气)

???????:
6?: Originally posted by jpchou at 2013-11-04 01:44:10
???
??????}????????
????f?????????????y????N???
INCAR ????N???
K?c?x?
?c VASP ?汾

??????
ENCUT ????y??c????????????????P?S
?????F?????? ENCUT ????y?
?Л]?п??]?????? ...

VASP??5.2
8楼2013-11-04 10:17:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

identation

金虫 (著名写手)

★ ★
franch: 金币+2, 谢谢回帖交流 2013-11-04 17:27:19
引用回帖:
5楼: Originally posted by losenq at 2013-11-03 22:09:10
我在捉摸是不是我没考虑到自旋极化之类的问题。。。。有机物的话是需要考虑这些么??迷惑中。。。。...

这个和无机有机没关系 如果含有磁性离子,就需要考虑了
9楼2013-11-04 10:29:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

identation

金虫 (著名写手)

★ ★
franch: 金币+2, 谢谢回帖交流 2013-11-04 17:27:27
引用回帖:
7楼: Originally posted by losenq at 2013-11-04 10:16:57
INCAR确实没啥新鲜的。。。就那么几个参数啊。
ISTATE=0
ICHRANG=2
ISMEAR=-5
EDIFF=1E-6
PREC=ACCURATE
就这么几个,静态计算也就这样了吧。
K点自动生成,4*4*4

当时做K点测试的时候,选的ENCUT=550, ...

这个是有可能的,我怀疑你的体系里有稀土元素,不知是不是这样?
稀土元素体系收敛不容易,因为f电子的存在
10楼2013-11-04 10:30:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 losenq 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 科研孤儿太难了 +14 我4大白菜 2026-08-20 15/750 2026-08-20 22:33 by 袁向阳007
[基金申请] 让我中一个面上吧! +8 大萍1987 2026-08-20 10/500 2026-08-20 21:48 by 大萍1987
[基金申请] 时间戳又变了 +12 wuchongjun 2026-08-20 17/850 2026-08-20 20:05 by spzhao
[基金申请] 今天系统多次维护,明天很可能放榜! +10 zju2000 2026-08-16 11/550 2026-08-20 20:02 by cl479861084
[基金申请] 时间戳变了,能看出什么问题? +18 基诺咪客 2026-08-17 23/1150 2026-08-20 17:19 by Godzela
[基金申请] 放榜前的不淡定 40+3 snowwithsea 2026-08-19 12/600 2026-08-20 15:56 by lfy8008
[基金申请] 重要消息,中午系统在维护 +11 yuleib84 2026-08-18 12/600 2026-08-20 11:09 by xskun
[基金申请] 今天基金会出结果吗?20260819 +14 kkkl_v 2026-08-19 15/750 2026-08-20 08:21 by sunzitan
[基金申请] filecode,4个jtjc了 +11 ziyangfang 2026-08-19 13/650 2026-08-19 21:28 by aasahr
[教师之家] 为什么余额宝的年化利率越来越低?主要原因有哪些? +5 瞬息宇宙 2026-08-15 5/250 2026-08-19 20:23 by super2002521
[基金申请] 今天放榜吗? +14 布布和一二 2026-08-19 15/750 2026-08-19 18:07 by gltch
[基金申请] filecode=后面第一个是大写字母 +10 wangze12014 2026-08-14 12/600 2026-08-19 16:56 by 苦难博士
[基金申请] 朋友圈看到的 +6 wangzilk 2026-08-18 8/400 2026-08-19 10:55 by Haru815
[基金申请] 什么时候开奖? +6 CrisMessi 2026-08-18 6/300 2026-08-19 08:06 by Equinoxhua
[基金申请] 明天放榜? +5 Shxjjxjkx 2026-08-18 5/250 2026-08-18 18:14 by -大大大大大-
[论文投稿] 投稿咨询 +4 wwm09 2026-08-17 6/300 2026-08-18 15:36 by wwm09
[基金申请] 今天维护系统维护 祝所有人 高中 +8 gjjjzhong 2026-08-18 9/450 2026-08-18 13:01 by 家与远方
[基金申请] 93BebMhtakh前后11位开头都是大写 +4 且听虎啸 2026-08-17 5/250 2026-08-18 00:49 by 蔡棒棒菂
[基金申请] 感觉是下周放榜了 +6 angus9576 2026-08-17 11/550 2026-08-17 23:57 by angus9576
[精细化工] 招聘 金属平磨液,抛光液研发工程师 +3 小天0311 2026-08-14 3/150 2026-08-16 07:31 by H9PLUS
信息提示
请填处理意见