VASP结构优化计算输入文件，程序执行和输出文件

VASP输入文件

1.INCAR文件

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

# I/O设置(读入、读出)
ISTART = 0  # 0代表一个全新的计算
ICHARG = 2
LWAVE = .FALSE.  
LCHARG = .FALSE.
LVOT = .FALSE.
LELF = .FALSE.  #以上四个文件一般在需要计算准确电子结构（自洽运算）时才会出现

# 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   #加快收敛

# 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文件中存在上述字符串，则表明计算成功收敛；若没有，则未收敛