| 查看: 7315 | 回復(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達(dá)到1E-5就認(rèn)為收斂 。 EDIFFG是離子步收斂的標(biāo)準(zhǔn),為負(fù)數(shù)時(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. 是否有這個(gè)提示就表示收斂了? 2. OUTCAR中為什么有兩列d Force, 都表示什么意思?d Force d Energy表示的都是一個(gè)范圍? 3. 有時(shí)在每一輪的離子步計(jì)算時(shí),其電子步dE沒達(dá)到也顯示該輪收斂。 4. 優(yōu)化好的構(gòu)型,以相同的設(shè)置再重新優(yōu)化一次,我發(fā)現(xiàn)其能量會(huì)發(fā)生變化,有時(shí)會(huì)大于1eV,但構(gòu)型變化倒不是很大?難道我們要一直優(yōu)化,直至能量不變,是否這樣才能得到可靠的計(jì)算結(jié)果? |
木蟲 (正式寫手)
|
難得的干貨,這些問題我也都遇到了……關(guān)于收斂標(biāo)準(zhǔn),我學(xué)的是在total force后有個(gè)total drift后的三列都小于力的收斂,另外出現(xiàn)計(jì)算時(shí)間,說明整個(gè)計(jì)算過程是正確運(yùn)行的! 關(guān)于多次優(yōu)化的問題,我也會(huì)這樣,尤其當(dāng)你的構(gòu)形越接近最優(yōu)構(gòu)形,你優(yōu)化次數(shù)越多,偏差越大,甚至構(gòu)形還會(huì)散或不收斂……一直也沒解決,我理解的計(jì)算原理這樣解釋:軟件在計(jì)算時(shí)會(huì)先進(jìn)行離子步移動(dòng),當(dāng)你接近最優(yōu)構(gòu)形時(shí),如果離子步過大,超過了收斂標(biāo)準(zhǔn)的閾值,而離子步移動(dòng)具有一定隨機(jī)性,因此可能導(dǎo)致更大的偏離最優(yōu)構(gòu)形! vasp本身最優(yōu)構(gòu)形就是相對的,精確的用聲子譜。優(yōu)化一段后最好減小離子步幅…… 當(dāng)然,以上都是我的個(gè)人理解,權(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ù)計(jì)是應(yīng)該構(gòu)型1穩(wěn)定的。所以我在想到底是什么問題造成的差別,思考如何得到可靠結(jié)果的問題。 再次優(yōu)化一般又會(huì)優(yōu)化10F至20F左右. |
木蟲 (正式寫手)

| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[教師之家] 焦慮 +8 | 水冰月月野兔 2026-03-13 | 11/550 |
|
|---|---|---|---|---|
|
[考研] 302求調(diào)劑 +10 | 呼呼呼。。。。 2026-03-17 | 10/500 |
|
|
[考研] 266求調(diào)劑 +4 | 陽陽哇塞 2026-03-14 | 8/400 |
|
|
[考研] 288求調(diào)劑,一志愿華南理工大學(xué)071005 +4 | ioodiiij 2026-03-17 | 4/200 |
|
|
[考研] 0703化學(xué)調(diào)劑 ,六級已過,有科研經(jīng)歷 +9 | 曦熙兮 2026-03-15 | 9/450 |
|
|
[考研] 332求調(diào)劑 +6 | Zz版 2026-03-13 | 6/300 |
|
|
[考研] 085600材料與化工求調(diào)劑 +5 | 緒幸與子 2026-03-17 | 5/250 |
|
|
[考研] 一志愿蘇州大學(xué)材料工程(085601)專碩有科研經(jīng)歷三項(xiàng)國獎(jiǎng)兩個(gè)實(shí)用型專利一項(xiàng)省級立項(xiàng) +6 | 大火山小火山 2026-03-16 | 8/400 |
|
|
[考博] 26申博 +4 | 八旬速覽 2026-03-16 | 4/200 |
|
|
[考研] 材料專碩274一志愿陜西師范大學(xué)求調(diào)劑 +5 | 薛云鵬 2026-03-13 | 5/250 |
|
|
[考研] 一志愿,福州大學(xué)材料專碩339分求調(diào)劑 +3 | 木子momo青爭 2026-03-15 | 3/150 |
|
|
[考研] 333求調(diào)劑 +3 | 文思客 2026-03-16 | 7/350 |
|
|
[考研] 304求調(diào)劑 +5 | 素年祭語 2026-03-15 | 5/250 |
|
|
[考研] 326求調(diào)劑 +3 | mlpqaz03 2026-03-15 | 3/150 |
|
|
[考研] 327求調(diào)劑 +6 | 拾光任染 2026-03-15 | 11/550 |
|
|
[考研] 288求調(diào)劑 +4 | 奇點(diǎn)0314 2026-03-14 | 4/200 |
|
|
[考研] 308 085701 四六級已過求調(diào)劑 +7 | 溫喬喬喬喬 2026-03-12 | 14/700 |
|
|
[考研] 材料與化工求調(diào)劑一志愿 985 總分 295 +8 | dream…… 2026-03-12 | 8/400 |
|
|
[考研] 四川大學(xué)085601材料工程專碩 初試294求調(diào)劑 +4 | 祝我們好在冬天 2026-03-11 | 4/200 |
|
|
[考研] 材料專碩350 求調(diào)劑 +4 | 王金科 2026-03-12 | 4/200 |
|