应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Siesta&Smeagol&Atk » 请教，为什么会出现这样的情况？谢谢大家
北京石油化工学院2026年研究生招生接收调剂公告
51/1返回列表
查看: 1234  |  回复: 6
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

identation

金虫 (著名写手)

[求助] 请教，为什么会出现这样的情况？谢谢大家

我在测试mesh cutoff的时候，一直很顺利，mesh cutoff从200Ry一直增加到750ry，200到700ry的测试中都很顺利完成了scf运算，但是当mesh cutoff增加到750ry时，出现了这样的输出情况

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
之前的计算中并无问题，请问怎么回事呢？谢谢
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2011-09-06 17:47:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhangguangping

木虫 (著名写手)

【答案】应助回帖

引用回帖:
4楼: Originally posted by identation at 2011-09-06 19:20:23:
测试meshcutoff时，收敛标准是怎么样的呢？Mesh cutoff大体上用多大数值呢？我的meshcutoff从600变到650RY，能量差大约是0.02eV，force和stress一般小数点后第三位开始变化 那取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的时候，为什么能量出现了很大的数，导致超出了输出长度呢？
赞一下 回复此楼
弘德明志博学笃行
5楼2011-09-06 20:13:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 7 个回答

zhangguangping

木虫 (著名写手)

【答案】应助回帖

identation(金币+30): 2011-09-10 21:13:38
这个是因为这些打着星号的地方的数值超出了输出字符格式输出的宽度。你看看能量分解的那地方，打着星号的地方的数，然后看看这个数是不是随着你的mesh cutoff的增大而增大。
赞一下 回复此楼
弘德明志博学笃行
2楼2011-09-06 19:10:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhangguangping

木虫 (著名写手)

【答案】应助回帖

identation(金币+20): 谢谢哈 2011-09-06 19:18:31
补充一下：我看到你的收敛输出好像是一直走向收敛的。虽然输出不正常了，但是数据在计算机内存里面是正常的，所以结果也是可信的。只不过如果你关心的数据恰好被打上了星号，那就无奈了。这个说明你关心的这个数据随着你的mesh cutoff没有一个收敛的趋势。另外你要想让输出的这些星号变成数值，简单的办法就是到程序里面找到输出这些数的代码，然后把格式输出的宽度增加就行。
赞一下 回复此楼
弘德明志博学笃行
3楼2011-09-06 19:13:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

identation

金虫 (著名写手)
引用回帖:
3楼: Originally posted by zhangguangping at 2011-09-06 19:13:34:
补充一下：我看到你的收敛输出好像是一直走向收敛的。虽然输出不正常了，但是数据在计算机内存里面是正常的，所以结果也是可信的。只不过如果你关心的数据恰好被打上了星号，那就无奈了。这个说明你关心的这个数据 ...

测试meshcutoff时，收敛标准是怎么样的呢？Mesh cutoff大体上用多大数值呢？我的meshcutoff从600变到650RY，能量差大约是0.02eV，force和stress一般小数点后第三位开始变化 那取mesh cutoff为600ry，结果可信么？谢谢您
赞一下 回复此楼
4楼2011-09-06 19:20:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 7 个回答
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 0854求调剂 +4 assdll 2026-04-05 4/200 2026-04-06 12:29 by 中飞院空管学院�
[考研] 一志愿南航，数一英一学硕317求调剂！！ +6 Acaciad 2026-04-04 6/300 2026-04-06 12:13 by 考研学校招点人
[考研] 0854电子信息319求调剂（接受跨专业调剂） +3 星星不眨眼喽 2026-04-05 3/150 2026-04-05 20:20 by 啵啵啵0119
[考研] 材料专硕(0856) 339分求调剂 +10 哈哈哈鹅哈哈哈 2026-04-04 10/500 2026-04-05 18:51 by 蓝云思雨
[考研] 08专硕275调剂 +5 AaAa7420 2026-04-05 5/250 2026-04-05 18:01 by jkddd
[考研] 0832食品科学与工程学硕282调剂 +6 鱼在水中游a 2026-04-02 9/450 2026-04-05 11:45 by flysky1234
[考研] 312求调剂 +3 Say Never 2026-04-04 3/150 2026-04-05 11:19 by guoweigw
[考研] 085602调剂 初试总分335 +12 19123253302 2026-04-04 12/600 2026-04-05 08:08 by 544594351
[考研] 求生物学学硕调剂——364分 +7 云朵遛弯指南 2026-04-04 7/350 2026-04-04 22:49 by zhyzzh
[考研] 11408,335分，本科211，求调剂，可转专业 +5 鳄梨大鳄鱼 2026-04-03 5/250 2026-04-04 22:49 by chongya
[考研] 283求调剂 +4 mcbbc 2026-04-03 5/250 2026-04-04 20:51 by imissbao
[考研] 400分求调剂 +3 尴尬且挠头 2026-04-04 3/150 2026-04-04 08:41 by jp9609
[考研] 调剂0855-288 +5 x熊二a 2026-04-03 5/250 2026-04-04 00:19 by 猪会飞
[考研] 材料科学与工程考研 +10 拯救皮特托先生 2026-04-02 10/500 2026-04-03 23:57 by userper
[考研] 266求调剂 +3 08电气工程 2026-04-03 3/150 2026-04-03 14:05 by 1753564080
[考研] 315分 085602 求调剂 +15 26考研上岸版26 2026-04-02 15/750 2026-04-03 12:45 by xingguangj
[考研] 求调剂 +3 心想事成可 2026-04-03 3/150 2026-04-03 11:22 by wangjy2002
[考研] 274求调剂 +10 薛定谔的虎。 2026-04-01 10/500 2026-04-03 10:13 by tianyyysss
[考研] 261求B区调剂 +5 明仔· 2026-04-01 7/350 2026-04-02 11:17 by 邹尉尉
[考研] 085600，321分求调剂 +13 大馋小子 2026-03-31 13/650 2026-04-01 12:35 by chemdavid
信息提示
关闭
请填处理意见
关闭