當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖

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

不知道這種情況到底算不算收斂？

回復(fù)此樓

» 搶金幣啦！回帖就可以得到:

查看全部散金貼

1樓 2011-03-15 09:32:02

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

sunyang1988

至尊木蟲 (著名寫手)

引用回帖:

Originally posted by bingmou at 2011-03-15 15:43:21:
考慮兩個(gè)氫原子，VASP通過調(diào)節(jié)他們的距離來獲得最低能量的結(jié)構(gòu)，即使在最后的時(shí)候能量精度也不一定能達(dá)到1E-5,而force已經(jīng)是0了，VASP也沒有辦法調(diào)節(jié)了，所以會(huì)自己停止。
其實(shí)force的判據(jù)比total energy要高得多 ...

能量很難收斂到1E-5，不過這里也收斂的太差了
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ì)稱的體系會(huì)出現(xiàn)這樣的問題，這里我已經(jīng)關(guān)掉了對(duì)稱性，ISYM=0，問題依舊

贊一下

回復(fù)此樓

9樓2011-03-15 16:19:59

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

查看全部 16 個(gè)回答

賀儀

鐵桿木蟲 (著名寫手)

★
sunyang1988(金幣+2): 2011-03-15 10:28:31
gzqdyouxia(金幣+1): 謝謝交流 2011-03-15 11:15:27

引用回帖:

Originally posted by sunyang1988 at 2011-03-15 09:32:02:
計(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
...

這些原子的受力基本上是每一個(gè)都有，是不是會(huì)存在某些原子受力大但是能量更低的點(diǎn)呢？如果有的話，我個(gè)人認(rèn)為還是以能量差值很小而且基本是最低的為好吧。樓主的EDIFFG = 0.0001要求太高了吧？我計(jì)算一般很難達(dá)到這個(gè)要求，樓主你資源豐富呀！羨慕……

贊一下(1人)

回復(fù)此樓

2樓2011-03-15 09:54:10

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

sunyang1988

至尊木蟲 (著名寫手)

引用回帖:

Originally posted by 賀儀 at 2011-03-15 09:54:10:
這些原子的受力基本上是每一個(gè)都有，是不是會(huì)存在某些原子受力大但是能量更低的點(diǎn)呢？如果有的話，我個(gè)人認(rèn)為還是以能量差值很小而且基本是最低的為好吧。樓主的EDIFFG = 0.0001要求太高了吧？我計(jì)算一般很難達(dá) ...

開始我標(biāo)準(zhǔn)設(shè)的沒這么高，發(fā)現(xiàn)有這個(gè)問題，于是提高了標(biāo)準(zhǔn)還是這樣。
按理說能量差值會(huì)收斂到一個(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

贊一下

回復(fù)此樓

3樓2011-03-15 10:28:19

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

leila3409

捐助貴賓 (正式寫手)

應(yīng)助: 0 (幼兒園)
金幣: 616.2
帖子: 310
在線: 71.5小時(shí)
蟲號(hào): 985885

★ ★
gzqdyouxia(金幣+1): 謝謝交流 2011-03-15 11:15:35
zzy870720z(金幣+1): 謝謝提示 2011-03-15 12:01:49
sunyang1988(金幣+2): 2011-03-15 15:16:55

可以認(rèn)為是收斂了。原子受力小于小數(shù)點(diǎn)以后第三位一般就可以停了。

贊一下(2人)

回復(fù)此樓

4樓2011-03-15 10:47:00

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

最具人氣熱帖推薦 [查看全部]	作者	回/看	最后發(fā)表

[考研] 能源材料化學(xué)課題組招收碩士研究生8-10名 +3	脫穎而出 2026-03-16	6/300	2026-03-17 21:19 by z1z2z3879
[考研] 材料與化工求調(diào)劑 +6	為學(xué)666 2026-03-16	6/300	2026-03-17 20:15 by peike
[考研] 328求調(diào)劑，英語(yǔ)六級(jí)551，有科研經(jīng)歷 +3	生物工程調(diào)劑 2026-03-16	8/400	2026-03-17 19:03 by Wangjingyue
[考研] 070300化學(xué)319求調(diào)劑 +4	錦鯉0909 2026-03-17	4/200	2026-03-17 18:21 by 重科小霸王
[考研] 085600材料與化工 +4	安全上岸！ 2026-03-16	4/200	2026-03-17 14:02 by 勇敢太監(jiān)王公公
[論文投稿] 有沒有大佬發(fā)小論文能帶我個(gè)二作 +3	增銳漏人 2026-03-17	4/200	2026-03-17 09:26 by xs74101122
[考研] 一志愿，福州大學(xué)材料專碩339分求調(diào)劑 +3	木子momo青爭(zhēng) 2026-03-15	3/150	2026-03-17 07:52 by laoshidan
[考研] 318求調(diào)劑 +3	Yanyali 2026-03-15	3/150	2026-03-16 16:41 by houyaoxu
[考研] 311求調(diào)劑 +5	26研0 2026-03-15	5/250	2026-03-16 16:21 by a不易
[考研] 283求調(diào)劑 +10	小樓。 2026-03-12	14/700	2026-03-16 16:08 by 13811244083
[考研] 326求調(diào)劑 +3	mlpqaz03 2026-03-15	3/150	2026-03-16 07:33 by Iveryant
[考博] 東華理工大學(xué)化材專業(yè)26屆碩士博士申請(qǐng) +6	zlingli 2026-03-13	6/300	2026-03-15 20:00 by ryzcf
[考研] 297一志愿上交085600求調(diào)劑 +5	指尖八千里 2026-03-14	5/250	2026-03-14 17:26 by a不易
[考研] 331求調(diào)劑（0703有機(jī)化學(xué) +5	ZY-05 2026-03-13	6/300	2026-03-14 10:51 by Jy?
[考研] 招收0805（材料）調(diào)劑 +3	18595523086 2026-03-13	3/150	2026-03-14 00:33 by 123%、
[考研] 0703，333分求調(diào)劑一志愿鄭州大學(xué)-物理化學(xué) +3	李魔女斗篷 2026-03-11	3/150	2026-03-13 22:24 by JourneyLucky
[考研] 304求調(diào)劑 +6	Mochaaaa 2026-03-12	7/350	2026-03-13 22:18 by 星空星月
[考研] 材料與化工求調(diào)劑一志愿 985 總分 295 +8	dream…… 2026-03-12	8/400	2026-03-13 22:17 by 星空星月
[考研] 329求調(diào)劑 +3	miaodesi 2026-03-12	4/200	2026-03-13 20:53 by 18595523086
[考研] 求調(diào)劑 +3	程雨杭 2026-03-12	3/150	2026-03-13 15:06 by JourneyLucky