求助：vasp计算DOS，积分电子数（IPDOS）为什么不等于原子的电子数

liheng

最近在算某一个体系的DOS（在Mo表面上沉积H,O和Cs，然后发现好像有氢氧化铯产生，就想着计算一下H,O和Cs的DOS），由于是第一次计算，遇到一些问题，具体如下：
1，从图中可以看出，我画了H6，O1，Cs2，Cs3，Cs4和Cs6的PDOS（元素后面的数字表示编号，而非主量子数）。问题是为什么不同的铯元素之间的DOS纵轴大小差距这么大？
2，我用vaspkit输出的PDOS，同时也输出了IPDOS。我的体系有80个Mo，8个H，8个Cs和2个O，所以应该有80*14+8*1+8*9+2*6=1212个电子，可是ITDOS的结果却是1463.7（我自己对TDOS积分也是一样的结果），请问这是怎么一回事呢？
3，不光是ITDOS，所有的IPDOS的结果都不等于原子的电子总数，同类原子间的IPDOS差距也很明显？（例如，IPDOSH6的结果为0.27，不等于1；IPDOSO1的结果为2.82不等于6；IPDOSCs2的结果为12.28，大于9；IPDOSCs3的结果为0.97，远小于9；IPDOSCs4的结果为6.92，小于9；IPDOSCs6的结果为0.90，远小于9；不同Cs之间的IPDOS的差距也印证了图中DOS的巨大差距；还有其它的结果就不一一列举了，反正就是不相等，而且差距还挺大）

卡开发发

和划分方式有关，VASP的PDOS用的WS划分使用的WS半径不一定能保证归一化。除去展宽的影响，实际上PDOS定积分到费米能级应当等于LOBRIT=11后面打印出来的值。

liheng

卡开发发 发表于 2023-12-18 11:43
和划分方式有关，VASP的PDOS用的WS划分使用的WS半径不一定能保证归一化。除去展宽的影响，实际上PDOS定积分 ...

那请问卡开发发老师要如何解决这个问题呢？这样的结果应该不能写到文章里面吧？用cp2k计算会不会好一点

卡开发发

liheng 发表于 2023-12-18 12:13
那请问卡开发发老师要如何解决这个问题呢？这样的结果应该不能写到文章里面吧？用cp2k计算会不会好一点

你可能需要自己手动调整WS半径来保证覆盖的体积为100%，你可以参考关键字RWIGS里面的具体说法。
理论上用VASP的话WS半径调整到合理应该也可以，当然你要用其他程序也行，大部分固体计算程序可能是用Lowdin和Mulliken的比较多，各自有各自的问题。

liheng

卡开发发 发表于 2023-12-18 12:55
你可能需要自己手动调整WS半径来保证覆盖的体积为100%，你可以参考关键字RWIGS里面的具体说法。
理论上 ...

好的，谢谢卡开发发老师

Sheeemong

liheng 发表于 2023-12-18 13:10
好的，谢谢卡开发发老师

请问问题解决了吗，还想问一下如何正确选择WS半径

zzyy3

卡开发发 发表于 2023-12-18 11:43
和划分方式有关，VASP的PDOS用的WS划分使用的WS半径不一定能保证归一化。除去展宽的影响，实际上PDOS定积分 ...

老师您好，我也有类似的问题，但我查看DOSCAR发现第三列（积分态密度）到费米能级的数值恰好等于价电子数目，官网说DOSCAR里DOS一栏是通过积分态密度差分计算得到的，所以有没有可能是vaspkit和vasp对态密度的积分方式不同或者是小数点位数带来的影响呢？DOSCAR中积分态密度通常写成0.2589E+03相当于只记录到第一位小数，后面的小数位数全都丢掉了