VASP参数设置详解软件主要功能: 采用周期性边界条件〔或超原胞模型〕处理原子、分子、团簇、纳米线〔或管〕、薄膜、晶体、准晶和无定性材料,以与外表体系和固体l 计算材料的结构参数〔键长、键角、晶格常数、原子位置等〕和构型l 计算材料的状态方程和力学性质〔体弹性模量和弹性常数〕l 计算材料的电子结构〔能级、电荷密度分布、能带、电子态密度和ELF〕l 计算材料的光学性质l 计算材料的磁学性质l 计算材料的晶格动力学性质〔声子谱等〕l 外表体系的模拟〔重构、外表态和STM模拟〕l 从头分子动力学模拟l 计算材料的激发态〔GW准粒子修正〕计算主要的四个参数文件:INCAR ,POSCAR,POTCAR ,KPOINTS,下面简要介绍,详细权威的请参照手册INCAR文件:该文件控制VASP进展何种性质的计算,并设置了计算方法中一些重要的参数,这些参数主要包括以下几类:l对所计算的体系进展注释:SYSTEMl定义如何输入或构造初始的电荷密度和波函数:ISTART,ICHARG,INIWAVl定义电子的优化–平面波切断动能和缀加电荷时的切断值:ENCUT,ENAUG–电子局部优化的方法:ALGO,IALGO,LDIAG–电荷密度混合的方法:IMIX,AMIX,AMIN,BMIX,AMIX_MAG,BMIX_MAG,WC,INIMIX,MIXPRE,MAXMIX–自洽迭代步数和收敛标准:NELM,NELMIN,NELMDL,EDIFFl定义离子或原子的优化–原子位置优化的方法、移动的步长和步数:IBRION,NFREE,POTIM,NSW–分子动力学相关参数:SMASS,TEBEG,TEEND,POMASS,NBLOCK,KBLOCK,PSTRESS–离子弛豫收敛标准:EDIFFGl定义态密度积分的方法和参数–smearing方法和参数:ISMEAR,SIGMA–计算态密度时能量围和点数:EMIN,EMAX,NEDOS–计算分波态密度的参数:RWIGS,LORBITl其它–计算精度控制:PREC–磁性计算:ISPIN,MAGMOM,NUPDOWN–交换关联函数:GGA,VOSKOWN–计算ELF和总的局域势:LELF,LVTOT–结构优化参数:ISIF–等等。
主要参数说明如下:?SYSTEM:该输入文件所要执行的任务的名字取值:字符串,缺省值:SYSTEM?NWRITE:输出容详细程度取值:0~4,缺省值:2如果是做长时间动力学计算的话,最好选0或1(首末步/每步核运动输出),短时运算用2,选3那么会在出错的时候给出说明信息ISTART:决定是否读取WAVECAR文件取值:0~2,缺省0/1 for无/有前次计算的WAVECAR〔波函数〕0:begin 'from scratch',根据INIWAV初始化波函数1:restart with constant energy cut-off,从WAVECAR读取波函数〔重定义平面波集〕2:restart with constant basis set,从WAVECAR读取波函数〔平面波集不变〕?ICHARG:决定如何建立初始电荷密度取值:0~2,缺省值: if ISTART=0 2 else 00:由初始波函数计算电荷密度1:从CHGCAR文件读取电荷密度2:使用原子电荷密度的叠加+10非自洽计算?ISPIN:是否进展spin polarized calculation取值:1,2〔1-no,2-yes〕,缺省值:2?MAGMOM:在ICHARG=2或在CHGCAR中未包含磁化密度〔ICHARG=1〕时,指定每个原子的初始磁化时刻。
取值:实数数组,缺省值:对ISPIN=2 NIONS*1.0,对非共线型磁化体系3*NIONS*1.0?INIWAV如何设置初始波函数,只在ISTART=0时使用取值:0,1〔0-最低动能的平面波,1-随机数〕,缺省值:1IDIPOL控制计算单极、偶极和四极修正取值:1~41~3只计算第一/二/三晶矢方向,适于厚板〔slab〕的计算4所有方向都计算,适于计算孤立分子?PREC进动〔precession〕取值:low/medium/high/normal/ accurate/single),缺省值VASP4.5+采用了优化的accurate来替代high,所以一般不推荐使用high不过high可以确保“绝对收敛〞,作为参考值有时也是必要的同样受推荐的是normal,适于作为日常计算选项受PREC影响的参数有四类:ENCUT;NGX,NGY,NGZ;NGXF,NGYF,NGZF;ROPT如果设置了PREC,这些参数就都不需要出现了,当然直接设置相应的参数也有同样效果具体影响效果见p53~54ENCUT平面波基组的截断能量〔eV〕取值:实数,缺省值:受PREC设置影响,从POTCAR文件中找出相应的ENMAX/ENMIN值来设置。
PREC = Low Medium Accurate NormalENCUT = ENMIN ENMAX ENMAX ENMAX Single HighENMAX ENMAX*1.3对于多个元素不同的ENMAX/ENMIN,都取最大值该参数非常重要,最好不要手工去设置,除非文献告诉你要用多少,或者经过结果可靠性的验证当然,为了测试一下提交的任务,也不妨先设个较小的值NGX,NGY,NGZ:控制FFT网格在三个晶矢方向上的格点数量NGFX,NGFY,NGFZ:控制第二次更准确的FFT网格的格点数量也是两类重要的最好不要去动的参数在未指定的情况下将根据PREC的设置从POTCAR中自动读取PREC=High/Accurate,基组中向量的2倍值,用来防止wrap around errors,得到准确解PREC=Low/Medium/Normal,基组中向量的3/4倍值〔已足够准确到1 meV/atom〕LREAL:决定投射是在实空间还是倒易空间进展取值:.TRUE.〔实空间〕/.FALSE.〔倒易空间〕,缺省值:.FALSE.用于求解赝势的非局域局部用到的一个积分,在倒格空间里采用平面波基组求解,在实空间里那么采用积分球求解。
其他还有两个选项:O or On,A or AutoROPT:在LREAL=Auto or On时,优化控制每个核周围的积分球的格点数取值:实数数组For LREAL=On PREC=Low, 700 points in the real space sphere (ROPT=0.67) PREC=Med, 1000 points in the real space sphere (ROPT=1.0) PREC=High, 1500 points in the real space sphere (ROPT=1.5)For LREAL=Auto PREC=Low, accuracy 10-2 (ROPT=0.01) PREC=Med, accuracy 2*10-3 (ROPT=0.002) PREC=High accuracy 2*10-4 (ROPT=2E-4)?NELM, NELMIN and NELMDL:控制电子自洽循环步数取值:整数NELM:电子自洽循环最大次数缺省值:60NELMIN:电子自洽循环最小次数缺省值:2NELMDL:弛豫次数缺省值:if ISTART=0, INIWAV=1, and IALGO=8,-5,if ISTART=0, INIWAV=1, and IALGO=48,-12,else 0NELMDL可以取负值。
如果初始波函数采用随机赋值,即ISTART=0,INIWAV=1,那么很可能开场的值比拟离谱,那么在第一步核运动循环之前采用NELMDL〔负值〕步的非自洽〔保存初始的H〕步计算将减少计算所需的时间EDIFF:指定电子自洽循环的全局中断条件,用于控制收敛精度取值:实数,缺省值:10-4注意,即使EDIFF=0,NELM步电子自洽循环也会执行EDIFFG:指定离子弛豫循环的中断条件,用于控制核运动的收敛精度取值,实数,缺省值:10*EDIFFEDIFFG>0在两个离子步的总自由能之差小于EDIFFG时停止EDIFFG<0在所有的力都小于EDIFFG时停止EDIFFG=0在NSW步弛豫后停止此参数不支持MD,仅用于弛豫NSW:给出最大离子步数取值:整数,缺省值:0NBLOCK,KBLOCK:取值:整数,缺省值:NBLOCK = 1,KBLOCK = NSW在NBLOCK离子步后对成对相关函数和DOS进展计算,并且把离子配置写入XDATCAR文件在KBLOCK*NBLOCK步主循环后平均的成对相关函数和DOS被写入PCDAT和DOSCAR文件IBRION:决定离子怎样更新和运动取值:-1~3,5~8〔-1-无更新,0-MD,1-RMM-DIIS,2-共轭梯度算法,3-Damped MD,5,6:有限差分,7,8:密度函数扰动理论〕,缺省值:if NSW=0/1,-1,else 0这个参数是和ISIF, IALGO/ALGO一起决定怎么算的最重要的参数。
1~3是三种弛豫的方法,根据ISIF决定是否固定离子位置、晶胞大小和形状,在INCAR中必须设置参数POTIM0是标准的ab-initio MD,不受ISIF影响,即不改变晶胞大小和形状5~8支持Hessian Matrix和phonon frequency的计算以与局部固定的MDPOTIM:IBRION=0时,给出MD每步步长〔fs〕,IBRION=1~3时,给出最小化的度量常量取值:实数,缺省值:IBRION=0无缺省,必须指定,IBRION=1,2,3 0.5?ISIF:决定是否计算应力量以与弛豫中晶胞变化的自由度取值:0~6,缺省值:if IBRION=0(MD) 0 else 2ISIF│calculate│ calculate │relax│ change │ change│ force │stress tensor│ions │cell shape│cell volume----┼-------┼-----------┼------┼---------┼---------0│ yes │ no │yes │ no │ no1│ yes │ trace only │yes │ no │ no2 │ yes │ yes │yes │ no │ no3 │ yes │ yes │yes │ yes │ yes4 │ yes │ yes │yes │ yes │ no5 │ yes │ yes │ no │ yes │ no6 │ yes │ yes │ no │ yes 。