VASP输入INCAR文件

计算的体系名称
SYSTEM = CaCO3

初始I/O设置(读入、读出)

ISYM = 0 #0代表不使用对称性，推荐表面计算使用，晶体计算默认使用2，杂化泛函默认使用3
ISTART = 0 # 0随机生成波函数，1从WAVECAR读取波函数，失败随机生成(一般续算会选择1)，初始计算推荐参数0
ICHARG = 2 #2初猜，原子电荷叠加生成电荷密度，初始计算推荐
LWAVE = .FALSE. #是否输出波函数，如果硬盘内存空间小，使用 .FALSE. 不输出波函数
LCHARG = .FALSE. #是否输出电子密度文件，如果内存不够，不建议输出
LVTOT = .TRUE. #loacl potential文件，局部的总电位，TRUE写入LOCPOT文件
LELF = .TRUE. #输出电子局域化函数，决定是否创建 ELFCAR文件，此时，应先设置ISPIN=2,计算精度PREC=High
LORBIT = 11 #输出材料的分波态密度，10输出每个原子角动量即能级s p d f…参数11输出每个原子的磁角动量s px py pz…
#以上四个文件一般在需要计算准确电子结构（自洽运算）时才会出现

Electronic Relaxation(电子步)

ENCUT = 520平面阶段能，由赝势决定（POTCAR中的ENMAX*1.0~1.3）（数值越大，精度越高，计算的基组多点）
NELM = 100 #最大电子步数量，默认60步（难收敛体系，如过渡金属，可设置为200、300~）
ALGO = Fast #自洽循环算法（F为N和V结合,N为DAV收敛性好速度慢；V 是RMM收敛差，速度快)
PREC = Accurate #精度设置，A高等精度，M中等，N普通精度
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 = 1000 #最大离子步的个数
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文件中存在上述字符串，则表明计算成功收敛；若没有，则未收敛