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文件中存在上述字符串,则表明计算成功收敛;若没有,则未收敛
标签:体系,文件,FALSE,程序执行,弛豫,VASP,收敛 来源: https://www.cnblogs.com/fengxiaoqiang/p/14089143.html