| 16 | 1/1 | 返回列表 |
| 查看: 3931 | 回復(fù): 15 | ||||||
sunyang1988至尊木蟲 (著名寫手)
|
[交流]
【求助】VASP收斂問題
|
|||||
|
計算中選用的收斂標(biāo)準(zhǔn)是能量 EDIFF = 1E-5 EDIFFG = 0.0001 但運(yùn)算正常結(jié)束后,能量差值并沒有達(dá)到這個要求 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 不知道這種情況到底算不算收斂? |
材料計算模擬實(shí)用技巧 | VASP | vasp計算資料 |
» 搶金幣啦!回帖就可以得到:
+1/84
+1/84
+1/79
+1/71
+1/67
+1/50
+2/38
+1/35
+1/32
+1/18
+1/17
+1/14
+1/13
+1/12
+1/7
+1/6
+1/6
+1/6
+1/2
+1/2
至尊木蟲 (著名寫手)
至尊木蟲 (著名寫手)
|
開始我標(biāo)準(zhǔn)設(shè)的沒這么高,發(fā)現(xiàn)有這個問題,于是提高了標(biāo)準(zhǔn)還是這樣。 按理說能量差值會收斂到一個比較小的變化,可這里并沒有達(dá)到那個要求,計算就終止了,并且是正常停止 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 |
至尊木蟲 (著名寫手)
金蟲 (著名寫手)
|
考慮兩個氫原子,VASP通過調(diào)節(jié)他們的距離來獲得最低能量的結(jié)構(gòu),即使在最后的時候能量精度也不一定能達(dá)到1E-5,而force已經(jīng)是0了,VASP也沒有辦法調(diào)節(jié)了,所以會自己停止。 其實(shí)force的判據(jù)比total energy要高得多,而且準(zhǔn)確的多。因?yàn)閒orce是要求每個原子都達(dá)到精度要求,而total energy只是總體的判據(jù)而已。而且根據(jù)F=-\frac{\partial{E}}{\partial{x}},只要某個原子受力,沿著那個力的方向運(yùn)動,體系的能量一定會降低。所以不存在受力而是最低能量結(jié)構(gòu)的情況。 PS:最低能量結(jié)構(gòu)不是基態(tài),基態(tài)是在給定結(jié)構(gòu)時求解電子波函數(shù)的最低能量狀態(tài)才是基態(tài)。 |
至尊木蟲 (著名寫手)
至尊木蟲 (著名寫手)
金蟲 (著名寫手)
至尊木蟲 (著名寫手)
| 16 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 0703化學(xué)336分求調(diào)劑 +5 | zbzihdhd 2026-03-15 | 6/300 |
|
|---|---|---|---|---|
|
[考研] 304求調(diào)劑 +9 | 小熊joy 2026-03-14 | 9/450 |
|
|
[考研] 301求調(diào)劑 +6 | yy要上岸呀 2026-03-17 | 6/300 |
|
|
[考研] 268求調(diào)劑 +7 | 好運(yùn)連綿不絕 2026-03-12 | 8/400 |
|
|
[考研] 268求調(diào)劑 +8 | 一定有學(xué)上- 2026-03-14 | 9/450 |
|
|
[考研] 290求調(diào)劑 +3 | p asserby. 2026-03-15 | 4/200 |
|
|
[考研] 材料與化工專碩調(diào)劑 +5 | heming3743 2026-03-16 | 5/250 |
|
|
[碩博家園] 深圳大學(xué)碩士招生(2026秋,傳感器方向,僅錄取第一志愿) +4 | xujiaoszu 2026-03-11 | 9/450 |
|
|
[考研] 化學(xué)調(diào)劑0703 +8 | 啊我我的 2026-03-11 | 8/400 |
|
|
[考研] 085600調(diào)劑 +5 | 漾漾123sun 2026-03-12 | 6/300 |
|
|
[考研] 085600材料與化工 求調(diào)劑 +13 | enenenhui 2026-03-13 | 14/700 |
|
|
[考研] 0703化學(xué)調(diào)劑 290分有科研經(jīng)歷,論文在投 +7 | 膩膩gk 2026-03-14 | 7/350 |
|
|
[考研] 26考研一志愿中國石油大學(xué)(華東)305分求調(diào)劑 +3 | 嘉年新程 2026-03-15 | 3/150 |
|
|
[考研] 297一志愿上交085600求調(diào)劑 +5 | 指尖八千里 2026-03-14 | 5/250 |
|
|
[考研] 328求調(diào)劑 +3 | 5201314Lsy! 2026-03-13 | 6/300 |
|
|
[考研] 【0703化學(xué)調(diào)劑】-一志愿華中師范大學(xué)-六級475 +5 | Becho359 2026-03-11 | 5/250 |
|
|
[考研] 材料工程調(diào)劑 +9 | 咪咪空空 2026-03-12 | 9/450 |
|
|
[考研] [0860]321分求調(diào)劑,ab區(qū)皆可 +4 | 寶貴熱 2026-03-13 | 4/200 |
|
|
[考研] 308求調(diào)劑 +3 | 是Lupa啊 2026-03-12 | 3/150 |
|
|
[考研] 283求調(diào)劑,材料、化工皆可 +8 | 蘇打水7777 2026-03-11 | 10/500 |
|