1） Wannier function的构造及其物理含义
2) Wannier function的应用： GW计算能带结构绘制
Thanks to Xiegang Zhu
认识1：第一性原理计算是为了获得wannier函数，有了wannier函数后就可以不需要第一性原理，直接通过wannier函数得到电子结构，尤其是表面态信息。
认识2： num_wann，num_band, Projections参数设置
认识3：画能带时高对称性k点位置的确定，可在wannier90_band.gnu文件中找到
vi wannier90_band.gnu
set data style dots
set nokey
set xrange [0: 5.73583]
set yrange [-49.75490 :211.32533]
set arrow from 0.51430, -49.75490 to 0.51430, 211.32533 nohead
set arrow from 1.76040, -49.75490 to 1.76040, 211.32533 nohead
set arrow from 2.27470, -49.75490 to 2.27470, 211.32533 nohead
set arrow from 3.52080, -49.75490 to 3.52080, 211.32533 nohead
set arrow from 4.59941, -49.75490 to 4.59941, 211.32533 nohead
set arrow from 5.11371, -49.75490 to 5.11371, 211.32533 nohead
set xtics (" G " 0.00000," A " 0.51430," H " 1.76040," K " 2.27470," G " 3.52080," M " 4.59941," L " 5.11371," H " 5.73583)
plot "wannier90_band.dat"
#num_wann = 64 ! set to NBANDS by VASP
num_wann=8
【nm_bands数量必须大于num_wann数量，从计算的64条能带中选取8条(num_bands)，对这8条轨道采用wannier函数投影（wannier函数是实空间的函数）；exclude_bands不是必须的；Nbands必须大于num_wann，否则会存在误投影情形】
num_bands=8
# for GW uncomment
exclude_bands 9-64
Begin Projections
Si:sp3
End Projections
【 这是初始猜测轨道，必须有；对Si来说，1个硅sp3的共4条，两个硅原子8条；对于ScH3，Sc以spd计算，则一个Sc原子需1+3+5=9条轨道，一个氢原子1条轨道，所以一个分子式的ScH3需要9+3=12条轨道，一个单胞有两个分子式，则需要24条轨道，因此，
num_wann=24,
num_bands=30 ！(大于24即可)；
如下方式设置投影
Begin Projections
Sc:l=0;l=1;l=2
H: l=-0
End Projections】
dis_froz_max=9
dis_num_iter=1000
guiding_centres=true
网络摘录1：
/html/11/6573747.html
0.1 下载wannier90-1.2，并编译安装
在/download.html下载
wannier90-1.2
修改make.sys
如果装了ifort编译器
LIBS = -lmkl_intel_lp64 -lmkl_sequential -lmkl_core
然后 make all
0.2 重新编译VASP
之后修改vasp
LIB = -L../vasp.5.lib -ldmy \
../vasp.5.lib/linpack_double.o \
../wannier90-1.2/libwannier.a $(SCA)(LAPACK)(LAPACK)(BLAS)
CPP = $(CPP_) -DMPI -DHOST=\"LinuxIFC\" -DIFC \
-DCACHE_SIZE=4000 -DPGF90 -Davoidalloc \
-DMPI_BLOCK=8000 -Duse_collective -DscaLAPACK -DNGZhalf -DVASP2WANNIER90
之后重新make
分五步计算
1. step1：静态计算
INCAR
ISMEAR = 0
SIGMA = 0.05
2.得到虚轨道（导带上方的空轨道）
INCAR
ISMEAR = 0
SIGMA = 0.05
ALGO = Exact
NBANDS = 64
LOPTICS=.TRUE.
NEDOS=2000
3. GW计算
INCAR
ALGO = GW0 ; LSPECTRAL = .TRUE. ; NOMEGA = 50
NBANDS=64
LRPA = .FALSE.
LWANNIER90=.TRUE.
前三步运行vasp 【第三步运行VASP前，需写好wannier90.win参数（可先不写wannier90.win文件，试运行在第三步运行VASP，自动产生符合POSCAR和KPOINTS的wannier90.win文件，再修改该文件，运行一遍VASP），以使得运行VASP后产生 wannier90.amn, wannier90.mmn, wannier90.eig三个文件，供后续第四步的wannier90.x程序调用】
4. step4：运行 wannier90.x wannier90
wannier90.win 【第三步运行完成之后得到的参数文件】
num_wann=8 【参数含义？】
num_bands=8
# for GW uncomment
exclude_bands 9-64
Begin Projections 【开始轨道投影？】
Si:sp3
End Projections
dis_froz_max=9
dis_num_iter=1000
guiding_centres=true
begin unit_cell_cart
2.7150000 2.7150000 0.0000000
0.0000000 2.7150000 2.7150000
2.7150000 0.0000000 2.7150000
end unit_cell_cart
begin atoms_cart
Si 0.0000000 0.0000000 0.0000000
Si 1.3575000 1.3575000 1.3575000
end atoms_cart
mp_grid = 4 4 4
begin kpoints
0.0000000 0.0000000 0.0000000
0.2500000 0.0000000 0.0000000
0.5000000 0.0000000 0.0000000
........
0.2500000 -0.5000000 -0.2500000
-0.2500000 0.2500000 -0.5000000
end kpoints
运行 wannier90.x wannier90
5.---step5:修改wannier90.win参数，再次运行 wannier90.x wannier90；能带即为wannier90_band.dat
num_wann=8
num_bands=8
# for GW uncomment
exclude_bands 9-64
Begin Projections
Si:sp3
End Projections
dis_froz_max=9
dis_num_iter=1000
guiding_centres=true
# Bandstructure plot
restart= plot
bands_plot = true
begin kpoint_path
L 0.50000 0.50000 0.5000 G 0.00000 0.00000 0.0000
G 0.00000 0.00000 0.0000 X 0.50000 0.00000 0.5000
X 0.50000 0.00000 0.5000 K 0.37500 -0.37500 0.0000
K 0.37500 -0.37500 0.0000 G 0.00000 0.00000 0.0000
end kpoint_path
bands_num_points 40
bands_plot_format gnuplot xmgrace
begin unit_cell_cart
2.7150000 2.7150000 0.0000000
0.0000000 2.7150000 2.7150000
2.7150000 0.0000000 2.7150000
end unit_cell_cart
begin atoms_cart
Si 0.0000000 0.0000000 0.0000000
Si 1.3575000 1.3575000 1.3575000
end atoms_cart
mp_grid = 4 4 4
begin kpoints
0.0000000 0.0000000 0.0000000
0.2500000 0.0000000 0.0000000
0.5000000 0.0000000 0.0000000
0.2500000 0.2500000 0.0000000
0.5000000 0.2500000 0.0000000
..........
0.2500000 -0.5000000 -0.2500000
-0.2500000 0.2500000 -0.5000000
end kpoints
再次运行 wannier90.x wannier90；能带即为wannier90_band.dat
如果你有gnuplot，可以直接执行wannier90_band.gnu。
