| 5 | 1/1 | 返回列表 |
| 查看: 1650 | 回復(fù): 12 | ||||
| 當(dāng)前主題已經(jīng)存檔。 | ||||
| 當(dāng)前只顯示滿足指定條件的回帖,點(diǎn)擊這里查看本話題的所有回帖 | ||||
nkleof銀蟲 (正式寫手)
|
[交流]
【求助】晶格優(yōu)化問題 已有7人參與
|
|||
|
我用vasp優(yōu)化一個體系,所得到的結(jié)果8.36A遠(yuǎn)小于實(shí)驗(yàn)值8.51A。 所采用的贗勢是PAWPBE的,過渡金屬元素Fe和Mn分別試過標(biāo)準(zhǔn)勢和pv勢。采用的優(yōu)化方法分別試過先優(yōu)化原胞體積(ISIF=7),后馳豫離子位置(ISIF=2);以及全優(yōu)化(ISIF=3)。馳豫方式分別試過CG法(IBRION=2)和準(zhǔn)牛頓法(IBRION=1)。為了消除pulay stress,每種方法都以得到的CONTCAR作為POSCAR反復(fù)優(yōu)化了3次。 但是以上所有方法所得到的結(jié)果都是8.36附近,遠(yuǎn)小于實(shí)驗(yàn)值,也小于其它文獻(xiàn)所給出的GGA計(jì)算值8.49。 請經(jīng)驗(yàn)豐富的蟲友幫忙想想可能是在什么地方出問題了? |
VASP |
木蟲 (小有名氣)
| 你在優(yōu)化之前要進(jìn)行K-point和ENCUT的測試。得到一個可靠的K點(diǎn)和ENCUT,然后晶格常數(shù)的求法也有相關(guān)方面的介紹,就是先求出體系能量與晶格之間的關(guān)系然后通過一個方程可以擬合出一條曲線。這個能量最點(diǎn)對應(yīng)的晶格長度就是晶格常數(shù)。還有,我也用過PAWPBE計(jì)算過晶格常數(shù),是要比試驗(yàn)值稍微小一點(diǎn)。希望我的回答對你有幫助 |
銀蟲 (正式寫手)
|
我用的k網(wǎng)格疏密程度是k點(diǎn)之間距離0.09 /A,PREC=Accurate因此ENCUT是贗勢內(nèi)默認(rèn)值的1.3倍,所以這兩個優(yōu)化條件我覺得精度還可以,如果有問題請經(jīng)驗(yàn)豐富的高人再給點(diǎn)建議。 擬合的優(yōu)化正打算做,但是我覺得vasp的自動優(yōu)化功能誤差應(yīng)該不會太大,因此用ISIF=2或3的設(shè)置算出來的值和靜態(tài)計(jì)算擬合優(yōu)化出來的值應(yīng)該不會有太大差異,是不是有可能我的什么設(shè)置出錯? 如下是我做優(yōu)化用的INCAR PREC = Accurate LREAL = A IALGO = 48 NPAR = 8 NSIM = 8 NELMIN = 4 NSW = 200 LORBIT = 12 EDIFF = 1E-5 IBRION = 2 ISIF = 3 ISPIN = 2 MAGMOM = 8*-6 16*6 32*0 [ Last edited by nkleof on 2010-4-19 at 21:27 ] |
木蟲 (著名寫手)
新人菜鳥
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 308求調(diào)劑 +7 | 墨墨漠 2026-03-27 | 7/350 |
|
|---|---|---|---|---|
|
[考研] 295求調(diào)劑 +5 | 1428151015 2026-03-27 | 6/300 |
|
|
[考研] 286求調(diào)劑 +3 | 丟掉懶惰 2026-03-27 | 6/300 |
|
|
[考研] 315調(diào)劑 +4 | 0860求調(diào)劑 2026-03-26 | 5/250 |
|
|
[考研]
|
nnnnnnn5 2026-03-25 | 5/250 |
|
|
[考研] 調(diào)劑 +4 | 柚柚yoyo 2026-03-26 | 4/200 |
|
|
[考研] 求調(diào)劑 +8 | Auroracx 2026-03-22 | 8/400 |
|
|
[考研] 286求調(diào)劑 +13 | Faune 2026-03-21 | 13/650 |
|
|
[考研] 中國科學(xué)院深圳先進(jìn)技術(shù)研究院-光纖傳感課題組招生-中國科學(xué)院大學(xué)、深圳理工大學(xué)聯(lián)培 +5 | YangTyu1 2026-03-26 | 5/250 |
|
|
[考研] 材料與化工304求B區(qū)調(diào)劑 +3 | 邱gl 2026-03-26 | 6/300 |
|
|
[考研] 機(jī)械學(xué)碩310分,數(shù)一英一,一志愿211本科雙非找調(diào)劑信息 +3 | @357 2026-03-25 | 3/150 |
|
|
[考研] 334分 一志愿武理 材料求調(diào)劑 +4 | 李李不服輸 2026-03-26 | 4/200 |
|
|
[考研] 尋找調(diào)劑 +5 | 倔強(qiáng)芒? 2026-03-21 | 8/400 |
|
|
[考研] 302求調(diào)劑 +4 | 錦衣衛(wèi)藤椒 2026-03-25 | 4/200 |
|
|
[考研] 306求0703調(diào)劑一志愿華中師范 +10 | 紙魚ly 2026-03-21 | 11/550 |
|
|
[考研] 一志愿國科過程所081700,274求調(diào)劑 +3 | 三水研0水立方 2026-03-23 | 3/150 |
|
|
[考研] 333求調(diào)劑 +3 | ALULU4408 2026-03-23 | 3/150 |
|
|
[考研] 276求調(diào)劑 +3 | YNRYG 2026-03-21 | 4/200 |
|
|
[考研] 求調(diào)劑院校信息 +6 | CX 330 2026-03-21 | 6/300 |
|
|
[考研] 求調(diào)劑 +4 | 要好好無聊 2026-03-21 | 4/200 |
|