 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 2498 | 回復(fù): 4 | |||||
[求助]
第一性原理計算 EOS擬合優(yōu)化晶胞 已有1人參與
|
|
對于EOS擬合優(yōu)化晶胞一直存在著很多疑問： 1. 大師兄的教程是計算不同體積晶胞的單點能(ISIF=2, NSW=0)，將獲得的能量和體積進行BM方程擬合來獲得平衡晶格常數(shù)； 2. 侯柱峰老師的方法是對不同體積的晶胞先進行一次離子弛豫（ISIF=2, 4或5），之后再進行靜態(tài)計算；最后靜態(tài)計算得到的能量再和體積進行BM方程擬合； ************* 這兩種方法的區(qū)別主要在于原子弛豫的順序，那在計算當中我們應(yīng)該選擇哪種方法呢？ |
回復(fù)此樓» 猜你喜歡
085600材料調(diào)劑 總分330
已經(jīng)有6人回復(fù)
-
歡迎211本科同學，過A區(qū)國家線，A區(qū)非偏遠一本，交叉學科課題組
已經(jīng)有30人回復(fù)
-
高分子化學與物理調(diào)劑
已經(jīng)有17人回復(fù)
-
化學專業(yè)調(diào)劑
已經(jīng)有6人回復(fù)
-
0703 東華大學 理學化學均過 a 線，總分 281 求調(diào)劑
已經(jīng)有11人回復(fù)
-
材料工程調(diào)劑
已經(jīng)有6人回復(fù)
-
材料學碩318求調(diào)劑
已經(jīng)有13人回復(fù)
-
08工科求調(diào)劑
已經(jīng)有4人回復(fù)
-
A區(qū)一本交叉課題組，低分調(diào)劑，招收機械電子信息通信等交叉方向
已經(jīng)有40人回復(fù)
-
070300調(diào)劑一志愿南京大學化學專業(yè)，初試273分，有論文有專利有競賽，求調(diào)劑B區(qū)211
已經(jīng)有11人回復(fù)
卡開發(fā)發(fā)
專家顧問 (著名寫手)
Ab Initio Amateur
-
專家經(jīng)驗: +224 - 1ST強帖: 2
- 應(yīng)助: 47 (小學生)
- 金幣: 3503.8
- 散金: 14
- 紅花: 223
- 帖子: 2477
- 在線: 1254小時
- 蟲號: 1369423
- 注冊: 2011-08-16
- 性別: GG
- 專業(yè): 金屬材料的磨損與磨蝕
- 管轄: 第一性原理
|
原理上應(yīng)該這么考慮，如果晶體處于穩(wěn)定結(jié)構(gòu)： <1>能量E對每個原子的核坐標xN有(∂E/∂xN)=0； <2>stress=0； <3>如果僅考慮<2>取trace的情況，能得到體積Ω滿足p=-(∂E/∂Ω)=0； 討論： 1、NSW=0，EOS找極小的點只是能夠滿足條件<3>，無法滿足條件<1>、<2>。唯一成立的條件是通過對稱性能把晶格參數(shù)和3N個原子坐標降低到只有一個獨立變量，比如Fe這種BCC或者Cu這樣的FCC的簡單格子，約束對稱性的情況只有晶格長度a可以自由變化。 2、isif=2，NSW>0，只能滿足條件<1>、<3>，此時對于三個長度參數(shù)任意或三個夾角任意的情況就不是很合理。成立的條件是能將晶格參數(shù)約化到一個獨立變量，比如一般的立方晶格應(yīng)該都行。 3、isif=4，我認為還算原理上合理。 如果只要最優(yōu)晶格的情況，我覺得這些擬合方法也不算是很好的方法。比如我們?nèi)?yōu)化H2分子來得到平衡鍵長，你可以計算不同鍵長的H2然后擬合勢能曲線來得到極小點，也可以通過給定某個結(jié)構(gòu)的能量和力來通過優(yōu)化算法確定極小點，通常后者肯定效率要更高一些。對于尋找最優(yōu)晶格的情況，擬合EOS就是做E-V曲線找極小，但直接擬合需要計算很多點才能完成；這個問題本身就是個一維極小問題，通過一維線搜索算法可以更高效進行解決，而且不需要假定晶格的E-V曲線長什么樣。如果晶格不是很復(fù)雜，直接用isif=3穩(wěn)定的情況下應(yīng)該更合適一些。 |
2樓2019-07-23 22:53:23
3樓2019-07-24 08:48:49
卡開發(fā)發(fā)
專家顧問 (著名寫手)
Ab Initio Amateur
-
專家經(jīng)驗: +224 - 1ST強帖: 2
- 應(yīng)助: 47 (小學生)
- 金幣: 3503.8
- 散金: 14
- 紅花: 223
- 帖子: 2477
- 在線: 1254小時
- 蟲號: 1369423
- 注冊: 2011-08-16
- 性別: GG
- 專業(yè): 金屬材料的磨損與磨蝕
- 管轄: 第一性原理
|
1、原理上一般復(fù)雜的情況其實isif=3最方便，如果把H2分子的掃描和優(yōu)化作為例子來看，isif=3相當于是那個優(yōu)化，isif=4+EOS就相當于掃描；同樣，做isif=5和做NSW=0的掃描也是這樣的關(guān)系，但兩者均缺乏原子弛豫，而isif=2與這兩個正好相反，僅考慮了原子弛豫而不考慮晶格形狀的變化。 2、NSW先優(yōu)化晶格，可能得到的晶格對于當前離子分數(shù)坐標固定的情況是最優(yōu)，在進行離子弛豫之后晶格就又不是最優(yōu)的了，要找的點其實應(yīng)該是兩者同時都是最優(yōu)。對于isif=4掃描體積的情況其實是一直約束原子坐標和晶格形狀是最優(yōu)的，所以才比較合理。 |
4樓2019-07-24 15:11:06
waveflying
新蟲 (著名寫手)
- 應(yīng)助: 12 (小學生)
- 金幣: 182
- 散金: 3611
- 紅花: 23
- 沙發(fā): 51
- 帖子: 2945
- 在線: 1696.5小時
- 蟲號: 1965238
- 注冊: 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
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 Cythia_cd4 的主題更新
| 5 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研]
 材料工程調(diào)劑
