vasp 输入文件中的ISMEAR参数指的是波函数占据数目，但是这个到底是什么意思？可以浅显一点讲吗？

就是说侍带电子在费米面附近占据数从0突变到1，这是个deta函数，为了计算方便，用一个平滑点的函数在费米面附近代替这个deta函数，这样计算就不容易振荡，易于收敛。和纳ismear就是唤谈没控制这种平滑函数的。

本来一开始打算用QE来算的，结果QE新版本bandup都不支持，banduppy例子都报错无语子，还是VASP吧。

他这里给了四步作为参考，但是实际上主要就是跑一个均匀网格点的自洽，和一个高对称点上的非自洽，然后做反折叠

首先第一步得到超胞CHGCAR

输入文件参考BANDUP的example的第一步，提交任务之后得到 CHGCAR 。

然后第二步主要是为了得到超胞非自洽计算的K点路径，这一步需要分别给出单胞和超胞的晶格坐标文件 prim_cell_lattice.in 与 supercell_lattice.in 。还有单胞的高对称K点路径

命令是 bandup kpts-sc-get -no_symm_avg

这一步会生成 KPOINTS_supercell.out 文件后面碧搜山会用到。

-no_symm_avg 这个标签是为了减少需要计算的k点个数，减小计算量和WAVECAR大小。

然后是第三步

这一步利用前面第一步生成的电荷密度来计算高对称点的波函数

INCAR为

ln -s ../step_1/CHGCAR CHGCAR

POSCAR POTCAR与第一步的一样，K点从第二步那个复制过来

cp ../step_2_*/KPOINTS_supercell.out KPOINTS

然后提交任务得到波函数WAVECAR

第四步就要开始做能带反折叠了，这一步其实就是一句命令，但是例子没有给它是怎么设置前边得到的波函数的信息的。这句命令是

会生成一个 unfolded_EBS_not-symmetry_averaged.dat 文件，画出来就行了。

Quantum Espresso 其实步骤也是差不多的，不同的地方主要在于首先bandup读取波函数悔中的方式还是QE老的格式，QE在6.4开始采用了全新的存储方式，在此之前的版本也必须编译的时候加上一个标签

configure --disable-xml

而之后的版本有一个人开发了 banduppy 声称能够解决问题，我还没试过。

另外不同的地方就在于第四步需要加一个指定prefix和outdir的标漏握签

---