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

返回列表

xianglannihao

銅蟲 (小有名氣)

應(yīng)助: 1 (幼兒園)
金幣: 1368.5
散金: 8
紅花: 4
帖子: 169
在線: 392.7小時(shí)
蟲號(hào): 267091
注冊(cè): 2006-07-22
性別: MM
專業(yè): 理論和計(jì)算化學(xué)

[交流] 關(guān)于VASP收斂判斷

關(guān)于VASP收斂判斷：
假定設(shè)置如下：
EDIFF = 1E-05    (SCF convergence criteria)
EDIFFG = -0.03 (Ionic convergence criteria, eV/AA)

我的理解收斂條件：EDIFF 是電子步收斂的條件，電子步dE達(dá)到1E-5就認(rèn)為收斂。
EDIFFG是離子步收斂的標(biāo)準(zhǔn)，為負(fù)數(shù)時(shí)表示力的收斂判據(jù)。是否收斂看OUTCAR中的 reached required accuracy - stopping structural energy minimisation提示和最后一輪的如下兩列：
d Force = 0.1423097E-03[ 0.108E-03, 0.177E-03]  d Energy = 0.1531149E-03-0.108E-04
d Force =-0.6844216E+00[-0.684E+00,-0.684E+00]  d Ewald  =-0.6844216E+00-0.177E-07

問題：1. 是否有這個(gè)提示就表示收斂了？
2. OUTCAR中為什么有兩列d Force, 都表示什么意思？d Force d Energy表示的都是一個(gè)范圍？
3. 有時(shí)在每一輪的離子步計(jì)算時(shí)，其電子步dE沒達(dá)到也顯示該輪收斂。
4. 優(yōu)化好的構(gòu)型，以相同的設(shè)置再重新優(yōu)化一次，我發(fā)現(xiàn)其能量會(huì)發(fā)生變化，有時(shí)會(huì)大于1eV，但構(gòu)型變化倒不是很大？難道我們要一直優(yōu)化，直至能量不變，是否這樣才能得到可靠的計(jì)算結(jié)果？

回復(fù)此樓

» 猜你喜歡

一志愿武理材料305分求調(diào)劑已經(jīng)有5人回復(fù)
材料專碩英一數(shù)二306 已經(jīng)有3人回復(fù)
085600材料與化工調(diào)劑 324分已經(jīng)有7人回復(fù)
311求調(diào)劑已經(jīng)有3人回復(fù)
26調(diào)劑/材料/英一數(shù)二/總分289/已過A區(qū)線已經(jīng)有7人回復(fù)
295求調(diào)劑已經(jīng)有5人回復(fù)
354求調(diào)劑已經(jīng)有5人回復(fù)
材料專業(yè)求調(diào)劑已經(jīng)有3人回復(fù)
求材料調(diào)劑已經(jīng)有7人回復(fù)
一志愿天大材料與化工（085600）總分338 已經(jīng)有3人回復(fù)

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

關(guān)于VASP計(jì)算收斂的判斷已經(jīng)有25人回復(fù)
關(guān)于能帶自洽計(jì)算是否完成的判斷已經(jīng)有4人回復(fù)
【原創(chuàng)】關(guān)于結(jié)構(gòu)弛豫設(shè)置及收斂判據(jù)選擇的一點(diǎn)體會(huì) 已經(jīng)有140人回復(fù)
【求助】Dmol3 K 收斂性判斷的問題已經(jīng)有6人回復(fù)

1樓 2015-10-05 11:19:44

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

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

shi398914439

木蟲 (正式寫手)

應(yīng)助: 15 (小學(xué)生)
金幣: 2290.3
散金: 100
紅花: 8
帖子: 666
在線: 228.3小時(shí)
蟲號(hào): 2335618
注冊(cè): 2013-03-10
性別: GG
專業(yè): 金屬材料的制備科學(xué)與跨學(xué)

★
小木蟲: 金幣+0.5, 給個(gè)紅包，謝謝回帖

