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

24小时热门版块排行榜    

Znn3bq.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » VASP计算中KPOINTS文件的作用？
51/1返回列表
查看: 4012  |  回复: 8
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

仲_______

金虫 (小有名气)

[求助] VASP计算中KPOINTS文件的作用？ 已有1人参与

大家好，本人初学vasp，对KPOINTS文件的作用不太理解。有人能比较形象地解释一下吗？

此外，今天做了两次bulk Fe的静态自洽计算，四个输入文件中只有KPOINTS文件改变了，然后得出的OSZICAR文件中的能量就变化了。
第一次情况：
KPOINTS为：
Gamma-point only
1          # one k-point
rec        #  in units of the reciprocal lattice vector
0 0 0 1    # 3 coordinates and weight
得到的OSZICAR文件为：
       N       E                     dE             d eps       ncg     rms          rms(c)
DAV:   1     0.379178474714E+02    0.37918E+02   -0.57325E+03    24   0.105E+03
DAV:   2    -0.222250175980E+02   -0.60143E+02   -0.60137E+02    24   0.207E+02
DAV:   3    -0.229294693776E+02   -0.70445E+00   -0.70445E+00    24   0.254E+01
DAV:   4    -0.229411027224E+02   -0.11633E-01   -0.11633E-01    24   0.362E+00
DAV:   5    -0.229412975069E+02   -0.19478E-03   -0.19478E-03    24   0.365E-01    0.122E+01
DAV:   6    -0.220321364331E+02    0.90916E+00   -0.37594E+00    24   0.310E+01    0.440E+00
DAV:   7    -0.220457337624E+02   -0.13597E-01   -0.23265E-02    24   0.204E+00    0.202E+00
DAV:   8    -0.220600462699E+02   -0.14313E-01   -0.66212E-03    24   0.107E+00    0.145E-01
DAV:   9    -0.220600639976E+02   -0.17728E-04   -0.39336E-06    12   0.301E-02
   1 F= -.22060064E+02 E0= -.22056723E+02  d E =-.668197E-02


第二次的情况：
KPOINTS为：
K-Points
0
Monkhorst Pack
11 11 11
0  0  0
得到的OSZICAR文件为：
       N       E                     dE             d eps       ncg     rms          rms(c)
DAV:   1     0.308870426411E+02    0.30887E+02   -0.54665E+03  1368   0.975E+02
DAV:   2    -0.179773985631E+02   -0.48864E+02   -0.48455E+02  1452   0.191E+02
DAV:   3    -0.184450276628E+02   -0.46763E+00   -0.46465E+00  1416   0.199E+01
DAV:   4    -0.184499709257E+02   -0.49433E-02   -0.49387E-02  1620   0.216E+00
DAV:   5    -0.184500084488E+02   -0.37523E-04   -0.37522E-04   948   0.140E-01    0.150E+01
DAV:   6    -0.157246299600E+02    0.27254E+01   -0.34328E+01  1344   0.829E+01    0.415E+00
DAV:   7    -0.154153110592E+02    0.30932E+00   -0.26867E+00  1344   0.219E+01    0.111E+00
DAV:   8    -0.153894943981E+02    0.25817E-01   -0.34199E-02  1344   0.195E+00    0.420E-01
DAV:   9    -0.153823563493E+02    0.71380E-02   -0.46595E-03  1020   0.511E-01    0.836E-02
DAV:  10    -0.153823954050E+02   -0.39056E-04   -0.15197E-04   672   0.138E-01
   1 F= -.15382395E+02 E0= -.15302896E+02  d E =-.158999E+00

能看出最后的能量E0相差还是挺大的。请问应该这是怎么回事？谢谢！
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2014-07-18 11:45:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

仲_______

金虫 (小有名气)
送红花一朵
引用回帖:
4楼: Originally posted by valenhou001 at 2014-07-18 16:10:16
这是对这样的分子计算，采用超原胞，而且取的很大，那对应的布里渊区就很小，取Gamma点（用一个k点）取样就足够了。

你可能需要看看固体物理理解一下相关的基本知识（晶体、基矢、布里渊区、波矢、原胞等等）。...

好的，非常感谢。
我本科学的是通信，然后现在研一导师让我用vasp计算，所以一些基础不太懂。目前正在学习固体物理学和固体能带理论，希望慢慢讲基础补上去。
非常感谢您的回复！
赞一下 回复此楼
5楼2014-07-18 16:23:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 9 个回答

valenhou001

至尊木虫 (职业作家)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
仲_______: 金币+10, ★★★★★最佳答案 2014-07-18 15:41:41
你计算的体材料，而且采用的初基原胞，当然不能用单k点（Gamma点）来进行计算。两次计算得到的总能当然差别很大。

