Vasp 石墨烯能带计算

这里暂时不做结构优化，使用数据库结构

首先做自洽计算产生正确的基态电子密度

INCAR

SYSTEM = C
PREC = accutat
ISTART = 0
ICHARG = 2
EDIFF = 0.1E-03
NSW = 0
IBRION = 2
ISIF = 2

POSCAR

C
2.45944
1.000 0.000 0.000
0.500 0.866 0.000
0.000 0.000 15.000
2
Cartesian
0.000000 0.000000 0.000000
0.500000 0.288667 0.000000

KPOINTS

Auto
0
Monkhorst Pack
5 5 5
0 0 0

POTCAR文件使用GGA泛函，从赝势文件拷过来。

执行计算脚本

 

第二步，固定电子密度，得到选取给定点的能量

修改INCAR

SYSTEM = C
PREC = accutat
ISTART = 1
ICHARG = 11
EDIFF = 0.1E-03
NSW = 0
IBRION = 2
ISIF = 2

修改KPOINTS

Auto
40
Line-mode
cart
0.0 0.0 0.0
1.000000 0.0 0.0

再次提交任务计算，计算的结果已经写在vasprun.xml文档中

 

后续借助pymatgen工具绘制能带图

导