| 3 | 1/1 | 返回列表 |
| 查看: 1128 | 回復(fù): 2 | ||||
fangyongxinxi新蟲 (正式寫手)
|
[求助]
為什么vasp能一直向高能量方向優(yōu)化? 已有1人參與
|
|
為什么vasp能一直向高能量方向優(yōu)化? 自己打算對(duì)一個(gè)體系以同樣參數(shù),進(jìn)行結(jié)構(gòu)多次優(yōu)化,ISIF=3,其中晶格常數(shù)abc如下,后面的數(shù)字是經(jīng)過(guò)優(yōu)化得到的,第一列是初始數(shù)據(jù)(別人文獻(xiàn)中查找的),第二列是第二次優(yōu)化,第三列是第三次優(yōu)化結(jié)果,依次類推。能量:第一次優(yōu)化的能量 -55.64ev, 第二次,-53.53ev ,第三次,-52.22ev,第四次,-51.44ev。能量隨著優(yōu)化次數(shù)逐漸升高。 a 5.75 5.62 5.56 5.52 5.51 b 4.06 3.97 3.92 3.89 3.88 c 4.06 3.97 3.92 3.89 3.88 求高手解答: 1. 為什么會(huì)出現(xiàn)優(yōu)化向高能量方向發(fā)生這樣的情況?(一般情況是如果一個(gè)體系達(dá)到能量和受力標(biāo)準(zhǔn),就會(huì)停止優(yōu)化,比如MS,為什么vaso能一直優(yōu)化?一種解釋是這個(gè)體系有亞穩(wěn)態(tài)存在,但能量趨勢(shì)如何解釋?自己將嘗試用同樣的方法優(yōu)化一些簡(jiǎn)單體系,比如cu ,si 等,看是否有同樣結(jié)果) 2. 就這個(gè)體系來(lái)說(shuō),自己應(yīng)該選哪個(gè)作為性質(zhì)計(jì)算的初始結(jié)構(gòu)?第一次優(yōu)化的,還是后面的,因?yàn)榈谝淮蝺?yōu)化的能量最小。 后面自己又進(jìn)行了嘗試:1. 首先設(shè)定ISIF=2,優(yōu)化,然后再用isif=3優(yōu)化,得到結(jié)果與上面直接用ISIF=3結(jié)果一樣;2. 首先用ISIF=4,再用ISIF=3,得到結(jié)果和上面直接用ISIF=3一樣。所以后面得出結(jié)論:對(duì)于本體系,直接用ISIF=3優(yōu)化一次就行了。這樣做對(duì)嗎? 謝謝。 |

