 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 3 | 1/1 | 返回列表 |
| 查看: 865 | 回復(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ì)算每輪次末尾都會有 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)收斂。 請問各位高手我這是什么問題呢？影響max grad (force/stress)值的因素有什么呢？ 我的認(rèn)為不管是內(nèi)部原子弛豫或是后來的幾何優(yōu)化，都應(yīng)該是滿足自洽計(jì)算,即 max grad (force/stress) < tolmxf的。不知我這樣的理解對嗎？ 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)還貌似有些下降的趨勢，可是后面就毫無規(guī)律。 2013-08-20 21:06:34的屏幕截圖.png |
回復(fù)此樓» 本帖附件資源列表
-
歡迎監(jiān)督和反饋：小木蟲僅提供交流平臺，不對該內(nèi)容負(fù)責(zé)。
本內(nèi)容由用戶自主發(fā)布，如果其內(nèi)容涉及到知識產(chǎn)權(quán)問題，其責(zé)任在于用戶本人，如對版權(quán)有異議，請聯(lián)系郵箱：xiaomuchong@tal.com - 附件 1 : zctc.out
2013-08-20 21:35:14, 250.75 K
» 猜你喜歡
材料科學(xué)與工程求調(diào)劑
已經(jīng)有8人回復(fù)
-
一志愿大連理工大學(xué)材料求調(diào)劑
已經(jīng)有3人回復(fù)
-
一志愿北化085600材料專碩275|有文章專利｜求調(diào)劑
已經(jīng)有9人回復(fù)
-
材料求調(diào)劑 一志愿哈工大總分298分，前三科223分
已經(jīng)有8人回復(fù)
-
274求調(diào)劑
已經(jīng)有3人回復(fù)
-
合肥區(qū)域性重點(diǎn)一本招收調(diào)劑
已經(jīng)有4人回復(fù)
-
085600材料與化工調(diào)劑
已經(jīng)有10人回復(fù)
-
22408 359分調(diào)劑
已經(jīng)有7人回復(fù)
-
抱歉
已經(jīng)有4人回復(fù)
-
085602化工求調(diào)劑（331分）
已經(jīng)有8人回復(fù)
charles0116
金蟲 (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 1062.5
- 紅花: 1
- 帖子: 43
- 在線: 112.5小時(shí)
- 蟲號: 1684499
- 注冊: 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í)
- 蟲號: 1684499
- 注冊: 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,再沒往下調(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ā)表 | |
|---|---|---|---|---|
|
[考研] 085701環(huán)境工程求調(diào)劑 +11 | 多久上課 2026-03-27 | 12/600 |
|
|---|---|---|---|---|
|
|
[考研] 309求調(diào)劑 +15 | 誰不是少年 2026-03-29 | 15/750 |
|
|
[考研] 295材料工程專碩求調(diào)劑 +10 | 1428151015 2026-03-27 | 10/500 |
|
|
|
[考研] 0703 化學(xué) 求調(diào)劑，一志愿山東大學(xué) 342 分 +7 | Shern—- 2026-03-28 | 7/350 |
|
|
|
[考研] 材料科學(xué)與工程調(diào)劑 +8 | 深V宿舍吧 2026-03-30 | 8/400 |
|
|
|
[考研] 298求調(diào)劑 +3 | 什么是胖頭魚 2026-03-30 | 5/250 |
|
|
|
[考研] 南京大學(xué)化學(xué)調(diào)劑 +10 | 景隨風(fēng) 2026-03-29 | 15/750 |
|
|
|
[考研] 求調(diào)劑 +10 | 張zz111 2026-03-27 | 11/550 |
|
|
|
[考研] 一志愿北京工業(yè)大學(xué)，324分求調(diào)劑 +6 | 零八# 2026-03-28 | 6/300 |
|
|
|
[考研] 求調(diào)劑 +4 | QiMing7 2026-03-25 | 5/250 |
|
|
|
[考研] 327求調(diào)劑 +6 | 汲亦昊 2026-03-29 | 6/300 |
|
|
|
[考研] 349求調(diào)劑 +6 | 李木子啊哈哈 2026-03-25 | 6/300 |
|
|
|
[考研] 298求調(diào)劑 +4 | 種圣賜 2026-03-28 | 4/200 |
|
|
|
[考研] 本科新能源科學(xué)與工程，一志愿華理能動285求調(diào)劑 +7 | AZMK 2026-03-28 | 11/550 |
|
|
|
[考研] 311（085601）求調(diào)劑 +4 | liziyeyeye 2026-03-28 | 4/200 |
|
|
|
[考研] 材料277求調(diào)劑 +7 | min3 2026-03-24 | 7/350 |
|
|
|
[考研] 292求調(diào)劑 +14 | 鵝鵝鵝額額額額?/a> 2026-03-25 | 15/750 |
|
|
|
[考研] 272求調(diào)劑 +7 | 腳滑的守法公民 2026-03-27 | 7/350 |
|
|
|
[考研] 266分求材料化工冶金礦業(yè)等專業(yè)的調(diào)劑 +4 | 哇呼哼呼哼 2026-03-26 | 4/200 |
|
|
|
[考研] 308求調(diào)劑 +7 | 墨墨漠 2026-03-25 | 7/350 |
|