计算化学公社

标题: 求助：vasp计算DOS，积分电子数（IPDOS）为什么不等于原子的电子数 [打印本页]

作者
Author: liheng 时间: 2023-12-18 11:11
标题: 求助：vasp计算DOS，积分电子数（IPDOS）为什么不等于原子的电子数
本帖最后由 liheng 于 2023-12-18 11:12 编辑

最近在算某一个体系的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的巨大差距；还有其它的结果就不一一列举了，反正就是不相等，而且差距还挺大）

作者
Author: 卡开发发 时间: 2023-12-18 11:43
和划分方式有关，VASP的PDOS用的WS划分使用的WS半径不一定能保证归一化。除去展宽的影响，实际上PDOS定积分到费米能级应当等于LOBRIT=11后面打印出来的值。

作者
Author: liheng 时间: 2023-12-18 12:13

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

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

作者
Author: 卡开发发 时间: 2023-12-18 12:55

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

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

作者
Author: liheng 时间: 2023-12-18 13:10

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

好的，谢谢卡开发发老师

作者
Author: Sheeemong 时间: 2025-7-18 02:27

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

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

作者
Author: zzyy3 时间: 2026-3-30 21:14

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

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

作者
Author: 卡开发发 时间: 2026-4-18 12:43

zzyy3 发表于 2026-3-30 21:14
老师您好，我也有类似的问题，但我查看DOSCAR发现第三列（积分态密度）到费米能级的数值恰好等于价电子数 ...

Total DOS不涉及到划分问题，所以进行以Fermi能级为上限的求积过程应该主要受到展宽、求积步长以及包括你提到的浮点数精度等影响。

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