wannier插值能带拟合2
1） Wannier function的构造及其物理含义
2) Wannier function的应用： GW计算能带结构绘制
Thanks to Xiegang Zhu
认识1：第一性原理计算是为了获得wannier函数，有了wannier函数后就可以不需要第一性原理，直接通过wannier函数得到电子结构，尤其是表面态信息。
认识2： num_wann，num_band, Projections参数设置
认识3：画能带时高对称性k点位置的确定，可在wannier90_band.gnu文件中找到
vi wannier90_band.gnu
set data style dots
set nokey
set xrange [0: 5.73583]
set yrange [-49.75490 :211.32533]
set arrow from 0.51430, -49.75490 to 0.51430, 211.32533 nohead
set arrow from 1.76040, -49.75490 to 1.76040, 211.32533 nohead
set arrow from 2.27470, -49.75490 to 2.27470, 211.32533 nohead
set arrow from 3.52080, -49.75490 to 3.52080, 211.32533 nohead
set arrow from 4.59941, -49.75490 to 4.59941, 211.32533 nohead
set arrow from 5.11371, -49.75490 to 5.11371, 211.32533 nohead
set xtics (" G " 0.00000," A " 0.51430," H " 1.76040," K " 2.27470," G " 3.52080," M " 4.59941," L " 5.11371," H " 5.73583)
plot "wannier90_band.dat"
#num_wann = 64 ! set to NBANDS by VASP
num_wann=8
【nm_bands数量必须大于num_wann数量，从计算的64条能带中选取8条(num_bands)，对这8条轨道采用wannier函数投影（wannier函数是实空间的函数）；exclude_bands不是必须的；Nbands必须大于num_wann，否则会存在误投影情形】
num_bands=8
# for GW uncomment
exclude_bands 9-64
Begin Projections
Si:sp3
End Projections
【 这是初始猜测轨道，必须有；对Si来说，1个硅sp3的共4条，两个硅原子8条；对于ScH3，Sc以spd计算，则一个Sc原子需1+3+5=9条轨道，一个氢原子1条轨道，所以一个分子式的ScH3需要9+3=12条轨道，一个单胞有两个分子式，则需要24条轨道，因此，
num_wann=24,
num_bands=30 ！(大于24即可)；
如下方式设置投影
Begin Projections
Sc:l=0;l=1;l=2
H: l=-0
End Projections】
dis_froz_max=9
dis_num_iter=1000
guiding_centres=true
网络摘录1：
http://emuch.net/html/201311/6573747.html
0.1 下载wannier90-1.2，并编译安装
在 http://www.wannier.org/download.html 下载
wannier90-1.2
修改make.sys
如果装了ifort编译器
LIBS = -lmkl_intel_lp64 -lmkl_sequential -lmkl_core
然后 make all
0.2 重新编译VASP
之后修改vasp
LIB = -L../vasp.5.lib -ldmy \
../vasp.5.lib/linpack_double.o \
../wannier90-1.2/libwannier.a $(SCA) (LAPACK)(LAPACK)(BLAS)
CPP = $(CPP_) -DMPI -DHOST=\"LinuxIFC\" -DIFC \
-DCACHE_SIZE=4000 -DPGF90 -Davoidalloc \
-DMPI_BLOCK=8000 -Duse_collective -DscaLAPACK -DNGZhalf -DVASP2WANNIER90
之后重新make
分五步计算
1. step1：静态计算
INCAR
ISMEAR = 0
SIGMA = 0.05
2.得到虚轨道（导带上方的空轨道）
INCAR
ISMEAR = 0
SIGMA = 0.05
ALGO = Exact
NBANDS = 64
LOPTICS=.TRUE.
NEDOS=2000
3. GW计算
INCAR
ALGO = GW0 ; LSPECTRAL = .TRUE. ; NOMEGA = 50
NBANDS=64
LRPA = .FALSE.
LWANNIER90=.TRUE.
前三步运行vasp 【第三步运行VASP前，需写好wannier90.win参数（可先不写wannier90.win文件，试运行在第三步运行VASP，自动产生符合POSCAR和KPOINTS的wannier90.win文件，再修改该文件，运行一遍VASP），以使得运行VASP后产生 wannier90.amn, wannier90.mmn, wannier90.eig三个文件，供后续第四步的wannier90.x程序调用】
4. step4：运行 wannier90.x wannier90
wannier90.win 【第三步运行完成之后得到的参数文件】
num_wann=8 【参数含义？】
num_bands=8
# for GW uncomment
exclude_bands 9-64
Begin Projections 【开始轨道投影？】
Si:sp3
End Projections
dis_froz_max=9
dis_num_iter=1000
guiding_centres=true
begin unit_cell_cart
2.7150000 2.7150000 0.0000000
0.0000000 2.7150000 2.7150000
2.7150000 0.0000000 2.7150000
end unit_cell_cart
begin atoms_cart
Si 0.0000000 0.0000000 0.0000000
Si 1.3575000 1.3575000 1.3575000
end atoms_cart
mp_grid = 4 4 4
begin kpoints
0.0000000 0.0000000 0.0000000
0.2500000 0.0000000 0.0000000
0.5000000 0.0000000 0.0000000
........
0.2500000 -0.5000000 -0.2500000
-0.2500000 0.2500000 -0.5000000
end kpoints
运行 wannier90.x wannier90
5.---step5:修改wannier90.win参数，再次运行 wannier90.x wannier90；能带即为wannier90_band.dat
num_wann=8
num_bands=8
# for GW uncomment
exclude_bands 9-64
Begin Projections
Si:sp3
End Projections
dis_froz_max=9
dis_num_iter=1000
guiding_centres=true
# Bandstructure plot
restart = plot
bands_plot = true
begin kpoint_path
L 0.50000 0.50000 0.5000 G 0.00000 0.00000 0.0000
G 0.00000 0.00000 0.0000 X 0.50000 0.00000 0.5000
X 0.50000 0.00000 0.5000 K 0.37500 -0.37500 0.0000
K 0.37500 -0.37500 0.0000 G 0.00000 0.00000 0.0000
end kpoint_path
bands_num_points 40
bands_plot_format gnuplot xmgrace
begin unit_cell_cart
2.7150000 2.7150000 0.0000000
0.0000000 2.7150000 2.7150000
2.7150000 0.0000000 2.7150000
end unit_cell_cart
begin atoms_cart
Si 0.0000000 0.0000000 0.0000000
Si 1.3575000 1.3575000 1.3575000
end atoms_cart
mp_grid = 4 4 4
begin kpoints
0.0000000 0.0000000 0.0000000
0.2500000 0.0000000 0.0000000
0.5000000 0.0000000 0.0000000
0.2500000 0.2500000 0.0000000
0.5000000 0.2500000 0.0000000
..........
0.2500000 -0.5000000 -0.2500000
-0.2500000 0.2500000 -0.5000000
end kpoints
再次运行 wannier90.x wannier90；能带即为wannier90_band.dat
如果你有gnuplot，可以直接执行wannier90_band.gnu。
wannier插值能带拟合2相关推荐
- wannier插值能带拟合5
Wannier90是一个构建局域轨道基组的有力工具,一般用来计算紧束缚模型中的占位能和跳跃积分,再结合杂化泛函或GW方法,插值得到更为准确的能带结构,或者计算极化.输运.拓扑不变量等物理性质.在其输入 ...
- wannier插值能带拟合3
今天老王来介绍一下VASP与wannier90结合构建基于第一性原理的紧束缚模型. 首先登录wannier90官网下载wannier90 2.1版本 Download - Wannier90www. ...
- wannier插值能带拟合1
构造Wannier90函数的要点 URL (以下内容最初来自"材料基因论坛"的Brook斑竹,本人只是内容进行了排版,如有不妥可以联系我进行删帖,具体可以参考帖子URL Ⅰ. 构造 ...
- 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数学建模·数据插值与拟合
本篇算是开始数据处理了.首先明确俩概念:插值和拟合.两者最根本的共同点都是基于现有数据进行预测.推演,比如根据现有的天文观测数据预测天体位置.插值问题是数学史上的经典问题,拉格朗日.高斯.牛顿等著名数 ...
最新文章
- 作为程序员，我是怎么通过一个方法记住7000个单词的！
- 在网页中加入百度搜索条
- 40个亿非负整数中找到未出现的数
- 打开电话Android系统调用
- ieee5节点系统图_存储系统的的关键技术之自精简技术
- [另类应用]让SQL Profiler拦不到你的SQL
- 想学PHP来兄弟连是正确的选择 初识兄弟连三周
- win10pe命令打开计算机,手把手教你如何制作win10PE启动盘
- 数电第三章 布尔代数与逻辑函数化简
- 配置库管理及版本管理规范
- 计算机房维修保养记录表,机房设备系统运行及维护记录学习表格.docx
- 通信专业顶刊_通信类期刊排名_2016通信期刊排名_2016通信前沿新技术
- OpenCV之图像锐化
- 在数据迁移过程中如何迁移holo（memory）表？
- VS code + Java 配置与使用
- php如何识别flac文件,flac是什么格式
- APP设计说明书概述
- OPA：open policy agent简介
- Python爬取王一博高清图片，竟然如此简单！
- 性能测试-工具篇：jmeter - Cookie管理和默认请求