版塊導(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í)熱門版塊排行榜

北京石油化工學(xué)院2026年研究生招生接收調(diào)劑公告

返回列表

下一頁(yè)

cenwanglai

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

老和山猥瑣派九段

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

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

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

但是TOTAL-FORCE和total drift都沒(méi)有達(dá)到EDIFFG的要求.PREC=Accurate,關(guān)閉實(shí)空間后,每個(gè)原子每個(gè)方向上的力和total drift還是沒(méi)有達(dá)到要求,計(jì)算就退出來(lái)了,運(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)過(guò)一段時(shí)間觀察發(fā)現(xiàn):EDIFFG<0是對(duì)自由原子(T T T)受力來(lái)做判斷的,對(duì)固定原子(F F F)不做為判斷依據(jù).只要自由原子達(dá)到力的收斂標(biāo)準(zhǔn)就可以了.

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

回復(fù)此樓

» 猜你喜歡

環(huán)境285分，過(guò)六級(jí)，求調(diào)劑已經(jīng)有4人回復(fù)
學(xué)碩化學(xué)工程與技術(shù)，一志愿中國(guó)海洋大學(xué)320+求調(diào)劑已經(jīng)有3人回復(fù)
能源動(dòng)力調(diào)劑已經(jīng)有3人回復(fù)
348求調(diào)劑已經(jīng)有4人回復(fù)
材料化工340求調(diào)劑已經(jīng)有5人回復(fù)
08工科求調(diào)劑290分已經(jīng)有6人回復(fù)
材料調(diào)劑已經(jīng)有3人回復(fù)
372分材料與化工（085600）一志愿湖南大學(xué)求調(diào)劑已經(jīng)有4人回復(fù)
283分材料與化工求調(diào)劑已經(jīng)有4人回復(fù)
一志愿華東理工大學(xué)，080500學(xué)碩，317分，求調(diào)劑已經(jīng)有12人回復(fù)

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

