VASP输入文件

1.INCAR文件

# 计算的体系的名称
SYSTEM = Hybrid

# I/O设置(读入、读出)
ISTART = 0  # 0代表一个全新的计算
ICHARG = 2
LWAVE = .FALSE.  
LCHARG = .TRUE.  #输出电子密度文件
LVOT = .TRUE.    #loacl potential文件
LELF = .TRUE.  #输出电子局域化函数
LORBIT = 11    #输出材料的分波态密度
#以上四个文件一般在需要计算准确电子结构（自洽运算）时才会出现

# Electronic Relaxation(电子步)
ENCUT = 600    #平面阶段能，由赝势决定（POTCAR中的ENMAX*1.25~1.50）（数值越大，精度越高）
NELM = 100    #最大电子步数量，默认60步（难收敛体系，如过渡金属，可设置为200、300~）
ALGO = Fast    #自洽循环算法（Fast,Normal,VaryFast)
PREC = Accurate    #精度设置
ISMEAR = 0    # 金属0+，非金属0-（金属不可以设置为负数）
SIGMA = 0.05  # 划分布里渊区,展宽
EDIFF = 1E-5    #电子步收敛精度
AMIX = 0.1    #
BMIX = 0.01   #加快收敛
NEDOT = 2000
EMIN = -10.0  #两个
EMAX = 10.0  #细化能级

# Ionic Relaxation (离子步)
IBRION = 2     #离子弛豫的算法（2较稳定）
NSW = 300      #最大离子步的个数
EDIFFG = -1E-2   #离子步收敛标准（默认是EDIFF*10）
ISIF = 3    # 对晶胞的弛豫方法(3:全弛豫 2:固定体积的弛豫 4:固定体积但允许形状改变)

# Polarization (含磁性体系)
ISPIN = 2       # 自旋极化   
LOSRBIT = .False.  #考虑非线性磁性（默认不考虑）

# Parallization(并行计算的相关设定)
LREAL = Auto   #计算进行在实空间或是导空间
NPAR = 1    #取节点数量或核心数量的根号值
LPLANE = .FALSE. #针对平面波的设置
  • ‘#’ 表示注释符，后面的内容不会被执行

2.KPOINTS文件

# 划分布里渊区网格
Automatic mesh
0
Gamma        #划分网格
8    8    8
0.0  0.0  0.0
# 另一种KPOINTS文件形式
Automatic mesh
0
Auto
25.0

 

3.POSCAR文件

# 通过VESTA建模得到
# 包括体系名称，晶胞基矢信息，原子元素种类、数量和具体坐标

4.POTCAR文件

# 从赝势库中得到的赝势信息，一般不需要进行操作

 

VASP程序的执行

mpirun -np x vasp_std < INCAR
# x是指CPU的线程数（核心数）

 

VASP输出文件

1.OSZICAR

#第一行内容
F   #体系的自由能
E0  #最后体系的总能量
dE  #体系能量的改变量

# 使用linux命令得到文件中我们需要的数据
grep E0 OSZICAR  #打印出OSZICAR文件中所有含有'E0'的行

2.OUTCAR文件

# 包含VASP全部输入信息的总结(日志文件)
# 最后一行：计算成本，磁性信息，电荷信息
grep 'reached required accurary' OUTCAR
# 若OUTCAR文件中存在上述字符串，则表明计算成功收敛；若没有，则未收敛

3.DOSCAR文件

# 材料的态密度信息（材料在不同能级上各自的占据情况/概率）（5列、3列）
# 若在电子步设置了LORBIT=11，则会输出分波态密度（19列）更大

4.CHGCAR文件

#电子在空间位置上的分布，可直接通过VESTA文件画图

5.LOCALCAR文件

#电子在空间中的能级分布

6.ELFCAR文件

#定义的电子局域化函数

7.EIGENVAL文件

#电子在导空间上的分布状态