版塊導(dǎo)航: 正在加載中...

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

24小時(shí)熱門(mén)版塊排行榜

返回列表

identation

金蟲(chóng) (著名寫(xiě)手)

1ST強(qiáng)帖: 16
應(yīng)助: 46 (小學(xué)生)
貴賓: 2.516
金幣: 253.4
散金: 2603
紅花: 44
帖子: 1283
在線: 500.5小時(shí)
蟲(chóng)號(hào): 842155
注冊(cè): 2009-09-06
性別: GG
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

[求助] 請(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)問(wèn)是為什么呢？
我的輸入文件如下：

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ì)算中并無(wú)問(wèn)題，請(qǐng)問(wèn)怎么回事呢？謝謝

回復(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ù)
268求調(diào)劑已經(jīng)有3人回復(fù)
311求調(diào)劑已經(jīng)有10人回復(fù)
被我言中：新模板不強(qiáng)調(diào)格式了，假專家開(kāi)始管格式了已經(jīng)有4人回復(fù)

» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:

請(qǐng)教各位前輩，我這樣的條件可以申請(qǐng)推免嗎？有沒(méi)有老師會(huì)要�。恐x謝~~~ 已經(jīng)有11人回復(fù)
【求助】有關(guān)標(biāo)準(zhǔn)曲線的疑問(wèn)，請(qǐng)大家?guī)蛶兔Γx謝！已經(jīng)有4人回復(fù)
【請(qǐng)教】簇合物晶體已經(jīng)有14人回復(fù)
【請(qǐng)教】問(wèn)個(gè)空心金屬球的TEM表征問(wèn)題已經(jīng)有21人回復(fù)

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

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

zhangguangping

木蟲(chóng) (著名寫(xiě)手)

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

【答案】應(yīng)助回帖

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

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

贊一下

回復(fù)此樓

弘德明志博學(xué)篤行

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

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

zhangguangping

木蟲(chóng) (著名寫(xiě)手)

1ST強(qiáng)帖: 23
應(yīng)助: 71 (初中生)
貴賓: 0.031
金幣: 1973.1
散金: 9059
紅花: 79
帖子: 2681
在線: 2128.8小時(shí)
蟲(chóng)號(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ò)如果你關(guān)心的數(shù)據(jù)恰好被打上了星號(hào)，那就無(wú)奈了。這個(gè)說(shuō)明你關(guān)心的這個(gè)數(shù)據(jù)隨著你的mesh cutoff沒(méi)有一個(gè)收斂的趨勢(shì)。另外你要想讓輸出的這些星號(hào)變成數(shù)值，簡(jiǎn)單的辦法就是到程序里面找到輸出這些數(shù)的代碼，然后把格式輸出的寬度增加就行。

贊一下

回復(fù)此樓

弘德明志博學(xué)篤行

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

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

identation

金蟲(chóng) (著名寫(xiě)手)

1ST強(qiáng)帖: 16
應(yīng)助: 46 (小學(xué)生)
貴賓: 2.516
金幣: 253.4
散金: 2603
紅花: 44
帖子: 1283
在線: 500.5小時(shí)
蟲(chóng)號(hào): 842155
注冊(cè): 2009-09-06
性別: GG
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

引用回帖:

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

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

贊一下

回復(fù)此樓

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

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

zhangguangping

木蟲(chóng) (著名寫(xiě)手)

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

【答案】應(yīng)助回帖

引用回帖:

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

