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

登錄注冊

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

24小時熱門版塊排行榜

返回列表

cenwanglai

榮譽版主 (知名作家)

老和山猥瑣派九段

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

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

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

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

這是什么原因呢?編譯錯了嗎?

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)過一段時間觀察發(fā)現(xiàn):EDIFFG<0是對自由原子(T T T)受力來做判斷的,對固定原子(F F F)不做為判斷依據(jù).只要自由原子達到力的收斂標準就可以了.

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

回復(fù)此樓

» 猜你喜歡

293求調(diào)劑已經(jīng)有14人回復(fù)
0703化學(xué)336分求調(diào)劑已經(jīng)有6人回復(fù)
304求調(diào)劑已經(jīng)有9人回復(fù)
281求調(diào)劑（0805）已經(jīng)有8人回復(fù)
環(huán)境領(lǐng)域全國重點實驗室招收博士1-2名已經(jīng)有3人回復(fù)
材料專碩306英一數(shù)二已經(jīng)有10人回復(fù)
301求調(diào)劑已經(jīng)有6人回復(fù)
一志愿天津大學(xué)化學(xué)工藝專業(yè)（081702）315分求調(diào)劑已經(jīng)有7人回復(fù)
302求調(diào)劑已經(jīng)有6人回復(fù)
26博士申請已經(jīng)有3人回復(fù)

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

