一、VASP
全称Vienna Ab-initio Simulation Package,是维也纳大学Hafner小组开发的进行电子结构计算和量子力学-分子动力学模拟软件包。
它是目前材料模拟和计算物质科学研究中最流行的商用软件之一。
二、输入文件
VASP一般来说最简单的输入文件有:INCAR,POSCAR,POTCAR,KPOINTS;如果是在服务器上计算的,那么还需要一个提交任务的脚本。
三、INCAR
SYSTEM--name of System 任务的名字
ISTART--startjob: restart选项 0-3
缺省0/1 for 无/有 前次计算的WAVECAR(波函数)
1)'restart with constant energy cut-off'
2)'restart with constant basis set'
3)'full restart including wave function and charge prediction'
一般不用写这个参数,除非你要继续算之前中断的任务。
ICHARG--charge:
1-file
2-atom
10-const
Default:if ISTART=0 2 else 0
如果想省空间,设置为2,想省时间就设置为1;2是读原子的数据,1是读之前算过的数据(如果要继续算的话这个比较好用)。
ISPIN--spin polarized calculation (2-yes 1-no),default 1,默认为1,选择是否自旋极化,1是不要自旋极化,2是选择自旋极化;
PREC--precession: medium, high or low(VASP.4.5+ also: normal, accurate)计算精度的选择:有3个选择medium,high or low;
ENCUT--energy cutoff in eV : default taken from POTCAR-file,截断能的选择,如400eV;
NELM--最大电子自洽循环次数 一般设置为NELM = 60;
EDIFF--电子SC循环的收敛精度,缺省:1e-4 注意,即使EDIFF=0,NELM步也会执行,一般为EDIFF=1e-4;
EDIFFG--核运动的收敛精度,缺省:EDIFF*10 (总能量),EDIFFG<0 则在所有的力都小于EDIFFG时停止;
EDIFFG=0 则在NSW步后停止 此参数不支持MD,仅用于Relax,一般EDIFFG=10*EDIFF就是EDIFF的10倍;
NSW--指定核运动步数 缺省: 0,一般NSW=60,就是让原子跑60步;
IBRION--ionic relaxation: -1-Fixed
0-MD
1-quasi-New
2-CG
3-Damp
5-freq
Default if NSW=0 or 1 IBRION=-1 else IBRION=0 这个参数是和ISIF;IALGO/ALGO一起决定怎么算的最重要的参数;
1-3 是三种Relax的方法,受ISIF决定是否固定核位置、晶胞大小和形状 0 是标准的ab-initio MD,不受ISIF影响,即不改变晶胞大小和形状;
一般做优化IBRION=2;做静态计算(原子坐标不动,只优化离子)IBRION=-1。
四、输出文件
OUTCAR:主要的输出文件,里面有计算的参数和每一个循环步骤的输出结果,含有每一步骤的能量结果,位置,力收敛标准,能量收敛,电荷,磁矩等信息;
CHG和CHGCAR:其中CHGCAR中的数据可以用来处理电荷密度图,电荷差分密度图ELF等信息;
OSICAR:为计算进程数据,通过这个文件可以看运行情况和能量收敛情况;通过 tail -100 OSICAR,可以查询最后100行的OSICAR文件,可以实时看运行进度。