24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

lijuanshu

铜虫 (正式写手)

应助: 1 (幼儿园)
金币: 330.2
帖子: 360
在线: 150.7小时
虫号: 1415189

[交流] 静态计算时出现警告

静态计算时的INCAR如下：
SYSTEM =
ENCUT =
ISTART = 0
ICHARG = 2
ISMEAR = -5
ISPIN = 2
NELM = 2000
EDIFF = 1E-4
IALGO = 48
LREAL = Auto
GGA = 91
VOSKOWN = 1
但是在计算中出现警告如下：
RMM:  22 -0.310967594429E+03 0.21293E+00 -0.24248E+00 10244 0.439E+00 0.422E+01
  WARNING: DENTET: can't reach specified precision
  Number of Electrons is NELECT = 476.996768664236
RMM:  23 -0.310376643978E+03 0.59095E+00 -0.83636E-01 10268 0.266E+00 0.375E+01

RMM:  43 -0.309736141363E+03 0.53076E-02 -0.50929E-02 10165 0.574E-01 0.356E+00
RMM:  44 -0.309811626871E+03 -0.75486E-01 -0.48571E-02 10210 0.597E-01 0.369E+00
  WARNING: DENTET: can't reach specified precision
  Number of Electrons is NELECT = 477.041023713897
RMM:  45 -0.309698348515E+03 0.11328E+00 -0.13902E-02 10260 0.310E-01BRMIX: very serious problems
the old and the new charge density differ
old charge density: 476.99930 new  477.04032
0.319E+00
  WARNING: DENTET: can't reach specified precision
  Number of Electrons is NELECT = 476.985853390457
RMM:  46 -0.310016610388E+03 -0.31826E+00 -0.42439E-02 10215 0.447E-01BRMIX: very serious problems
the old and the new charge density differ
old charge density: 477.04032 new  476.98515
0.296E+00
RMM:  47 -0.309997772396E+03 0.18838E-01 -0.27267E-02 10259 0.413E-01BRMIX: very serious problems
the old and the new charge density differ
old charge density: 476.98515 new  476.99930
0.198E+00
RMM:  48 -0.310064999937E+03 -0.67228E-01 -0.44514E-02 10200 0.397E-01 0.171E+00
这些警告对电荷密度的结果有影响吗？应该如何消除？

回复此楼

» 抢金币啦！回帖就可以得到:
查看全部散金贴

1楼2013-03-05 16:19:38

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

hypnos

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 11.8
帖子: 7
在线: 2小时
虫号: 1513520

★
小木虫: 金币+0.5, 给个红包，谢谢回帖

引用回帖:

11楼: Originally posted by lijuanshu at 2013-03-06 14:29:38
这三种方法都试了，只有第一种方法把ISMEAR改为0才可以消除这些警告？但是在静态计算的时候可以设置ISMEAR=0，SIGMA=0.2吗？期待你的回复，谢谢。...

静态计算时候我个人觉得可以用ISMEAR=0
但是不同系统对于ISMEAR的设置是不同的，=0的设置对大多数系统都是有效的。SIGMA取值需要测试，用多个SIGMA计算单点能，取一个合理的entropy。具体做法我不太清楚了，不妨查查论坛里类似的帖子。

比如你是小系统绝缘体或半导体，用-5和较多的KPOINTS
如果是大系统绝缘体，用=0并取小的SIGMA值，如0.05。
比如你是金属体系，用>0，调试合适的SIGMA值