| 9 | 1/1 | 返回列表 |
| 查看: 3966 | 回復(fù): 8 | |||
[求助]
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)問應(yīng)該這是怎么回事?謝謝! |
至尊木蟲 (職業(yè)作家)
|
你計(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的簡(jiǎn)單圖示。 https://www.vasp.at/vasp-workshop/slides/k-points.pdf http://www.jncasr.ac.in/ccms/fps ... day1/waghmare_1.pdf |
至尊木蟲 (職業(yè)作家)
至尊木蟲 (職業(yè)作家)
至尊木蟲 (職業(yè)作家)
鐵蟲 (初入文壇)
| 9 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 材料與化工求調(diào)劑 +7 | 為學(xué)666 2026-03-16 | 7/350 |
|
|---|---|---|---|---|
|
[考研] 287求調(diào)劑 +3 | 晨昏線與星海 2026-03-19 | 4/200 |
|
|
[考研] 一志愿天大材料與化工(085600)總分338 +5 | 蔡大美女 2026-03-13 | 5/250 |
|
|
[考研] 0703化學(xué)調(diào)劑,求各位老師收留 +10 | 秋有木北 2026-03-14 | 10/500 |
|
|
[考研] 304求調(diào)劑 +6 | 司空. 2026-03-18 | 6/300 |
|
|
[考研] 085410人工智能專碩317求調(diào)劑(0854都可以) +3 | xbxudjdn 2026-03-18 | 3/150 |
|
|
[考研] 材料專業(yè)求調(diào)劑 +5 | hanamiko 2026-03-18 | 5/250 |
|
|
[考研] 一志愿中國(guó)海洋大學(xué),生物學(xué),301分,求調(diào)劑 +4 | 1孫悟空 2026-03-17 | 4/200 |
|
|
[考研] 26調(diào)劑/材料/英一數(shù)二/總分289/已過A區(qū)線 +7 | 步川酷紫123 2026-03-13 | 7/350 |
|
|
[考研] 0854可跨調(diào)劑,一作一項(xiàng)核心論文五項(xiàng)專利,省、國(guó)級(jí)證書40+數(shù)一英一287 +8 | 小李0854 2026-03-16 | 8/400 |
|
|
[考研] 070300化學(xué)319求調(diào)劑 +6 | 錦鯉0909 2026-03-17 | 6/300 |
|
|
[考研] 303求調(diào)劑 +4 | 睿08 2026-03-17 | 6/300 |
|
|
[碩博家園] 湖北工業(yè)大學(xué) 生命科學(xué)與健康學(xué)院-課題組招收2026級(jí)食品/生物方向碩士 +3 | 1喜春8 2026-03-17 | 5/250 |
|
|
[考研] 302求調(diào)劑 +4 | 小賈同學(xué)123 2026-03-15 | 8/400 |
|
|
[考研] 機(jī)械專碩325,尋找調(diào)劑院校 +3 | y9999 2026-03-15 | 5/250 |
|
|
[考研] 一志愿211 0703方向310分求調(diào)劑 +3 | 努力奮斗112 2026-03-15 | 3/150 |
|
|
[考研] 277材料科學(xué)與工程080500求調(diào)劑 +3 | 自由煎餅果子 2026-03-16 | 3/150 |
|
|
[考研] 26考研一志愿中國(guó)石油大學(xué)(華東)305分求調(diào)劑 +3 | 嘉年新程 2026-03-15 | 3/150 |
|
|
[考研] 中科大材料專碩319求調(diào)劑 +3 | 孟鑫材料 2026-03-13 | 3/150 |
|
|
[考研] 招收0805(材料)調(diào)劑 +3 | 18595523086 2026-03-13 | 3/150 |
|