Wannier90 使用教程

第一步，计算能带 （正常DFT计算）。因为之后要跟拟合的能带对比。

第二步，首先做一次自洽，加上NBANDS（正常DFT自洽）, 通过使用wannier90拟合能带。

第二步的详细过程：（注意： 以下2 ,3步是合并在一步的。）

1. DFT自洽（并在INCAR中添加NBANDS，可以不加，之后运行与wannier90接口的vasp需要加上NBANDS，且需要与此步的NBANDS一致，从OUTCAR中读取即可）

2.运行与wannier90接口的vasp，INCAR中需要添加参数：LWANNIER90 = .TRUE.，读取第1步的CHGCAR文件，ICHARG = 11。（使用的k点和自洽的时候一样，注意不是能带计算的KPOINTS, 是自洽的）

3. 在第2步的同时构建wannier90.win 文件，需要设置如下参数：

num_wann=24  （看态密度在费米能级附近的贡献，之后确定所要拟合的轨道：比如算碳材料，费米能级处都是p轨道的贡献，原胞中有8个原子，因此，带数目=3*8=24. 3代表px,py,pz轨道。如s轨道也有贡献，则设为 4*8=32）

num_bands=30 （手册中说的比较清楚，可以参考一下）

exclude_bands 31-80 

Begin Projections

C:px;py;pz （对应num_wann你所要投影的轨道）

End Projections

dis_froz_max=6

dis_num_iter=1000

guiding_centres=true

运行vasp程序，计算。

计算完成之后运行命令：wannier90.x  wannier90

4. 上步完成之后，修改wannier90.win文件：添加高对称点信息，和一些想要输出的信息。（此步即为拟合能带）

#Bandstructure plot

restart   = plot

hr_plot  = true          #要输出的信息 wannier90_hr.dat （计算哈密顿量, 可用于之后的表面态， Fermic Arc， 等的计算，具体参考wannier tool教程。）

bands_plot = true  (输出能带数据 wannier90_band.dat)

begin kpoint_path

G 0.000   0.000   0.000         X  0.000   0.000   0.500

X 0.000   0.000   0.500         M -0.500   0.500   0.500   

end kpoint_path

注意：高对称点与DFT计算能带的高对称点相同。

上步计算完成之后再次运行命令：wannier90.x  wannier90 即可得到wannier90_band.dat, 跟PBE的结果对比，如果费米能级附近符合的很完美，即为拟合成功。（注意，这里的数据没有减费米能级，比较的时候要减一下）

如下图：费米能级附近与PBE拟合的很好。

以上内容为本人自制，禁止转载，欢迎转发：如有错误，欢迎指正。有问题可以下面评论，大家一起讨论。