 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 3318 | 回復(fù): 9 | ||||||||
| 當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖 | ||||||||
奧古拉斯無蟲 (小有名氣)
|
[交流]
【原創(chuàng)】晶格優(yōu)化中ISIF取值的探究（六方MgO為例） 已有6人參與
|
|||||||
|
一、目的：搞清楚用vasp進(jìn)行晶格優(yōu)化中ISIF值的重要性 二、方法：采用分步優(yōu)化和一次性優(yōu)化兩種方法探究 三、重點(diǎn)參數(shù)說明： ISIF：2－>relax ions?/Y; || change cell shape?/N; || change cell volume?/N 3－>relax ions?/Y; || change cell shape?/Y; || change cell volume?/Y 4－>relax ions?/Y; || change cell shape?/Y; || change cell volume?/N 5－>relax ions?/N; || change cell shape?/Y; || change cell volume?/N 6－>relax ions?/N; || change cell shape?/Y; || change cell volume?/Y 7－>relax ions?/N; || change cell shape?/N; || change cell volume?/Y POSCAR： a 0.0 -1.0 0.0 0.8660254037844 0.5 0.0 0.0 0.0 1.6230529595 2 Direct 0.6666666666666667 0.3333333333333333 0.750 0.3333333333333333 0.6666666666666667 0.250 四、實(shí)驗(yàn)步驟： A、采用分步優(yōu)化法： 算法原理： a.先進(jìn)行一次馳豫優(yōu)化計(jì)算 b.將CONTCAR拷貝成POSCAR c.再進(jìn)行一次靜態(tài)計(jì)算，求出總能和體積，擬合狀態(tài)方程，得出晶格常 數(shù) I、取ISIF=2時(shí)：（＃上下對應(yīng)。上行為系列縮放系數(shù)；下行為相應(yīng)的總能；所取CONTAR＆volume為縮放系數(shù)取3.71時(shí)的結(jié)構(gòu)情況） 2.91 3.01 3.11 3.21 3.31 3.41 3.51 3.61 3.71 -0.09768706 -2.60381322 -2.62609278 -2.62612787 -2.62612791 -2.60064488 -2.59725747 -2.59725252 -2.597253 volume of cell : 71.78 CONTCAR： 3.710000000000000 0.0000000000000000 -1.0000000000000000 0.0000000000000000 0.8660254037844000 0.5000000000000000 0.0000000000000000 0.0000000000000000 0.0000000000000000 1.6230529595000000 II、取ISIF＝3時(shí)： 2.91 3.01 3.11 3.21 3.31 3.41 3.51 3.61 3.71 -1.47880875 -3.56749270 -3.57829654 -3.57830700 -3.57830701 -3.55534231 -3.55344206 -3.55344601 -3.553446 volume of cell : 43.12 CONTCAR： 3.710000000000000 0.0000000000000000 -0.8480309609621116 0.0000000000000000 0.7344163553888973 0.4240154804810558 0.0000000000000000 0.0000000000000000 0.0000000000000000 1.3556928126902820 III、取ISIF＝4時(shí)： 2.91 3.01 3.11 3.21 3.31 3.41 3.51 3.61 3.71 -0.24918897 -2.60258815 -2.63108804 -2.63118786 -2.63118806 -2.60762910 -2.60458644 -2.60459724 -2.604597 volume of cell : 71.78 CONTCAR： 3.710000000000000 0.0000000000000001 -1.0553127846071310 0.0000000000000000 0.9139276804085778 0.5276563923035652 0.0000000000000000 0.0000000000000001 0.0000000000000000 1.4573715661012910 IV、取ISIF＝5時(shí)： 2.91 3.01 3.11 3.21 3.31 3.41 3.51 3.61 3.71 -0.24918897 -2.60258815 -2.63108804 -2.63118786 -2.63118806 -2.60762910 -2.60458644 -2.60459724 -2.604597 volume of cell : 71.78 CONTCAR： 3.710000000000000 0.0000000000000001 -1.0553127846071310 0.0000000000000000 0.9139276804085778 0.5276563923035652 0.0000000000000000 0.0000000000000001 0.0000000000000000 1.4573715661012910 V、取ISIF＝6時(shí)： 2.91 3.01 3.11 3.21 3.31 3.41 3.51 3.61 3.71 -1.47880875 -3.56749270 -3.57829654 -3.57830700 -3.57830701 -3.55534231 -3.55344206 -3.55344601 -3.553446 volume of cell : 43.12 CONTCAR: 3.710000000000000 0.0000000000000000 -0.8480309609621116 0.0000000000000000 0.7344163553888973 0.4240154804810558 0.0000000000000000 0.0000000000000000 0.0000000000000000 1.3556928126902820 VI、取ISIF＝7時(shí)： 2.91 3.01 3.11 3.21 3.31 3.41 3.51 3.61 3.71 -1.28296288 -3.56234825 -3.57418319 -3.57419641 -3.57419642 -3.55112710 -3.54922372 -3.54922856 -3.549229 volume of cell : 43.07 CONTCAR： 3.710000000000000 0.0000000000000000 -0.8434465875014347 0.0000000000000000 0.7304461715115043 0.4217232937507173 0.0000000000000000 0.0000000000000000 0.0000000000000000 1.3689584800243790 B：采用一次優(yōu)化法： 與上面方法不同的是INCAR中取ISIF＝3，取定的實(shí)驗(yàn)縮放系數(shù)a＝3.21，計(jì)算的精度選得更高，且僅做一次優(yōu)化。 結(jié)果為： TOTEN = -3.549568 eV volume of cell : 43.12 CONTCAR： Mg-hex 3.210000000000000 0.0000000000000000 -0.9801232359458133 0.0000000000000000 0.8488116211684474 0.4900616179729067 0.0000000000000000 0.0000000000000000 0.0000000000000000 1.5669819270266370 2 Direct 0.6666666666666643 0.3333333333333357 0.7500000000000000 0.3333333333333357 0.6666666666666643 0.2500000000000000 結(jié)論： 1、特殊位置的粒子隨晶格形狀的改變而變化，不論允許馳豫與否。 2、不同ISIF的取值對晶格優(yōu)化影響不大，本例中縮放系數(shù)均在3.21~3.31之范圍內(nèi) 3、采用分步優(yōu)化法得出的優(yōu)化晶格比較可靠 |
回復(fù)此樓» 收錄本帖的淘帖專輯推薦
 VASP | 計(jì)算-vasp | VASP相關(guān)主題 | 科研 |