難得的干貨，這些問題我也都遇到了……關(guān)于收斂標(biāo)準(zhǔn)，我學(xué)的是在total force后有個(gè)total drift后的三列都小于力的收斂，另外出現(xiàn)計(jì)算時(shí)間，說明整個(gè)計(jì)算過程是正確運(yùn)行的！
  關(guān)于多次優(yōu)化的問題，我也會(huì)這樣，尤其當(dāng)你的構(gòu)形越接近最優(yōu)構(gòu)形，你優(yōu)化次數(shù)越多，偏差越大，甚至構(gòu)形還會(huì)散或不收斂……一直也沒解決，我理解的計(jì)算原理這樣解釋：軟件在計(jì)算時(shí)會(huì)先進(jìn)行離子步移動(dòng)，當(dāng)你接近最優(yōu)構(gòu)形時(shí)，如果離子步過大，超過了收斂標(biāo)準(zhǔn)的閾值，而離子步移動(dòng)具有一定隨機(jī)性，因此可能導(dǎo)致更大的偏離最優(yōu)構(gòu)形！
  vasp本身最優(yōu)構(gòu)形就是相對(duì)的，精確的用聲子譜。優(yōu)化一段后最好減小離子步幅……
  當(dāng)然，以上都是我的個(gè)人理解，權(quán)當(dāng)拋磚引玉，坐等權(quán)威解釋……

發(fā)自小木蟲Android客戶端

贊一下(2人)

回復(fù)此樓

10樓2015-10-07 11:18:31

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

wmy8802217

木蟲 (正式寫手)

應(yīng)助: 89 (初中生)
金幣: 2877
紅花: 12
帖子: 350
在線: 461.9小時(shí)
蟲號(hào): 2935682
注冊(cè): 2014-01-16
性別: GG
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

★ ★ ★ ★
小木蟲: 金幣+0.5, 給個(gè)紅包，謝謝回帖
xianglannihao: 金幣+3 2015-10-05 12:19:58

EDIFFG應(yīng)該說的是每個(gè)原子的受力收斂吧，應(yīng)該看OUTCAR里的Total-Force；
EDIFF我覺得應(yīng)該是D eps和dE的收斂；
至于第四個(gè)問題，你應(yīng)該加大NSW，增大離子步，因?yàn)槟阒暗慕Y(jié)果應(yīng)該是沒收斂的。

贊一下(1人)

回復(fù)此樓

凝聚態(tài)理論與計(jì)算物理

2樓2015-10-05 12:02:52

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

普通回帖

xianglannihao

銅蟲 (小有名氣)

應(yīng)助: 1 (幼兒園)
金幣: 1368.5
散金: 8
紅花: 4
帖子: 169
在線: 392.7小時(shí)
蟲號(hào): 267091
注冊(cè): 2006-07-22
性別: MM
專業(yè): 理論和計(jì)算化學(xué)

引用回帖:

2樓: Originally posted by wmy8802217 at 2015-10-05 12:02:52
EDIFFG應(yīng)該說的是每個(gè)原子的受力收斂吧，應(yīng)該看OUTCAR里的Total-Force；
EDIFF我覺得應(yīng)該是D eps和dE的收斂；
至于第四個(gè)問題，你應(yīng)該加大NSW，增大離子步，因?yàn)槟阒暗慕Y(jié)果應(yīng)該是沒收斂的。

謝謝回復(fù)。
1. OUTCAR里的每個(gè)原子在三個(gè)文獻(xiàn)上Total-Force我看了，除了設(shè)置了固定原子位置的原子，其它原子的力都小于0.03，已達(dá)到力收斂標(biāo)準(zhǔn)。
2. EDIFF只表示其中一個(gè)吧？
3. NSW = 500,已經(jīng)足夠大了吧，之前的結(jié)果都提示已達(dá)收斂。除了設(shè)置了固定原子位置的原子，其它原子的力都小于0.03，已達(dá)到力收斂標(biāo)準(zhǔn)。

贊一下

回復(fù)此樓

3樓2015-10-05 12:19:47

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

wmy8802217

木蟲 (正式寫手)

應(yīng)助: 89 (初中生)
金幣: 2877
紅花: 12
帖子: 350
在線: 461.9小時(shí)
蟲號(hào): 2935682
注冊(cè): 2014-01-16
性別: GG
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

★ ★ ★ ★
小木蟲: 金幣+0.5, 給個(gè)紅包，謝謝回帖
xianglannihao: 金幣+3 2015-10-05 16:19:06

引用回帖:

3樓: Originally posted by xianglannihao at 2015-10-05 12:19:47
謝謝回復(fù)。
1. OUTCAR里的每個(gè)原子在三個(gè)文獻(xiàn)上Total-Force我看了，除了設(shè)置了固定原子位置的原子，其它原子的力都小于0.03，已達(dá)到力收斂標(biāo)準(zhǔn)。
2. EDIFF只表示其中一個(gè)吧？
3. NSW = 500,已經(jīng)足夠大 ...

