计算化学公社

标题: 求助：恳请指教、金红石二氧化钛的bader电荷得到的价态非常不合理 [打印本页]

作者
Author: 2955778771 时间: 2024-10-28 17:05
标题: 求助：恳请指教、金红石二氧化钛的bader电荷得到的价态非常不合理
结构是正常的二维金红石氧化钛结构，比例为TiO2，其中一个钛被替换成镍

(, 下载次数 Times of downloads: 1)
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和对吸附物的电子转移数的计算有影响吗
恳请大佬指教

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

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

作者
Author: 2955778771 时间: 2024-10-28 19:19
感谢两位大佬，我大概理解了

欢迎光临计算化学公社 (http://bbs.keinsci.com/)