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]。为了解纠缠,还需多预留4条价带 [2],因此需要关注第31-64条能带。如果用QE生成波函数,在第二步计算时(calculation="bands")需指定nbnd = 64。如果用VASP生成波函数,非自洽计算时需指定NBANDS = 64。因为我们只关心第31-64条能带,所以需要把第1-30条能带排除。对应参数设置为:num_wann = 30
num_bands = 34
exclude_bands = 1-30
生成波函数后,通过输出文件或EIGENVAL,确定各能带能量范围如下(以区间表示,单位为eV,均未作任何平移):
第64条能带能量范围为[4.58, 5.31],第35条能带(对应num_wann下限)能量范围为[-3.70, -3.09],第31条能带(对应num_bands下限)能量范围为[-4.91, -3.94]。
解纠缠时涉及两个能量窗口:outer_window和frozen_window。若不考虑其物理意义,frozen_window取值范围无限制,outer_window内至少包含num_wann条能带。在此例中,dis_win_min要低于第35条能带能量最小值(-3.70 eV),否则报如下错误:
实际上一般将outer_window设置为用exclude_bands排除无关能带后,剩下的能带能量范围,两端再留些余量。此例中为第31-64条能带,对应参数设置为:dis_win_min = -5.0
dis_win_max = 5.4
frozen_window则设定为num_wann对应能带的能量范围,两端留些余量。此例中为第35-64条能带,对应参数设置为:dis_froz_min = -3.8
dis_froz_max = 5.4
如此设置,一般能得到较好的拟合结果。
补充说明
[1] 费米能级上下各取15条能带,是根据fatband大概数出来的。取14条价带和16条导带,或18条价带和12条导带均可,不必十分严格。
[2] 关于什么情况下需要解纠缠,Wannier90的帮助文档和MLWF的理论文献中有十分详细的描述。简单来说,一旦num_wann所指定的能带范围把一组完整的能带一分为二,就需要解纠缠。
以上图为例,当num_wann指定的能带范围包含完整的A组和(或)B组价带时,不需要解纠缠。一旦切到了某组价带,或者切到了导带C,就需要解纠缠。实际上,只有在拟合完整的半导体或绝缘体价带时,才不需要解纠缠。对于半导体的导带或金属能带,很难找到一组和其它能带完全分离的能带。这种情况下就需要解纠缠。
解纠缠需要num_bands > num_wann,即额外预留一部分能带。本例中预留了4条价带,这没有什么严格的标准,预留4条导带或2条价带+2条导带均可,也可预留更多能带。
转载本文请联系原作者获取授权,同时请注明本文来自李云海科学网博客。
链接地址:http://wap.sciencenet.cn/blog-2909108-1154273.html
上一篇:QE使用经验总结:(一)编译
下一篇:MobaXterm及OpenSSH配置说明
wannier拟合能带总是拟合不上_科学网—Wannier90输入文件中num_wann, num_bands, 和energy window等参数设置规则 - 李云海的博文...相关推荐
- wannier拟合能带总是拟合不上_科学网-Wannier90输入文件中num_wann, num_bands, 和energy window等参数设置规则-李云海的博文...
基本规则 Wannier90是一个构建局域轨道基组的有力工具,一般用来计算紧束缚模型中的占位能和跳跃积分,再结合杂化泛函或GW方法,插值得到更为准确的能带结构,或者计算极化.输运.拓扑不变量等物理性质 ...
- python斐波那契数列第四十项_科学网—不死神兔的繁衍生息——神奇的斐波那契数列 - 霍开拓的博文...
不死神兔的繁衍生息--神奇的斐波那契数列 • 故事得从西元1202年说起,话说有一位意大利青年,名叫斐波那契.在他的一部著作中提出了一个有趣的问题:假设一对刚出生的小兔一个月后就能长成大兔,再过一个月 ...
- python课程结课感悟_科学网—《互联网+引论与Python》课堂感想(七) - 张忆文的博文...
课程感想 文/2018级社会学 方向阳 在本学期的课程中,我受益良多,这种收获不仅只是课程的知识,更多的是关于学习的态度.下面我将从三个方面讲述对课程的感想感悟. 一.课堂模式 互联网+课程的模式是让 ...
- python服务器怎么获取时间_科学网—python 笔记及示例(三)从时间服务器读取时间 - 钱磊的博文...
时间同步是一个观测站必须考虑的问题.实际的时间同步可以通过GPS进行,不过考虑一下用时间服务器同步也是很有意思的. 附件中的例子最初是 Paul Rubin写的(http://www.nightson ...
- wannier拟合能带总是拟合不上_Wannier 90使用教程
Wannier90 使用教程 第一步,计算能带 (正常DFT计算).因为之后要跟拟合的能带对比. 第二步,首先做一次自洽,加上NBANDS(正常DFT自洽), 通过使用wannier90拟合能带. 第 ...
- 将serversocket 写在按钮事件中连接不上_看 Netty 在 Dubbo 中如何应用
作者:莫那 鲁道http://cnblogs.com/stateis0/p/9062171.html 目录 dubbo 的 Consumer 消费者如何使用 Netty dubbo 的 Provide ...
- ftl不存在为真_科学网—世界上第一个真正的超光速实验 - 张操的博文
世界上第一个真正的超光速实验 2015年2月初,我们发表了第一篇关于交流电可以超光速实验的论文[1].近3年来,在科学网上引起了很多学者的讨论和质疑,也引发了一些理论解释的新猜想.我在这里向参加讨论的 ...
- ei会议论文录用但不参加会议_科学网—推荐一个录用论文全部发表在EI期刊上的会议 - 徐庆征的博文...
会议网址 http://www1.tyust.edu.cn/yuanxi/yjjg/iwsis2011/iwsis2011.htm 重要日期: Submission Date: March 10, 2 ...
- matlab mic系数_科学网—最大信息系数 (Maximal Information Coefficient, MIC)详解(1) - 彭勇的博文...
最大信息系数 (Maximal Information Coefficient, MIC)详解(1) 四年前看过的一篇论文,当时还在组会上报告过,很确信当时把它弄懂了,由于当时是用机器学习的方法来做预 ...
最新文章
- pl/sql dev连接报错Access violation at address 67614F04 in module 'oranls11.dll'
- Django从理论到实战(part32)--外键删除
- Java后端开发需要具备哪些知识结构
- Spring2..5整合Ehacahe
- java 反射 类名_java – 从反射中获取字段的类名
- python自动发送邮件_Python自动发送邮件
- html让图片自动旋转360,html5 canvas 360图片旋转制作抽奖转盘代码
- Vue安装 devTool 时报错的解决办法
- 红米 刷机 android7.1,有没有红米note3全网通的android7.1刷机包
- SQL基础篇 (增 删 查 改)
- java使用代理请求https
- 资产配置那些事-标准普尔家庭资产象限图
- Vivado IP核fifo使用指南
- 什么是Autorun病毒?它的运作原理是什么?如何手工清除?
- html导出excel换行,ASP.NET 导出到Excel时保留换行的代码
- 学会python,妈妈再也不用担心我乱花钱了!
- Python程序员爬出百套美女写真集,同样是爬虫,他为何如此突出?
- 2021年湖南省高考体考成绩查询,2021年湖南体育专业考试成绩查询网址:http://jyt.hunan.gov.cn/...
- 我开发了一个女朋友陪聊系统【全天24小时推送问候,自动回复女友的微信消息】
- 使用three.js + geojson 完成广西地图的绘制(上篇)
热门文章
- cad隐藏图层命令快捷键_CAD隐藏当前图层和只显示当前图层的快捷键是什么快捷键?...
- 鸢(yuan)尾花分类任务
- 计算机二级考试公网入口和教育网入口,考试入口
- 本以为能躺着进华为,结果陆续收到京东/滴滴/爱奇艺offer的我迷茫了
- 隐马尔可夫模型HMM
- Maven学习(基础部分)
- 全球及中国辉光放电质谱仪(GDMS)市场商业模式与十四五投资战略规划研究报告2022-2028年
- 1336A - Linova and Kingdom
- linux 7 开启远程桌面,CentOS 7 开启VNC远程桌面
- 十字军之王3 mac中文版