 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 6 | 1/1 | 返回列表 |
| 查看: 1120 | 回復(fù): 5 | ||||
[求助]
菜鳥求助關(guān)于加壓優(yōu)化的問(wèn)題。。。
|
|
最近在做關(guān)于加壓優(yōu)化的東西， 第一個(gè)問(wèn)題：我在toturial 中沒(méi)找見具體的加壓優(yōu)化的例子（只有l(wèi)esson3中有不加壓優(yōu)化的例子），請(qǐng)問(wèn)有沒(méi)有具體關(guān)于加壓優(yōu)化的例子？有的話在那一節(jié)？ 第二個(gè)問(wèn)題：Out 文件中：壓力張量Cartesian components of stress tensor (hartree/bohr^3) sigma(1 1)= -4.76514276E-08 sigma(3 2)= 0.00000000E+00 sigma(2 2)= -4.76514276E-08 sigma(3 1)= 0.00000000E+00 sigma(3 3)= -4.76514276E-08 sigma(2 1)= 0.00000000E+00 和梯度收斂標(biāo)準(zhǔn)：At Broyd/MD step 3, gradients are converged : max grad (force/stress) = 4.7651E-06 < tolmxf= 5.0000E-05 ha/bohr (free atoms) 之間為什么呈現(xiàn)sigma(1 1)+22+33整體除以三 再乘以10的2次方的關(guān)系？按道理說(shuō)max grad的單位是 ha/boh，而壓力張量單位是hartree/bohr^3。這是不同量綱的東西�。� 三 請(qǐng)問(wèn)為什么立方晶系sigma(1 1)=22=33而四方晶系（比如金紅石）sigma(1 1)=22不等于33呢，11 22 和33之間的關(guān)系式什么呢？如果我壓價(jià)優(yōu)化金紅石晶胞時(shí)應(yīng)該怎眼設(shè)置11 22 33呢？ 最后 跪求高手了。。。。 |
回復(fù)此樓» 猜你喜歡
278求調(diào)劑
已經(jīng)有7人回復(fù)
-
26考研一志愿中國(guó)石油大學(xué)(華東)305分求調(diào)劑
已經(jīng)有5人回復(fù)
-
280求調(diào)劑
已經(jīng)有9人回復(fù)
-
279求調(diào)劑
已經(jīng)有5人回復(fù)
-
286求調(diào)劑
已經(jīng)有5人回復(fù)
-
材料與化工（0856）304求 B區(qū) 調(diào)劑
已經(jīng)有3人回復(fù)
-
268求調(diào)劑
已經(jīng)有8人回復(fù)
-
265求調(diào)劑
已經(jīng)有14人回復(fù)
-
22 350 本科985求調(diào)劑，求老登收留
已經(jīng)有3人回復(fù)
-
【考研調(diào)劑】化學(xué)專業(yè) 281分，一志愿四川大學(xué)，誠(chéng)心求調(diào)劑
已經(jīng)有9人回復(fù)
souledge
專家顧問(wèn) (著名寫手)
-
專家經(jīng)驗(yàn): +174 - 1ST強(qiáng)帖: 19
- 應(yīng)助: 303 (大學(xué)生)
- 金幣: 11306.4
- 散金: 1123
- 紅花: 108
- 帖子: 2814
- 在線: 916.1小時(shí)
- 蟲號(hào): 941486
- 注冊(cè): 2010-01-12
- 專業(yè): 結(jié)構(gòu)陶瓷
- 管轄: 第一性原理
【答案】應(yīng)助回帖
★ ★ ★ ★ ★
感謝參與，應(yīng)助指數(shù) +1
fzx2008: 金幣+5, 專家考核, 謝謝指教！ 2012-12-04 18:38:56
感謝參與，應(yīng)助指數(shù) +1
fzx2008: 金幣+5, 專家考核, 謝謝指教！ 2012-12-04 18:38:56
|
1，沒(méi)有加壓優(yōu)化的例子，因?yàn)樗^的不加壓優(yōu)化，就是為strtarget = 0的優(yōu)化，而加壓優(yōu)化，是strtarget為一個(gè)finite value的優(yōu)化，造個(gè)例子就屬于重復(fù)建設(shè)了。 2，Hartree是能量單位，可以看作是F*l的結(jié)果，Bohr是長(zhǎng)度單位，也就是l。于是Hartree/Bohr^3就可以看作是F*l/l^3=F/l^2，正好是壓強(qiáng)的單位。 3，理論上也應(yīng)該相等，但是由于四方晶的c軸與a、b非等效軸，所以難以優(yōu)化后完全一致，但是該問(wèn)題顯然在對(duì)稱性極高的立方晶上不存在。其實(shí)可以看到，如果優(yōu)化的收斂條件很高，33方向的壓強(qiáng)與11，22的差別很小。至于加壓，一般都是加等靜壓，當(dāng)然是11=22=33了。 |
2樓2012-12-04 17:43:08
 |
