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

vkituser

nature protocol的文章里有写对于这些体系，设置GET_IRREDUCIBLE_KPOINTS = .FALSE.即可。这样可行，但是因为是对全布里渊区采样，相当于auto kmesh的ISYM=-1关掉对称性，这样k点数量会很多，计算时间很长而且导致可能内存不够算不动。直接用vasp生成的IBZKPT应该更快，因为vasp已经根据对称性生成了不可约k点，但这样的话继续用vaspkit 263 处理出来的费米面是错的。不知道vaspkit以后有没有这方面的新功能，对soc和磁性材料情况下也能用到对称性，不用全关掉。比如我下面这个例子，体系有soc，k点是553。
这个是IBZKPT的结果。

1. Automatically generated mesh
   
2.       21
   
3. Reciprocal lattice
   
4.     0.00000000000000    0.00000000000000    0.00000000000000             1
   
5.     0.20000000000000   -0.00000000000000    0.00000000000000             3
   
6.    -0.20000000000000    0.00000000000000    0.00000000000000             3
   
7.     0.40000000000000    0.00000000000000    0.00000000000000             3
   
8.    -0.40000000000000   -0.00000000000000    0.00000000000000             3
   
9.     0.20000000000000    0.20000000000000    0.00000000000000             6
   
10.     0.40000000000000    0.20000000000000    0.00000000000000             6
    
11.     0.00000000000000    0.00000000000000    0.33333333333334             1
    
12.     0.00000000000000    0.00000000000000   -0.33333333333334             1
    
13.     0.20000000000000   -0.00000000000000    0.33333333333334             3
    
14.    -0.20000000000000    0.00000000000000   -0.33333333333334             3
    
15.     0.00000000000000    0.20000000000000    0.33333333333334             3
    
16.     0.00000000000000   -0.20000000000000   -0.33333333333334             3
    
17.     0.40000000000000    0.00000000000000    0.33333333333334             3
    
18.    -0.40000000000000   -0.00000000000000   -0.33333333333334             3
    
19.     0.00000000000000    0.40000000000000    0.33333333333334             3
    
20.     0.00000000000000   -0.40000000000000   -0.33333333333334             3
    
21.     0.20000000000000    0.20000000000000    0.33333333333334             6
    
22.    -0.20000000000000   -0.20000000000000   -0.33333333333334             6
    
23.     0.40000000000000    0.20000000000000    0.33333333333334             6
    
24.    -0.40000000000000   -0.20000000000000   -0.33333333333334             6
    

复制代码

这个是vaspkit生成结果，相当于vasp关掉soc的情况的IBZKPT。

1. G    5    5    3    0.050     # Parameters to Generate KPOINTS (Don't Edit This Line)
   
2.       12
   
3. Reciprocal lattice
   
4.     0.000000000000    0.000000000000    0.000000000000       1
   
5.     0.200000000000    0.000000000000    0.000000000000       6
   
6.     0.400000000000    0.000000000000    0.000000000000       6
   
7.     0.200000000000    0.200000000000    0.000000000000       6
   
8.     0.400000000000    0.200000000000    0.000000000000       6
   
9.     0.000000000000    0.000000000000    0.333333333333       2
   
10.     0.200000000000    0.000000000000    0.333333333333       6
    
11.     0.400000000000    0.000000000000    0.333333333333       6
    
12.    -0.400000000000    0.000000000000    0.333333333333       6
    
13.    -0.200000000000    0.000000000000    0.333333333333       6
    
14.     0.200000000000    0.200000000000    0.333333333333      12
    
15.     0.400000000000    0.200000000000    0.333333333333      12
    

复制代码

如果设置GET_IRREDUCIBLE_KPOINTS = .FALSE.那就是全部75个k点。

vaspkit

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

vkituser

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

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