各种计算中INCAR设置汇总

1、ICHARG=11ENCUT=600PREC=Accuracyprecision:Medium/High/Lowdefaultmedium.thesametime,ROPTaccuracyofapproximately1meV/atom(10e-3).recommendtousePRECNormalforeverydaycalculationsinVASP.4.5LWAVE=.FALSE.LCHARG=.FALSE.Band计算要读取静态(精度提高后比如提高PREC,ENCUT,)的CHGCAR,除了ICHARG=11以及KPOINTS改为线性模式以外，其他的设置与静态计算的设置一致。bade

2、rISTART=1ICHARG=0LCHARG=.TURE.LAECHG=.TURE.NGX=132NGY=188NGZ=268Bader计算读取优化后或静态(与优化一致)的WAVECAR(ISTART=1)，要保存Bader计算后的CHGCAR(LCHARG=.TRUE.),要设置LAECHG=.TRUE.,NGX/Y/Z的设置是优化计算中OUTCAR里面推荐的NGX/Y/Z值的2倍(比如，优化时OUTCAR中WARNING:aliasingerrorsmustbeexpectedsetNGXto54toavoidthem,此时bader计算中应设置NGX=108).其他设置与优化时一致。G

3、GA+ULMAXMIX=4!fordelementsincreaseLMAXMIXto4,f:LMAXMIX=6LDAU=.TURE.LDAUTYPE=2LDAUU=300LDAUL=2-1-1GGA+U设置，可从OUTCAR中看出，LDA+Uisselected,typeissettoLDAUTYPE=2angularmomentumforeachspeciesLDAUL=2-1-1U(eV)foreachspeciesLDAUU=3.00.00.0J(eV)foreachspeciesLDAUJ=0.00.00.0HSEHFISTART=1ICHARG=1GGA=PERPstandsfor

4、RPBE,PEforPBE,91forPW91LHFCALC=.TRUE.HartreeFockissetto；Theflagspecifies,whetherHartree-FocktypecalculationsareperformedHFSCREEN=0.207screeninglength(eitherq_TFor0.3A-1)ALGO=Damped;TIME=0.3timestepforELMPRECFOCK=FastNormal,FastorAccurate(LoworMediumforcompatibility)NKREDZ=3reducek-pointgridbyOutputC

5、ontrol:LWAVE=.FALSE.LCHARG=.FALSE.LORBIT=11EMIN=-10.0EMAX=10.0NEDOS=2001HSE计算时由于计算量较大，一般只做单点校正，顺便会一起做DOS计算。做HSE计算时由于只做单点优化，因此最好读取PBE或RPBE优化后的WAVECAWCHGCAR.具体设置只需在PBE优化中加入HF部分即可，HF的计算一般设置LHFCALC=.TRUE.,ALGO=Damped,TIME=0.3（一般推荐设为0.4,但不易收敛时可将其调小），具体设置详见VASP5.3.5手册119-121,123页，也可参见网页Related:ISTART=1ICH

6、ARG=0PREC=Accurateprecision:Medium/High/Lowdefaultmedium.thesametime,ROPTaccuracyofapproximatelyISMEAR=0SIGMA=0.02LORBIT=11NEDOS=2101EMAX=6EMIN=-15DOS计算要读取静态或优化后的WAVECAR,然后提高精度ENCUT,PREC=Accurate,KPOINTS(abc)设置比优化的要高即a+2,b+2,c+2,加2是为了使取点为基数，也可以+4/6/8等，但计算量会很大，有时候服务器承受不了，计算不动，因此不宜取太大；或者读取提高精度后(ENCUT,

7、PREC,KPOINTS)静态的CHARG,然后将SIGMA设小一点，之后就是能窗的设置，一般是优化后的E-Fermi值加减10就可以，比如上面所设的值，优化后的fermi能级应该在4-5之间，然后EMAX设为6,EMIN设为-15,此时NEDOS=(EMAX-EMIN)*100+1,力口1是由于加上本身，比如0-100之间总共有100+1个个数。其他的设置与一般优化设置相同，最后不必保留计算后的WAVECAR,CHGCAR,我们需要的是DOSCARPARCHGLPARD=.TURE.Evaluatepartial(bandand/ork-point)decomposedchargedensi

8、ty.NBMOD=-3!>0NumberofvaluesinthearrayIBAND.IfIBANDisspecified,NBMODissetautomaticallytothecorrectvalue(inthatcaseNBMODshouldnotbesetmanuallyintheINCARfile)0Takeallbandstocalculatethechargedensity,evenunoccupiedbandsaretakenintoaccount.- 1Calculatethetotalchargedensityasusual.Thisisthedefaultvalu

9、eifnothingelseisgiven.- 2CalculatethepartialchargedensityforelectronswiththereeigenvaluesintherangespecifiedbyEINT.- 3Thesameasbefore,buttheenergyrangeisgivenvs.theFermienergyEINT=-1.2LSEPB=.FALSE.SpecifieswhetherthechargedensityiscalculatedforeverybandseparatelyandwrittentoafilePARCHG.nb.?(TRUE)orw

10、hetherchargedensityismergedforallselectedbandsandwritetothefilePARCHG.ALLB.?orPARCHG.DefaultisFALSE.LSEPK=.FALSE.Specifieswhetherthechargedensityofeveryk-pointiswritetothefilesPARCHG.?.nk(TRUE)orwhetheritismerged(FALSE)toasinglefile.Ifthemergedfileiswritten,thentheweightofeachk-pointisdeterminedfrom

