 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 查看: 2476 | 回復(fù): 11 | ||||||
| 本帖產(chǎn)生 2 個(gè) 1ST強(qiáng)帖 ，點(diǎn)擊這里進(jìn)行查看 | ||||||
navyzhang木蟲 (正式寫手)
|
[求助]
最復(fù)雜情況下的參數(shù)優(yōu)化
|
|||||
|
看了侯老師的VASP教程，關(guān)于晶格參數(shù)優(yōu)化存在一些疑問！ 對(duì)于a,b,c都不相等的體系，是不是可以這樣 1、先優(yōu)化形狀(即c/a,b/a)和位置，(將該計(jì)算得到的CONTCAR作為下一步的POSCAR) 2、僅優(yōu)化位置，得到該情況下的能量E1 3、改變a參數(shù)，按上面兩步，先優(yōu)化形狀和位置，再單優(yōu)化位置，得到該情況下的能量E2 4、如此往復(fù)共改變十次a參數(shù)，得到是個(gè)能量E1~E10，擬合E-V曲線，得到能量最低時(shí)的V0 那么下面參數(shù)a到底該怎么求呢？將V0/(b*c)？如果是這樣，b和c從哪里來？是能量最低情況下第一步對(duì)應(yīng)的結(jié)果？ |
回復(fù)此樓» 收錄本帖的淘帖專輯推薦
 第一原理資料匯編 | 計(jì)算-vasp | 分子模擬專題 | vasp從結(jié)構(gòu)優(yōu)化到光學(xué)性質(zhì) |
