| 5 | 1/1 | 返回列表 |
| 查看: 3965 | 回復(fù): 8 | ||
| 當(dāng)前只顯示滿足指定條件的回帖,點擊這里查看本話題的所有回帖 | ||
[求助]
VASP計算中KPOINTS文件的作用? 已有1人參與
|
||
|
大家好,本人初學(xué)vasp,對KPOINTS文件的作用不太理解。有人能比較形象地解釋一下嗎? 此外,今天做了兩次bulk Fe的靜態(tài)自洽計算,四個輸入文件中只有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)該這是怎么回事?謝謝! |
至尊木蟲 (職業(yè)作家)
|
你計算的體材料,而且采用的初基原胞,當(dāng)然不能用單k點(Gamma點)來進(jìn)行計算。兩次計算得到的總能當(dāng)然差別很大。 強(qiáng)烈建議你在做計算之前讀讀計算方法相關(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 |
至尊木蟲 (職業(yè)作家)
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
|
百度網(wǎng)盤 |
360云盤 |
千易網(wǎng)盤 |
華為網(wǎng)盤
在新窗口頁面中打開自己喜歡的網(wǎng)盤網(wǎng)站,將文件上傳后,然后將下載鏈接復(fù)制到帖子內(nèi)容中就可以了。 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 招收調(diào)劑碩士 +3 | lidianxing 2026-03-19 | 9/450 |
|
|---|---|---|---|---|
|
[考研] 一志愿 西北大學(xué) ,070300化學(xué)學(xué)碩,總分287,雙非一本,求調(diào)劑。 +3 | 晨昏線與星海 2026-03-19 | 3/150 |
|
|
[考研] 0817調(diào)劑 +3 | 沒有答案_ 2026-03-14 | 3/150 |
|
|
[考研] 能源材料化學(xué)課題組招收碩士研究生8-10名 +4 | 脫穎而出 2026-03-16 | 11/550 |
|
|
[考研] 0703化學(xué)調(diào)劑,求各位老師收留 +10 | 秋有木北 2026-03-14 | 10/500 |
|
|
[考研] 344求調(diào)劑 +6 | knight344 2026-03-16 | 7/350 |
|
|
[考研] 295求調(diào)劑 +3 | 一志愿京區(qū)211 2026-03-18 | 5/250 |
|
|
[考研] 311求調(diào)劑 +11 | 冬十三 2026-03-15 | 12/600 |
|
|
[考研] 304求調(diào)劑 +12 | 小熊joy 2026-03-14 | 13/650 |
|
|
[考研] 工科材料085601 279求調(diào)劑 +6 | 困于星晨 2026-03-17 | 6/300 |
|
|
[考研] 26考研求調(diào)劑 +6 | 丶宏Sir 2026-03-13 | 6/300 |
|
|
[考研] 283求調(diào)劑 +3 | 聽風(fēng)就是雨; 2026-03-16 | 3/150 |
|
|
[考研] 085600調(diào)劑 +5 | 漾漾123sun 2026-03-12 | 6/300 |
|
|
[考研] 0856求調(diào)劑 +3 | 劉夢微 2026-03-15 | 3/150 |
|
|
[考研] 求老師收留調(diào)劑 +4 | jiang姜66 2026-03-14 | 5/250 |
|
|
[考研] 080500,材料學(xué)碩302分求調(diào)劑學(xué)校 +4 | 初識可樂 2026-03-14 | 5/250 |
|
|
[考研] 297求調(diào)劑 +4 | 學(xué)海漂泊 2026-03-13 | 4/200 |
|
|
[考研] 招收0805(材料)調(diào)劑 +3 | 18595523086 2026-03-13 | 3/150 |
|
|
[考研] [0860]321分求調(diào)劑,ab區(qū)皆可 +4 | 寶貴熱 2026-03-13 | 4/200 |
|
|
[考研] 一志愿山大07化學(xué) 332分 四六級已過 本科山東雙非 求調(diào)劑! +3 | 不想理你 2026-03-12 | 3/150 |
|