24小时热门版块排行榜

返回列表

syding1531

银虫 (初入文坛)

应助: 1 (幼儿园)
金币: 30.7
散金: 300
红花: 4
帖子: 17
在线: 8小时
虫号: 571452
注册: 2008-06-10
专业: 电化学

[求助] K点设置报错：VERY BAD NEWS! internal error in subroutine IBZKPT已有1人参与

我是vasp的菜鸟，想就下面的问题请教各位大虾：

我想研究一下硫酸根在Au111表面上的吸附问题，单胞取为 sqrt(3) x sqrt(7).
短长轴之比约为0.6547，于是我设置K点为 8,5,1  （5/8=0.625）
KPOINT文件为：
--------------------------
Automomatic mesh
0
Monkhorst-Pack
8  5  1
0  0  0
--------------------------

于是就报出下面的错，并且任务彻底停掉了。
VERY BAD NEWS! internal error in subroutine IBZKPT:
Reciprocal lattice and k-lattice belong to different class of lattices. Often results are still useful...    96
VERY BAD NEWS! internal error in subroutine SGRCON:
Found some non-integer element in rotation matrix    5

（2）我把对称性关掉，
还是出现这样的报错

（3）我把K点改成7,5,1，(5/7约为0.7143)
--------------------------
Automomatic mesh
0
Monkhorst-Pack
7  5  1
0  0  0
--------------------------
结果还是会报
VERY BAD NEWS! internal error in subroutine IBZKPT:
Reciprocal lattice and k-lattice belong to different class of lattices. Often results are still useful...    96
但是没有
VERY BAD NEWS! internal error in subroutine SGRCON:
Found some non-integer element in rotation matrix    5
而且vasp程序一直还在run，没有跳出来。

我的INCAR为
-----------------------------
System=SO4onAulll
ISPIN=1
ISTART=0
ICHARG=2

IBRION=2
ISIF=2
POTIM=0.5
ISYM=0

LREAL=Auto
PREC=A

GGA=PE
ENCUT=400

EDIFF=1.0E-05
EDIFFG=-0.01
NELM=120
NSW=240
#NPAR=8

#LVTOT=.TURE.
#LELF=.TURE.

ISMEAR=1
SIGMA=0.2
IALGO=48
---------------------

POSCAR为：
--------------------------------------------------------------------
SO4onAulll
1.00
  5.11208509 0.00000000 0.00000000
  2.55604255 7.37865926 0.00000000
  0.00000000 0.00000000 28.91832000
  25 1 4
Selective Dynamic
Cartesian
4.68607809 5.16506132 0.00000000 F    F    F
2.13003555 3.68932947 0.00000000 F    F    F
7.24212064 6.64079318 0.00000000 F    F    F
4.68607800 2.21359778 0.00000000 F    F    F
2.13003546 0.73786592 0.00000000 F    F    F
3.83406382 3.68932963 2.40986002 F    F    F
1.27802127 2.21359778 2.40986002 F    F    F
6.39010637 5.16506148 2.40986002 F    F    F
3.83406382 0.73786593 2.40986002 F    F    F
3.83406382 6.64079333 2.40986002 F    F    F
2.98204955 2.21359794 4.81972004 T    T    T
0.42600700 0.73786608 4.81972004 T    T    T
5.53809209 3.68932979 4.81972004 T    T    T
5.53809218 6.64079334 4.81972004 T    T    T
2.98204964 5.16506148 4.81972004 T    T    T
4.68607809 5.16506132 7.22958007 T    T    T
2.13003555 3.68932947 7.22958007 T    T    T
7.24212064 6.64079318 7.22958007 T    T    T
4.68607800 2.21359778 7.22958007 T    T    T
2.13003546 0.73786592 7.22958007 T    T    T
3.83406382 3.68932963 9.63944009 T    T    T
1.27802127 2.21359778 9.63944009 T    T    T
6.39010637 5.16506148 9.63944009 T    T    T
3.83406382 0.73786593 9.63944009 T    T    T
3.83406382 6.64079333 9.63944009 T    T    T
2.98204955 2.21359794 12.44497618 T    T    T
1.27802127 2.21359778 11.82944009 T    T    T
3.83406382 0.73786593 11.82944009 T    T    T
3.83406382 3.68932963 11.82944009 T    T    T
2.98204955 2.21359794 13.97178357 T    T    T
---------------------------------------------------------------------------

