版塊導(dǎo)航: 正在加載中...

應(yīng)《網(wǎng)絡(luò)安全法》要求，自2017年10月1日起，未進(jìn)行實(shí)名認(rèn)證將不得使用互聯(lián)網(wǎng)跟帖服務(wù)。為保障您的帳號(hào)能夠正常使用，請(qǐng)盡快對(duì)帳號(hào)進(jìn)行手機(jī)號(hào)驗(yàn)證，感謝您的理解與支持！

24小時(shí)熱門版塊排行榜

返回列表

cenwanglai

榮譽(yù)版主 (知名作家)

老和山猥瑣派九段

1ST強(qiáng)帖: 5
應(yīng)助: 46 (小學(xué)生)
貴賓: 8.842
金幣: 7440.4
散金: 9654
紅花: 118
沙發(fā): 23
帖子: 5306
在線: 1961.4小時(shí)
蟲號(hào): 537452
注冊(cè): 2008-04-01
性別: GG
專業(yè): 理論和計(jì)算化學(xué)
管轄: 仿真模擬

[交流] 【求助成功】為什么total drift沒有達(dá)到EDIFFG的要求,卻說reached required accuracy? 已有8人參與

EDIFFG=-0.001,呵呵,比較猛吧.
OUTCAR后面報(bào)告"reached required accuracy-stoping structural energy minimisation".

但是TOTAL-FORCE和total drift都沒有達(dá)到EDIFFG的要求.PREC=Accurate,關(guān)閉實(shí)空間后,每個(gè)原子每個(gè)方向上的力和total drift還是沒有達(dá)到要求,計(jì)算就退出來了,運(yùn)行離子步小于NSW.

這是什么原因呢?編譯錯(cuò)了嗎?

POSITION                                     TOTAL-FORCE (eV/Angst)
-----------------------------------------------------------------------------------
   2.20903    0.00000    6.13890       0.000000    0.000000    0.055890
   0.00000    0.00000    0.00000       0.000000    0.000000    0.401501
   1.10452    1.91308    3.08095       0.000000    0.000000    -0.139539
   2.20903    0.00000    1.56202       0.000000    0.000000    0.013240
   0.00000    0.00000    4.65030       0.000000    0.000000    0.044110
   1.10452    1.91308    7.62370       0.000000    0.000000    0.036490
   5.52258    -1.91308    6.13890       0.000000    0.000000    0.055890
   3.31355    -1.91308    0.00000       0.000000    0.000000    0.401501
   4.41806    0.00000    3.08095       0.000000    0.000000    -0.139539
   5.52258    -1.91308    1.56202       0.000000    0.000000    0.013240
   3.31355    -1.91308    4.65030       0.000000    0.000000    0.044110
   4.41806    0.00000    7.62370       0.000000    0.000000    0.036490
   2.20903    3.82615    6.13890       0.000000    0.000000    0.055890
   0.00000    3.82615    0.00000       0.000000    0.000000    0.401501
   1.10452    5.73923    3.08095       0.000000    0.000000    -0.139539
   2.20903    3.82615    1.56202       0.000000    0.000000    0.013240
   0.00000    3.82615    4.65030       0.000000    0.000000    0.044110
   1.10452    5.73923    7.62370       0.000000    0.000000    0.036490
   5.52258    1.91308    6.13890       0.000000    0.000000    0.055890
   3.31355    1.91308    0.00000       0.000000    0.000000    0.401501
   4.41806    3.82615    3.08095       0.000000    0.000000    -0.139539
   5.52258    1.91308    1.56202       0.000000    0.000000    0.013240
   3.31355    1.91308    4.65030       0.000000    0.000000    0.044110
   4.41806    3.82615    7.62370       0.000000    0.000000    0.036490
   0.00000    0.00000    6.92021       0.000000    0.000000    -0.138665
   2.20903    0.00000    3.86577       0.000000    0.000000    -0.147815
   1.10452    1.91308    0.78101       0.000000    0.000000    -0.125211
   3.31355    -1.91308    6.92021       0.000000    0.000000    -0.138665
   5.52258    -1.91308    3.86577       0.000000    0.000000    -0.147815
   4.41806    0.00000    0.78101       0.000000    0.000000    -0.125211
   0.00000    3.82615    6.92021       0.000000    0.000000    -0.138665
   2.20903    3.82615    3.86577       0.000000    0.000000    -0.147815
   1.10452    5.73923    0.78101       0.000000    0.000000    -0.125211
   3.31355    1.91308    6.92021       0.000000    0.000000    -0.138665
   5.52258    1.91308    3.86577       0.000000    0.000000    -0.147815
   4.41806    3.82615    0.78101       0.000000    0.000000    -0.125211
