Example4: Copper的费米面和能带轨道特征

Detail: 超软赝势；k点4×4×4；5个以Cu为中心的d轨道，2个位于四面体间隙的s轨道

Input (见Wannier90官网):

copper.win        输入文件

copper.mmn        积分矩阵

copper.amn        Bloch态在轨道上的尝试投影

copper.eig        Bloch态的本征值 

运行过程：

1 执行 wannier90.x copper

2 画费米面，费米能级为12.2123 eV，操作见example2

3 画差值能带结构，操作见example3

begin kpoint_path 

G 0.00 0.00 0.00 X 0.50 0.50 0.00 

X 0.50 0.50 0.00 W 0.50 0.75 0.25 

W 0.50 0.75 0.25 L 0.00 0.50 0.00 

L 0.00 0.50 0.00 G 0.00 0.00 0.00 

G 0.00 0.00 0.00 K 0.00 0.50 -0.50 

end kpoint_path

4 添加参数：bands_plot_project = 6,7  #在能带图上投影两个s轨道

执行 wannier90.x copper

5 对copper_band_proj.gnu进行画图处理，得到投影轨道能带图。

用生成的文件画图会得到和教程中不一样的结果，图像会比较扁，可以将gnu文件中的

#set term postscript enhanced eps color lw 0.5 dl 0.5

取消注释，放在文件开头 （果然set term pdf很草率哈哈哈）。

修改： set pointsize 0.5

在splot质量前添加行：set output "out.pdf"

除此以外高对称点的位置会和实际线的位置发生偏移，只需要在gnu文件中删除多余的空格就行了。

结果：

由于我们在copper.win里面设置的num_wann=7，所以这边一共有7个Wannier function，括号里对应其中心坐标，最右边对应的是spread（延展范围？）。

根据坐标，可以看出6,7对应的是四面体空隙中的s轨道，1-5对应的是以Cu原子为中心的d轨道。5个d轨道中，1和4的范围相同，235的范围相同，对应的分别是t2g轨道和eg轨道。

放一张费米面的图，只有band6穿过了费米能级，所以只有band6有费米面：

下面3张图是投影轨道能带图，可以看到不同的轨道对能带的贡献是不同的。

当然，我们也可以改变能量窗口的位置，来进行MLWFs的构造，需要注意的是，要确保能量窗口中包含的能带数要大于等于我们设置的num_wann，更改能量窗口的同时得到的能带结构会发生一些变化，操作相同，就不重复了。放一张教程里的结果图：

新手上路，敬请指正。