專家顧問(wèn) (著名寫手)
![]() |
專家經(jīng)驗(yàn): +20 |
|
1. 為什么會(huì)出現(xiàn)優(yōu)化向高能量方向發(fā)生這樣的情況?(一般情況是如果一個(gè)體系達(dá)到能量和受力標(biāo)準(zhǔn),就會(huì)停止優(yōu)化,比如MS,為什么vaso能一直優(yōu)化?一種解釋是這個(gè)體系有亞穩(wěn)態(tài)存在,但能量趨勢(shì)如何解釋?自己將嘗試用同樣的方法優(yōu)化一些簡(jiǎn)單體系,比如cu ,si 等,看是否有同樣結(jié)果) > VASP 就如你想的,就像一般的軟體 (比如 MS) 一樣,結(jié)構(gòu)優(yōu)化是會(huì)往能量低的方向前進(jìn) ISIF=3 通常在 1~2 次就會(huì)找到最低能量,之後再怎麼重覆做,結(jié)果都會(huì)是一樣,會(huì)發(fā)生像你遇到的情況,我能想到以下幾點(diǎn) (1) 你的計(jì)算沒(méi)有 "收斂" VASP 裡面有個(gè) electron loop convergence criterion,叫 NELM。default value 是 60 如果你electron loop 到達(dá)這次數(shù),VASP 會(huì)?娭仆V,?dǎo)致計(jì)算不收斂,其結(jié)果是不可信的 (2) 你的 VASP 沒(méi)有正確編輯 有可能你使用的 math kernel library 有 bug,或是編輯上有參數(shù)沒(méi)有給正確,導(dǎo)致計(jì)算結(jié)果是不可信 (3) 你 INCAR 內(nèi)可能有加了些特殊的參數(shù),導(dǎo)致計(jì)算不收斂,或是收斂到不正確的地方 在沒(méi)有看到 INCAR, POSCAR, POTCAR, KPOINTS 前,很難判斷發(fā)生什麼事 但能肯定的是,正常的 VASP 計(jì)算,是不會(huì)往能量高的地方去收斂 2. 就這個(gè)體系來(lái)說(shuō),自己應(yīng)該選哪個(gè)作為性質(zhì)計(jì)算的初始結(jié)構(gòu)?第一次優(yōu)化的,還是后面的,因?yàn)榈谝淮蝺?yōu)化的能量最小。 > 要先確認(rèn)計(jì)算是"正確" |
| 3 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 調(diào)劑 +4 | 好好讀書。 2026-03-28 | 6/300 |
|
|---|---|---|---|---|
|
[考研] 考研調(diào)劑 +9 | 小蠟新筆 2026-03-29 | 10/500 |
|
|
[考研] 復(fù)試調(diào)劑 +7 | 雙馬尾痞老板2 2026-03-31 | 7/350 |
|
|
[考研] 一志愿華東師范大學(xué)有機(jī)化學(xué)專業(yè),初試351分,復(fù)試被刷求調(diào)劑! +9 | 真名有冰 2026-03-29 | 10/500 |
|
|
[考研] 322求調(diào)劑:一志愿湖南大學(xué) 材料與化工(085600),已過(guò)六級(jí)。 +10 | XX小鄧 2026-03-29 | 10/500 |
|
|
[考研] 求化學(xué)調(diào)劑 +12 | wulanna 2026-03-28 | 12/600 |
|
|
[考研] 322求調(diào)劑 +3 | 熹僖XX 2026-03-31 | 3/150 |
|
|
[考研] 08工科,295,接受跨專業(yè)調(diào)劑 +4 | lmnlzy 2026-03-31 | 4/200 |
|
|
[考研] 一志愿哈爾濱工業(yè)大學(xué)材料與化工方向336分 +13 | 辰沐5211314 2026-03-26 | 13/650 |
|
|
[考研] 食品工程專碩一志愿中海洋309求調(diào)劑 +5 | 小張zxy張 2026-03-26 | 10/500 |
|
|
[考研] 343求調(diào)劑085601 +3 | 要努力學(xué)習(xí)x 2026-03-29 | 3/150 |
|
|
[考研] 085600,專業(yè)課化工原理,321分求調(diào)劑 +5 | 大饞小子 2026-03-28 | 5/250 |
|
|
[考研] 085600,材料與化工321分調(diào)劑 +4 | 大饞小子 2026-03-27 | 6/300 |
|
|
[考研] 315調(diào)劑 +4 | 0860求調(diào)劑 2026-03-26 | 5/250 |
|
|
[考研] 292求調(diào)劑 +4 | 求求了收下我吧?/a> 2026-03-26 | 4/200 |
|
|
[考研] 286求調(diào)劑 +4 | lim0922 2026-03-26 | 4/200 |
|
|
[考研] 081200-11408-276學(xué)碩求調(diào)劑 +4 | 崔wj 2026-03-26 | 4/200 |
|
|
[考研] 打過(guò)很多競(jìng)賽,085406控制工程300分,求調(diào)劑 +3 | askeladz 2026-03-26 | 3/150 |
|
|
[考研] 網(wǎng)絡(luò)空間安全0839招調(diào)劑 +4 | w320357296 2026-03-25 | 6/300 |
|
|
[考研] 302求調(diào)劑 +4 | 錦衣衛(wèi)藤椒 2026-03-25 | 4/200 |
|