| 5 | 1/1 | 返回列表 |
| 查看: 2362 | 回復(fù): 12 | |||
| 當(dāng)前只顯示滿足指定條件的回帖,點(diǎn)擊這里查看本話題的所有回帖 | |||
cenwanglai榮譽(yù)版主 (知名作家)
老和山猥瑣派九段
|
[交流]
【求助成功】從原胞到表面吸附結(jié)構(gòu),一套完整的計(jì)算,怎么規(guī)劃優(yōu)化? 已有5人參與
|
||
|
在各位蟲友的幫助啟發(fā)下,現(xiàn)在VASP總能計(jì)算結(jié)果基本上都能與文獻(xiàn)結(jié)果保持?jǐn)?shù)量級(jí)上的一致.現(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)化需要做嗎?在哪個(gè)階段做比較合適? (2)對(duì)于含有f電子的模型,帶U優(yōu)化結(jié)構(gòu)有必要嗎?先不帶U優(yōu)化結(jié)構(gòu),計(jì)算總能和電子結(jié)構(gòu)的時(shí)候再把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)過了一段時(shí)間,計(jì)算任務(wù)積累的越來越多了,各位怎么管理和保存這些結(jié)算結(jié)果呢? 這些還是入門性的問題,很多可以自己慢慢摸索.但是希望可以和大家交流一下,以節(jié)約時(shí)間.也希望這個(gè)帖子對(duì)一些人有幫助. [ Last edited by cenwanglai on 2010-6-17 at 14:46 ] |
榮譽(yù)版主 (知名作家)
老和山猥瑣派九段
|
候博士的指南中給出了SIGMA優(yōu)化的一個(gè)例子.去找來看看.意思是在0.1-0.3之間設(shè)置幾個(gè)SIGMA,分別做靜態(tài)計(jì)算,然后看EENTRO和SIGMA的關(guān)系,一般要控制EENTRO較小,比如0.01(這個(gè)我不肯定,請(qǐng)高手給予判斷)就可以了. 不過很多文獻(xiàn)直接給出所有計(jì)算SIGMA=0.2不變化,或者算某個(gè)特別的值是指定SIGMA=0.05或者0.07這樣很小的數(shù)值. 但是可以看到,SIGMA對(duì)OUTCAR中的三個(gè)能量數(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)?吹接胓rep "TOTEN" OUTCAR |tail -1 來看能量的. 希望高手可以給個(gè)結(jié)論. |
榮譽(yù)版主 (知名作家)
老和山猥瑣派九段
|
發(fā)現(xiàn)這個(gè)貼字不錯(cuò),大家傳閱. http://www.gaoyang168.com/bbs/viewthread.php?tid=2075970&fpage=1 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿 西北大學(xué) ,070300化學(xué)學(xué)碩,總分287,雙非一本,求調(diào)劑。 +3 | 晨昏線與星海 2026-03-20 | 3/150 |
|
|---|---|---|---|---|
|
[考研] 319求調(diào)劑 +4 | 小力氣珂珂 2026-03-20 | 4/200 |
|
|
[考研] 291 求調(diào)劑 +3 | 化工2026屆畢業(yè)?/a> 2026-03-21 | 3/150 |
|
|
[考研] 286求調(diào)劑 +10 | Faune 2026-03-21 | 10/500 |
|
|
[考研] 333求調(diào)劑 +5 | 87639 2026-03-21 | 7/350 |
|
|
[考研] 求調(diào)劑 +4 | 要好好無聊 2026-03-21 | 4/200 |
|
|
[考研] 材料工程專碩 348分求調(diào)劑 +3 | 冬辭. 2026-03-17 | 5/250 |
|
|
[考研] 材料 271求調(diào)劑 +5 | 展信悅_ 2026-03-21 | 5/250 |
|
|
[考研] 266求調(diào)劑 +3 | 哇呼哼呼哼 2026-03-20 | 3/150 |
|
|
[考研] 22 350 本科985求調(diào)劑,求老登收留 +3 | 李軼男003 2026-03-20 | 3/150 |
|
|
[考研] 一志愿南昌大學(xué),327分,材料與化工085600 +9 | Ncdx123456 2026-03-19 | 9/450 |
|
|
[考研] 330求調(diào)劑 +4 | 小材化本科 2026-03-18 | 4/200 |
|
|
[考研]
|
然11 2026-03-19 | 4/200 |
|
|
[考研] 中南大學(xué)化學(xué)學(xué)碩337求調(diào)劑 +3 | niko- 2026-03-19 | 6/300 |
|
|
[考研] 廣西大學(xué)家禽遺傳育種課題組2026年碩士招生(接收計(jì)算機(jī)專業(yè)調(diào)劑) +3 | 123阿標(biāo) 2026-03-17 | 3/150 |
|
|
[考研] 298-一志愿中國(guó)農(nóng)業(yè)大學(xué)-求調(diào)劑 +9 | 手機(jī)用戶 2026-03-17 | 9/450 |
|
|
[考研] 081700化工學(xué)碩調(diào)劑 +3 | 【1】 2026-03-16 | 3/150 |
|
|
[考研] 收復(fù)試調(diào)劑生 +4 | 雨后秋荷 2026-03-18 | 4/200 |
|
|
[考研] 312求調(diào)劑 +8 | 陌宸希 2026-03-16 | 9/450 |
|
|
[考研] 308求調(diào)劑 +4 | 是Lupa啊 2026-03-16 | 4/200 |
|