应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » ENCUT收敛测试时，能量差为什么总不能达到1meV呢？
211/3返回列表123下一页
查看: 3024  |  回复: 20
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

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。。。。
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2013-10-31 22:49:25
已阅   回复此楼   关注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???????о??
赞一下(1人) 回复此楼
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
赞一下(1人) 回复此楼
请各位高手不吝赐教
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 的主题更新
211/3返回列表123下一页
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