收斂標(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
這個(gè)不就說(shuō)明你的計(jì)算在走向收斂嗎？
如果你的能量差別在0.002eV的話，當(dāng)取750Ry的時(shí)候，為什么能量出現(xiàn)了很大的數(shù)，導(dǎo)致超出了輸出長(zhǎng)度呢？

贊一下

回復(fù)此樓

弘德明志博學(xué)篤行

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

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

identation

金蟲(chóng) (著名寫(xiě)手)

1ST強(qiáng)帖: 16
應(yīng)助: 46 (小學(xué)生)
貴賓: 2.516
金幣: 253.4
散金: 2603
紅花: 44
帖子: 1283
在線: 500.5小時(shí)
蟲(chóng)號(hào): 842155
注冊(cè): 2009-09-06
性別: GG
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

引用回帖:

5樓: Originally posted by zhangguangping at 2011-09-06 20:13:22:
收斂標(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
...

我指的收斂是相對(duì)于mesh cutoff來(lái)講能量的收斂啊
另外
MeshCutoff 750. Ry
FilterCutoff 750. Ry

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

贊一下

回復(fù)此樓

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

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

zhangguangping

木蟲(chóng) (著名寫(xiě)手)

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

【答案】應(yīng)助回帖

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

引用回帖:

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

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

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

贊一下(1人)

回復(fù)此樓

弘德明志博學(xué)篤行

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

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

相關(guān)版塊跳轉(zhuǎn) 我要訂閱樓主 identation 的主題更新

返回列表

普通表情龍兔虎貓

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

[考博] 26博士申請(qǐng) +3	1042136743 2026-03-17	3/150	2026-03-17 23:30 by 輕松不少隨
[考研] 311求調(diào)劑 +9	冬十三 2026-03-15	10/500	2026-03-17 22:53 by lbsjt
[考研] 材料與化工求調(diào)劑 +6	為學(xué)666 2026-03-16	6/300	2026-03-17 20:15 by peike
[考研] 268求調(diào)劑 +8	一定有學(xué)上- 2026-03-14	9/450	2026-03-17 17:47 by laoshidan
[考研] 293求調(diào)劑 +6	世界首富 2026-03-11	6/300	2026-03-17 17:04 by ruiyingmiao
[考研] 332求調(diào)劑 +6	Zz版 2026-03-13	6/300	2026-03-17 17:03 by ruiyingmiao
[考研] 08工科 320總分求調(diào)劑 +4	梨花珞晚風(fēng) 2026-03-17	4/200	2026-03-17 13:38 by houyaoxu
[考研] 275求調(diào)劑 +4	太陽(yáng)花天天開(kāi)心 2026-03-16	4/200	2026-03-17 10:53 by 功夫瘋狂
[考研] 材料與化工304求B區(qū)調(diào)劑 +7	邱gl 2026-03-11	8/400	2026-03-17 09:36 by 努力學(xué)習(xí)賺彩禮
[考研] 267一志愿南京工業(yè)大學(xué)0817化工求調(diào)劑 +6	SUICHILD 2026-03-12	6/300	2026-03-17 09:24 by 霧散后相遇lc
[考研] [導(dǎo)師推薦]西南科技大學(xué)國(guó)防/材料導(dǎo)師推薦 +3	尖角小荷 2026-03-16	6/300	2026-03-16 23:21 by 尖角小荷
[考研] 070305求調(diào)劑 +3	mlpqaz03 2026-03-14	4/200	2026-03-15 11:04 by peike
[考研] 308 085701 四六級(jí)已過(guò)求調(diào)劑 +7	溫喬喬喬喬 2026-03-12	14/700	2026-03-14 10:49 by JourneyLucky
[考研] 材料與化工（0856）304求B區(qū)調(diào)劑 +6	邱gl 2026-03-12	7/350	2026-03-13 23:24 by 邱gl
[考研] 337一志愿華南理工0805材料求調(diào)劑 +7	mysdl 2026-03-11	9/450	2026-03-13 22:43 by JourneyLucky
[考研] ［0860］321分求調(diào)劑，ab區(qū)皆可 +4	寶貴熱 2026-03-13	4/200	2026-03-13 22:01 by 星空星月
[考研] 一志愿西南交大，材料專碩317求調(diào)劑 +5	lx8568 2026-03-11	5/250	2026-03-13 21:43 by peike
[考研] 333求調(diào)劑 +3	球球古力 2026-03-11	3/150	2026-03-13 21:27 by JourneyLucky
[考研] （081700）化學(xué)工程與技術(shù)-298分求調(diào)劑 +12	11啦啦啦 2026-03-11	35/1750	2026-03-13 21:25 by JourneyLucky
[考研] 270求調(diào)劑 085600材料與化工專碩 +3	YXCT 2026-03-11	3/150	2026-03-13 10:13 by houyaoxu

24小時(shí)熱門(mén)版塊排行榜

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

» 猜你喜歡

» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:

【答案】應(yīng)助回帖

【答案】應(yīng)助回帖

【答案】應(yīng)助回帖

【答案】應(yīng)助回帖

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

» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助: