Ⅰ. 构造Wannier函数通常要仔细调节以下三大参数：

总能带数Num_bands，你需要有足够多的DFT能带数，以至于能够覆盖住你所需要研究的能带以及覆盖你投影子具有投影成分的能带。
投影子projector。选择投影子之前，请彻底做一次能带成分分析，也就是画Fatband图。请确认费米面附件的轨道成分，把这些轨道作为投影子。
解纠缠窗口和Frozen窗口，参考说明书，多试!

Ⅱ. 质量判断：

最初的判断就是对比DFT的能带和Wannier函数插值后的能带。看看会不会相等，如果你看到了一些曲折的条纹，那么说明你的Wannier函数不够局域。
建议一开始不做最局域化计算，也就是说设置num_iter=0。检查*.wout文件中Wannier函数的spread。如果Spread太大，比如超过了两个晶格常数。就说明你构造过程中的三大参数没有选好，请好好调节三大参数。
- 改变投影子，有时候增加投影子可以使得轨道更加局域；很多人构造wannier函数的时候把core level也加进来，这个实在不是很有必要。
- 增大解纠缠的能量空间也会有帮助，有时候你需要关注的投影子在你设置的解纠缠能量空间之外有投影成分，这时候就需要增大解纠缠空间；如果解纠缠空间过于大，也不利于wannier函数的构造。
- 慎重调节Frozen窗口，一定要确保这个窗口内除了有投影子的成分外，没有其它轨道的成分。
对称性粗略判断，可以通过观察Wannier函数的Spread，粗略地判断轨道是否简并。这个对非磁soc体系比较重要，如果你发现自旋上轨道和自旋下轨道的spread不同，那么就表明你构造的Wannier函数没有时间反演对称性。因此在时间反演的k点没有简并，同时计算出来的Z2拓扑数也会是错的。

Ⅲ. Tips：

增加kmesh。增加kmesh的意义在于可以得到更远距离的hopping，对于金属体系比较重要。
不要轻易使用num_iter，一旦使用了它，你就没法控制Wannier函数会怎么演化了。
对于同一个原子，尽量不要设置两个相同角量子数但是不同主量子数的轨道，比如1s和2s同时。对于金刚石这种三维结构，使用sp3轨道，对于平面的那种，使用sp2轨道。
能不能覆盖住所有的能带你需要做能带成分分析，也就是分析能带轨道波函数的投影。对于vasp，在计算能带的INCAR中添加LORBIT=11，分析fat-band。这一步最重要，能帮你确定你所需要设置的能量窗口。
dis_froz_min的设置有时候也可以设置小点，这样可以包括更多的能带，以至于收敛更快。但是还是要使用前面所说的原则。
建议一开始设置num_iter=0，看看你构造的Wannier函数的spread怎么样。如果spread太大，那就调节这些参数，使得这个spread减小。如果你实在没辙了，那就使用num_iter吧。

举个Bi2Se3的例子，从wannier90.wout中截取出来的（只要搜索"Final State"就可以找到）。解释一下数据，括号内的是center，最后一列是spread，也就是WF的展宽。

FinalStateWF centre and spread    1  (-0.000040, -1.194745,  6.638646 )     3.00135581WF centre and spread    2 (  0.000038, -1.196699,  6.640059 )    3.05024704WF centre and spread    3  (-0.000032, -1.192363,  6.640243 )     3.03256302WF centre and spread    4  (-0.000086, -3.583414,  2.908040 )     3.00131291WF centre and spread    5 (  0.000047, -3.581457,  2.906587 )    3.05019607WF centre and spread    6  (-0.000033, -3.585864,  2.906443 )     3.03267256WF centre and spread    7  (-0.000001,  1.194527,  4.773338 )    2.09798868WF centre and spread    8 (  0.000003,  1.194538, 4.773336 )     2.05011802WF centre and spread    9  (-0.000037,  1.194536,  4.773327 )    2.04466361WF centre and spread   10 (  0.000006, -1.194384,  1.130261 )    2.00572846WF centre and spread   11  (-0.000018, -1.216986,  1.140267 )     1.96783297WF centre and spread   12 (  0.000007, -1.172216,  1.140684 )    1.95999088WF centre and spread   13 (  0.000011, -3.583770,  8.416406 )    2.00573957WF centre and spread   14  (-0.000002, -3.561169,  8.406398 )     1.96786947WF centre and spread   15  (-0.000007, -3.605960,  8.405979 )     1.95994538WF centre and spread   16 (  0.000086, -1.194737,  6.638626 )    3.00130962WF centre and spread   17  (-0.000047, -1.196693,  6.640080 )     3.05019277WF centre and spread   18 (  0.000033, -1.192286,  6.640223 )    3.03266922WF centre and spread   19 (  0.000040, -3.583406,  2.908021 )    3.00135911WF centre and spread   20  (-0.000038, -3.581452,  2.906608 )     3.05025038WF centre and spread   21 (  0.000032, -3.585788,  2.906424 )    3.03256640WF centre and spread   22 (  0.000001,  1.194548, 4.773330 )     2.09798869WF centre and spread   23  (-0.000003,  1.194537,  4.773332 )    2.05011802WF centre and spread   24 (  0.000037,  1.194539, 4.773340 )     2.04466360WF centre and spread   25  (-0.000011, -1.194381,  1.130260 )     2.00573073WF centre and spread   26 (  0.000002, -1.216981,  1.140268 )    1.96786075WF centre and spread   27 (  0.000007, -1.172191,  1.140687 )    1.95993647WF centre and spread   28  (-0.000006, -3.583766,  8.416405 )     2.00573730WF centre and spread   29 (  0.000018, -3.561165,  8.406400 )    1.96784170WF centre and spread   30  (-0.000007, -3.605935,  8.405982 )     1.95999981Sum of centres and spreads (  0.000000,-50.170584,143.200000 )  72.45644902

