VASP全参数设置详解

1、实用标准文档VAS哆数设置详解软件主要功能：采用周期性边界条件（或超原胞模型）处理原子、分子、团簇、纳米线（或管）、薄膜、晶体、准晶和无定性材料，以及表面体系和固体i计算材料的结构参数（键长、键角、晶格常数、原子位置等）和构型i计算材料的状态方程和力学性质（体弹性模量和弹性常数）l计算材料的电子结构（能级、电荷密度分布、能带、电子态密度和ELF）l计算材料的光学性质l计算材料的磁学性质l计算材料的晶格动力学性质（声子谱等）l表面体系的模拟（重构、表面态和STM模拟）l从头分子动力学模拟l计算材料的激发态（GW隹粒子修正）计算主要的四个参数文件：INCAR,POSCAR,POTCAR,KPOIN

2、TS下面简要介绍，详细权威的请参照手册INCAR文件：该文件控制VASP®行何种性质的计算，并设置了计算方法中一些重要的参数，这些参数主要包括以下几类：对所计算的体系进行注释：SYSTEM定义如何输入或构造初始的电荷密度和波函数：ISTART,ICHARGINIWAV定义电子的优化平面波切断动能和缀加电荷时的切断值：ENCUTENAUG电子部分优化的方法：ALGOIALGO,LDIAG电荷密度混合的方法：IMIX,AMIX,AMIN,BMIX,AMIX_MAGBMIX_MAG文案大全实用标准文档WCINIMIX,MIXPREMAXMIX自洽迭代步数和收敛标准：NELMNELMINNE

3、LMDLEDIFF定义离子或原子的优化原子位置优化的方法、移动的步长和步数：BRION,NFREEPOTIM,NSW分子动力学相关参数：SMASSTEBEGTEENDPOMASSNBLOCKKBLOCKPSTRESS离子弛豫收敛标准：EDIFFG定义态密度积分的方法和参数smearing方法和参数：ISMEARSIGMA计算态密度时能量范围和点数：EMIN,EMAXNEDOS计算分波态密度的参数：RWIGSLORBIT其它计算精度控制：PREC磁性计算：ISPIN,MAGMOMNUPDOWN交换关联函数：GGAVOSKOWN计算ELF和总的局域势：LELF,LVTOT结构优化参数：ISIF一等

4、等。主要参数说明如下：SYSTEM该输入文件所要执行的任务的名字。取值：字符串，缺省值：SYSTEMNWRITE输出内容详细程度。取值：04,缺省值：2如果是做长时间动力学计算的话，最好选0或1（首末步/每步核运动输出），短时运算用2,选3则会在出错的时候给出说明信文案大全实用标准文档ISTART:决定是否读取WAVECAR件。取值：02,缺省0/1for无/有前次计算的WAVECAR波函数）0:begin'fromscratch',根据INIWAV初始化波函数1:restartwithconstantenergycut-off,从WAVECAR取波函数（重定义平面波集）2：r

5、estartwithconstantbasisset,从WAVECA犊取波函数（平面波集不变）ICHARG决定如何建立初始电荷密度。取值：02,缺省值：ifISTART=02else00:由初始波函数计算电荷密度1:从CHGCAR：件读取电荷密度2:使用原子电荷密度的叠加+10非自洽计算ISPIN:是否进行spinpolarizedcalculation。取值：1,2（1-no,2-yes）,缺省值：2MAGMOMBICHARG=2£在CHGCA即未包含磁化密度（ICHARG=1时，指定每个原子的初始磁化时刻。取值：实数数组，缺省值：对ISPIN=2NIONS*1.0,对非共线型磁化

6、体系3*NIONS*1.0INIWAV如何设置初始波函数，只在ISTART=0时使用。取值：0,1（0-最低动能的平面波，1-随机数），缺省值：1OIDIPOL控制计算单极、偶极和四极修正。取值：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。不

