版塊導航: 正在加載中...

登錄注冊

應《網(wǎng)絡安全法》要求，自2017年10月1日起，未進行實名認證將不得使用互聯(lián)網(wǎng)跟帖服務。為保障您的帳號能夠正常使用，請盡快對帳號進行手機號驗證，感謝您的理解與支持！

24小時熱門版塊排行榜

返回列表

identation

金蟲 (著名寫手)

1ST強帖: 16
應助: 46 (小學生)
貴賓: 2.516
金幣: 253.4
散金: 2603
紅花: 44
帖子: 1283
在線: 500.5小時
蟲號: 842155
注冊: 2009-09-06
性別: GG
專業(yè): 凝聚態(tài)物性 II ：電子結構

[求助] 請教，為什么會出現(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
之前的計算中并無問題，請問怎么回事呢？謝謝

回復此樓

» 猜你喜歡

293求調劑已經有14人回復
0703化學336分求調劑已經有6人回復
304求調劑已經有9人回復
281求調劑（0805）已經有8人回復
環(huán)境領域全國重點實驗室招收博士1-2名已經有3人回復
材料專碩306英一數(shù)二已經有10人回復
301求調劑已經有6人回復
一志愿天津大學化學工藝專業(yè)（081702）315分求調劑已經有7人回復
302求調劑已經有6人回復
26博士申請已經有3人回復

» 本主題相關價值貼推薦，對您同樣有幫助:

請教各位前輩，我這樣的條件可以申請推免嗎？有沒有老師會要�。恐x謝~~~ 已經有11人回復
【求助】有關標準曲線的疑問，請大家?guī)蛶兔Γx謝！已經有4人回復
【請教】簇合物晶體已經有14人回復
【請教】問個空心金屬球的TEM表征問題已經有21人回復

1樓 2011-09-06 17:47:58

已閱回復此樓關注TA 給TA發(fā)消息送TA紅花 TA的回帖

zhangguangping

木蟲 (著名寫手)

1ST強帖: 23
應助: 71 (初中生)
貴賓: 0.031
金幣: 1973.1
散金: 9059
紅花: 79
帖子: 2681
在線: 2128.8小時
蟲號: 529624
注冊: 2008-03-20
性別: GG
專業(yè): 原子和分子物理

【答案】應助回帖

identation(金幣+30): 2011-09-10 21:13:38

這個是因為這些打著星號的地方的數(shù)值超出了輸出字符格式輸出的寬度。你看看能量分解的那地方，打著星號的地方的數(shù)，然后看看這個數(shù)是不是隨著你的mesh cutoff的增大而增大。

贊一下

回復此樓

弘德明志博學篤行

2樓2011-09-06 19:10:07

已閱回復此樓關注TA 給TA發(fā)消息送TA紅花 TA的回帖

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ù)在計算機內存里面是正常的，所以結果也是可信的。只不過如果你關心的數(shù)據(jù)恰好被打上了星號，那就無奈了。這個說明你關心的這個數(shù)據(jù)隨著你的mesh cutoff沒有一個收斂的趨勢。另外你要想讓輸出的這些星號變成數(shù)值，簡單的辦法就是到程序里面找到輸出這些數(shù)的代碼，然后把格式輸出的寬度增加就行。

贊一下

回復此樓

弘德明志博學篤行

3樓2011-09-06 19:13:34

已閱回復此樓關注TA 給TA發(fā)消息送TA紅花 TA的回帖

identation

金蟲 (著名寫手)

1ST強帖: 16
應助: 46 (小學生)
貴賓: 2.516
金幣: 253.4
散金: 2603
紅花: 44
帖子: 1283
在線: 500.5小時
蟲號: 842155
注冊: 2009-09-06
性別: GG
專業(yè): 凝聚態(tài)物性 II ：電子結構

引用回帖:

