VASP-中INCAR设置

INCAR（参数的首字母代表的参数性质：I初始化L逻辑开关E能量N数目T温度）一般要设置的关键词为：SYSTEM,ENCUT,ISTART,ICHARG,PREC,ISMEAR,SIGMA。缺省值就是默认值。是指一个属性、参数在被修改前的初始值。SYSTEM任务的名字ENCUT平面波的切断能（energycutoffineV）:defaulttakenfromPOTCAR-fileimportant!重要到几乎最好不要手工去设置除非文献告诉你要用多少，或者经过结果可靠性的验证当然，为了测试一下提交的任务，也不妨先设个较小的值附加说明：当且仅当POTCAR里头没有设置ENCUT时(其实貌似没有才是常态)，才受PREC设置影响从POTCAR里找出相应的ENMAX/ENMIN值来设置。PREC=LowMediumAccurateHighENCUT=ENMINENMAXENMAX130%ENMAX对于多个元素的POTCAR不同的ENMAX/ENMIN，都取最大值ISTART决定是否读取WAVECAR文件。取值：0~2，缺省0/1for无/有前次计算的WAVECAR（波函数）0：begin'fromscratch'，根据INIWAV初始化波函数1：restartwithconstantenergycut-off，从WAVECAR读取波函数（重定义平面波集）2：restartwithconstantbasisset，从WAVECAR读取波函数（平面波集不变）在进行能带结构、电子态密度等性质计算时设置为1，其他情况一般都设置为0。ICHARG决定如何建立初始电荷密度。取值：0~2，缺省值:ifISTART=02else00：由初始波函数计算电荷密度1：从CHGCAR文件读取电荷密度2：使用原子电荷密度的叠加+10非自洽计算PRECprecession:medium,highorlow(VASP.4.5+also:normal,accurate)Default:MediumVASP4.5+采用了优化的accurate来替代high,所以一般不推荐使用high。不过high可以确保'绝对收敛'，作为参考值有时也是必要的。同样受推荐的是normal，作为日常计算选项，可惜的是说明文档提供的信息不足。受PREC影响的参数有四类：ENCUT;NGX,NGY,NGZ;NGXF,NGYF,NGZF;ROPT;如果设置了PREC，这些参数就都不需要出现了当然直接设置相应的参数也是同样效果的，这里不展开了，随后详释ISMEAR决定每个波函数的部分占位fnk如何设置。取值：-5|-4|-3|-2|0|N-5-带有Blochl修正的四面体方法,-4-不带Blochl修正的四面体方法,-3-根据INCAR文件中提供的smearing参数执行循环，-2-从WAVECAR文件中读取,-1-Fermi-smearing,0-Gaussiansmearing,0-methodofMethfessel-PaxtonorderN，缺省值：1采用部分占位波函数，用一个函数来平滑积分，尤其是对于金属体系可减少k点。采用所谓部分占有波函数，用一个函数来平滑积分，尤其是对于金属体系可减少k点SIGMA决定smearing的宽度（eV）。取值：实数，缺省值：0.2其他参数：NWRITE输出内容详细程度。取值：0~4，缺省值：2如果是做长时间动力学计算的话，最好选0或1(首末步/每步核运动输出)，短时运算用2，选3则会在出错的时候给出说明信息。ISPIN是否进行spinpolarizedcalculation。取值：1，2（1-no，2-yes），缺省值：2MAGMOM在ICHARG=2或在CHGCAR中未包含磁化密度（ICHARG=1）时，指定每个原子的初始磁化时刻。取值：实数数组，缺省值:对ISPIN=2NIONS*1.0，对非共线型磁化体系3*NIONS*1.0INIWAV如何设置初始波函数，只在ISTART=0时使用。取值：0，1（0-最低动能的平面波，1-随机数），缺省值：1。IDIPOL控制计算单极、偶极和四极修正。取值：1~4。1~3只计算第一/二/三晶矢方向，适于厚板（slab）的计算4所有方向都计算，适于计算孤立分子LREAL决定投射是在实空间还是倒易空间进行。取值：.TRUE.（实空间）/.FALSE.（倒易空间），缺省值：.FALSE.用于求解赝势的非局域部分用到的一个积分，在倒格空间里采用平面波基组求解，在实空间里则采用积分球求解。其他还有两个选项：OorOn，AorAuto。On和.TRUE.的差别在于是否使用King-Smith算法优化，设为Auto则进行自动选择，推荐使用。ROPT在LREAL=AutoorOn时，优化控制每个核周围的积分球内的格点数。取值：实数数组ForLREAL=OnPREC=Low,700pointsintherealspacesphere(ROPT=0.67)PREC=Med,1000pointsintherealspacesphere(ROPT=1.0)PREC=High,1500pointsintherealspacesphere(ROPT=1.5)ForLREAL=AutoPREC=Low,accuracy10-2(ROPT=0.01)PREC=Med,accuracy2*10-3(ROPT=0.002)PREC=Highaccuracy2*10-4(ROPT=2E-4)EDIFF指定电子自洽循环的全局中断条件，用于控制收敛精度。取值：实数，缺省值：10-4注意，即使EDIFF=0，NELM步电子自洽循环也会执行。是电子结构部分自洽迭代循环是，判断是否自洽了的条件，上次和当前两次迭代中总能和本征值的变化都小于ＥＤＩＦＦ，则电子构件部分迭代循环停止。如果ＥＤＩＦＦ＝０，则进行ＮＥＬＭ步迭代后停止迭代。EDIFFG指定离子弛豫循环的中断条件，用于控制核运动的收敛精度。取值，实数，缺省值：10*EDIFFEDIFFG0在两个离子步的总自由能之差小于EDIFFG时停止EDIFFG0在所有的力都小于EDIFFG时停止。EDIFFG=0在NSW步弛豫后停止此参数不支持MD，仅用于弛豫。ISIF=3calculatestressandwhattorelaxNSW指定核运动步数（最大离子步数）缺省：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-DampedMD，5,6：有限差分，7,8：密度函数扰动理论），缺省值：ifNSW=0/1，-1，else0这个参数是和ISIF,IALGO/ALGO一起决定怎么算的最重要的参数。~3是三种弛豫的方法，根据ISIF决定是否固定离子位置、晶胞大小和形状，在INCAR中必须设置参数POTIM。0是标准的ab-initioMD，不受ISIF影响，即不改变晶胞大小和形状。5~8支持HessianMatrix和phononfrequency的计算以及部分固定的MD。NGX；NGY；NGZ控制FFT网格在三个晶矢方向上的格点数量。NGFX，NGFY，NGFZ控制第二次更精确的FFT网格的格点数量。也是两类重要的最好不要去动的参数。在未指定的情况下将根据PREC的设置从POTCAR中自动读取。PREC=High/Accurate，基组中向量的2倍值，用来避免wraparounderrors，得到精确解。PREC=Low/Medium/Normal，基组中向量的3/4倍值（已足够精确到1meV/atom）。POTIMBRION=0时，给出MD每步步长（fs），IBRION=1~3时，给出最小化的度量常量。取值：实数，缺省值：IBRION=0无缺省，必须指定，IBRION=1,2,30.5RWIGSWigner-Seitz半径DOS计算用。用在计算分波态密度以及每根能带对应的波函数按ｓｐｄ和位置投影NPARNPAR并行计算band的节点数，每一个节点计算一个band当然可以提高并行效率，减少通讯量，不过貌似现在硬件的主要限制还是内存，而这个选项的使用可能会大幅增加内存的需求ISYM=0决定是否计算应力张量以及弛豫中晶胞变化的自由度。取值：0~6，缺省值：ifIBRION=0(MD)0else2；0-nonsym1-usesym是否使用对称性ISIF│calculate│calculate│relax│change│change│force│stresstensor│ions│cellshape│cellvolume──┼─────┼───────┼───┼──────┼──────0│yes│no│yes│no│no1│yes│traceonly│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│yes7│yes│yes│no│no│yes