各种计算中INCAR设置汇总

--------------------------------------------Band--------------------------------------------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改为线性模式以外，其他的设置与静态计算的设置一致。--------------------------------------------bader--------------------------------------------ISTART=1ICHARG=0LCHARG=.TURE.LAECHG=.TURE.NGX=132NGY=188NGZ=268~Bader计算读取优化后或静态（与优化一致）的WAVECAR（ISTART=1）,要保存Bader计算后的CHGCAR(LCHARG=.TRUE.),要设置LAECHG=.TRUE.，NGX/Y/Z的设置是优化计算中OUTCAR里面推荐的NGX/Y/Z值的2倍（比如，优化时OUTCAR中WARNING:aliasingerrorsmustbeexpectedsetNGXto54toavoidthem，此时bader计算中应设置NGX=108）.其他设置与优化时一致。-------------------------------------------GGA+U-------------------------------------------LMAXMIX=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.0---------------------------------------------HSEHF----------------------------------------------ISTART=1ICHARG=1GGA=PERPstandsforRPBE,PEforPBE,91forPW91LHFCALC=.TRUE.HartreeFockissetto；Theflagspecifies,whetherHartree-FocktypecalculationsareperformedHFSCREEN=0.207screeninglength(eitherq_TFor0.3A-1)ALGO=Damped;TIME=0.3timestepforELMPRECFOCK=FastNormal,FastorAccurate(LoworMediumforcompatibility)NKREDZ=3reducek-pointgridby--------------------------------------------OutputControl:--------------------------------------------LWAVE=.FALSE.LCHARG=.FALSE.LORBIT=11EMIN=-10.0EMAX=10.0NEDOS=2001HSE计算时由于计算量较大，一般只做单点校正，顺便会一起做DOS计算。做HSE计算时由于只做单点优化，因此最好读取PBE或RPBE优化后的WAVECAR/CHGCAR.具体设置只需在PBE优化中加入HF部分即可，HF的计算一般设置LHFCALC=.TRUE.,ALGO=Damped,TIME=0.3(一般推荐设为0.4，但不易收敛时可将其调小),具体设置详见VASP5.3.5手册119-121，123页，也可参见网页:--------------------------------------------ISTART=1ICHARG=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,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,我们需要的是DOSCAR-------------------------------------------PARCHG-------------------------------------------LPARD=.TURE.Evaluatepartial(bandand/ork-point)decomposedchargedensity.NBMOD=-3!0NumberofvaluesinthearrayIBAND.IfIBANDisspecified,NBMODissetautomaticallytothecorrectvalue(inthatcaseNBMODshouldnotbesetmanuallyintheINCARfile)0Takeallbandstocalculatethechargedensity,evenunoccupiedbandsaretakenintoaccount.-1Calculatethetotalchargedensityasusual.Thisisthedefaultvalueifnothingelseisgiven.-2CalculatethepartialchargedensityforelectronswiththereeigenvaluesintherangespecifiedbyEINT.-3Thesameasbefore,buttheenergyrangeisgivenvs.theFermienergyEINT=-1.2LSEPB=.FALSE.SpecifieswhetherthechargedensityiscalculatedforeverybandseparatelyandwrittentoafilePARCHG.nb.⋆(TRUE)orwhetherchargedensityismergedforallselectedbandsandwritetothefilePARCHG.ALLB.⋆orPARCHG.DefaultisFALSE.LSEPK=.FALSE.Specifieswhetherthechargedensityofeveryk-pointiswritetothefilesPARCHG.⋆.nk(TRUE)orwhetheritismerged(FALSE)toasinglefile.Ifthemergedfileiswritten,thentheweightofeachk-pointisdeterminedfromtheKPOINTSfile,otherwisethekpointsweightsofonearechosen.PARCHG分析主要是可以分析DOS图中具体的轨道，从而看其组成。PARCHG计算主要是基于DOS计算之后选取某一能区（能量区间）范围内的轨道进行分析，相当于是从DOSCAR中取出单独的信息，屏蔽不需要的轨道信息，因此不需要另外单独的计算，只是分解DOSCAR而已，孤儿其他的设置应该与DOS计算的设置一致，有时候KPOINTS的设置也可以小点，主要是精度会上不去，但都是想要单独分析某一能区的信息。ENIT设置区间。详见VASP5.3.5手册105页--------------------------------------------TSopt--------------------------------------------POTIM=0.05defaultis0.5forfrequencycalculationsasmallervaluesshouldbeset,0.02example.time-stepforionic-motionEDIFFG=-0.1Convergecriterion(positiveforenergy,negetiveforforce)NSW=300numberofstepsforIOMIBRION=2-1:noupdate0:MD1:quasi-New2:CGIMAGES=4time-stepforionic-motionSPRING=-5recommendedsetting;ThenudgedelasticbandmethodisappliedwhenSPRINGissettoanegativevalueLCLIMB=.FALSE.NEB方法计算过渡态，在00,05（XX）中间插入4个点（IMAGES=4，XX=numberofimages+1）,通过SPRING设置力常数(InadditiontotheIMAGEStag,aspringconstantcanbesuppliedintheSPRINGtag)，SPRING为负数时表示采用NEB方法计算力常数，一般默认为-5，LCLIMB=.FALSE.表示采用NEB方法计算，如果LCLIMB=.TRUE.则表示采用CINEB方法计算过渡态，一般最好不要用CINEB，有时候结果会比较诡异，一般难以收敛，除非过渡态很简单；过渡态优化时只能用POTIM=0.05，无论哪种计算方法，QN或CG;刚开始计算TS时可取精度低一点比如力收敛标准为0.1，此时可以做单点（能量最高的点）优化，或者继续NEB优化，刚开始切不可取过大的NSW，一般300步以内力不到0.1左右，则需改方法或改变初始猜测。也可参见网页=.TURE.LCHARG=.TURE.