| 查看: 7314 | 回復(fù): 15 | |||
xianglannihao銅蟲 (小有名氣)
|
[交流]
關(guān)于VASP收斂判斷
|
||
|
關(guān)于VASP收斂判斷: 假定設(shè)置如下: EDIFF = 1E-05 (SCF convergence criteria) EDIFFG = -0.03 (Ionic convergence criteria, eV/AA) 我的理解收斂條件:EDIFF 是電子步收斂的條件,電子步dE達到1E-5就認(rèn)為收斂 。 EDIFFG是離子步收斂的標(biāo)準(zhǔn),為負數(shù)時表示力的收斂判據(jù)。是否收斂看OUTCAR中的 reached required accuracy - stopping structural energy minimisation提示和最后一輪的如下兩列: d Force = 0.1423097E-03[ 0.108E-03, 0.177E-03] d Energy = 0.1531149E-03-0.108E-04 d Force =-0.6844216E+00[-0.684E+00,-0.684E+00] d Ewald =-0.6844216E+00-0.177E-07 問題:1. 是否有這個提示就表示收斂了? 2. OUTCAR中為什么有兩列d Force, 都表示什么意思?d Force d Energy表示的都是一個范圍? 3. 有時在每一輪的離子步計算時,其電子步dE沒達到也顯示該輪收斂。 4. 優(yōu)化好的構(gòu)型,以相同的設(shè)置再重新優(yōu)化一次,我發(fā)現(xiàn)其能量會發(fā)生變化,有時會大于1eV,但構(gòu)型變化倒不是很大?難道我們要一直優(yōu)化,直至能量不變,是否這樣才能得到可靠的計算結(jié)果? |
木蟲 (正式寫手)
|
難得的干貨,這些問題我也都遇到了……關(guān)于收斂標(biāo)準(zhǔn),我學(xué)的是在total force后有個total drift后的三列都小于力的收斂,另外出現(xiàn)計算時間,說明整個計算過程是正確運行的! 關(guān)于多次優(yōu)化的問題,我也會這樣,尤其當(dāng)你的構(gòu)形越接近最優(yōu)構(gòu)形,你優(yōu)化次數(shù)越多,偏差越大,甚至構(gòu)形還會散或不收斂……一直也沒解決,我理解的計算原理這樣解釋:軟件在計算時會先進行離子步移動,當(dāng)你接近最優(yōu)構(gòu)形時,如果離子步過大,超過了收斂標(biāo)準(zhǔn)的閾值,而離子步移動具有一定隨機性,因此可能導(dǎo)致更大的偏離最優(yōu)構(gòu)形! vasp本身最優(yōu)構(gòu)形就是相對的,精確的用聲子譜。優(yōu)化一段后最好減小離子步幅…… 當(dāng)然,以上都是我的個人理解,權(quán)當(dāng)拋磚引玉,坐等權(quán)威解釋…… 發(fā)自小木蟲Android客戶端 |
木蟲 (正式寫手)

銅蟲 (小有名氣)
木蟲 (正式寫手)

銅蟲 (小有名氣)
木蟲 (正式寫手)

|
銅蟲 (小有名氣)
|
謝謝您的討論。 我的測試結(jié)果如下 : 構(gòu)型1, 優(yōu)化一次,能量為-219.837378; 以第一次優(yōu)化得到的構(gòu)型,以相同設(shè)置再次優(yōu)化,能量為-220.206788; 依此再優(yōu)化能量為-219.12072; 構(gòu)型2, 優(yōu)化一次,能量為-220.292638; 以第一次優(yōu)化得到的構(gòu)型,以相同設(shè)置再次優(yōu)化,能量為-219.157869; 依此再優(yōu)化能量為-218.285171; 構(gòu)型1和模型分別是同種分子不同吸附構(gòu)型,相對穩(wěn)定是如果按第一次優(yōu)化的結(jié)果是構(gòu)型2更穩(wěn)定,而按后兩次優(yōu)化的結(jié)果是構(gòu)型1穩(wěn)定。而按理論的預(yù)計是應(yīng)該構(gòu)型1穩(wěn)定的。所以我在想到底是什么問題造成的差別,思考如何得到可靠結(jié)果的問題。 再次優(yōu)化一般又會優(yōu)化10F至20F左右. |
木蟲 (正式寫手)

| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 能源材料化學(xué)課題組招收碩士研究生8-10名 +4 | 脫穎而出 2026-03-16 | 8/400 |
|
|---|---|---|---|---|
|
[考研] 一志愿中國海洋大學(xué),生物學(xué),301分,求調(diào)劑 +3 | 1孫悟空 2026-03-17 | 3/150 |
|
|
[考研] 工科材料085601 279求調(diào)劑 +6 | 困于星晨 2026-03-17 | 6/300 |
|
|
[考研] 生物學(xué)071000 329分求調(diào)劑 +3 | 我愛生物生物愛?/a> 2026-03-17 | 3/150 |
|
|
[考研] 化學(xué)工程321分求調(diào)劑 +13 | 大米飯! 2026-03-15 | 16/800 |
|
|
[考研] 293求調(diào)劑 +11 | zjl的號 2026-03-16 | 16/800 |
|
|
[考研] 材料與化工求調(diào)劑 +6 | 為學(xué)666 2026-03-16 | 6/300 |
|
|
[考研] 328求調(diào)劑,英語六級551,有科研經(jīng)歷 +3 | 生物工程調(diào)劑 2026-03-16 | 8/400 |
|
|
[考研] 26考研求調(diào)劑 +6 | 丶宏Sir 2026-03-13 | 6/300 |
|
|
[考研] 085600材料與化工 +4 | 安全上岸! 2026-03-16 | 4/200 |
|
|
[考研] 一志愿南京大學(xué),080500材料科學(xué)與工程,調(diào)劑 +4 | Jy? 2026-03-16 | 4/200 |
|
|
[考研] 277材料科學(xué)與工程080500求調(diào)劑 +3 | 自由煎餅果子 2026-03-16 | 3/150 |
|
|
[考研] 0856專碩279求調(diào)劑 +5 | 加油加油!? 2026-03-15 | 5/250 |
|
|
[考研] 080500,材料學(xué)碩302分求調(diào)劑學(xué)校 +4 | 初識可樂 2026-03-14 | 5/250 |
|
|
[考研] 255求調(diào)劑 +3 | 李嘉慧, 2026-03-12 | 4/200 |
|
|
[考研] 304求調(diào)劑 +7 | 7712b 2026-03-13 | 7/350 |
|
|
[考研] 求調(diào)劑 +5 | 一定有學(xué)上- 2026-03-12 | 5/250 |
|
|
[考研] 0703化學(xué)求調(diào)劑 +7 | 綠豆芹菜湯 2026-03-12 | 7/350 |
|
|
[考研] 求b區(qū)學(xué)校調(diào)劑 +3 | 周56 2026-03-11 | 3/150 |
|
|
[考研] 一志愿山大07化學(xué) 332分 四六級已過 本科山東雙非 求調(diào)劑! +3 | 不想理你 2026-03-12 | 3/150 |
|