-----------------------------------------------------------------------------------
total drift:                               0.000296    -0.000609    -0.021872

INCAR:

ENCUT = 400
EDIFF = 1E-4
EDIFFG = -0.001
IBRION = 2
ISMEAR = 0
SIGMA = 0.05
NSW = 10
ALGO = Fast
#LREAR = Auto
POTIM = 0.2

呵呵,經(jīng)過一段時(shí)間觀察發(fā)現(xiàn):EDIFFG<0是對(duì)自由原子(T T T)受力來做判斷的,對(duì)固定原子(F F F)不做為判斷依據(jù).只要自由原子達(dá)到力的收斂標(biāo)準(zhǔn)就可以了.

[ Last edited by cenwanglai on 2010-5-21 at 13:37 ]

回復(fù)此樓

» 猜你喜歡

268求調(diào)劑已經(jīng)有3人回復(fù)
302求調(diào)劑已經(jīng)有5人回復(fù)
311求調(diào)劑已經(jīng)有10人回復(fù)
被我言中：新模板不強(qiáng)調(diào)格式了，假專家開始管格式了已經(jīng)有4人回復(fù)
303求調(diào)劑已經(jīng)有5人回復(fù)
考研求調(diào)劑已經(jīng)有4人回復(fù)
085601專碩，總分342求調(diào)劑，地區(qū)不限已經(jīng)有4人回復(fù)
085601材料工程專碩求調(diào)劑已經(jīng)有5人回復(fù)
一志愿天津大學(xué)化學(xué)工藝專業(yè)（081702）315分求調(diào)劑已經(jīng)有6人回復(fù)
能源材料化學(xué)課題組招收碩士研究生8-10名已經(jīng)有6人回復(fù)

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

為什么 Johnson wiley 中chirality的數(shù)據(jù)庫有部分文章是可以免費(fèi)下的呢？已經(jīng)有4人回復(fù)
求助：PLGA合成中要求壓力較低，1000pa，怎么控制壓力不變呢？已經(jīng)有10人回復(fù)
紫外分光光度計(jì)和原子吸收分光光度計(jì)單色器位置為什么不同？已經(jīng)有7人回復(fù)
為什么液相色譜壓力一直不穩(wěn)？已經(jīng)有6人回復(fù)
Required Reviews Completed 變?yōu)閑nder review 又變?yōu)镽equired Reviews Completed 已經(jīng)有12人回復(fù)
做上環(huán)手術(shù)弄得大出血，醫(yī)院卻說正常已經(jīng)有4人回復(fù)
★★★教你圖片如何處理才能達(dá)到SCI投稿要求已經(jīng)有61人回復(fù)
請(qǐng)問這個(gè)數(shù)據(jù)是否符合要求已經(jīng)有5人回復(fù)
required reviews completed已經(jīng)一個(gè)月已經(jīng)有19人回復(fù)
等待Required Reviews Completed 已經(jīng)有27人回復(fù)
SBB 一個(gè)月Required Reviews Completed 已經(jīng)有12人回復(fù)

1樓 2010-05-07 09:18:42

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

回帖支持 ( 顯示支持度最高的前 50 名 )

