| 5 | 1/1 | 返回列表 |
| 查看: 3972 | 回復: 8 | ||
| 當前只顯示滿足指定條件的回帖,點擊這里查看本話題的所有回帖 | ||
[求助]
VASP計算中KPOINTS文件的作用? 已有1人參與
|
||
|
大家好,本人初學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è)作家)
至尊木蟲 (職業(yè)作家)
|
你計算的體材料,而且采用的初基原胞,當然不能用單k點(Gamma點)來進行計算。兩次計算得到的總能當然差別很大。 強烈建議你在做計算之前讀讀計算方法相關的文獻,比如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è)作家)
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 求調劑 +5 | 十三加油 2026-03-21 | 5/250 |
|
|---|---|---|---|---|
|
[考研] 一志愿深大,0703化學,總分302,求調劑 +4 | 七月-七七 2026-03-21 | 4/200 |
|
|
[考研] 280求調劑 +10 | 咕嚕曉曉 2026-03-18 | 11/550 |
|
|
[考研] 求調劑 +3 | .m.. 2026-03-21 | 4/200 |
|
|
[考研] 306求調劑 +4 | chuanzhu川燭 2026-03-18 | 4/200 |
|
|
[考研] 南昌大學材料專碩311分求調劑 +6 | 77chaselx 2026-03-20 | 6/300 |
|
|
[考研] 265求調劑 +3 | Jack?k?y 2026-03-17 | 3/150 |
|
|
[考研] 328求調劑,英語六級551,有科研經(jīng)歷 +4 | 生物工程調劑 2026-03-17 | 8/400 |
|
|
[考研] 材料 336 求調劑 +3 | An@. 2026-03-18 | 4/200 |
|
|
[考研] 一志愿華中科技大學,080502,354分求調劑 +5 | 守候夕陽CF 2026-03-18 | 5/250 |
|
|
[考研] 藥學383 求調劑 +3 | 藥學chy 2026-03-15 | 5/250 |
|
|
[考研] 260求調劑 +3 | 朱芷琳 2026-03-20 | 3/150 |
|
|
[考研] 材料學求調劑 +4 | Stella_Yao 2026-03-20 | 4/200 |
|
|
[考研] 求調劑 +3 | eation27 2026-03-20 | 3/150 |
|
|
[考研] 086500 325 求調劑 +3 | 領帶小熊 2026-03-19 | 3/150 |
|
|
[考研] 工科材料085601 279求調劑 +7 | 困于星晨 2026-03-17 | 9/450 |
|
|
[考研] 廣西大學家禽遺傳育種課題組2026年碩士招生(接收計算機專業(yè)調劑) +3 | 123阿標 2026-03-17 | 3/150 |
|
|
[考研] 招收調劑碩士 +4 | lidianxing 2026-03-19 | 12/600 |
|
|
[考研] 085601專碩,總分342求調劑,地區(qū)不限 +5 | share_joy 2026-03-16 | 5/250 |
|
|
[考研] 有沒有道鐵/土木的想調劑南林,給自己招師弟中~ +3 | TqlXswl 2026-03-16 | 7/350 |
|