 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 1201 | 回復: 6 | |||
| 當前只顯示滿足指定條件的回帖，點擊這里查看本話題的所有回帖 | |||
[求助]
請教，為什么會出現(xiàn)這樣的情況？謝謝大家
|
|||
|
我在測試mesh cutoff的時候，一直很順利，mesh cutoff從200Ry一直增加到750ry，200到700ry的測試中都很順利完成了scf運算，但是當mesh cutoff增加到750ry時，出現(xiàn)了這樣的輸出情況 InitMesh: MESH = 80 x 80 x 80 = 512000 InitMesh: Mesh cutoff (required, used) = 750.000 870.761 Ry * Maximum dynamic memory allocated = 168 MB stepf: Fermi-Dirac step function siesta: Program's energy decomposition (eV): siesta: Ebs = -48733249.309373 siesta: Eions = 10539.985266 siesta: Ena =***************** siesta: Ekin = 5039.821056 siesta: Enl = -2320.791946 siesta: DEna = -0.148201 siesta: DUscf = 0.000000 siesta: DUext = 0.000000 siesta: Exc = -1407.884017 siesta: eta*DQ = 0.000000 siesta: Emadel = 0.000000 siesta: Emeta = 0.000000 siesta: Emolmec = 0.000000 siesta: Ekinion = 0.000000 siesta: Eharris =***************** siesta: Etot =***************** siesta: FreeEng =***************** siesta: iscf Eharris(eV) E_KS(eV) FreeEng(eV) dDmax Ef(eV) siesta: 1******************************************390.3300******** timer: Routine,Calls,Time,% = IterSCF 1 2189.282 98.32 elaps: Routine,Calls,Wall,% = IterSCF 1 327.087 98.11 siesta: 2****************************************** 0.0553******** siesta: 3****************************************** 0.0453******** siesta: 4****************************************** 0.0027******** siesta: 5****************************************** 0.0030******** siesta: 6****************************************** 0.0009******** siesta: 7****************************************** 0.0007******** siesta: 8****************************************** 0.0003******** siesta: 9****************************************** 0.0002******** siesta: 10****************************************** 0.0002******** siesta: 11****************************************** 0.0002******** siesta: 12****************************************** 0.0002******** siesta: 13****************************************** 0.0002******** siesta: 14****************************************** 0.0001******** siesta: 15****************************************** 0.0001******** siesta: 16****************************************** 0.0002******** siesta: 17****************************************** 0.0002******** siesta: 18****************************************** 0.0002******** siesta: 19****************************************** 0.0002******** siesta: 20****************************************** 0.0002******** siesta: 21****************************************** 0.0001******** 請問是為什么呢？ 我的輸入文件如下： SystemName SystemLabel NumberOfAtoms 10 NumberOfSpecies 4 %block ChemicalSpeciesLabel 1 * * # Species index, atomic number, species label 2 * * 3 * * 4 * * %endblock ChemicalSpeciesLabel LatticeConstant 1.889727 Bohr %block LatticeVectors * * * * * * * * * %endblock LatticeVectors AtomicCoordinatesFormat ScaledByLatticeVectors %block AtomicCoordinatesAndAtomicSpecies * * * 1 * * * 1 * * * 2 * * * 3 * * * 4 * * * 4 * * * 0 4 * * * 4 * * * 4 * * * 4 %endblock AtomicCoordinatesAndAtomicSpecies %block AtomicCoordinatesOrigin 0.000 0.000 0.000 %endblock AtomicCoordinatesOrigin PAO.BasisType split PAO.SplitNorm 0.15 PAO.EnergyShift 0.01 Ry %block PAO.Basis # Define Basis set * 5 2.23376 n=5 0 1 E 75.06727 3.90515 4.47365 1.00000 n=6 0 2 E 71.19116 3.69937 5.35635 4.45480 1.00000 1.00000 n=5 1 1 E 9.59098 0.01349 4.68524 1.00000 n=6 1 1 E 87.28316 3.65766 5.49527 1.00000 n=5 2 1 E 14.20169 4.25960 5.28396 1.00000 * 5 2.25014 n=3 0 1 E 0.64399 2.65247 3.03150 1.00000 n=4 0 2 E 18.52214 3.13305 4.15184 3.42190 1.00000 1.00000 n=3 1 1 E 8.12894 2.27894 3.02759 1.00000 n=4 1 1 E 16.26911 2.00889 4.26349 1.00000 n=3 2 2 E 14.23042 3.38062 3.78964 2.62819 1.00000 1.00000 * 3 0.09 n=2 0 2 E 24.13 2.59 4.601 2.082 1.000 1.000 n=2 1 2 E 47.15 2.61 5.813 2.127 1.000 1.000 n=3 2 1 E 41.49 4.31 7.798 1.000 %endblock PAO.Basis %block ProjectedDensityOfStates -20.00 10.00 0.050 2000 eV %endblock ProjectedDensityOfStates MeshCutoff 750. Ry FilterCutoff 750. Ry ElectronicTemperature 0.001 eV DM.MixingWeight 0.1 DM.NumberPulay 4 MaxSCFiterations 400 SCFMustConverge .true. WriteMullikenPop 1 COOP.Write .true. SpinPolarized yes NonCollinearSpin no %block kgrid_Monkhorst_Pack 6 0 0 0.5 0 6 0 0.5 0 0 6 0.5 %endblock kgrid_Monkhorst_Pack DM.UseSaveDM yes %block DM.InitSpin * +* 0. 0. # atom,spin,theta,phi (deg) * +* 0. 0. %endblock DM.InitSpin 之前的計算中并無問題，請問怎么回事呢？謝謝 |
回復此樓» 猜你喜歡
0703化學336分求調(diào)劑
已經(jīng)有6人回復
-
304求調(diào)劑
已經(jīng)有9人回復
-
293求調(diào)劑
已經(jīng)有13人回復
-
281求調(diào)劑（0805）
已經(jīng)有8人回復
-
環(huán)境領(lǐng)域全國重點實驗室招收博士1-2名
已經(jīng)有3人回復
-
材料專碩306英一數(shù)二
已經(jīng)有10人回復
-
301求調(diào)劑
已經(jīng)有6人回復
-
一志愿天津大學化學工藝專業(yè)（081702）315分求調(diào)劑
已經(jīng)有7人回復
-
302求調(diào)劑
已經(jīng)有6人回復
-
26博士申請
已經(jīng)有3人回復
» 本主題相關(guān)價值貼推薦，對您同樣有幫助:
- 請教各位前輩，我這樣的條件可以申請推免嗎？有沒有老師會要��？謝謝~~~
已經(jīng)有11人回復
- 【求助】有關(guān)標準曲線的疑問，請大家?guī)蛶兔�，謝謝！
已經(jīng)有4人回復
- 【請教】簇合物晶體
已經(jīng)有14人回復
- 【請教】問個空心金屬球的TEM表征問題
已經(jīng)有21人回復
zhangguangping
木蟲 (著名寫手)
- 1ST強帖: 23
- 應助: 71 (初中生)
- 貴賓: 0.031
- 金幣: 1973.1
- 散金: 9059
- 紅花: 79
- 帖子: 2681
- 在線: 2128.8小時
- 蟲號: 529624
- 注冊: 2008-03-20
- 性別: GG
- 專業(yè): 原子和分子物理
【答案】應助回帖
identation(金幣+20): 謝謝哈 2011-09-06 19:18:31
| 補充一下：我看到你的收斂輸出好像是一直走向收斂的。雖然輸出不正常了，但是數(shù)據(jù)在計算機內(nèi)存里面是正常的，所以結(jié)果也是可信的。只不過如果你關(guān)心的數(shù)據(jù)恰好被打上了星號，那就無奈了。這個說明你關(guān)心的這個數(shù)據(jù)隨著你的mesh cutoff沒有一個收斂的趨勢。另外你要想讓輸出的這些星號變成數(shù)值，簡單的辦法就是到程序里面找到輸出這些數(shù)的代碼，然后把格式輸出的寬度增加就行。 |
3樓2011-09-06 19:13:34
zhangguangping
木蟲 (著名寫手)
- 1ST強帖: 23
- 應助: 71 (初中生)
- 貴賓: 0.031
- 金幣: 1973.1
- 散金: 9059
- 紅花: 79
- 帖子: 2681
- 在線: 2128.8小時
- 蟲號: 529624
- 注冊: 2008-03-20
- 性別: GG
- 專業(yè): 原子和分子物理
2樓2011-09-06 19:10:07
4樓2011-09-06 19:20:23
zhangguangping
木蟲 (著名寫手)
- 1ST強帖: 23
- 應助: 71 (初中生)
- 貴賓: 0.031
- 金幣: 1973.1
- 散金: 9059
- 紅花: 79
- 帖子: 2681
- 在線: 2128.8小時
- 蟲號: 529624
- 注冊: 2008-03-20
- 性別: GG
- 專業(yè): 原子和分子物理
5樓2011-09-06 20:13:22
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 303求調(diào)劑 +3 | 睿08 2026-03-17 | 5/250 |
|
|---|---|---|---|---|
|
|
[考研] 277調(diào)劑 +5 | 自由煎餅果子 2026-03-16 | 6/300 |
|
|
|
[考研] 化學工程321分求調(diào)劑 +11 | 大米飯！ 2026-03-15 | 14/700 |
|
|
[考研] 085601求調(diào)劑 +4 | Du.11 2026-03-16 | 4/200 |
|
|
|
[考研] 271求調(diào)劑 +12 | 生如夏花… 2026-03-11 | 14/700 |
|
|
|
[考研] 材料與化工304求B區(qū)調(diào)劑 +7 | 邱gl 2026-03-11 | 8/400 |
|
|
|
[考研] 考研調(diào)劑 +3 | 淇ya_~ 2026-03-17 | 5/250 |
|
|
|
[考研] 環(huán)境工程調(diào)劑 +6 | 大可digkids 2026-03-16 | 6/300 |
|
|
|
[考研] 318求調(diào)劑 +3 | Yanyali 2026-03-15 | 3/150 |
|
|
|
[考研] 0856求調(diào)劑 +3 | 劉夢微 2026-03-15 | 3/150 |
|
|
|
[考研] 求老師收留調(diào)劑 +4 | jiang姜66 2026-03-14 | 5/250 |
|
|
|
[考研] 255求調(diào)劑 +3 | 李嘉慧， 2026-03-12 | 4/200 |
|
|
|
[考研] 266求調(diào)劑 +4 | 學員97LZgn 2026-03-13 | 4/200 |
|
|
|
[考研] 336求調(diào)劑 +6 | Iuruoh 2026-03-11 | 6/300 |
|
|
|
[考研] 四川大學085601材料工程專碩 初試294求調(diào)劑 +4 | 祝我們好在冬天 2026-03-11 | 4/200 |
|
|
|
[考研] 一志愿211化學學碩310分求調(diào)劑 +8 | 努力奮斗112 2026-03-12 | 9/450 |
|
|
|
[考研] 308求調(diào)劑 +3 | 是Lupa啊 2026-03-12 | 3/150 |
|
|
|
[考研] 0856化學工程280分求調(diào)劑 +4 | shenzxsn 2026-03-11 | 4/200 |
|
|
|
[考研] 283求調(diào)劑，材料、化工皆可 +8 | 蘇打水7777 2026-03-11 | 10/500 |
|
|
|
[考博] 2026年博士申請 +3 | QwQwQW10 2026-03-11 | 3/150 |
|