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

上一篇:基于CPU加速的VASP安装


下一篇:VASP服务器第一次安装各种软件(中)