11、theKPOINTSfile,otherwisethekpointsweightsofonearechosen.PARCHG分析主要是可以分析DOS图中具体的轨道，从而看其组成。PARCHG计算主要是基于DOS计算之后选取某一能区(能量区间)范围内的轨道进行分析，相当于是从DOSCAR中取出单独的信息，屏蔽不需要的轨道信息，因此不需要另外单独的计算，只是分解DOSCAR而已，孤儿其他的设置应该与DOS计算的设置一致，有时候KPOINTS的设置也可以小点，主要是精度会上不去，但都是想要单独分析某一能区的信息。ENIT设置区间。详见VASP5.3.5手册105页TSoptPOTIM=0.05def

12、aultis0.5forfrequencycalculationsasmallervaluesshouldbeset,0.02example.time-stepforionic-motionEDIFFG=-0.1Convergecriterion(positiveforenergy,negetiveforforce)NSW=300numberofstepsforIOMIBRION=2-1:noupdate0:MD1:quasi-New2:CGIMAGES=4time-stepforionic-motionSPRING=-5recommendedsetting;Thenudgedelasticb

13、andmethodisappliedwhenSPRINGissettoanegativevalueLCLIMB=.FALSE.NEB方法计算过渡态，在00,05(XX)中间插入4个点(IMAGES=4,XX=numberofimages+1),通过SPRING设置力常数(InadditiontotheIMAGEStag,aspringconstantcanbesuppliedintheSPRINGtag),SPRING为负数时表示采用NEB方法计算力常数，一般默认为-5,LCLIMB=.FALSE.表示采用NEB方法计算，如果LCLIMB=.TRUE.则表示采用CINEB方法计算过渡态，一般最

14、好不要用CINEB,有时候结果会比较诡异，一般难以收敛，除非过渡态很简单；过渡态优化时只能用POTIM=0.05，无论哪种计算方法，QN或CG;刚开始计算TS时可取精度低一点比如力收敛标准为0.1,此时可以做单点(能量最高的点)优化，或者继续NEB优化，刚开始切不可取过大的NSW,一般300步以内力不到0.1左右，则需改方法或改变初始猜测。也可参见网页/vtsttools/optmizationforfreqLWAVE=.TURE.LCHARG=.TURE.POTIM=0.05set,0.02example.IBRION=1defaultis

15、0.5forfrequencycalculationsasmallervaluesshouldbe-1:noupdate0:MD1:quasi-New2:CGEDIFFG=-0.03NSW=1000Convergecriterion(positiveforenergy,negetiveforforce)numberofstepsforIOMNEB插入点做过渡态后取能量最高的点然后只能用QN的方法优化得到过渡态，此时POTIM可取范围0.050.2,最好刚开始用0.05,当力接近0.1但一直在0.1附近徘徊时，可以重新建POTIM取大一点比如0.1或0.2,之后再继续优化，但当力趋近0.05又一

16、直不能更小时，这时可以将能量的收敛精度调小比如EDIFF=1e-5,然后继续优化，一般过渡态收敛到力<0.05,能量标准1e-5就可以啦(有时候过渡态比较难收敛时做到0.05就可以啦)。最好保留WAVECAR或CHARG,下一步频率计算需要用到。ISTART=1ICHARG=1POTIM=0.02defaultis0.5forfrequencycalculationsasmallervaluesshouldbeset,0.02example.NSW=1numberofstepsforIOMIBRION=5-1:noupdate0:MD1:quasi-New2:CG#NPAR=4NFREE

17、=2频率计算时要读上一步优化的WAVECAR或CHGCAR,或2者同时读取也可以，有时候也可以不用读，但一般是磁性体系或是自选变化比较太的体系其态不易收敛的必须要读，主要是设置POTIM=0.02或0.01,NFREE=2或4,NSW=1,IBRION=5,这4个参数是计算频率的必设，其中NFREE=2表示是在我们所取点的左右2边各取一个点来拟合曲线做其二阶导求力常数k,继而求得频率，详见VASP5.3.5手册60-61页(NFREE=2usescentraldifference,i.e.eachionisdisplacedineachdirectionbyasmallpositiveandn

18、egativedisplacement,±POTIMx?x,±POTIMx?y,土POTIM；xForzNFREE=4,fourdisplacementareused,±POTIMx?xand±2POTIMx?x±POTIMx?yand±2POTIM不FoyNFREE=1,onlyasingledisplacementisapplied(itisstronglyrecommendtoavoidNFREE=1).)；ForIBRION=5,allatomsaredisplacedinallthreeCartesiandirections

19、,resultinginasignificantcomputationaleffortevenformoderatelysizedhighsymmetrysystems.ForIBRION=6,howeveronlysymmetryinequivalentdisplacementsareconsidered,andtheremainderoftheHessianmatrixisf川edusingsymmetryconsiderations.SelectivedynamicsarepresentlyonlysupportedforIBRION=5;TheparameterNFREEdetermi

20、neshowmanydisplacementsareusedforeachdirectionandion,andPOTIMdeterminesthestepsize.Thestepsizeisdefaultedto0.015?A,我们推荐的使用0.01或0.02其中对于我们的服务器cluster,不知道出于什么原因要注释掉NPAR,不然设为1或4都会报错，然后死掉。IonicRelaxation:EDIFFG=-0.03Convergecriterion(positiveforenergy,negetiveforforce)POTIM=2defaultis0.5forfrequencycalculationsasmallervaluesshouldbeset,0.02example.timestepinfsNSW=2000IBRION=0-1:noupdate0:MD1:quasi-New2:CGSMASS=3controlsthevelocitiesduringanab-initiomoleculardynamics.TE