fanchen021

鐵桿木蟲 (正式寫手)

應(yīng)助: 7 (幼兒園)
金幣: 7994.9
散金: 1844
紅花: 4
帖子: 776
在線: 569小時(shí)
蟲號(hào): 773188
注冊(cè): 2009-05-18
專業(yè): 理論和計(jì)算化學(xué)

★
cenwanglai(金幣+3): 2010-05-10 18:11:52
qasd(金幣+1):謝謝 2010-05-10 22:58:57

那你試試 grep 'FORCES: max atom, RMS' OUTCAR

贊一下(2人)

回復(fù)此樓

11樓2010-05-10 13:54:04

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

普通回帖

jinhuan123

金蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 529.3
散金: 21
紅花: 1
帖子: 106
在線: 51.7小時(shí)
蟲號(hào): 805643
注冊(cè): 2009-07-08
專業(yè): 催化化學(xué)

★
cenwanglai(金幣+1):你用什么?呵呵 2010-05-07 10:05:10
gavinliu7390(金幣+1):謝謝交流！ 2010-05-10 08:47:44

這個(gè)貌似沒有用到過哦

贊一下(1人)

回復(fù)此樓

2樓2010-05-07 09:50:38

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

jinhuan123

金蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 529.3
散金: 21
紅花: 1
帖子: 106
在線: 51.7小時(shí)
蟲號(hào): 805643
注冊(cè): 2009-07-08
專業(yè): 催化化學(xué)

★
gavinliu7390(金幣+1):感謝交流！ 2010-05-07 10:46:46

OUTCAR 里確實(shí)有EDIFFG，但是具體怎么體現(xiàn)出來，就像你列出的這個(gè)，現(xiàn)在還沒有研究過呢，期待......

贊一下(1人)

回復(fù)此樓

3樓2010-05-07 10:10:01

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

fanchen021

鐵桿木蟲 (正式寫手)

應(yīng)助: 7 (幼兒園)
金幣: 7994.9
散金: 1844
紅花: 4
帖子: 776
在線: 569小時(shí)
蟲號(hào): 773188
注冊(cè): 2009-05-18
專業(yè): 理論和計(jì)算化學(xué)

★
cenwanglai(金幣+2):呵呵,怎么運(yùn)行呢?我在運(yùn)行結(jié)束后含有OUTCAR的目錄下運(yùn)行$vef.pl,生成的fe.dat是空的呢. 2010-05-07 20:56:45
gavinliu7390(金幣+1):謝謝交流！ 2010-05-10 08:48:24

用vef.pl 看一下啊

贊一下(1人)

回復(fù)此樓

4樓2010-05-07 10:11:28

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

citrine

鐵桿木蟲 (正式寫手)

應(yīng)助: 0 (幼兒園)
金幣: 8823.3
散金: 1
紅花: 7
帖子: 460
在線: 299.9小時(shí)
蟲號(hào): 427638
注冊(cè): 2007-07-29
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

★ ★
cenwanglai(金幣+2):謝謝參與.按你說的意思,是要判斷每一個(gè)原子的受力的哈.但是不知道為什么我的任務(wù)會(huì)在力沒有達(dá)到收斂精度的情況下報(bào)告"達(dá)到收斂精度,退停止結(jié)構(gòu)能量最小化" 2010-05-07 20:59:32
gavinliu7390(金幣+2):謝謝交流！ 2010-05-10 08:48:43

POSITION TOTAL-FORCE (eV/Angst)
-----------------------------------------------------------------------------------
2.20903 0.00000 6.13890 0.000000 0.000000 0.055890

最后三個(gè)數(shù)是原子收到的力，顯然沒有達(dá)到收斂標(biāo)準(zhǔn)。
我注意到NSW=10，這個(gè)參數(shù)應(yīng)該變大，比如說100

贊一下(1人)

回復(fù)此樓

5樓2010-05-07 11:27:40

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

stractor

