wannier拟合能带总是拟合不上_科学网-Wannier90输入文件中num_wann, num_bands, 和energy window等参数设置规则-李云海的博文...

基本规则

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]。为了解纠缠,

  • 0
    点赞
  • 5
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
平面波展开法是计算周期性结构(如晶体)能结构的一种常用方法。对于声子晶体,可以将其看作具有周期性结构的介质,因此也可以使用平面波展开法来计算其隙特性。 在二维声子晶体,声波的传播方向被限制在平面内,因此只需要考虑平面波在二维平面内的展开。假设晶体的基元周期为$a$,则可以将平面波表示为: $$ e^{i\vec{k}\cdot\vec{r}}=e^{ik_xx}e^{ik_yy} $$ 其$\vec{k}$为波矢,$x$和$y$为晶体平面内的坐标。将平面波代入声子晶体的动力学方程,可以得到一个本征值问题,其解给出了声子晶体的能结构。 在实际计算,需要对波矢$\vec{k}$进行离散化,即将其分解为$k_x=2\pi n_x/L_x$和$k_y=2\pi n_y/L_y$,其$L_x$和$L_y$为晶体的尺寸,$n_x$和$n_y$为整数。然后,可以将平面波展开为: $$ e^{i\vec{k}\cdot\vec{r}}=\sum_{n_x,n_y}c_{n_x,n_y}e^{i\frac{2\pi}{L_x}n_xx}e^{i\frac{2\pi}{L_y}n_yy} $$ 其$c_{n_x,n_y}$为系数,需要通过求解本征值问题来确定。将展开后的平面波代入动力学方程,可以得到一个矩阵本征值问题,其解给出了声子晶体的能结构和隙特性。 需要注意的是,由于离散化导致的误差和计算量的增加,平面波展开法在计算大尺寸的声子晶体时可能会遇到困难。因此,一些改进方法,如Wannier函数和投影算子方法,已经被提出来用于加速计算。

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值