 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 2499 | 回復(fù): 4 | ||||
| 當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖 | ||||
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ù)此樓» 猜你喜歡
26申博自薦
已經(jīng)有7人回復(fù)
-
進(jìn)入個(gè)人成果庫(kù)好難，一下午都沒(méi)進(jìn)去
已經(jīng)有5人回復(fù)
-
接受26屆調(diào)劑生
已經(jīng)有12人回復(fù)
-
求調(diào)劑
已經(jīng)有5人回復(fù)
-
求調(diào)劑 一志愿蘇州大學(xué)，0856化工323分 | 本科應(yīng)化 | 有專利/競(jìng)賽/科研助手經(jīng)歷 |
已經(jīng)有6人回復(fù)
-
求調(diào)劑
已經(jīng)有4人回復(fù)
-
308求調(diào)劑
已經(jīng)有9人回復(fù)
-
303求調(diào)劑
已經(jīng)有7人回復(fù)
-
085600材料與化工 292分 求調(diào)劑
已經(jīng)有4人回復(fù)
-
一志愿中石油，化學(xué)工程，301求調(diào)劑
已經(jīng)有4人回復(fù)
Cythia_cd4
新蟲(chóng) (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 帖子: 16
- 在線: 17.8小時(shí)
- 蟲(chóng)號(hào): 15043045
- 注冊(cè): 2019-05-07
- 性別: MM
|
非常感謝從原理出發(fā)進(jìn)行了如此詳細(xì)的解釋！ 不知道自己是否理解正確： 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ā)
專家顧問(wèn) (著名寫(xiě)手)
Ab Initio Amateur
-
專家經(jīng)驗(yàn): +224 - 1ST強(qiáng)帖: 2
- 應(yīng)助: 47 (小學(xué)生)
- 金幣: 3503.8
- 散金: 14
- 紅花: 223
- 帖子: 2477
- 在線: 1254小時(shí)
- 蟲(chóng)號(hào): 1369423
- 注冊(cè): 2011-08-16
- 性別: GG
- 專業(yè): 金屬材料的磨損與磨蝕
- 管轄: 第一性原理
|
原理上應(yīng)該這么考慮，如果晶體處于穩(wěn)定結(jié)構(gòu)： <1>能量E對(duì)每個(gè)原子的核坐標(biāo)xN有(∂E/∂xN)=0； <2>stress=0； <3>如果僅考慮<2>取trace的情況，能得到體積Ω滿足p=-(∂E/∂Ω)=0； 討論： 1、NSW=0，EOS找極小的點(diǎn)只是能夠滿足條件<3>，無(wú)法滿足條件<1>、<2>。唯一成立的條件是通過(guò)對(duì)稱性能把晶格參數(shù)和3N個(gè)原子坐標(biāo)降低到只有一個(gè)獨(dú)立變量，比如Fe這種BCC或者Cu這樣的FCC的簡(jiǎn)單格子，約束對(duì)稱性的情況只有晶格長(zhǎng)度a可以自由變化。 2、isif=2，NSW>0，只能滿足條件<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ì)能曲線來(lái)得到極小點(diǎn)，也可以通過(guò)給定某個(gè)結(jié)構(gòu)的能量和力來(lái)通過(guò)優(yōu)化算法確定極小點(diǎn)，通常后者肯定效率要更高一些。對(duì)于尋找最優(yōu)晶格的情況，擬合EOS就是做E-V曲線找極小，但直接擬合需要計(jì)算很多點(diǎn)才能完成；這個(gè)問(wèn)題本身就是個(gè)一維極小問(wèn)題，通過(guò)一維線搜索算法可以更高效進(jìn)行解決，而且不需要假定晶格的E-V曲線長(zhǎng)什么樣。如果晶格不是很復(fù)雜，直接用isif=3穩(wěn)定的情況下應(yīng)該更合適一些。 |
2樓2019-07-23 22:53:23
卡開(kāi)發(fā)發(fā)
專家顧問(wèn) (著名寫(xiě)手)
Ab Initio Amateur
-
專家經(jīng)驗(yàn): +224 - 1ST強(qiáng)帖: 2
- 應(yīng)助: 47 (小學(xué)生)
- 金幣: 3503.8
- 散金: 14
- 紅花: 223
- 帖子: 2477
- 在線: 1254小時(shí)
- 蟲(chóng)號(hào): 1369423
- 注冊(cè): 2011-08-16
- 性別: GG
- 專業(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
waveflying
新蟲(chóng) (著名寫(xiě)手)
- 應(yīng)助: 12 (小學(xué)生)
- 金幣: 182
- 散金: 3611
- 紅花: 23
- 沙發(fā): 51
- 帖子: 2945
- 在線: 1696.5小時(shí)
- 蟲(chóng)號(hào): 1965238
- 注冊(cè): 2012-08-31
- 性別: GG
- 專業(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
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 306求調(diào)劑 +7 | Bahati 2026-03-05 | 7/350 |
|
|---|---|---|---|---|
|
|
[考研] 0805總分292，求調(diào)劑 +13 | 幻想之殤 2026-03-01 | 13/650 |
|
|
[考研] 材料與化工328求調(diào)劑 +4 | 。，。，。，。i 2026-03-02 | 4/200 |
|
|
|
[考研] 材料277分求調(diào)劑 +13 | 飯飯星球 2026-03-04 | 14/700 |
|
|
|
[考研] 材料與化工354調(diào)劑 +4 | Lucy-xiao 2026-03-06 | 7/350 |
|
|
|
[考研] 0856材料與化工求調(diào)劑！ +5 | 化工考生111 2026-03-04 | 11/550 |
|
|
|
[考研] 080500材料科學(xué)與工程 +13 | 202114020319 2026-03-03 | 13/650 |
|
|
|
[考研] 環(huán)境調(diào)劑 +7 | 柒槿levana 2026-03-01 | 7/350 |
|
|
|
[考研] 295求調(diào)劑 +4 | 小賽不吃香菜 2026-03-04 | 4/200 |
|
|
|
[考研] 312求調(diào)劑 +8 | 醋精華了一下發(fā)?/a> 2026-03-03 | 10/500 |
|
|
|
[考研] 331求調(diào)劑 +3 | zzZ&zZ 2026-03-03 | 3/150 |
|
|
|
[考研] 0703 總分319求調(diào)劑 +4 | Xinyuu 2026-03-02 | 4/200 |
|
|
|
[考研] 求調(diào)劑 +4 | Guo_yuxuan 2026-03-02 | 5/250 |
|
|
|
[考研]
 計(jì)算機(jī)學(xué)碩 分?jǐn)?shù)285求調(diào)劑
+4
|
glwshine 2026-03-02 | 5/250 |
|
|
|
[考研] 0857調(diào)劑 +6 | 一ll半 2026-02-28 | 8/400 |
|
|
|
[考研] 0856求調(diào)劑285 +11 | 呂仔龍 2026-02-28 | 11/550 |
|
|
|
[考研] 材料085601調(diào)劑 +5 | 多多子. 2026-03-02 | 5/250 |
|
|
|
[考研] 一志愿東北大學(xué)材料專碩328，求調(diào)劑 +3 | shs1083 2026-03-02 | 3/150 |
|
|
|
[考研] 317一志愿華南理工電氣工程求調(diào)劑 +6 | Soliloquy_Q 2026-02-28 | 11/550 |
|
|
|
[考研]
尋找調(diào)劑
+4
|
LYidhsjabdj 2026-02-28 | 4/200 |
|