VASP计算DOS相关问题求助

小螃蟹

各位研友，近来我用VASP进行态密度的计算分析，我看到有的帖子说ISPIN参数=1/ISPIN=2是一样的，我自己也算了一下，当我不考虑自旋时，将ISPIN参数设置为1得出的数据与考虑自旋ISPIN=2的数据进行对比，发现它们的横坐标能量几乎没有变化，纵坐标数据呈现2倍的关系，即ISPIN=1是ISPIN=2的2倍。这是否是因为我不考虑自旋时，所得的贡献度是自旋向上和自旋向下的叠加。

卡开发发

哪个帖子这么说的？闭壳层的体系ISPIN=1和ISPIN=2才是一样的，此时spin DOS=alpha电子DOS-beta电子DOS=0，total DOS=alpha DOS+beta DOS=2*alpha DOS，这是没问题的。开壳层的体系ISPIN=1会强制按照闭壳层处理，得到的根本都不是同一个结果，总态密度怎么可能一样？

小螃蟹

卡开发发 发表于 2019-7-29 11:47
哪个帖子这么说的？闭壳层的体系ISPIN=1和ISPIN=2才是一样的，此时spin DOS=alpha电子DOS-beta电子DOS=0，t ...

嗯 谢谢您的详细解答。 我是在百度上搜ISPIN=1和ISPIN=2的区别时，看到说两者是没有区别的，可能是我没有看全或者没理解吧。还有一个问题，就是您说的开壳层与闭壳层体系：有一个以上亚层未充满的电子组态称为开壳层组态，简称开壳层;每个亚层都充满的组态为闭壳层组态，简称闭壳层，我的计算的是染料分子吸附TO2周期性晶胞上体系的DOS，我不会判断这个是开壳层还是闭壳层。（知道定义但不会判断）

卡开发发

小螃蟹 发表于 2019-7-29 13:06
嗯 谢谢您的详细解答。 我是在百度上搜ISPIN=1和ISPIN=2的区别时，看到说两者是没有区别的，可能是我没有 ...

简单点说就是看磁序。TiO2没磁性，一般情况下吸附的有机分子，这种情况就是闭壳层的，此时用ISPIN=1能够节约计算时间。

小螃蟹

卡开发发 发表于 2019-7-29 14:03
简单点说就是看磁序。TiO2没磁性，一般情况下吸附的有机分子，这种情况就是闭壳层的，此时用ISPIN=1能够 ...

非常感谢您我就是计算的这种闭壳层体系，然后纵坐标的数据ISPIN=1是ISPIN=2时的二倍，那这就是您说的：此时spin DOS=alpha电子DOS-beta电子DOS=0，total DOS=alpha DOS+beta DOS=2*alpha DOS这种情况吧？那我做图时，用不用将ISPIN=1得到的数据除以2。

卡开发发

小螃蟹 发表于 2019-7-29 14:30
非常感谢您我就是计算的这种闭壳层体系，然后纵坐标的数据ISPIN=1是ISPIN=2时的二倍，那这就是您说的 ...

不用，DOS的纵坐标单位本来就是a.u.（arbitary unit）。

小螃蟹

卡开发发 发表于 2019-7-29 18:45
不用，DOS的纵坐标单位本来就是a.u.（arbitary unit）。

我计算了孤立状态下的染料分子及染料分子-TiO2周期性晶胞吸附体系的DOS及PDOS6，发现孤立状态下染料分子的DOS在0附近没有费米能级，这可能是什么原因造成的呢？（我计算了3个不同有机分子在孤立状态下的DOS，只有一个出现了上述情况，所以想请教一下。）

卡开发发

小螃蟹 发表于 2019-8-14 15:52
我计算了孤立状态下的染料分子及染料分子-TiO2周期性晶胞吸附体系的DOS及PDOS6，发现孤立状态下染料分子 ...

这个影响应该不大。