基本规则
Wannier90是一个构建局域轨道基组的有力工具,一般用来计算紧束缚模型中的占位能和跳跃积分,再结合杂化泛函或GW方法,插值得到更为准确的能带结构,或者计算极化、输运、拓扑不变量等物理性质。在其输入文件中,num_wann、num_bands、projections和energy window是几类极其重要的参数,如果设置不当,容易导致Wannier函数不够局域,计算结果误差很大。对于比较复杂的体系,这些参数的默认值一般不符合要求,需要手动指定。在这里总结一下个人关于如何设定这些参数的经验,欢迎大家批评指正。
在使用最大局域化Wannier函数(MLWF)理论及Wannier90程序之前,一定要先确定打算拟合哪一部分能带。拟合全部能带一来非常耗时,二来容易拟合失败。一般来说,如果做能带插值,取费米能级附近几条能带就足够了。如果计算输运、拓扑性质等物理量,可能有额外的判定标准。我没有算过这些性质,不敢妄下结论。
明确拟合哪一部分能带后,还需要结合fatband或者PDOS,确定这部分能带由哪些单原子态贡献。在这里以某化合物为例,其原胞中包含28个A原子和1个B原子,最高的价带为第49条能带,最低的导带为第50条能带。通过fatband得知,价带顶和导带底附近的能带主要由A原子的pz轨道和B原子的dz2和dx2-y2轨道贡献。有了这些信息,便可指定projections了:begin projections
A:pz
B:dz2;dx2-y2
end projections
Wannier函数个数为28*1+1*2=30。均分之后,各考虑15条价带和15条导带 [1]。为了解纠缠,