金蟲 (著名寫手)

1ST強(qiáng)帖: 6
應(yīng)助: 91 (初中生)
金幣: 1487.6
散金: 666
紅花: 30
沙發(fā): 1
帖子: 2224
在線: 426.8小時(shí)
蟲號(hào): 394087
注冊(cè): 2007-06-06
專業(yè): 半導(dǎo)體材料

★
cenwanglai(金幣+2):呵呵,即使設(shè)置成EDIFFG=-0.001,即使實(shí)際上力沒有達(dá)到要求的收斂精度,但是還是會(huì)以"reached required accuracy"結(jié)束.不知道是不是編譯的原因. 2010-05-07 20:53:45
gavinliu7390(金幣+1):謝謝交流！ 2010-05-10 08:48:58

EDIFFG=-0.001,收斂標(biāo)準(zhǔn)確實(shí)挺高，建議改為-0.01至-0.03.

贊一下(1人)

回復(fù)此樓

6樓2010-05-07 12:10:18

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

zengxiangm

木蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 1981.2
帖子: 221
在線: 96.9小時(shí)
蟲號(hào): 676900
注冊(cè): 2008-12-18
性別: GG
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

老師說是每一個(gè)原子每一個(gè)方向上的受力都要小于|EDIFFG|，但我的情況也和樓主一樣，希望有高人來解答

贊一下

回復(fù)此樓

7樓2010-05-07 12:18:40

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

valenhou001

至尊木蟲 (職業(yè)作家)

1ST強(qiáng)帖: 13
應(yīng)助: 241 (大學(xué)生)
金幣: 25701.7
散金: 602
紅花: 166
帖子: 3782
在線: 873.8小時(shí)
蟲號(hào): 1007127
注冊(cè): 2010-04-27
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

★
cenwanglai(金幣+2):謝謝參與.grep "F" OSZICAR|tail -1給出的F就是總能? 2010-05-07 20:55:08
gavinliu7390(金幣+1):感謝交流！ 2010-05-07 23:27:39

EDIFF為負(fù)數(shù)時(shí)，以原子所受的力來判斷收斂；如果為正數(shù)，以前后兩次總能的變化為依據(jù)來判斷收斂。

贊一下(1人)

回復(fù)此樓

8樓2010-05-07 14:11:45

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

stractor

金蟲 (著名寫手)

1ST強(qiáng)帖: 6
應(yīng)助: 91 (初中生)
金幣: 1487.6
散金: 666
紅花: 30
沙發(fā): 1
帖子: 2224
在線: 426.8小時(shí)
蟲號(hào): 394087
注冊(cè): 2007-06-06
專業(yè): 半導(dǎo)體材料

★ ★ ★ ★
cenwanglai(金幣+1):謝謝參與
cenwanglai(金幣+1): 設(shè)置為-0.001純粹是跑一下看看會(huì)出現(xiàn)什么情況.呵呵 2010-05-08 13:43:52
gavinliu7390(金幣+2):謝謝交流！ 2010-05-10 08:49:15
zzy870720z(金幣+1):感謝交流 2010-05-10 08:49:44

此外，比如你優(yōu)化的步數(shù)如果設(shè)置為20步，NSW=20，到第20步結(jié)束時(shí)，哪怕力沒有達(dá)到標(biāo)準(zhǔn)，VASP也會(huì)reached required accuracy而結(jié)束。
你為什么要設(shè)置-0.001？這個(gè)值太小了，費(fèi)時(shí)。-0.01至-0.03夠了。

贊一下(3人)

回復(fù)此樓

9樓2010-05-08 13:27:57

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

cenwanglai

榮譽(yù)版主 (知名作家)

老和山猥瑣派九段

1ST強(qiáng)帖: 5
應(yīng)助: 46 (小學(xué)生)
貴賓: 8.842
金幣: 7440.4
散金: 9654
紅花: 118
沙發(fā): 23
帖子: 5306
在線: 1961.4小時(shí)
蟲號(hào): 537452
注冊(cè): 2008-04-01
性別: GG
專業(yè): 理論和計(jì)算化學(xué)
管轄: 仿真模擬

