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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 求助，bader电荷计算出现真空电荷不为零
61/1返回列表
查看: 1928  |  回复: 5
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

mingmbm

铜虫 (小有名气)

[求助] 求助，bader电荷计算出现真空电荷不为零

如果在算bader电荷时，VACUUM CHARGE和VACUUM VOLUME不为零，出现价电荷总和不等于总价电子数，请问这个情况该怎么处理呢
INCAR:
SYSTEM=BN
PREC=Accurate
ENCUT=520
ISTART=0
ICHARG=2
EDIFF=1e-5
EDIFFG=-0.01
ISMEAR=-5
SIGMA=0.2
LAECHG=.True.
ISPIN=2


POSCAR:
B1 N1
   1.00000000000000
     3.6253608600271261   -0.0000000000000000    0.0000000000000000
     0.0000000000000000    3.6253608600271261    0.0000000000000000
    -0.0000000000000000    0.0000000000000000    3.6253608600271261
   B    N
     4     4
Direct
-0.0000000000000000  0.0000000000000000 -0.0000000000000000
-0.0000000000000000  0.5000000000000000  0.5000000000000000
  0.5000000000000000  0.0000000000000000  0.5000000000000000
  0.5000000000000000  0.5000000000000000 -0.0000000000000000
  0.2500000000000000  0.2500000000000000  0.2500000000000000
  0.7500000000000000  0.7500000000000000  0.2500000000000000
  0.7500000000000000  0.2500000000000000  0.7500000000000000
  0.2500000000000000  0.7500000000000000  0.7500000000000000

KPOINTS file
0
Gamma-centered
10  10  10
0 0 0

算bader电荷并处理后的结果：

    #         X           Y           Z       CHARGE      MIN DIST   ATOMIC VOL
--------------------------------------------------------------------------------
    1    0.000000    0.000000    0.000000    1.266260     0.424519     2.809570
    2    0.000000    1.812680    1.812680    1.229569     0.424519     2.700226
    3    1.812680    0.000000    1.812680    1.211632     0.424519     2.651387
    4    1.812680    1.812680    0.000000    1.211632     0.424519     2.651387
    5    0.906340    0.906340    0.906340    5.913683     0.929174     7.448141
    6    2.719021    2.719021    0.906340    5.914086     0.915542     7.389263
    7    2.719021    0.906340    2.719021    5.916489     0.915542     7.404457
    8    0.906340    2.719021    2.719021    5.909063     0.908649     7.386821
--------------------------------------------------------------------------------
    VACUUM CHARGE:               3.4276
    VACUUM VOLUME:               7.2077
    NUMBER OF ELECTRONS:        32.0000
从结果中可以看出，CHARGE的总和不等于电子总数，差值为真空电荷值，而我查看以往算例真空电荷和真空体积应该为零。请问这是什么原因造成的呢？我该怎么修改相关设置呢？
回复此楼

» 猜你喜欢
1楼 2017-11-03 11:38:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

有何不可琼

银虫 (著名写手)
我觉得，CHGCAR之和不等于总量，是因为没有收敛，应该在INCAR中增加NGX等参数，不过这个问题，我也没有完全解决，你可以看下我发的帖子，

发自小木虫IOS客户端
赞一下 回复此楼
优秀的研究员是自己学出来的
2楼2017-11-03 12:25:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mingmbm

铜虫 (小有名气)
引用回帖:
2楼: Originally posted by 有何不可琼 at 2017-11-03 06:25:44
我觉得，CHGCAR之和不等于总量，是因为没有收敛，应该在INCAR中增加NGX等参数，不过这个问题，我也没有完全解决，你可以看下我发的帖子，

我设置过NGX和NGXF等参数，但结果还是一样的，现在相当于NGX等参数是默认值
赞一下 回复此楼
3楼2017-11-03 15:20:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

有何不可琼

银虫 (著名写手)
引用回帖:
3楼: Originally posted by mingmbm at 2017-11-03 15:20:32
我设置过NGX和NGXF等参数，但结果还是一样的，现在相当于NGX等参数是默认值...

