當(dāng)前只顯示滿(mǎn)足指定條件的回帖，點(diǎn)擊這里查看本話(huà)題的所有回帖

仲_______

應(yīng)助: 1 (幼兒園)
金幣: 831.2
散金: 36
紅花: 71
帖子: 295
在線(xiàn): 689小時(shí)
蟲(chóng)號(hào): 2866770
注冊(cè): 2013-12-12
性別: GG
專(zhuān)業(yè): 新型信息器件

[求助] VASP計(jì)算中KPOINTS文件的作用？已有1人參與

大家好，本人初學(xué)vasp，對(duì)KPOINTS文件的作用不太理解。有人能比較形象地解釋一下嗎？

此外，今天做了兩次bulk Fe的靜態(tài)自洽計(jì)算，四個(gè)輸入文件中只有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相差還是挺大的。請(qǐng)問(wèn)應(yīng)該這是怎么回事？謝謝！

回復(fù)此樓

» 猜你喜歡

【考研調(diào)劑】化學(xué)專(zhuān)業(yè) 281分，一志愿四川大學(xué)，誠(chéng)心求調(diào)劑已經(jīng)有4人回復(fù)
317求調(diào)劑已經(jīng)有5人回復(fù)
0817 化學(xué)工程 299分求調(diào)劑有科研經(jīng)歷有二區(qū)文章已經(jīng)有10人回復(fù)
一志愿西北大學(xué) ，070300化學(xué)學(xué)碩，總分287，雙非一本，求調(diào)劑。已經(jīng)有3人回復(fù)
一志愿福大288有機(jī)化學(xué)，求調(diào)劑已經(jīng)有3人回復(fù)
085600材料與化工求調(diào)劑已經(jīng)有6人回復(fù)
一志愿西安交通大學(xué)材料工程專(zhuān)業(yè) 282分求調(diào)劑已經(jīng)有6人回復(fù)
085600材料與化工調(diào)劑 324分已經(jīng)有10人回復(fù)
294求調(diào)劑材料與化工專(zhuān)碩已經(jīng)有11人回復(fù)
一志愿南昌大學(xué)，327分，材料與化工085600 已經(jīng)有3人回復(fù)

» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:

Vasp在計(jì)算二維材料如石墨烯的時(shí)候怎么設(shè)置POSCAR和其他輸入文件？已經(jīng)有6人回復(fù)
vasp KPOINTS值如何選擇已經(jīng)有5人回復(fù)
vasp計(jì)算能帶，能量范圍的選取問(wèn)題已經(jīng)有5人回復(fù)
VASP計(jì)算態(tài)密度時(shí)K點(diǎn)選擇已經(jīng)有6人回復(fù)
用VASP計(jì)算單點(diǎn)能的時(shí)候總是出錯(cuò)，怎么辦�。看髠b來(lái)給知道啊？已經(jīng)有8人回復(fù)
【求助】vasp中計(jì)算石墨烯能帶結(jié)構(gòu)K點(diǎn)路徑該如何取已經(jīng)有11人回復(fù)
vasp 計(jì)算 partial charge 時(shí)k點(diǎn)選取的問(wèn)題已經(jīng)有19人回復(fù)
VASP5.3 k-points parallelism計(jì)算, KPAR的設(shè)置，KPAR=2 已經(jīng)有11人回復(fù)
關(guān)于VASP計(jì)算得到CHGCAR文件的數(shù)據(jù)處理以及數(shù)據(jù)單位的已經(jīng)有16人回復(fù)
vasp結(jié)構(gòu)優(yōu)化中EDIFF以及KPOINTS設(shè)置，以及氫吸附結(jié)合能問(wèn)題！已經(jīng)有17人回復(fù)
VASP中如何創(chuàng)建 INCAR , POSCAR, KPOINTS和POTCAR文件已經(jīng)有7人回復(fù)
VASP計(jì)算能帶時(shí)是否讀取WAVECAR 已經(jīng)有16人回復(fù)
vasp計(jì)算過(guò)程中總得把CONTCAR改成POSCAR再次計(jì)算已經(jīng)有7人回復(fù)
VASP 5.2.12 以上可以直接在INCAR中設(shè)置K點(diǎn) 已經(jīng)有9人回復(fù)
vasp計(jì)算中，如何選擇正確的K點(diǎn) 已經(jīng)有7人回復(fù)
VASP計(jì)算：LORBIT設(shè)置出現(xiàn)問(wèn)題已經(jīng)有5人回復(fù)
能帶計(jì)算 vasp 已經(jīng)有11人回復(fù)
VASP計(jì)算時(shí)的幾個(gè)細(xì)節(jié)設(shè)置問(wèn)題，希望有經(jīng)驗(yàn)的人過(guò)來(lái)指點(diǎn)一下！已經(jīng)有23人回復(fù)
【求助】VASP計(jì)算能帶時(shí)KPOINTS設(shè)置出錯(cuò) 已經(jīng)有14人回復(fù)
【求助】VASP計(jì)算的KPOINTS 設(shè)置，求助已經(jīng)有10人回復(fù)
【求助】VASP中K點(diǎn)設(shè)置已經(jīng)有14人回復(fù)

