| 8 | 1/1 | 返回列表 |
| 查看: 5262 | 回復(fù): 7 | ||||
Prettyswan木蟲 (小有名氣)
|
[交流]
【求助】結(jié)構(gòu)優(yōu)化中,如何選擇POTIM? 已有6人參與
|
|
說明手冊(cè)上說: 這里推薦在結(jié)構(gòu)優(yōu)化中使用IBRION=2和POTIM=0.5. 在侯柱峰編寫的《vasp軟件包使用入門指南》中《§6.2結(jié)構(gòu)參數(shù)的優(yōu)化》中“一次成型”一節(jié)寫道: 那么應(yīng)該如何選擇POTIM的值呢? |
VASP |
金蟲 (著名寫手)
木蟲 (小有名氣)
新蟲 (小有名氣)
|
我現(xiàn)在也在優(yōu)化晶體的結(jié)構(gòu)參數(shù)?戳撕蟛┦康闹形娜腴T手冊(cè),雖然講的很詳細(xì),當(dāng)有個(gè)地方還不是很明白。就他例子的關(guān)于Mg的晶格參數(shù)的優(yōu)化(P25)。 我取10個(gè)點(diǎn),循環(huán)完后,擬合得到最佳的a值。他后面說“在該體積下,重復(fù)上面第2和3步,得到平衡狀態(tài)下的其他晶胞參數(shù)。這一步就是:在得到了E(V)曲線后,通過擬合狀態(tài)方程得到平衡狀態(tài)下的體積,計(jì)算出上面腳本中的變量$i的值,并改變$i的循環(huán)值,再運(yùn)行run_cell計(jì)算一次,得到其他的結(jié)構(gòu)參數(shù)c和位置u.! 對(duì)這句話我不太理解,安我自己的意思,通過擬合上面E(V)曲線得到平衡狀態(tài)下的體積,計(jì)算出變量$i,然后在該值下遲豫,就可以得到b,c。而不是候博士上面說的“改變$i的循環(huán)值,再運(yùn)行run_cell計(jì)算一次”。 不知道大家怎么理解他這句話的。 |
至尊木蟲 (職業(yè)作家)
|
“改變$i的循環(huán)值,再運(yùn)行run_cell計(jì)算一次” 此時(shí)$i只附了一個(gè)值,就是前面由EOS擬合得到的平衡體積V_0計(jì)算(V_0 = 0.5*sqrt(3.0) a^3 * c/a, c/a的值取自run_cell所設(shè)置的,即c/a的值與前一步中run_cell中的相同。這樣得到了一個(gè)平衡體積下初始的a, c/a)。這一次run_cell就是固定體積(將體積固定在所擬合出來的V_0),根據(jù)run_cell所設(shè)置的初始a, c/a進(jìn)行固定體積的優(yōu)化。 E(a,c) -----> E(V, a)------->E(V_0,a)------> V_0下最優(yōu)的a, c的值。 |
|
謝謝斑竹的回答,我理解了這個(gè)意思。當(dāng)我還有點(diǎn)疑問,vasp 好像是在開始十步優(yōu)化的時(shí)候,都是固定角度和c/a的值進(jìn)行優(yōu)化的。萬一我優(yōu)化的是一個(gè)三斜晶系,我該怎么辦。我還是安你的這個(gè)步驟:E(a,b,c)->E(v,a)->E(V_0,a)->V_0下的最佳a(bǔ),b,c,alpha,beta,gama嗎?不知道這個(gè)有通用步驟沒有? 因?yàn)槲疫想優(yōu)化結(jié)構(gòu):a。絙!=c,三個(gè)角度都為90度。如果還是按Mg那種方法的話,b/a的值將會(huì)始終不變。 [ Last edited by swx0789 on 2010-12-25 at 15:22 ] |
銀蟲 (正式寫手)
| 8 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 293求調(diào)劑 +9 | zjl的號(hào) 2026-03-16 | 14/700 |
|
|---|---|---|---|---|
|
[考研] 304求調(diào)劑 +9 | 小熊joy 2026-03-14 | 9/450 |
|
|
[考研] 301求調(diào)劑 +4 | A_JiXing 2026-03-16 | 4/200 |
|
|
[考研] 344求調(diào)劑 +4 | knight344 2026-03-16 | 4/200 |
|
|
[碩博家園] 湖北工業(yè)大學(xué) 生命科學(xué)與健康學(xué)院-課題組招收2026級(jí)食品/生物方向碩士 +3 | 1喜春8 2026-03-17 | 5/250 |
|
|
[考研] 考研調(diào)劑 +3 | 淇ya_~ 2026-03-17 | 5/250 |
|
|
[考研] 286求調(diào)劑 +3 | lemonzzn 2026-03-16 | 5/250 |
|
|
[考研] 0854控制工程 359求調(diào)劑 可跨專業(yè) +3 | 626776879 2026-03-14 | 9/450 |
|
|
[考研] 0703化學(xué)調(diào)劑,求各位老師收留 +8 | 秋有木北 2026-03-14 | 8/400 |
|
|
[考博] 東華理工大學(xué)化材專業(yè)26屆碩士博士申請(qǐng) +6 | zlingli 2026-03-13 | 6/300 |
|
|
[考研] 中科大材料專碩319求調(diào)劑 +3 | 孟鑫材料 2026-03-13 | 3/150 |
|
|
[考研] 0703,333分求調(diào)劑 一志愿鄭州大學(xué)-物理化學(xué) +3 | 李魔女斗篷 2026-03-11 | 3/150 |
|
|
[考研] 336求調(diào)劑 +6 | Iuruoh 2026-03-11 | 6/300 |
|
|
[考研] 材料工程調(diào)劑 +9 | 咪咪空空 2026-03-12 | 9/450 |
|
|
[考研] 0856材料與化工301求調(diào)劑 +5 | 奕束光 2026-03-13 | 5/250 |
|
|
[考研] 333求調(diào)劑 +3 | 球球古力 2026-03-11 | 3/150 |
|
|
[考研] 329求調(diào)劑 +3 | miaodesi 2026-03-12 | 4/200 |
|
|
[考研] 310求調(diào)劑 +3 | 【上上簽】 2026-03-11 | 3/150 |
|
|
[考研] 274求調(diào)劑 +3 | S.H1 2026-03-12 | 3/150 |
|
|
[考研] 289求調(diào)劑 +3 | 李政瑩 2026-03-12 | 3/150 |
|