计算Wilson loop

准备好wannier计算后的wannier90_hr.dat, wt.in,运行以下命令：
nohup wt.x &
wt.in文件如下：

&TB_FILE
Hrfile = 'wannier90_hr.dat'
Package = 'VASP'
/

&CONTROL
Z2_3D_calc            = T
/

&SYSTEM
NumOccupied =24
SOC = 1                 ! soc
/

&PARAMETERS
Nk1 = 51          ! number k points
Nk2 = 501          ! number k points
/

LATTICE
Angstrom
     4.1960747270310463    0.2582250830477404    0.0000000000000000
    -1.8744078817023773    3.7630198512971482    0.0000000000000000
     0.0000000000000000    0.0000000000000000   20.0000000000000000

ATOM_POSITIONS
4                              ! number of atoms for projectors
Direct                          ! Direct or Cartisen coordinate
Sb  0.3339682286449707  0.6650317843550384  0.5413060536727565
Sb  0.6650317843550286  0.3339682286449621  0.4586939463272437
F   0.4953343945142382  0.5036656184857289  0.3792253454505144
F   0.5036656184857613  0.4953343945142714  0.6207746245494831

PROJECTORS
4 4 4 4        ! number of projectors
Sb  s px py pz
Sb  s px py pz
F   s px py pz
F   s px py pz


MILLER_INDICES
1 0 0


MILLER_INDEX        ! this is equal to the SURFACE card
1 0 0