版塊導(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)版塊排行榜

返回列表

當(dāng)前只顯示滿足指定條件的回帖，點(diǎ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)有9人回復(fù)
293求調(diào)劑已經(jīng)有14人回復(fù)
0703化學(xué)336分求調(diào)劑已經(jīng)有6人回復(fù)
304求調(diào)劑已經(jīng)有9人回復(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)推免嗎？有沒(méi)有老師會(huì)要啊？謝謝~~~ 已經(jīng)有11人回復(fù)
【求助】有關(guān)標(biāo)準(zhǔn)曲線的疑問(wèn)，請(qǐng)大家?guī)蛶兔�，謝謝！已經(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)助回帖

引用回帖:

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的回帖

查看全部 7 個(gè)回答

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的回帖

查看全部 7 個(gè)回答

普通表情龍兔虎貓

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

[考研] 281求調(diào)劑（0805） +3	煙汐憶海 2026-03-16	9/450	2026-03-18 07:23 by 煙汐憶海
[考研] 268求調(diào)劑 +3	簡(jiǎn)單點(diǎn)0 2026-03-17	3/150	2026-03-17 23:15 by wangjihu
[考研] 材料專碩326求調(diào)劑 +6	墨煜姒莘 2026-03-15	7/350	2026-03-17 17:10 by ruiyingmiao
[考研] 085601求調(diào)劑 +4	Du.11 2026-03-16	4/200	2026-03-17 17:08 by ruiyingmiao
[考研] 材料專碩274一志愿陜西師范大學(xué)求調(diào)劑 +5	薛云鵬 2026-03-13	5/250	2026-03-17 10:15 by Sammy2
[考研] 278求調(diào)劑 +3	Yy7400 2026-03-13	3/150	2026-03-17 08:24 by laoshidan
[考研] 283求調(diào)劑 +3	聽(tīng)風(fēng)就是雨； 2026-03-16	3/150	2026-03-17 07:41 by 熱情沙漠
[基金申請(qǐng)] 國(guó)自科面上基金字體 +6	iwuli 2026-03-12	7/350	2026-03-16 21:18 by sculhf
[考研] 277材料科學(xué)與工程080500求調(diào)劑 +3	自由煎餅果子 2026-03-16	3/150	2026-03-16 14:10 by 運(yùn)氣yunqi
[考研] 機(jī)械專碩調(diào)劑 +3	笨笨兔子 2026-03-12	3/150	2026-03-15 20:02 by 栗子粥?
[考研] 070305求調(diào)劑 +3	mlpqaz03 2026-03-14	4/200	2026-03-15 11:04 by peike
[考研] 復(fù)試調(diào)劑 +3	呼呼？~+123456 2026-03-14	3/150	2026-03-14 16:53 by WTUChen
[考研] 材料工程調(diào)劑 +9	咪咪空空 2026-03-12	9/450	2026-03-13 22:05 by 星空星月
[考研] 290求調(diào)劑 +9	ADT 2026-03-11	9/450	2026-03-13 21:55 by JourneyLucky
[考研] 一志愿西南交大，材料專碩317求調(diào)劑 +5	lx8568 2026-03-11	5/250	2026-03-13 21:43 by peike
[考研] 【考研調(diào)劑求收留】 +3	Ceciilia 2026-03-11	3/150	2026-03-13 20:18 by JourneyLucky
[考研] 材料301分求調(diào)劑 +5	Liyouyumairs 2026-03-12	5/250	2026-03-13 14:42 by JourneyLucky
[考研] 070303一志愿西北大學(xué)學(xué)碩310找調(diào)劑 +3	d如愿上岸 2026-03-13	3/150	2026-03-13 10:43 by houyaoxu
[考研] 290求調(diào)劑 +3	ADT 2026-03-13	3/150	2026-03-13 10:19 by peike
[考研] 化工學(xué)碩306求調(diào)劑 +9	42838695 2026-03-12	9/450	2026-03-13 10:16 by houyaoxu