為什么 Johnson wiley 中chirality的數(shù)據(jù)庫有部分文章是可以免費下的呢？已經(jīng)有4人回復(fù)
求助：PLGA合成中要求壓力較低，1000pa，怎么控制壓力不變呢？已經(jīng)有10人回復(fù)
紫外分光光度計和原子吸收分光光度計單色器位置為什么不同？已經(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ù)
★★★教你圖片如何處理才能達到SCI投稿要求已經(jīng)有61人回復(fù)
請問這個數(shù)據(jù)是否符合要求已經(jīng)有5人回復(fù)
required reviews completed已經(jīng)一個月已經(jīng)有19人回復(fù)
等待Required Reviews Completed 已經(jīng)有27人回復(fù)
SBB 一個月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小時
蟲號: 773188
注冊: 2009-05-18
專業(yè): 理論和計算化學(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小時
蟲號: 805643
注冊: 2009-07-08
專業(yè): 催化化學(xué)

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

這個貌似沒有用到過哦

贊一下(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小時
蟲號: 805643
注冊: 2009-07-08
專業(yè): 催化化學(xué)

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

OUTCAR 里確實有EDIFFG，但是具體怎么體現(xiàn)出來，就像你列出的這個，現(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小時
蟲號: 773188
注冊: 2009-05-18
專業(yè): 理論和計算化學(xué)

★
cenwanglai(金幣+2):呵呵,怎么運行呢?我在運行結(jié)束后含有OUTCAR的目錄下運行$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小時
蟲號: 427638
注冊: 2007-07-29
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

★ ★
cenwanglai(金幣+2):謝謝參與.按你說的意思,是要判斷每一個原子的受力的哈.但是不知道為什么我的任務(wù)會在力沒有達到收斂精度的情況下報告"達到收斂精度,退停止結(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

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

贊一下(1人)

回復(fù)此樓

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

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

stractor

金蟲 (著名寫手)

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

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

EDIFFG=-0.001,收斂標準確實挺高，建議改為-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小時
蟲號: 676900
注冊: 2008-12-18
性別: GG
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

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

贊一下

回復(fù)此樓

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

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

valenhou001

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

1ST強帖: 13
應(yīng)助: 241 (大學(xué)生)
金幣: 25701.7
散金: 602
紅花: 166
帖子: 3782
在線: 873.8小時
蟲號: 1007127
注冊: 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為負數(shù)時，以原子所受的力來判斷收斂；如果為正數(shù)，以前后兩次總能的變化為依據(jù)來判斷收斂。

贊一下(1人)

回復(fù)此樓

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

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

stractor

金蟲 (著名寫手)

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

★ ★ ★ ★
cenwanglai(金幣+1):謝謝參與
cenwanglai(金幣+1): 設(shè)置為-0.001純粹是跑一下看看會出現(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é)束時，哪怕力沒有達到標準，VASP也會reached required accuracy而結(jié)束。
你為什么要設(shè)置-0.001？這個值太小了，費時。-0.01至-0.03夠了。

贊一下(3人)

回復(fù)此樓

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

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

cenwanglai

榮譽版主 (知名作家)

老和山猥瑣派九段

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

這個體系收斂沒有問題了,問題是收斂的結(jié)果不對.
vef.pl 怎么運行呢?運行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 的主題更新

返回列表

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

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

[考研] 0703化學(xué)調(diào)劑，六級已過，有科研經(jīng)歷 +8	曦熙兮 2026-03-15	8/400	2026-03-17 20:31 by xilongliang
[考研] 277調(diào)劑 +5	自由煎餅果子 2026-03-16	6/300	2026-03-17 19:26 by 李leezz
[考研] 341求調(diào)劑 +5	搗蛋豬豬 2026-03-11	7/350	2026-03-17 19:09 by 搗蛋豬豬
[考研] 0854可跨調(diào)劑，一作一項核心論文五項專利，省、國級證書40+數(shù)一英一287 +3	小李0854 2026-03-16	3/150	2026-03-17 13:40 by 熱情沙漠
[考研] 211本，11408一志愿中科院277分，曾在中科院自動化所實習(xí) +6	Losir 2026-03-12	7/350	2026-03-17 12:09 by danranxie
[論文投稿] 有沒有大佬發(fā)小論文能帶我個二作 +3	增銳漏人 2026-03-17	4/200	2026-03-17 09:26 by xs74101122
[考研] 278求調(diào)劑 +3	Yy7400 2026-03-13	3/150	2026-03-17 08:24 by laoshidan
[考研] 070303 總分349求調(diào)劑 +3	LJY9966 2026-03-15	5/250	2026-03-16 14:24 by xwxstudy
[考研] 機械專碩調(diào)劑 +3	笨笨兔子 2026-03-12	3/150	2026-03-15 20:02 by 栗子粥?
[考研] 26考研一志愿中國石油大學(xué)(華東)305分求調(diào)劑 +3	嘉年新程 2026-03-15	3/150	2026-03-15 13:58 by 哈哈哈哈嘿嘿嘿
[考研] 0856專碩279求調(diào)劑 +5	加油加油！? 2026-03-15	5/250	2026-03-15 11:58 by 2020015
[考研] 080500，材料學(xué)碩302分求調(diào)劑學(xué)校 +4	初識可樂 2026-03-14	5/250	2026-03-14 21:08 by peike
[考研] 289求調(diào)劑 +4	這么名字咋樣 2026-03-14	6/300	2026-03-14 18:58 by userper
[考研] 297一志愿上交085600求調(diào)劑 +5	指尖八千里 2026-03-14	5/250	2026-03-14 17:26 by a不易
[考研] 0703化學(xué)調(diào)劑 +4	快樂的香蕉 2026-03-11	4/200	2026-03-13 22:41 by JourneyLucky
[考研] 求材料調(diào)劑 +5	隔壁陳先生 2026-03-12	5/250	2026-03-13 22:03 by 星空星月
[考研] 一志愿西南交大，材料專碩317求調(diào)劑 +5	lx8568 2026-03-11	5/250	2026-03-13 21:43 by peike
[考研] 26調(diào)劑/材料科學(xué)與工程/總分295/求收留 +9	2026調(diào)劑俠 2026-03-12	9/450	2026-03-13 20:46 by 18595523086
[考研] 307求調(diào)劑 +5	超級伊昂大王 2026-03-12	5/250	2026-03-13 15:56 by 棒棒球手
[考研] 295求調(diào)劑 +3	小匕仔汁 2026-03-12	3/150	2026-03-13 15:17 by vgtyfty

24小時熱門版塊排行榜

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

» 猜你喜歡

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

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