求助VASP表面结构优化自旋极化ISPIN及MAGMOM设置问题

Adair

各位老师好，刚开始接触VASP计算不久，关于VASP计算过程中是否考虑自旋极化问题一直不是特别清楚，想请教各位老师，感谢各位老师的指导。现在在优化 BiOBr的（001）面，不清楚到底该加不加自旋极化，就先添加了自旋极化计算，优化完成后在OSZICAR最后输出的结果如下（文末最后一张图，不知道为什么调整不过来）：

文末最后一张图：考虑自旋极化时OSZICAR输出的结果
不知道对不对，就在大师兄VASP学习网站上查找了相关的资料，看到大师兄在帖子中说可以做一个简单的测试来判断是否应该添加自旋极化。于是就进行了相关的测试：
测试计算1：不开自旋极化 ISPIN=1
测试计算2：打开自旋极化 ISPIN=2，并且在INCAR中添加 LORBIT=11
根据大师兄所说，查看测试计算2的OUTCAR文件，如果末尾输出的体系原子磁矩不为０，那么就需要考虑自旋极化（如果体系中所有原子磁矩都为零，那么就不需要考虑自旋极化）。我的输出结果明显不是０（如图２），也就意味着我在计算的时候是要考虑自旋极化的。

图２：测试计算２中ISPIN=2 LORBIT=11计算输出的OUTCAR中原子磁矩




在Vaspwiki上看了添加自旋的时候需要同时设置MAGMOM，Vaspwiki上给出的MAGMOM设置方法如下：


图３：Vaspwiki中MAGMOM设置的格式




这里有个问题想请教各位老师：
问题１：根据上面写的MAGMOM＝各元素原子数＊原子磁矩，我想知道各原子的磁矩（上图红框里设置的数值）是如何确定的？是通过计算得来的吗？比如我这个体系中Bi，O，Br的磁矩该如何计算出来？
另外还有两个问题：
问题２：如果在计算的过程中打开自旋极化，那么我们如何判断计算结果的正确行？
问题３：对于一个给定的体系，我们如何判断计算的时候是否考虑添加自旋极化？在不知道的情况下是不是按上面提到的大师兄的方法可以加以判断？

loovfnd

同样有疑惑，楼主找到答案了吗？

oneseven

同样有疑问，希望楼主如果找到解决方法了可以更新一下，谢谢~

successjiang

您的问题我也有同样的疑惑，除此之外，我对自己的体系也有关于磁矩上面的困扰，我是在V2O5表面掺杂很多过渡金属元素，其中一些是典型带有磁性的原子比如Fe，Co等，但是他们自己不成一相，请问这个需要考虑磁矩吗？又该怎么确定磁矩是否正常？针对这样的体系我只是统一设置ISPIN=2，而没有关于初始磁矩有多余的设置，这样做可以吗？有没有大神能顺便解答一下小弟的疑问，谢谢。

山里面的和尚

请卡卡发发老师来解答一下@卡开发发

卡开发发

山里面的和尚 发表于 2021-8-7 15:40
请卡卡发发老师来解答一下@卡开发发

LZ的问题其实私下讨论过了，这里不再重复讨论。连同楼上的问题写在这吧。

通常情况下铁磁性体系如果磁矩本身也不大的不猜也就无所谓，或者无磁性的索性就关了自旋极化。但是如果是反铁磁、亚铁磁的很可能不猜磁矩会导致算出来是铁磁性的，或者落到其他的反铁磁的状态上（反铁磁的磁序可能性就比较多了）。目前来说，磁序特别复杂的情况，确实没什么好的办法去猜，

wprison

标记收藏一下这贴

jiyuanyuan

successjiang 发表于 2021-7-31 11:25
您的问题我也有同样的疑惑，除此之外，我对自己的体系也有关于磁矩上面的困扰，我是在V2O5表面掺杂很多过渡 ...

我也是掺杂了只加了ISPIN=2, 算出来不知道对不对

15162142191

successjiang 发表于 2021-7-31 11:25
您的问题我也有同样的疑惑，除此之外，我对自己的体系也有关于磁矩上面的困扰，我是在V2O5表面掺杂很多过渡 ...

请问你最后打开自旋了吗