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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Siesta&Smeagol&Atk » 请教，为什么会出现这样的情况？谢谢大家
北京石油化工学院2026年研究生招生接收调剂公告
51/1返回列表
查看: 1233  |  回复: 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

木虫 (著名写手)

【答案】应助回帖

★
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发消息 送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 个回答
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 307求调剂 +3 所念及所望 2026-04-06 3/150 2026-04-06 17:30 by 土木硕士招生
[考研] 319分085702安全工程求调剂 +4 rious 2026-04-05 4/200 2026-04-06 11:08 by 静静静静静静静�
[考研] 化学0703-一志愿211-338分求调剂 +8 vants 2026-04-05 8/400 2026-04-06 06:17 by houyaoxu
[考研] 295求调剂 +8 FZAC123 2026-04-03 8/400 2026-04-05 17:46 by 蓝云思雨
[考研] 本科211 分数293请求调剂 +4 莲菜就是藕吧 2026-04-01 4/200 2026-04-04 22:32 by hemengdong
[考研] 材料与化工306分找调剂 +23 沧海轻舟e 2026-04-02 27/1350 2026-04-04 21:52 by laoshidan
[考研] 278求调剂 +3 依旧！ 2026-04-02 4/200 2026-04-04 20:27 by 蓝云思雨
[考研] 怎么删帖子啊 +3 缝曦1000 2026-04-04 3/150 2026-04-04 14:20 by 土木硕士招生
[考研] 0710生物学336分求调剂 +6 kiyy 2026-04-01 8/400 2026-04-04 10:10 by kiyy
[考研] 材料科学与工程考研 +10 拯救皮特托先生 2026-04-02 10/500 2026-04-03 23:57 by userper
[考研] 317分 一志愿江南大学 化学工程学硕 求调剂 +6 YinTai 2026-04-03 6/300 2026-04-03 22:30 by 无际的草原
[考研] 338求调剂 +7 晟功? 2026-04-03 7/350 2026-04-03 16:46 by wxiongid
[考研] 289-求调剂 +4 这里是_ 2026-04-03 4/200 2026-04-03 14:23 by 1753564080
[考研] 求调剂 +3 usbdndj 2026-04-03 3/150 2026-04-03 14:10 by dxiaoxin
[考研] 求调剂 +9 akdhjs 2026-03-31 11/550 2026-04-03 13:32 by akdhjs
[考研] 302求调剂 +9 zyx上岸！ 2026-04-02 9/450 2026-04-02 23:07 by 马儿快快地跑
[考研] 土木304求调剂 +4 兔突突突， 2026-04-02 5/250 2026-04-02 21:16 by 兔突突突，
[考研] 282求调剂 +13 呼吸都是减肥 2026-04-01 13/650 2026-04-02 14:10 by baoball
[考研] 266分，一志愿电气工程，本科材料，求材料专业调剂 +10 哇呼哼呼哼 2026-04-01 11/550 2026-04-02 11:31 by lnilvy
[考研] 270调剂 +7 maxjxbsk 2026-04-02 7/350 2026-04-02 09:50 by yulian1987
信息提示
关闭
请填处理意见
关闭