24小时热门版块排行榜    

Znn3bq.jpeg
北京石油化工学院2026年研究生招生接收调剂公告
查看: 4003  |  回复: 8

仲_______

金虫 (小有名气)

[求助] 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相差还是挺大的。请问应该这是怎么回事?谢谢!
回复此楼

» 猜你喜欢

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

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

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的回帖

仲_______

金虫 (小有名气)

送红花一朵
引用回帖:
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的回帖

valenhou001

至尊木虫 (职业作家)

引用回帖:
5楼: Originally posted by 仲_______ at 2014-07-18 17:23:04
好的,非常感谢。
我本科学的是通信,然后现在研一导师让我用vasp计算,所以一些基础不太懂。目前正在学习固体物理学和固体能带理论,希望慢慢讲基础补上去。
非常感谢您的回复!...

本科学通信的,来做这个材料计算,好像有点不值得。继续做通信的话,毕业后找工作要比这做材料计算容易得多。
6楼2014-07-18 18:55:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

仲_______

金虫 (小有名气)

引用回帖:
6楼: Originally posted by valenhou001 at 2014-07-18 18:55:55
本科学通信的,来做这个材料计算,好像有点不值得。继续做通信的话,毕业后找工作要比这做材料计算容易得多。...

导师跳到了自旋电子的方向上,所以实验室都做自旋电子相关的工作。有做相关的电路的、系统的、纠错的,导师让我做材料计算,好多东西没学,也只能硬着头皮上了。
7楼2014-07-18 19:14:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

valenhou001

至尊木虫 (职业作家)

引用回帖:
7楼: Originally posted by 仲_______ at 2014-07-18 20:14:03
导师跳到了自旋电子的方向上,所以实验室都做自旋电子相关的工作。有做相关的电路的、系统的、纠错的,导师让我做材料计算,好多东西没学,也只能硬着头皮上了。...

个人的看法,既然你这样的专业和背景,能逃掉材料计算,尽量避免。。。,选做其他的课题。
自旋电子学方面的材料计算现在已经不是2003年左右那么热了。计算能做的也大部分,也大部分做的是半导体或氧化物掺过渡金属,简单的做些磁性的计算,离指导实验还差很多。
8楼2014-07-18 20:18:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

绿色心情5977

铁虫 (初入文坛)

同问?
9楼2017-11-30 08:17:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 仲_______ 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 有没有学校材料专业收跨调(一志愿085410) +3 momo(上岸版) 2026-04-06 6/300 2026-04-08 14:53 by 化学521
[考研] 278求调剂 +17 范婷娜 2026-04-07 18/900 2026-04-08 14:43 by lijunpoly
[考研] 372求调剂 +4 jj涌77 2026-04-02 4/200 2026-04-07 09:31 by 白云朵朵飞
[考研] 269电子信息求调剂,可转专业 +5 独酌wl 2026-04-06 5/250 2026-04-06 17:23 by 土木硕士招生
[考研] 324求调剂 +3 k可乐 2026-04-05 4/200 2026-04-06 09:54 by 蓝云思雨
[考研] 一志愿哈尔滨工业大学085600英一数二337分求调剂 +12 lyz0427 2026-04-03 12/600 2026-04-06 06:37 by houyaoxu
[考研] 求调剂 +10 Hll胡 2026-04-04 10/500 2026-04-05 20:09 by nepu_uu
[考研] 315求调剂 +5 &123456789 2026-04-05 5/250 2026-04-05 19:55 by nepu_uu
[考研] 313求调剂 +5 海日海日 2026-04-04 7/350 2026-04-05 13:58 by imissbao
[考研] 一志愿江南大学085501机械工程专硕326分,本科佳木斯大学 +5 顾若浮生 2026-04-03 9/450 2026-04-05 09:57 by 1753564080
[考研] 材料与化工306分找调剂 +23 沧海轻舟e 2026-04-02 27/1350 2026-04-04 21:52 by laoshidan
[考研] 一志愿沪9,求生物学调剂,326分 +6 刘墨墨 2026-04-04 6/300 2026-04-04 19:44 by 唐沐儿
[考研] 368求调剂 +5 今华习 2026-04-03 7/350 2026-04-04 18:47 by imissbao
[考研] 复试调剂 +6 范根培 2026-04-04 6/300 2026-04-04 14:27 by 土木硕士招生
[考研] 一志愿北京科技大学材料工程085601,求调剂 +17 cdyw 2026-04-02 18/900 2026-04-04 11:14 by w_xuqing
[考研] 考研调剂 +8 不爱喝饮料 2026-04-03 8/400 2026-04-03 16:40 by Mistake-J
[考研] 材料考研调剂 +10 Gs大王 2026-04-02 10/500 2026-04-03 09:47 by 遗忘消失的灆
[考研] 26考研调剂 +4 Wnz.20030617 2026-04-01 5/250 2026-04-02 16:11 by 1939136013狗壮
[考研] 材料求调剂 +10 呢呢妮妮 2026-04-01 13/650 2026-04-02 09:17 by olim
[考研] 353求调剂 +4 拉钩不许变 2026-04-01 4/200 2026-04-01 18:10 by 记事本2026
信息提示
请填处理意见