| vasp |
» 猜你喜歡
291求調(diào)劑
已經(jīng)有4人回復(fù)
-
生物學(xué)一志愿985，分?jǐn)?shù)349求調(diào)劑
已經(jīng)有8人回復(fù)
-
328求調(diào)劑，英語六級551，有科研經(jīng)歷
已經(jīng)有12人回復(fù)
-
070300化學(xué)求調(diào)劑
已經(jīng)有6人回復(fù)
-
考研化學(xué)308分求調(diào)劑
已經(jīng)有7人回復(fù)
-
一志愿上海交大生物與醫(yī)藥專碩324分，求調(diào)劑
已經(jīng)有5人回復(fù)
-
工科0856求調(diào)劑
已經(jīng)有5人回復(fù)
-
一志愿中國石油大學(xué)（華東） 本科齊魯工業(yè)大學(xué)
已經(jīng)有4人回復(fù)
-
279分求調(diào)劑 一志愿211
已經(jīng)有16人回復(fù)
-
求調(diào)劑一志愿武漢理工大學(xué)材料工程（085601）
已經(jīng)有5人回復(fù)
» 本主題相關(guān)價(jià)值貼推薦，對您同樣有幫助:
- 求助KH晶格參數(shù)，急!
已經(jīng)有3人回復(fù)
- 能帶圖圖中k點(diǎn)取值
已經(jīng)有6人回復(fù)
- 晶界處的氧氣和晶格中的氧空位如何表征
已經(jīng)有8人回復(fù)
- 在做CASTP模擬時(shí)，建立晶格需要晶格常數(shù)，如何獲取呢？
已經(jīng)有10人回復(fù)
- 用LDA贗勢優(yōu)化出來的晶格常數(shù)一定比實(shí)驗(yàn)值小嗎？謝謝！
已經(jīng)有10人回復(fù)
- 求問用X-ray diffraction 能定量的測定出晶體的晶格參數(shù)嗎？
已經(jīng)有5人回復(fù)
- 晶格常數(shù)優(yōu)化
已經(jīng)有7人回復(fù)
- 超晶格和量子阱的區(qū)別？
已經(jīng)有7人回復(fù)
- 【求助】如何證明摻雜離子進(jìn)入晶格
已經(jīng)有8人回復(fù)
- 【求助】vasp優(yōu)化時(shí)選擇了isif=3 怎么限制晶格的一個(gè)維度呢？
已經(jīng)有16人回復(fù)
- 【求助】蟲友們對ISIF取值有什么高見呢
已經(jīng)有10人回復(fù)
- 【求助成功】晶格優(yōu)化中，murn.f的使用方法
已經(jīng)有7人回復(fù)
- 【求助】si晶格常數(shù)
已經(jīng)有7人回復(fù)
奧古拉斯
無蟲 (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: -8.7
- 散金: 4
- 帖子: 93
- 在線: 25.9小時(shí)
- 蟲號: 714435
- 注冊: 2009-03-04
- 專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)
★ ★
xiaohunhun(金幣+2):謝謝 2010-06-08 23:39:26
xiaohunhun(金幣+2):謝謝 2010-06-08 23:39:26
|
謝謝鼓勵(lì)  1. Q: “II、取ISIF＝3”，和“B：采用一次優(yōu)化法” 都是ISIF=3.這兩個(gè)有什么區(qū)別呢？ A:“II、取ISIF＝3”，這樣取定ISIF后，會(huì)對所有的縮放系數(shù)試探值下的結(jié)構(gòu)全優(yōu)化，也就是在每次優(yōu)化過程中離子可以馳豫，晶格形狀允許改變。這樣就得到了一系列的V～E數(shù)對，然后進(jìn)行狀態(tài)方程的擬合，從而得到最穩(wěn)定的結(jié)構(gòu)情況；而“B：采用一次優(yōu)化法”，僅僅是取定了實(shí)驗(yàn)上得到的縮放系數(shù)，在這個(gè)基礎(chǔ)上進(jìn)行結(jié)構(gòu)的全優(yōu)化，而且僅作這一次運(yùn)算，所得到的穩(wěn)定結(jié)構(gòu)情況在輸出文件CONTCAR中。 2.Q:，為什么有的volume of cell 是七十多，有的四十多呢? A:看I和II這兩種情況，I是ISIF＝2，II是ISIF＝3，然后你看再前面的ISIF參數(shù)說明，就會(huì)知道ISIF＝2在優(yōu)化過程中要保持c/a，以及volume不變，而ISIF＝3則相反；同時(shí)我已經(jīng)說明了我列出的volume of cell都是在縮放系數(shù)為3.71時(shí)的情況。據(jù)此可以這樣看，ISIF＝2求出的結(jié)構(gòu)是ISIF＝3的初始結(jié)構(gòu)，經(jīng)過可變晶格的優(yōu)化計(jì)算，當(dāng)然會(huì)得到與初始結(jié)構(gòu)完全不同的結(jié)構(gòu)。 呵呵，希望我解釋清楚了 |
3樓2010-06-08 17:16:19
xirainbow
木蟲 (正式寫手)
- 1ST強(qiáng)帖: 2
- 應(yīng)助: 0 (幼兒園)
- 貴賓: 0.03
- 金幣: 2132.6
- 紅花: 2
- 帖子: 775
- 在線: 6.1小時(shí)
- 蟲號: 431211
- 注冊: 2007-08-11
- 專業(yè): 光譜學(xué)
2樓2010-06-08 14:47:21
xirainbow
木蟲 (正式寫手)
- 1ST強(qiáng)帖: 2
- 應(yīng)助: 0 (幼兒園)
- 貴賓: 0.03
- 金幣: 2132.6
- 紅花: 2
- 帖子: 775
- 在線: 6.1小時(shí)
- 蟲號: 431211
- 注冊: 2007-08-11
- 專業(yè): 光譜學(xué)
★ ★ ★
小木蟲(金幣+0.5):給個(gè)紅包，謝謝回帖交流
xiaohunhun(金幣+2):謝謝 2010-06-08 23:39:36
小木蟲(金幣+0.5):給個(gè)紅包，謝謝回帖交流
xiaohunhun(金幣+2):謝謝 2010-06-08 23:39:36
|
非常非常感謝您給出如此詳細(xì)的說明；） 但是對于第一個(gè)小問題，我依然有些疑問：） 既然ISIF都是=3，那么vasp將對所有的原子位置和晶包參數(shù)進(jìn)行優(yōu)化。 那么A和B中的ISIF=3的時(shí)，在相同的輸入結(jié)構(gòu)的前提下，他們的計(jì)算過程和結(jié)果有差別么？ 另：2.91 3.01 3.11 3.21 3.31 3.41 3.51 3.61 3.71 以上的這幾個(gè)參數(shù)，是您人工指定的，還是vasp在優(yōu)化時(shí)，自己產(chǎn)生的中間結(jié)構(gòu)？ 如果是vasp自己產(chǎn)生的，您如何保證A中的ISIF=不同數(shù)值的時(shí)候，都取這幾個(gè)參數(shù)？ 再：將V～E進(jìn)行狀態(tài)方程的擬合。vasp有提供這個(gè)功能么？還請指教；） [ Last edited by xirainbow on 2010-6-8 at 20:36 ] |
4樓2010-06-08 20:27:41
奧古拉斯
無蟲 (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: -8.7
- 散金: 4
- 帖子: 93
- 在線: 25.9小時(shí)
- 蟲號: 714435
- 注冊: 2009-03-04
- 專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)
5樓2010-06-08 21:48:24
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 279分求調(diào)劑 一志愿211 +14 | chaojifeixia 2026-03-19 | 16/800 |
|
|---|---|---|---|---|
|
|
[考研] 316求調(diào)劑 +7 | 梁茜雯 2026-03-19 | 7/350 |
|
|
[考研]
 0854 考研調(diào)劑 招生了！AI 方向