1樓 2014-07-18 11:45:29

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

valenhou001

至尊木蟲(chóng) (職業(yè)作家)

1ST強(qiáng)帖: 13
應(yīng)助: 241 (大學(xué)生)
金幣: 25701.7
散金: 602
紅花: 166
帖子: 3782
在線(xiàn): 873.8小時(shí)
蟲(chóng)號(hào): 1007127
注冊(cè): 2010-04-27
專(zhuān)業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

引用回帖:

3樓: Originally posted by 仲_______ at 2014-07-18 16:44:57
您好，非常感謝您的回答和文獻(xiàn)！
我在學(xué)習(xí)vasp官網(wǎng)給的教程“hand on session I”,其中提到“for molecules and atoms only a single k-point is required”. 然后在例子中，對(duì)于O原子、O2分子、CO分子、H2O分子都 ...

這是對(duì)這樣的分子計(jì)算，采用超原胞，而且取的很大，那對(duì)應(yīng)的布里淵區(qū)就很小，取Gamma點(diǎn)（用一個(gè)k點(diǎn)）取樣就足夠了。

你可能需要看看固體物理理解一下相關(guān)的基本知識(shí)（晶體、基矢、布里淵區(qū)、波矢、原胞等等）。

贊一下

回復(fù)此樓

» 本帖已獲得的紅花（最新10朵）

仲_______

4樓2014-07-18 16:10:16

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

查看全部 9 個(gè)回答

valenhou001

至尊木蟲(chóng) (職業(yè)作家)

1ST強(qiáng)帖: 13
應(yīng)助: 241 (大學(xué)生)
金幣: 25701.7
散金: 602
紅花: 166
帖子: 3782
在線(xiàn): 873.8小時(shí)
蟲(chóng)號(hào): 1007127
注冊(cè): 2010-04-27
專(zhuān)業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

【答案】應(yīng)助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感謝參與，應(yīng)助指數(shù) +1
仲_______: 金幣+10, ★★★★★最佳答案 2014-07-18 15:41:41

你計(jì)算的體材料，而且采用的初基原胞，當(dāng)然不能用單k點(diǎn)（Gamma點(diǎn)）來(lái)進(jìn)行計(jì)算。兩次計(jì)算得到的總能當(dāng)然差別很大。

強(qiáng)烈建議你在做計(jì)算之前讀讀計(jì)算方法相關(guān)的文獻(xiàn)，比如k-point sampling之類(lèi)的。另外理解布里淵區(qū)的概念。有些教程里給出了k-point sampling的簡(jiǎn)單圖示。
https://www.vasp.at/vasp-workshop/slides/k-points.pdf
http://www.jncasr.ac.in/ccms/fps ... day1/waghmare_1.pdf

贊一下

回復(fù)此樓

2樓2014-07-18 13:00:54

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

仲_______

金蟲(chóng) (小有名氣)

應(yīng)助: 1 (幼兒園)
金幣: 831.2
散金: 36
紅花: 71
帖子: 295
在線(xiàn): 689小時(shí)
蟲(chóng)號(hào): 2866770
注冊(cè): 2013-12-12
性別: GG
專(zhuān)業(yè): 新型信息器件

引用回帖:

2樓: Originally posted by valenhou001 at 2014-07-18 13:00:54
你計(jì)算的體材料，而且采用的初基原胞，當(dāng)然不能用單k點(diǎn)（Gamma點(diǎn)）來(lái)進(jìn)行計(jì)算。兩次計(jì)算得到的總能當(dāng)然差別很大。

強(qiáng)烈建議你在做計(jì)算之前讀讀計(jì)算方法相關(guān)的文獻(xiàn)，比如k-point sampling之類(lèi)的。另外理解布里淵區(qū) ...

您好，非常感謝您的回答和文獻(xiàn)！
我在學(xué)習(xí)vasp官網(wǎng)給的教程“hand on session I”,其中提到“for molecules and atoms only a single k-point is required”. 然后在例子中，對(duì)于O原子、O2分子、CO分子、H2O分子都采用的是我上文中提到的單K點(diǎn)的KPOINTS。請(qǐng)問(wèn)這個(gè)是為什么呢？謝謝！

贊一下

回復(fù)此樓

3樓2014-07-18 15:44:57

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