强烈建议你在做计算之前读读计算方法相关的文献，比如k-point sampling之类的。另外理解布里渊区的概念。有些教程里给出了k-point sampling的简单图示。
https://www.vasp.at/vasp-workshop/slides/k-points.pdf
http://www.jncasr.ac.in/ccms/fps ... day1/waghmare_1.pdf
赞一下 回复此楼
2楼2014-07-18 13:00:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

仲_______

金虫 (小有名气)
引用回帖:
2楼: Originally posted by valenhou001 at 2014-07-18 13:00:54
你计算的体材料，而且采用的初基原胞，当然不能用单k点（Gamma点）来进行计算。两次计算得到的总能当然差别很大。

强烈建议你在做计算之前读读计算方法相关的文献，比如k-point sampling之类的。另外理解布里渊区 ...

您好，非常感谢您的回答和文献！
我在学习vasp官网给的教程“hand on session I”,其中提到“for molecules and atoms only a single k-point is required”. 然后在例子中，对于O原子、O2分子、CO分子、H2O分子都采用的是我上文中提到的单K点的KPOINTS。请问这个是为什么呢？谢谢！
赞一下 回复此楼
3楼2014-07-18 15:44:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

valenhou001

至尊木虫 (职业作家)
引用回帖:
3楼: Originally posted by 仲_______ at 2014-07-18 16:44:57
您好，非常感谢您的回答和文献！
我在学习vasp官网给的教程“hand on session I”,其中提到“for molecules and atoms only a single k-point is required”. 然后在例子中，对于O原子、O2分子、CO分子、H2O分子都 ...

这是对这样的分子计算，采用超原胞，而且取的很大，那对应的布里渊区就很小，取Gamma点（用一个k点）取样就足够了。

你可能需要看看固体物理理解一下相关的基本知识（晶体、基矢、布里渊区、波矢、原胞等等）。
赞一下 回复此楼

» 本帖已获得的红花（最新10朵）

仲_______
4楼2014-07-18 16:10:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 9 个回答
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 材料化工总分334求调剂 +14 Riot2025 2026-04-08 14/700 2026-04-09 17:36 by yongzhesheng
[考研] 材料调剂 +10 18815505510 2026-04-09 11/550 2026-04-09 17:07 by 544594351
[考研] 一志愿郑州大学 22408 305分求调剂 +3 安小满zzz 2026-04-08 3/150 2026-04-09 16:16 by wp06
[考研] 085501机械英二77总分294求调剂，接受跨专业学习 +6 守法公民亓纪 2026-04-08 6/300 2026-04-09 15:55 by wp06
[考研] 26考研调剂0710 0860 +10 补补不补 2026-04-03 15/750 2026-04-09 11:13 by ditto77778
[基金申请] 有爆料，一个青年教师卖房得400万，然后换了一个四青帽子 +6 babu2015 2026-04-08 6/300 2026-04-09 08:14 by kudofaye
[考研] 材料考研求调剂总分280 +30 mkjlz1 2026-04-06 35/1750 2026-04-08 21:25 by cyh—315
[考研] 283分求调剂 +14 试试看呗 2026-04-04 14/700 2026-04-08 07:03 by lijunpoly
[考研] 计算机408｜在校多次国家级竞赛获奖｜申请调剂 +4 东山大白鹅 2026-04-05 4/200 2026-04-08 00:18 by chongya
[考研] 277、学硕，求调剂 数一104， +11 瓶子PZ 2026-04-07 12/600 2026-04-07 23:30 by 一只好果子?
[考研] 307求调剂 +3 Youth@@ 2026-04-07 3/150 2026-04-07 22:00 by hemengdong
[考研] 081700学硕，323分，一志愿中国海洋大学求调剂学校 +19 披星河 2026-04-04 19/950 2026-04-07 15:00 by 上岸快快
[考研] 工科 22408 267求推荐 +4 wanwan00 2026-04-05 5/250 2026-04-06 22:47 by chenzhimin
[考研] 工科370求调剂 +3 溏心煎鸡蛋 2026-04-05 3/150 2026-04-06 10:55 by 这是一个无聊的�
[考研] 302分 085601求调剂推荐 +11 zyx上岸！ 2026-04-05 11/550 2026-04-05 22:13 by dongzh2009
[考研] 271分求调剂学校 +12 zph158488！ 2026-04-02 13/650 2026-04-05 10:13 by lqwchd
[考研] 考研调剂 +11 小sun要好运 2026-04-04 11/550 2026-04-05 08:02 by qlm5820
[考研] 306求调剂 +3 hyb上名工 2026-04-02 3/150 2026-04-04 18:12 by 热情沙漠
[考研] 11408 一志愿西电，277分求调剂 +4 zhouzhen654 2026-04-03 4/200 2026-04-04 18:10 by 猪会飞
[考研] 考研调剂 +3 Draa 2026-04-03 3/150 2026-04-03 17:37 by hgwz7468
信息提示
关闭
请填处理意见
关闭