+3
|
pk3725069 2026-03-19 | 14/700 |
|
|
|
[考研] 336化工調(diào)劑 +3 | 王大坦1 2026-03-23 | 4/200 |
|
|
|
[考研] 317求調(diào)劑 +12 | 申子申申 2026-03-19 | 18/900 |
|
|
|
[考研] 307求調(diào)劑 +11 | 冷笙123 2026-03-17 | 11/550 |
|
|
|
[考研] 308求調(diào)劑 +3 | 墨墨漠 2026-03-21 | 3/150 |
|
|
|
[考研] 275求調(diào)劑 +6 | shansx 2026-03-22 | 8/400 |
|
|
|
[考研] 考研調(diào)劑 +4 | 來好運(yùn)來來來 2026-03-21 | 4/200 |
|
|
|
[基金申請] 山東省面上項(xiàng)目限額評審 +4 | 石瑞0426 2026-03-19 | 4/200 |
|
|
|
[考研] 初試 317 +7 | 半拉月丙 2026-03-20 | 7/350 |
|
|
|
[考研] 一志愿南大，0703化學(xué)，分?jǐn)?shù)336，求調(diào)劑 +3 | 收到VS 2026-03-21 | 3/150 |
|
|
|
[考研] 336求調(diào)劑 +5 | rmc8866 2026-03-21 | 5/250 |
|
|
|
[考研] 0805材料320求調(diào)劑 +3 | 深海物語 2026-03-20 | 3/150 |
|
|
|
[考研] 265求調(diào)劑 +9 | 梁梁校校 2026-03-17 | 9/450 |
|
|
|
[考研] 085700資源與環(huán)境308求調(diào)劑 +12 | 墨墨漠 2026-03-18 | 13/650 |
|
|
|
[考研] 材料 336 求調(diào)劑 +3 | An@. 2026-03-18 | 4/200 |
|
|
|
[考研] 材料學(xué)碩297已過四六級求調(diào)劑推薦 +11 | adaie 2026-03-19 | 11/550 |
|
|
|
[考研] 353求調(diào)劑 +3 | 拉鉤不許變 2026-03-20 | 3/150 |
|
|
|
[考研] 本科鄭州大學(xué)物理學(xué)院，一志愿華科070200學(xué)碩，346求調(diào)劑 +4 | 我不是一根蔥 2026-03-18 | 4/200 |
|