| 16 | 1/1 | 返回列表 |
| 查看: 3930 | 回復(fù): 15 | ||||||
sunyang1988至尊木蟲 (著名寫手)
|
[交流]
【求助】VASP收斂問題
|
|||||
|
計(jì)算中選用的收斂標(biāo)準(zhǔn)是能量 EDIFF = 1E-5 EDIFFG = 0.0001 但運(yùn)算正常結(jié)束后,能量差值并沒有達(dá)到這個(gè)要求 free energy TOTEN = -462.008990 eV free energy TOTEN = -462.024537 eV free energy TOTEN = -462.031260 eV free energy TOTEN = -462.037495 eV 我查看了下原子受力情況,都已經(jīng)非常小了 POSITION TOTAL-FORCE (eV/Angst) ----------------------------------------------------------------------------------- 1.15865 1.15865 1.15865 0.000002 -0.000001 -0.000001 3.47596 3.47596 1.15865 0.000000 -0.000002 -0.000002 3.47596 1.15865 3.47596 0.000003 -0.000001 0.000001 1.15865 3.47596 3.47596 0.000001 0.000000 0.000000 1.15865 1.15865 3.47596 0.000001 -0.000001 0.000000 3.47596 3.47596 3.47596 0.000002 -0.000001 0.000002 3.47596 1.15865 1.15865 0.000000 0.000000 -0.000001 1.15865 3.47596 1.15865 0.000002 -0.000001 -0.000001 5.79326 1.15865 1.15865 -0.000002 -0.000002 0.000001 8.11056 3.47596 1.15865 0.000001 0.000001 0.000000 8.11056 1.15865 3.47596 -0.000001 -0.000002 -0.000001 5.79326 3.47596 3.47596 -0.000002 -0.000001 0.000001 5.79326 1.15865 3.47596 -0.000003 -0.000002 0.000001 8.11056 3.47596 3.47596 -0.000001 0.000000 0.000000 8.11056 1.15865 1.15865 0.000000 -0.000002 0.000000 5.79326 3.47596 1.15865 -0.000002 -0.000001 0.000000 1.15865 5.79326 1.15865 0.000001 0.000001 -0.000002 3.47596 8.11056 1.15865 0.000001 0.000001 0.000000 3.47596 5.79326 3.47596 0.000002 0.000003 0.000002 1.15865 8.11056 3.47596 0.000002 0.000000 0.000000 1.15865 5.79326 3.47596 0.000001 0.000001 0.000001 3.47596 8.11056 3.47596 0.000002 0.000000 0.000001 3.47596 5.79326 1.15865 0.000000 0.000002 0.000000 1.15865 8.11056 1.15865 0.000000 0.000001 -0.000001 5.79326 5.79326 1.15865 0.000000 0.000002 0.000001 8.11056 8.11056 1.15865 -0.000001 0.000000 -0.000001 8.11056 5.79326 3.47596 -0.000002 0.000002 0.000001 5.79326 8.11056 3.47596 -0.000002 0.000000 0.000000 5.79326 5.79326 3.47596 -0.000001 0.000003 0.000001 8.11056 8.11056 3.47596 -0.000002 0.000000 -0.000001 8.11056 5.79326 1.15865 -0.000001 0.000001 0.000000 不知道這種情況到底算不算收斂? |
材料計(jì)算模擬實(shí)用技巧 | VASP | vasp計(jì)算資料 |
» 搶金幣啦!回帖就可以得到:
+1/167
+2/114
+2/52
+1/40
+1/38
+1/37
+1/35
+1/32
+1/19
+1/7
+1/7
+1/6
+1/6
+1/5
+1/5
+1/3
+1/2
+1/2
+1/1
+1/1
至尊木蟲 (著名寫手)
鐵桿木蟲 (著名寫手)
至尊木蟲 (著名寫手)
|
開始我標(biāo)準(zhǔn)設(shè)的沒這么高,發(fā)現(xiàn)有這個(gè)問題,于是提高了標(biāo)準(zhǔn)還是這樣。 按理說能量差值會收斂到一個(gè)比較小的變化,可這里并沒有達(dá)到那個(gè)要求,計(jì)算就終止了,并且是正常停止 curvature: -1.53 expect dE=-0.765E-07 dE for cont linesearch -0.765E-07 trial: gam=-0.00171 g(F)= 0.296E-10 g(S)= 0.499E-07 ort = 0.291E-04 (trialstep = 0.138E+01) search vector abs. value= 0.131E-09 reached required accuracy - stopping structural energy minimisation |
至尊木蟲 (著名寫手)
金蟲 (著名寫手)
|
考慮兩個(gè)氫原子,VASP通過調(diào)節(jié)他們的距離來獲得最低能量的結(jié)構(gòu),即使在最后的時(shí)候能量精度也不一定能達(dá)到1E-5,而force已經(jīng)是0了,VASP也沒有辦法調(diào)節(jié)了,所以會自己停止。 其實(shí)force的判據(jù)比total energy要高得多,而且準(zhǔn)確的多。因?yàn)閒orce是要求每個(gè)原子都達(dá)到精度要求,而total energy只是總體的判據(jù)而已。而且根據(jù)F=-\frac{\partial{E}}{\partial{x}},只要某個(gè)原子受力,沿著那個(gè)力的方向運(yùn)動,體系的能量一定會降低。所以不存在受力而是最低能量結(jié)構(gòu)的情況。 PS:最低能量結(jié)構(gòu)不是基態(tài),基態(tài)是在給定結(jié)構(gòu)時(shí)求解電子波函數(shù)的最低能量狀態(tài)才是基態(tài)。 |
至尊木蟲 (著名寫手)
至尊木蟲 (著名寫手)
金蟲 (著名寫手)
至尊木蟲 (著名寫手)
| 16 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 材料專碩306英一數(shù)二 +8 | z1z2z3879 2026-03-16 | 10/500 |
|
|---|---|---|---|---|
|
[考研] 085601材料工程專碩求調(diào)劑 +5 | 慕寒mio 2026-03-16 | 5/250 |
|
|
[考研] 341求調(diào)劑 +5 | 搗蛋豬豬 2026-03-11 | 7/350 |
|
|
[考研] 268求調(diào)劑 +8 | 一定有學(xué)上- 2026-03-14 | 9/450 |
|
|
[考研] 290求調(diào)劑 +6 | 孔志浩 2026-03-12 | 11/550 |
|
|
[考研] 275求調(diào)劑 +4 | 太陽花天天開心 2026-03-16 | 4/200 |
|
|
[論文投稿] 有沒有大佬發(fā)小論文能帶我個(gè)二作 +3 | 增銳漏人 2026-03-17 | 4/200 |
|
|
[考研] 333求調(diào)劑 +3 | 文思客 2026-03-16 | 7/350 |
|
|
[考研] 304求調(diào)劑 +5 | 素年祭語 2026-03-15 | 5/250 |
|
|
[考研] 277材料科學(xué)與工程080500求調(diào)劑 +3 | 自由煎餅果子 2026-03-16 | 3/150 |
|
|
[教師之家] 焦慮 +7 | 水冰月月野兔 2026-03-13 | 9/450 |
|
|
[考研] 070305求調(diào)劑 +3 | mlpqaz03 2026-03-14 | 4/200 |
|
|
[考研] 一志愿哈工大材料324分求調(diào)劑 +5 | 閆旭東 2026-03-14 | 5/250 |
|
|
[考研] 308 085701 四六級已過求調(diào)劑 +7 | 溫喬喬喬喬 2026-03-12 | 14/700 |
|
|
[考研] 材料與化工(0856)304求B區(qū)調(diào)劑 +6 | 邱gl 2026-03-12 | 7/350 |
|
|
[考研] 304求調(diào)劑 +6 | Mochaaaa 2026-03-12 | 7/350 |
|
|
[考研] 求調(diào)劑(材料與化工327) +4 | 愛吃香菜啦 2026-03-11 | 4/200 |
|
|
[考研] 314求調(diào)劑 +7 | 無懈可擊的巨人 2026-03-12 | 7/350 |
|
|
[考研] 277求調(diào)劑 +4 | anchor17 2026-03-12 | 4/200 |
|
|
[考研] 283求調(diào)劑,材料、化工皆可 +8 | 蘇打水7777 2026-03-11 | 10/500 |
|