《VASP参数设置详解.docx》由会员分享，可在线阅读，更多相关《VASP参数设置详解.docx（15页珍藏版）》请在三一办公上搜索。
1、VASP参数设置详解软件主要功能：采用周期性边界条件（或超原胞模型）处理原子、分子、团簇、纳米线（或管）、薄膜、晶体、准晶和无定性材料，以及外表体系和固体I计算材料的构造参数（键长、键角、晶格常数、原子位置等）和构型I计算材料的状态方程和力学性质体弹性模量和弹性常数）I计算材料的电子构造（能级、电荷密度分布、能带、电子态密度和ELF）I计算材料的光学性质I计算材料的磁学性质I计算材料的晶格动力学性质（声子谱等）I外表体系的模拟（重构、外表态和STM模拟）I从头分子动力学模拟I计算材料的激发态（GW准粒子修正）计算主要的四个参数文件：INCAR,poscar,potcar,kpoints,下面简
2、要介绍，详细权威的请参照手册INCAR文件：该文件控制VASP进展何种性质的计算，并设置了计算方法中一些重要的参数，这些参数主要包括以下几类： 对所计算的体系进展注释：SYSTEM 定义若何输入或构造初始的电荷密度和波函数：ISTART,ICHARG,INIWAV 定义电子的优化- 平面波切断动能和缀加电荷时的切断值：ENCUT,ENAUG- 电子局部优化的方法：ALGO,1ALG0,LDIAG一电荷密度混合的方法：IMIX,AMIX,AMIN,BMIX,AMIXMAG,BMIXMAG,WC,INlMIX,MIXPRE,MAXMlX- 自治迭代步数和收敛标准：NELM,NELMIN,NELMD
3、L,EDIFF定义离子或原子的优化- 原子位置优化的方法、移动的步长和步数：IBRION,NFREE,POTIM,NSW- 分子动力学相关参数:SMASS，TEBEG,TEEND,POMASS，NBLOCK,KBLOCK,PSTRESS- 离子弛豫收敛标准：EDIFFG定义态密度积分的方法和参数Tmearing方法和参数：ISMEAR,SIGMA-计算态密度时能量范围和点数：EMIN,EMAX,NEDOS-计算分波态密度的参数：RWIGS,LORBIT其它- 计算精度控制：PREC- 三性计算：1SPIN,MAGMOM,NUPD0WN- 交换关联函数：GGA,VOSKOWN- 计算ELF和总的
4、局域势：LELF,LVTOT- 构造优化参数：ISIF- 等等。主要参数说明如下：/SYSTEM：该输入文件所要执行的任务的名字。取值：字符串，缺省值：SYSTEM拈NWRITE：输出内容详细程度。取值：04,缺省值：2如果是做长时间动力学计算的话，最好选0或1（首末步/每步核运动输出），短时运算用2,选3那么会在出错的时候给出说明信息。nISTART：决定是否读取WAVECAR文件。取值：02,缺省0/1for无/有前次计算的WAVECAR（波函数）0：begin,fromscratch,根据INlWAV初始化波函数1：restartwithconstantenergycut-off,从WA
5、VECAR读取波函数（重定义平面波集）2：restartwithconstantbasisset,从WAVECAR读取波函数（平面波集不变）必ICHARG：决定若何建设初始电荷密度。取值：()2,缺省值:ifISTART=O2else。0：由初始波函数计算电荷密度1：从CHGCAR文件读取电荷密度2：使用原子电荷密度的叠加+10非自洽计算gISPIN:是否进展SPinPoIariZedCaICUIation。取值：1,2(l-no,2-yes),缺省值：2MAGM0M：在ICHARG=2或在CHGCAR中未包含磁化密度(ICHARG=I)时，指定每个原子的初始磁化时刻。取值：实数数组，缺省值:
6、对ISPIN=2NIONS*1.0,对非共线型磁化体系3*NIONS*LONINlWAV若何设置初始波函数，只在ISTART=O时使用。取值：0,110-最低动能的平面波，1-随机数)，缺省值：1。与IDlPOL控制计算单极、偶极和四极修正。取值：14。13只计算第一/二/三晶矢方向，适于厚板(SIab)的计算4所有方向都计算，适于计算孤立分子NPREC进动(precession)。取值：IowZmediumZhighZnormal/accurate/single),缺省值：Normal(VASP.4.X)；Medium(VASR5.X)VASP4.5+采用了优化的accurate来替代hig
7、h,所以一般不推荐使用high。不过high可以确保“绝对收敛，作为参考值有时也是必要的。同样受推荐的是normal,适于作为日常计算选项。受PREC影响的参数有四类：ENCUT；NGX,NGY,NGZ；NGXF,NGYF,NGZF；ROPTo如果设置了PREC,这些参数就都不需要出现了，当然直接设置相应的参数也有同样效果。具体影响效果见p5354NENeUT平面波基组的截断能量（eV）o取值：实数，缺省值：受PREC设置影响，从potcar文件中找出相应的Enmaxzenmin值来设置。PREC=LowMediumAccurateNormalENCUT=ENMINENMAXENMAXENMA
8、XSingleHighENMAXENMAX*1.3对于多个元素不同的ENMAXENMIN,都取最大值。该参数非常重要，最好不要手工去设置，除非文献告诉你要用多少，或者经过结果可靠性的验证。当然，为了测试一下提交的任务，也不妨先设个较小的值。NNGX,NGY,NGZ：控制FFT网格在三个晶矢方向上的格点数量。MNGFX,NGFY,NGFZ：控制第二次更准确的FFT网格的格点数量。也是两类重要的最好不要去动的参数。在未指定的情况下将根据PREC的设置从POTCAR中自动读取。PREC=HighZAccurate,基组中向量的2倍值，用来防止WraParounderrors,得到准确解。PREc=L
9、owZMediumZNormal,基组中向量的3/4倍值（已足够准确到1meV/atom）oBLREAL:决定投射是在实空间还是倒易空间进展。取值：.TRUE.（实空间）/.FALSE.（倒易空间），缺省值：.FALSE.用于求解惯势的非局域局部用到的一个积分，在倒格空间里采用平面波基组求解，在实空间里那么采用积分球求解。其他还有两个选项：OorOn,AorAutOoOn和.TRUE.的差异在于是否使用King-Smith算法优化，设为Auto那么进展自动选择，推荐使用。/ROPT:在LREAL=AUtOoron时，优化控制每个核周围的积分球内的格点数。取值：实数数组ForLREAL=OnPR
10、EC=Low,700pointsintherealspacesphere(ROPT=0.67)PREC=Med,I(X)Opointsintherealspacesphere(ROPT=LO)PREC=High,1500pointsintherealspacesphere(ROPT=1.5)ForLREAL=AutoPREC=Low,accuracyIO2(ROPT=O.Ol)PREC=Med,accuracy2*IO3(ROPT=0.002)PREC=Highaccuracy2*104(ROPT=2E-4)与NELM,NELMINandNELMDL：控制电子自治循环步数。取值：整数NELM：
11、电子自洽循环最大次数。缺省值：60NELMIN：电子自洽循环最小次数。缺省值：2NELMDL：弛豫次数。缺省值：ifISTART=O,INIWAV=LandIALGO=8,-5,ifISTART=O,INIWAV=1,andIALGO=48,-12,else0NELMDL可以取负值。如果初始波函数采用随机赋值，即ISTART=O,INIWAV=I,那么很可能开场的值比较离谱，那么在第一步核运动循环之前采用NELMDL(负值)步的非自洽(保存初始的H)步计算将减少计算所需的时间。NEDIFF:指定电子自洽循环的全局中断条件，用于控制收敛精度。取值：实数，缺省值：1(尸注意，即使EDlFF=0,N
12、ELM步电子自洽循环也会执行。幺EDIFFG：指定离子弛豫循环的中断条件，用于控制核运动的收敛精度。取值，实数,缺省值：10*EDIFFEDIFFG()在两个离子步的总自由能之差小于EDIFFG时停顿EDIFFGVO在所有的力都小于EDIFFG时停顿。EDlFFG=O在NSW步弛豫后停顿此参数不支持MD,仅用于弛豫。NSW：给出最大离子步数。取值：整数，缺省值：OoUNBLOCK,KBLOCK:取值：整数，缺省值：NBLOCK=I,KBLOCK=NSW在NBLOCK离子步后对成对相关函数和DOS进展计算，并且把离子配置写入XDATCAR文件。在KBLOCK*NBLOCK步主循环后平均的成对相关
13、函数和DOS被写入PCDAT和DOSCAR文件。MBRION：决定离子若何更新和运动。取值：-13,58(-l-无更新，()-MD,I-RMM-DIIS,2.共痈梯度算法，3-DampedMD,5,6：有限差分，7,8：密度函数扰动理论)，缺省值：ifNSW=O1,-1,elseO这个参数是和ISIF,IALGO/ALGO起决定若何算的最重要的参数。13是三种弛豫的方法，根据ISIF决定是否固定离子位置、晶胞大小和形状，在INCAR中必须设置参数POTIMoO是标准的ab-initioMD,不受ISlF影响，即不改变晶胞大小和形状。58支持HessianMatrix和phononfrequen
14、cy的计算以及局部固定的MDo幺POTlM：IBRloN=O时，给出MD每步步长(fs),IBRlON=13时，给出最小化的度量常量。取值：实数，缺省值：IBRION=O无缺省，必须指定，IBRION=I,2,30.5NlSlF：决定是否计算应力张量以及弛豫中晶胞变化的自由度。取值：()6,缺省值：ifIBRIOn=O(MD)0else2ISIFIcalculatecalculaterelaxchangechangeiforceIstresstensorionscellshapecellvolume0IyesnoIyesnono1IyestraceonlyIyesInono2IyesIyesI
15、yesnono3IyesIyesIyesIyesIyes4IyesIyesIyesIyesno5IyesIyesnoIyesno6IyesIyesnoIyesIyes7IyesIyesnonoIyestraceonlymeansthatonlythetotalpressureiscorrectNlWAVPR:决定波函数和/或电荷密度若何从一个离子配置向下一个离子配置进展推测。取值：03,KM3无推测，1,11-用原子电荷密度进展简单推测,2,12-二阶推测,3,13混合前两种方法)，缺省值：iflBRION=0(MD),l,2(relaxation)2else(静态计算)0推测结果保存在外部文
16、件TMPCAR中，取值+10那么全部使用内存，不保存此文件。NlSYM:决定是否使用对称性。取值：13(-l,0不使用，1,2,3-使用),缺省值：if使用US-PP1,if使用PAW2ISYM=2使用种效率更高也更节省内存的电荷密度对称性，1SYM=3时仅考虑力和应力张量的对称性，而电荷密度是非对称的。幺SYMPREC：决定POSCAR文件中给出的位置的精度。取值：实数，缺省值：104NLCORR:决定是否对非完全自洽计算中的力进展HarriS修正。取值：.TRUE.FALSE.,缺省值：.TRUE.NTEBEG,TEEND：控制从头分子动力学计算中的起始温度和最终温度(MD有效)。取值：实
17、数，缺省值：TEBEG=0TEEND=TEBEG注意VASP的温度定义与实际温度有细微的差异，所以TEBEG=TX(N-I)/N,T为实际要求的温度，N为原子数。SMASS：控制从头MD中的速度。取值：-3(),缺省值：-3微正那么系综(总自由能不变），2保持初速度不变，1每NBLOCK步调整速度，来保证动能连续，=0Nose算法模拟正那么系综WNPACO：成对相关函数的槽数。取值：整数，缺省值：256APACO：成对相关函数求值中的最大距离取值：整数，缺省值：16简单说就是在不超过APACO的NPACO个距离上求成对相关函数PCFt./RW1GS：给出Wigner-SeitZ半径，DOS计算
18、用。取值：实数数组，缺省值：从POTCAR文件中读取MNELECT:总电子数，如果系统不是电中性的就必须设置，所带电荷作为均一的背景电子气考虑。取值：实数，缺省值：-（价电子数），由POSCAR和Pe）TCAR文件自动决定（通常不必给出）。/NUPDOWN上下自旋成分间的电子数之差。取值：整数，缺省值：未设置（此时将进展完全弛豫）MEMIN,EMAX：DOS求值的最小/最大能量。取值：实数，缺省值：EMIN=-（lowestKS-eigenvalue-）,EMAX=（highestKS-eigenvalue-）NISMEAR：决定每个波函数的局部占位加若何设置。取值：-54-3-2IOlN（-
19、5-带有Blochl修正的四面体方法,-4.不带BloChl修正的四面体方法,3-根据INCAR文件中提供的smearing参数执行循环，-2-从WAVECAR文件中读取， 1-Fermi-Smearing,O-Gaussiansmearing,0-methodofMethfessel-PaxtonOrderM缺省值：1采用局部占位波函数，用一个函数来平滑积分，尤其是对于金属体系可减少k点。NSIGMA:决定Smearing的宽度（eV）。取值：实数，缺省值：0.2KALGO：指定电子最小化算法。取值：Normal（blockedDavidsonblockiterationSCheme）/Ve
20、ryFast（RMM-DHS）/FaSt（前两个算法的混合）A11（波函数的所有带同时更新）/Damped（dampedvelocityfrictionalgorithm）,缺省值:NormalIALGO:指定主算法（整数选择算法）。取值：8（共挽梯度算法）/38（DaVidSOnbIoCkiterationscheme）48（RMM-DIIS）,缺省值：38算法是最重要的参数之一。一般VASP推荐使用的是以上三种算法，一般来说8/38是初期比较快收敛，在接近平衡时采用48较快，在初期或MD时使用48可能会遇到不收敛的情况。也可以使用ALGO参数来替代IALG0,设置Fasi,VASP会先用3
21、8,再自动切换到48。各种算法只要收敛，结果应该一致。另一个可能有用的选项是-1。不进展实际的计算，只对重要的步骤做计算测试,并将测试得到的各局部耗时输出到OUTPUT中。幺VOSKOWN：决定是否使用VWN插值算法。取值：0（不使用）/1（使用），缺省值：0如果使用了PW91泛函或需要计算磁性质时可以设为1使用。必MiXing-tags： IMIX：混合的类型，取值：整数，缺省值：4 AMIX：线性混合参数。取值：实数，缺省值：0.8IUS-PP）,0.4（PAW） AMIN：最小混合参数。取值：实数，缺省值：0.1 BMIX：Kerker混合方案的截断波向量。取值：实数，缺省值：1.0 A
22、MIX_MAG：磁化过程的线性混合参数。取值：实数，缺省值：1.6 BMIX_MAG：磁化过程的Kerker混合方案的截断波向量。取值：实数，缺省值：1.0 WC：BrOyden混合方案中每步的加权因子。取值：实数，缺省值：1000.0 INIMIX：Broyden混合方案中的初始混合类型。取值：整数，缺省值：1 MIXPRE：Broyden混合方案中的预处理类型。取值：整数，缺省值：1 MAXMIX：BrOyden混合器中存储的最大步数。取值：整数，缺省值：-45值得注意的是，在MD或者弛豫的时候，设置MAXMIX（0,一般约3倍的电子SC步数）可能会大大减少核运动步数，但同时也会增加对内存
23、的要求。NLWAVE,LCHARG：决定是否把波函数（或电荷密度）写入外部文件WAVECAR（或CHGCAR和CHG）中。取值：.TRUE./.FALSE缺省值：.TRUE.LVTOT：决定是否把总局域势写入外部文件LoCPoT中。取值：.TRUE./.FALSE.,缺省值：.FALSE.NLELF：决定是否创立ELFCAR文件。取值：.TRUE./.FALSE缺省值：.FALSE.ELFCAR用于保存ELF（electronlocalizationfunction）LORBIT:和适当的RWIGS一起决定是否创立PROCAR或PROOUT文件。取值：0|1251011I12,缺省值：0（创立
24、DOSCAR和PROCAR文件）/NPAR：用于控制VASP切换到带间并行模式。取值：整数，缺省值：总节点数。每一个节点计算一个带，可以提高并行效率，减少通讯量，但可能会大幅增加内存的需求。/NBANDS：给出计算中实际的总能带数。取值：整数，缺省值：NELECT/2+NIONS/2(非自旋)，0.6*NELECT+NMAG(自旋)NBANDS的取值应使计算中包含相当数量的空带。因为计算需要大量的空带，至少要求1个空带(否那么VASP会给出警告NBANDS对于解决内存需求非常重要。一般NBANDS在NELECT/2+NIONS/2以上可以得到较为准确的结果，但如果内存不够就只能减少NBANDS
25、,在牺牲精度和体系大小之间平衡了。POTCAR文件鹰势文件，最重要的输入文件之一。可以理解为分子力学模拟中的力场文件，但包含的信息更多。VASP将各元素优化的INCAR里的参数也包含在该文件中，作为支持PREC的缺省选择。通常各元素的POTCAR文件已经包含在VASP软件包中的度势库里了，用户需要做的是确定自己具体需要哪几种腰势，然后按照POSCAR文件里的顺序，将所选择的POTCAR文件按顺序连接起来就可以了。如以下命令：catfilelfile2file3POTCAR燃势库中的鹰势文件可以分为以下几类：根据产生方法的不同有Ultra-SOfl底势(USPP)和投影扩大波震势(PAW)； 根
26、据交换关联函数的不同有LDA(localdensityapproximation)和GGA(generalizedgradientapproximation又可以再分为PW91和PBE)； 根据半芯态的不同元素A的鹰势文件还可以分为A,A_sv和A_pv： 根据选取的不同截断能量(ENMAX)元素A的鹰势文件还可以分为A,A_s和A_h。VASP的鹰势文件放在目录vaspPOlentiaIS下，该目录又包含五个子目录：pot,pot_GGA,potpaw,potpaw_GGA,POtPaW_PBE,其中每个子目录对应一种势形式：pot=PRLDA：pot.GGA=PP,GGA；potpaw=P
27、AW,LDA；potpaw_GGA=PAW,GGA,PW91；potpaw_PBE=PAW,GGA,PBEo每个目录中每种元素根据截断能量和半芯态的不同还会有多个对应的鹰势文件存在。在具体选取时可以参考各版本同目录下的V_RHFIN和PSCTR文件,这两个文件说明了该版本的腰势是若何生成的。选好所要使用的鹰势之后，进入对应的目录，会看到里面有四个文件：POTCAR.Z,PSCTR.Z,V_RHFIN.Z和WS_FTP.LOG。现在需要用到的是第一个，把各元素的该文件解压(zcatPOTCAR.Zfile),然后用CP或mv命令把这些文件移到工作目录里，再用Cat命令把它们合并至UPoTCAR中
28、，就得到了我们需要的PoTCAR。注意要记住这里元素的排列顺序，以后在POSCAR文件中各元素的排列就是按照这个顺序。POSCAR文件位置文件。描述所计算体系的晶胞参数、原子个数及晶胞中原子的位置，以及分子动力学计算时原子的初始速度。POSCAR文件例如：Si-fcc5.430.00.500.500.500.000.50().5()().50().(M)2Direct0.00.000.000.250.250.25POSCAR文件的构造：第1行：任意文字注释。第2行：晶格常数(?),也是晶矢的缩放系数，后面所有长度值得自原值除以此值。a=b=c时取a即可，否那么一般取三者最大值，假设取负值，那么
29、为晶胞体积(?3)。第35行:定义晶矢。第6行：每种元素的原子个数，特别注意顺序,要与下面的坐标顺序以及Pc)TCAR中的顺序一致。第7行：可省略，无需空行。做动力学时，用于指定是否需要固定局部原子的坐标。假设是，此行以S或者s作为首字母即可。第8行开场为原子的坐标位置，格式为optionlinecoordinate!ofdementicoordinate?ofdementiCoordinaleNofelement1optionlinecoordinate!ofelement2coordinate?ofelement2CoordinaleMofelement2其中，optionline指定输入
30、坐标的格式，除了第一个以外，如果后面的输入格式同前，那么都可以无空行省略。optionIine可指定的输入坐标格式有两种:DordfordirectmodeCorcorKorkforcartesianmode顾名思义，前者是定义在三个晶矢方向上的坐标：R=Rlx+R2y+R3z,R1,R2,R3为前面的晶矢，x,y,z为输入的三个坐标，R为原子坐标位矢，是相对于晶矢给出的，而后者那么是以笛卡尔坐标系来给出原子的绝对坐标，实际就是将直角坐标除以前面第二行定义的晶胞常数后得到的坐标值。如果第7行设定为S(SelectiveDynamic),那么可以用以下形式定义各坐标是否可以移动：Selectiv
31、edynamicsCartesian().(X)().(X)().()OTTF0.250.250.25FFFKPOINTS文件设置布里渊区k点网格取样大小或能带构造计算时沿高对称方向的k点。一般来说，k点越密越多，计算精度也就越高，但计算成本也会越高。因此，对于原子数较多的体系的计算，需要慎重的尝试k点数目，在防止或者预先评估wrap-arounderror的前提下尽量减少k点数目。KPOINTS文件例如：k-pointsalonghighsymmetrylines11Reciprocal0.000.000.001.000.05().00().051.000.050.00().5()1.00K
32、POINTS文件构造：第行：注释行。第2行：k点总数，或者0(自动生成网格(Automatick-meshgeneration)(如果是前者，给出k点总数，又分为两种情况：M：全手动(Enteringallk-pointsexplicitly),手动输入即自定义各个k点的坐标和权重，推荐只在能带计算时使用。此时：第3行：输入格式标识。直角坐标(Cartesian)或倒易坐标(Reciprocal)o同样的，c、C、k、*K,forCartesian,其他首字母那么自动切换到ReCiPrOCa1。第4n行:逐个k点的坐标描述。格式为“xyzW。其中“xyz”是该点的三个坐标，W是权重。所有k点的
33、权重相互之间的比例对了就行，VASP会自动归一的。注意C坐标和R坐标的定义:C：k=(2a)(xyz)R：k=x*bl+y*b2+z*b3(bl-b3为倒易晶矢)这里“xyz”只是代表了坐标的顺序，与坐标轴无关。一般如非必要，可以先用自动模式生成k点，VASP会自动生成一个简约化后的k点矩阵，存于IBZKPT文件，可以直接复制里面的数据到KPoINTS文件中使用，这也是该输入法的主要用途，可以减少重复自动生成格点的时间。另一个用途是为了做准确的DoS(Densityofstatus)计算，由于这类计算所需的k点数极大，通过全手动尽可能的优化k点也就必需了。1.：半手动/线形模式(Line-mo
34、de)：在计算能带构造时使用，此时需要准确地选取k点，在指定的高对称性方向上生成指定数目的k点。此时：第2行：指定两点间生成的k点数，不同于全自动的总k点数。第2.5行：L表示是线形模式。第3行：输入格式标识，同前。第4n行：每行描述一个点，格式为“xyz。每两行的点连成一条直线，在两点间生成指定数目的k点。每隔两行间以空行区分，例如：k-pointsalonghighsymmetrylines101.ine-modeReciprocalOOO0.50.500.50.50().50.750.25最常用的方法是让程序自动生成k点网格，用户只需定义网格取样大小即可。此时：第2行：0第3行：生成k点
35、的方法。“A：全自动模式，可以看作以点为圆心以1为半径做圆，当然各晶矢不同时，相应的圆就被拉成了椭圆，以确保三个倒格晶矢方向上覆盖的k点数都为1“G：17Gamma方法M：MOnkhorSI-PaCk方法其它首字母：高级模式，就是用C坐标或R坐标直接输入新的晶矢，一般不使用第4行：在各个晶矢方向上分割各晶矢的点数。一般为奇数，使得产生的k点是以Gamma点为中心的。第5行：是否移动网格点以及移动多少，也即偏移原点的位矢，一般设成“000”,表示不移动。例如：Automaticgeneration0Monhkorst-Pack9990.00.0().0VASP的输出文件主要包括：OUTCAR：包
36、含运行VASP作业后得到的绝大局部计算结果以及每步迭代的详细情况;CHG,CHGCAR：给出体系的电荷密度。两个文件内容一样，但前者的数据精度比后者略低一些；DOSCAR,EIGENVAL：分别给出所计算体系的电子态密度和本征值。两个文件中的能量值都是绝对的；WAVECAR：给出所计算体系的电子波函数，二进制文件，不可编辑；OSZICAR：包含每次迭代或离子移动情况的简单汇总；C0NTCAR：给出离子进展弛豫时，每次移动后体系的晶格参数，与POSCAR内容一样；PCDAT,XDATCAR：给出有关分子动力学模拟中的一些结果，如配对相关函数；4推荐以I分享到：阅读(197)评论(0)|引用(0)|举报
