VASP+wannier90 构建紧束缚模型(补充)
欢迎新来的朋友们,
昨天发文以后微信群里收到很多好的反馈,老王在这里谢谢大家的建议。
昨天的文章:
一文搞定VASP+wannier90构建紧束缚模型
首先,老王更正一个错误。
patch -p0 < mlwf.patch
上面这个命令是在vasp编译目录下运行的,而不是在src目录下运行。
另外,老王对昨天的文章做两个补充。
第一个对于初学者的。
wannier90.win 文件是用于构造wannier函数的初始输入文件。
在这个文件最开始只要设置好num_bands, num_wann,frozen窗口 以及投影轨道部分就好。老王习惯把这些都写上。
num_bands = 48 #num_bands = NBANDS of INCAR >= num_wann
num_wann = 23 # number of wannier bands
dis_win_min = -20.00d0
dis_win_max = 30d0
dis_froz_min = -11.50d0 # energy window
dis_froz_max = 4.0d0
#dis_num_iter = 1000
num_iter = 200
num_print_cycles = 20
#spinors = .true.
# projections
begin projections
Mo : d
Si : p
N : p
end projections
#Bandstructure plot
#restart = plot
bands_plot = true
#fermi_energy = 7.9099
#fermi_surface_plot = true
begin kpoint_path
K -0.3333 0.6666 0.0 G 0.0 0.0 0.0
G 0.0 0.0 0.0 M 0.0 0.5 0.0
M 0.0 0.5 0.0 K -0.3333 0.6666 0.0
end kpoint_path
bands_num_points 101
bands_plot_format gnuplot xmgrace
运行完VASP以后这个文件会自动补全。包括结构信息,k点信息。
然后检查是否生成wannier90.amm,wannier90.mmn 以及wannier90.eig文件。
有了这三个文件就可以进一步调节参数,构造wannier函数了。
运行命令为:
wannier90.x wannier90 &
这里后面一个winner90代表wannier90.win的文件.win前面的部分。如果你的win文件是name.win运行的地方就可以是
wannier90.x name &
运行过程中当num_bands > num_wann 有个迭代过程,控制迭代步数可以通过
dis_num_iter = 200 # 迭代1,默认200步
num_iter = 200 # 迭代2,默认200步
当num_bands = num_wann时只有后面一个迭代过程。
最后会生成wannier90.chk文件,restart计算会读取这个文件。
加入能带部分的话第一次运行就会出来,wannier90_band.dat文件,这个可以用来和DFT结果对比。
#Bandstructure plot
#restart = plot
bands_plot = true
#fermi_energy = 7.9099
#fermi_surface_plot = true
begin kpoint_path
K -0.3333 0.6666 0.0 G 0.0 0.0 0.0
G 0.0 0.0 0.0 M 0.0 0.5 0.0
M 0.0 0.5 0.0 K -0.3333 0.6666 0.0
end kpoint_path
bands_num_points 101
bands_plot_format gnuplot xmgrace
好消息:
昨天一直没找到的一篇拓扑绝缘体Bi2Se3练习教程,今天找到了。
https://zhuanlan.zhihu.com/p/25432728
第二个问题如下:
可以结合一些设置非原子点位轨道初始位置猜测的案例,解释下轨道初始位置猜测的重要性和技巧
初始猜测设置对了 基本上就稳了 必然收敛好结果 物理图像还清晰得狠
但是老王没有接触过相关案例。有懂这方面的可以和老王说一说。
重要性不言而喻,对于这样的体系,能够猜测出费原子点位的初始位置,收敛会很快,物理图像也很清晰。这个一般比较难,考验物理功底,是我们追求的方向。做科研要习惯于谋定而后动,想清楚了一般都比较顺利,节省时间。
老王要强调的是初始轨道设置问题。这部分要结合一些化学轨道方面的知识。例如,石墨烯,sp2杂化,金刚石sp3杂化。投影轨道用杂化轨道代替。
评论
我爱学习爱我

干货推荐

提取文献数据点,你只知道GetData?
木若林溪
558 1

模态分析到底在算啥?搞懂这个,结构为什么会自己振动就明白了
莫名其妙
64 4

lammps中文教程12
春风得意马蹄疾
87 2

计算机材料设计Materials-Studio教程5
活着
33 2

小胶质细胞
Mr弘🔬
63 8

这些真的是棕色脂肪吗?
Mr弘🔬
164 2