7、过high可以确保“绝对收敛”，作为参考值有时也是必要的。同样受推荐的是normal,适于作为日常计算选项。受PRE澎响的参数有四类：ENCUTNGXNGYNGZNGXFNGYFNGZFROPT如果设置了PREC这些参数就都不需要出现了，当然直接设置相应的参数也有同样效果。具体影响效果见p5354。ENCU"面波基组的截断能量(eV)。取值：实数，缺省值：受PRE彼置影响，从POTCA就件中找出相应的ENMAX/ENMIN1来设置。PREC=LowMediumAccurateNormalENCUT=ENMINENMAXENMAXENMAXSingleHighENMAXENMAX*1.

8、3对于多个元素不同的ENMAX/ENMJNTB取最大值。该参数非常重要，最好不要手工去设置，除非文献告诉你要用多少，或者经过结果可靠性的验证。当然，为了测试一下提交的任务，也不妨先设个较小的值。NGXNGYNGZ控制FFT网格在三个晶矢方向上的格点数量。NGFXNGFYNGFZ控制第二次更精确的FFT网格的格点数量文案大全实用标准文档也是两类重要的最好不要去动的参数。在未指定的情况下将根据PREC勺设置从POTCA中自动读取。PREC=High/Accurate,基组中向量的2倍值，用来避免wraparounderrors,得至U精确解。PREC=Low/Medium/Normal基组中向量的

9、3/4倍值(已足够精确到1meV/atom)。LREAL:决定投射是在实空间还是倒易空间进行。取值：.TRUE.(实空间)/.FALSE.(倒易空间)，缺省值：.FALSE.用于求解震势的非局域部分用到的一个积分，在倒格空间里采用平面波基组求解，在实空间里则采用积分球求解。其他还有两个选项：OorOn,AorAuto。On和.TRUE.的差别在于是否使用King-Smith算法优化，设为Auto则进行自动选择，推荐使用。ROPT在LREAL=AutoorOn时，优化控制每个核周围的积分球内的格点数。取值：实数数组ForLREAL=OnPREC=Low,700pointsintherealspa

10、cesphere(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)文案大全实用标准文档NELM,NELMINandNELMDL:控制电子自洽循环步数。取值：整数NELM电子自洽循环最大次数。缺省值：60

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

12、子弛豫循环的中断条件，用于控制核运动的收敛精度。取值，实数，缺省值：10*EDIFFEDIFFG驻两个离子步的总自由能之差小于EDIFFG时停止EDIFFG驻所有的力都小于EDIFFG时停止。EDIFFG=ONSWH也豫后停止此参数不支持MD仅用于弛豫。NSW给出最大离子步数。取值：整数，缺省值：0。文案大全实用标准文档NBLOCKKBLOCK取值：整数，缺省值：NBLOCK=1,KBLOCK=NSW在NBLOCK子步后对成对相关函数和DOS进行计算，并且把离子配置写入XDATCA或件。在KBLOCK*NBLOC发主循环后平均的成对相关函数和DOSa写入PCDATADOSCARC件。IBRIO

13、N:决定离子怎样更新和运动。取值：-13,58（-1-无更新，0-MD,1-RMM-DIIS,2-共朝梯度算法，3-DampedMD,5,6:有限差分，7,8:密度函数扰动理论），缺省值：ifNSW=0/1,-1,else0这个参数是和ISIF,IALGO/ALGO一起决定怎么算的最重要的参数。13是三种弛豫的方法，根据ISIF决定是否固定离子位置、晶胞大小和形状，在INCAR中必须设置参数POTIM0是标准的ab-initioMD,不受ISIF影响，即不改变晶胞大小和形状。58支持HessianMatrix和phononfrequency的计算以及部分固定的MDPOTIM:旧RION=0时，

14、给出MD每步步长（fs）,IBRION=13时，给出最小化的度量常量。取值：实数，缺省值：旧RION=0无缺省，必须指定，IBRION=1,2,30.5ISIF:决定是否计算应力张量以及弛豫中晶胞变化的自由度。取值：06,缺省值：ifIBRION=0（MD）0else2relaxchangeionscellshapecellISIFcalculatecalculatechangeforcestresstensorvolume文案大全实用标准文档no0yesnoyesno1yestraceonlyyesnono2yesyesyesnono3yesyesyesyesyes4yesyesyesyes