再探索探索吧。

发自小木虫IOS客户端
回复此楼
优秀的研究员是自己学出来的
4楼2017-11-03 21:05:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zheng_meng

新虫 (正式写手)
请问楼主这个问题解决了吗，我的真空电荷也不是0
赞一下 回复此楼
5楼2017-11-27 08:54:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mingmbm

铜虫 (小有名气)
引用回帖:
5楼: Originally posted by zheng_meng at 2017-11-27 02:54:23
请问楼主这个问题解决了吗，我的真空电荷也不是0

bader 电荷处理程序的问题，用老版本的就没问题
赞一下 回复此楼
6楼2017-12-22 16:08:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 mingmbm 的主题更新
61/1返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[博后之家] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +3 ASdOkHsho7FD 2026-08-28 3/150 2026-08-29 03:23 by gy1nBQXYQJqL
[硕博家园] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +3 ASdOkHsho7FD 2026-08-28 3/150 2026-08-29 03:03 by gy1nBQXYQJqL
[基金申请] 我就是申请一个面上项目而已，这评审意见是按照杰青的条件评的吧？ +3 gouxfjh 2026-08-28 3/150 2026-08-29 01:01 by jurkat.1640
[基金申请] 为什么资助数各大高校都创新高，自己申请怎么就这么难 +12 Kittylucky 2026-08-27 13/650 2026-08-29 00:04 by superceng
[基金申请] 国自然面上复盘~欢迎讨论 (金币+15) +13 晴天加油 2026-08-26 14/700 2026-08-28 18:36 by huagongfeihu
[基金申请] 有没有仍没收到信息的 +3 德尚中行 2026-08-27 3/150 2026-08-28 15:24 by rongshuxia
[教师之家] 导师吐槽：我怎么摊上了这么个极品研究生！ +8 苏东坡二世 2026-08-23 8/400 2026-08-28 14:41 by lkforward
[基金申请] 基金不中，共勉 +11 eulota 2026-08-26 11/550 2026-08-28 14:22 by 火星超人xi
[基金申请] 麻烦专家们看看评委们的意见（F口面上） +4 gdd2018 2026-08-28 9/450 2026-08-28 14:22 by jklily
[基金申请] 怎么查啊 +6 huang1991js 2026-08-26 6/300 2026-08-28 08:42 by winsaint
[基金申请] 申请删除本帖 +6 lyz123lyz 2026-08-27 7/350 2026-08-27 17:31 by 宁静致远sy
[基金申请] 怎么看青基中了没有啊 +5 叶九微 2026-08-26 5/250 2026-08-27 10:35 by l_zh2008
[基金申请] 为什么 国际(地区)合作与交流项目 没有放榜？ 10+3 majunge000 2026-08-26 11/550 2026-08-27 08:42 by 北京莱茵编辑
[基金申请] 范进中举一文的中心思想 +9 炎黄贵胄 2026-08-22 10/500 2026-08-26 15:40 by semaglutide
[基金申请] 为什么国自然不能直接公布 +4 bjdxyxy 2026-08-26 4/200 2026-08-26 13:12 by qingmu1201
[基金申请] 系统进不去 +4 yanglien 2026-08-26 5/250 2026-08-26 11:10 by wenfengw83
[基金申请] 项目信息和经费信息在系统里都可以看到了 +6 wittyboy 2026-08-26 14/700 2026-08-26 10:55 by wittyboy
[基金申请] 今天务委会开完了，明天出结果吗 +19 angus9576 2026-08-25 23/1150 2026-08-26 10:03 by zp519
[职场人生] 学生撞见辅导员送外卖，第二天全班都沉默了 +3 鲸鱼融金_浙江_� 2026-08-22 3/150 2026-08-26 09:01 by zzuzxg
[基金申请] 国合现在查不到了吗？ +10 chengyan1220 2026-08-24 20/1000 2026-08-26 08:57 by peasantsprig
信息提示
关闭
请填处理意见
关闭