| 5 | 1/1 | 返回列表 |
| 查看: 1773 | 回復(fù): 5 | |||||
| 當(dāng)前只顯示滿足指定條件的回帖,點(diǎn)擊這里查看本話題的所有回帖 | |||||
小小鑫000鐵蟲 (小有名氣)
|
[求助]
關(guān)于VASP中對(duì)ENCUT收斂性測(cè)試的疑問(wèn) 已有1人參與
|
||||
|
各位高手大家好: 有個(gè)問(wèn)題一直困擾我很多天了,是這樣的。我近期用vasp對(duì)Co進(jìn)行模擬計(jì)算。首先要確定Co的ENCUT和KPOINT大小,但是文獻(xiàn)上沒(méi)告訴,所以得自己先進(jìn)行收斂性測(cè)試。資料上說(shuō)先要假設(shè)一個(gè)ENCUT值,然后進(jìn)行KPOINT的收斂性測(cè)試,疑問(wèn)出來(lái)了,那就是假設(shè)的ENCUT應(yīng)該取多大?這是第一個(gè)疑問(wèn)。 我當(dāng)時(shí)就取了PONTCAR中的ENMAX=270,接著進(jìn)行了K點(diǎn)的收斂性測(cè)試,最后確定為8x8x8,然后用這個(gè)KPOINT再進(jìn)行ENCUT的收斂性測(cè)試,得到這個(gè)結(jié)果: 270 -27.073083 300 -27.112801 330 -27.106534 360 -27.089872 390 -27.080460 420 -27.077593 沒(méi)有出現(xiàn)總能變化在0.001,而且420這已經(jīng)是一個(gè)大的不能再大的值了,所以420后面的值更不可能收斂了,我的腳本也應(yīng)該沒(méi)問(wèn)題呀! #!/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 到底哪兒出問(wèn)題了,求高手解釋呀!不勝感激! [ Last edited by 小小鑫000 on 2013-11-27 at 10:32 ] |

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

專家顧問(wèn) (著名寫手)
![]() |
專家經(jīng)驗(yàn): +20 |


| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 313求調(diào)劑 +3 | 賣個(gè)關(guān)子吧 2026-03-31 | 3/150 |
|
|---|---|---|---|---|
|
[考研] 材料科學(xué)與工程調(diào)劑 +10 | 深V宿舍吧 2026-03-30 | 11/550 |
|
|
[考研] 一志愿 南京航空航天大學(xué) ,080500材料科學(xué)與工程學(xué)碩 +4 | @taotao 2026-03-30 | 4/200 |
|
|
[考研] 理學(xué)07化學(xué) 303求調(diào)劑 +4 | 睿08 2026-03-27 | 4/200 |
|
|
[考研] 317求調(diào)劑 +8 | 十閑wx 2026-03-24 | 8/400 |
|
|
[考研] 調(diào)劑 +4 | GK72 2026-03-30 | 4/200 |
|
|
[考研] 332求調(diào)劑 +14 | 032500 2026-03-25 | 14/700 |
|
|
[考研] 342求調(diào)劑 +4 | 加油a李zs 2026-03-26 | 4/200 |
|
|
[考研] 298求調(diào)劑 +3 | 什么是胖頭魚 2026-03-30 | 5/250 |
|
|
[考研] 284求調(diào)劑 +14 | junqihahaha 2026-03-26 | 15/750 |
|
|
[考研] 一志愿南航 335分 | 0856 | GPA 4.07 | 有科研經(jīng)歷 +8 | cccchenso 2026-03-29 | 8/400 |
|
|
[考研] 356求調(diào)劑 +3 | gysy?s?a 2026-03-28 | 3/150 |
|
|
[考研] 283求調(diào)劑 +3 | A child 2026-03-28 | 3/150 |
|
|
[考研] 347求調(diào)劑 +3 | 山頂見α 2026-03-25 | 3/150 |
|
|
[考研] 299求調(diào)劑 +7 | 嗯嗯嗯嗯2 2026-03-27 | 7/350 |
|
|
[考研] 0856調(diào)劑 +5 | 求求讓我有書讀?/a> 2026-03-26 | 6/300 |
|
|
[考研] 考研調(diào)劑 +9 | 小蠟新筆 2026-03-26 | 9/450 |
|
|
[考研] 一志愿天津大學(xué)339材料與化工求調(diào)劑 +3 | 江往賣魚 2026-03-26 | 3/150 |
|
|
[考研] 334分 一志愿武理-080500 材料求調(diào)劑 +4 | 李李不服輸 2026-03-25 | 4/200 |
|
|
[考研] 347求調(diào)劑 +4 | L when 2026-03-25 | 4/200 |
|