15、no5yesyesnoyesno6yesyesnoyesyes7yesyesnonoyestraceonlymeansthatonlythetotalpressureiscorrectIWAVPR决定波函数和/或电荷密度怎样从一个离子配置向下一个离子配置进行推测。取值：03,1013（0-无推测，1,11-用原子电荷密度进行简单推测，2,12-二阶推测，3,13-混合前两种方法），缺省值：iHBRION=0（MD）,1,2（relaxation）2else（静态计算）0推测结果保存在外部文件TMPCARP,取彳t+10则全部使用内存，不保存此文件。ISYM:决定是否使用对称性。取值：-13（-

16、1,0-不使用，1,2,3-使用），缺省值:if使用US-PP1,if使用PAW2ISYM=2使用一种效率更高也更节省内存的电荷密度对称性，ISYM=3寸仅考虑力和应力张量的对称性，而电荷密度是非对称的。文案大全实用标准文档.、-5SYMPREC决定POSCARC件中给出的位置的精度。取值：实数，缺省值：10LCORR决定是否对非完全自洽计算中的力进行Harris修正。取值：.TRUE./.FALSE.,缺省值：.TRUE.TEBEG,TEEND控制从头分子动力学计算中的起始温度和最终温度（MDf效）。取值：实数，缺省值：TEBEG=0TEEND=TEBEG注意VASP的温度定义与实际温度有细

17、微的差别，所以TEBEG=T（N-1）/N,T为实际要求的温度，N为原子数。SMASS控制从头M叶的速度。取值：-30,缺省值：-3微正则系综（总自由能不变），-2保持初速度不变，-1每NBLOC处调整速度，来保证动能连续，>=0Nose算法模拟正则系综NPACO成对相关函数的槽数。取值：整数，缺省值：256APACO成对相关函数求值中的最大距离（？）。取值：整数，缺省值：16简单说就是在不超过APACO勺NPAC价距离上求成又t相关函数PCRRWIGS给出Wigner-Seitz半径，DOStf算用。取值：实数数组，缺省值：从POTCAR文件中读取NELECT总电子数，如果系统不是电中

