24小时热门版块排行榜    

查看: 1690  |  回复: 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 的主题更新
信息提示
请填处理意见