 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 1200 | 回復(fù): 6 | ||
| 當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖 | ||
identation金蟲 (著名寫手)
|
[求助]
請教，為什么會出現(xiàn)這樣的情況？謝謝大家
|
|
|
我在測試mesh cutoff的時候，一直很順利，mesh cutoff從200Ry一直增加到750ry，200到700ry的測試中都很順利完成了scf運(yùn)算，但是當(dāng)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ù)此樓» 猜你喜歡
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)域全國重點(diǎ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博士申請
已經(jīng)有3人回復(fù)
» 本主題相關(guān)價值貼推薦，對您同樣有幫助:
- 請教各位前輩，我這樣的條件可以申請推免嗎？有沒有老師會要��？謝謝~~~
已經(jīng)有11人回復(fù)
- 【求助】有關(guān)標(biāo)準(zhǔn)曲線的疑問，請大家?guī)蛶兔�，謝謝！
已經(jīng)有4人回復(fù)
- 【請教】簇合物晶體
已經(jīng)有14人回復(fù)
- 【請教】問個空心金屬球的TEM表征問題
已經(jīng)有21人回復(fù)
identation
金蟲 (著名寫手)
- 1ST強(qiáng)帖: 16
- 應(yīng)助: 46 (小學(xué)生)
- 貴賓: 2.516
- 金幣: 253.4
- 散金: 2603
- 紅花: 44
- 帖子: 1283
- 在線: 500.5小時
- 蟲號: 842155
- 注冊: 2009-09-06
- 性別: GG
- 專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)
4樓2011-09-06 19:20:23
zhangguangping
木蟲 (著名寫手)
- 1ST強(qiáng)帖: 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強(qiáng)帖: 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
| 補(bǔ)充一下：我看到你的收斂輸出好像是一直走向收斂的。雖然輸出不正常了，但是數(shù)據(jù)在計算機(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強(qiáng)帖: 23
- 應(yīng)助: 71 (初中生)
- 貴賓: 0.031
- 金幣: 1973.1
- 散金: 9059
- 紅花: 79
- 帖子: 2681
- 在線: 2128.8小時
- 蟲號: 529624
- 注冊: 2008-03-20
- 性別: GG
- 專業(yè): 原子和分子物理
【答案】應(yīng)助回帖
|
收斂標(biāo)準(zhǔn)不是你自己定的嗎？ 我看到好像你定義的是DM收斂1.d-4 0.0553 0.0453 0.0027 0.0030 0.0009 0.0007 0.0003 0.0002 0.0002 0.0002 0.0002 0.0002 0.0001 0.0001 0.0002 0.0002 0.0002 0.0002 0.0002 0.0001 這個不就說明你的計算在走向收斂嗎？ 如果你的能量差別在0.002eV的話，當(dāng)取750Ry的時候，為什么能量出現(xiàn)了很大的數(shù)，導(dǎo)致超出了輸出長度呢？ |
5樓2011-09-06 20:13:22
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 302求調(diào)劑 +6 | 呼呼呼。。。。 2026-03-17 | 6/300 |
|
|---|---|---|---|---|
|
[考研] 能源材料化學(xué)課題組招收碩士研究生8-10名 +3 | 脫穎而出 2026-03-16 | 6/300 |
|
|
|
[考研] 材料與化工求調(diào)劑 +6 | 為學(xué)666 2026-03-16 | 6/300 |
|
|
|
[考研] 268求調(diào)劑 +8 | 一定有學(xué)上- 2026-03-14 | 9/450 |
|
|
|
[考研] 材料與化工專碩調(diào)劑 +5 | heming3743 2026-03-16 | 5/250 |
|
|
|
[考研] 302求調(diào)劑 +4 | 小賈同學(xué)123 2026-03-15 | 8/400 |
|
|
|
[考研] 267一志愿南京工業(yè)大學(xué)0817化工求調(diào)劑 +6 | SUICHILD 2026-03-12 | 6/300 |
|
|
|
[考研] 一志愿，福州大學(xué)材料專碩339分求調(diào)劑 +3 | 木子momo青爭 2026-03-15 | 3/150 |
|
|
|
[考研] [導(dǎo)師推薦]西南科技大學(xué)國防/材料導(dǎo)師推薦 +3 | 尖角小荷 2026-03-16 | 6/300 |
|
|
|
[考研] 333求調(diào)劑 +3 | 文思客 2026-03-16 | 7/350 |
|
|
|
[考研] 環(huán)境工程調(diào)劑 +6 | 大可digkids 2026-03-16 | 6/300 |
|
|
|
[考研] 304求調(diào)劑 +4 | ahbd 2026-03-14 | 4/200 |
|
|
|
[考研] 327求調(diào)劑 +6 | 拾光任染 2026-03-15 | 11/550 |
|
|
|
[考研] 080500，材料學(xué)碩302分求調(diào)劑學(xué)校 +4 | 初識可樂 2026-03-14 | 5/250 |
|
|
|
[考研] 265求調(diào)劑 +4 | 威化餅07 2026-03-12 | 4/200 |
|
|
|
[考研] 一志愿中科院，化學(xué)方向，295求調(diào)劑 +4 | 一氧二氮 2026-03-11 | 4/200 |
|
|
|
[考研] 材料與化工求調(diào)劑一志愿 985 總分 295 +8 | dream…… 2026-03-12 | 8/400 |
|
|
|
[考研] 307求調(diào)劑 +5 | 超級伊昂大王 2026-03-12 | 5/250 |
|
|
|
[考研] 328化工專碩求調(diào)劑 +4 | 。，。，。，。i 2026-03-12 | 4/200 |
|
|
|
[考博] 26讀博 +4 | Rui135246 2026-03-12 | 10/500 |
|