wannier拟合能带总是拟合不上_Wannier 90使用教程

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拟合的很好。

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