1、计算磁性体系,首先要知道它的稳定态是什么,一般第一步要进行优化
2、分两步来看,第一,做非磁的优化和静态计算,参数同前面小编的设置
3、按照上面的设置做一个非磁的计算,先看一下能量。接下来做磁性的计算:INCAR中磁性的计算要用ISPIN和MAGMOM两个tag来设置:ISPIN=2 表示计算磁性MAGMOM 是每个原子的磁性,这个要手动去写
4、磁性第一步优化完成可以跟非磁的对比能量,能量低的为基态,如下图,结构
5、小编的例子是CrI3的结构,2x2的超胞,也就是有4个Cr原子,这里Cr是有磁性的,所以小编给加了4个磁性 3 3 3 3 如下图:12*0是指其余12个磁矩都是0
6、这就是铁磁态的计算,至于反铁磁,自己要先找好反铁磁结构,才能手动的在MAGMOM中设定磁矩,来计算不同的反铁磁情况
7、如上图,蓝色是Cr原子,可以看出它的排列是石墨烯结构,这里小编设定了四种不同的磁耦合情况,包括FM铁磁,和3种反铁磁
8、如上图右边,列出了四中磁耦合态
9、所以小编的四种计算中FMMAGMOM= 3 3 3 3 12*0AFMMAGMOM= 3 -3 -3 3 12*02AFMMAGMOM= 3 -3 3 -3 12*03AFMMAGMOM= 3 3 -3 -3 12*0
10、四种磁性态算完后,可以比较能量来看哪一种是基态。