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

仲_______

金蟲 (小有名氣)

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

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

大家好，本人初學(xué)vasp，對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相差還是挺大的。請問應(yīng)該這是怎么回事？謝謝！

回復(fù)此樓

» 猜你喜歡

08工科，295，接受跨專業(yè)調(diào)劑已經(jīng)有4人回復(fù)
318求調(diào)劑已經(jīng)有8人回復(fù)
307分求調(diào)劑已經(jīng)有4人回復(fù)
288資源與環(huán)境專碩求調(diào)劑，不限專業(yè)，有學(xué)上就行已經(jīng)有17人回復(fù)
085601 329分調(diào)劑已經(jīng)有6人回復(fù)
343求調(diào)劑已經(jīng)有6人回復(fù)
070300化學(xué)求調(diào)劑已經(jīng)有10人回復(fù)
學(xué)碩274求調(diào)劑已經(jīng)有17人回復(fù)
346求調(diào)劑一志愿070303有機(jī)化學(xué) 已經(jīng)有11人回復(fù)
085601 材料工程 313分求調(diào)劑已經(jīng)有6人回復(fù)

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

Vasp在計(jì)算二維材料如石墨烯的時(shí)候怎么設(shè)置POSCAR和其他輸入文件？已經(jīng)有6人回復(fù)
vasp KPOINTS值如何選擇已經(jīng)有5人回復(fù)
vasp計(jì)算能帶，能量范圍的選取問題已經(jīng)有5人回復(fù)
VASP計(jì)算態(tài)密度時(shí)K點(diǎn)選擇已經(jīng)有6人回復(fù)
用VASP計(jì)算單點(diǎn)能的時(shí)候總是出錯(cuò)，怎么辦啊？大俠來給知道��？已經(jīng)有8人回復(fù)
【求助】vasp中計(jì)算石墨烯能帶結(jié)構(gòu)K點(diǎn)路徑該如何取已經(jīng)有11人回復(fù)
vasp 計(jì)算 partial charge 時(shí)k點(diǎ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é)合能問題！已經(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ì)算過程中總得把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)問題已經(jīng)有5人回復(fù)
能帶計(jì)算 vasp 已經(jīng)有11人回復(fù)
VASP計(jì)算時(shí)的幾個(gè)細(xì)節(jié)設(shè)置問題，希望有經(jīng)驗(yàn)的人過來指點(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

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

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

引用回帖:

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

本科學(xué)通信的，來做這個(gè)材料計(jì)算，好像有點(diǎn)不值得。繼續(xù)做通信的話，畢業(yè)后找工作要比這做材料計(jì)算容易得多。

贊一下

回復(fù)此樓

6樓2014-07-18 18:55:55

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

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

valenhou001

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

1ST強(qiáng)帖: 13
應(yīng)助: 241 (大學(xué)生)
金幣: 25701.7
散金: 602
紅花: 166
帖子: 3782
在線: 873.8小時(shí)
蟲號: 1007127
注冊: 2010-04-27
專業(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)）來進(jìn)行計(jì)算。兩次計(jì)算得到的總能當(dāng)然差別很大。

強(qiáng)烈建議你在做計(jì)算之前讀讀計(jì)算方法相關(guān)的文獻(xiàn)，比如k-point sampling之類的。另外理解布里淵區(qū)的概念。有些教程里給出了k-point sampling的簡單圖示。
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的回帖

仲_______

金蟲 (小有名氣)

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

引用回帖:

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

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

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

贊一下

回復(fù)此樓

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

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

valenhou001

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

1ST強(qiáng)帖: 13
應(yīng)助: 241 (大學(xué)生)
金幣: 25701.7
散金: 602
紅花: 166
帖子: 3782
在線: 873.8小時(shí)
蟲號: 1007127
注冊: 2010-04-27
專業(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”. 然后在例子中，對于O原子、O2分子、CO分子、H2O分子都 ...

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

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

贊一下

回復(fù)此樓

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

仲_______

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

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

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

[考研] 085601 329分調(diào)劑 +6	yzsa12 2026-03-31	6/300	2026-03-31 15:23 by yanflower7133
[考研] 289求調(diào)劑 +3	Acesczlo 2026-03-29	4/200	2026-03-31 14:48 by 熱情沙漠
[考研] 本2一志愿C9-333分，材料科學(xué)與工程，求調(diào)劑 +6	升升不降 2026-03-31	6/300	2026-03-31 14:33 by 記事本2026
[考研] 考研調(diào)劑求助 +5	13287130938 2026-03-31	5/250	2026-03-31 14:16 by limeifeng
[考研] 290分調(diào)劑求助 +9	吉祥止止陳 2026-03-25	9/450	2026-03-31 14:16 by 六翼TMR
[考研] 環(huán)境工程297分求調(diào)劑一志愿杭高院 +5	GENJIOW 2026-03-31	5/250	2026-03-31 13:58 by limeifeng
[考研] 材料專碩調(diào)劑 +12	椰椰。 2026-03-29	12/600	2026-03-31 10:52 by foria
[考研] 南京大學(xué)化學(xué)調(diào)劑 +11	景隨風(fēng) 2026-03-29	16/800	2026-03-31 10:14 by herarysara
[考研] 297 地理學(xué)070500 復(fù)試求調(diào)劑 +3	小圓圈圈ooo 2026-03-30	3/150	2026-03-30 21:05 by 余震yz
[考研] 一志愿：西北大學(xué)，英一數(shù)一408-284分求調(diào)劑 +5	12.27 2026-03-27	5/250	2026-03-30 19:31 by 源_2020
[考研] 求化學(xué)調(diào)劑 +11	wulanna 2026-03-28	11/550	2026-03-30 10:59 by 探123
[考研] 318一志愿吉林大學(xué)生物與醫(yī)藥求調(diào)劑 +5	篤行致遠(yuǎn). 2026-03-28	5/250	2026-03-30 06:56 by ilovexiaobin
[考研] 290求調(diào)劑 +3	dfffsar 2026-03-29	3/150	2026-03-29 22:38 by 毛毛毛阿莫2
[考研] 【求調(diào)劑】085601材料工程專碩 \| 總分272 \| +7	腳滑的守法公民 2026-03-27	7/350	2026-03-29 20:21 by dophin1985
[考研] 329求調(diào)劑 +10	鈕恩雪 2026-03-25	10/500	2026-03-29 13:32 by peike
[考研] 279求調(diào)劑 +4	蝶舞輕繞 2026-03-29	4/200	2026-03-29 09:45 by laoshidan
[考研] 305求調(diào)劑 +8	RuiFairyrui 2026-03-28	8/400	2026-03-29 08:22 by fmesaito
[考研] 材料求調(diào)劑一志愿哈工大324 +7	閆旭東 2026-03-28	9/450	2026-03-28 08:51 by Xu de nuo
[考研] 復(fù)試調(diào)劑，一志愿南農(nóng)083200食品科學(xué)與工程 +5	XQTJZ 2026-03-26	5/250	2026-03-27 14:49 by 狂炫麥當(dāng)當(dāng)
[考研] 調(diào)劑推薦 +5	清酒714 2026-03-26	6/300	2026-03-27 11:12 by 不吃魚的貓