成果关键词详解.pptxVIP

VASP的关键词详解侯柱锋复旦大学物理系2005,12,3北京，宏剑

@HouZhufeng2内容描述体系：SYSTEM设置如何输入或构造初始的电荷密度和波函数：ISTART,ICHARG,INIWAV设置电子的优化：平面波切断动能和缀加电荷时的切断值：ENCUT,ENAUG电子部分优化的方法：ALGO,IALGO,LDIAG电荷密度混合的方法：IMIX,AMIX,AMIN,BMIX,AMIX_MAG,BMIX_MAG,WC,INIMIX,MIXPRE,MAXMIX自洽迭代步数和收敛标准：NELM,NELMIN,NELMDL,EDIFF设置原子的驰豫：原子如何移动以及步长和步数：IBRION,NFREE,POTIM,NSW分子动力学相关参数：SMASS,TEBEG,TEEND,POMASS,NBLOCK,KBLOCK,PSTRESS原子驰豫收敛标准:EDIFFG定义态密度积分的方法和参数:smearing方法和参数：ISMEAR,SIGMA计算态密度时能量范围和点数:EMIN,EMAX,NEDOS计算分波态密度的参数：RWIGS,LORBIT其他：计算精度控制:PREC磁性计算:ISPIN,MAGMOM,NUPDOWN交换关联函数:GGA,VOSKOWN计算ELF和总的局域势:LELF,LVTOT结构优化参数:ISIF

@HouZhufeng3SYSTEM:注释所计算的体系，以示说明。NWRITE:默认值为2，可赋予值为0|1|2|3|4。决定OUTCAR中输出内容的详细程度ISTART:如果计算目录中有WAVECAR文件，则默认值为1，否则为0。可赋予值为0|1|2|3。决定是否读入WAVECAR：0:开始新的计算，按INIWAV初始化波函数1:接着计算，通常用在测试ENCUT的收敛性以及计算结合能曲线(也就是体积和总能的关系)2:接着计算，通常用在希望保持基矢不变的计算中3:接着计算，读入上一次计算得到的电荷密度和波函数，不推荐用ICHARG:如果ISTART=0，则默认值为2，否则为0。可赋予值为0|1|2。决定了如何构造初始的电荷密度0:从初始的波函数构造1：读入CHGCAR读入，并同原子密度进行线性插值2:构造原子密度11：读入自洽的CHGCAR，并进行能带计算或态密度的非自洽计算12:非自洽的原子密度计算初始化

@HouZhufeng4INIWAV:默认值为1，可赋予值为0|1。只在开始新的计算(也就是ISTART=0)中有效。决定了如何初始化波函数0:采用”jellium波函数”1:对波函数赋予随机数NELM:允许电子自洽迭代的最大步数。默认值为60。如果超过了40步还没有收敛的话，推荐对IALGO、LDIAG和混合参数进行手动设置到合理的值NELMIN:电子自洽迭代的最小步数。默认值为2。通常不需手动设置，在有些情况可以赋予更大数NELMDL:在一开始计算时初始化的过程中电子非自洽迭代的步数。当ISTART=0,INIWAV=1和IALGO=8时默认值为-5。当ISTART=0,INIWAV=1和IALGO=48(vasp4.4版本中)时默认值为-12。其他情况下为0EDIFF:在电子自洽迭代的循环中，总能收敛的标准。默认值为10-4电子自洽收敛

@HouZhufeng5EDIFFG:原子迟豫收敛的标准。默认值为EDIFF*10。如果它的值为正的，则表示前后两次的总自由能之差小于EDIFFG，原子停止迟豫。如果为负的，则原子所受的最大的力小于EDIFFG的绝对值，原子停止迟豫。NSW:原子迟豫的最大步数和分子动力学的步数。默认值为0。在每一步内，电子进行自洽计算，并精确计算原子所受的H-F力和和应力NBLOCK:默认值为1。原子每迟豫NBLOCK步后，将计算对关联函数、态密度以及输出原子迟豫后的位置(XDATCAR)。如果SMASS=-1，它将控制原子的动能在每多少步后进行缩放。KBLOCK：默认值为NSW。每KBLOCK*NBLOCK步后，将输出对关联函数和态密度。原子迟豫

@HouZhufeng6原子迟豫IBRION:决定了原子如何移动或迟豫。如果NSW=0或1，则默认值为－1，否则为0。可赋予值为－1|0|1|2|3|5。IBRION=0表示进行分子动力学模拟。为-1表示原子不移动。为其他值设置结构优化的方法。-1：原子位置不移动0:标准的分子动力学模拟。采用Verlet算法来积分原子的牛顿运动方程。通过POTIM来控制时间步长(单位是fs)。SMASS控制系综的设置1:采用准牛顿算法来优化原子的位置2