+6
|
zbcm_zbcm 2026-03-05 | 6/300 |
|
|---|---|---|---|---|
|
[考研] A區(qū)一本交叉課題組，低分調(diào)劑，招收機械電子信息通信等交叉方向 +28 | lisimayy 2026-03-04 | 40/2000 |
|
|
|
[考研] 085602 293分求調(diào)劑 +3 | SivanNano. 2026-03-05 | 3/150 |
|
|
|
[考研] 接收調(diào)劑 +19 | 津萌津萌 2026-03-02 | 29/1450 |
|
|
|
[考研] 求調(diào)劑 +3 | 泡了個椒 2026-03-04 | 4/200 |
|
|
|
[考研] 304求調(diào)劑 +3 | 曼殊2266 2026-03-04 | 3/150 |
|
|
|
[考研] 中國科學技術(shù)大學材料與化工281求調(diào)劑，有科研和獲獎經(jīng)歷 +11 | wsxw 2026-03-02 | 12/600 |
|
|
|
[考研] 材料學碩080500復(fù)試調(diào)劑294 +3 | 四葉zjz 2026-03-04 | 3/150 |
|
|
|
[考研] 成績276，專業(yè)代碼0856求調(diào)劑 +10 | 小陳朵 2026-03-03 | 10/500 |
|
|
|
[考研] 0855機械工程324求調(diào)劑 +3 | 鄭SK 2026-03-02 | 7/350 |
|
|
|
[考研] 材料化工調(diào)劑 +15 | 今夏不夏 2026-03-01 | 18/900 |
|
|
|
[考研] 一志愿西交大材料學碩 346 求調(diào)劑 +3 | zju51 2026-03-04 | 3/150 |
|
|
|
[考研]
 0703化學學碩 288分
+4
|
旅行中的紫葡萄 2026-03-03 | 4/200 |
|
|
|
[考研] 276求調(diào)劑 +8 | 路lyh123 2026-02-28 | 10/500 |
|
|
|
[考研] 298求調(diào)劑 +3 | 人間唯你是清歡 2026-03-03 | 4/200 |
|
|
|
[考研] 清華大學 材料與化工 353分求調(diào)劑 +5 | awaystay 2026-03-02 | 6/300 |
|
|
|
[考研] 求調(diào)劑 +11 | yunziaaaaa 2026-03-01 | 13/650 |
|
|
|
[考研] 一志愿華南理工大學材料與化工326分，求調(diào)劑 +3 | wujinrui1 2026-02-28 | 3/150 |
|
|
|
[考研] 303求調(diào)劑 +5 | 今夏不夏 2026-03-01 | 5/250 |
|
|
|
[考研] 化工專碩342，一志愿大連理工大學，求調(diào)劑 +6 | kyf化工 2026-02-28 | 7/350 |
|