 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 9 | 1/1 | 返回列表 |
| 查看: 2152 | 回復(fù): 8 | |||
qphll金蟲 (正式寫手)
|
[交流]
【求助】關(guān)于結(jié)構(gòu)參數(shù)的優(yōu)化
|
|
一直有個問題不是很清楚,請有經(jīng)驗的大家解答一下. 在做結(jié)構(gòu)參數(shù)優(yōu)化的時候,如果需要優(yōu)化的參數(shù)不止一個,那么應(yīng)該怎樣? 比如,如果需要對Orthorhombic(正交)盒子進行優(yōu)化, 那么就會涉及到需要優(yōu)化的參數(shù)有三個. 如果是在pwscf中,那么就是a, b/a和c/a. 對于我而已,我是想研究小分子在表面的分解, 那么一般c取到一定大以后,就不去考慮c/a的優(yōu)化了. 那么在實際的操作中怎么分別優(yōu)化a和b/a? 我的理解和做法是這樣的: 分兩步來 (1) 固定b/a的值,看a變化的過程中,體系總能量的變化,從而找到a的優(yōu)化值; (2) 固定a值為(1)中找到的最優(yōu)值,然后調(diào)節(jié)b/a的值, (事實上這里就是等同于改變b的值), 然后找到b的最優(yōu)值. 請大家評價一下,是不是就是這樣做的? 還有,如果我是這樣兩步做的: (1) 固定b的值,看a變化的過程中,體系總能量的變化,從而找到a的優(yōu)化值; (2) 固定a值為(1)中找到的最優(yōu)值,然后改變b的值, 然后找到b的最優(yōu)值. 如果這樣做,那么和前面的那種優(yōu)化方法有啥本質(zhì)的區(qū)別嗎? 搬個板凳,帶著筆記,等大家拍磚. |
回復(fù)此樓» 猜你喜歡
材料考研調(diào)劑
已經(jīng)有3人回復(fù)
-
材料調(diào)劑
已經(jīng)有12人回復(fù)
-
英一數(shù)一408，總分284，二戰(zhàn)真誠求調(diào)劑
已經(jīng)有14人回復(fù)
-
085410 一志愿211 22408分數(shù)359求調(diào)劑
已經(jīng)有4人回復(fù)
-
271求調(diào)劑
已經(jīng)有19人回復(fù)
-
385分 生物學（071000）求調(diào)劑
已經(jīng)有3人回復(fù)
-
一志愿安徽大學計算機科學與技術(shù)學碩，331分求調(diào)劑
已經(jīng)有3人回復(fù)
-
318求調(diào)劑，計算材料方向
已經(jīng)有8人回復(fù)
-
291求調(diào)劑
已經(jīng)有25人回復(fù)
-
一志愿北京科技大學085601材料工程英一數(shù)二初試總分335求調(diào)劑
已經(jīng)有6人回復(fù)
» 本主題相關(guān)商家推薦: (我也要在這里推廣)
» 本主題相關(guān)價值貼推薦，對您同樣有幫助:
- MS 結(jié)構(gòu)優(yōu)化的結(jié)果不能理解 求高人指點
已經(jīng)有10人回復(fù)
- 關(guān)于VASP一些瑣事的求助（K點，吸附優(yōu)化等）
已經(jīng)有12人回復(fù)
- 結(jié)構(gòu)優(yōu)化參數(shù)設(shè)置
已經(jīng)有6人回復(fù)
- 關(guān)于建模結(jié)構(gòu)優(yōu)化的一些問題
已經(jīng)有11人回復(fù)
- 【活動or資源】vasp晶體結(jié)構(gòu)優(yōu)化
已經(jīng)有142人回復(fù)
- 請大家?guī)蛶兔�，看看這CASTEP結(jié)構(gòu)優(yōu)化怎么失敗了
已經(jīng)有11人回復(fù)
- 【原創(chuàng)】關(guān)于結(jié)構(gòu)優(yōu)化的總結(jié)，體會
已經(jīng)有33人回復(fù)
- 【求助】結(jié)構(gòu)優(yōu)化中，如何選擇POTIM？
已經(jīng)有7人回復(fù)
- 【求助】VASP怎么做晶格參數(shù)的結(jié)構(gòu)優(yōu)化
已經(jīng)有17人回復(fù)
- 【求助】關(guān)于優(yōu)化時對稱性是否加入的問題 能量的選擇問題 贗勢的選擇問題
已經(jīng)有5人回復(fù)
- 【求助】關(guān)于優(yōu)化結(jié)果
已經(jīng)有17人回復(fù)
- 【求助】簡單問題----關(guān)于vasp晶格優(yōu)化
已經(jīng)有7人回復(fù)
- 【求助】關(guān)于復(fù)雜的結(jié)構(gòu)參數(shù)的優(yōu)化
已經(jīng)有5人回復(fù)
- 【求助】siesta中關(guān)于幾何結(jié)構(gòu)優(yōu)化的問題【完結(jié)】
已經(jīng)有17人回復(fù)
- 【求助】DMol3優(yōu)化參數(shù)的選擇
已經(jīng)有17人回復(fù)
zxzj05
榮譽版主 (著名寫手)
- 1ST強帖: 1
- 應(yīng)助: 0 (幼兒園)
- 貴賓: 1.242
- 金幣: 5912.5
- 紅花: 2
- 帖子: 1528
- 在線: 144.5小時
- 蟲號: 406147
- 注冊: 2007-06-18
- 性別: GG
- 專業(yè): 金屬功能材料
- 管轄: 第一性原理
2樓2009-09-23 09:32:08
qphll
金蟲 (正式寫手)
- 應(yīng)助: 18 (小學生)
- 金幣: 2708.7
- 散金: 2294
- 紅花: 8
- 帖子: 554
- 在線: 227.1小時
- 蟲號: 89654
- 注冊: 2005-08-29
- 性別: GG
- 專業(yè): 碳素材料與超硬材料
3樓2009-09-23 09:37:32
y1ding
鐵桿木蟲 (著名寫手)
- 1ST強帖: 1
- 應(yīng)助: 61 (初中生)
- 貴賓: 0.33
- 金幣: 5959.3
- 散金: 1
- 紅花: 21
- 帖子: 1884
- 在線: 491.1小時
- 蟲號: 142265
- 注冊: 2005-12-21
- 專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)
4樓2009-09-23 10:52:22
xirainbow
木蟲 (正式寫手)
- 1ST強帖: 2
- 應(yīng)助: 0 (幼兒園)
- 貴賓: 0.03
- 金幣: 2132.6
- 紅花: 2
- 帖子: 775
- 在線: 6.1小時
- 蟲號: 431211
- 注冊: 2007-08-11
- 專業(yè): 光譜學
★ ★ ★ ★
wuchenwf(金幣+2,VIP+0):謝謝 9-24 22:00
qphll(金幣+2,VIP+0): 9-25 11:42
wuchenwf(金幣+2,VIP+0):謝謝 9-24 22:00
qphll(金幣+2,VIP+0): 9-25 11:42
|
PWscf的參數(shù)“cell_dofree”專門處理你提出的問題，只優(yōu)化a,b而固定c 但是這個方法bfgs好像是暫時還沒有實現(xiàn)，你需使用共軛梯度法 你可以看看這個：“http://blog.sina.com.cn/s/blog_5f15ead20100e1lu.html###”（Some fundamental concept and parameters of pwscf） 里面有關(guān)于此問題PWscf forum的網(wǎng)址 |
5樓2009-09-24 09:47:17
qphll
金蟲 (正式寫手)
- 應(yīng)助: 18 (小學生)
- 金幣: 2708.7
- 散金: 2294
- 紅花: 8
- 帖子: 554
- 在線: 227.1小時
- 蟲號: 89654
- 注冊: 2005-08-29
- 性別: GG
- 專業(yè): 碳素材料與超硬材料
|
多謝樓上三位的參與. 對于我自己的特定體系而已, 我發(fā)帖的里面提到的結(jié)構(gòu)參數(shù)的優(yōu)化問題已經(jīng)想明白了. 對于graphene而言,因為我自己編程生成graphene的原因,我的盒子是Orthorhombic的,a,b,c三個值各不相等,alpha, beta, gamma 都是90度. 如果想想結(jié)構(gòu)優(yōu)化的目的,我想對于我的體系而已,無非是因為在生成graphene的時候, C-C的鍵長使用了實驗值1.42 那么這里結(jié)構(gòu)優(yōu)化的目的,就是通過盒子的'擴張'或者'收縮', 找到一個合適的C-C鍵長而已. 如果keep that in mind, 那么雖然我這里a和b的數(shù)值是不相等的,但是其實它們也不是互相獨立的.如果我將a (X軸)拉伸,那么為了維持graphene中C-C鍵長的一致性, b(Y軸)必定得要按比例拉伸. 這里c (Z軸)在取到一定值以后,保證graphene面和它的Z軸鏡像沒有相互左右以后,就不需要優(yōu)化啦. 所以我想,對于我的這個長方體的模擬盒子,放在pwscf里面考慮lattice parameter優(yōu)化的話,是可以這么一步做的: (在in文件中,有celldm(1)=a, celldm(2)=b/a,celldm(3)=c/a 這三個參數(shù)) 固定b/a的值,固定c的值, 看a變化的過程中,體系總能量的變化,從而找到a的優(yōu)化值 (事實上b/a固定的話,那么a的最優(yōu)值找到了,那么b的值也相應(yīng)找到了). 超級新手, pw里面的vc-relax沒有嘗試,看manual也似乎沒說啥,例子也沒有找到. y1ding可否提供一個in文件給我學習學習? |
6樓2009-09-25 12:02:36
7樓2009-09-28 15:30:47
zxzj05
榮譽版主 (著名寫手)
- 1ST強帖: 1
- 應(yīng)助: 0 (幼兒園)
- 貴賓: 1.242
- 金幣: 5912.5
- 紅花: 2
- 帖子: 1528
- 在線: 144.5小時
- 蟲號: 406147
- 注冊: 2007-06-18
- 性別: GG
- 專業(yè): 金屬功能材料
- 管轄: 第一性原理
8樓2009-09-28 17:20:17
改王之王
金蟲 (正式寫手)
- 1ST強帖: 1
- 應(yīng)助: 0 (幼兒園)
- 金幣: 1120.7
- 散金: 61
- 紅花: 3
- 帖子: 817
- 在線: 113.5小時
- 蟲號: 968169
- 注冊: 2010-03-11
- 專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)
9樓2011-03-31 12:35:32
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 qphll 的主題更新
| 9 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿北京科技，085601總分305求調(diào)劑 +9 | 半生瓜！ 2026-04-01 | 10/500 |
|
|---|---|---|---|---|
|
[考研] 282求調(diào)劑 +17 | ycy1201 2026-04-01 | 19/950 |
|
|
|
[考研] 289求調(diào)劑 +21 | 新時代材料 2026-03-27 | 23/1150 |
|
|
|
[考研] 0710生物學求調(diào)劑 +6 | manman511 2026-04-01 | 6/300 |
|
|
|
[考研] 0817化工學碩調(diào)劑 +11 | 努力上岸中！ 2026-03-31 | 11/550 |
|
|
|
[考研] 一志愿西安交大材料學碩（英一數(shù)二）347，求調(diào)劑到高分子/材料相關(guān)專業(yè) +7 | zju51 2026-03-31 | 9/450 |
|
|
|
[考研] 一志愿鄭大085600，310分求調(diào)劑 +6 | 李瀟可 2026-03-26 | 6/300 |
|
|
|
[考研] 322求調(diào)劑 +8 | 三水sss 2026-04-01 | 8/400 |
|
|
|
[考研] 318一志愿吉林大學生物與醫(yī)藥 求調(diào)劑 +6 | 篤行致遠. 2026-03-28 | 6/300 |
|
|
|
[考研] 085701環(huán)境工程，267求調(diào)劑 +17 | minht 2026-03-26 | 17/850 |
|
|
|
[考研] 考研生物與醫(yī)藥調(diào)劑 +7 | 鐵憨憨123425 2026-03-31 | 7/350 |
|
|
|
[考研] 322求調(diào)劑 +4 | 熹僖XX 2026-03-31 | 4/200 |
|
|
|
[考研] 318求調(diào)劑 +3 | 篤行致遠. 2026-03-31 | 3/150 |
|
|
|
[考研] 086000調(diào)劑 +5 | 7901117076 2026-03-26 | 5/250 |
|
|
|
[考研] 吉大生物學326分求調(diào)劑 +3 | sunnyupup 2026-03-31 | 3/150 |
|
|
|
[考研] 一志愿北京理工大學本科211材料工程294求調(diào)劑 +8 | mikasa的圍巾 2026-03-28 | 8/400 |
|
|
|
[考研] 315分求調(diào)劑 +7 | 26考研上岸版26 2026-03-26 | 7/350 |
|
|
|
[考研] 070300化學求調(diào)劑 +4 | 起個名咋這么難 2026-03-27 | 4/200 |
|
|
|
[考研] 0856調(diào)劑 +5 | 求求讓我有書讀?/a> 2026-03-26 | 6/300 |
|
|
|
[考研] 081200-11408-276學碩求調(diào)劑 +4 | 崔wj 2026-03-26 | 4/200 |
|