shinkon 发表于 2023-12-29 17:32
卡老师，请教一下，我有个稀土掺杂的系统，某种情况下做了自洽计算之后，磁矩变为负的（绝对值没问题）但 ...

负的无所谓吧。NUPDOWN设置完的话，磁矩就会被固定在这个值上，如果非要得到总磁矩为正，我建议调整初猜以及密度混合参数上面来变成正的，但这不是必要的。

卡开发发 发表于 2021-9-25 18:14
1、其实直接ISPIN=2和MAGMOM一般就够了，结果一样其实也正常，就是谁效率高谁效率低的问题。PBE优化结构的 ...

卡老师，请教一下，我有个稀土掺杂的系统，某种情况下做了自洽计算之后，磁矩变为负的（绝对值没问题）但能量是对的，是否可以通过设置NUPDOWN强行把它定在自旋向上呢

ymeng 发表于 2021-12-4 11:11
发发老师，这个outcar中的magnetization是不是只代表了该原子上单电子数，那如果想知道原子上的总电子数 ...

outcar给的magnetization和charge应该都是基于Wigner-Sietz半径划分的，Bader划分是基于zero-flux面，主要看你要讨论什么问题，粗糙看看两者都可以，要求高的话可以考虑做DDEC charge。

卡开发发 发表于 2021-9-25 18:14
1、其实直接ISPIN=2和MAGMOM一般就够了，结果一样其实也正常，就是谁效率高谁效率低的问题。PBE优化结构的 ...

发发老师，这个outcar中的magnetization是不是只代表了该原子上单电子数，那如果想知道原子上的总电子数要怎么办呢？我总觉得bader有点不能说明情况。

卡开发发 发表于 2021-9-25 18:14
1、其实直接ISPIN=2和MAGMOM一般就够了，结果一样其实也正常，就是谁效率高谁效率低的问题。PBE优化结构的 ...

谢谢老师！根据您的建议我再试一试。

1、其实直接ISPIN=2和MAGMOM一般就够了，结果一样其实也正常，就是谁效率高谁效率低的问题。PBE优化结构的阶段，尤其是做变晶格计算，最好多拷贝几次CONTCAR到POSCAR进行优化。
2、参考4.
3、恰恰相反，平滑才是因为ISMEAR=0，可能还是要足够密集的k点和一定的展宽，展宽小了峰就显得比较尖。
4、和WS半径很有关系，不只是total charge，包括你的磁矩偏低很可能也是类似问题导致的。