本文最后更新于 3 分钟前,文中所描述的信息可能已发生改变。
Example
text
num_bands = # 设置能带数量
exclude_bands = # 设置排除的能带数量
num_wann = # 设置用于计算的Wannier函数数量
dis_win_min = # 设置解耦合过程的能量窗口的最小值
dis_win_max = # 设置解耦合过程的能量窗口的最大值
dis_froz_min = # 设置冻结窗口的最小值
dis_froz_max = # 设置冻结窗口的最大值
num_iter = 200 # 设置计算的迭代次数
num_print_cycles = 20 # 设置打印输出的周期
# 通过fatband分析判断主要产生贡献的轨道
# projection
begin projections
Co: d # Co的d轨道作为Wannier函数的初始猜测
Sn: p # Sn的p轨道作为Wannier函数的初始猜测
end projections
# 设置是否绘制费米弧和费米面,生成wannier90_hr.dot
#restart = plot
#fermi_energy =
#fermi_surface_plot = true
#write_hr = true
mp_grid = 4 4 4 # 设置k点网格
# band structure plot
bands_plot = true # 设置是否绘制能带结构图
# kpoint_path
begin kpoint_path # 定义能带结构计算的k点路径
G 0.0 0.0 0.0 M 0.5 0.0 0.0
M 0.5 0.0 0.0 K 0.3333333333 0.3333333333 0.0
K 0.3333333333 0.3333333333 0.0 G 0.0 0.0 0.0
G 0.0 0.0 0.0 A 0.0 0.0 0.5
A 0.0 0.0 0.5 L 0.5 0.0 0.5
L 0.5 0.0 0.5 H 0.3333333333 0.3333333333 0.5
H 0.3333333333 0.3333333333 0.5 A 0.0 0.0 0.5
L 0.5 0.0 0.5 M 0.5 0.0 0.0
H 0.3333333333 0.3333333333 0.5 K 0.3333333333 0.3333333333 0.0
end kpoint_path
num_bands
解纠缠需要 num_bands > num_wann,即额外预留一部分能带。预留 4 条导带或 2 条价带+2 条导带均可,也可预留更多能带。
根据第一性计算中的能带数进行设置。
num_wann
根据 projection 中选择的轨道进行计算,s 贡献 1 条轨道,p 贡献 3 条轨道,d 贡献 5 条轨道。
根据原子数,轨道数,自旋状态进行计算。
dis_win_min | dis_win_max | dis_froz_min | dis_froz_max
projections
选择投影轨道之前,做一次能带成分分析,也就是画 Fatband 图和 PDOS 图。确认费米面附件的轨道成分,把这些轨道作为投影轨道。
kpoint_path
从 KPOINTS(KPATH.in) 文件中复制路径。
Template
dis_win_min =
dis_win_max =
dis_froz_min =
dis_froz_max =
num_wann =
num_bands =
exclude_bands =
num_iter = 200
num_print_cycles = 20
begin projections
···
end projections
begin unit_cell_cart
... ... ...
... ... ...
... ... ...
end unit_cell_cart
begin atoms_cart
... ... ...
... ... ...
... ... ...
... ... ...
end atoms_cart
bands_plot = true
#restart = plot
#fermi_energy =
#fermi_surface_plot = true
#write_hr = true
mp_grid = .. .. ..
begin kpoints
... ... ...
... ... ...
... ... ...
... ... ...
end kpoints