 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 4 | 1/1 | 返回列表 |
| 查看: 835 | 回復: 3 | ||
xiongxiong5712新蟲 (小有名氣)
|
[求助]
原子受力沒有達到收斂標準計算就結束了 已有1人參與
|
|
用siesta進行結構優(yōu)化，由于是六方晶系，首先通過擬合得到了能量最低對應的體積，MD.TypeOfRun設置為CG，然后再保持這個體積不變，優(yōu)化晶格參數(shù)，MD.TypeOfRun 設置為Broyden，具體的fdf文件如下： SystemName bulk SystemLabel bulk NumberOfAtoms 14 NumberOfSpecies 2 %block ChemicalSpeciesLabel 1 14 Si 2 7 N %endblock ChemicalSpeciesLabel LatticeConstant 1 Ang %block LatticeVectors 6.571864087507342 -3.794267500000001 0.000000000000000 0.000000000000000 7.588534999999999 0.000000000000000 0.000000000000000 0.000000000000000 2.906249000000000 %endblock LatticeVectors %block kgrid_Monkhorst_Pack 5 0 0 0 0 5 0 0 0 0 12 0 %endblock kgrid_Monkhorst_Pack AtomicCoordinatesFormat Fractional %block AtomicCoordinatesAndAtomicSpecies 0.3323000000000000 0.0314000000000000 0.2500000000000000 2 -0.0314000000000000 0.3009000000000000 0.2500000000000000 2 -0.3009000000000000 -0.3323000000000000 0.2500000000000000 2 -0.3323000000000000 -0.0314000000000000 0.7500000000000000 2 0.0314000000000000 -0.3009000000000000 0.7500000000000000 2 0.3009000000000000 0.3323000000000000 0.7500000000000000 2 0.3333333333333330 0.6666666666666670 0.2500000000000000 2 -0.3333333333333330 -0.6666666666666670 0.7500000000000000 2 0.1733000000000000 -0.2306000000000000 0.2500000000000000 1 0.2306000000000000 0.4039000000000000 0.2500000000000000 1 -0.4039000000000000 -0.1733000000000000 0.2500000000000000 1 -0.1733000000000000 0.2306000000000000 0.7500000000000000 1 -0.2306000000000000 -0.4039000000000000 0.7500000000000000 1 0.4039000000000000 0.1733000000000000 0.7500000000000000 1 %endblock AtomicCoordinatesAndAtomicSpecies xc.functional LDA xc.authors CA MeshCutoff 100 Ry PAO.BasisSize DZP PAO.EnergyShift 75 mev MD.MaxForceTol 0.05 eV/Ang DM.MixingWeight 0.01 DM.NumberPulay 10 MaxSCFIterations 200 SolutionMethod diagon MD.TypeOfRun Broyden MD.NumCGsteps 500 MD.VariableCell .true. MD.ConstantVolume .true. MD.RelaxCellOnly .true. 優(yōu)化很快， Begin CG move = 1的時候就結束了，dMax達到了收斂標準，但是原子受力沒有收斂，結果是： siesta: Atomic forces (eV/Ang): ---------------------------------------- Tot -0.000014 -0.000041 -0.000144 ---------------------------------------- Max 0.722244 Res 0.314921 sqrt( Sum f_i^2 / 3N ) ---------------------------------------- Max 0.722244 constrained 可見原子受力沒有達到0.05以下，而計算就正常結束了。請教大家，這是什么原因引起的�。渴遣皇俏襇D.TypeOfRun 設置為Broyden的原因呢？用這種方法優(yōu)化的時候應該注意什么？有勞各位幫幫忙，多謝了。 [ Last edited by xiongxiong5712 on 2013-5-31 at 21:46 ] |
回復此樓» 猜你喜歡
323求調劑
已經有6人回復
-
一志愿北京化工大學 070300 學碩 336分 求調劑
已經有4人回復
-
352求調劑
已經有3人回復
-
一志愿東華大學化學070300，求調劑
已經有8人回復
-
277材料科學與工程080500求調劑
已經有7人回復
-
317求調劑
已經有18人回復
-
293求調劑
已經有5人回復
-
280分求調劑 一志愿085802
已經有7人回復
-
0854電子信息求調劑
已經有3人回復
-
263求調劑
已經有4人回復
» 本主題相關價值貼推薦，對您同樣有幫助:
- 燃燒計算，先收斂后發(fā)散
已經有12人回復
- 關于回流與計算結果收斂的討論
已經有11人回復
- ansys接觸計算不收斂問題
已經有16人回復
- pwscf計算integrated charge出錯及不收斂
已經有13人回復
- 怎樣判定fluent計算是否收斂
已經有4人回復
- NEB計算擴散路徑的力收斂精度設置問題
已經有16人回復
- Castep計算無法收斂，請高人指教
已經有17人回復
- 加溶劑后計算，不收斂，請高手指點
已經有5人回復
- CASTEP計算摻雜不收斂
已經有5人回復
- 計算時候總是不能收斂如何解決？
已經有5人回復
- 靜態(tài)計算能收斂，同樣的參數(shù)弛豫不能收斂，這是為什么？
已經有4人回復
- ELK計算嚴重不收斂～
已經有10人回復
- vasp計算：優(yōu)化能收斂，但同樣的設置靜態(tài)計算不收斂
已經有3人回復
- IRC計算沒收斂
已經有9人回復
- DMol3計算頻率不收斂
已經有4人回復
- 【求助】用Orca計算一個片簇的單點能，不收斂。
已經有7人回復
- 【求助】計算不收斂怎么辦？
已經有4人回復
- 【求助】請教PWSCF晶體建模和自洽計算，不收斂
已經有21人回復
xiongxiong5712
新蟲 (小有名氣)
- 應助: 1 (幼兒園)
- 金幣: 543
- 散金: 12
- 紅花: 2
- 帖子: 257
- 在線: 267.7小時
- 蟲號: 706467
- 注冊: 2009-02-23
- 專業(yè): 半導體材料
|
自己頂一下。麻煩各位大俠百忙中幫我看看，或者提示一下也行。我先說說我的看法：是不是因為MD.TypeOfRun 不是CG，計算不會進行CGmove的循環(huán)，所以在Begin CG move = 1的時候就結束了，是否收斂也不是看force里的Max，而是有一個地方給出了每個原子的受力： siesta: Atomic forces (eV/Ang): siesta: 1 -0.368915 -0.160954 -0.000062 siesta: 2 0.323778 -0.238986 -0.000043 siesta: 3 0.045044 0.399972 0.000016 siesta: 4 0.368938 0.160918 0.000009 siesta: 5 -0.323802 0.238974 -0.000020 siesta: 6 -0.045044 -0.399962 -0.000040 siesta: 7 0.000000 -0.000020 -0.000041 siesta: 8 -0.000020 0.000009 0.000023 siesta: 9 0.450554 -0.573922 -0.000001 siesta: 10 0.271759 0.677123 0.000003 siesta: 11 -0.722307 -0.103215 -0.000009 siesta: 12 -0.450542 0.573917 0.000002 siesta: 13 -0.271775 -0.677130 0.000005 siesta: 14 0.722304 0.103217 0.000011 siesta: ---------------------------------------- siesta: Tot -0.000027 -0.000058 -0.000148 最后一行是總的受力，是不是看這個地方收斂即可？（CG方法沒有給出這一結果）。不知道我猜測的對不對，請大家指正。 |
2樓2013-06-01 15:10:48
newton3915
鐵桿木蟲 (正式寫手)
- 應助: 4 (幼兒園)
- 金幣: 4742.3
- 散金: 40
- 紅花: 4
- 帖子: 418
- 在線: 391.3小時
- 蟲號: 609189
- 注冊: 2008-09-22
- 專業(yè): 凝聚態(tài)物性 II ：電子結構
3樓2014-12-15 04:22:57
4樓2016-05-10 22:10:45
相關版塊跳轉
我要訂閱樓主 xiongxiong5712 的主題更新
| 4 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 0854電子信息求調劑 +3 | α____ 2026-03-22 | 3/150 |
|
|---|---|---|---|---|
|
|
[考研] 一志愿中南化學（0703）總分337求調劑 +9 | niko- 2026-03-19 | 10/500 |
|
|
|
[考研] 材料學碩301分求調劑 +7 | Liyouyumairs 2026-03-21 | 7/350 |
|
|
[考研] 初試 317 +7 | 半拉月丙 2026-03-20 | 7/350 |
|
|
|
[考研] 333求調劑 +5 | 87639 2026-03-21 | 7/350 |
|
|
|
[考研]
 一志愿東華大學控制學碩320求調劑
