 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 4 | 1/1 | 返回列表 |
| 查看: 841 | 回復: 3 | |||
xiongxiong5712新蟲 (小有名氣)
|
[求助]
原子受力沒有達到收斂標準計算就結(jié)束了 已有1人參與
|
|
用siesta進行結(jié)構(gòu)優(yōu)化，由于是六方晶系，首先通過擬合得到了能量最低對應(yīng)的體積，MD.TypeOfRun設(shè)置為CG，然后再保持這個體積不變，優(yōu)化晶格參數(shù)，MD.TypeOfRun 設(shè)置為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的時候就結(jié)束了，dMax達到了收斂標準，但是原子受力沒有收斂，結(jié)果是： 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以下，而計算就正常結(jié)束了。請教大家，這是什么原因引起的��？是不是我MD.TypeOfRun 設(shè)置為Broyden的原因呢？用這種方法優(yōu)化的時候應(yīng)該注意什么？有勞各位幫幫忙，多謝了。 [ Last edited by xiongxiong5712 on 2013-5-31 at 21:46 ] |
回復此樓» 猜你喜歡
招08考數(shù)學
已經(jīng)有11人回復
-
291求調(diào)劑
已經(jīng)有4人回復
-
276求調(diào)劑
已經(jīng)有4人回復
-
一志愿070300浙大化學358分，求調(diào)劑！
已經(jīng)有4人回復
-
一志愿中南大學化學學碩0703總分337求調(diào)劑
已經(jīng)有4人回復
-
一志愿北京化工大學 070300 學碩 336分 求調(diào)劑
已經(jīng)有5人回復
-
323求調(diào)劑
已經(jīng)有6人回復
-
352求調(diào)劑
已經(jīng)有3人回復
-
一志愿東華大學化學070300，求調(diào)劑
已經(jīng)有8人回復
-
277材料科學與工程080500求調(diào)劑
已經(jīng)有7人回復
» 本主題相關(guān)價值貼推薦，對您同樣有幫助:
- 燃燒計算，先收斂后發(fā)散
已經(jīng)有12人回復
- 關(guān)于回流與計算結(jié)果收斂的討論
已經(jīng)有11人回復
- ansys接觸計算不收斂問題
已經(jīng)有16人回復
- pwscf計算integrated charge出錯及不收斂
已經(jīng)有13人回復
- 怎樣判定fluent計算是否收斂
已經(jīng)有4人回復
- NEB計算擴散路徑的力收斂精度設(shè)置問題
已經(jīng)有16人回復
- Castep計算無法收斂，請高人指教
已經(jīng)有17人回復
- 加溶劑后計算，不收斂，請高手指點
已經(jīng)有5人回復
- CASTEP計算摻雜不收斂
已經(jīng)有5人回復
- 計算時候總是不能收斂如何解決？
已經(jīng)有5人回復
- 靜態(tài)計算能收斂，同樣的參數(shù)弛豫不能收斂，這是為什么？
已經(jīng)有4人回復
- ELK計算嚴重不收斂～
已經(jīng)有10人回復
- vasp計算：優(yōu)化能收斂，但同樣的設(shè)置靜態(tài)計算不收斂
已經(jīng)有3人回復
- IRC計算沒收斂
已經(jīng)有9人回復
- DMol3計算頻率不收斂
已經(jīng)有4人回復
- 【求助】用Orca計算一個片簇的單點能，不收斂。
已經(jīng)有7人回復
- 【求助】計算不收斂怎么辦？
已經(jīng)有4人回復
- 【求助】請教PWSCF晶體建模和自洽計算，不收斂
已經(jīng)有21人回復
xiongxiong5712
新蟲 (小有名氣)
- 應(yīng)助: 1 (幼兒園)
- 金幣: 543
- 散金: 12
- 紅花: 2
- 帖子: 257
- 在線: 267.7小時
- 蟲號: 706467
- 注冊: 2009-02-23
- 專業(yè): 半導體材料
|
自己頂一下。麻煩各位大俠百忙中幫我看看，或者提示一下也行。我先說說我的看法：是不是因為MD.TypeOfRun 不是CG，計算不會進行CGmove的循環(huán)，所以在Begin CG move = 1的時候就結(jié)束了，是否收斂也不是看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方法沒有給出這一結(jié)果）。不知道我猜測的對不對，請大家指正。 |
2樓2013-06-01 15:10:48
newton3915
鐵桿木蟲 (正式寫手)
- 應(yīng)助: 4 (幼兒園)
- 金幣: 4742.3
- 散金: 40
- 紅花: 4
- 帖子: 418
- 在線: 391.3小時
- 蟲號: 609189
- 注冊: 2008-09-22
- 專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)
3樓2014-12-15 04:22:57
4樓2016-05-10 22:10:45
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 xiongxiong5712 的主題更新
| 4 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 323求調(diào)劑 +6 | 洼小桶 2026-03-18 | 6/300 |
|
|---|---|---|---|---|
|
|
[考研] 277材料科學與工程080500求調(diào)劑 +7 | 自由煎餅果子 2026-03-16 | 7/350 |
|
|
|
[考研] 310求調(diào)劑 +4 | baibai1314 2026-03-16 | 4/200 |
|
|
[考研] 求調(diào)劑一志愿海大，0703化學學碩304分，有大創(chuàng)項目，四級已過 +6 | 幸運哩哩 2026-03-22 | 10/500 |
|
|
|
[考研] 275求調(diào)劑 +6 | shansx 2026-03-22 | 8/400 |
|
|
|
[考研] 材料求調(diào)劑 +5 | @taotao 2026-03-21 | 5/250 |
|
|
|
[考研] 一志愿南大，0703化學，分數(shù)336，求調(diào)劑 +3 | 收到VS 2026-03-21 | 3/150 |
|
|
|
[基金申請]
學校已經(jīng)提交到NSFC，還能修改嗎？
 40+4
|
babangida 2026-03-19 | 9/450 |
|
|
|
[考研] 316求調(diào)劑 +6 | 梁茜雯 2026-03-19 | 6/300 |
|
|
|
[考研] 301求調(diào)劑 +10 | yy要上岸呀 2026-03-17 | 10/500 |
|
|
|
[考研] 求調(diào)劑 +3 | Ma_xt 2026-03-17 | 3/150 |
|
|
|
[考研] 南京大學化學376求調(diào)劑 +3 | hisfailed 2026-03-19 | 6/300 |
|
|
|
[考研] 一志愿中海洋材料工程專碩330分求調(diào)劑 +8 | 小材化本科 2026-03-18 | 8/400 |
|
|
|
[考研] 一志愿 西北大學 ，070300化學學碩，總分287，雙非一本，求調(diào)劑。 +4 | 晨昏線與星海 2026-03-19 | 4/200 |
|
|
|
[考研] 中南大學化學學碩337求調(diào)劑 +3 | niko- 2026-03-19 | 6/300 |
|
|
|
[考研] 求調(diào)劑一志愿南京航空航天大學289分 +3 | @taotao 2026-03-19 | 3/150 |
|
|
|
[考研] A區(qū)線材料學調(diào)劑 +5 | 周周無極 2026-03-20 | 5/250 |
|
|
|
[考研] 085601材料工程專碩求調(diào)劑 +10 | 慕寒mio 2026-03-16 | 10/500 |
|
|
|
[考研] 收復試調(diào)劑生 +4 | 雨后秋荷 2026-03-18 | 4/200 |
|
|
|
[考研] 312求調(diào)劑 +8 | 陌宸希 2026-03-16 | 9/450 |
|