①首先看center，看看你的center是不是束缚在原子本身的位置上。这点可以和原子坐标做对比。原子坐标也可以在wannier90.wout中找到。在构造的Wannier函数中，可以发现Wannier函数还是很好的坐落在原子位置上。如果你的Wannier函数能够很好的坐落在原子位置上，那么也能反应出你的Wannier函数的对称性很好。

*----------------------------------------------------------------------------*
|  Site       Fractional Coordinate          Cartesian Coordinate (Ang)     |
+----------------------------------------------------------------------------+
|Bi   1  0.39900   0.39900   0.69700  |    0.00000  -1.19454  6.65403    |
|Bi   2  0.60100   0.60100   0.30300  |    0.00000  -3.58361  2.89264    |
|Se   1 -0.00000  -0.00000   0.50000  |    0.00000   1.19454  4.77333    |
|Se   2  0.20600   0.20600   0.11800  |    0.00000  -1.19454  1.12651    |
|Se   3  0.79400   0.79400   0.88200  |    0.00000 -3.58361   8.42016    |
*----------------------------------------------------------------------------*

②检查轨道的简并性问题，在Bi2Se3里头，px, py, pz能够近似简并。从spread中能够看出我们的Wannier函数符合这一条件。

③检查自旋的简并。也就是对比第n+1条和第n+15条的center和spread。这个直接决定了你的Wannier函数是否具有时间反演不变性。这里的Wannier函数的顺序决定于不同的软件包。VASP的顺序和QE的顺序是不同的，这个需要自己去看相关的说明书。

作者：水田水田
链接：https://www.jianshu.com/p/ed0c5fd0d14f
来源：简书
著作权归作者所有。商业转载请联系作者获得授权，非商业转载请注明出处。

wannier插值能带拟合1相关推荐

wannier插值能带拟合5
Wannier90是一个构建局域轨道基组的有力工具,一般用来计算紧束缚模型中的占位能和跳跃积分,再结合杂化泛函或GW方法,插值得到更为准确的能带结构,或者计算极化.输运.拓扑不变量等物理性质.在其输入 ...
wannier插值能带拟合3
今天老王来介绍一下VASP与wannier90结合构建基于第一性原理的紧束缚模型. 首先登录wannier90官网下载wannier90 2.1版本 Download - Wannier90www. ...
wannier插值能带拟合2
1) Wannier function的构造及其物理含义 2) Wannier function的应用: GW计算能带结构绘制 Thanks to Xiegang Zhu 认识1:第一性原理计算 ...
wannier插值能带拟合4
这里利用Wannier90得到Wannier轨道后, 通过其来计算一下对应的材料能带, 这个通常用来和vasp计算得到的能带来做对比, 二者在费米面附近通常是要吻合的非常好, 这样才能说明在Wanni ...
MATLAB-wannier90_hr.dat后处理以及紧束缚能带拟合（晶体规范）
初次接触vasp,在学习的过程中发现没有完整的关于wannier画能带的Matlab程序,经过三周的学习编写了后处理wannier90_hr.dat的Matlab代码.在这里参考了Supriyo Da ...
【机器学习入门到精通系列】插值与拟合
文章目录 1 插值! 2 一维插值拉格朗日插值法分段线性插值法例题 3 二维插值网络结点插值法例题散乱数据插值 4 拟合! 5 区别与联系 1 插值! 在离散数据的基础上补插连续函数,使得 ...
python 曲线拟合_Python实现数学模型(插值、拟合和微分方程)
这篇文章主要介绍了python实现数学模型(插值.拟合和微分方程),本文给大家介绍的非常详细,对大家的学习或工作具有一定的参考借鉴价值,需要的朋友可以参考下问题1 车辆数量估计题目描述交通管理部门 ...
数学建模-预测模型总结(适用范围、优缺点)【灰色预测模型、插值与拟合、时间序列预测法、马尔科夫预测、差分方程、微分方程模型、神经元网络】
基于数学建模的预测方法种类繁多,从经典的单耗法.弹性系数法.统计分析法,到目前的灰色预测法.当在使用相应的预测方法建立预测模型时,我们需要知道主要的一些预测方法的研究特点,优缺点和适用范围.下面就当下 ...
matlab插值与拟合例题_菜鸟进阶系列:MATLAB数学建模·数据插值与拟合
本篇算是开始数据处理了.首先明确俩概念:插值和拟合.两者最根本的共同点都是基于现有数据进行预测.推演,比如根据现有的天文观测数据预测天体位置.插值问题是数学史上的经典问题,拉格朗日.高斯.牛顿等著名数 ...

wannier插值能带拟合1

Ⅰ. 构造Wannier函数通常要仔细调节以下三大参数：

Ⅱ. 质量判断：

Ⅲ. Tips：

wannier插值能带拟合1相关推荐

最新文章

热门文章