| 6 | 1/1 | 返回列表 |
| 查看: 1738 | 回復(fù): 5 | ||||
小小鑫000鐵蟲 (小有名氣)
|
[求助]
關(guān)于VASP中對ENCUT收斂性測試的疑問 已有1人參與
|
|
各位高手大家好: 有個(gè)問題一直困擾我很多天了,是這樣的。我近期用vasp對Co進(jìn)行模擬計(jì)算。首先要確定Co的ENCUT和KPOINT大小,但是文獻(xiàn)上沒告訴,所以得自己先進(jìn)行收斂性測試。資料上說先要假設(shè)一個(gè)ENCUT值,然后進(jìn)行KPOINT的收斂性測試,疑問出來了,那就是假設(shè)的ENCUT應(yīng)該取多大?這是第一個(gè)疑問。 我當(dāng)時(shí)就取了PONTCAR中的ENMAX=270,接著進(jìn)行了K點(diǎn)的收斂性測試,最后確定為8x8x8,然后用這個(gè)KPOINT再進(jìn)行ENCUT的收斂性測試,得到這個(gè)結(jié)果: 270 -27.073083 300 -27.112801 330 -27.106534 360 -27.089872 390 -27.080460 420 -27.077593 沒有出現(xiàn)總能變化在0.001,而且420這已經(jīng)是一個(gè)大的不能再大的值了,所以420后面的值更不可能收斂了,我的腳本也應(yīng)該沒問題呀! #!/bin/sh rm WAVECAR for i in 15 270 300 330 360 390 420 do cat > INCAR <<! SYSTEM = Co-Diamond ENCUT = $i ISTART = 0 ; ICHARG = 2 ISMEAR = -5 PREC = Accurate ! echo "ENCUT = $i eV" ; time vasp E=`grep "TOTEN" OUTCAR | tail -1 | awk '{printf "%12.6f \n", $5 }'` echo $i $E >>comment done 到底哪兒出問題了,求高手解釋呀!不勝感激! [ Last edited by 小小鑫000 on 2013-11-27 at 10:32 ] |

專家顧問 (著名寫手)
![]() |
專家經(jīng)驗(yàn): +20 |
專家顧問 (著名寫手)
![]() |
專家經(jīng)驗(yàn): +174 |
|
雖然一般可能更多是先假設(shè)K-mesh進(jìn)行ECUT的收斂測試,不過反過來也沒啥原則性錯(cuò)誤。 不過問題就在于,如果要假設(shè)ECUT,那就必然要假設(shè)個(gè)“一定能夠滿足收斂”的ECUT,否則得到的K-mesh又能有什么意義? 所以,這就需要把ECUT設(shè)置得比一般別人文獻(xiàn)報(bào)告的還要略高,才比較安全。如果我是LZ,可能ECUT至少要從500eV或者更高來假設(shè),然后對K-mesh進(jìn)行測試。之后在固定合理的k-mesh的情況下,回頭對ecut進(jìn)行假設(shè)。 話說,為什么"420這已經(jīng)是一個(gè)大得不能再大的值了"呢?很多用VASP計(jì)算的文獻(xiàn),給出的cut-off energy都是500eV左右的啊~ |



| 6 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 07化學(xué)280分求調(diào)劑 +4 | 722865 2026-03-23 | 4/200 |
|
|---|---|---|---|---|
|
[考研]
|
pk3725069 2026-03-19 | 16/800 |
|
|
[考研] 材料專碩英一數(shù)二306 +8 | z1z2z3879 2026-03-18 | 8/400 |
|
|
[考研] 工科材料085601 279求調(diào)劑 +8 | 困于星晨 2026-03-17 | 10/500 |
|
|
[考研] 333求調(diào)劑 +6 | 87639 2026-03-21 | 10/500 |
|
|
[考研] 352求調(diào)劑 +3 | 大米飯! 2026-03-22 | 3/150 |
|
|
[考研] 0854電子信息求調(diào)劑 +3 | α____ 2026-03-22 | 3/150 |
|
|
[考研] 石河子大學(xué)(211、雙一流)碩博研究生長期招生公告 +3 | 李子目 2026-03-22 | 3/150 |
|
|
[考研] 306求調(diào)劑 +5 | 來好運(yùn)來來來 2026-03-22 | 5/250 |
|
|
[考研] 289求調(diào)劑 +7 | 懷瑾握瑜l 2026-03-20 | 7/350 |
|
|
[考博] 招收博士1-2人 +3 | QGZDSYS 2026-03-18 | 4/200 |
|
|
[考研] 297求調(diào)劑 +3 | 喜歡還是不甘心 2026-03-20 | 3/150 |
|
|
[考研] 268求調(diào)劑 +9 | 簡單點(diǎn)0 2026-03-17 | 9/450 |
|
|
[考研] 332求調(diào)劑 +3 | 鳳凰院丁真 2026-03-20 | 3/150 |
|
|
[考研] 301求調(diào)劑 +10 | yy要上岸呀 2026-03-17 | 10/500 |
|
|
[考研] 一志愿 南京航空航天大學(xué)大學(xué) ,080500材料科學(xué)與工程學(xué)碩 +5 | @taotao 2026-03-20 | 5/250 |
|
|
[考研] 261求B區(qū)調(diào)劑,科研經(jīng)歷豐富 +3 | 牛奶很忙 2026-03-20 | 4/200 |
|
|
[考研] 材料學(xué)碩318求調(diào)劑 +5 | February_Feb 2026-03-19 | 5/250 |
|
|
[考研] 收復(fù)試調(diào)劑生 +4 | 雨后秋荷 2026-03-18 | 4/200 |
|
|
[考研] 考研求調(diào)劑 +3 | 橘頌. 2026-03-17 | 4/200 |
|