计算化学公社

标题: 求助：vasp计算开soc的MAE(磁各向异性能)以及能带的参数疑惑 [打印本页]

作者
Author: kkk99999 时间: 2025-8-2 09:30
标题: 求助：vasp计算开soc的MAE(磁各向异性能)以及能带的参数疑惑
本帖最后由 kkk99999 于 2025-8-2 09:33 编辑

最近老师让计算MAE以及能带，两者均开soc，体系是单层MnBr,我们老师让计算MAE不是采用vasp官方的自洽完读非自洽，而是直接读取自洽能量即可，我遇到的问题就是（1）在进行非共线计算时需不需要加上ISYM = 0/-1,这个参数如果加是-1还是0，以及原因。（2）我的INCAR参数计算能带以及MAE有没有要改进的地方基于我前面提到的MAE计算方法，因为目前此参数容易震荡。（3）要不要用GGA_COMPAT = .F.参数

图1为计算能带自洽步骤INCAR,图二为计算MAE INCAR

作者
Author: 山有辰桉 时间: 2025-9-14 16:03
本帖最后由山有辰桉于 2025-9-14 16:09 编辑

感觉可以把磁性原子个数在扩展一下，可以讨论下更多的磁基态

作者
Author: Ivy_cao 时间: 2025-11-18 23:31
请问，第二张图noncolinnear的计算不需要加SAXIS指定自旋旋转方向吗？

作者
Author: kkk99999 时间: 2025-11-20 17:15

Ivy_cao 发表于 2025-11-18 23:31
请问，第二张图noncolinnear的计算不需要加SAXIS指定自旋旋转方向吗？

我是直接更改Magmom的大小自己应该可以计算夹角

作者
Author: Ivy_cao 时间: 2025-11-21 20:53

kkk99999 发表于 2025-11-20 17:15
我是直接更改Magmom的大小自己应该可以计算夹角

比如我想就指定正交的x，y和z三个方向，算这三个方向的能量，我用magmom指定了每个原子的初猜磁矩，指定x方向用SAXIS = 1 0 0， y方向 SAXIS = 0 1 0， z方向 SAXIS = 0 0 1，具体的x方向的磁设置如下：
# noncollinear + SOC
LORBMOM = .TRUE.
LNONCOLLINEAR = .TRUE.
LSORBIT = .TRUE.
MAGMOM = 0 0 2 \
      0 0 0.5 \
      0 0 0.5 \
      0 0 0.5 \
      0 0 0.5 \
      0 0 0.5 \
      0 0 0.5 \
      0 0 0.5 \
      0 0 0.5 \
      0 0 0.5 \
      0 0 0.5 \
SAXIS = 1 0 0
GGA_COMPAT = .FALSE.

我发现SAXIS并没有起到作用，OUTCAR里的magnetization还是都在z方向上。请教下你知道是什么问题吗？或者我都做错了，请教下应该怎么设置这几个参数呢？期待回复，谢谢！

作者
Author: yurainingak 时间: 2025-11-21 21:15
本帖最后由 yurainingak 于 2025-11-21 21:17 编辑

https://www.vasp.at/wiki/SAXIS 建议看VASP官方的Wiki，里面有清晰的说明：“There are two methods to rotate the initial magnetization in an arbitrary direction: either by changing the initial magnetic moments MAGMOM or by changing SAXIS. Both methods should, in principle, yield exactly the same energy, but for implementation reasons, the second method might be more precise.”
还有官方的例子：https://vasp.at/wiki/Determining_the_Magnetic_Anisotropy

作者
Author: chenxc 时间: 2026-4-2 13:38

Ivy_cao 发表于 2025-11-21 20:53
比如我想就指定正交的x，y和z三个方向，算这三个方向的能量，我用magmom指定了每个原子的初猜磁矩，指定 ...

同学/老师，您好！我在计算MAE的时候也碰到了SAXIS不生效的情况，无论怎么设置SAXIS的方向，计算结果中只能得到z方向的磁矩，请问您后来有解决吗？如果有这方面经验还请指导回复一下，十分感谢！

欢迎光临计算化学公社 (http://bbs.keinsci.com/)