当用VASP 计算态密度时，采用LORBIT=1是，需要手动设置RWIGS，在INCAR里面的RWIGS是单位为A，
添写适当的RWIGS半径，使得用RWIGS计算得到的体积和整个体积近似相等，这个时候RWIGS半径设置最为正确。
V = 4/3 * pi *(sum [N_i * RWIG_i^3])
根据上面的公式，可以计算出得到的采用球近似时体系的体积，这个体积和体系真实体积之比可以从OUTCAR里面找到，采用命令 grep % OURCAR,得可以得到这个百分比；
所以我们知道体积V，原子数N，根据上面的公式就可以直接求出RWIGS，填入INCAR。
如转载请注明出处。
若是你要算PDOS，或是需要各个原子的电荷量和磁矩等，
可以设置
1. PAW情况下，LORBIT >10，此时无需RWIGS，这里NPAR不等于1，应该也能得到PDOS等。
2. 若LORBIT小于10，要得到PDOS，就计算时把 NPAR 设为 1。
LORBIT和RWIGS的设置，详见VASP说明书
阅读(7058) | 评论(2) | 转发(2) |
