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

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

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

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

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

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

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

好的，谢谢卡开发发老师

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

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

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

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

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