求助：恳请指教、金红石二氧化钛的bader电荷得到的价态非常不合理

2955778771

结构是正常的二维金红石氧化钛结构，比例为TiO2，其中一个钛被替换成镍


Bader电荷计算得到的ACF文件，钛为正二价，氧为负一价
（钛的价电子是9.803644，potcar里面是12。氧是的价电子是7.1，镍的价电子是8.6，+1.4价）
正常氧化钛，钛应该是正四价，氧是负二价，明显不符。

一些参数：（发现加u对电子数没有影响）
EDIFFG = -0.02
ENCUT = 450 eV
EDIFF = 1E-05
ISMEAR = 0
SIGMA = 0.2
LDAU = .T.   #开启加U
LDAUTYPE = 2   #d方法
LDAUL = 3 3 -1   #加u的轨道
LDAUU = 5 6 0    #u值
LMAXMIX = 4     #混合参数


GGA-PBE，k点G 2 1 1，OPT和SCF已经收敛
第一次电子步在算到100步的时候都无法收敛RMM: 200    -0.932994010508E+03   -0.20204E-02   -0.37063E-02  2314   0.163E-01
第二次电子步收敛RMM:  39    -0.952331171963E+03   -0.67018E-05   -0.96136E-05  1807   0.147E-02

请问为什么会不合理，应该怎么改正，这对后续进行吸附能，PDOS和对吸附物的电子转移数的计算有影响吗
恳请大佬指教

ChemG

严格区分“原子电荷”和“价态”，根本不是一回事儿。好好看看以下的博文
一篇深入浅出、完整全面介绍原子电荷的综述文章已发表！http://sobereva.com/714
使用Multiwfn通过LOBA方法计算氧化态 http://sobereva.com/362
况且直接把Ni塞到TiO2里，你掺杂了个四价镍啊？四价镍可是能直接氧化水的。

sobereva

除了楼上提到的外，下文也应该看了
使用Multiwfn结合CP2K计算晶体中原子的氧化态
http://sobereva.com/711（http://bbs.keinsci.com/thread-45876-1-1.html）

2955778771

感谢两位大佬，我大概理解了