求助：vasp计算开soc的MAE(磁各向异性能)以及能带的参数疑惑

kkk99999

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

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

山有辰桉

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

Ivy_cao

请问，第二张图noncolinnear的计算不需要加SAXIS指定自旋旋转方向吗？

kkk99999

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

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

Ivy_cao

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方向上。请教下你知道是什么问题吗？或者我都做错了，请教下应该怎么设置这几个参数呢？期待回复，谢谢！

yurainingak

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