临界应力欧拉公式适用条件，四态法计算磁耦合强度
临界应力欧拉公式适用条件，四态法计算磁耦合强度LWAVE =.TRUE.LCHARG = .TRUE.步骤2. 非线性磁的静态计算（IBRION=-1; NSW=0; NELM =300），固定磁矩方向INCAR中开启:python2.7 super.py POSCAR 3 3 1 >POSCAR_super然后复制超胞成POSCAR：cp POSCAR_super POSCAR步骤1. 非磁计算在进行非线性计算之前，将超胞作为输入结构文件，做一步非磁性的静态计算，输出WAVECAR 和 CHGCAR
上一节（此处可链接）介绍了四态法的原理，这一节介绍一下如何设置上一节中提到的四个态。
重要参数：1. 非线性磁计算；2. 固定磁矩方向
以CrI3为例，如下图，我们计算红色标出的两个Cr原子的磁耦合强度Jz(最近邻z方向自旋磁耦合)。为了保证这两个原子与超胞外的3 4之间的影响可忽略不计，我们构造3*3的超胞。
对于原包CrI3结构文件，我们可以用小编的脚本super.py来快速处理：
python2.7 super.py POSCAR 3 3 1 >POSCAR_super
然后复制超胞成POSCAR：
cp POSCAR_super POSCAR
步骤1. 非磁计算
在进行非线性计算之前，将超胞作为输入结构文件，做一步非磁性的静态计算，输出WAVECAR 和 CHGCAR
INCAR中开启:
LWAVE =.TRUE.LCHARG = .TRUE.
步骤2. 非线性磁的静态计算（IBRION=-1; NSW=0; NELM =300），固定磁矩方向
假设3*3超胞中的CrI3中， 我们考虑的两个近邻原子序号为1 2。共18个Cr，54个I
读取上一步的WAVECAR 和CHGCAR。分四个文件夹分别计算这四个态。
4个态的INCAR中，关于磁性的内容如下：
注意由于下面显示问题MAGMOM自动换行了，MAGMOM= 后面的数字要写在同一行里的。
态1 ：
LNONCOLLINEAR = .TRUE.
I_CONSTRAINED_M = 1
ISPIN = 2
MAGMOM = 0 0 3 0 0 3 0 0 3 0 0 3 0 0 3... ... 0 0 3 162*0
LDAU = .TRUE.
LDAUTYPE = 2
LDAUL = 2 -1
LDAUU = 2 0.0
LDAUJ = 0.0 0.0
态2.
LNONCOLLINEAR = .TRUE.
I_CONSTRAINED_M = 1
ISPIN = 2
MAGMOM = 0 0 3 0 0 -3 0 0 3 0 0 3 0 0 3... ... 0 0 3 162*0
LDAU = .TRUE.
LDAUTYPE = 2
LDAUL = 2 -1
LDAUU = 2 0.0
LDAUJ = 0.0 0.0
态3.
LNONCOLLINEAR = .TRUE.
I_CONSTRAINED_M = 1
ISPIN = 2
MAGMOM = 0 0 -3 0 0 3 0 0 3 0 0 3 0 0 3... ... 0 0 3 162*0
LDAU = .TRUE.
LDAUTYPE = 2
LDAUL = 2 -1
LDAUU = 2 0.0
LDAUJ = 0.0 0.0
态4.
LNONCOLLINEAR = .TRUE.
I_CONSTRAINED_M = 1
ISPIN = 2
MAGMOM = 0 0 -3 0 0 -3 0 0 3 0 0 3 0 0 3... ... 0 0 3 162*0
LDAU = .TRUE.
LDAUTYPE = 2
LDAUL = 2 -1
LDAUU = 2 0.0
LDAUJ = 0.0 0.0
L
抽奖
关注计算凝聚态物理公众号，回复：高量抽奖
邮寄地址仅限高校跟研究所，非高校研究所读者请者私信，上传自己的教师或者研究生证明。
M
预告
抱歉，本该把vasp算例分享给读者，但是小编的服务器最近很紧张，后面有时间会放出来。后续列出其他磁耦合参数的求解公式以及四态设置。对于磁耦合参数J会有对角项以及非对角项 (如DMI作用)，所以我们需要推导每一个J分量的四态设置（大致有9个分量）。
N
超胞脚本
链接:
https://pan.baidu.com/s/1ebOxyG80QO5Kit-RyodNXw
提取码: psu2
