 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 2513 | 回復(fù): 4 | ||||
| 當(dāng)前只顯示滿(mǎn)足指定條件的回帖，點(diǎn)擊這里查看本話(huà)題的所有回帖 | ||||
Cythia_cd4新蟲(chóng) (初入文壇)
|
[求助]
第一性原理計(jì)算 EOS擬合優(yōu)化晶胞 已有1人參與
|
|||
|
對(duì)于EOS擬合優(yōu)化晶胞一直存在著很多疑問(wèn)： 1. 大師兄的教程是計(jì)算不同體積晶胞的單點(diǎn)能(ISIF=2, NSW=0)，將獲得的能量和體積進(jìn)行BM方程擬合來(lái)獲得平衡晶格常數(shù)； 2. 侯柱峰老師的方法是對(duì)不同體積的晶胞先進(jìn)行一次離子弛豫（ISIF=2, 4或5），之后再進(jìn)行靜態(tài)計(jì)算；最后靜態(tài)計(jì)算得到的能量再和體積進(jìn)行BM方程擬合； ************* 這兩種方法的區(qū)別主要在于原子弛豫的順序，那在計(jì)算當(dāng)中我們應(yīng)該選擇哪種方法呢？ |
回復(fù)此樓» 猜你喜歡
307求調(diào)劑
已經(jīng)有4人回復(fù)
-
收調(diào)劑
已經(jīng)有3人回復(fù)
-
315食品工程求調(diào)劑
已經(jīng)有6人回復(fù)
-
研究生招生
已經(jīng)有6人回復(fù)
-
材料與化工，291，求調(diào)劑
已經(jīng)有12人回復(fù)
-
0817一志愿蘇州大學(xué)280復(fù)試調(diào)劑
已經(jīng)有8人回復(fù)
-
求調(diào)劑 一志愿蘇州大學(xué)，0856化工323分 | 本科應(yīng)化 | 有專(zhuān)利/競(jìng)賽/科研助手經(jīng)歷 |
已經(jīng)有9人回復(fù)
-
085602化工求調(diào)劑
已經(jīng)有4人回復(fù)
-
083000環(huán)境科學(xué)與工程調(diào)劑
已經(jīng)有6人回復(fù)
-
求材料讀博院校
已經(jīng)有9人回復(fù)
waveflying
新蟲(chóng) (著名寫(xiě)手)
- 應(yīng)助: 12 (小學(xué)生)
- 金幣: 182
- 散金: 3611
- 紅花: 23
- 沙發(fā): 51
- 帖子: 2945
- 在線(xiàn): 1696.5小時(shí)
- 蟲(chóng)號(hào): 1965238
- 注冊(cè): 2012-08-31
- 性別: GG
- 專(zhuān)業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)
【答案】應(yīng)助回帖
感謝參與，應(yīng)助指數(shù) +1
|
Pseudo-code For (set of Volumes: equilibrium volume ±5%){ Step 1 : Fixed Volume relaxation (IBRION = 2, ISIF=4, ENCUT = 1.3x ENMAX, LCHARG=.TRUE., NSW=100) Step 2→n-1: Second and following fixed Volume relaxation (until a threshold is crossed and the structure is relaxed in fewer than N ionic steps) (IBRION = 2, ISIF=4, ENCUT = 1.3x ENMAX, ICHARG=1, LCHARG=.TRUE., NSW=100) Step n : Static calculation (IBRION = -1, no ISIF parameter, ICHARG=1, ENCUT = 1.3x ENMAX, ICHARG=1, LCHARG=.TRUE., NSW=0) } Fit Volume-Energy to Equation of State. Fixed volume relaxation at equilibrium volume. (With continuations if too many ionic steps are required.) Static calculation at equilibrium volume |
5樓2019-07-24 16:22:32
卡開(kāi)發(fā)發(fā)
專(zhuān)家顧問(wèn) (著名寫(xiě)手)
Ab Initio Amateur
-
專(zhuān)家經(jīng)驗(yàn): +224 - 1ST強(qiáng)帖: 2
- 應(yīng)助: 47 (小學(xué)生)
- 金幣: 3505.8
- 散金: 14
- 紅花: 223
- 帖子: 2477
- 在線(xiàn): 1254小時(shí)
- 蟲(chóng)號(hào): 1369423
- 注冊(cè): 2011-08-16
- 性別: GG
- 專(zhuān)業(yè): 金屬材料的磨損與磨蝕
- 管轄: 第一性原理
|
原理上應(yīng)該這么考慮，如果晶體處于穩(wěn)定結(jié)構(gòu)： <1>能量E對(duì)每個(gè)原子的核坐標(biāo)xN有(∂E/∂xN)=0； <2>stress=0； <3>如果僅考慮<2>取trace的情況，能得到體積Ω滿(mǎn)足p=-(∂E/∂Ω)=0； 討論： 1、NSW=0，EOS找極小的點(diǎn)只是能夠滿(mǎn)足條件<3>，無(wú)法滿(mǎn)足條件<1>、<2>。唯一成立的條件是通過(guò)對(duì)稱(chēng)性能把晶格參數(shù)和3N個(gè)原子坐標(biāo)降低到只有一個(gè)獨(dú)立變量，比如Fe這種BCC或者Cu這樣的FCC的簡(jiǎn)單格子，約束對(duì)稱(chēng)性的情況只有晶格長(zhǎng)度a可以自由變化。 2、isif=2，NSW>0，只能滿(mǎn)足條件<1>、<3>，此時(shí)對(duì)于三個(gè)長(zhǎng)度參數(shù)任意或三個(gè)夾角任意的情況就不是很合理。成立的條件是能將晶格參數(shù)約化到一個(gè)獨(dú)立變量，比如一般的立方晶格應(yīng)該都行。 3、isif=4，我認(rèn)為還算原理上合理。 如果只要最優(yōu)晶格的情況，我覺(jué)得這些擬合方法也不算是很好的方法。比如我們?nèi)?yōu)化H2分子來(lái)得到平衡鍵長(zhǎng)，你可以計(jì)算不同鍵長(zhǎng)的H2然后擬合勢(shì)能曲線(xiàn)來(lái)得到極小點(diǎn)，也可以通過(guò)給定某個(gè)結(jié)構(gòu)的能量和力來(lái)通過(guò)優(yōu)化算法確定極小點(diǎn)，通常后者肯定效率要更高一些。對(duì)于尋找最優(yōu)晶格的情況，擬合EOS就是做E-V曲線(xiàn)找極小，但直接擬合需要計(jì)算很多點(diǎn)才能完成；這個(gè)問(wèn)題本身就是個(gè)一維極小問(wèn)題，通過(guò)一維線(xiàn)搜索算法可以更高效進(jìn)行解決，而且不需要假定晶格的E-V曲線(xiàn)長(zhǎng)什么樣。如果晶格不是很復(fù)雜，直接用isif=3穩(wěn)定的情況下應(yīng)該更合適一些。 |
2樓2019-07-23 22:53:23
Cythia_cd4
新蟲(chóng) (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 帖子: 16
- 在線(xiàn): 17.8小時(shí)
- 蟲(chóng)號(hào): 15043045
- 注冊(cè): 2019-05-07
- 性別: MM
|
非常感謝從原理出發(fā)進(jìn)行了如此詳細(xì)的解釋?zhuān)?br /> 不知道自己是否理解正確： 1. 對(duì)于ISIF的取值需要我們根據(jù)我們體系而定：對(duì)于只有一個(gè)參數(shù)（比如BCC,FCC簡(jiǎn)單格子只需要優(yōu)化晶格長(zhǎng)度a）我們可以用ISIF=2；對(duì)于有多個(gè)參數(shù)需要優(yōu)化，ISIF=4或者5更為合理； 2. 對(duì)于NSW的取值還是有些許疑問(wèn)：我們是否可以在優(yōu)化晶格參數(shù)的時(shí)候先用NSW=0來(lái)固定原子位置；在用BM方程確定晶格參數(shù)后再進(jìn)行一次離子弛豫； ***** 覺(jué)得自己對(duì)第一性原理計(jì)算理解得太淺顯了，向您學(xué)習(xí)！ |
3樓2019-07-24 08:48:49
卡開(kāi)發(fā)發(fā)
專(zhuān)家顧問(wèn) (著名寫(xiě)手)
Ab Initio Amateur
-
專(zhuān)家經(jīng)驗(yàn): +224 - 1ST強(qiáng)帖: 2
- 應(yīng)助: 47 (小學(xué)生)
- 金幣: 3505.8
- 散金: 14
- 紅花: 223
- 帖子: 2477
- 在線(xiàn): 1254小時(shí)
- 蟲(chóng)號(hào): 1369423
- 注冊(cè): 2011-08-16
- 性別: GG
- 專(zhuān)業(yè): 金屬材料的磨損與磨蝕
- 管轄: 第一性原理
|
1、原理上一般復(fù)雜的情況其實(shí)isif=3最方便，如果把H2分子的掃描和優(yōu)化作為例子來(lái)看，isif=3相當(dāng)于是那個(gè)優(yōu)化，isif=4+EOS就相當(dāng)于掃描；同樣，做isif=5和做NSW=0的掃描也是這樣的關(guān)系，但兩者均缺乏原子弛豫，而isif=2與這兩個(gè)正好相反，僅考慮了原子弛豫而不考慮晶格形狀的變化。 2、NSW先優(yōu)化晶格，可能得到的晶格對(duì)于當(dāng)前離子分?jǐn)?shù)坐標(biāo)固定的情況是最優(yōu)，在進(jìn)行離子弛豫之后晶格就又不是最優(yōu)的了，要找的點(diǎn)其實(shí)應(yīng)該是兩者同時(shí)都是最優(yōu)。對(duì)于isif=4掃描體積的情況其實(shí)是一直約束原子坐標(biāo)和晶格形狀是最優(yōu)的，所以才比較合理。 |
4樓2019-07-24 15:11:06
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 288求調(diào)劑 +10 | 王曉陽(yáng)- 2026-03-09 | 15/750 |
|
|---|---|---|---|---|
|
|
[考研] 265求調(diào)劑 +4 | 小木蟲(chóng)085600 2026-03-09 | 7/350 |
|
|
|
[考研] 材料與化工求調(diào)劑 +10 | 與冬清寧 2026-03-07 | 13/650 |
|
|
[考研] 家人們 調(diào)劑不迷路 看這里 +8 | likeihood 2026-03-09 | 13/650 |
|
|
|
[考研] 歡迎211本科同學(xué)，過(guò)A區(qū)國(guó)家線(xiàn)，A區(qū)非偏遠(yuǎn)一本，交叉學(xué)科課題組 +22 | lisimayy 2026-03-04 | 35/1750 |
|
|
|
[考研] 一志愿天津大學(xué)，英一數(shù)二305分求調(diào)劑，四六級(jí)已過(guò) +5 | 小小番的茄 2026-03-09 | 5/250 |
|
|
|
[考研] 0703化學(xué) 學(xué)碩 理工科均可 不區(qū)分研究方向 總分279求調(diào)劑 +9 | 1一11 2026-03-03 | 9/450 |
|
|
|
[考研] 求調(diào)劑，一志愿江南大學(xué)，食品科學(xué)與工程，總分，320 +3 | yyyyyukino 2026-03-07 | 3/150 |
|
|
|
[考研] 招調(diào)劑 +20 | 帆船哥 2026-03-04 | 23/1150 |
|
|
|
[考研] 0817化學(xué)工程與技術(shù)312分求調(diào)劑 +7 | T123 tt 2026-03-04 | 7/350 |
|
|
|
[考研] 材料調(diào)劑 +7 | ounce. 2026-03-05 | 15/750 |
|
|
|
[考研] 08工科 +4 | li李樂(lè)成 2026-03-06 | 4/200 |
|
|
|
[考研] 304求調(diào)劑 +4 | 52hz~~ 2026-03-05 | 5/250 |
|
|
|
[考研] 材料調(diào)劑 +10 | ounce. 2026-03-04 | 12/600 |
|
|
|
[考研] 0856材料與化工求調(diào)劑！ +5 | 化工考生111 2026-03-04 | 11/550 |
|
|
|
[考研] 332材料求調(diào)劑 +6 | zjy101327 2026-03-05 | 7/350 |
|
|
|
[考研] 085602 293分求調(diào)劑 +3 | SivanNano. 2026-03-05 | 3/150 |
|
|
|
[考研] 一志愿西交化工專(zhuān)碩288專(zhuān)業(yè)課93求調(diào)劑 +6 | 好運(yùn)好運(yùn)接接 2026-03-04 | 6/300 |
|
|
|
[考研] 293求調(diào)劑 +3 | 是樂(lè)渝哇 2026-03-04 | 3/150 |
|
|
|
[考研] 085601 材料305分求助 +4 | 泡泡郵件 2026-03-03 | 6/300 |
|