 |
 |
 |
| 24小時(shí)熱門版塊排行榜 |
- >論壇更新日志 (409)
- >考研 (29)
- >導(dǎo)師招生 (21)
- >考博 (13)
- >蟲(chóng)友互識(shí) (13)
- >無(wú)機(jī)/物化 (7)
- >文獻(xiàn)求助 (5)
- >綠色求助(高懸賞) (3)
- >碩博家園 (3)
- >教師之家 (2)
- >論文投稿 (2)
- >休閑灌水 (2)
- >留學(xué)生活 (1)
- >論文道賀祈福 (1)
- >基金申請(qǐng) (1)
- >檢索知識(shí) (1)
 |
 |
 |
||
 |
 |
 |
 |
 |
| 10 | 1/1 | 返回列表 |
| 查看: 3289 | 回復(fù): 9 | ||||||||
奧古拉斯無(wú)蟲(chóng) (小有名氣)
|
[交流]
【原創(chuàng)】晶格優(yōu)化中ISIF取值的探究（六方MgO為例） 已有6人參與
|
|
一、目的：搞清楚用vasp進(jìn)行晶格優(yōu)化中ISIF值的重要性 二、方法：采用分步優(yōu)化和一次性優(yōu)化兩種方法探究 三、重點(diǎn)參數(shù)說(shuō)明： 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í)：（＃上下對(duì)應(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的取值對(duì)晶格優(yōu)化影響不大，本例中縮放系數(shù)均在3.21~3.31之范圍內(nèi) 3、采用分步優(yōu)化法得出的優(yōu)化晶格比較可靠 |
回復(fù)此樓» 收錄本帖的淘帖專輯推薦
 VASP | 計(jì)算-vasp | VASP相關(guān)主題 | 科研 |
| vasp |
» 猜你喜歡
環(huán)境領(lǐng)域全國(guó)重點(diǎn)實(shí)驗(yàn)室招收博士1-2名
已經(jīng)有3人回復(fù)
-
材料專碩306英一數(shù)二
已經(jīng)有10人回復(fù)
-
301求調(diào)劑
已經(jīng)有6人回復(fù)
-
一志愿天津大學(xué)化學(xué)工藝專業(yè)（081702）315分求調(diào)劑
已經(jīng)有7人回復(fù)
-
302求調(diào)劑
已經(jīng)有6人回復(fù)
-
26博士申請(qǐng)
已經(jīng)有3人回復(fù)
-
268求調(diào)劑
已經(jīng)有3人回復(fù)
-
311求調(diào)劑
已經(jīng)有10人回復(fù)
-
被我言中：新模板不強(qiáng)調(diào)格式了，假專家開(kāi)始管格式了
已經(jīng)有4人回復(fù)
-
303求調(diào)劑
已經(jīng)有5人回復(fù)
» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:
- 求助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贗勢(shì)優(yōu)化出來(lái)的晶格常數(shù)一定比實(shí)驗(yàn)值小嗎？謝謝！
已經(jīng)有10人回復(fù)
- 求問(wèn)用X-ray diffraction 能定量的測(cè)定出晶體的晶格參數(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ù)
- 【求助】蟲(chóng)友們對(duì)ISIF取值有什么高見(jiàn)呢
已經(jīng)有10人回復(fù)
- 【求助成功】晶格優(yōu)化中，murn.f的使用方法
已經(jīng)有7人回復(fù)
- 【求助】si晶格常數(shù)
已經(jīng)有7人回復(fù)
xirainbow
木蟲(chóng) (正式寫手)
- 1ST強(qiáng)帖: 2
- 應(yīng)助: 0 (幼兒園)
- 貴賓: 0.03
- 金幣: 2132.6
- 紅花: 2
- 帖子: 775
- 在線: 6.1小時(shí)
- 蟲(chóng)號(hào): 431211
- 注冊(cè): 2007-08-11
- 專業(yè): 光譜學(xué)
2樓2010-06-08 14:47:21
奧古拉斯
無(wú)蟲(chóng) (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: -8.7
- 散金: 4
- 帖子: 93
- 在線: 25.9小時(shí)
- 蟲(chóng)號(hào): 714435
- 注冊(cè): 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ì)對(duì)所有的縮放系數(shù)試探值下的結(jié)構(gòu)全優(yōu)化，也就是在每次優(yōu)化過(guò)程中離子可以馳豫，晶格形狀允許改變。這樣就得到了一系列的V～E數(shù)對(duì)，然后進(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ù)說(shuō)明，就會(huì)知道ISIF＝2在優(yōu)化過(guò)程中要保持c/a，以及volume不變，而ISIF＝3則相反；同時(shí)我已經(jīng)說(shuō)明了我列出的volume of cell都是在縮放系數(shù)為3.71時(shí)的情況。據(jù)此可以這樣看，ISIF＝2求出的結(jié)構(gòu)是ISIF＝3的初始結(jié)構(gòu)，經(jīng)過(guò)可變晶格的優(yōu)化計(jì)算，當(dāng)然會(huì)得到與初始結(jié)構(gòu)完全不同的結(jié)構(gòu)。 呵呵，希望我解釋清楚了 |
3樓2010-06-08 17:16:19
xirainbow
木蟲(chóng) (正式寫手)
- 1ST強(qiáng)帖: 2
- 應(yīng)助: 0 (幼兒園)
- 貴賓: 0.03
- 金幣: 2132.6
- 紅花: 2
- 帖子: 775
- 在線: 6.1小時(shí)
- 蟲(chóng)號(hào): 431211
- 注冊(cè): 2007-08-11
- 專業(yè): 光譜學(xué)
★ ★ ★
小木蟲(chóng)(金幣+0.5):給個(gè)紅包，謝謝回帖交流
xiaohunhun(金幣+2):謝謝 2010-06-08 23:39:36
小木蟲(chóng)(金幣+0.5):給個(gè)紅包，謝謝回帖交流
xiaohunhun(金幣+2):謝謝 2010-06-08 23:39:36
|
非常非常感謝您給出如此詳細(xì)的說(shuō)明；） 但是對(duì)于第一個(gè)小問(wèn)題，我依然有些疑問(wèn)：） 既然ISIF都是=3，那么vasp將對(duì)所有的原子位置和晶包參數(shù)進(jìn)行優(yōu)化。 那么A和B中的ISIF=3的時(shí)，在相同的輸入結(jié)構(gòu)的前提下，他們的計(jì)算過(guò)程和結(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è)功能么？還請(qǐng)指教；） [ Last edited by xirainbow on 2010-6-8 at 20:36 ] |
4樓2010-06-08 20:27:41
奧古拉斯
無(wú)蟲(chóng) (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: -8.7
- 散金: 4
- 帖子: 93
- 在線: 25.9小時(shí)
- 蟲(chóng)號(hào): 714435
- 注冊(cè): 2009-03-04
- 專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)
★ ★ ★
xiaohunhun(金幣+3):謝謝 2010-06-08 23:38:41
xiaohunhun(金幣+3):謝謝 2010-06-08 23:38:41
|
呵呵，你的問(wèn)題提的很好！  限于時(shí)間，我寫的簡(jiǎn)略些，呵呵 1. 從過(guò)程上看更清楚：在馳豫階段，A和B只有精度上的差別，接下來(lái)，A還進(jìn)行了靜態(tài)計(jì)算，而B(niǎo)則沒(méi)有 2.這幾個(gè)值是在實(shí)驗(yàn)值基礎(chǔ)上左右等間隔人為取定的 3.擬合已經(jīng)有現(xiàn)成程序，不是vasp自己的，你搜搜'狀態(tài)方程的擬合'就有相關(guān)結(jié)果 |
5樓2010-06-08 21:48:24
xirainbow
木蟲(chóng) (正式寫手)
- 1ST強(qiáng)帖: 2
- 應(yīng)助: 0 (幼兒園)
- 貴賓: 0.03
- 金幣: 2132.6
- 紅花: 2
- 帖子: 775
- 在線: 6.1小時(shí)
- 蟲(chóng)號(hào): 431211
- 注冊(cè): 2007-08-11
- 專業(yè): 光譜學(xué)
6樓2010-06-09 09:28:15
zhoudq
新蟲(chóng) (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 69.5
- 帖子: 30
- 在線: 19.3小時(shí)
- 蟲(chóng)號(hào): 998260
- 注冊(cè): 2010-04-16
- 性別: GG
- 專業(yè): 極端條件下使用的金屬材料
7樓2010-06-09 10:20:46
awmc2008
至尊木蟲(chóng) (文壇精英)
學(xué)習(xí)者
- 應(yīng)助: 62 (初中生)
- 金幣: 13659.6
- 紅花: 10
- 帖子: 12257
- 在線: 736.8小時(shí)
- 蟲(chóng)號(hào): 600677
- 注冊(cè): 2008-09-12
- 性別: GG
- 專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)
8樓2010-06-09 16:33:53
ufo991
鐵蟲(chóng) (正式寫手)
- 1ST強(qiáng)帖: 1
- 應(yīng)助: 0 (幼兒園)
- 金幣: 1148.3
- 散金: 14
- 紅花: 1
- 帖子: 358
- 在線: 104.6小時(shí)
- 蟲(chóng)號(hào): 671678
- 注冊(cè): 2008-12-10
- 性別: MM
- 專業(yè): 凝聚態(tài)物性I:結(jié)構(gòu)、力學(xué)和
9樓2011-07-19 09:56:04
10樓2013-12-04 11:23:17
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 奧古拉斯 的主題更新
| 10 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 301求調(diào)劑 +4 | A_JiXing 2026-03-16 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 材料工程專碩調(diào)劑 +5 | 204818@lcx 2026-03-17 | 5/250 |
|
|
|
[考研] 302求調(diào)劑 +9 | 負(fù)心者當(dāng)誅 2026-03-11 | 9/450 |
|
|
|
[考研] 材料與化工專碩調(diào)劑 +5 | heming3743 2026-03-16 | 5/250 |
|
|
|
[考博] 26申博 +4 | 八旬速覽 2026-03-16 | 4/200 |
|
|
|
[考研] [導(dǎo)師推薦]西南科技大學(xué)國(guó)防/材料導(dǎo)師推薦 +3 | 尖角小荷 2026-03-16 | 6/300 |
|
|
|
[考研] 0856求調(diào)劑 +3 | 劉夢(mèng)微 2026-03-15 | 3/150 |
|
|
|
[考博] 東華理工大學(xué)化材專業(yè)26屆碩士博士申請(qǐng) +6 | zlingli 2026-03-13 | 6/300 |
|
|
|
[考研] 289求調(diào)劑 +4 | 這么名字咋樣 2026-03-14 | 6/300 |
|
|
|
[考研] 308 085701 四六級(jí)已過(guò)求調(diào)劑 +7 | 溫喬喬喬喬 2026-03-12 | 14/700 |
|
|
|
[考研] 337一志愿華南理工0805材料求調(diào)劑 +7 | mysdl 2026-03-11 | 9/450 |
|
|
|
[考研] 0703，333分求調(diào)劑 一志愿鄭州大學(xué)-物理化學(xué) +3 | 李魔女斗篷 2026-03-11 | 3/150 |
|
|
|
[考研] ［0860］321分求調(diào)劑，ab區(qū)皆可 +4 | 寶貴熱 2026-03-13 | 4/200 |
|
|
|
[考研] 材料工程調(diào)劑 +4 | 咪咪空空 2026-03-11 | 4/200 |
|
|
|
[碩博家園] 085600 260分求調(diào)劑 +3 | 天空還下雨么 2026-03-13 | 5/250 |
|
|
|
[考研] 材料專碩350 求調(diào)劑 +4 | 王金科 2026-03-12 | 4/200 |
|
|
|
[考研] 一志愿211化學(xué)學(xué)碩310分求調(diào)劑 +8 | 努力奮斗112 2026-03-12 | 9/450 |
|
|
|
[考研] 308求調(diào)劑 +3 | 是Lupa啊 2026-03-12 | 3/150 |
|
|
|
[考博] 福州大學(xué)楊黃浩課題組招收2026年專業(yè)學(xué)位博士研究生，2026.03.20截止 +3 | Xiangyu_ou 2026-03-12 | 3/150 |
|
|
|
[考研] 321求調(diào)劑（食品/專碩） +3 | xc321 2026-03-12 | 6/300 |
|