VASP参数设置详解报告

1、VASP 参数设置详解计算材料2010-11-30 20:11:32阅读 197 评论 0字号： 大中小 订阅转自小木虫，略有增减软件主要功能：采用周期性边界条件（或超原胞模型）处理原子 、分子 、团簇、 纳米线 （或管）、薄膜、晶体 、准晶和无定性材料，以及表面体系和固体l 计算材料的结构参数（键长、键角、晶格常数、原子位置等）和构型l 计算材料的状态方程和力学性质（体弹性模量和弹性常数）l计算材料的电子结构（能级、电荷密度分布、能带、电子态密度和ELF ）l 计算材料的光学性质l 计算材料的磁学性质l 计算材料的晶格动力学性质（声子谱等）l表面体系的模拟（重构、表面态和STM 模拟）l 从

2、头分子动力学模拟l计算材料的激发态（GW 准粒子修正）计算主要的四个参数文件：INCAR ,POSCAR,POTCAR ,KPOINTS,下面简要介绍，详细权威的请参照手册INCAR 文件：该文件控制VASP 进行何种性质的计算，并设置了计算方法中一些重要的参数，这些参数主要包括以下几类：对所计算的体系进行注释：SYSTEM定义如何输入或构造初始的电荷密度和波函数：ISTART ， ICHARG ， INIWA V定义电子的优化平面波切断动能和缀加电荷时的切断值：ENCUT ， ENAUG电子部分优化的方法：ALGO ， IALGO ， LDIAG电荷密度混合的方法：IMIX ， AMIX ，

3、 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

4、， EMAX ， NEDOS计算分波态密度的参数：RWIGS ， LORBIT其它计算精度控制： PREC磁性计算： ISPIN ， MAGMOM， NUPDOWN交换关联函数： GGA ， VOSKOWN计算 ELF 和总的局域势：LELF ， LVTOT结构优化参数： ISIF等等。主要参数说明如下：SYSTEM ：该输入文件所要执行的任务的名字。取值：字符串，缺省值：SYSTEMNWRITE ：输出内容详细程度。取值：04，缺省值：2如果是做长时间动力学计算的话，最好选0或1(首末步 /每步核运动输出 )，短时运算用 2，选3则会在出错的时候给出说明信息。ISTART ：决定是否读取WA

5、VECAR文件。取值：02，缺省 0/1 for 无 /有前次计算的 WAVECAR （波函数）0： begin 'from scratch' ，根据 INIWA V 初始化波函数1：restart with constant energy cut-off ，从 WAVECAR读取波函数 （重定义平面波集）2： restart with constant basis set ，从 WAVECAR读取波函数（平面波集不变）ICHARG ：决定如何建立初始电荷密度。取值：02，缺省值 : if ISTART=0 2 else 00：由初始波函数计算电荷密度1：从 CHGCAR 文件

6、读取电荷密度2：使用原子电荷密度的叠加+10非自洽计算ISPIN ：是否进行spin polarized calculation 。取值： 1， 2（ 1-no ， 2-yes ），缺省值：2MAGMOM：在 ICHARG=2 或在 CHGCAR 中未包含磁化密度（ICHARG=1 ）时，指定每个原子的初始磁化时刻。取值：实数数组，缺省值对非共线型磁化体系3*NIONS*1.0:对ISPIN=2 NIONS*1.0，INIWA V 如何设置初始波函数，只在ISTART=0 时使用。取值：0， 1（ 0-最低动能的平面波，1-随机数），缺省值： 1。IDIPOL控制计算单极、偶极和四极修正。取值

7、：14。13只计算第一 /二/三晶矢方向，适于厚板（slab）的计算4 所有方向都计算，适于计算孤立分子PREC 进动（ precession）。取值： low/medium/high/normal/ accurate/single)，缺省值 :Normal （ VASP.4.X ）；Medium （ VASP.5.X ）VASP4.5+采用了优化的 accurate来替代 high，所以一般不推荐使用 high。不过 high 可以确保“绝对收敛”，作为参考值有时也是必要的。同样受推荐的是normal，适于作为日常计算选项。受 PREC 影响的参数有四类： ENCUT ；NGX ，NGY ，

8、NGZ；NGXF ，NGYF ，NGZF；ROPT。如果设置了 PREC，这些参数就都不需要出现了，当然直接设置相应的参数也有同样效果。具体影响效果见p5354。ENCUT 平面波基组的截断能量（eV ）。取值：实数，缺省值：受PREC 设置影响，从 POTCAR 文件中找出相应的 ENMAX/ENMIN 值来设置。PREC =LowMediumAccurateNormalENCUT =ENMINENMAXENMAXENMAXSingleHighENMAXENMAX*1.3对于多个元素不同的ENMAX/ENMIN ，都取最大值。该参数非常重要， 最好不要手工去设置， 除非文献告诉你要用多少，或

