24小时热门版块排行榜    

CyRhmU.jpeg
查看: 3376  |  回复: 5

losenq

新虫 (小有名气)

[求助] VASP直接计算弹性常数时报错:Very bad news! internal error in subroutine IBZKPT:

用VASP直接计算弹性常数时报错:

Very bad news! internal error in subroutine IBZKPT: Reciprocal lattice and k-lattice belong to different class of lattice. Often results are still useful....

先说一下我在网上搜索的一些解决办法:
1。改变K点设置。对于一些低对称性的体系,偶数K点有时会导致这种情况的发生。将偶数K点改为奇数。
2.设置ISYM=0.
3.设置SYMPREC=1E-8.

以上方法我都试过了,但仍旧未解决问题。期待高手!多谢了!!

体系是NiTi的B2相,体心立方。

以下是我用的INCAR文件:
ENCUT=700
ISTART=0
ICHARG=2
ISMEAR=1
IBRION=6
NSW=1
ISIF=3
POTIM=0.01
NFREE=4
ISIM=0
EDIFF=1E-6
EDIFFG=-0.01
PREC=ACCURATE

另外我的K点设置是
15 15 15 自动产生
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xwangan

金虫 (初入文坛)

【答案】应助回帖


sunyang1988: 金币+1, 谢谢提示 2013-01-03 12:55:35
改变K点个数试试
Tomorrowisanotherday
2楼2013-01-01 21:44:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

匿名

用户注销 (小有名气)

饼干队长

★ ★
franch: 金币+2, 谢谢回帖交流,, 2013-01-07 15:35:17
本帖仅楼主可见
3楼2013-01-07 00:18:28
已阅   申请1ST强帖   回复此楼   编辑   查看我的主页

fzuvivi

金虫 (小有名气)

【答案】应助回帖

★ ★
franch: 金币+2, 谢谢回帖交流,, 2013-01-07 15:35:28
这个是因为你的K点设置不合理造成的。比如你的胞是六方的,你给的却是立方的Kgrid。可以尝试改变KPOTINS里的K点情况来解决
4楼2013-01-07 08:42:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fanarsenal

木虫 (正式写手)

【答案】应助回帖

★ ★
franch: 金币+2, 谢谢回帖交流,, 2013-01-07 15:35:39
之前我计算B2结构的时候也遇到这个问题,我之前是因为精度设置太高了,计算能力达不到。个人认为,体系精度要求太高,在对称性略有破坏的情况下vasp就会认为正倒空间对称性不一致。
建议lz像我一样增加cpu数目重新算一下,或者降低一下精度。这个结果如报错信息,肯定是有意义的。
5楼2013-01-07 09:42:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

losenq

新虫 (小有名气)

引用回帖:
3楼: Originally posted by 旭旭宝宝 at 2013-01-07 00:18:28
1.你的材料有NI  所以要考虑自旋   ,2 KPOINTS 取的有点大  我56个原子才取 4*4*4  , K点 最好取偶数!!3.你的B2结构是否是 立方?

K点是收敛测试结果,只有取15以上,能量差才会在1meV内。取奇数是因为手册上有说:K<8时最好取偶数,K>8时最好取奇数。
6楼2013-11-06 09:00:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 losenq 的主题更新
信息提示
请填处理意见