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

登錄注冊

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

24小時熱門版塊排行榜

返回列表

cenwanglai

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

老和山猥瑣派九段

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

[交流] 【求助成功】從原胞到表面吸附結(jié)構(gòu),一套完整的計(jì)算,怎么規(guī)劃優(yōu)化? 已有5人參與

在各位蟲友的幫助啟發(fā)下,現(xiàn)在VASP總能計(jì)算結(jié)果基本上都能與文獻(xiàn)結(jié)果保持?jǐn)?shù)量級上的一致.現(xiàn)在迫切需要提高計(jì)算精度,希望可以把偏差控制到0.02eV以下.我把我的計(jì)算流程列出來,希望大家提供改善或糾正的意見.

我做計(jì)算的流程如下:
(1)原胞優(yōu)化,建模.
(2)ISTART=0,ICHARG=2,EDIFF=1E-3,EDIFFG=-0.1,IBRION=2,LREAL=Auto粗結(jié)構(gòu)優(yōu)化;
(3)rm WAVECAR.EDIFF=1E-4,EDIFFG=-0.01,IBRION=2,LREAL=Auto結(jié)構(gòu)優(yōu)化;
(4)靜態(tài)計(jì)算,優(yōu)化ENCUT,K-points,SIGMA
(5)rm WAVECAR,用優(yōu)化好的參數(shù)和結(jié)構(gòu),做靜態(tài)自恰計(jì)算,得到總能和電荷密度數(shù)據(jù).
(6)能帶等其它性質(zhì)的計(jì)算.
(7)數(shù)據(jù)處理.

問題:
(1)FFT的網(wǎng)格優(yōu)化需要做嗎?在哪個階段做比較合適?
(2)對于含有f電子的模型,帶U優(yōu)化結(jié)構(gòu)有必要嗎?先不帶U優(yōu)化結(jié)構(gòu),計(jì)算總能和電子結(jié)構(gòu)的時候再把U加上滿足精度要求嗎?
(3)同種材料不同結(jié)構(gòu)的模型,SIGMA要分別優(yōu)化嗎?比如原胞總能計(jì)算,需要優(yōu)化SIGMA;做出表面后SIGMA又要優(yōu)化一次;做出表面吸附結(jié)構(gòu)后,SIGMA需要再優(yōu)化一次?提出這樣的問題是因?yàn)槲覚z查同樣的SIGMA設(shè)置,不同模型中EENTRO差別較大,只有少數(shù)情況是較小的(<0.01)
(4)過了一段時間,計(jì)算任務(wù)積累的越來越多了,各位怎么管理和保存這些結(jié)算結(jié)果呢?

這些還是入門性的問題,很多可以自己慢慢摸索.但是希望可以和大家交流一下,以節(jié)約時間.也希望這個帖子對一些人有幫助.

[ Last edited by cenwanglai on 2010-6-17 at 14:46 ]

回復(fù)此樓

» 猜你喜歡

294求調(diào)劑材料與化工專碩已經(jīng)有10人回復(fù)
288求調(diào)劑已經(jīng)有6人回復(fù)
287求調(diào)劑已經(jīng)有4人回復(fù)
招收調(diào)劑碩士已經(jīng)有7人回復(fù)
085601材料工程專碩求調(diào)劑已經(jīng)有8人回復(fù)
一志愿武漢理工材料工程專碩調(diào)劑已經(jīng)有3人回復(fù)
【考研調(diào)劑】化學(xué)專業(yè) 281分，一志愿四川大學(xué)，誠心求調(diào)劑已經(jīng)有3人回復(fù)
281求調(diào)劑（0805）已經(jīng)有19人回復(fù)
材料考研調(diào)劑已經(jīng)有3人回復(fù)
346求調(diào)劑[0856] 已經(jīng)有6人回復(fù)

» 本主題相關(guān)商家推薦: (我也要在這里推廣)

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

