 |
 |
 |
| 24小時(shí)熱門(mén)版塊排行榜 |
- >論壇更新日志 (989)
- >考研 (185)
- >導(dǎo)師招生 (132)
- >基金申請(qǐng) (19)
- >考博 (14)
- >蟲(chóng)友互識(shí) (10)
- >碩博家園 (8)
- >休閑灌水 (7)
- >博后之家 (5)
- >招聘信息布告欄 (4)
- >論文道賀祈福 (4)
- >復(fù)合材料 (2)
- >學(xué)術(shù)會(huì)議 (2)
- >教師之家 (2)
- >找工作 (2)
- >專(zhuān)業(yè)外語(yǔ) (1)
 |
 |
 |
||
 |
 |
 |
 |
 |
| 6 | 1/1 | 返回列表 |
| 查看: 1119 | 回復(fù): 5 | ||||
[求助]
菜鳥(niǎo)求助關(guān)于加壓優(yōu)化的問(wèn)題。。。
|
|
最近在做關(guān)于加壓優(yōu)化的東西， 第一個(gè)問(wèn)題：我在toturial 中沒(méi)找見(jiàn)具體的加壓優(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ù)此樓» 猜你喜歡
070300化學(xué)319求調(diào)劑
已經(jīng)有7人回復(fù)
-
08工科 320總分 求調(diào)劑
已經(jīng)有6人回復(fù)
-
一志愿天津大學(xué)化學(xué)工藝專(zhuān)業(yè)（081702）315分求調(diào)劑
已經(jīng)有12人回復(fù)
-
307求調(diào)劑
已經(jīng)有3人回復(fù)
-
265求調(diào)劑
已經(jīng)有3人回復(fù)
-
301求調(diào)劑
已經(jīng)有10人回復(fù)
-
303求調(diào)劑
已經(jīng)有7人回復(fù)
-
材料工程（專(zhuān)）一志愿985 初試335求調(diào)劑
已經(jīng)有4人回復(fù)
-
299求調(diào)劑
已經(jīng)有6人回復(fù)
-
083200學(xué)碩321分一志愿暨南大學(xué)求調(diào)劑
已經(jīng)有3人回復(fù)
souledge
專(zhuān)家顧問(wèn) (著名寫(xiě)手)
-
專(zhuān)家經(jīng)驗(yàn): +174 - 1ST強(qiáng)帖: 19
- 應(yīng)助: 303 (大學(xué)生)
- 金幣: 11306.4
- 散金: 1123
- 紅花: 108
- 帖子: 2814
- 在線: 916.1小時(shí)
- 蟲(chóng)號(hào): 941486
- 注冊(cè): 2010-01-12
- 專(zhuān)業(yè): 結(jié)構(gòu)陶瓷
- 管轄: 第一性原理
【答案】應(yīng)助回帖
★ ★ ★ ★ ★
感謝參與，應(yīng)助指數(shù) +1
fzx2008: 金幣+5, 專(zhuān)家考核, 謝謝指教！ 2012-12-04 18:38:56
感謝參與，應(yīng)助指數(shù) +1
fzx2008: 金幣+5, 專(zhuān)家考核, 謝謝指教！ 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非等效軸，所以難以?xún)?yōu)化后完全一致，但是該問(wèn)題顯然在對(duì)稱(chēng)性極高的立方晶上不存在。其實(shí)可以看到，如果優(yōu)化的收斂條件很高，33方向的壓強(qiáng)與11，22的差別很小。至于加壓，一般都是加等靜壓，當(dāng)然是11=22=33了。 |
» 本帖已獲得的紅花（最新10朵）
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)化！）也就是說(shuō)我應(yīng)該看壓力張量分量呢還是以整體殘余壓力為標(biāo)準(zhǔn)？ 最后，我懸賞的10個(gè)金幣應(yīng)該怎樣支付給你？給你兩朵小紅花嗎？ |
4樓2012-12-05 10:30:33
souledge
專(zhuān)家顧問(wèn) (著名寫(xiě)手)
-
專(zhuān)家經(jīng)驗(yàn): +174 - 1ST強(qiáng)帖: 19
- 應(yīng)助: 303 (大學(xué)生)
- 金幣: 11306.4
- 散金: 1123
- 紅花: 108
- 帖子: 2814
- 在線: 916.1小時(shí)
- 蟲(chóng)號(hào): 941486
- 注冊(cè): 2010-01-12
- 專(zhuān)業(yè): 結(jié)構(gòu)陶瓷
- 管轄: 第一性原理
【答案】應(yīng)助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★
若博特: 金幣+5, ★★★很有幫助, 很熱心的專(zhuān)家 2012-12-14 09:03:28
若博特: 金幣+5, ★★★很有幫助, 很好的寫(xiě)手，謝謝了，希望能及時(shí)回復(fù)下邊的問(wèn)題。 2012-12-15 16:20:41
若博特: 金幣+5, ★★★很有幫助, 很熱心的專(zhuān)家 2012-12-14 09:03:28
若博特: 金幣+5, ★★★很有幫助, 很好的寫(xiě)手，謝謝了，希望能及時(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)嗄苤?lè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ā)表 | |
|---|---|---|---|---|
|
[考研]
 271材料工程求調(diào)劑