vasp论坛上也有人提这样的问题，http://cms.mpi.univie.ac.at/vasp-forum/forum_viewtopic.php?3.489
但是也没有说的很具体，看不懂。

道理上 8 5 1是最合适的，如何实现用8 5 1的K点？我知道文献中用 SeqQuest软件是可以这么设置的。
请大家帮帮我看看到底是什么原因好吗？谢谢
附件为 POTCAR[ Last edited by syding1531 on 2011-9-23 at 22:07 ]

回复此楼

» 猜你喜欢

寻求一种能扛住强氧化性腐蚀性的容器密封件已经有7人回复
到新单位后，换了新的研究方向，没有团队，持续积累2区以上论文，能申请到面上吗已经有8人回复
申请2026年博士已经有6人回复
请问哪里可以有青B申请的本子可以借鉴一下。已经有5人回复
天津工业大学郑柳春团队欢迎化学化工、高分子化学或有机合成方向的博士生和硕士生加入已经有5人回复
2025冷门绝学什么时候出结果已经有7人回复
请问有评职称，把科研教学业绩算分排序的高校吗已经有6人回复
Bioresource Technology期刊，第一次返修的时候被退回好几次了已经有7人回复
请问下大家为什么这个铃木偶联几乎不反应呢已经有5人回复
康复大学泰山学者周祺惠团队招收博士研究生已经有6人回复

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

origin这样的坐标刻度怎么设置已经有3人回复
求助关于vasp磁矩设置的说明书的详细理解已经有12人回复
Autodock时 gridbox的参数怎么设置?? 已经有18人回复
VASP报错：VERY BAD NEWS! 已经有5人回复
interface的设置已经有6人回复
【求助】gaussian03 error2070 已经有10人回复
【求助】VERY BAD NEWS! internal error in subroutine IBZKPT 已经有13人回复
【求助】Error termination request processed by link 9999是什么错误？程序在下面已经有13人回复
【求助成功】大家帮忙改改错 internal error in subroutine INVGRP: 已经有6人回复
【求助】求VASP计算ZnO能带和态密度的详细输入设置已经有15人回复
【求助】FLUENT中的内热源怎么设置？在哪里设置？已经有16人回复

1楼2011-09-23 21:51:33

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

utdwang

银虫 (初入文坛)

应助: 0 (幼儿园)
金币: 455.7
红花: 1
帖子: 49
在线: 46.6小时
虫号: 837149
注册: 2009-08-31
性别: GG
专业: 催化化学

【答案】应助回帖

★
zzy870720z(金币+1): 谢谢建议 2011-09-24 17:13:46
syding1531(金币+2): 这个我也知道，但是为什么要用Gamma呢？ 2011-09-26 01:52:23

Not sure. you can try to switch your k-point scheme 'Monkhorst-Pack' to 'Gamma'.

赞一下(1人)

回复此楼

Thekeytosuccessinnearlyeverythingisdailyaction

2楼2011-09-24 00:13:25

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

jinhua1983

金虫 (小有名气)

应助: 2 (幼儿园)
金币: 737.4
散金: 1010
红花: 3
帖子: 161
在线: 60.6小时
虫号: 713045
注册: 2009-03-02
性别: MM
专业: 理论和计算化学

【答案】应助回帖

★
uuv2010(金币+1): 多谢提示 2011-09-27 09:43:58

我之前貌似也出过类似的错误，建议你可以改一下计算方法试试，用38算，不用48
。还有我算表面的时候一直ISIF用的是0，很少用2。一般是表面用0，块才用3吧，还有就是ISMEAR这个关键词对于表面是用0的，用1好像是用在纯金属的体系吧，具体的不是很记得了~
也期待听到你的好声音，弄清楚如何解决后也麻烦告知一下，大家相互学习！

赞一下(2人)

回复此楼

3楼2011-09-27 08:32:11

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