9、者经过结果可靠性的验证。当然，为了测试一下提交的任务，也不妨先设个较小的值。NGX ， NGY ， NGZ ：控制FFT 网格在三个晶矢方向上的格点数量。NGFX ， NGFY ， NGFZ ：控制第二次更精确的FFT 网格的格点数量。也是两类重要的最好不要去动的参数。 在未指定的情况下将根据 PREC 的设置从 POTCAR 中自动读取。PREC=High/Accurate，基组中向量的 2倍值，用来避免wraparound errors，得到精确解。PREC=Low/Medium/Normal ，基组中向量的 3/4倍值（已足够精确到 1 meV/atom）。LREAL: 决定投射是在实空

10、间还是倒易空间进行。取值： .TRUE.（实空间） /.FALSE.（倒易空间） ，缺省值：.FALSE.用于求解赝势的非局域部分用到的一个积分，在倒格空间里采用平面波基组求解，在实空间里则采用积分球求解。其他还有两个选项： O or On，A or Auto。On 和.TRUE.的差别在于是否使用King-Smith 算法优化，设为 Auto 则进行自动选择，推荐使用。ROPT: 在 LREAL=Auto or On时，优化控制每个核周围的积分球内的格点数。取值：实数数组For LREAL=OnPREC=Low, 700 points in the real space sphere (RO

11、PT=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=AutoPREC=Low, accuracy 10 -2 (ROPT=0.01)-3PREC=Med, accuracy 2*10(ROPT=0.002)PREC=High accuracy 2*10 -4 (ROPT=2E-4)NELM, NELMIN and NELMDL：控制电子自洽循环步数。取值：整数NELM ：电子自洽循环

12、最大次数。缺省值：60NELMIN ：电子自洽循环最小次数。缺省值：2NELMDL ：弛豫次数。缺省值：if ISTART=0, INIW AV=1, andIALGO=8 ，-5，if ISTART=0, INIW AV=1, and IALGO=48 ，-12，else 0NELMDL 可以取负值。如果初始波函数采用随机赋值，即ISTART=0，INIW AV=1，那么很可能开始的值比较离谱，那么在第一步核运动循环之前采用 NELMDL（负值）步的非自洽（保留初始的 H）步计算将减少计算所需的时间。EDIFF ：指定电子自洽循环的全局中断条件，用于控制收敛精度。取值：实数，缺省值： 10

13、-4注意，即使 EDIFF=0 ，NELM 步电子自洽循环也会执行。EDIFFG ：指定离子弛豫循环的中断条件，用于控制核运动的收敛精度。取值，实数，缺省值：10*EDIFFEDIFFG>0 在两个离子步的总自由能之差小于EDIFFG 时停止EDIFFG<0 在所有的力都小于EDIFFG 时停止。EDIFFG=0 在 NSW 步弛豫后停止此参数不支持MD ，仅用于弛豫。NSW ：给出最大离子步数。取值：整数，缺省值：0。NBLOCK， KBLOCK：取值：整数，缺省值：NBLOCK = 1， KBLOCK = NSW在NBLOCK离子步后对成对相关函数和DOS进行计算， 并且把离子

14、配置写入XDATCAR文件。在 KBLOCK*NBLOCK步主循环后平均的成对相关函数和DOS 被写入 PCDAT和 DOSCAR 文件。IBRION ：决定离子怎样更新和运动。取值：-13 ， 58（ -1- 无更新， 0-MD ，1-RMM-DIIS， 2-共轭梯度算法，3-Damped MD ，5,6：有限差分，7,8：密度函数扰动理论），缺省值：if NSW=0/1 ， -1， else 0这个参数是和ISIF, IALGO/ALGO一起决定怎么算的最重要的参数。13是三种弛豫的方法，根据ISIF 决定是否固定离子位置、晶胞大小和形状，在 INCAR 中必须设置参数 POTIM 。0是

15、标准的ab-initio MD ，不受ISIF 影响，即不改变晶胞大小和形状。58支持 Hessian Matrix 和 phonon frequency 的计算以及部分固定的MD 。POTIM：IBRION=0时，给出MD每步步长（fs）， IBRION=13时，给出最小化的度量常量。取值：实数，缺省值：IBRION=0无缺省，必须指定，IBRION=1,2,3 0.5ISIF ：决定是否计算应力张量以及弛豫中晶胞变化的自由度。取值：06，缺省值：if IBRION=0(MD) 0 else 2ISIF calculatecalculate relaxchangechange force s

16、tress tensor ions cell shapecell volume- -0yesno yesnono1yes trace only yes no no2yesyes yes nono3yesyes yes yesyes4yesyes yes yesno5yesyes no yesno6yesyes no yesyes7yesyes no noyestrace only means that only the total pressure is correctIWAVPR ：决定波函数和/或电荷密度怎样从一个离子配置向下一个离子配置进行推测。取值：03， 1013（ 0-无推测， 1

17、,11-用原子电荷密度进行简单推测，2,12- 二阶推测， 3,13-混合前两种方法） ，缺省值： ifIBRION=0(MD),1,2(relaxation) 2 else(静态计算 )0推测结果保存在外部文件TMPCAR中，取值 +10则全部使用内存，不保存此文件。ISYM ：决定是否使用对称性。取值：-13 （ -1,0-不使用， 1,2,3-使用），缺省值：if 使用 US-PP 1 ， if 使用 PAW 2ISYM=2 使用一种效率更高也更节省内存的电荷密度对称性，ISYM=3 时仅考虑力和应力张量的对称性，而电荷密度是非对称的。SYMPREC ：决定 POSCAR 文件中给出的位

18、置的精度。取值：实数，缺省值：10-5LCORR ：决定是否对非完全自洽计算中的力进行Harris 修正。取值： .TRUE./.FALSE. ，缺省值：.TRUE.TEBEG, TEEND ：控制从头分子动力学计算中的起始温度和最终温度（ MD 有效）。取值：实数，缺省值：TEBEG = 0 TEEND = TEBEG注意 VASP 的温度定义与实际温度有细微的差别，所以 TEBEG=T × (N-1)/N ，T为实际要求的温度，N 为原子数。SMASS ：控制从头MD 中的速度。取值：-30 ，缺省值：-3微正则系综（总自由能不变），-2 保持初速度不变， -1每 NBLOCK步

19、调整速度， 来保证动能连续， >=0 Nos é算法模拟正则系综NPACO ：成对相关函数的槽数。取值：整数，缺省值：256APACO ：成对相关函数求值中的最大距离（?）。取值：整数，缺省值：16简单说就是在不超过APACO 的 NPACO 个距离上求成对相关函数PCF。RWIGS ：给出 Wigner-Seitz 半径，DOS 计算用。取值：实数数组，缺省值：从 POTCAR文件中读取NELECT ：总电子数， 如果系统不是电中性的就必须设置，所带电荷作为均一的背景电子气考虑。取值：实数，缺省值：- （价电子数） ，由 POSCAR 和 POTCAR 文件自动决定（通常不必

20、给出）。NUPDOWN上下自旋成分间的电子数之差。取值：整数，缺省值：未设置（此时将进行完全弛豫）EMIN, EMAX：DOS 求值的最小 /最大能量。 取值：实数，缺省值： EMIN = - (lowestKS-eigenvalue -)， EMAX =- (highest KS-eigenvalue -)ISMEAR ：决定每个波函数的部分占位fnk 如何设置。 取值： -5 | -4 | -3 | -2 | 0 | N（ -5-带有 Blochl 修正的四面体方法, -4-不带 Blochl 修正的四面体方法, -3- 根据 INCAR 文件中提供的smearing 参数执行循环，-2-

21、 从 WAVECAR 文件中读取,-1-Fermi-smearing, 0-Gaussian smearing, >0-method of Methfessel-Paxton orderN，缺省值： 1采用部分占位波函数，用一个函数来平滑积分，尤其是对于金属体系可减少k点。SIGMA ：决定 smearing 的宽度（ eV）。取值：实数，缺省值：0.2ALGO ：指定电子最小化算法。取值：Normal(blocked Davidson block iterationscheme)/VeryFast(RMM-DIIS)/Fast (前两个算法的混合)/All( 波函数的所有带同时更新 )

22、/ Damped(damped velocity friction algorithm)，缺省值 : NormalIALGO ：指定主算法（整数选择算法）。取值： 8( 共轭梯度算法)/38(Davidson blockiteration scheme)/48(RMM-DIIS)，缺省值：38算法是最重要的参数之一。一般VASP 推荐使用的是以上三种算法，一般来说8/38是初期比较快收敛，在接近平衡时采用48较快，在初期或MD 时使用 48可能会遇到不收敛的情况。也可以使用ALGO 参数来替代IALGO ，设置 Fast，VASP 会先用 38，再自动切换到 48 。各种算法只要收敛，结果应该

23、一致。另一个可能有用的选项是 -1 。不进行实际的计算， 只对重要的步骤做计算测试，并将测试得到的各部分耗时输出到OUTPUT 中。VOSKOWN ：决定是否使用VWN 插值算法。取值：0( 不使用 )/1( 使用 )，缺省值：0如果使用了PW91 泛函或需要计算磁性质时可以设为1使用。Mixing-tags ：IMIX ：混合的类型，取值：整数，缺省值：4AMIX ：线性混合参数。取值：实数，缺省值：0.8（ US-PP）， 0.4（ PAW）AMIN ：最小混合参数。取值：实数，缺省值：0.1BMIX ： Kerker 混合方案的截断波向量。取值：实数，缺省值：1.0AMIX_MAG ：磁

24、化过程的线性混合参数。取值：实数，缺省值：1.6BMIX_MAG：磁化过程的Kerker 混合方案的截断波向量。取值：实数，缺省值： 1.0WC ：Broyden 混合方案中每步的加权因子。取值：实数，缺省值：1000.0INIMIX ： Broyden 混合方案中的初始混合类型。取值：整数，缺省值：1MIXPRE ： Broyden 混合方案中的预处理类型。取值：整数，缺省值：1MAXMIX ： Broyden 混合器中存储的最大步数。取值：整数，缺省值：-45值得注意的是，在MD 或者弛豫的时候，设置MAXMIX （ >0,一般约 3倍的电子 SC 步数）可能会大大减少核运动步数，但

25、同时也会增加对内存的要求。LWAVE,LCHARG：决定是否把波函数（或电荷密度）写入外部文件WAVECAR（或 CHGCAR和 CHG ）中。取值：.TRUE./.FALSE. ，缺省值：.TRUE.LVTOT ：决定是否把总局域势写入外部文件LOCPOT 中。取值： .TRUE./.FALSE. ，缺省值： . FALSE.LELF ：决定是否创建ELFCAR 文件。取值：.TRUE./.FALSE. ，缺省值：. FALSE.ELFCAR用于保存ELF(electron localization function)。LORBIT ：和适当的RWIGS 一起决定是否创建PROCAR 或 P

26、ROOUT 文件。取值：0 | 1 | 2 | 5 | 10 | 11 | 12，缺省值：0(创建 DOSCAR 和 PROCAR 文件 )NPAR ：用于控制VASP 切换到带间并行模式。取值：整数，缺省值：总节点数。每一个节点计算一个带，可以提高并行效率，减少通讯量，但可能会大幅增加内存的需求。NBANDS ：给出计算中实际的总能带数。取值：整数，缺省值： NELECT/2+NIONS/2（非自旋）， 0.6*NELECT + NMAG（自旋）NBANDS的取值应使计算中包含相当数量的空带。因为计算需要大量的空带，至少要求 1个空带（否则VASP 会给出警告） 。NBANDS对于解决内存需

27、求非常重要。一般 NBANDS在 NELECT/2+NIONS/2以上可以得到较为精确的结果，但如果内存不够就只能减少NBANDS ，在牺牲精度和体系大小之间平衡了。POTCAR 文件赝势文件，最重要的输入文件之一。可以理解为分子力学模拟中的力场文件，但包含的信息更多。 VASP 将各元素优化的INCAR里的参数也包含在该文件中，作为支持PREC 的缺省选择。通常各元素的POTCAR 文件已经包含在VASP 软件包中的赝势库里了，用户需要做的是确定自己具体需要哪几种赝势，然后按照POSCAR 文件里的顺序，将所选择的POTCAR 文件按顺序连接起来就可以了。如以下命令：cat file1 fi

28、le2 file3 > POTCAR赝势库中的赝势文件可以分为以下几类：根据产生方法的不同有Ultra-soft 赝势（ USPP ）和投影扩充波赝势（PAW）；根据交换关联函数的不同有LDA（ local density approximation）和 GGA（ generalizedgradient approximation ，又可以再分为PW91 和 PBE ）；根据半芯态的不同元素A 的赝势文件还可以分为A ， A_sv 和 A_pv ；根据选取的不同截断能量（ENMAX ）元素A 的赝势文件还可以分为A ， A_s 和A_h 。VASP 的赝势文件放在目录/vasp/pote

29、ntials 下，该目录又包含五个子目录：pot ，pot_GGA ，potpaw ，potpaw_GGA ，potpaw_PBE ，其中每个子目录对应一种赝势形式：pot => PP, LDA ；pot_GGA => PP, GGA ； potpaw => PAW, LDA ； potpaw_GGA => PAW, GGA, PW91；potpaw_PBE => PAW, GGA, PBE。每个目录中每种元素根据截断能量和半芯态的不同还会有多个对应的赝势文件存在。在具体选取时可以参考各版本同目录下的V_RHFIN和 PSCTR 文件， 这两个文件说明了该版本的赝

30、势是如何生成的。选好所要使用的赝势之后，进入对应的目录，会看到里面有四个文件：POTCAR.Z ，PSCTR.Z ，V_RHFIN.Z和 WS_FTP.LOG 。现在需要用到的是第一个，把各元素的该文件解压（ zcat POTCAR.Z > file），然后用 cp 或 mv 命令把这些文件移到工作目录里，再用 cat 命令把它们合并到POTCAR 中，就得到了我们需要的POTCAR 。注意要记住这里元素的排列顺序，以后在 POSCAR 文件中各元素的排列就是按照这个顺序。POSCAR 文件位置文件。描述所计算体系的晶胞参数、原子个数及晶胞中原子的位置，以及分子动力学计算时原子的初始速度

31、。POSCAR 文件示例：Si-fcc5.430.00.50 0.500.50 0.00 0.500.50 0.50 0.002Direct0.00.00 0.000.25 0.25 0.25POSCAR 文件的结构：第1行：任意文字注释。第2行：晶格常数（ ?），也是晶矢的缩放系数，后面所有长度值得自原值除以此值。a=b=c时取 a 即可，否则一般取三者最大值，若取负值，则为晶胞体积（?3）。第3-5行：定义晶矢。第6行：每种元素的原子个数，特别注意顺序,要与下面的坐标顺序以及POTCAR 中的顺序一致。第7行：可省略，无需空行。做动力学时，用于指定是否需要固定部分原子的坐标。若是，此行以S

32、或者 s作为首字母即可。第8行开始为原子的坐标位置，格式为option linecoordinate1 of element1coordinate2 of element1.coordinateN of element1option linecoordinate1 of element2coordinate2 of element2.coordinateM of element2.其中， option line 指定输入坐标的格式，除了第一个以外，如果后面的输入格式同前，则都可以无空行省略。option line 可指定的输入坐标格式有两种：D or d for direct modeC or

33、 cor K or k for cartesian mode顾名思义，前者是定义在三个晶矢方向上的坐标：R=R1× x+ R2×y +R3×z，R1,R2,R3 为前面的晶矢， x,y,z 为输入的三个坐标，R 为原子坐标位矢，是相对于晶矢给出的，而后者则是以笛卡尔坐标系来给出原子的绝对坐标，实际就是将直角坐标除以前面第二行定义的晶胞常数后得到的坐标值。如果第 7行设定为S(Selective Dynamic) ，则可以用以下形式定义各坐标是否可以移动：Selective dynamicsCartesian0.00 0.00 0.00 T T F0.25 0.25

34、 0.25 F F FKPOINTS 文件设置布里渊区k 点网格取样大小或能带结构计算时沿高对称方向的k 点。一般来说， k 点越密越多，计算精度也就越高，但计算成本也会越高。因此，对于原子数较多的体系的计算，需要谨慎的尝试k 点数目，在避免或者预先评估wrap-around error 的前提下尽量减少k 点数目。KPOINTS 文件示例：k-points along high symmetry lines11Reciprocal0.00 0.00 0.00 1.000.05 0.00 0.05 1.000.05 0.00 0.50 1.00KPOINTS 文件结构：第1行：注释行。第2行：

35、 k 点总数，或者0（自动生成网格（Automatic k-mesh generation ）。如果是前者，给出k 点总数，又分为两种情况：M ：全手动（ Entering all k-points explicitly），手动输入即自定义各个k 点的坐标和权重，推荐只在能带计算时使用。此时：第3行：输入格式标识。直角坐标（Cartesian）或倒易坐标（Reciprocal ）。同样的，c、 C、 k、 K for Cartesian ，其他首字母则自动切换到Reciprocal 。第4n 行 :逐个 k 点的坐标描述。格式为“x y z W ”。其中“ x y z ”是该点的三个坐标，W是权重。所有k 点的权重相互之间的比例对了就行，VASP 会自动归一的。注意 C 坐标和R 坐标的定义：C： k=(2 /a)(x y z)R： k=x*b1+y*b2+z*b3（ b1-b3 为倒易晶矢）这里“ x y z ”只是代表了坐标的顺序，与