+8
|
.6lL 2026-03-18 | 8/400 |
|
|---|---|---|---|---|
|
[考研] 317求調(diào)劑 +8 | 申子申申 2026-03-19 | 13/650 |
|
|
|
[考研] 296求調(diào)劑 +6 | www_q 2026-03-18 | 10/500 |
|
|
|
[考研] 288求調(diào)劑 +16 | 于海海海海 2026-03-19 | 16/800 |
|
|
|
[考研] 317求調(diào)劑 +5 | 申子申申 2026-03-19 | 9/450 |
|
|
|
[考研] 一志愿西南交通 專(zhuān)碩 材料355 本科雙非 求調(diào)劑 +5 | 西南交通專(zhuān)材355 2026-03-19 | 5/250 |
|
|
|
[考研]
 295復(fù)試調(diào)劑
+8
|
簡(jiǎn)木ChuFront 2026-03-19 | 8/400 |
|
|
|
[考研] 一志愿西安交通大學(xué) 學(xué)碩 354求調(diào)劑 +3 | 我想要讀研究生 2026-03-20 | 3/150 |
|
|
|
[考研] 招收調(diào)劑碩士 +4 | lidianxing 2026-03-19 | 12/600 |
|
|
|
[考研] 能源材料化學(xué)課題組招收碩士研究生8-10名 +5 | 脫穎而出 2026-03-16 | 14/700 |
|
|
|
[考研] 生物學(xué)調(diào)劑招人！�。� +3 | 山海天嵐 2026-03-17 | 4/200 |
|
|
|
[考研] 本科鄭州大學(xué)物理學(xué)院，一志愿華科070200學(xué)碩，346求調(diào)劑 +4 | 我不是一根蔥 2026-03-18 | 4/200 |
|
|
|
[考研] 【同濟(jì)軟件】軟件（085405）考研求調(diào)劑 +3 | 2026eternal 2026-03-18 | 3/150 |
|
|
|
[考研] 化學(xué)工程321分求調(diào)劑 +15 | 大米飯！ 2026-03-15 | 18/900 |
|
|
|
[考研] 生物學(xué)071000 329分求調(diào)劑 +3 | 我愛(ài)生物生物愛(ài)?/a> 2026-03-17 | 3/150 |
|
|
|
[考研] 290求調(diào)劑 +3 | p asserby. 2026-03-15 | 4/200 |
|
|
|
[考博] 26申博 +4 | 八6八68 2026-03-16 | 4/200 |
|
|
|
[考研]
11408 一志愿西電，277分求調(diào)劑
+3
|
zhouzhen654 2026-03-16 | 3/150 |
|
|
|
[考研] 機(jī)械專(zhuān)碩325，尋找調(diào)劑院校 +3 | y9999 2026-03-15 | 5/250 |
|
|
|
[考研] 0854控制工程 359求調(diào)劑 可跨專(zhuān)業(yè) +3 | 626776879 2026-03-14 | 9/450 |
|