我覺得你可以把標(biāo)準(zhǔn)再調(diào)低一點(diǎn)看看，既然已經(jīng)到了收斂標(biāo)準(zhǔn)，沒必要再來一次；（也許對(duì)于你這個(gè)力收斂精度，對(duì)于總能的影響可能就在1eV左右）
比較能量你可以用更大的K點(diǎn)和截?cái)嗄苡?jì)算更精確的總能來比較。

贊一下

回復(fù)此樓

凝聚態(tài)理論與計(jì)算物理

4樓2015-10-05 13:11:00

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

xianglannihao

銅蟲 (小有名氣)

應(yīng)助: 1 (幼兒園)
金幣: 1368.5
散金: 8
紅花: 4
帖子: 169
在線: 392.7小時(shí)
蟲號(hào): 267091
注冊(cè): 2006-07-22
性別: MM
專業(yè): 理論和計(jì)算化學(xué)

引用回帖:

4樓: Originally posted by wmy8802217 at 2015-10-05 13:11:00
我覺得你可以把標(biāo)準(zhǔn)再調(diào)低一點(diǎn)看看，既然已經(jīng)到了收斂標(biāo)準(zhǔn)，沒必要再來一次；（也許對(duì)于你這個(gè)力收斂精度，對(duì)于總能的影響可能就在1eV左右）
比較能量你可以用更大的K點(diǎn)和截?cái)嗄苡?jì)算更精確的總能來比較。...

只是很好奇，為什么優(yōu)化后的構(gòu)型，再次以相同設(shè)置優(yōu)化按道理應(yīng)該不變化才對(duì)吧？你試過拿著優(yōu)化好的構(gòu)型再優(yōu)化嗎？

贊一下

回復(fù)此樓

5樓2015-10-05 16:21:40

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

wmy8802217

木蟲 (正式寫手)

應(yīng)助: 89 (初中生)
金幣: 2877
紅花: 12
帖子: 350
在線: 461.9小時(shí)
蟲號(hào): 2935682
注冊(cè): 2014-01-16
性別: GG
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

★ ★ ★
小木蟲: 金幣+0.5, 給個(gè)紅包，謝謝回帖
xianglannihao: 金幣+2 2015-10-07 10:30:44

引用回帖:

5樓: Originally posted by xianglannihao at 2015-10-05 16:21:40
只是很好奇，為什么優(yōu)化后的構(gòu)型，再次以相同設(shè)置優(yōu)化按道理應(yīng)該不變化才對(duì)吧？你試過拿著優(yōu)化好的構(gòu)型再優(yōu)化嗎？...

再優(yōu)化跑了幾個(gè)離子步？

贊一下

回復(fù)此樓

凝聚態(tài)理論與計(jì)算物理

6樓2015-10-05 19:20:27

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

stl409

7樓2015-10-05 19:56:40

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

xianglannihao

銅蟲 (小有名氣)

應(yīng)助: 1 (幼兒園)
金幣: 1368.5
散金: 8
紅花: 4
帖子: 169
在線: 392.7小時(shí)
蟲號(hào): 267091
注冊(cè): 2006-07-22
性別: MM
專業(yè): 理論和計(jì)算化學(xué)

引用回帖:

6樓: Originally posted by wmy8802217 at 2015-10-05 19:20:27
再優(yōu)化跑了幾個(gè)離子步？...

謝謝您的討論。
我的測(cè)試結(jié)果如下：
構(gòu)型1, 優(yōu)化一次,能量為-219.837378; 以第一次優(yōu)化得到的構(gòu)型,以相同設(shè)置再次優(yōu)化,能量為-220.206788; 依此再優(yōu)化能量為-219.12072;

構(gòu)型2, 優(yōu)化一次,能量為-220.292638; 以第一次優(yōu)化得到的構(gòu)型,以相同設(shè)置再次優(yōu)化,能量為-219.157869; 依此再優(yōu)化能量為-218.285171;

構(gòu)型1和模型分別是同種分子不同吸附構(gòu)型，相對(duì)穩(wěn)定是如果按第一次優(yōu)化的結(jié)果是構(gòu)型2更穩(wěn)定，而按后兩次優(yōu)化的結(jié)果是構(gòu)型1穩(wěn)定。而按理論的預(yù)計(jì)是應(yīng)該構(gòu)型1穩(wěn)定的。所以我在想到底是什么問題造成的差別，思考如何得到可靠結(jié)果的問題。

再次優(yōu)化一般又會(huì)優(yōu)化10F至20F左右.

贊一下

回復(fù)此樓

8樓2015-10-07 10:31:19

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

wmy8802217

木蟲 (正式寫手)

