當(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ù)此樓

» 收錄本帖的淘帖專輯推薦

材料計(jì)算模擬實(shí)用技巧

VASP

vasp計(jì)算資料

» 猜你喜歡

深圳大學(xué)碩士招生（2026秋，傳感器方向，僅錄取第一志愿）已經(jīng)有10人回復(fù)
278求調(diào)劑已經(jīng)有5人回復(fù)
化學(xué)工程321分求調(diào)劑已經(jīng)有15人回復(fù)
材料專碩306英一數(shù)二已經(jīng)有11人回復(fù)
334求調(diào)劑已經(jīng)有3人回復(fù)
268求調(diào)劑已經(jīng)有5人回復(fù)
293求調(diào)劑已經(jīng)有16人回復(fù)
302求調(diào)劑已經(jīng)有7人回復(fù)
312求調(diào)劑已經(jīng)有6人回復(fù)
281求調(diào)劑（0805）已經(jīng)有9人回復(fù)

» 本主題相關(guān)價(jià)值貼推薦，對您同樣有幫助:

vasp計(jì)算過渡態(tài)問題已經(jīng)有13人回復(fù)
VASP結(jié)構(gòu)優(yōu)化時(shí)出現(xiàn)的POTIM問題已經(jīng)有4人回復(fù)
關(guān)于vasp收斂和優(yōu)化已經(jīng)有24人回復(fù)
vasp5.2計(jì)算誤差太大已經(jīng)有11人回復(fù)
求助：vasp結(jié)構(gòu)優(yōu)化不能收斂已經(jīng)有18人回復(fù)
vasp計(jì)算：優(yōu)化能收斂，但同樣的設(shè)置靜態(tài)計(jì)算不收斂已經(jīng)有3人回復(fù)
VASP馳豫自動(dòng)停止問題。。已經(jīng)有9人回復(fù)
vasp優(yōu)化結(jié)構(gòu)中出現(xiàn)的問題已經(jīng)有10人回復(fù)
VASP計(jì)算時(shí)的幾個(gè)細(xì)節(jié)設(shè)置問題，希望有經(jīng)驗(yàn)的人過來指點(diǎn)一下！已經(jīng)有23人回復(fù)
【求助成功】vasp測試K點(diǎn)收斂性的腳本問題已經(jīng)有20人回復(fù)
【求助】vasp靜態(tài)計(jì)算能量值錯(cuò)誤已經(jīng)有21人回復(fù)
【求助】vasp不同磁結(jié)構(gòu)最后的收斂問題已經(jīng)有6人回復(fù)

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

查看全部散金貼

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

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

sunyang1988

至尊木蟲 (著名寫手)

引用回帖:

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

開始我標(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

贊一下

回復(fù)此樓

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

已閱回復(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è)都有，是不是會存在某些原子受力大但是能量更低的點(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的回帖

leila3409

捐助貴賓 (正式寫手)

應(yīng)助: 0 (幼兒園)
金幣: 616.2
帖子: 310
在線: 71.5小時(shí)
蟲號: 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的回帖

sunyang1988

至尊木蟲 (著名寫手)

引用回帖:

Originally posted by leila3409 at 2011-03-15 10:47:00:
可以認(rèn)為是收斂了。原子受力小于小數(shù)點(diǎn)以后第三位一般就可以停了。

可是還沒有到收斂到0.001eV啊

贊一下

回復(fù)此樓

5樓2011-03-15 15:17:19

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

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

[考研] 考研求調(diào)劑 +3	橘頌. 2026-03-17	4/200	2026-03-17 21:43 by 有只貍奴
[考研] 296求調(diào)劑 +5	大口吃飯身體健 2026-03-13	5/250	2026-03-17 21:05 by 不惑可樂
[考研] 一志愿蘇州大學(xué)材料工程（085601）專碩有科研經(jīng)歷三項(xiàng)國獎(jiǎng)兩個(gè)實(shí)用型專利一項(xiàng)省級立項(xiàng) +6	大火山小火山 2026-03-16	8/400	2026-03-17 15:05 by 無懈可擊111
[考研] 211本，11408一志愿中科院277分，曾在中科院自動(dòng)化所實(shí)習(xí) +6	Losir 2026-03-12	7/350	2026-03-17 12:09 by danranxie
[考研] 275求調(diào)劑 +4	太陽花天天開心 2026-03-16	4/200	2026-03-17 10:53 by 功夫瘋狂
[考研] 267一志愿南京工業(yè)大學(xué)0817化工求調(diào)劑 +6	SUICHILD 2026-03-12	6/300	2026-03-17 09:24 by 霧散后相遇lc
[考研] 東南大學(xué)364求調(diào)劑 +5	JasonYuiui 2026-03-15	5/250	2026-03-16 21:28 by 木瓜膏
[考研] 321求調(diào)劑 +5	大米飯！ 2026-03-15	5/250	2026-03-16 16:33 by houyaoxu
[考研] 277材料科學(xué)與工程080500求調(diào)劑 +3	自由煎餅果子 2026-03-16	3/150	2026-03-16 14:10 by 運(yùn)氣yunqi
[考研] 326求調(diào)劑 +3	mlpqaz03 2026-03-15	3/150	2026-03-16 07:33 by Iveryant
[考研] 材料與化工 323 英一+數(shù)二+物化，一志愿：哈工大本人本科雙一流 +4	自由的_飛翔 2026-03-13	5/250	2026-03-14 19:39 by hmn_wj
[考研] 329求調(diào)劑 +3	miaodesi 2026-03-12	4/200	2026-03-13 20:53 by 18595523086
[考研] 求調(diào)劑 +5	一定有學(xué)上- 2026-03-12	5/250	2026-03-13 18:31 by ms629
[考研] 314求調(diào)劑 +7	無懈可擊的巨人 2026-03-12	7/350	2026-03-13 15:40 by JourneyLucky
[考研] 土木第一志愿276求調(diào)劑，科研和技能十分豐富，求新興方向的導(dǎo)師收留 +3	土木小天才 2026-03-12	3/150	2026-03-13 15:01 by JourneyLucky
[考研] 0856化學(xué)工程280分求調(diào)劑 +4	shenzxsn 2026-03-11	4/200	2026-03-13 11:55 by ymwdoctor
[考研] 289求調(diào)劑 +3	李政瑩 2026-03-12	3/150	2026-03-13 11:02 by 求調(diào)劑zz
[考研] 290求調(diào)劑 +3	ADT 2026-03-13	3/150	2026-03-13 10:19 by peike
[考研] 283求調(diào)劑，材料、化工皆可 +8	蘇打水7777 2026-03-11	10/500	2026-03-13 09:06 by Linda Hu
[考研] 333求調(diào)劑 +3	152697 2026-03-12	4/200	2026-03-13 07:08 by Iveryant