计算石墨烯能带_六，能带计算实例

跌停

4320人浏览 · 2021-01-13 02:15:39

跌停 · 2021-01-13 02:15:39 发布

我们以第四节中的石墨烯C为例，来进行一个简单的能带计算。分为三个步骤：结构优化，自洽计算，非自洽计算。

1，结构优化

INCAR文件：

SYSTEM = Graphene
ISTART = 0 
ICHARG = 2  
PREC   = Accurate
ENCUT  = 600
EDIFF  = 0.1E-07
NSW    = 100
IBRION = 2
ISIF   = 2  
ISMEAR = -5

POSCAR文件：

Graphene
1.0
 2.4600000000    0.0000000000    0.0000000000
-1.2300000000    2.1300000000    0.0000000000
 0.0000000000    0.0000000000    10.0000000000
C
2
Direct
0.666666667    0.666666667    0.000000000
0.000000000    0.000000000    0.000000000

POTCAR文件自行从赝势库里载入。

KPOINTS文件：

Automatic generation
0
Gamma
5 5 5
0 0 0

把上述文件放在一个目录下运行vasp，得到一系列的文件，把产生的CONTCAR命名为POSCAR供下一步计算。

2，自洽计算

INCAR文件：

SYSTEM = Graphene
ISTART = 0 
ICHARG = 2  
PREC   = Accurate
ENCUT  = 600
EDIFF  = 0.1E-07
#NSW    = 100
#IBRION = 2
#ISIF   = 2  
ISMEAR = -5

POSCAR文件：（上一步的CONTCAR文件）

Graphene
   1.00000000000000
     2.4684164524000001    0.0000000000000000    0.0000000000000000
    -1.2342082262000000    2.1377113549000000    0.0000000000000000
     0.0000000000000000    0.0000000000000000    9.9990577698000003
   C
     2
Direct
  0.6666669769999984  0.3333329880000022  0.0000000000000000
  0.0000000000000000  0.0000000000000000  0.0000000000000000

  0.00000000E+00  0.00000000E+00  0.00000000E+00
  0.00000000E+00  0.00000000E+00  0.00000000E+00

POTCAR文件同上。

KPOINTS文件：

Automatic generation
0
Gamma
5 5 5
0 0 0

把上述文件放在一个目录下运行vasp，保存产生的文件。

3，非自洽运算（能带）

将第二步产生的CHG，CHGCAR文件保存复制到新的目录。

INCAR文件：

SYSTEM = Graphene
ISTART = 1 
ICHARG = 11  
PREC   = Accurate
ENCUT  = 600
EDIFF  = 0.1E-07
#NSW    = 100
#IBRION = 2
#ISIF   = 2  
ISMEAR = 1
LORBIT = 11

POSCAR文件：

Graphene
   1.00000000000000
     2.4684164524000001    0.0000000000000000    0.0000000000000000
    -1.2342082262000000    2.1377113549000000    0.0000000000000000
     0.0000000000000000    0.0000000000000000    9.9990577698000003
   C
     2
Direct
  0.6666669769999984  0.3333329880000022  0.0000000000000000
  0.0000000000000000  0.0000000000000000  0.0000000000000000

  0.00000000E+00  0.00000000E+00  0.00000000E+00
  0.00000000E+00  0.00000000E+00  0.00000000E+00

POTCAR文件同上。

KPOINTS文件：（此时k点要手动选取高的对称点）

k-points along high symmetry lines
11
Line-mode
Reciprocal
0.00 0.00 0.00 
0.50 0.00 0.00

0.50 0.00 0.00
0.33 0.33 0.00

0.33 0.33 0.00
0.00 0.00 0.00

运行vasp，得到vasprun.xml文件，用p4vasp打开该文件即可查看能带。

第一性计算的能带(与下面kpoint路径不一致)

TB模型计算的能带

POSCAR文件：网上有一些现成的POSCAR库，下载cif文件，然后以文本形式打开即可相关参数。

DAMO开发者矩阵

DAMO开发者矩阵，由阿里巴巴达摩院和中国互联网协会联合发起，致力于探讨最前沿的技术趋势与应用成果，搭建高质量的交流与分享平台，推动技术创新与产业应用链接，围绕“人工智能与新型计算”构建开放共享的开发者生态。

更多推荐

光感日常，电懂人心：消费电子与生活的光电融合

DAMO开发者矩阵

山东工厂西门子6RA直流调速器频繁跳闸？6RA24/6RA70/6RA80五大报警代码故障定位维修指南

维修品类覆盖大功率进口变频器、高压变频器、直流调速器、大功率伺服驱动器、伺服电机、进口触摸屏、功率单元、软启动器、工业机器人，并可承接空压机变频器控制系统维修。跳闸报警时，不同系列的代码体系并不完全相同：6RA24与6RA70采用一致的F0xx代码，6RA80改为F60xxx代码，6RA23/6RA28/6RA27则是更早的代码体系。缺相 F004（电枢缺相） F60004（电枢缺相）、F6000