VASP磁性计算总结篇

以下是从 VASP 在线说明书整理出来的非线性磁矩和自旋轨道耦合的计算说明.非线性磁矩计算：1）计算非磁性基态产生 WAVECAR 和 CHGCAR 文件。2）然后 INCAR 中加上ISPIN=2ICHARG=1 或 11 ！读取 WAVECAR 和 CHGCAR 文件LNONCOLLINEAR=。TRUE. MAGMOM= 注意：①对于非线性磁矩计算,要在 x, y 和 z 方向分别加上磁矩，如MAGMOM = 1 0 0 0 1 0 !表示第一个原子在 x 方向，第二个原子的 y 方向有磁矩② 在任何时候，指定 MAGMOM 值的前提是 ICHARG=2(没有 WAVECAR 和 CHGCAR 文件）或者ICHARG=1 或 11（有 WAVECAR 和 CHGCAR 文件），但是前一步的计算是非磁性的(ISPIN=1）。 磁各向异性能（自旋轨道耦合)计算：注意： LSORBIT=.TRUE。 会自动打开 LNONCOLLINEAR= 。TRUE。选项,且自旋轨道计算只适用于 PAW 赝势，不适于超软赝势。自旋轨道耦合效应就意味着能量对磁矩的方向存在依赖，即存在磁各向异性能（MAE），所以要定义初始磁矩的方向。如下：LSORBIT = 。TRUE。SAXIS = s_x s_y s_z （quantisation axis for spin)默认值： SAXIS=(0+,0，1)，即 x 方向有正的无限小的磁矩,Z 方向有磁矩。 要使初始的磁矩方向平行于选定方向,有以下两种方法：MAGMOM = x y z ！ local magnetic moment in x,y,zSAXIS = 0 0 1 ！ quantisation axis parallel to zorMAGMOM = 0 0 total_magnetic_moment ！ local magnetic moment parallel to SAXIS (注意每个原子分别指定）SAXIS = x y z ！ quantisation axis parallel to vector （x,y,z)，如 0 0 1 两种方法原则上应该是等价的,但是实际上第二种方法更精确。第二种方法允许读取已存在的WAVECAR（来自线性或者非磁性计算)文件，并且继续另一个自旋方向的计算（改变 SAXIS 值而MAGMOM 保持不变）。当读取一个非线性磁矩计算的 WAVECAR 时,自旋方向会指定平行于 SAXIS。 计算磁各向异性的推举步骤是：1）首先计算线性磁矩以产生 WAVECAR 和 CHGCAR 文件（注意加入 LMAXMIX）。2）然后 INCAR 中加入：LSORBIT = 。TRUE。ICHARG = 11 ! non selfconsistent run, read CHGCAR！或 ICHARG ==1 优化到易磁化轴,但此时应提高 EDIFF 的精度LMAXMIX = 4 ! for d elements increase LMAXMIX to 4， f： LMAXMIX = 6！ you need to set LMAXMIX alr...