VASP参数设置详解
计算材料 2010-11-30 20:11:32 阅读197 评论0 字号:大中小 订阅
转自小木虫,略有增减
软件主要功能:
采用周期性边界条件(或超原胞模型)处理原子、分子、团簇、纳米线(或管)、薄膜、晶体、准晶和无定性材料,以及表面体系和固体
l 计算材料的结构参数(键长、键角、晶格常数、原子位置等)和构型 l 计算材料的状态方程和力学性质(体弹性模量和弹性常数)
l 计算材料的电子结构(能级、电荷密度分布、能带、电子态密度和ELF) l 计算材料的光学性质 l 计算材料的磁学性质
l 计算材料的晶格动力学性质(声子谱等) l 表面体系的模拟(重构、表面态和STM模拟) l 从头分子动力学模拟
l 计算材料的激发态(GW准粒子修正)
计算主要的四个参数文件:INCAR ,POSCAR,POTCAR ,KPOINTS,下面简要介绍,详细权威的请参照手册 INCAR文件:
该文件控制VASP进行何种性质的计算,并设置了计算方法中一些重要的参数,这些参数主要包括以下几类:
? 对所计算的体系进行注释: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 – 等等。
主要参数说明如下:
? 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 0
0:由初始波函数计算电荷密度 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-随机数),缺省值:1。