24小时热门版块排行榜    

查看: 1712  |  回复: 5

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的总和不等于电子总数,差值为真空电荷值,而我查看以往算例真空电荷和真空体积应该为零。请问这是什么原因造成的呢?我该怎么修改相关设置呢?
回复此楼

» 猜你喜欢

已阅   回复此楼   关注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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 286求调剂 +6 lemonzzn 2026-03-16 9/450 2026-03-18 15:31 by vgtyfty
[考研] 08工科 320总分 求调剂 +5 梨花珞晚风 2026-03-17 5/250 2026-03-18 14:49 by haxia
[考研] 311求调剂 +6 26研0 2026-03-15 6/300 2026-03-18 14:43 by haxia
[考研] 298-一志愿中国农业大学-求调剂 +7 手机用户 2026-03-17 7/350 2026-03-18 14:34 by vgtyfty
[考研] 材料专硕274一志愿陕西师范大学求调剂 +6 薛云鹏 2026-03-13 6/300 2026-03-18 14:14 by 脱颖而出
[考研] 070300化学319求调剂 +6 锦鲤0909 2026-03-17 6/300 2026-03-18 13:22 by Iveryant
[考研] 302求调剂 +10 呼呼呼。。。。 2026-03-17 10/500 2026-03-18 12:45 by Linda Hu
[考研] 0703化学调剂 +3 妮妮ninicgb 2026-03-17 3/150 2026-03-18 10:29 by macy2011
[考研] 材料工程专硕调剂 +5 204818@lcx 2026-03-17 5/250 2026-03-17 17:27 by Little-xue
[考博] 26申博 +4 八旬速览 2026-03-16 4/200 2026-03-17 13:00 by 轻松不少随
[考研] 302求调剂 +4 小贾同学123 2026-03-15 8/400 2026-03-17 10:33 by 小贾同学123
[考研] 11408 一志愿西电,277分求调剂 +3 zhouzhen654 2026-03-16 3/150 2026-03-17 07:03 by laoshidan
[考研] [导师推荐]西南科技大学国防/材料导师推荐 +3 尖角小荷 2026-03-16 6/300 2026-03-16 23:21 by 尖角小荷
[考研] 0703化学调剂 +6 妮妮ninicgb 2026-03-15 9/450 2026-03-16 16:40 by houyaoxu
[考博] 东华理工大学化材专业26届硕士博士申请 +6 zlingli 2026-03-13 6/300 2026-03-15 20:00 by ryzcf
[考研] 22408总分284求调剂 +3 InAspic 2026-03-13 3/150 2026-03-15 11:10 by zhq0425
[考研] 080500,材料学硕302分求调剂学校 +4 初识可乐 2026-03-14 5/250 2026-03-14 21:08 by peike
[考研] 297一志愿上交085600求调剂 +5 指尖八千里 2026-03-14 5/250 2026-03-14 17:26 by a不易
[考研] 304求调剂 +6 Mochaaaa 2026-03-12 7/350 2026-03-13 22:18 by 星空星月
[论文投稿] 投稿问题 5+4 星光灿烂xt 2026-03-12 6/300 2026-03-13 14:17 by god_tian
信息提示
请填处理意见