vasp在计算磁性的实例和讨论(6)

最小原胞。比如对铁晶体的铁磁状态，你可以采用bcc的原胞来计算，但是在进行反铁磁的Fe计算，这是你需要采用sc的结构来计算，计算的晶胞中包括两个原子，你要设置一个原子的MAGMOM为正的，另一个原子的MAGMOM设置为负，但是它们的绝对值一样。因此在进行反铁磁的计算时，应该确定好反铁磁的磁胞，以及磁序，要判断哪种磁序和磁胞是最可能的反铁磁状态，那只能是先做好各种可能的排列组合，然后分别计算这些可能组合的情况，最后比较它们的总能，总能最低的就是可能的磁序。同样也可以与它们同铁磁或顺磁的进行比较。了解到该材料究竟是铁磁的、还是顺磁或反铁磁的。

NUPDOWN 设置为一个为非0而且非-1的数，就是相当是将在自洽迭代循环的计算过程加了一个限制，即限制自旋向上和自旋向下的电子数之差为一个固定的值，对应的就是给体系设置一个固定的总磁矩。具体在程序里面实现起来，量子化学的程序(比如Gaussian)和广泛用于固体物理的第一原理程序(比如WIEN2K)略有不同。

亚铁, 结构, 原子, 磁矩

磁结构，亚铁磁结构。ISPIN=2，

MAGMOM=......，这样来说，计算铁磁结构只需要把各个原子的磁矩方向设为一致，如，都为正值就可以了；计算反铁磁结构只需要把不同子子晶格的原子磁矩方向设为相反，如，一正一负就可以了；计算亚铁磁结构或许还要用到NUPDOWN。

但是有个问题，我做了一个亚铁磁结构的计算，初始MAGMOM设置了有正有负的原子磁矩，没有设置NUPDOWN值，而优化出来的结果里各个原子的磁矩却是按同一方向排列的！也就是说没有得到我想要的亚铁磁结构。如果要设置NUPDOWN进行强制总磁矩的话，又不知道这个值该设成多大？

非常希望各位多给意见，只要是有相关内容的回帖都送金币。 收藏 分享 评分 回复 引用

TOP

定义: 自旋多重度（spin multiplicity）＝2S+1, S=n*1/2,n为单电子数。所以，关键是单电子的数目是多少。

当有偶数个电子时，例如 O2，共有16个电子，那么单电子数目可能是0，即8个alpha和