VASP参数设置详解报告 下载本文

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。