VSe2CDW态能带计算

第一步&＃xff0c;得到单胞的POSCAR和POTCAR

第二部&＃xff0c;推导CDW超胞与单胞的晶格矢量关系&＃xff0c;写TRANSMAT.in文件&＃xff0c;如下

Read transformation matrix from the TRANSMAT.in file if it exists. 第一行是注释行
2 0 0 # must be three integers
0 3 0 # must be three integers
0 0 3 # must be three integers

使用vaspkit 401生成超包的POSCAR

并写一个标准的KPOINTS

Kpoints from vasprun.xml
0
Gamma
13 13 1


第三步&＃xff1a;开始进行结构优化&＃xff0c;ISIF&＃61;2, 固定晶格矢量&＃xff0c;INCAR如下

SYSTEM&＃61;2D_VSe2
ISTART&＃61;0
NWRITE&＃61;2
ENCUT&＃61;500
GGA&＃61;PE
NSW&＃61;200
ISIF&＃61;2
ISYM&＃61;2
NBLOCK&＃61;1
KBLOCK&＃61;1
IBRION&＃61;2
NELM &＃61; 80
EDIFF &＃61; 1E-05
EDIFFG &＃61; -0.01
ALGO &＃61; Normal
LDIAG&＃61;.True.
LREAL&＃61;.FALSE.
ISMEAR&＃61;0
SIGMA&＃61;0.02
ICHARG&＃61;2
LPLANE&＃61;.TRUE.
NPAR&＃61;4
LSCALU &＃61; .FALSE.
NSIM&＃61;4
LWAVE&＃61;.FALSE.
LCHARG&＃61;.FALSE.
ICORELEVEL&＃61;1


从输出文件查看计算是否达到精度要求

grep "reached required accuracy" output.419046


查看体系总能量

grep TOTEN OUTCAR | tail -n 1


比较优化前和优化后晶体结构差别

vim -O CONTCAR POSCAR

方法一&＃xff1a;由于4x4均匀体系&＃xff0c;布里渊区只有缩放&＃xff0c;手动缩小倍数布里渊区计算&＃xff0c;没有能带折叠

所以进行静态计算

SYSTEM &＃61; 2D_VSe2
ISTART &＃61; 0
NWRITE&＃61;2
PREC&＃61;Accurate
ENCUR&＃61;500
GGA&＃61;PE
NSW&＃61;0
ISIF&＃61;2
ISYM&＃61;2
NBLOCK&＃61;1
KBLOCK&＃61;1
IBRION&＃61;-1
NELM&＃61;80
EDIFF&＃61;1E-05
EDIFFG&＃61;-0.01
ALGO&＃61; Normal
LDIAG&＃61;.TRUE.
LREAL&＃61;.FALSE.
ISMEAR&＃61;0
SIGMA&＃61;0.02
ICHARG&＃61;2
LPLANE&＃61;.TRUE.
NPAR &＃61; 4
LSCALU &＃61; .FALSE.
NSIM &＃61; 4
LWAVE &＃61; .TRUE.
LCHARG &＃61; .TRUE.
ICORELEVEL&＃61;1


KPOINTS用标准

Kpoints from vasprun.xml
0
Gamma
16 16 1


得到自洽的电荷密度和波函数文件&＃xff0c;再进行一次能带计算

SYSTEM &＃61; 2D_VSe2
ISTART &＃61; 1
NWRITE&＃61;2
PREC&＃61;Accurate
ENCUR&＃61;500
GGA&＃61;PE
NSW&＃61;0
ISIF&＃61;2
ISYM&＃61;2
NBLOCK&＃61;1
KBLOCK&＃61;1
IBRION&＃61;-1
NELM&＃61;80
EDIFF&＃61;1E-05
EDIFFG&＃61;-0.01
ALGO&＃61; Normal
LDIAG&＃61;.TRUE.
LREAL&＃61;.FALSE.
ISMEAR&＃61;0
SIGMA&＃61;0.02
ICHARG&＃61;11
LPLANE&＃61;.TRUE.
NPAR &＃61; 4
LSCALU &＃61; .FALSE.
NSIM &＃61; 4
LWAVE &＃61; .FALSE.
LCHARG &＃61; .FALSE.
ICORELEVEL&＃61;1
LORBIT&＃61;11


布里渊区用缩小4倍的

K-Path Generated by VASPKIT.20
Line-Mode
Reciprocal0.0000000000 0.0000000000 0.0000000000 GAMMA 0.1250000000 0.0000000000 0.0000000000 M 0.1250000000 0.0000000000 0.0000000000 M 0.0833333333 0.0833333333 0.0000000000 K 0.0833333333 0.0833333333 0.0000000000 K 0.0000000000 0.0000000000 0.0000000000 GAMMA

画图显示

方法二

结构优化后&＃xff0c;进行一次能带反折叠计算

INCAR

SYSTEM &＃61; VSe2
ISTART &＃61; 1
LREAL &＃61; F
PREC &＃61; N
LWAVE &＃61; .TRUE.
LCHARG &＃61; F
ISMEAR &＃61; 0
SIGMA &＃61; 0.05
NELM &＃61; 60
NELMIN &＃61; 6
EDIFF &＃61; 1E-08


由于之前结构优化已经产生了超包&＃xff0c;这里需要准备以下几个文件

用于产生超胞的变化矩阵TRANSMAT.in文件&＃xff08;可选文件&＃xff09;。我们先看看TRANSMAT.in的格式

Read transformation matrix from the TRANSMAT.in file if it exists. 第一行是注释行
4 0 0   # must be three integers
0 4 0   # must be three integers
0 0 0   # must be three integers

原始晶包的KPOINTS文件

KPATH for MoS2
20
Line mode
Rec
0.0000000 0.5000000 0.0000000 # M
0.3333333 0.3333333 0.0000000 # K0.3333333 0.3333333 0.0000000 # K
0.0000000 0.0000000 0.0000000 # GAMMA

第三步&＃xff1a;cp SUPERCELL.vasp POSCAR和cp TRANSMAT TRANSMAT.in&＃xff08;当TRANSMAT.in文件不存在时&＃xff09;, 并运行vaspkit-281生成KPOINTS文件&＃xff1b;

这里生成的KPOINTS就是新晶胞的KPOINTS文件

第四步&＃xff1a;提交vasp作业&＃xff1b;

第五步&＃xff1a;运行vaspkit-282&＃xff1b;

第六步&＃xff1a;最后&＃xff0c;可以利用vaspkit/examples/band_unfolding/eps_plot.py脚本画图