» 猜你喜歡
材料專業(yè)求調(diào)劑
已經(jīng)有5人回復(fù)
-
一志愿天津大學(xué)化學(xué)工藝專業(yè)（081702）315分求調(diào)劑
已經(jīng)有10人回復(fù)
-
344求調(diào)劑
已經(jīng)有7人回復(fù)
-
一志愿985，本科211，0817化學(xué)工程與技術(shù)319求調(diào)劑
已經(jīng)有7人回復(fù)
-
271材料工程求調(diào)劑
已經(jīng)有4人回復(fù)
-
材料專碩英一數(shù)二306
已經(jīng)有4人回復(fù)
-
294求調(diào)劑材料與化工專碩
已經(jīng)有4人回復(fù)
-
一志愿中海洋材料工程專碩330分求調(diào)劑
已經(jīng)有4人回復(fù)
-
070303一志愿西北大學(xué)學(xué)碩310找調(diào)劑
已經(jīng)有9人回復(fù)
-
0703化學(xué)調(diào)劑，求各位老師收留
已經(jīng)有9人回復(fù)
» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:
- 結(jié)構(gòu)優(yōu)化參數(shù)設(shè)置
已經(jīng)有6人回復(fù)
- 求回歸方程參數(shù)估計(jì)的最優(yōu)化算法matlab代碼
已經(jīng)有9人回復(fù)
- LC-MS 優(yōu)化方法時(shí)用到的一些參數(shù)DP, EP, FP, CXP等分別起什么作用？
已經(jīng)有7人回復(fù)
- 結(jié)構(gòu)優(yōu)化參數(shù)設(shè)置
已經(jīng)有18人回復(fù)
- 【求助】用matlab最優(yōu)化方法進(jìn)行參數(shù)擬合
已經(jīng)有17人回復(fù)
- 【求助】castep做幾何優(yōu)化時(shí)其參數(shù)如何設(shè)置？
已經(jīng)有19人回復(fù)
- 【求助】vasp優(yōu)化晶格參數(shù)a
已經(jīng)有13人回復(fù)
- 【求助】關(guān)于復(fù)雜的結(jié)構(gòu)參數(shù)的優(yōu)化
已經(jīng)有5人回復(fù)
- 【求助】Castep中結(jié)構(gòu)優(yōu)化后晶胞參數(shù)差別較大 vasp高手也幫我看一下
已經(jīng)有22人回復(fù)
- 【求助】DMol3優(yōu)化參數(shù)的選擇
已經(jīng)有17人回復(fù)
qiqi2926
銅蟲 (小有名氣)
- 1ST強(qiáng)帖: 1
- 應(yīng)助: 5 (幼兒園)
- 金幣: 401.2
- 散金: 15
- 紅花: 3
- 帖子: 264
- 在線: 87.8小時(shí)
- 蟲號(hào): 1004137
- 注冊(cè): 2010-04-23
- 性別: GG
- 專業(yè): 原子和分子物理
【答案】應(yīng)助回帖
★ ★ ★
youzhizhe(金幣+3): 謝謝交流。 2011-05-26 12:27:49
youzhizhe(1ST強(qiáng)帖+1): 謝謝交流。 2011-05-26 16:41:46
youzhizhe(金幣+3): 謝謝交流。 2011-05-26 12:27:49
youzhizhe(1ST強(qiáng)帖+1): 謝謝交流。 2011-05-26 16:41:46
|
首先，結(jié)構(gòu)優(yōu)化是為了尋找晶體的基態(tài)。 在侯博士的那片文章中提到的對(duì)復(fù)雜晶體的優(yōu)化，再腳本中的a首先是在實(shí)驗(yàn)值的左右各取若干個(gè)，然后對(duì)晶格采用固定體積的方法進(jìn)行優(yōu)化，這樣就會(huì)得到一系列的a-E圖像，做出曲線，這樣的話，你就能看到這個(gè)曲線呈現(xiàn)U型。我們知道基態(tài)的能量應(yīng)該是最低的，所以a的取值應(yīng)該是對(duì)應(yīng)于E最小的a，但是一次并不能精確的確定a的值，只能錯(cuò)略的估計(jì)a的范圍，但是通過多次的縮小這個(gè)區(qū)間，就能得到較為準(zhǔn)確的a值，而由于是允許晶格形可變，且原子位置馳豫的，所以就可以得到基態(tài)結(jié)構(gòu)。 |
2樓2011-05-26 11:37:46
navyzhang
木蟲 (正式寫手)
- 應(yīng)助: 20 (小學(xué)生)
- 金幣: 1411.8
- 散金: 1197
- 紅花: 5
- 帖子: 640
- 在線: 187.1小時(shí)
- 蟲號(hào): 883307
- 注冊(cè): 2009-10-25
- 性別: GG
- 專業(yè): 理論和計(jì)算化學(xué)
3樓2011-05-26 15:17:14
qiqi2926
銅蟲 (小有名氣)
- 1ST強(qiáng)帖: 1
- 應(yīng)助: 5 (幼兒園)
- 金幣: 401.2
- 散金: 15
- 紅花: 3
- 帖子: 264
- 在線: 87.8小時(shí)
- 蟲號(hào): 1004137
- 注冊(cè): 2010-04-23
- 性別: GG
- 專業(yè): 原子和分子物理
4樓2011-05-26 15:28:44
tider
木蟲 (正式寫手)
- 應(yīng)助: 8 (幼兒園)
- 金幣: 1950.4
- 散金: 40
- 紅花: 5
- 帖子: 313
- 在線: 87.2小時(shí)
- 蟲號(hào): 1190614
- 注冊(cè): 2011-01-16
- 性別: GG
- 專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)
5樓2011-05-26 18:47:07
nicolemikeyi
禁蟲 (著名寫手)
- 1ST強(qiáng)帖: 1
- 應(yīng)助: 6 (幼兒園)
- 金幣: 5725.1
- 散金: 94
- 紅花: 15
- 帖子: 2539
- 在線: 616.8小時(shí)
- 蟲號(hào): 1241122
- 注冊(cè): 2011-03-22
- 性別: GG
- 專業(yè): 固體力學(xué)
【答案】應(yīng)助回帖
★ ★ ★
gzqdyouxia(金幣+3, 1ST強(qiáng)帖+1): 謝謝詳細(xì)指教 2011-05-27 10:43:47
gzqdyouxia(金幣+3, 1ST強(qiáng)帖+1): 謝謝詳細(xì)指教 2011-05-27 10:43:47
|
那里頭的腳本分為兩部分，一部分是固定晶胞體積的弛豫，得到CONTCAR后復(fù)制到POSCAR中進(jìn)行另一部分靜態(tài)計(jì)算，得出能量。所以侯的優(yōu)化思路是這樣的： 1.取不同的a值，首先進(jìn)行ISIF=5的固定體積但不固定形狀的弛豫，然后靜態(tài)計(jì)算出對(duì)應(yīng)的能量。 2.每一個(gè)a值對(duì)應(yīng)一個(gè)體積和能量，多個(gè)a值可以得到一條E~V曲線，通過擬合可以得到最小的那個(gè)E對(duì)應(yīng)的體積和a值。 3.注意，此時(shí)如果是立方晶格，a值即為優(yōu)化結(jié)果。但例如是四方或其他，此時(shí)你只是得到了最佳的形狀（a/c值），你再取不同的a進(jìn)行下一次循環(huán)，以得到最佳的a值與c值。 |
» 本帖已獲得的紅花（最新10朵）
6樓2011-05-26 20:21:07
nicolemikeyi
禁蟲 (著名寫手)
navyzhang(金幣+10): 謝謝 2011-06-02 15:58:59
|
本帖內(nèi)容被屏蔽 |
7樓2011-05-26 20:28:24
navyzhang
木蟲 (正式寫手)
- 應(yīng)助: 20 (小學(xué)生)
- 金幣: 1411.8
- 散金: 1197
- 紅花: 5
- 帖子: 640
- 在線: 187.1小時(shí)
- 蟲號(hào): 883307
- 注冊(cè): 2009-10-25
- 性別: GG
- 專業(yè): 理論和計(jì)算化學(xué)
8樓2011-06-02 15:58:41
liliangfang(金幣+1): 多謝交流 2012-02-26 08:52:35
|
本帖內(nèi)容被屏蔽 |
9樓2012-02-26 00:16:07
洋洋1984
新蟲 (正式寫手)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 1791.1
- 散金: 143
- 紅花: 2
- 帖子: 474
- 在線: 113.1小時(shí)
- 蟲號(hào): 1303656
- 注冊(cè): 2011-05-23
- 性別: GG
- 專業(yè): 數(shù)值預(yù)報(bào)與數(shù)值模擬
送鮮花一朵 |
nicolemikeyi 你好，在下有幾個(gè)問題想請(qǐng)教一下，麻煩了。 1.對(duì)于在您的回復(fù)：2.每一個(gè)a值對(duì)應(yīng)一個(gè)體積和能量，多個(gè)a值可以得到一條E~V曲線，通過擬合可以得到最小的那個(gè)E對(duì)應(yīng)的體積和a值。此處您用的是Birch-Murnaghan方程還是二次曲線擬合的？ 2.在結(jié)果中得到了a值，但為何是“此時(shí)你只是得到了最佳的形狀（a/c值）”？如果直接用a/(a/c)是否可以認(rèn)為得到了c值？你再取不同的a進(jìn)行下一次循環(huán)，以得到最佳的a值與c值？ 麻煩指教一下，謝謝！  |
10樓2012-05-09 11:24:59
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 navyzhang 的主題更新
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿天津大學(xué)化學(xué)工藝專業(yè)（081702）315分求調(diào)劑 +10 | yangfz 2026-03-17 | 10/500 |
|
|---|---|---|---|---|
|
[考研] 一志愿中海洋材料工程專碩330分求調(diào)劑 +4 | 小材化本科 2026-03-18 | 4/200 |
|
|
|
[考研] 070303一志愿西北大學(xué)學(xué)碩310找調(diào)劑 +6 | d如愿上岸 2026-03-12 | 9/450 |
|
|
|
[考研] 求材料調(diào)劑 +7 | 隔壁陳先生 2026-03-12 | 7/350 |
|
|
|
[考研] 288求調(diào)劑，一志愿華南理工大學(xué)071005 +4 | ioodiiij 2026-03-17 | 4/200 |
|
|
|
[考研] 生物學(xué)071000 329分求調(diào)劑 +3 | 我愛生物生物愛?/a> 2026-03-17 | 3/150 |
|
|
|
[考研] 334求調(diào)劑 +3 | 志存高遠(yuǎn)意在機(jī)?/a> 2026-03-16 | 3/150 |
|
|
|
[考研] 301求調(diào)劑 +4 | A_JiXing 2026-03-16 | 4/200 |
|
|
|
[考研] 梁成偉老師課題組歡迎你的加入 +8 | 一鴨鴨喲 2026-03-14 | 10/500 |
|
|
|
[考研] 一志愿，福州大學(xué)材料專碩339分求調(diào)劑 +3 | 木子momo青爭(zhēng) 2026-03-15 | 3/150 |
|
|
|
[考研]
 11408 一志愿西電，277分求調(diào)劑
+3
|
zhouzhen654 2026-03-16 | 3/150 |
|
|
|
[考研] 東南大學(xué)364求調(diào)劑 +5 | JasonYuiui 2026-03-15 | 5/250 |
|
|
|
[考研] 333求調(diào)劑 +3 | 文思客 2026-03-16 | 7/350 |
|
|
|
[考研] 085600材料與化工 求調(diào)劑 +13 | enenenhui 2026-03-13 | 14/700 |
|
|
|
[考研] 294求調(diào)劑 +3 | Zys010410@ 2026-03-13 | 4/200 |
|
|
|
[碩博家園] 085600 260分求調(diào)劑 +3 | 天空還下雨么 2026-03-13 | 5/250 |
|
|
|
[考研] 290求調(diào)劑 +7 | ADT 2026-03-12 | 7/350 |
|
|
|
[考研] 274求調(diào)劑 +3 | S.H1 2026-03-12 | 3/150 |
|
|
|
[考研] 材料301分求調(diào)劑 +5 | Liyouyumairs 2026-03-12 | 5/250 |
|
|
|
[論文投稿]
投稿問題
 5+4
|
星光燦爛xt 2026-03-12 | 6/300 |
|