| 9 | 1/1 | 返回列表 |
| 查看: 2155 | 回復(fù): 8 | |||
qphll金蟲 (正式寫手)
|
[交流]
【求助】關(guān)于結(jié)構(gòu)參數(shù)的優(yōu)化
|
|
一直有個(gè)問題不是很清楚,請(qǐng)有經(jīng)驗(yàn)的大家解答一下. 在做結(jié)構(gòu)參數(shù)優(yōu)化的時(shí)候,如果需要優(yōu)化的參數(shù)不止一個(gè),那么應(yīng)該怎樣? 比如,如果需要對(duì)Orthorhombic(正交)盒子進(jìn)行優(yōu)化, 那么就會(huì)涉及到需要優(yōu)化的參數(shù)有三個(gè). 如果是在pwscf中,那么就是a, b/a和c/a. 對(duì)于我而已,我是想研究小分子在表面的分解, 那么一般c取到一定大以后,就不去考慮c/a的優(yōu)化了. 那么在實(shí)際的操作中怎么分別優(yōu)化a和b/a? 我的理解和做法是這樣的: 分兩步來 (1) 固定b/a的值,看a變化的過程中,體系總能量的變化,從而找到a的優(yōu)化值; (2) 固定a值為(1)中找到的最優(yōu)值,然后調(diào)節(jié)b/a的值, (事實(shí)上這里就是等同于改變b的值), 然后找到b的最優(yōu)值. 請(qǐng)大家評(píng)價(jià)一下,是不是就是這樣做的? 還有,如果我是這樣兩步做的: (1) 固定b的值,看a變化的過程中,體系總能量的變化,從而找到a的優(yōu)化值; (2) 固定a值為(1)中找到的最優(yōu)值,然后改變b的值, 然后找到b的最優(yōu)值. 如果這樣做,那么和前面的那種優(yōu)化方法有啥本質(zhì)的區(qū)別嗎? 搬個(gè)板凳,帶著筆記,等大家拍磚. |

榮譽(yù)版主 (著名寫手)

金蟲 (正式寫手)

鐵桿木蟲 (著名寫手)
木蟲 (正式寫手)
|
PWscf的參數(shù)“cell_dofree”專門處理你提出的問題,只優(yōu)化a,b而固定c 但是這個(gè)方法bfgs好像是暫時(shí)還沒有實(shí)現(xiàn),你需使用共軛梯度法 你可以看看這個(gè):“http://blog.sina.com.cn/s/blog_5f15ead20100e1lu.html###”(Some fundamental concept and parameters of pwscf) 里面有關(guān)于此問題PWscf forum的網(wǎng)址 |
金蟲 (正式寫手)
|
多謝樓上三位的參與. 對(duì)于我自己的特定體系而已, 我發(fā)帖的里面提到的結(jié)構(gòu)參數(shù)的優(yōu)化問題已經(jīng)想明白了. 對(duì)于graphene而言,因?yàn)槲易约壕幊躺蒰raphene的原因,我的盒子是Orthorhombic的,a,b,c三個(gè)值各不相等,alpha, beta, gamma 都是90度. 如果想想結(jié)構(gòu)優(yōu)化的目的,我想對(duì)于我的體系而已,無非是因?yàn)樵谏蒰raphene的時(shí)候, C-C的鍵長使用了實(shí)驗(yàn)值1.42 那么這里結(jié)構(gòu)優(yōu)化的目的,就是通過盒子的'擴(kuò)張'或者'收縮', 找到一個(gè)合適的C-C鍵長而已. 如果keep that in mind, 那么雖然我這里a和b的數(shù)值是不相等的,但是其實(shí)它們也不是互相獨(dú)立的.如果我將a (X軸)拉伸,那么為了維持graphene中C-C鍵長的一致性, b(Y軸)必定得要按比例拉伸. 這里c (Z軸)在取到一定值以后,保證graphene面和它的Z軸鏡像沒有相互左右以后,就不需要優(yōu)化啦. 所以我想,對(duì)于我的這個(gè)長方體的模擬盒子,放在pwscf里面考慮lattice parameter優(yōu)化的話,是可以這么一步做的: (在in文件中,有celldm(1)=a, celldm(2)=b/a,celldm(3)=c/a 這三個(gè)參數(shù)) 固定b/a的值,固定c的值, 看a變化的過程中,體系總能量的變化,從而找到a的優(yōu)化值 (事實(shí)上b/a固定的話,那么a的最優(yōu)值找到了,那么b的值也相應(yīng)找到了). 超級(jí)新手, pw里面的vc-relax沒有嘗試,看manual也似乎沒說啥,例子也沒有找到. y1ding可否提供一個(gè)in文件給我學(xué)習(xí)學(xué)習(xí)? |

榮譽(yù)版主 (著名寫手)

金蟲 (正式寫手)
| 9 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 08生物與醫(yī)藥專碩初試346找調(diào)劑 +6 | dianeeee 2026-04-01 | 7/350 |
|
|---|---|---|---|---|
|
[考研] 271求調(diào)劑 +14 | 勒布朗@ 2026-03-31 | 19/950 |
|
|
[考研] 081200-11408-276學(xué)碩求調(diào)劑 +5 | 崔wj 2026-03-26 | 5/250 |
|
|
[考研] 070300化學(xué)求調(diào)劑 +13 | 小黃鴨寶 2026-03-30 | 13/650 |
|
|
[考研] 一志愿085600中科院寧波所276分求調(diào)劑 +22 | 材料學(xué)257求調(diào)劑 2026-03-28 | 23/1150 |
|
|
[考研] 材料科學(xué)與工程339求調(diào)劑 +11 | hyz0119 2026-03-31 | 12/600 |
|
|
[考研] 調(diào)劑 +3 | 好好讀書。 2026-04-01 | 3/150 |
|
|
[考研] 一志愿 南京航空航天大學(xué) ,080500材料科學(xué)與工程學(xué)碩 +7 | @taotao 2026-03-30 | 7/350 |
|
|
[考研] 318一志愿吉林大學(xué)生物與醫(yī)藥 求調(diào)劑 +6 | 篤行致遠(yuǎn). 2026-03-28 | 6/300 |
|
|
[考研] 318求調(diào)劑 +3 | 篤行致遠(yuǎn). 2026-03-31 | 3/150 |
|
|
[考研] 一志愿華東理工大學(xué),080500學(xué)碩,317分,求調(diào)劑 +9 | s1145 2026-03-31 | 9/450 |
|
|
[考研] 一志愿a區(qū)211,085601-307分求調(diào)劑 +10 | 黨嘉豪 2026-03-31 | 23/1150 |
|
|
[考研] 本2一志愿C9-333分,材料科學(xué)與工程,求調(diào)劑 +9 | 升升不降 2026-03-31 | 9/450 |
|
|
[考研] 生醫(yī)工0831調(diào)劑求推薦 +4 | 小熊睿睿_s 2026-03-27 | 6/300 |
|
|
[考研] 考研調(diào)劑求助 +7 | 13287130938 2026-03-31 | 7/350 |
|
|
[考研] 353求調(diào)劑 +3 | 江上楓_26 2026-03-28 | 3/150 |
|
|
[考研]
|
鐘llll 2026-03-26 | 5/250 |
|
|
[考研] 340求調(diào)劑 +6 | Amber00 2026-03-26 | 6/300 |
|
|
[考研] 070300化學(xué)求調(diào)劑 +4 | 起個(gè)名咋這么難 2026-03-27 | 4/200 |
|
|
[考研] 0856調(diào)劑 +5 | 求求讓我有書讀?/a> 2026-03-26 | 6/300 |
|