 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 1205 | 回復(fù): 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 之前的計算中并無問題，請問怎么回事呢？謝謝 |
回復(fù)此樓» 猜你喜歡
一志愿中國海洋大學(xué)，生物學(xué)，301分，求調(diào)劑
已經(jīng)有4人回復(fù)
-
一志愿武理材料305分求調(diào)劑
已經(jīng)有6人回復(fù)
-
材料專碩英一數(shù)二306
已經(jīng)有3人回復(fù)
-
085600材料與化工調(diào)劑 324分
已經(jīng)有7人回復(fù)
-
311求調(diào)劑
已經(jīng)有3人回復(fù)
-
26調(diào)劑/材料/英一數(shù)二/總分289/已過A區(qū)線
已經(jīng)有7人回復(fù)
-
295求調(diào)劑
已經(jīng)有5人回復(fù)
-
354求調(diào)劑
已經(jīng)有5人回復(fù)
-
材料專業(yè)求調(diào)劑
已經(jīng)有3人回復(fù)
-
求材料調(diào)劑
已經(jīng)有7人回復(fù)
» 本主題相關(guān)價值貼推薦，對您同樣有幫助:
- 請教各位前輩，我這樣的條件可以申請推免嗎？有沒有老師會要�。恐x謝~~~
已經(jīng)有11人回復(fù)
- 【求助】有關(guān)標準曲線的疑問，請大家?guī)蛶兔�，謝謝！
已經(jīng)有4人回復(fù)
- 【請教】簇合物晶體
已經(jīng)有14人回復(fù)
- 【請教】問個空心金屬球的TEM表征問題
已經(jīng)有21人回復(fù)
6樓2011-09-06 20:19:36
zhangguangping
木蟲 (著名寫手)
- 1ST強帖: 23
- 應(yīng)助: 71 (初中生)
- 貴賓: 0.031
- 金幣: 1973.1
- 散金: 9059
- 紅花: 79
- 帖子: 2681
- 在線: 2128.8小時
- 蟲號: 529624
- 注冊: 2008-03-20
- 性別: GG
- 專業(yè): 原子和分子物理
2樓2011-09-06 19:10:07
zhangguangping
木蟲 (著名寫手)
- 1ST強帖: 23
- 應(yīng)助: 71 (初中生)
- 貴賓: 0.031
- 金幣: 1973.1
- 散金: 9059
- 紅花: 79
- 帖子: 2681
- 在線: 2128.8小時
- 蟲號: 529624
- 注冊: 2008-03-20
- 性別: GG
- 專業(yè): 原子和分子物理
【答案】應(yīng)助回帖
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
4樓2011-09-06 19:20:23
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿中國海洋大學(xué)，生物學(xué)，301分，求調(diào)劑 +4 | 1孫悟空 2026-03-17 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 一志愿武理材料305分求調(diào)劑 +5 | 想上岸的鯉魚 2026-03-18 | 6/300 |
|
|
|
[考研] 一志愿天大材料與化工（085600）總分338 +3 | 蔡大美女 2026-03-13 | 3/150 |
|
|
|
[考研] 344求調(diào)劑 +5 | knight344 2026-03-16 | 6/300 |
|
|
|
[考研] 0854可跨調(diào)劑，一作一項核心論文五項專利，省、國級證書40+數(shù)一英一287 +8 | 小李0854 2026-03-16 | 8/400 |
|
|
|
[考研] 331求調(diào)劑（0703有機化學(xué) +7 | ZY-05 2026-03-13 | 8/400 |
|
|
|
[考研] 0703化學(xué)336分求調(diào)劑 +6 | zbzihdhd 2026-03-15 | 7/350 |
|
|
|
[考研] 301求調(diào)劑 +9 | yy要上岸呀 2026-03-17 | 9/450 |
|
|
|
[考研] 296求調(diào)劑 +5 | 大口吃飯 身體健 2026-03-13 | 5/250 |
|
|
|
[考研] 326求調(diào)劑 +5 | 上岸的小葡 2026-03-15 | 6/300 |
|
|
|
[考研] 288求調(diào)劑 +4 | 奇點0314 2026-03-14 | 4/200 |
|
|
|
[考研] 304求調(diào)劑 +6 | Mochaaaa 2026-03-12 | 7/350 |
|
|
|
[考研] 材料與化工求調(diào)劑一志愿 985 總分 295 +8 | dream…… 2026-03-12 | 8/400 |
|
|
|
[考研] 求調(diào)劑 +7 | 18880831720 2026-03-11 | 7/350 |
|
|
|
[考研] 工科材料085601 279求調(diào)劑 +8 | 困于星晨 2026-03-12 | 10/500 |
|
|
|
[考研] 295求調(diào)劑 +3 | 小匕仔汁 2026-03-12 | 3/150 |
|
|
|
[考研] 328化工專碩求調(diào)劑 +4 | 。，。，。，。i 2026-03-12 | 4/200 |
|
|
|
[論文投稿]
投稿問題
 5+4
|
星光燦爛xt 2026-03-12 | 6/300 |
|
|
|
[考研] 290求調(diào)劑 +3 | ADT 2026-03-13 | 3/150 |
|
|
|
[考研] 081200-11408-276學(xué)碩求調(diào)劑 +3 | 崔wj 2026-03-12 | 4/200 |
|