3樓: Originally posted by zhangguangping at 2011-09-06 19:13:34:
補充一下：我看到你的收斂輸出好像是一直走向收斂的。雖然輸出不正常了，但是數(shù)據(jù)在計算機內存里面是正常的，所以結果也是可信的。只不過如果你關心的數(shù)據(jù)恰好被打上了星號，那就無奈了。這個說明你關心的這個數(shù)據(jù) ...

測試meshcutoff時，收斂標準是怎么樣的呢？Mesh cutoff大體上用多大數(shù)值呢？我的meshcutoff從600變到650RY，能量差大約是0.02eV，force和stress一般小數(shù)點后第三位開始變化那取mesh cutoff為600ry，結果可信么？謝謝您

贊一下

回復此樓

4樓2011-09-06 19:20:23

已閱回復此樓關注TA 給TA發(fā)消息送TA紅花 TA的回帖

zhangguangping

木蟲 (著名寫手)

1ST強帖: 23
應助: 71 (初中生)
貴賓: 0.031
金幣: 1973.1
散金: 9059
紅花: 79
帖子: 2681
在線: 2128.8小時
蟲號: 529624
注冊: 2008-03-20
性別: GG
專業(yè): 原子和分子物理

【答案】應助回帖

引用回帖:

4樓: Originally posted by identation at 2011-09-06 19:20:23:
測試meshcutoff時，收斂標準是怎么樣的呢？Mesh cutoff大體上用多大數(shù)值呢？我的meshcutoff從600變到650RY，能量差大約是0.02eV，force和stress一般小數(shù)點后第三位開始變化那取mesh cutoff為600ry，結果可信么？ ...

收斂標準不是你自己定的嗎？
我看到好像你定義的是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的話，當取750Ry的時候，為什么能量出現(xiàn)了很大的數(shù)，導致超出了輸出長度呢？

贊一下

回復此樓

弘德明志博學篤行

5樓2011-09-06 20:13:22

已閱回復此樓關注TA 給TA發(fā)消息送TA紅花 TA的回帖

identation

金蟲 (著名寫手)

1ST強帖: 16
應助: 46 (小學生)
貴賓: 2.516
金幣: 253.4
散金: 2603
紅花: 44
帖子: 1283
在線: 500.5小時
蟲號: 842155
注冊: 2009-09-06
性別: GG
專業(yè): 凝聚態(tài)物性 II ：電子結構

引用回帖:

5樓: Originally posted by zhangguangping at 2011-09-06 20:13:22:
收斂標準不是你自己定的嗎？
我看到好像你定義的是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
...

我指的收斂是相對于mesh cutoff來講能量的收斂啊
另外
MeshCutoff 750. Ry
FilterCutoff 750. Ry

這里FilterCutoff 怎么指定呢？如果僅制定mesh cutoff，其余的用默認值可以么？
謝謝

贊一下

回復此樓

6樓2011-09-06 20:19:36

已閱回復此樓關注TA 給TA發(fā)消息送TA紅花 TA的回帖

zhangguangping

木蟲 (著名寫手)

1ST強帖: 23
應助: 71 (初中生)
貴賓: 0.031
金幣: 1973.1
散金: 9059
紅花: 79
帖子: 2681
在線: 2128.8小時
蟲號: 529624
注冊: 2008-03-20
性別: GG
專業(yè): 原子和分子物理

【答案】應助回帖

★
uuv2010(金幣+1): 多謝提示 2011-09-06 22:46:07

引用回帖:

6樓: Originally posted by identation at 2011-09-06 20:19:36:
我指的收斂是相對于mesh cutoff來講能量的收斂啊
另外
MeshCutoff 750. Ry
FilterCutoff 750. Ry

這里FilterCutoff 怎么指定呢？如果僅制定mesh cutoff，其余的用默認值可以么？
謝謝

這個我不是很清楚。但是我看到手冊上說要是用了FilterCutoff，并且跟著改變的話，就不要嘗試著做mesh cutoff的收斂測試。你正是做的這個測試。
所以看看手冊。把這個Filter的去掉或者固定一個值。

