初学VASP最重要的INCAR参数资料下载.pdf
- 文档编号:5972461
- 上传时间:2023-05-05
- 格式:PDF
- 页数:6
- 大小:138.71KB
初学VASP最重要的INCAR参数资料下载.pdf
《初学VASP最重要的INCAR参数资料下载.pdf》由会员分享,可在线阅读,更多相关《初学VASP最重要的INCAR参数资料下载.pdf(6页珍藏版)》请在冰点文库上搜索。
ifISTART=02else0ISPINspinpolarizedcalculation(2-yes1-no)default2MAGMOMinitialmagmoment/atomDefaultNIONS*1INIWAVinitialelectrwf.:
0-lowe1-randDefault1onlyusedforstartjobs(ISTART=0)IDIPOLcalculatemonopole/dipoleandquadrupolecorrections1-3只计算第一/二/三晶矢方向适于slab的计算4全部计算尤其适于就算孤立分子PRECprecession:
medium,highorlow(VASP.4.5+also:
normal,accurate)Default:
MediumVASP4.5+采用了优化的accurate来替代high,所以一般不推荐使用high。
不过high可以确保绝对收敛,作为参考值有时也是必要的。
同样受推荐的是normal,作为日常计算选项,可惜的是说明文档提供的信息不足。
受PREC影响的参数有四类:
ENCUT;
NGX,NGY,NGZ;
NGXF,NGYF,NGZF;
ROPT如果设置了PREC,这些参数就都不需要出现了当然直接设置相应的参数也是同样效果的,这里不展开了,随后详释ENCUTenergycutoffineV:
defaulttakenfromPOTCAR-fileimportant!
重要到几乎最好不要手工去设置除非文献告诉你要用多少,或者经过结果可靠性的验证当然,为了测试一下提交的任务,也不妨先设个较小的值附加说明:
当且仅当POTCAR里头没有设置ENCUT时(其实貌似没有才是常态),才受PREC设置影响从POTCAR里找出相应的ENMAX/ENMIN值来设置。
PREC=LowMediumAccurateHighENCUT=ENMINENMAXENMAX130%ENMAX对于多个元素的POTCAR不同的ENMAX/ENMIN,都取最大值NGX,NGY,NGZ:
FFTmeshforwavefunctionsNGFX,NGFY,NGFZ:
FFTmeshforcharges也是两类重要的最好不要去动的参数,PREC设置将从POTCAR中自动读取。
PREC=High,Accurate2倍值,用来避免wraparounderrors得到精确解PREC=Low,Medium,Normal3/4也已经足够精确到1meV/atomLREAL:
Default=.FALSE.赝势的非局域部分用到的一个积分在倒格空间或者实空间都可以求值。
这个选项就是决定是在哪个空间里求。
在倒格空间里,采用平面波基组求解,在实空间里,采用积分球求解。
缺省是.FALSE,即不在实空间求。
但效率会低一些。
其他选项是OorOn,AorAuto和.True.。
On和.TRUE.的差别在于是否使用King-Smith算法优化,Auto则自动选择,推荐。
ROPT:
优化控制每个核周围的积分球内的格点数,LREAL=AutoorOnForLREAL=OnPREC=Low700pointsintherealspacesphere(ROPT=0.67)PREC=Med1000pointsintherealspacesphere(ROPT=1.0)PREC=High1500pointsintherealspacesphere(ROPT=1.5)ForLREAL=AutoPREC=Lowaccuracy1e-2(ROPT=0.01)PREC=Medaccuracy2e-3(ROPT=0.002)PREC=Highaccuracy2e-4(ROPT=2E-4)NELM,NELMINandNELMDLnr.ofelectronicstepsDefault最大电子自洽循环次数NELM=60最小次数NELMIN=2弛豫次数NELMDL=-5ifISTART=0,INIWAV=1,andIALGO=8-12ifISTART=0,INIWAV=1,andIALGO=480else如果初始的波函数采取随机赋值,即ISTART=0,INIWAV=1,那么很可能开始的值比较离谱,那么在第一步核运动循环之前采用NELMDL(负值)步的非自洽(保留初始的H)步计算将减少计算所需的时间。
如果NELMDL取正值,将在每次核运动之后附加指定次数的弛豫步,目前不知道可以干嘛EDIFF电子SC循环的收敛精度缺省:
1e-4注意,即使EDIFF=0,NELM步也会执行EDIFFG核运动的收敛精度缺省:
EDIFF*10(总能量)EDIFFGNSW指定核运动步数缺省:
0NBLOCKandKBLOCKinnerblock;
outerblockDefaultNBLOCK=1KBLOCK=NSWIBRIONionicrelaxation:
-1-Fixed0-MD1-quasi-New2-CG3-Damp5-freqDefaultifNSW=0or1IBRION=-1elseIBRION=0这个参数是和ISIF;
IALGO/ALGO一起决定怎么算的最重要的参数1-3是三种Relax的方法,受ISIF决定是否固定核位置、晶胞大小和形状0是标准的ab-initioMD,不受ISIF影响,即不改变晶胞大小和形状5大概是和0差不多吧?
支持Hessian和Freq(仅点)的计算以及部分固定的MD详细的要在示例中具体情况具体分析了。
ISIFcalculatestressandwhattorelaxDefaultifIBRION=0(MD)0else2ISIFcalculatecalculaterelaxchangechangeforcestresstensorionscellshapecellvolume0yesnoyesnono1yestraceonlyyesnono2yesyesyesnono3yesyesyesyesyes4yesyesyesyesno5yesyesnoyesno6yesyesnoyesyes7yesyesnonoyesTraceonlymeansthatonlythetotalpressureIWAVPRpredictionofwf.:
0-non1-charg2-wave3-combDefaultifIBRION=0(MD)2ifIBRION=1,2(relaxation)1else(staticcalculation)0以上选项保存TMPCAR+10则全部使用内存,不保存此文件IWAVPRdetermineshowwavefunctionsand/orchargedensityareextrapolatedfromoneionicconfigurationtothenextconfiguration.ISYMsymmetry:
0-nonsym1-usesym是否使用对称性Default1SYMPRECdeterminesprecisionofthepositionsinPOSCARfile.Default1e-5LCORRHarris-correctiontoforces.Default.TRUE.POTIMtime-stepforion-motion(fs)DefaultIBRION=0(MD)nodefault,必须指定,MD每步步长IBRION=1,2,3(relaxation)0.5最小化的scalingconstant,尤其是IBRION=1TEBEG,TEENDtemperatureduringrun(MD有效)Default:
TEBEG=0TEEND=TEBEG注意VASP的温度定义与实际温度有细微的差别,所以TEBEG=T(N-1)/NT为实际温度,N为原子数SMASS控制MD中的速度模拟方法default-3微正则系综(总自由能不变)-2保持初速度不变-1每NBLOCK步调整速度,来保证动能连续=0Nos算法模拟正则系综,(不懂-,-)NPACOandAPACONPACO:
numberofslotsforpaircorrelation(PC)function.Default256APACO:
maximumdistancefortheevaluationofPCfunctioninA.Default16简单说就是在不超过APACO的NPACO个距离上求成对相关函PCFRWIGSWigner-Seitz半径DOS计算用NELECT总电子数如果系统不是电中性的就必须设置,所带电荷作为均一的背景电子气考虑NUPDOWNdefault不考虑电子自旋态改变的可能EMIN,EMAXenergy-rangeforDOSCARfileISMEARpart.occupancies:
-5tetwithBlochl-4-tet-1-fermi0-gaus0MP采用所谓部分占有波函数,用一个函数来平滑积分,尤其是对于金属体系可减少k点DefaultISMEAR=1如果在KPOINTS里使用了tetrahedra方法推荐ISMEAR=5SIGMAdeterminesthewidthofthesmearingineVDefaultSIGMA=0.2ALGOalgorithm:
Normal(Davidson)|Fast(mixed)|Very_Fast(RMM-DIIS)IALGOalgorithm:
useonly48(RMM-DIIS)or38(Davidson)or8(CG)DefaultIALGO=38forVASP4.5算法是最重要的参数之一。
一般VASP推荐使用的是以上三种算法,一般来说8/38是初期比较快收敛,在接近平衡时采用48较快,在初期或MD时使用48可能会遇到不收敛的情况。
也可以使用ALGO参数来替代IALGO,设置Fast,VASP会先用38,再自动切换到48。
各种算法只要收敛,结果应该一致。
另一个可能有用的选项是-1。
不进行实际的计算,只对重要的步骤做计算测试,并将测试得到的各部分耗时输出在OUTPUT里。
VOSKOWNuseVWNinterpolation算法,default0不用,如果使用了PW91或者需要计算磁性质,不妨设为1用mixingtag&
MAXMIXIMIX=typeofmixingAMIX=linearmixingparameterAMIN=minimalmixingparameterBMIX=cutoffwavevectorforKerkermixingschemeAMIX_MAG=linearmixingparameterformagnetizationBMIX_MAG=cutoffwavevectorforKerkermixingschemeformag.WC=weightfactorforeachstepinBroydenmixingschemeINIMIX=typeofinitialmixinginBroydenmixingschemeMIXPRE=typeofpreconditioninginBroydenmixingschemeMAXMIX=maximumnumberstepsstoredinBroydenmixerDefaultforUS-PPnon-magneticIMIX=4AMIX=0.8BMIX=1.0WC=1000.INIMIX=1MIXPRE=1MAXMIX=-45值得注意的是,在MD或者Relax的时候,设置MAXMIX(0,一般约3倍的电子SC步数),可能会大大减少核运动步数。
但是同时会增加对内存的要求。
LWAVE,LCHARGandLVTOTcreateWAVECAR/CHGCAR/LOCPOTLELFcreateELFCARLORBITcreatePROOUT输出文件的选项NPAR并行计算band的节点数,每一个节点计算一个band当然可以提高并行效率,减少通讯量,不过貌似现在硬件的主要限制还是内存,而这个选项的使用可能会大幅增加内存的需求NBANDS总能带数。
之所以把它放在最后,是因为它对于解决内存需求的重要性。
计算需要大量的能带(空带),至少要1个空带(否则VASP会给出警告)。
一般NBANDS=NELECT/2+NIONS/2以上可得到较精确的结果,如果内存不够就只好减少NBANDS,在牺牲精度和体系大小之间平衡了。
最后提示一下大多数参数的首字母代表了参数的性质I初始化L逻辑开关E能量N数目T温度便于记忆
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 初学 VASP 最重要 INCAR 参数