表面吸附體系的結(jié)構(gòu)優(yōu)化（弛豫）無法收斂求分析求指導(dǎo)！已經(jīng)有5人回復(fù)
離子吸附結(jié)構(gòu)優(yōu)化后太奇怪了��？？？已經(jīng)有9人回復(fù)
vasp表面和吸附分子的結(jié)構(gòu)分別優(yōu)化后，如何將它們放到一起已經(jīng)有11人回復(fù)
fluent 計(jì)算收斂后，看不到結(jié)果已經(jīng)有4人回復(fù)
關(guān)于在構(gòu)建表面吸附結(jié)構(gòu)時吸附物的位置問題已經(jīng)有5人回復(fù)
EDIFF EDIFFG之間設(shè)置關(guān)系已經(jīng)有12人回復(fù)
vasp結(jié)構(gòu)優(yōu)化中EDIFF以及KPOINTS設(shè)置，以及氫吸附結(jié)合能問題！已經(jīng)有17人回復(fù)
文獻(xiàn)中通常所說的結(jié)構(gòu)優(yōu)化是單胞的優(yōu)化還是對超胞進(jìn)行的優(yōu)化呢？已經(jīng)有5人回復(fù)
緊急求助vasp5.2 計(jì)算氧化鈰表面吸附遇到的若干問題已經(jīng)有4人回復(fù)
Encut和K點(diǎn)個數(shù)之間有沒有聯(lián)系？已經(jīng)有6人回復(fù)
【求助】在用MS搭建表面模型時，怎么看出來撘了幾層模型呢？？已經(jīng)有15人回復(fù)
不算完美的2010（我的年終總結(jié)和展望）已經(jīng)有10人回復(fù)
【求助完畢】vasp-slab-邊界附近原子頂位吸附小分子的結(jié)構(gòu)優(yōu)化很難收斂嗎？已經(jīng)有5人回復(fù)
【求助】結(jié)構(gòu)優(yōu)化時如何將原子位置固定在原胞內(nèi) 已經(jīng)有9人回復(fù)
【原創(chuàng)】關(guān)于結(jié)構(gòu)弛豫設(shè)置及收斂判據(jù)選擇的一點(diǎn)體會已經(jīng)有140人回復(fù)
【求助完畢】硫酸根離子的結(jié)構(gòu)優(yōu)化，偶極修正及勢能修正。已經(jīng)有20人回復(fù)
【求助】ENCUT的意思已經(jīng)有12人回復(fù)
【求助完畢】計(jì)算吸附能的時候要除以原胞的個數(shù),有什么物理意義呢? 已經(jīng)有10人回復(fù)

1樓 2010-05-23 09:48:51

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

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

kgdu

金蟲 (正式寫手)

應(yīng)助: 1 (幼兒園)
金幣: 475.4
紅花: 8
帖子: 396
在線: 295.8小時
蟲號: 893805
注冊: 2009-11-04
性別: GG
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

引用回帖:

Originally posted by lono75 at 2010-05-24 18:19:45:

幾個meV

請問這是對任何體系都是這樣嗎？如果原子數(shù)目較多，是不是應(yīng)該大點(diǎn)呢？這有沒有一個對不同體系都適用的標(biāo)準(zhǔn)？謝謝了！

贊一下(1人)

回復(fù)此樓

7樓2010-05-24 18:28:07

已閱回復(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小時
蟲號: 537452
注冊: 2008-04-01
性別: GG
專業(yè): 理論和計(jì)算化學(xué)
管轄: 仿真模擬

引用回帖:

Originally posted by kgdu at 2010-05-24 16:45:00:

請問優(yōu)化ENCUT，K-point時，一般能量變化小于多少時就可以了，這個值是不是還體系原子的個數(shù)有關(guān)？還請問一下sigma怎么優(yōu)化？謝謝了！！

候博士的指南中給出了SIGMA優(yōu)化的一個例子.去找來看看.意思是在0.1-0.3之間設(shè)置幾個SIGMA,分別做靜態(tài)計(jì)算,然后看EENTRO和SIGMA的關(guān)系,一般要控制EENTRO較小,比如0.01(這個我不肯定,請高手給予判斷)就可以了.

不過很多文獻(xiàn)直接給出所有計(jì)算SIGMA=0.2不變化,或者算某個特別的值是指定SIGMA=0.05或者0.07這樣很小的數(shù)值.

但是可以看到,SIGMA對OUTCAR中的三個能量數(shù)值是有影響的:
freen energy TOTEN
energy without entropy
energy (sigma->0)
到底選什么能量作為計(jì)算能量差計(jì)算基準(zhǔn)?不知道.因?yàn)橛形墨I(xiàn)交代所有能量都插值sigma->0.候博士的指南中好象用的是energy without entropy.但是經(jīng)常看到用grep "TOTEN" OUTCAR |tail -1 來看能量的.

希望高手可以給個結(jié)論.

贊一下(1人)

回復(fù)此樓

9樓2010-05-24 20:43:23

