 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 839 | 回復: 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 ] |
回復此樓» 猜你喜歡
一志愿070300浙大化學358分，求調劑！
已經有4人回復
-
一志愿中南大學化學學碩0703總分337求調劑
已經有4人回復
-
一志愿北京化工大學 070300 學碩 336分 求調劑
已經有5人回復
-
323求調劑
已經有6人回復
-
352求調劑
已經有3人回復
-
一志愿東華大學化學070300，求調劑
已經有8人回復
-
277材料科學與工程080500求調劑
已經有7人回復
-
317求調劑
已經有18人回復
-
293求調劑
已經有5人回復
-
280分求調劑 一志愿085802
已經有7人回復
» 本主題相關價值貼推薦，對您同樣有幫助:
- 燃燒計算，先收斂后發(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人回復
4樓2016-05-10 22:10:45
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
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 293求調劑 +3 | 濤濤Wjt 2026-03-22 | 5/250 |
|
|---|---|---|---|---|
|
|
[考研] 298求調劑 +6 | 上岸6666@ 2026-03-20 | 6/300 |
|
|
|
[考研] 311求調劑 +6 | 冬十三 2026-03-18 | 6/300 |
|
|
[考研] 291 求調劑 +3 | 化工2026屆畢業(yè)?/a> 2026-03-21 | 3/150 |
|
|
|
[考研] 318求調劑 +4 | plum李子 2026-03-21 | 7/350 |
|
|
|
[考研] 廣西大學材料導師推薦 +3 | 夏夏夏小正 2026-03-17 | 5/250 |
|
|
|
[考研] 297求調劑 +3 | 喜歡還是不甘心 2026-03-20 | 3/150 |
|
|
|
[考研] 一志愿深大，0703化學，總分302，求調劑 +4 | 七月-七七 2026-03-21 | 4/200 |
|
|
|
[考研] 311求調劑 +3 | 勇敢的小吳 2026-03-20 | 3/150 |
|
|
|
[考研] 299求調劑 +4 | 某某某某位 2026-03-21 | 4/200 |
|
|
|
[考研] 268求調劑 +9 | 簡單點0 2026-03-17 | 9/450 |
|
|
|
[考研] 材料專業(yè)求調劑 +6 | hanamiko 2026-03-18 | 6/300 |
|
|
|
[考研] 274求調劑 +10 | S.H1 2026-03-18 | 10/500 |
|
|
|
[考研] 22408 344分 求調劑 一志愿 華電計算機技術 +4 | solanXXX 2026-03-20 | 4/200 |
|
|
|
[考研] 中南大學化學學碩337求調劑 +3 | niko- 2026-03-19 | 6/300 |
|
|
|
[考研] 材料與化工專碩調劑 +7 | heming3743 2026-03-16 | 7/350 |
|
|
|
[碩博家園] 湖北工業(yè)大學 生命科學與健康學院-課題組招收2026級食品/生物方向碩士 +3 | 1喜春8 2026-03-17 | 5/250 |
|
|
|
[論文投稿] 有沒有大佬發(fā)小論文能帶我個二作 +3 | 增銳漏人 2026-03-17 | 4/200 |
|
|
|
[考研]
 11408 一志愿西電，277分求調劑
+3
|
zhouzhen654 2026-03-16 | 3/150 |
|
|
|
[考研] 070300化學學碩求調劑 +6 | 太想進步了0608 2026-03-16 | 6/300 |
|