| 8 | 1/1 | 返回列表 |
| 查看: 5256 | 回復(fù): 7 | |||
Prettyswan木蟲 (小有名氣)
|
[交流]
【求助】結(jié)構(gòu)優(yōu)化中,如何選擇POTIM? 已有6人參與
|
|
說明手冊上說: 這里推薦在結(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手冊,雖然講的很詳細,當(dāng)有個地方還不是很明白。就他例子的關(guān)于Mg的晶格參數(shù)的優(yōu)化(P25)。 我取10個點,循環(huán)完后,擬合得到最佳的a值。他后面說“在該體積下,重復(fù)上面第2和3步,得到平衡狀態(tài)下的其他晶胞參數(shù)。這一步就是:在得到了E(V)曲線后,通過擬合狀態(tài)方程得到平衡狀態(tài)下的體積,計算出上面腳本中的變量$i的值,并改變$i的循環(huán)值,再運行run_cell計算一次,得到其他的結(jié)構(gòu)參數(shù)c和位置u.! 對這句話我不太理解,安我自己的意思,通過擬合上面E(V)曲線得到平衡狀態(tài)下的體積,計算出變量$i,然后在該值下遲豫,就可以得到b,c。而不是候博士上面說的“改變$i的循環(huán)值,再運行run_cell計算一次”。 不知道大家怎么理解他這句話的。 |
至尊木蟲 (職業(yè)作家)
|
“改變$i的循環(huán)值,再運行run_cell計算一次” 此時$i只附了一個值,就是前面由EOS擬合得到的平衡體積V_0計算(V_0 = 0.5*sqrt(3.0) a^3 * c/a, c/a的值取自run_cell所設(shè)置的,即c/a的值與前一步中run_cell中的相同。這樣得到了一個平衡體積下初始的a, c/a)。這一次run_cell就是固定體積(將體積固定在所擬合出來的V_0),根據(jù)run_cell所設(shè)置的初始a, c/a進行固定體積的優(yōu)化。 E(a,c) -----> E(V, a)------->E(V_0,a)------> V_0下最優(yōu)的a, c的值。 |
|
謝謝斑竹的回答,我理解了這個意思。當(dāng)我還有點疑問,vasp 好像是在開始十步優(yōu)化的時候,都是固定角度和c/a的值進行優(yōu)化的。萬一我優(yōu)化的是一個三斜晶系,我該怎么辦。我還是安你的這個步驟:E(a,b,c)->E(v,a)->E(V_0,a)->V_0下的最佳a,b,c,alpha,beta,gama嗎?不知道這個有通用步驟沒有? 因為我還想優(yōu)化結(jié)構(gòu):a。絙。絚,三個角度都為90度。如果還是按Mg那種方法的話,b/a的值將會始終不變。 [ Last edited by swx0789 on 2010-12-25 at 15:22 ] |
銀蟲 (正式寫手)
| 8 | 1/1 | 返回列表 |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
|
百度網(wǎng)盤 |
360云盤 |
千易網(wǎng)盤 |
華為網(wǎng)盤
在新窗口頁面中打開自己喜歡的網(wǎng)盤網(wǎng)站,將文件上傳后,然后將下載鏈接復(fù)制到帖子內(nèi)容中就可以了。 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 070300化學(xué)319求調(diào)劑 +4 | 錦鯉0909 2026-03-17 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 344求調(diào)劑 +4 | knight344 2026-03-16 | 4/200 |
|
|
[考研] 308求調(diào)劑 +4 | 是Lupa啊 2026-03-16 | 4/200 |
|
|
[考研] 293求調(diào)劑 +6 | 世界首富 2026-03-11 | 6/300 |
|
|
[考研] 290求調(diào)劑 +3 | p asserby. 2026-03-15 | 4/200 |
|
|
[考研] 求調(diào)劑,總分315,考的生物醫(yī)藥,一志愿湖南師范大學(xué)。調(diào)劑到任何專業(yè)都可以 +4 | 小丁想進步 2026-03-11 | 5/250 |
|
|
[考博] 26申博 +4 | 八旬速覽 2026-03-16 | 4/200 |
|
|
[考研] 275求調(diào)劑 +4 | 太陽花天天開心 2026-03-16 | 4/200 |
|
|
[考研] 東南大學(xué)364求調(diào)劑 +5 | JasonYuiui 2026-03-15 | 5/250 |
|
|
[考研] 中科院材料273求調(diào)劑 +4 | yzydy 2026-03-15 | 4/200 |
|
|
[考研] 26考研一志愿中國石油大學(xué)(華東)305分求調(diào)劑 +3 | 嘉年新程 2026-03-15 | 3/150 |
|
|
[考研] 0856專碩279求調(diào)劑 +5 | 加油加油!? 2026-03-15 | 5/250 |
|
|
[考研] 085601材料工程315分求調(diào)劑 +3 | yang_0104 2026-03-15 | 3/150 |
|
|
[考研] 289求調(diào)劑 +4 | 這么名字咋樣 2026-03-14 | 6/300 |
|
|
[考研] 26調(diào)劑/材料/英一數(shù)二/總分289/已過A區(qū)線 +6 | 步川酷紫123 2026-03-13 | 6/300 |
|
|
[考研] 333求調(diào)劑 +3 | 球球古力 2026-03-11 | 3/150 |
|
|
[考研] 考研調(diào)劑 +4 | 芬達46 2026-03-12 | 4/200 |
|
|
[考研] 307求調(diào)劑 +5 | 超級伊昂大王 2026-03-12 | 5/250 |
|
|
[考研] 0817化學(xué)工程與技術(shù)考研312分調(diào)劑 +3 | T123 tt 2026-03-12 | 3/150 |
|
|
[考研] 081200-11408-276學(xué)碩求調(diào)劑 +3 | 崔wj 2026-03-12 | 4/200 |
|