 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 1199 | 回復(fù): 6 | ||
| 當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖 | ||
identation金蟲 (著名寫手)
|
[求助]
請(qǐng)教，為什么會(huì)出現(xiàn)這樣的情況？謝謝大家
|
|
|
我在測(cè)試mesh cutoff的時(shí)候，一直很順利，mesh cutoff從200Ry一直增加到750ry，200到700ry的測(cè)試中都很順利完成了scf運(yùn)算，但是當(dāng)mesh cutoff增加到750ry時(shí)，出現(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******** 請(qǐng)問是為什么呢？ 我的輸入文件如下： 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 之前的計(jì)算中并無問題，請(qǐng)問怎么回事呢？謝謝 |
回復(fù)此樓» 猜你喜歡
0703化學(xué)336分求調(diào)劑
已經(jīng)有6人回復(fù)
-
304求調(diào)劑
已經(jīng)有9人回復(fù)
-
293求調(diào)劑
已經(jīng)有13人回復(fù)
-
281求調(diào)劑（0805）
已經(jīng)有8人回復(fù)
-
環(huán)境領(lǐng)域全國(guó)重點(diǎn)實(shí)驗(yàn)室招收博士1-2名
已經(jīng)有3人回復(fù)
-
材料專碩306英一數(shù)二
已經(jīng)有10人回復(fù)
-
301求調(diào)劑
已經(jīng)有6人回復(fù)
-
一志愿天津大學(xué)化學(xué)工藝專業(yè)（081702）315分求調(diào)劑
已經(jīng)有7人回復(fù)
-
302求調(diào)劑
已經(jīng)有6人回復(fù)
-
26博士申請(qǐng)
已經(jīng)有3人回復(fù)
» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:
- 請(qǐng)教各位前輩，我這樣的條件可以申請(qǐng)推免嗎？有沒有老師會(huì)要啊？謝謝~~~
已經(jīng)有11人回復(fù)
- 【求助】有關(guān)標(biāo)準(zhǔn)曲線的疑問，請(qǐng)大家?guī)蛶兔�，謝謝！
已經(jīng)有4人回復(fù)
- 【請(qǐng)教】簇合物晶體
已經(jīng)有14人回復(fù)
- 【請(qǐng)教】問個(gè)空心金屬球的TEM表征問題
已經(jīng)有21人回復(fù)
zhangguangping
木蟲 (著名寫手)
- 1ST強(qiáng)帖: 23
- 應(yīng)助: 71 (初中生)
- 貴賓: 0.031
- 金幣: 1973.1
- 散金: 9059
- 紅花: 79
- 帖子: 2681
- 在線: 2128.8小時(shí)
- 蟲號(hào): 529624
- 注冊(cè): 2008-03-20
- 性別: GG
- 專業(yè): 原子和分子物理
7樓2011-09-06 22:02:04
zhangguangping
木蟲 (著名寫手)
- 1ST強(qiáng)帖: 23
- 應(yīng)助: 71 (初中生)
- 貴賓: 0.031
- 金幣: 1973.1
- 散金: 9059
- 紅花: 79
- 帖子: 2681
- 在線: 2128.8小時(shí)
- 蟲號(hào): 529624
- 注冊(cè): 2008-03-20
- 性別: GG
- 專業(yè): 原子和分子物理
2樓2011-09-06 19:10:07
zhangguangping
木蟲 (著名寫手)
- 1ST強(qiáng)帖: 23
- 應(yīng)助: 71 (初中生)
- 貴賓: 0.031
- 金幣: 1973.1
- 散金: 9059
- 紅花: 79
- 帖子: 2681
- 在線: 2128.8小時(shí)
- 蟲號(hào): 529624
- 注冊(cè): 2008-03-20
- 性別: GG
- 專業(yè): 原子和分子物理
【答案】應(yīng)助回帖
identation(金幣+20): 謝謝哈 2011-09-06 19:18:31
| 補(bǔ)充一下：我看到你的收斂輸出好像是一直走向收斂的。雖然輸出不正常了，但是數(shù)據(jù)在計(jì)算機(jī)內(nèi)存里面是正常的，所以結(jié)果也是可信的。只不過如果你關(guān)心的數(shù)據(jù)恰好被打上了星號(hào)，那就無奈了。這個(gè)說明你關(guān)心的這個(gè)數(shù)據(jù)隨著你的mesh cutoff沒有一個(gè)收斂的趨勢(shì)。另外你要想讓輸出的這些星號(hào)變成數(shù)值，簡(jiǎn)單的辦法就是到程序里面找到輸出這些數(shù)的代碼，然后把格式輸出的寬度增加就行。 |
3樓2011-09-06 19:13:34
identation
金蟲 (著名寫手)
- 1ST強(qiáng)帖: 16
- 應(yīng)助: 46 (小學(xué)生)
- 貴賓: 2.516
- 金幣: 253.4
- 散金: 2603
- 紅花: 44
- 帖子: 1283
- 在線: 500.5小時(shí)
- 蟲號(hào): 842155
- 注冊(cè): 2009-09-06
- 性別: GG
- 專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)
4樓2011-09-06 19:20:23
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 085601材料工程專碩求調(diào)劑 +5 | 慕寒mio 2026-03-16 | 5/250 |
|
|---|---|---|---|---|
|
|
[考研] 材料與化工求調(diào)劑 +6 | 為學(xué)666 2026-03-16 | 6/300 |
|
|
[考研] 070300化學(xué)319求調(diào)劑 +4 | 錦鯉0909 2026-03-17 | 4/200 |
|
|
|
[考研] 化學(xué)工程321分求調(diào)劑 +11 | 大米飯！ 2026-03-15 | 14/700 |
|
|
|
[考研] 332求調(diào)劑 +6 | Zz版 2026-03-13 | 6/300 |
|
|
|
[考研] 085600材料與化工 +4 | 安全上岸！ 2026-03-16 | 4/200 |
|
|
|
[考研] 材料專碩274一志愿陜西師范大學(xué)求調(diào)劑 +5 | 薛云鵬 2026-03-13 | 5/250 |
|
|
|
[考研] 材料與化工304求B區(qū)調(diào)劑 +7 | 邱gl 2026-03-11 | 8/400 |
|
|
|
[論文投稿] 有沒有大佬發(fā)小論文能帶我個(gè)二作 +3 | 增銳漏人 2026-03-17 | 4/200 |
|
|
|
[考研] 278求調(diào)劑 +3 | Yy7400 2026-03-13 | 3/150 |
|
|
|
[考研] [導(dǎo)師推薦]西南科技大學(xué)國(guó)防/材料導(dǎo)師推薦 +3 | 尖角小荷 2026-03-16 | 6/300 |
|
|
|
[考研] 318求調(diào)劑 +3 | Yanyali 2026-03-15 | 3/150 |
|
|
|
[考研] 297求調(diào)劑 +4 | 學(xué)海漂泊 2026-03-13 | 4/200 |
|
|
|
[考研] 學(xué)碩285求調(diào)劑 +13 | Wisjxn 2026-03-12 | 46/2300 |
|
|
|
[考研] 333求調(diào)劑 +3 | 球球古力 2026-03-11 | 3/150 |
|
|
|
[考研] 310求調(diào)劑 +3 | 【上上簽】 2026-03-11 | 3/150 |
|
|
|
[考研] 307求調(diào)劑 +5 | 超級(jí)伊昂大王 2026-03-12 | 5/250 |
|
|
|
[考研] 工科材料085601 279求調(diào)劑 +8 | 困于星晨 2026-03-12 | 10/500 |
|
|
|
[考研] 0817化學(xué)工程與技術(shù)考研312分調(diào)劑 +3 | T123 tt 2026-03-12 | 3/150 |
|
|
|
[考研] 290求調(diào)劑 +3 | ADT 2026-03-13 | 3/150 |
|