+3
|
Grand777 2026-03-21 | 3/150 |
|
|
|
[考研] 求調劑 +3 | .m.. 2026-03-21 | 4/200 |
|
|
|
[考研] 22 350 本科985求調劑，求老登收留 +3 | 李軼男003 2026-03-20 | 3/150 |
|
|
|
[考研] 346求調劑[0856] +4 | WayneLim327 2026-03-16 | 7/350 |
|
|
|
[考研] 332求調劑 +4 | ydfyh 2026-03-17 | 4/200 |
|
|
|
[考研] 085700資源與環(huán)境308求調劑 +12 | 墨墨漠 2026-03-18 | 13/650 |
|
|
|
[考研] 材料專碩英一數(shù)二306 +7 | z1z2z3879 2026-03-18 | 7/350 |
|
|
|
[考研] 321求調劑 +9 | 何潤采123 2026-03-18 | 11/550 |
|
|
|
[考研]
材料與化工 322求調劑
+4
|
然11 2026-03-19 | 4/200 |
|
|
|
[考研]
 295復試調劑
+8
|
簡木ChuFront 2026-03-19 | 8/400 |
|
|
|
[考研] 材料學求調劑 +4 | Stella_Yao 2026-03-20 | 4/200 |
|
|
|
[考研] 一志愿西安交通大學 學碩 354求調劑211或者雙一流 +3 | 我想要讀研究生 2026-03-20 | 3/150 |
|
|
|
[考研] 一志愿福大288有機化學，求調劑 +3 | 小木蟲200408204 2026-03-18 | 3/150 |
|
|
|
[考研] 0854可跨調劑，一作一項核心論文五項專利，省、國級證書40+數(shù)一英一287 +8 | 小李0854 2026-03-16 | 8/400 |
|
|
|
[考研] 333求調劑 +3 | 文思客 2026-03-16 | 7/350 |
|