贊一下(1人)

回復此樓

弘德明志博學篤行

7樓2011-09-06 22:02:04

已閱回復此樓關注TA 給TA發(fā)消息送TA紅花 TA的回帖

相關版塊跳轉我要訂閱樓主 identation 的主題更新

返回列表

最具人氣熱帖推薦 [查看全部]		作者	回/看	最后發(fā)表

[考研] 268求調劑 +3	簡單點0 2026-03-17	3/150	2026-03-17 23:15 by wangjihu
[考研] 085600材料與化工 +4	安全上岸！ 2026-03-16	4/200	2026-03-17 14:02 by 勇敢太監(jiān)王公公
[考研] 0854可跨調劑，一作一項核心論文五項專利，省、國級證書40+數(shù)一英一287 +3	小李0854 2026-03-16	3/150	2026-03-17 13:40 by 熱情沙漠
[考研] 一志愿，福州大學材料專碩339分求調劑 +3	木子momo青爭 2026-03-15	3/150	2026-03-17 07:52 by laoshidan
[基金申請] 國自科面上基金字體 +6	iwuli 2026-03-12	7/350	2026-03-16 21:18 by sculhf
[考研] 一志愿211 0703方向310分求調劑 +3	努力奮斗112 2026-03-15	3/150	2026-03-16 16:44 by houyaoxu
[考研] 0703一志愿211 285分求調劑 +5	ly3471z 2026-03-13	5/250	2026-03-16 16:16 by 哦哦123
[考研] 070303一志愿西北大學學碩310找調劑 +5	d如愿上岸 2026-03-12	8/400	2026-03-16 15:19 by peike
[考研] 277材料科學與工程080500求調劑 +3	自由煎餅果子 2026-03-16	3/150	2026-03-16 14:10 by 運氣yunqi
[考研] 中科大材料與化工319求調劑 +3	孟鑫材料 2026-03-14	3/150	2026-03-14 20:10 by ms629
[考研] 材料與化工 323 英一+數(shù)二+物化，一志愿：哈工大本人本科雙一流 +4	自由的_飛翔 2026-03-13	5/250	2026-03-14 19:39 by hmn_wj
[考研] 265求調劑 +4	威化餅07 2026-03-12	4/200	2026-03-14 17:23 by userper
[考研] 招收0805（材料）調劑 +3	18595523086 2026-03-13	3/150	2026-03-14 00:33 by 123%、
[考研] 304求調劑 +7	7712b 2026-03-13	7/350	2026-03-13 21:42 by peike
[考研] 求調劑 +5	一定有學上- 2026-03-12	5/250	2026-03-13 18:31 by ms629
[考研] 0703化學求調劑 +7	綠豆芹菜湯 2026-03-12	7/350	2026-03-13 17:25 by njzyff
[考研] 工科278分求調劑 +5	周慢熱啊 2026-03-12	7/350	2026-03-13 15:49 by JourneyLucky
[考研] 0817化學工程與技術考研312分調劑 +3	T123 tt 2026-03-12	3/150	2026-03-13 10:49 by houyaoxu
[考研] 290求調劑 +3	ADT 2026-03-13	3/150	2026-03-13 10:19 by peike
[考研] 化工學碩306求調劑 +9	42838695 2026-03-12	9/450	2026-03-13 10:16 by houyaoxu

24小時熱門版塊排行榜

[求助] 請教，為什么會出現(xiàn)這樣的情況？謝謝大家

» 猜你喜歡

» 本主題相關價值貼推薦，對您同樣有幫助:

【答案】應助回帖

【答案】應助回帖

【答案】應助回帖

【答案】應助回帖

[求助] 請教，為什么會出現(xiàn)這樣的情況？謝謝大家

» 本主題相關價值貼推薦，對您同樣有幫助: