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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 静态计算时出现警告
201/1返回列表
查看: 1758  |  回复: 19
只看楼主@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

lijuanshu

铜虫 (正式写手)

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

静态计算时的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
这些警告对电荷密度的结果有影响吗？应该如何消除？
回复此楼

» 收录本帖的淘贴专辑推荐

vasp讨论帖子

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:

» 抢金币啦！回帖就可以得到:

查看全部散金贴
1楼2013-03-05 16:19:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

hypnos

新虫 (初入文坛)
★
lijuanshu(金币+1): 谢谢参与
the reason for this error message (DENTET) is that the Fermi level
cannot be determined accurately enough by the tetrahedron method (i.e, the integrated DOS stays different from the number of valence electrons)
please try (one of) the following:
--choose a different BZ integration scheme (ISMEAR)
--if you want or have to keep the tetrahedron integration,
-- increase the number of k-points (mind that Gamma has to be included in the k-mesh)
-- increase NEDOS

以上是vasp论坛对此类问题的回复

大意是费米能级不能精确确定。解决办法有如下三条：
1. 改变ISMEAR（如使用ISMEAR=0）
2. 增加K点
3. 增加NEDOS（DOS的采点个数）
赞一下(1人) 回复此楼

» 本帖已获得的红花（最新10朵）

lijuanshu
9楼2013-03-06 11:18:54
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lijuanshu

铜虫 (正式写手)
引用回帖:
14楼: Originally posted by hypnos at 2013-03-09 11:08:33
静态计算时候我个人觉得可以用ISMEAR=0
但是不同系统对于ISMEAR的设置是不同的，=0的设置对大多数系统都是有效的。SIGMA取值需要测试，用多个SIGMA计算单点能，取一个合理的entropy。具体做法我不太清楚了，不妨查 ...

我的体系时金属间化合物，算的是有关晶界的性质，在静态计算的时候ISMEAR是取0还是应该取2呢？
赞一下(1人) 回复此楼
15楼2013-03-09 11:58:55
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

171713294

木虫 (正式写手)
★
lijuanshu(金币+1): 谢谢参与
你这个静态计算的INCAR设置不合适吧
赞一下 回复此楼
2楼2013-03-05 17:30:07
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lijuanshu

铜虫 (正式写手)
引用回帖:
2楼: Originally posted by 171713294 at 2013-03-05 17:30:07
你这个静态计算的INCAR设置不合适吧

那些参数不合适，还希望你详细指出
赞一下 回复此楼
3楼2013-03-05 17:46:48
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

171713294

木虫 (正式写手)
加上NSW=0 IBRION=-1 关掉ISIF
赞一下 回复此楼
4楼2013-03-05 17:49:04
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lijuanshu

铜虫 (正式写手)
引用回帖:
4楼: Originally posted by 171713294 at 2013-03-05 17:49:04
加上NSW=0 IBRION=-1 关掉ISIF

但是NSW的默认值就是0啊，但NSW=0时，IBRION的默认值就是-1啊
赞一下 回复此楼
5楼2013-03-05 17:55:33
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lijuanshu

铜虫 (正式写手)
引用回帖:
6楼: Originally posted by yingying1588 at 2013-03-05 23:01:11

还希望您指点下，不胜感激
赞一下 回复此楼
7楼2013-03-06 08:48:32
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lijuanshu

铜虫 (正式写手)
怎么没人回复啊？这种警告对结果有影响吗？可以不用管它吗？
赞一下 回复此楼
8楼2013-03-06 09:12:37
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lijuanshu

铜虫 (正式写手)
送鲜花一朵
引用回帖:
9楼: Originally posted by hypnos at 2013-03-06 11:18:54
the reason for this error message (DENTET) is that the Fermi level
cannot be determined accurately enough by the tetrahedron method (i.e, the integrated DOS stays different from the number of valen ...

这三种方法都试了，只有第一种方法把ISMEAR改为0才可以消除这些警告？但是在静态计算的时候可以设置ISMEAR=0，SIGMA=0.2吗？期待你的回复，谢谢。
赞一下 回复此楼
11楼2013-03-06 14:29:38
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

future_wl

木虫 (著名写手)
★
lijuanshu(金币+1): 谢谢参与
你是不是用的VASP.4.6,编译器版本是ifort10 ?
赞一下 回复此楼
12楼2013-03-07 09:02:08
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lijuanshu

铜虫 (正式写手)
引用回帖:
12楼: Originally posted by future_wl at 2013-03-07 09:02:08
你是不是用的VASP.4.6,编译器版本是ifort10 ?

我用的是VASP.4.6，编译器版本是ifort9
赞一下 回复此楼
13楼2013-03-07 09:41:24
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hypnos

新虫 (初入文坛)
★
小木虫: 金币+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值
赞一下 回复此楼
14楼2013-03-09 11:08:33
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hypnos

新虫 (初入文坛)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
金属化合物如果是绝缘体，有带隙，就用0
如果是导体，考虑用2
赞一下 回复此楼
16楼2013-03-10 10:45:54
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

umklapp

银虫 (小有名气)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
你算的金属吧。 当你有bands过ferim level， ISMear=-5 不work
赞一下 回复此楼
18楼2013-03-10 15:53:23
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lijuanshu

铜虫 (正式写手)
引用回帖:
18楼: Originally posted by umklapp at 2013-03-10 15:53:23
你算的金属吧。 当你有bands过ferim level， ISMear=-5 不work

我算的体系就具有金属性质的
赞一下 回复此楼
19楼2013-03-11 08:26:37
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

暖暖67

新虫 (初入文坛)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
引用回帖:
11楼: Originally posted by lijuanshu at 2013-03-06 14:29:38
这三种方法都试了，只有第一种方法把ISMEAR改为0才可以消除这些警告？但是在静态计算的时候可以设置ISMEAR=0，SIGMA=0.2吗？期待你的回复，谢谢。...

我的就是这样设置的还是有问题 而且不是针对所有的点
赞一下 回复此楼
20楼2013-06-19 19:45:42
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
简单回复
yingying15886楼
2013-03-05 23:01   回复  
lijuanshu(金币+1): 谢谢参与
yanwuliaorao10楼
2013-03-06 11:36   回复  
lijuanshu(金币+1): 谢谢参与
lijuanshu17楼
2013-03-10 14:03   回复  
引用回帖:
16楼: Originally posted by hypnos at 2013-03-10 10:45:54 金属化合物如果是绝缘体，有带隙，就用0 如果是导体，考虑用2

相关版块跳转 我要订阅楼主 lijuanshu 的主题更新
201/1返回列表
普通表情 龙 兔 虎 猫 高级回复(可上传附件) SCI英文润色翻译-EnPapers美团队-先修改后付款
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 审不上青基又非升即走的青椒 和 牢里踩缝纫机的犯人哪个活的更舒服一点？ +19 非非飞远了 2024-05-20 22/1100 2024-05-23 21:09 by chengmy19
[基金申请] 基金评审 +5 阿呆不呆 2024-05-20 5/250 2024-05-23 20:57 by 3001160025
[考博] 大扩招 +13 version211 2024-05-21 19/950 2024-05-23 20:46 by keyaner23
[基金申请] 化学B口多少分能上会呀 +7 WOWO159357 2024-05-22 15/750 2024-05-23 17:25 by shuguang8600
[有机交流] TsCl保护羟基为什么不反应 +7 853015158 2024-05-21 25/1250 2024-05-23 17:19 by 853015158
[论文投稿] 为什么有的影响因子高的期刊分区不高呢？ +7 安处一室 2024-05-21 7/350 2024-05-23 15:51 by 晓目崇
[基金申请] 国自然的面上项目，５个审稿人，５个B能上会吗？ 4+11 lancet0903 2024-05-20 33/1650 2024-05-23 14:07 by 天狼2019
[论文投稿] 因为邮箱故障，没有收到论文接受后进入生产的邮件 5+3 VCOISTNT 2024-05-22 5/250 2024-05-23 10:03 by bobvan
[论文投稿] word转成pdf之后公式里面的字体变了，正文字体没变。 +9 1255037206 2024-05-20 11/550 2024-05-23 05:54 by tjushede
[基金申请] bless bless bless bless bless bless +6 chenwenqnig 2024-05-19 7/350 2024-05-22 22:59 by 957083516
[论文投稿] 求能源类，故障诊断，机器学习结合的EI期刊。 36+3 wj0318 2024-05-19 5/250 2024-05-22 22:05 by 知识产权服务
[考博] 换导师 +16 是柠檬呀！ 2024-05-18 29/1450 2024-05-22 16:29 by oooooo?o
[考博] 想被211以上高校课题组接收 +11 风起沧澜 2024-05-16 13/650 2024-05-22 14:37 by keyaner23
[硕博家园] 博士复试，申请成绩复核，有机会翻盘吗？ +15 长海二声笑 2024-05-21 22/1100 2024-05-22 12:44 by 带甲三千
[基金申请] 又有两个中科院三、四区计算机领域SCI/SSCI期刊爆出问题了！ +4 瞬息宇宙 2024-05-19 6/300 2024-05-22 12:37 by yp_nupt
[复合材料] 关于导电胶 +3 wgx2000 2024-05-17 6/300 2024-05-22 09:42 by wgx2000
[基金申请] 这个模块怎么成了烧香拜佛的地方了 +7 shrz98 2024-05-18 7/350 2024-05-21 10:26 by lancet0903
[考博] 25年博士申请 +6 lixinmiao9 2024-05-18 6/300 2024-05-20 11:19 by 裴先生533
[考博] 考博跨专业可以跨多远 +5 竹外。。。 2024-05-18 6/300 2024-05-19 21:51 by 摩天思瑞
[硕博家园] 五氯化铌怎么溶解啊 +3 南南枝枝 2024-05-17 5/250 2024-05-17 11:37 by ad_fish
信息提示
关闭
请填处理意见
关闭