已閱回復(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小時
蟲號: 537452
注冊: 2008-04-01
性別: GG
專業(yè): 理論和計(jì)算化學(xué)
管轄: 仿真模擬

發(fā)現(xiàn)這個貼字不錯,大家傳閱.
http://www.gaoyang168.com/bbs/viewthread.php?tid=2075970&fpage=1

贊一下

回復(fù)此樓

2樓2010-05-24 14:19:12

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

kgdu

金蟲 (正式寫手)

應(yīng)助: 1 (幼兒園)
金幣: 475.4
紅花: 8
帖子: 396
在線: 295.8小時
蟲號: 893805
注冊: 2009-11-04
性別: GG
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

cenwanglai(金幣+2):謝謝參與討論 2010-06-17 14:44:56

引用回帖:

Originally posted by cenwanglai at 2010-05-23 09:48:51:
在各位蟲友的幫助啟發(fā)下,現(xiàn)在VASP總能計(jì)算結(jié)果基本上都能與文獻(xiàn)結(jié)果保持?jǐn)?shù)量級上的一致.現(xiàn)在迫切需要提高計(jì)算精度,希望可以把偏差控制到0.02eV以下.我把我的計(jì)算流程列出來,希望大家提供改善或糾正的意見.

我做 ...

請問優(yōu)化ENCUT，K-point時，一般能量變化小于多少時就可以了，這個值是不是還體系原子的個數(shù)有關(guān)？還請問一下sigma怎么優(yōu)化？謝謝了��！

贊一下

回復(fù)此樓

3樓2010-05-24 16:45:00

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

rdw00

木蟲 (著名寫手)

應(yīng)助: 43 (小學(xué)生)
金幣: 28
散金: 11042
紅花: 23
帖子: 2447
在線: 895.5小時
蟲號: 998186
注冊: 2010-04-16
性別: GG
專業(yè): 催化化學(xué)

cenwanglai(金幣+2):謝謝回答. 2010-05-24 20:25:27

引用回帖:

不能看能量小于多少，是要求收斂。取最佳收斂點(diǎn)就可以

贊一下

回復(fù)此樓

4樓2010-05-24 17:40:38

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

kgdu

金蟲 (正式寫手)

應(yīng)助: 1 (幼兒園)
金幣: 475.4
紅花: 8
帖子: 396
在線: 295.8小時
蟲號: 893805
注冊: 2009-11-04
性別: GG
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

cenwanglai(金幣+1):支持相互交流學(xué)習(xí). 2010-05-24 20:26:02

對于每個給定的ENCUT，計(jì)算收斂后，比較不同ENCUT下的系統(tǒng)能量，看其差值，小于某一值時，就可以取那個較小ENCUT。這個差值一般是多少呢？

贊一下

回復(fù)此樓

5樓2010-05-24 17:56:42

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

lono75

鐵桿木蟲 (著名寫手)

應(yīng)助: 0 (幼兒園)
金幣: 6963.1
散金: 100
紅花: 2
帖子: 2433
在線: 536.4小時
蟲號: 932458
注冊: 2009-12-23
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

cenwanglai(金幣+2):謝謝回復(fù). 2010-05-24 20:26:20

引用回帖:

Originally posted by kgdu at 2010-05-24 17:56:42:
對于每個給定的ENCUT，計(jì)算收斂后，比較不同ENCUT下的系統(tǒng)能量，看其差值，小于某一值時，就可以取那個較小ENCUT。這個差值一般是多少呢？

1meV/atom

[ Last edited by lono75 on 2010-5-24 at 20:13 ]

回復(fù)此樓

由天命,成人事！

6樓2010-05-24 18:19:45

