基本参考cnqiang的教程Intel编译WannierTools (cndaqiang.github.io)
，唯一不同是Makefile，不知道为什么使用新版本的wannier tools里面的Makefile.intel-mpi报错，

mpif90 -fpp -DMPI -fpe3 -O3 -DARPACK -DINTELMKL -O3 -static-intel -nogen-interface -I/opt/ohpc/pub/apps/oneapi2022/mkl/2022.1.0/include -c module.f90
gfortran: error: unrecognized command line option ‘-fpp’; did you mean ‘-cpp’?
gfortran: error: unrecognized command line option ‘-fpe3’; did you mean ‘-fdec’?
gfortran: error: unrecognized command line option ‘-static-intel’; did you mean ‘-static-pie’?
gfortran: error: unrecognized command line option ‘-nogen-interface’; did you mean ‘-Wimplicit-interface’?
make: *** [module.o] Error 1

换成之前安装的老版本的Makefile之后就顺利编译了，放在这里备忘。注意把OBJ换成新版的。

OBJ =  module.o sparse.o wt_aux.o math_lib.o symmetry.o \
                 readHmnR.o inverse.o proteus.o \
       eigen.o ham_qlayer2qlayer.o psi.o unfolding.o rand.o \
                 ham_slab.o ham_bulk.o ek_slab.o ek_bulk_polar.o ek_bulk.o \
       readinput.o fermisurface.o surfgreen.o surfstat.o \
                 mat_mul.o ham_ribbon.o ek_ribbon.o \
       fermiarc.o berrycurvature.o \
                 wanniercenter.o dos.o  orbital_momenta.o \
                 landau_level_sparse.o landau_level.o lanczos_sparse.o \
                 berry.o wanniercenter_adaptive.o \
                 effective_mass.o findnodes.o \
                 sigma_OHE.o sigma.o Boltz_transport_anomalous.o \
                 2D_TSC.o \
        main.o


#OBJ =  module.o wt_aux.o math_lib.o symmetry.o readHmnR.o inverse.o proteus.o \
#       eigen.o ham_qlayer2qlayer.o psi.o unfolding.o \
#                ham_slab.o ham_bulk.o ek_slab.o ek_bulk_polar.o ek_bulk.o \
#       readinput.o fermisurface.o surfgreen.o surfstat.o \
#                mat_mul.o ham_ribbon.o ek_ribbon.o \
#       fermiarc.o berrycurvature.o \
#                wanniercenter.o dos.o  orbital_momenta.o \
#                 berry.o wanniercenter_adaptive.o \
#                effective_mass.o findnodes.o sigma_AHC.o \
#       main.o

# compiler
F90  = mpiifort -fpp -DMPI

INCLUDE = -I${MKLROOT}/include
WFLAG = -nogen-interface
OFLAG = -O3 -g -traceback -static-intel
FFLAG = $(OFLAG) $(WFLAG)
LFLAG = $(OFLAG)

# blas and lapack libraries
# static linking
LIBS = -Wl,--start-group ${MKLROOT}/lib/intel64/libmkl_intel_lp64.a \
        ${MKLROOT}/lib/intel64/libmkl_sequential.a \
        ${MKLROOT}/lib/intel64/libmkl_core.a -Wl,--end-group -lpthread -lm -ldl

# dynamic linking
# LIBS = -L/${MKLROOT}/lib/intel64 -lmkl_core -lmkl_sequential -lmkl_intel_lp64 -lpthread

main : $(OBJ)
        $(F90) $(LFLAG) $(OBJ) -o wt.x $(LIBS)
        cp -f wt.x ../bin

.SUFFIXES: .o .f90

.f90.o :
        $(F90) $(FFLAG) $(INCLUDE) -c $*.f90

clean :
        rm -f *.o *.mod *~ wt.x

wannier tools计算声子拓扑

使用wannier_tools计算拓扑声子，需要将力常数文件FORCE_CONSTANTS转化为哈密顿量文件phonopyTB_hr.dat，这个可以通过wannier_tools里面的utility中的脚本实现，需要注意的是，它里面的脚本是python2的，而且适配的phonopy版本比较低，1.11.8，试了高版本的会不间断的报错。
相比于低版本的phonopy，高版本里面可以自动指定PRIMITIVE_AXES = AUTO
但是这只对1.14.0之后的版本适用，我的方法是先使用高版本的phonopy输出转移矩阵（在计算声子谱的命令里加一个 -v），然后在band.conf文件中使用PRIMITIVE_AXIS手动指定转移矩阵，注意老版本的这个参数是PRIMITIVE_AXIS而不是PRIMITIVE_AXES，指定后面这个不管用，识别不了。。。