 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 3 | 1/1 | 返回列表 |
| 查看: 864 | 回復(fù): 2 | ||
charles0116金蟲 (初入文壇)
|
[求助]
高手看過來：內(nèi)部原子弛豫時(shí)自洽計(jì)算max grad (force/stress）無法滿足收斂條件tolmxf
|
|
優(yōu)化結(jié)構(gòu)之前是要先進(jìn)行內(nèi)部原子弛豫吧？ 但在做內(nèi)部原子弛豫時(shí)（optcell=0 ionmov=2 or 3) 單次scf在開始的幾輪滿足條件退出，后面就時(shí)而converged,時(shí)而not converged,這是為什么呢？ e.g. scprqt: WARNING - nstep= 200 was not enough SCF cycles to converge; potential residual= 1.475E-05 exceeds tolvrs= 1.000E-10 nstep設(shè)為200己經(jīng)差不多了吧？再調(diào)大也只是浪費(fèi)時(shí)間？ 最讓人不解的是，在前十幾輪成功converged scf計(jì)算或是后面一直not enough SCF cycles to converge 的scf計(jì)算每輪次末尾都會(huì)有 max grad (force/stress) = ××××× > tolmxf= 5.0000E-05 ha/bohr (free atoms) 類似字樣，而開始下一輪scf計(jì)算，直到in文件中設(shè)置ntime的100次scf計(jì)算完成，整個(gè)計(jì)算才退出完成。 為此，我調(diào)整了occopt=1 or 7 ;iscf=7 or 5; 還調(diào)整了各種收斂參數(shù) tolmxf; tolvrs or toldff or tolrff or tolfe ; diemac 等值，也曾嘗試過直接進(jìn)行幾何優(yōu)化（optcell=/0 ionmov=2 or 3)但都無法在ntime之內(nèi)收斂。 請(qǐng)問各位高手我這是什么問題呢？影響max grad (force/stress)值的因素有什么呢？ 我的認(rèn)為不管是內(nèi)部原子弛豫或是后來的幾何優(yōu)化，都應(yīng)該是滿足自洽計(jì)算,即 max grad (force/stress) < tolmxf的。不知我這樣的理解對(duì)嗎？ in文件內(nèi)容如下： #Structural optimisation #************************************ nshiftk 1 shiftk 0 0 0 #********************************* ecut 26 ecutsm 0.5 ngkpt 3 3 4 occopt 7 tsmear 0.01 #********************************* #Common input data chksymbreak 0 chkprim 0 brvltt -1 #Definition of the unit cell #********************************** acell 2*11.135000 4.376000 Angstroms angdeg 90 90 90 #Definition of the atoms #********************* spgroup 82 ntypat 5 znucl 48 30 16 7 6 natom 14 natrd 5 typat 1 2 3 4 5 xred 0.000000 0.000000 0.000000 0.000000 0.500000 0.250000 -0.124500 0.152800 -0.302900 -0.048000 0.359600 0.006400 -0.078300 0.273400 -0.117700 #Definition of the scf procedure #****************************** iscf 7 nstep 300 ixc 1 toldfe 5.0d-6 diemac 10.0 #Structural relaxation #************************************ ionmov 2 optcell 0 tolmxf 5.0d-5 ntime 100 #************************************ out文件我附在了帖子中 圖片是我摘出的每輪次scf結(jié)束的收斂判定條件，感覺開始max grad (force/stress)還貌似有些下降的趨勢(shì)，可是后面就毫無規(guī)律。 2013-08-20 21:06:34的屏幕截圖.png |
回復(fù)此樓» 本帖附件資源列表
-
歡迎監(jiān)督和反饋：小木蟲僅提供交流平臺(tái)，不對(duì)該內(nèi)容負(fù)責(zé)。
本內(nèi)容由用戶自主發(fā)布，如果其內(nèi)容涉及到知識(shí)產(chǎn)權(quán)問題，其責(zé)任在于用戶本人，如對(duì)版權(quán)有異議，請(qǐng)聯(lián)系郵箱：xiaomuchong@tal.com - 附件 1 : zctc.out
2013-08-20 21:35:14, 250.75 K
» 猜你喜歡
284求調(diào)劑
已經(jīng)有13人回復(fù)
-
材料科學(xué)與工程調(diào)劑
已經(jīng)有5人回復(fù)
-
346求調(diào)劑 一志愿070303有機(jī)化學(xué)
已經(jīng)有4人回復(fù)
-
329求調(diào)劑
已經(jīng)有8人回復(fù)
-
一志愿南航 335分 | 0856材料化工 | GPA 4.07 | 有科研經(jīng)歷
已經(jīng)有7人回復(fù)
-
322求調(diào)劑
已經(jīng)有8人回復(fù)
-
考研調(diào)劑
已經(jīng)有5人回復(fù)
-
307求調(diào)劑
已經(jīng)有11人回復(fù)
-
292求調(diào)劑
已經(jīng)有9人回復(fù)
-
本科211生物醫(yī)學(xué)工程085409求調(diào)劑339分
已經(jīng)有3人回復(fù)
charles0116
金蟲 (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 1062.5
- 紅花: 1
- 帖子: 43
- 在線: 112.5小時(shí)
- 蟲號(hào): 1684499
- 注冊(cè): 2012-03-12
- 專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)
2樓2013-08-21 11:27:09
charles0116
金蟲 (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 1062.5
- 紅花: 1
- 帖子: 43
- 在線: 112.5小時(shí)
- 蟲號(hào): 1684499
- 注冊(cè): 2012-03-12
- 專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)
|
根本沒有人理我，這是為什么呢？ 看到一個(gè)帖子，他說他的結(jié)構(gòu)也是無法達(dá)到tolmxf收斂。于是他改到了tolmxf=1.0d-3。（Oh,my god,我怎么就沒想到。。。）其實(shí)我之前將tolmxf調(diào)到了-4,再?zèng)]往下調(diào)。 不出我所料，我這體系果然在將tolmxf設(shè)置為5.0d-3之后： max grad (force/stress) = ××××× < tolmxf= 5.0000E-03 ha/bohr (free atoms) 完成自洽計(jì)算退出程序 可是 tolmxf=5.0-3是不是太大了？ |
3樓2013-08-23 01:59:25
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 charles0116 的主題更新
| 3 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 化學(xué)308分調(diào)劑 +7 | 你好明天你好 2026-03-30 | 8/400 |
|
|---|---|---|---|---|
|
[考研] 289求調(diào)劑 +14 | 新時(shí)代材料 2026-03-27 | 14/700 |
|
|
|
[考研] 南京大學(xué)化學(xué)調(diào)劑 +10 | 景隨風(fēng) 2026-03-29 | 15/750 |
|
|
|
[考研] 求調(diào)劑 +10 | 張zz111 2026-03-27 | 11/550 |
|
|
|
[考研] 化學(xué)0703 調(diào)劑 306分 一志愿211 +7 | 26要上岸 2026-03-28 | 7/350 |
|
|
|
[考研] 298求調(diào)劑 +3 | 種圣賜 2026-03-29 | 3/150 |
|
|
|
[考研] 343求調(diào)劑 +6 | 愛羈絆 2026-03-29 | 6/300 |
|
|
|
[考研] 317求調(diào)劑 +6 | 十閑wx 2026-03-24 | 6/300 |
|
|
|
[考研] 求調(diào)劑 +6 | 蘆lty 2026-03-25 | 7/350 |
|
|
|
[考研] 292求調(diào)劑 +14 | 鵝鵝鵝額額額額?/a> 2026-03-25 | 15/750 |
|
|
|
[考研] 張芳銘-中國(guó)農(nóng)業(yè)大學(xué)-環(huán)境工程專碩-298 +4 | 手機(jī)用戶 2026-03-26 | 4/200 |
|
|
|
[考研] 086000調(diào)劑 +3 | 7901117076 2026-03-26 | 3/150 |
|
|
|
[考研] 一志愿211院校 344分 東北農(nóng)業(yè)大學(xué)生物學(xué)學(xué)碩，求調(diào)劑 +5 | 丶風(fēng)雪夜歸人丶 2026-03-26 | 8/400 |
|
|
|
[考研] 085601 材料工程 313分 求調(diào)劑 +5 | Ong3 2026-03-27 | 5/250 |
|
|
|
[考研]
 304材料求調(diào)劑
+4
|
鐘llll 2026-03-26 | 4/200 |
|
|
|
[考研] 336材料求調(diào)劑 +7 | 陳瀅瑩 2026-03-26 | 9/450 |
|
|
|
[考研] 321求調(diào)劑 +6 | wasdssaa 2026-03-26 | 6/300 |
|
|
|
[考研] 一志愿哈工大，085400，320，求調(diào)劑 +4 | gdlf9999 2026-03-24 | 4/200 |
|
|
|
[考研] 一志愿武理085500機(jī)械專業(yè)總分300求調(diào)劑 +3 | an10101 2026-03-24 | 7/350 |
|
|
|
[考研] 一志愿北化315 求調(diào)劑 +3 | akrrain 2026-03-24 | 3/150 |
|