已閱回復(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小時
蟲號: 537452
注冊: 2008-04-01
性別: GG
專業(yè): 理論和計(jì)算化學(xué)
管轄: 仿真模擬

引用回帖:

Originally posted by lono75 at 2010-05-24 18:19:45:

1meV/atom

[ Last edited by lono75 on 2010-5-24 at 20:13 ]

那如果50個原子,總能誤差限就是50atoms*1meV/atom=0.05eV?
對應(yīng)EDIFF=5E-2?

一般設(shè)置是EDIFF=1E-4.怎么解釋?

贊一下

回復(fù)此樓

8樓2010-05-24 20:23:49

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

lono75

鐵桿木蟲 (著名寫手)

應(yīng)助: 0 (幼兒園)
金幣: 6963.1
散金: 100
紅花: 2
帖子: 2433
在線: 536.4小時
蟲號: 932458
注冊: 2009-12-23
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

cenwanglai(金幣+2):明白.謝謝. 2010-05-24 21:18:51
cenwanglai(金幣+2):謝謝回復(fù),再鼓勵一下. 2010-05-25 09:21:52

兩個不是一回事

EDIFF是取某一個Ecut時,電子自洽迭代時,相鄰兩次的能量差,達(dá)到此值時,ION動一下,ION步對應(yīng)EDIFFG,.....

ΔEi<0.001eV/atom for different Ecut

電子迭代的ΔEi->EDIFF not the ΔEi of the later

贊一下

回復(fù)此樓

由天命,成人事！

10樓2010-05-24 21:08:58

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

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

返回列表

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

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

[考研] 能源材料化學(xué)課題組招收碩士研究生8-10名 +4	脫穎而出 2026-03-16	11/550	2026-03-19 08:11 by 脫穎而出
[考研] 材料專碩英一數(shù)二306 +5	z1z2z3879 2026-03-18	5/250	2026-03-19 07:43 by BruceLiu320
[考研] 0703化學(xué)調(diào)劑 +8	妮妮ninicgb 2026-03-15	12/600	2026-03-19 05:55 by anny19840123
[考研] 085700資源與環(huán)境308求調(diào)劑 +3	墨墨漠 2026-03-18	3/150	2026-03-18 22:35 by bingxueer79
[考研] 085600材料與化工 +5	安全上岸！ 2026-03-16	5/250	2026-03-18 15:33 by cmz0325
[考研] 286求調(diào)劑 +6	lemonzzn 2026-03-16	9/450	2026-03-18 15:31 by vgtyfty
[考研] 0703化學(xué)調(diào)劑，六級已過，有科研經(jīng)歷 +10	曦熙兮 2026-03-15	10/500	2026-03-18 14:19 by 007_lilei
[考研] 070300化學(xué)319求調(diào)劑 +6	錦鯉0909 2026-03-17	6/300	2026-03-18 13:22 by Iveryant
[考博] 環(huán)境領(lǐng)域全國重點(diǎn)實(shí)驗(yàn)室招收博士1-2名 +3	QGZDSYS 2026-03-13	5/250	2026-03-18 11:13 by QGZDSYS
[考研] 材料，紡織，生物（0856、0710），化學(xué)招生啦 +3	Eember. 2026-03-17	9/450	2026-03-18 10:28 by Eember.
[考研] 334求調(diào)劑 +3	志存高遠(yuǎn)意在機(jī)?/a> 2026-03-16	3/150	2026-03-18 08:34 by lm4875102
[碩博家園] 湖北工業(yè)大學(xué) 生命科學(xué)與健康學(xué)院-課題組招收2026級食品/生物方向碩士 +3	1喜春8 2026-03-17	5/250	2026-03-17 17:18 by ber川cool子
[考研] 一志愿蘇州大學(xué)材料工程（085601）專碩有科研經(jīng)歷三項(xiàng)國獎兩個實(shí)用型專利一項(xiàng)省級立項(xiàng) +6	大火山小火山 2026-03-16	8/400	2026-03-17 15:05 by 無懈可擊111
[考研] 278求調(diào)劑 +3	Yy7400 2026-03-13	3/150	2026-03-17 08:24 by laoshidan
[考研] 283求調(diào)劑 +3	聽風(fēng)就是雨； 2026-03-16	3/150	2026-03-17 07:41 by 熱情沙漠
[考研] 0854控制工程 359求調(diào)劑可跨專業(yè) +3	626776879 2026-03-14	9/450	2026-03-16 17:42 by 626776879
[考研] 321求調(diào)劑 +5	大米飯！ 2026-03-15	5/250	2026-03-16 16:33 by houyaoxu
[考研] 085600材料與化工求調(diào)劑 +13	enenenhui 2026-03-13	14/700	2026-03-16 15:19 by 了了了了。。
[考研] 中科大材料專碩319求調(diào)劑 +3	孟鑫材料 2026-03-13	3/150	2026-03-14 18:10 by houyaoxu
[考研] 26調(diào)劑/材料科學(xué)與工程/總分295/求收留 +9	2026調(diào)劑俠 2026-03-12	9/450	2026-03-13 20:46 by 18595523086

24小時熱門版塊排行榜

[交流] 【求助成功】從原胞到表面吸附結(jié)構(gòu),一套完整的計(jì)算,怎么規(guī)劃優(yōu)化? 已有5人參與

» 猜你喜歡

» 本主題相關(guān)商家推薦: (我也要在這里推廣)

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

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