[zz]限制磁矩方向的计算

2009-11-25 16:10

天应要求又测试了Constraining the direction of magnetic moments（限制磁矩方向的计算） 和fully non-collinear magnetic structure calculations（完全非共线磁结构的计算）的关系。也就是让I_CONSTRAINED_M = 1 ，而分别计算LNONCOLLINEAR= .TRUE. 和.FALSE.的情况。发现在LNONCOLLINEAR=.FALSE.时的计算，得不到任何有关磁矩的数据（这些计算中LAMBDA = 0）。只有在LNONCOLLINEAR= .TRUE.时，才能得到有关磁矩的信息。也就是说，要进行限制性磁矩计算时，须让LNONCOLLINEAR= .TRUE.。

今天做对BCC结构的Ni，用SC晶格计算了Ni的磁矩(VASP程序，PAW+GGA)，有几点需记下来的，初始磁矩大小（差别的不是很大的话）对最后程序收敛得到的磁矩没有多大的影响。k mesh大小对总磁矩（我参看的是OSIZCAR文件中给出的mag值，我想一般还是在计算时在INCAR中设置RWIGS值，看每个原子的局域磁矩）影响比较大。

奇怪的一点，我用bcc晶格计算得到的磁矩要比sc晶格计算得到的磁矩小（因为我采用的同样的k mesh，导致对不同格子得到的不可约布里渊区中的k点数目不同，前者的较少，而后者的较多，所以导致了磁矩的差别）

在打开了I_CONSTRAINED_M=1或2后，则开始进行限制磁矩方向的计算（当然还有好几个其他的相关参数需设置），发现I_CONSTRAINED_M = 1和I_CONSTRAINED_M = 2给出的结果是一样。但是从源代码来看两者似乎得到的结果有所不同的。奇怪！（不知道为什么！）

另外如果对LAMBDA = 0，则应该是类似通常的自旋极化的计算。但是得到的磁矩又与spin-polarized calculations得到的磁矩又差别。奇怪（还不知道为什么！）