仲_______

金蟲(chóng) (小有名氣)

應(yīng)助: 1 (幼兒園)
金幣: 831.2
散金: 36
紅花: 71
帖子: 295
在線(xiàn): 689小時(shí)
蟲(chóng)號(hào): 2866770
注冊(cè): 2013-12-12
性別: GG
專(zhuān)業(yè): 新型信息器件

送紅花一朵

引用回帖:

4樓: Originally posted by valenhou001 at 2014-07-18 16:10:16
這是對(duì)這樣的分子計(jì)算，采用超原胞，而且取的很大，那對(duì)應(yīng)的布里淵區(qū)就很小，取Gamma點(diǎn)（用一個(gè)k點(diǎn)）取樣就足夠了。

你可能需要看看固體物理理解一下相關(guān)的基本知識(shí)（晶體、基矢、布里淵區(qū)、波矢、原胞等等）。...

好的，非常感謝。
我本科學(xué)的是通信，然后現(xiàn)在研一導(dǎo)師讓我用vasp計(jì)算，所以一些基礎(chǔ)不太懂。目前正在學(xué)習(xí)固體物理學(xué)和固體能帶理論，希望慢慢講基礎(chǔ)補(bǔ)上去。
非常感謝您的回復(fù)！

贊一下

回復(fù)此樓

5樓2014-07-18 16:23:04

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

最具人氣熱帖推薦 [查看全部]	作者	回/看	最后發(fā)表

[考研] 一志愿南昌大學(xué)，327分，材料與化工085600 +3	Ncdx123456 2026-03-19	3/150	2026-03-19 13:18 by houyaoxu
[考研] 一志愿武漢理工材料工程專(zhuān)碩調(diào)劑 +3	Doleres 2026-03-19	3/150	2026-03-19 12:29 by peike
[考研] 281求調(diào)劑（0805） +9	煙汐憶海 2026-03-16	19/950	2026-03-19 11:42 by laoshidan
[考研] 085410人工智能專(zhuān)碩317求調(diào)劑（0854都可以） +3	xbxudjdn 2026-03-18	3/150	2026-03-18 22:14 by zhq0425
[考研] 材料專(zhuān)業(yè)求調(diào)劑 +5	hanamiko 2026-03-18	5/250	2026-03-18 20:19 by 星空星月
[考研] 344求調(diào)劑 +6	knight344 2026-03-16	7/350	2026-03-18 20:13 by walc
[考研] 一志愿中國(guó)海洋大學(xué)，生物學(xué)，301分，求調(diào)劑 +4	1孫悟空 2026-03-17	4/200	2026-03-18 17:59 by fivewind
[考研] 材料與化工一志愿南昌大學(xué)327求調(diào)劑推薦 +8	Ncdx123456 2026-03-13	9/450	2026-03-18 14:40 by haxia
[考研] 收復(fù)試調(diào)劑生 +4	雨后秋荷 2026-03-18	4/200	2026-03-18 14:16 by elevennnne
[考研] 331求調(diào)劑（0703有機(jī)化學(xué) +7	ZY-05 2026-03-13	8/400	2026-03-18 14:13 by 007_lilei
[考研] 293求調(diào)劑 +11	zjl的號(hào) 2026-03-16	16/800	2026-03-18 08:10 by zhukairuo
[考研] 一志愿，福州大學(xué)材料專(zhuān)碩339分求調(diào)劑 +3	木子momo青爭(zhēng) 2026-03-15	3/150	2026-03-17 07:52 by laoshidan
[考研] 277材料科學(xué)與工程080500求調(diào)劑 +3	自由煎餅果子 2026-03-16	3/150	2026-03-16 14:10 by 運(yùn)氣yunqi
[考研] 求老師收留調(diào)劑 +4	jiang姜66 2026-03-14	5/250	2026-03-15 20:11 by Winj1e
[考研] 070305求調(diào)劑 +3	mlpqaz03 2026-03-14	4/200	2026-03-15 11:04 by peike
[考研] 328求調(diào)劑 +3	5201314Lsy！ 2026-03-13	6/300	2026-03-14 15:31 by hyswxzs
[考研] 一志愿哈工大材料324分求調(diào)劑 +5	閆旭東 2026-03-14	5/250	2026-03-14 14:53 by 木瓜膏
[考研] 266求調(diào)劑 +4	學(xué)員97LZgn 2026-03-13	4/200	2026-03-14 08:37 by zhukairuo
[考研] 招收0805（材料）調(diào)劑 +3	18595523086 2026-03-13	3/150	2026-03-14 00:33 by 123%、
[考研] 290求調(diào)劑 +3	ADT 2026-03-13	3/150	2026-03-13 10:19 by peike