應(yīng)助: 89 (初中生)
金幣: 2877
紅花: 12
帖子: 350
在線: 461.9小時(shí)
蟲號(hào): 2935682
注冊(cè): 2014-01-16
性別: GG
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

★
小木蟲: 金幣+0.5, 給個(gè)紅包，謝謝回帖

引用回帖:

8樓: Originally posted by xianglannihao at 2015-10-07 10:31:19
謝謝您的討論。
我的測(cè)試結(jié)果如下：
構(gòu)型1, 優(yōu)化一次,能量為-219.837378; 以第一次優(yōu)化得到的構(gòu)型,以相同設(shè)置再次優(yōu)化,能量為-220.206788; 依此再優(yōu)化能量為-219.12072;

構(gòu)型2, 優(yōu)化一次,能量為-220.292638; ...

你應(yīng)該試試加大精度。。
我覺得有可能是因?yàn)樵谀隳壳?.03eV/A的精度下，就會(huì)照成1eV左右的誤差，加大精度試試。
打個(gè)比方：一個(gè)方程F(x)=0的精確解是5，但是F(4.5-5.5)的值都小于0.01，你如果把精度定到0.01，在數(shù)值求解的時(shí)候就會(huì)有誤差，辦法就是加大精度到0.001試試。

贊一下

回復(fù)此樓

凝聚態(tài)理論與計(jì)算物理

9樓2015-10-07 10:37:00

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

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

返回列表

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

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

[考研] 材料專碩英一數(shù)二306 +3	z1z2z3879 2026-03-18	3/150	2026-03-18 17:18 by yeahyou
[考研] 295求調(diào)劑 +3	一志愿京區(qū)211 2026-03-18	5/250	2026-03-18 17:03 by zhaoqian0518
[考研] 材料專業(yè)求調(diào)劑 +3	hanamiko 2026-03-18	3/150	2026-03-18 16:30 by 30660438
[考研] 344求調(diào)劑 +5	knight344 2026-03-16	6/300	2026-03-18 15:17 by knight344
[考研] 化學(xué)工程321分求調(diào)劑 +15	大米飯！ 2026-03-15	18/900	2026-03-18 14:52 by haxia
[考研] 材料與化工一志愿南昌大學(xué)327求調(diào)劑推薦 +8	Ncdx123456 2026-03-13	9/450	2026-03-18 14:40 by haxia
[考研] 312求調(diào)劑 +8	陌宸希 2026-03-16	9/450	2026-03-18 12:39 by Linda Hu
[考研] 280求調(diào)劑 +6	咕嚕曉曉 2026-03-18	7/350	2026-03-18 11:25 by 無際的草原
[考研] 0703化學(xué)調(diào)劑 +3	妮妮ninicgb 2026-03-17	3/150	2026-03-18 10:29 by macy2011
[考研] 268求調(diào)劑 +8	一定有學(xué)上- 2026-03-14	9/450	2026-03-17 17:47 by laoshidan
[考研] 302求調(diào)劑 +9	負(fù)心者當(dāng)誅 2026-03-11	9/450	2026-03-17 17:13 by ruiyingmiao
[考研] 085601求調(diào)劑 +4	Du.11 2026-03-16	4/200	2026-03-17 17:08 by ruiyingmiao
[考研] 290求調(diào)劑 +3	p asserby. 2026-03-15	4/200	2026-03-17 16:35 by wangkm
[考研] 化學(xué)調(diào)劑0703 +8	啊我我的 2026-03-11	8/400	2026-03-16 17:23 by 我的船我的海
[考研] 085600調(diào)劑 +5	漾漾123sun 2026-03-12	6/300	2026-03-16 15:58 by 漾漾123sun
[考研] 070303 總分349求調(diào)劑 +3	LJY9966 2026-03-15	5/250	2026-03-16 14:24 by xwxstudy
[考研] 0856求調(diào)劑 +3	劉夢(mèng)微 2026-03-15	3/150	2026-03-16 10:00 by houyaoxu
[考研] 26考研一志愿中國石油大學(xué)(華東)305分求調(diào)劑 +3	嘉年新程 2026-03-15	3/150	2026-03-15 13:58 by 哈哈哈哈嘿嘿嘿
[基金申請(qǐng)] 有必要更換申報(bào)口嗎 20+3	fannyamoy 2026-03-11	3/150	2026-03-14 00:52 by zhanghaozhu
[考研] 307求調(diào)劑 +5	超級(jí)伊昂大王 2026-03-12	5/250	2026-03-13 15:56 by 棒棒球手

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

[交流] 關(guān)于VASP收斂判斷

» 猜你喜歡

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

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