這個(gè)體系收斂沒有問題了,問題是收斂的結(jié)果不對(duì).
vef.pl 怎么運(yùn)行呢?運(yùn)行vasp后,直接在工作目錄下的命令提示符下輸入?
$vef.pl
有文件fe.dat生成,但是是空的文件.

贊一下

回復(fù)此樓

10樓2010-05-10 08:42:21

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

相關(guān)版塊跳轉(zhuǎn) 我要訂閱樓主 cenwanglai 的主題更新

返回列表

普通表情龍兔虎貓高級(jí)回復(fù) (可上傳附件)

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

[考研] 268求調(diào)劑 +7	好運(yùn)連綿不絕 2026-03-12	8/400	2026-03-17 20:28 by xilongliang
[考研] 材料與化工求調(diào)劑 +6	為學(xué)666 2026-03-16	6/300	2026-03-17 20:15 by peike
[考研] 材料工程專碩調(diào)劑 +5	204818@lcx 2026-03-17	5/250	2026-03-17 17:27 by Little-xue
[考研] 302求調(diào)劑 +9	負(fù)心者當(dāng)誅 2026-03-11	9/450	2026-03-17 17:13 by ruiyingmiao
[考研] 274求調(diào)劑 +5	時(shí)間點(diǎn) 2026-03-13	5/250	2026-03-17 07:34 by 熱情沙漠
[考研] 333求調(diào)劑 +3	文思客 2026-03-16	7/350	2026-03-16 18:21 by 文思客
[考研] 318求調(diào)劑 +3	Yanyali 2026-03-15	3/150	2026-03-16 16:41 by houyaoxu
[考研] 283求調(diào)劑 +10	小樓。 2026-03-12	14/700	2026-03-16 16:08 by 13811244083
[考研] 0703 物理化學(xué)調(diào)劑 +3	我可以上岸的對(duì)?/a> 2026-03-13	5/250	2026-03-16 10:50 by 我可以上岸的對(duì)?/a>
[考研] 0856求調(diào)劑 +3	劉夢(mèng)微 2026-03-15	3/150	2026-03-16 10:00 by houyaoxu
[考研] 中科大材料專碩319求調(diào)劑 +3	孟鑫材料 2026-03-13	3/150	2026-03-14 18:10 by houyaoxu
[考研] 297一志愿上交085600求調(diào)劑 +5	指尖八千里 2026-03-14	5/250	2026-03-14 17:26 by a不易
[考研] 復(fù)試調(diào)劑 +3	呼呼？~+123456 2026-03-14	3/150	2026-03-14 16:53 by WTUChen
[考研] 328求調(diào)劑 +3	5201314Lsy！ 2026-03-13	6/300	2026-03-14 15:31 by hyswxzs
[考研] 26調(diào)劑/材料/英一數(shù)二/總分289/已過A區(qū)線 +6	步川酷紫123 2026-03-13	6/300	2026-03-13 21:59 by 星空星月
[考研] 310求調(diào)劑 +3	【上上簽】 2026-03-11	3/150	2026-03-13 16:16 by JourneyLucky
[考研] 307求調(diào)劑 +5	超級(jí)伊昂大王 2026-03-12	5/250	2026-03-13 15:56 by 棒棒球手
[考研] 328化工專碩求調(diào)劑 +4	。，。，。，。i 2026-03-12	4/200	2026-03-13 14:44 by JourneyLucky
[考研] 085600材料與化工 309分請(qǐng)求調(diào)劑 +7	dtdxzxx 2026-03-12	8/400	2026-03-13 14:43 by jxchenghu
[考研] 0856化學(xué)工程280分求調(diào)劑 +4	shenzxsn 2026-03-11	4/200	2026-03-13 11:55 by ymwdoctor