3樓2012-12-04 18:14:05
送鮮花一朵 |
首先 謝謝了 關(guān)于第二個(gè)問(wèn)題：我想問(wèn)max grad（就是max grad (force/stress) = 4.7651E-06），和 sigma(1 1) 22 33 的關(guān)系，這三個(gè)是壓強(qiáng)，但max grad明顯是在和tolmxf作比較，所以應(yīng)該是力。（11+22+33）除以三等于 -4.76514276E-08，這是壓強(qiáng)，但max grad (force/stress) = 4.7651E-06 ，這是力，難道說(shuō)壓強(qiáng)乘以100就是力（二者之間相差100倍），這是為什么呢？ 另外一個(gè)問(wèn)題：在計(jì)算金紅石0壓力優(yōu)化時(shí)，殘余壓力為： -Cartesian components of stress tensor (GPa) [Pressure= 4.8974E-03 GPa] - sigma(1 1)= -1.70282756E+00 sigma(3 2)= 0.00000000E+00 - sigma(2 2)= -1.70282756E+00 sigma(3 1)= 0.00000000E+00 - sigma(3 3)= 3.39096300E+00 sigma(2 1)= 0.00000000E+00 可以看到雖然整體的 [Pressure= 4.8974E-03 GPa] 相對(duì)比較小，但是11 22 和33都超過(guò)了1GPa................，這合適嗎？（就連我設(shè)置1G壓力優(yōu)化時(shí)三個(gè)分量才設(shè)為-1G，而這個(gè)是0壓力優(yōu)化�。┮簿褪钦f(shuō)我應(yīng)該看壓力張量分量呢還是以整體殘余壓力為標(biāo)準(zhǔn)？ 最后，我懸賞的10個(gè)金幣應(yīng)該怎樣支付給你？給你兩朵小紅花嗎？ |
4樓2012-12-05 10:30:33
souledge
專家顧問(wèn) (著名寫手)
-
專家經(jīng)驗(yàn): +174 - 1ST強(qiáng)帖: 19
- 應(yīng)助: 303 (大學(xué)生)
- 金幣: 11306.4
- 散金: 1123
- 紅花: 108
- 帖子: 2814
- 在線: 916.1小時(shí)
- 蟲號(hào): 941486
- 注冊(cè): 2010-01-12
- 專業(yè): 結(jié)構(gòu)陶瓷
- 管轄: 第一性原理
【答案】應(yīng)助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★
若博特: 金幣+5, ★★★很有幫助, 很熱心的專家 2012-12-14 09:03:28
若博特: 金幣+5, ★★★很有幫助, 很好的寫手，謝謝了，希望能及時(shí)回復(fù)下邊的問(wèn)題。 2012-12-15 16:20:41
若博特: 金幣+5, ★★★很有幫助, 很熱心的專家 2012-12-14 09:03:28
若博特: 金幣+5, ★★★很有幫助, 很好的寫手，謝謝了，希望能及時(shí)回復(fù)下邊的問(wèn)題。 2012-12-15 16:20:41
|
1，根據(jù)fzx2008版主的提示，該問(wèn)題可以由下面幾個(gè)信息解答： a）At Broyd/MD step 14, gradients are converged : max grad (force/stress) = 5.5774E-06 < tolmxf= 1.0000E-04 ha/bohr (free atoms) 所以max grad表示最大梯度。 b） 所以在Abinit中，直接由stress來(lái)預(yù)測(cè)力。一般系數(shù)是100. 綜上，該方法是取stress矩陣中最大的值來(lái)乘100（當(dāng)然可以少乘，自己設(shè)置strfact即可），然后與tolmxf比較。 2，LZ的優(yōu)化有點(diǎn)粗糙啊，雖然總Pressure過(guò)得去，但是，各個(gè)方向的確很大，建議讓收斂條件更嚴(yán)苛點(diǎn)，順便看看截?cái)嗄苤惖膮?shù)是否收斂。 3，這個(gè)問(wèn)題應(yīng)該問(wèn)版主而不是我…… |
5樓2012-12-05 15:06:45
6樓2012-12-10 17:12:09
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 若博特 的主題更新
| 6 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 材料與化工（0856）304求 B區(qū) 調(diào)劑 +3 | 邱gl 2026-03-21 | 3/150 |
|
|---|---|---|---|---|
|
[考研] 能源材料化學(xué)課題組招收碩士研究生8-10名 +5 | 脫穎而出 2026-03-16 | 15/750 |
|
|
|
[考研] 一志愿山大07化學(xué) 332分 四六級(jí)已過(guò) 本科山東雙非 求調(diào)劑！ +3 | 不想理你 2026-03-16 | 3/150 |
|
|
|
[考研] 297求調(diào)劑 +9 | 戲精丹丹丹 2026-03-17 | 9/450 |
|
|
|
[考研] 華東師范大學(xué)-071000生物學(xué)-293分-求調(diào)劑 +3 | 研究生何瑤明 2026-03-18 | 3/150 |
|
|
|
[考研] 材料專業(yè)求調(diào)劑 +6 | hanamiko 2026-03-18 | 6/300 |
|
|
|
[考研] 308求調(diào)劑 +3 | 阿姐阿姐家啊 2026-03-18 | 3/150 |
|
|
|
[考研] 求調(diào)劑，一志愿:南京航空航天大學(xué)大學(xué) ，080500材料科學(xué)與工程學(xué)碩，總分289分 +4 | @taotao 2026-03-19 | 4/200 |
|
|
|
[考研] 一志愿蘇州大學(xué)材料求調(diào)劑，總分315（英一） +5 | sbdksD 2026-03-19 | 5/250 |
|
|
|
[考研] 261求B區(qū)調(diào)劑，科研經(jīng)歷豐富 +3 | 牛奶很忙 2026-03-20 | 4/200 |
|
|
|
[考研] 一志愿南理工085701環(huán)境302求調(diào)劑院校 +3 | 葵梓衛(wèi)隊(duì) 2026-03-20 | 3/150 |
|
|
|
[考研]
 320求調(diào)劑0856
+3
|
不想起名字112 2026-03-19 | 3/150 |
|
|
|
[考研] 085601材料工程專碩求調(diào)劑 +10 | 慕寒mio 2026-03-16 | 10/500 |
|
|
|
[考研] 085600材料與化工調(diào)劑 324分 +10 | llllkkkhh 2026-03-18 | 12/600 |
|
|
|
[考研] 一志愿福大288有機(jī)化學(xué)，求調(diào)劑 +3 | 小木蟲200408204 2026-03-18 | 3/150 |
|
|
|
[考研] 【同濟(jì)軟件】軟件（085405）考研求調(diào)劑 +3 | 2026eternal 2026-03-18 | 3/150 |
|
|
|
[考研] 化學(xué)工程321分求調(diào)劑 +15 | 大米飯！ 2026-03-15 | 18/900 |
|
|
|
[考研] 293求調(diào)劑 +11 | zjl的號(hào) 2026-03-16 | 16/800 |
|
|
|
[考研] 308求調(diào)劑 +4 | 是Lupa啊 2026-03-16 | 4/200 |
|
|
|
[考研] 材料工程專碩274一志愿211求調(diào)劑 +6 | 薛云鵬 2026-03-15 | 6/300 |
|