18、性的就必须设置，所带电荷作为均一的背景电子气考虑。取值：实数，缺省值：-（价电子数），由POSCAR口POTCA戍件自动决定（通常不必给出）。NUPDOWN下自旋成分间的电子数之差。取值：整数，缺省值：未设置（此时将进行完全弛豫）文案大全实用标准文档EMIN,EMAXDOSt值的最小/最大能量。取值：实数，缺省值：EMIN=-(lowestKS-eigenvalue-A),EMAX=-(highestKS-eigenvalue-A)ISMEAR决定每个波函数的部分占位fnk如何设置。取值：-5|-4|-3|-2|0|N(-5-带有Blochl修正的四面体方法,-4-不带Blochl修正的四面体

19、方法,-3-根据INCAR文件中提供的smearing参数执行循环，-2-从WAVECAR件中读取，-1-Fermi-smearing,0-Gaussiansmearing,>0-methodofMethfessel-PaxtonorderN,缺省值：1采用部分占位波函数，用一个函数来平滑积分，尤其是对于金属体系可减少k点。SIGMA决定smearing的宽度(eV)。取值：实数，缺省值：0.2ALGO指定电子最小化算法。取值：Normal(blockedDavidsonblockiterationscheme)/VeryFast(RMM-DIIS)/Fast(前两个算法的混合)/All

20、(波函数的所有带同时更新)/Damped(dampedvelocityfrictionalgorithm),缺省值:NormalIALGO:指定主算法(整数选择算法)。取值：8(共朝梯度算法)/38(Davidsonblockiterationscheme)/48(RMM-DIIS),缺省值：38算法是最重要的参数之一。一般VASP推荐使用的是以上三种算法，一般来说8/38是初期比较快收敛，在接近平衡时采用48较快，在初期或MD寸使用48可能会遇到不收敛的情况。也可以使用ALG3数来替代IALGO,设置Fast,VASP会先用38,再自动切换到48。各种算法只要收敛，结果应该一致。另一个可能有

21、用的选项是-1。不进行实际的计算，只对重要的步骤做计算测试，并将测试得到的各部分耗时输出到OUTPUTKVOSKOW做定是否使用VWN1值算法。取值：0(不使用)/1(使用)，缺省值：0如果使用了PW91乏函或需要计算磁性质时可以设为1使用。文案大全实用标准文档Mixing-tags:IMIX:混合的类型，取值：整数，缺省值：4AMIX:线性混合参数。取值：实数，缺省值：0.8（US-PP）,0.4（PAWAMIN:最小混合参数。取值：实数，缺省值：0.1BMIX：Kerker混合方案的截断波向量。取值：实数，缺省值：1.0AMIX_MAG磁化过程的线性混合参数。取值：实数，缺省值：1.6BM

22、IX_MAG磁化过程的Kerker混合方案的截断波向量。取值：实数，缺省值：1.0WCBroyden混合方案中每步的加权因子。取值：实数，缺省值：1000.0INIMIX:Broyden混合方案中的初始混合类型。取值：整数，缺省值：1MIXPREBroyden混合方案中的预处理类型。取值：整数，缺省值：1MAXMIXBroyden混合器中存储的最大步数。取值：整数，缺省值：-45值得注意的是，在MD或者弛豫的时候，设置MAXMIX（0,一般约3倍的电子SC步数）可能会大大减少核运动步数，但同时也会增加对内存的要求。LWAVE,LCHAR映定是否把波函数（或电荷密度）写入外部文件WAVECA碱C

23、HGCAR和CHG中。取值：.TRUE./.FALSE.,缺省值：.TRUE.LVTOT决定是否把总局域势写入外部文件LOCPO用。取值：.TRUE./.FALSE.,缺省值：.FALSE.LELF:决定是否创建ELFCAR5：件。取值：.TRUE./.FALSE.,缺省值：.FALSE.ELFCA观于保存ELF（electronlocalizationfunction）。LORBIT:和适当的RWIGS-起决定是否创建PROCARPPROOUTC件。取值：0|1|2|5|10|11|12,缺省值：0（创建DOSCAR口PROCARC件）文案大全实用标准文档NPAR用于控制VASP切换到带间并

24、行模式。取值：整数，缺省值：总节点数。每一个节点计算一个带，可以提高并行效率，减少通讯量，但可能会大幅增加内存的需求。NBANDS给出计算中实际的总能带数。取值：整数，缺省值：NELECT/2+NIONS/2（非自旋），0.6*NELECT+NMAG（自旋）NBANDS勺取值应使计算中包含相当数量的空带。因为计算需要大量的空带，至少要求1个空带（否则VASP会给出警告）。NBANDS寸于解决内存需求非常重要。一般NBAND院NELECT/2+NIONS/科上可以得到较为精确的结果，但如果内存不够就只能减少NBANDS在牺牲精度和体系大小之间平衡了。POTCA皮件鹰势文件，最重要的输入文件之一。

25、可以理解为分子力学模拟中的力场文件，但包含的信息更多。VASP将各元素彳t化的INCAR里的参数也包含在该文件中，作为支持PREC的缺省选择。通常各元素的POTCA欣件已经包含在VASP软件包中的鹰势库里了，用户需要做的是确定自己具体需要哪几种震势，然后按照POSCARt件里的顺序，将所选择的POTCARC件按顺序连接起来就可以了。如以下命令：catfile1file2file3>POTCAR震势库中的鹰势文件可以分为以下几类：根据产生方法的不同有Ultra-soft鹰势（USPP和投影扩充波震势（PAW；根据交换关联函数的不同有LDA（localdensityapproximation

26、）和GGA（generalizedgradientapproximation,又可以再分为PW9和PBE>；根据半芯态的不同元素A的鹰势文件还可以分为A,A_sv和A_pv；根据选取的不同截断能量（ENMAX元素A的震势文件还可以分为A,A_s和A_h。文案大全实用标准文档VASP勺鹰势文件放在目录/vasp/potentials下，该目录又包含五个子目录：pot,pot_GGA,potpaw,potpaw_GGApotpaw_PBE,其中每个子目录对应一种震势形式：pot=>PP,LDA；pot_GGA=>PP,GGApotpaw=>PAW,LDA；potpaw_GG

27、A=>PAW,GGA,PW91；potpaw_PBE=>PAW,GGA,PBE每个目录中每种元素根据截断能量和半芯态的不同还会有多个对应的震势文件存在。在具体选取时可以参考各版本同目录下的V_RHFIN和PSCTRC件，这两个文件说明了该版本的鹰势是如何生成的。选好所要使用的震势之后，进入对应的目录，会看到里面有四个文件：POTCAR.ZPSCTRZV_RHFIN.Z和WS_FTP.LOG现在需要用到的是第一个，把各元素的该文件解压(zcatPOTCAR.Z>file)，然后用cp或mv命令把这些文件移到工作目录里，再用cat命令把它们合并到POTCARK就得到了我们需要的P

28、OTCAR注意要记住这里元素的排列顺序，以后在POSCARt件中各元素的排列就是按照这个顺序。POSCAR件位置文件。描述所计算体系的晶胞参数、原子个数及晶胞中原子的位置，以及分子动力学计算时原子的初始速度。POSCARC件示例：Si-fcc5.430.00.500.500.500.000.500.500.500.002Direct0.00.000.00文案大全实用标准文档5POSCARC件的结构：第1行：任意文字注释。第2行：晶格常数（？），也是晶矢的缩放系数，后面所有长度值得自原值除以此值。a=b=c时取a即可，否则一般取三者最大值，若取负值，则为晶胞体积（？3）。第

29、3-5行：定义晶矢。第6行：每种元素的原子个数，特别注意顺序，要与下面的坐标顺序以及POTCA珅的顺序一致。第7行：可省略，无需空行。做动力学时，用于指定是否需要固定部分原子的坐标。若是，此行以S'或者's'作为首字母即可。第8行开始为原子的坐标位置，格式为optionlinecoordinate1ofelement1coordinate?ofelement1.coordinateNofelement1optionlinecoordinate1ofelement2coordinate2ofelement2.coordinateMofelement2.其中，optionl

30、ine指定输入坐标的格式，除了第一个以外，如果后面的输入格式同前，则都可以无空行省略。optionline可指定的输入坐标格式有两种：'D'or'd'fordirectmode文案大全实用标准文档'C'or'c'or'K'or'k'forcartesianmode顾名思义，前者是定义在三个晶矢方向上的坐标：R=R1Xx+R2Xy+R3Xz,R1,R2,R3为前面的晶矢，x,y,z为输入的三个坐标，R为原子坐标位矢，是相对于晶矢给出的，而后者则是以笛卡尔坐标系来给出原子的绝对坐标，实际就是将直角坐标除

31、以前面第二行定义的晶胞常数后得到的坐标值。如果第7行设定为S(SelectiveDynamic)，则可以用以下形式定义各坐标是否可以移动：SelectivedynamicsCartesian0.000.000.00TTF5FFFKPOINTSC件设置布里渊区k点网格取样大小或能带结构计算时沿高对称方向的k点。一般来说，k点越密越多，计算精度也就越高，但计算成本也会越高。因此，对于原子数较多的体系的计算，需要谨慎的尝试k点数目，在避免或者预先评估wrap-arounderror的前提下尽量减少k点数目。KPOINTSC件示例：k-pointsalonghighsymmetr

32、ylines11Reciprocal0.000.000.001.000.050.000.051.000.050.000.501.00KPOINTSC件结构:第1行：注释行。文案大全实用标准文档第2行：k点总数，或者0'(自动生成网格(Automatick-meshgeneration)。如果是前者，给出k点总数，又分为两种情况：M全手动(Enteringallk-pointsexplicitly),手动输入即自定义各个k点的坐标和权重，推荐只在能带计算时使用。此时：第3行：输入格式标识。直角坐标(Cartesian)或倒易坐标(Reciprocal)。同样的，c'、C、k、K'forCartesian,其他首字母则自动切换到Reciprocal。第4门行：逐个k点的坐标描述。格式为"xyzW其中"xyz”是该点的三个坐标，W是权重。所有k点的权重相互之间的比例对了就行，VASP会自动归一的。注意C坐标和R坐标的定义：C:k=(2Tt/a)(xyz)R:k=x*b1+y*b2+z*b3(b1-b3为倒易晶矢)这里“xyz”只是代表了坐标的顺序，与坐标轴无关。一般如非必要，可以先用自动模式生成k点，VASP会自动生成一个简约化后的k点矩阵，存于旧Z