当用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说明书
阅读(7234) | 评论(2) | 转发(2) |