论坛 › VASPKIT使用相关 › soc和磁性体系的费米面k点选择问题

vkituser 发表于 2025-12-6 12:27:59

soc和磁性体系的费米面k点选择问题

nature protocol的文章里有写对于这些体系，设置GET_IRREDUCIBLE_KPOINTS = .FALSE.即可。这样可行，但是因为是对全布里渊区采样，相当于auto kmesh的ISYM=-1关掉对称性，这样k点数量会很多，计算时间很长而且导致可能内存不够算不动。直接用vasp生成的IBZKPT应该更快，因为vasp已经根据对称性生成了不可约k点，但这样的话继续用vaspkit 263 处理出来的费米面是错的。不知道vaspkit以后有没有这方面的新功能，对soc和磁性材料情况下也能用到对称性，不用全关掉。比如我下面这个例子，体系有soc，k点是553。
这个是IBZKPT的结果。
Automatically generated mesh
      21
Reciprocal lattice
    0.00000000000000    0.00000000000000    0.00000000000000             1
    0.20000000000000   -0.00000000000000    0.00000000000000             3
   -0.20000000000000    0.00000000000000    0.00000000000000             3
    0.40000000000000    0.00000000000000    0.00000000000000             3
   -0.40000000000000   -0.00000000000000    0.00000000000000             3
    0.20000000000000    0.20000000000000    0.00000000000000             6
    0.40000000000000    0.20000000000000    0.00000000000000             6
    0.00000000000000    0.00000000000000    0.33333333333334             1
    0.00000000000000    0.00000000000000   -0.33333333333334             1
    0.20000000000000   -0.00000000000000    0.33333333333334             3
   -0.20000000000000    0.00000000000000   -0.33333333333334             3
    0.00000000000000    0.20000000000000    0.33333333333334             3
    0.00000000000000   -0.20000000000000   -0.33333333333334             3
    0.40000000000000    0.00000000000000    0.33333333333334             3
   -0.40000000000000   -0.00000000000000   -0.33333333333334             3
    0.00000000000000    0.40000000000000    0.33333333333334             3
    0.00000000000000   -0.40000000000000   -0.33333333333334             3
    0.20000000000000    0.20000000000000    0.33333333333334             6
   -0.20000000000000   -0.20000000000000   -0.33333333333334             6
    0.40000000000000    0.20000000000000    0.33333333333334             6
   -0.40000000000000   -0.20000000000000   -0.33333333333334             6
这个是vaspkit生成结果，相当于vasp关掉soc的情况的IBZKPT。
G    5    5    3    0.050   # Parameters to Generate KPOINTS (Don't Edit This Line)
      12
Reciprocal lattice
    0.000000000000    0.000000000000    0.000000000000       1
    0.200000000000    0.000000000000    0.000000000000       6
    0.400000000000    0.000000000000    0.000000000000       6
    0.200000000000    0.200000000000    0.000000000000       6
    0.400000000000    0.200000000000    0.000000000000       6
    0.000000000000    0.000000000000    0.333333333333       2
    0.200000000000    0.000000000000    0.333333333333       6
    0.400000000000    0.000000000000    0.333333333333       6
   -0.400000000000    0.000000000000    0.333333333333       6
   -0.200000000000    0.000000000000    0.333333333333       6
    0.200000000000    0.200000000000    0.333333333333      12
    0.400000000000    0.200000000000    0.333333333333      12
如果设置GET_IRREDUCIBLE_KPOINTS = .FALSE.那就是全部75个k点。

vaspkit 发表于 2025-12-23 21:22:41

考虑SOC需要在INCAR中设置ISYM=-1或0，等价于GET_IRREDUCIBLE_KPOINTS = .FALSE.。你提供的第一个IBZKPT并没有关闭对称性识别，计算SOC结果可能不合理。

vkituser 发表于 2025-12-31 22:08:40

vaspkit 发表于 2025-12-23 21:22
考虑SOC需要在INCAR中设置ISYM=-1或0，等价于GET_IRREDUCIBLE_KPOINTS = .FALSE.。你提供的第一个IBZKPT并 ...

明白了 感谢解答。官网看到说soc体系要测试比较关或者不关对称性的结果。以ISYM=-1为准。HSE用ISYM=-1的话会很慢，如果对比结果看不出区别那应该还是用对称性会快一点。

查看完整版本: soc和磁性体系的费米面k点选择问题