為什么 Johnson wiley 中chirality的數(shù)據(jù)庫(kù)有部分文章是可以免費(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ī)院卻說(shuō)正常已經(jīng)有4人回復(fù)
★★★教你圖片如何處理才能達(dá)到SCI投稿要求已經(jīng)有61人回復(fù)
請(qǐng)問(wèn)這個(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

鐵桿木蟲(chóng) (正式寫手)

應(yīng)助: 7 (幼兒園)
金幣: 7994.9
散金: 1844
紅花: 4
帖子: 776
在線: 569小時(shí)
蟲(chóng)號(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

金蟲(chóng) (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 529.3
散金: 21
紅花: 1
帖子: 106
在線: 51.7小時(shí)
蟲(chóng)號(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è)貌似沒(méi)有用到過(guò)哦

贊一下(1人)

回復(fù)此樓

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

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

jinhuan123

金蟲(chóng) (小有名氣)

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

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

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

贊一下(1人)

回復(fù)此樓

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

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

fanchen021

鐵桿木蟲(chóng) (正式寫手)

應(yīng)助: 7 (幼兒園)
金幣: 7994.9
散金: 1844
紅花: 4
帖子: 776
在線: 569小時(shí)
蟲(chóng)號(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

鐵桿木蟲(chóng) (正式寫手)

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

★ ★
cenwanglai(金幣+2):謝謝參與.按你說(shuō)的意思,是要判斷每一個(gè)原子的受力的哈.但是不知道為什么我的任務(wù)會(huì)在力沒(méi)有達(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ù)是原子收到的力，顯然沒(méi)有達(dá)到收斂標(biāo)準(zhǔn)。
我注意到NSW=10，這個(gè)參數(shù)應(yīng)該變大，比如說(shuō)100

贊一下(1人)

回復(fù)此樓

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

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

stractor

金蟲(chóng) (著名寫手)

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

★
cenwanglai(金幣+2):呵呵,即使設(shè)置成EDIFFG=-0.001,即使實(shí)際上力沒(méi)有達(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

木蟲(chóng) (小有名氣)

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

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

贊一下

回復(fù)此樓

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

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

valenhou001

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

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

贊一下(1人)

回復(fù)此樓

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

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

stractor

金蟲(chóng) (著名寫手)

1ST強(qiáng)帖: 6
應(yīng)助: 91 (初中生)
金幣: 1487.6
散金: 666
紅花: 30
沙發(fā): 1
帖子: 2224
在線: 426.8小時(shí)
蟲(chóng)號(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í)，哪怕力沒(méi)有達(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í)
蟲(chóng)號(hào): 537452
注冊(cè): 2008-04-01
性別: GG
專業(yè): 理論和計(jì)算化學(xué)
管轄: 仿真模擬

這個(gè)體系收斂沒(méi)有問(wèn)題了,問(wèn)題是收斂的結(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 的主題更新

返回列表

下一頁(yè)

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

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

[考研] 一志愿華東理工大學(xué)，080500學(xué)碩，317分，求調(diào)劑 +12	s1145 2026-03-31	12/600	2026-04-02 12:38 by smileboy2006
[考研] 初試301，代碼085701環(huán)境工程，本碩一致，四六級(jí)已過(guò)，有二區(qū)一作，共發(fā)表5篇論文 +5	axibli 2026-04-01	5/250	2026-04-02 12:08 by 1753564080
[考研] 材料工程322分 +8	哈哈哈吼吼吼哈 2026-04-01	8/400	2026-04-02 11:53 by 3041
[考研] 0805求調(diào)劑 +8	是水分 2026-03-31	8/400	2026-04-02 10:46 by guanxin1001
[考研] 材料調(diào)劑 +12	一樣YWY 2026-04-01	12/600	2026-04-02 09:15 by olim
[考研] 291求調(diào)劑 +20	Y-cap 2026-03-29	25/1250	2026-04-01 23:49 by 欣喜777
[考研] 070300化學(xué)279求調(diào)劑 +15	哈哈哈^_^ 2026-03-31	17/850	2026-04-01 21:37 by 給你你注意休息
[考研] 332求調(diào)劑 +8	Lyy930824@ 2026-03-29	8/400	2026-04-01 18:40 by 千鈞澄玉宇啊
[考研] 材料與化工（0856）304求B區(qū)調(diào)劑 +8	邱gl 2026-03-30	16/800	2026-04-01 17:58 by 邱gl
[考研] 294分080500材料科學(xué)與工程求調(diào)劑 +15	柳溪邊 2026-03-26	16/800	2026-04-01 17:07 by 七度不信任
[考研] 267求調(diào)劑 +13	uiybh 2026-03-31	13/650	2026-04-01 10:25 by 探123
[考研] 358求調(diào)劑 +3	王向陽(yáng)花 2026-03-31	3/150	2026-04-01 09:56 by zzchen2000
[考研] 考研調(diào)劑 +9	小蠟新筆 2026-03-29	10/500	2026-03-31 19:52 by Dyhoer
[考研] 0856求調(diào)劑 +9	楒桉 2026-03-28	9/450	2026-03-31 19:06 by 暮澤12
[考研] 本科211總分289，08工學(xué)真心求調(diào)劑 +3	utopiaE 2026-03-30	3/150	2026-03-30 23:42 by ms629
[有機(jī)交流] 甲基亞磺磺酸鈉和甲基磺酸酯反應(yīng)機(jī)理 10+3	kaobao456 2026-03-29	4/200	2026-03-30 23:16 by nBu鋰
[考研] 085701環(huán)境工程求調(diào)劑 +11	多久上課 2026-03-27	12/600	2026-03-30 21:21 by 研究僧導(dǎo)導(dǎo)
[考研] 340求調(diào)劑 +6	Amber00 2026-03-26	6/300	2026-03-29 12:06 by 無(wú)際的草原
[考研] 285求調(diào)劑 +4	AZMK 2026-03-27	7/350	2026-03-27 20:59 by AZMK
[考研] 298調(diào)劑 +3	jiyingjie123 2026-03-27	3/150	2026-03-27 11:57 by wxiongid

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

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

» 猜你喜歡

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

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