| 24小時(shí)熱門(mén)版塊排行榜 |
| 16 | 1/1 | 返回列表 |
| 查看: 3956 | 回復(fù): 15 | |||||
sunyang1988至尊木蟲(chóng) (著名寫(xiě)手)
|
[交流]
【求助】VASP收斂問(wèn)題
|
||||
|
計(jì)算中選用的收斂標(biāo)準(zhǔn)是能量 EDIFF = 1E-5 EDIFFG = 0.0001 但運(yùn)算正常結(jié)束后,能量差值并沒(méi)有達(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/183
+1/86
+1/81
+2/46
+1/37
+1/37
+1/36
+1/31
+3/25
+1/25
+1/10
+1/7
+1/7
+1/7
+1/6
+1/4
+1/3
+1/2
+1/2
+3/1
至尊木蟲(chóng) (著名寫(xiě)手)
鐵桿木蟲(chóng) (著名寫(xiě)手)
至尊木蟲(chóng) (著名寫(xiě)手)
|
開(kāi)始我標(biāo)準(zhǔn)設(shè)的沒(méi)這么高,發(fā)現(xiàn)有這個(gè)問(wèn)題,于是提高了標(biāo)準(zhǔn)還是這樣。 按理說(shuō)能量差值會(huì)收斂到一個(gè)比較小的變化,可這里并沒(méi)有達(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 |
捐助貴賓 (正式寫(xiě)手)
至尊木蟲(chóng) (著名寫(xiě)手)
金蟲(chóng) (著名寫(xiě)手)
|
考慮兩個(gè)氫原子,VASP通過(guò)調(diào)節(jié)他們的距離來(lái)獲得最低能量的結(jié)構(gòu),即使在最后的時(shí)候能量精度也不一定能達(dá)到1E-5,而force已經(jīng)是0了,VASP也沒(méi)有辦法調(diào)節(jié)了,所以會(huì)自己停止。 其實(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)動(dòng),體系的能量一定會(huì)降低。所以不存在受力而是最低能量結(jié)構(gòu)的情況。 PS:最低能量結(jié)構(gòu)不是基態(tài),基態(tài)是在給定結(jié)構(gòu)時(shí)求解電子波函數(shù)的最低能量狀態(tài)才是基態(tài)。 |
捐助貴賓 (正式寫(xiě)手)
至尊木蟲(chóng) (著名寫(xiě)手)
至尊木蟲(chóng) (著名寫(xiě)手)
|
能量很難收斂到1E-5,不過(guò)這里也收斂的太差了 free energy TOTEN = -462.008990 eV free energy TOTEN = -462.024537 eV free energy TOTEN = -462.031260 eV free energy TOTEN = -462.037495 eV 之收斂到1E-2 之前算的一些體系一般力收斂到1E-2能量都收斂到接近1E-4 總結(jié)了一下,高對(duì)稱(chēng)的體系會(huì)出現(xiàn)這樣的問(wèn)題,這里我已經(jīng)關(guān)掉了對(duì)稱(chēng)性,ISYM=0,問(wèn)題依舊 |
金蟲(chóng) (著名寫(xiě)手)
銀蟲(chóng) (小有名氣)
木蟲(chóng) (正式寫(xiě)手)
至尊木蟲(chóng) (著名寫(xiě)手)
| 16 | 1/1 | 返回列表 |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
|
百度網(wǎng)盤(pán) |
360云盤(pán) |
千易網(wǎng)盤(pán) |
華為網(wǎng)盤(pán)
在新窗口頁(yè)面中打開(kāi)自己喜歡的網(wǎng)盤(pán)網(wǎng)站,將文件上傳后,然后將下載鏈接復(fù)制到帖子內(nèi)容中就可以了。 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 085601 材料工程 313分 求調(diào)劑 +6 | Ong3 2026-03-27 | 6/300 |
|
|---|---|---|---|---|
|
[考研] 化學(xué)0703 調(diào)劑 306分 一志愿211 +9 | 26要上岸 2026-03-28 | 9/450 |
|
|
[考研] 一志愿中海洋材料357 +4 | 麥恩莉. 2026-03-30 | 4/200 |
|
|
[考研] 求調(diào)劑 生物學(xué) 377分 +5 | zzll03 2026-03-31 | 5/250 |
|
|
[考研] 英一數(shù)一總分334求調(diào)劑 +4 | 陳陽(yáng)坤 2026-03-31 | 4/200 |
|
|
[電化學(xué)] 070300化學(xué)調(diào)劑 +3 | 山頂見(jiàn)α 2026-03-25 | 3/150 |
|
|
[考研] 一志愿食品科學(xué)與工程083200求調(diào)劑 +4 | XQTJZ 2026-03-30 | 4/200 |
|
|
[考研] 083000環(huán)境科學(xué)與工程調(diào)劑,總分281 +4 | 橙子(勝意) 2026-03-30 | 4/200 |
|
|
[考研] 322求調(diào)劑 +10 | 宋明欣 2026-03-27 | 10/500 |
|
|
[考研] 一志愿廈門(mén)大學(xué)材料工程專(zhuān)碩354找調(diào)劑。! +5 | 貝唄鋇鋇 2026-03-30 | 5/250 |
|
|
[考研] 332求調(diào)劑 +6 | @MZB382400 2026-03-28 | 6/300 |
|
|
[碩博家園] 求調(diào)劑 有機(jī)化學(xué)考研356分 +10 | Nadiums 2026-03-25 | 11/550 |
|
|
[考研] 085602 化學(xué)工程專(zhuān)碩 340分求調(diào)劑 +4 | qianbai11 2026-03-29 | 4/200 |
|
|
[考研] 086000生物與醫(yī)藥調(diào)劑 +5 | Feisty。 2026-03-28 | 9/450 |
|
|
[考研] 085600,專(zhuān)業(yè)課化工原理,321分求調(diào)劑 +5 | 大饞小子 2026-03-28 | 5/250 |
|
|
[考研] 298求調(diào)劑 +4 | 種圣賜 2026-03-28 | 4/200 |
|
|
[考研] 299求調(diào)劑 +7 | 嗯嗯嗯嗯2 2026-03-27 | 7/350 |
|
|
[考研] 283求調(diào)劑 +7 | A child 2026-03-28 | 7/350 |
|
|
[考研] 265求調(diào)劑11408 +3 | 劉小鹿lu 2026-03-27 | 3/150 |
|
|
[考研] 296求調(diào)劑 +4 | 汪。! 2026-03-25 | 7/350 |
|