INCAR各个参数是什么意思.docx

INCAR各个参数是什么意思.docx

《INCAR各个参数是什么意思.docx》由会员分享，可在线阅读，更多相关《INCAR各个参数是什么意思.docx（10页珍藏版）》请在冰豆网上搜索。

INCAR各个参数是什么意思

INCAR是决定howtodo的文件

  限于能力，只对部份最大体的一些参数（>,没有那个标志的参数都是能够不显现的）

  详细说明，在那个地址只是简单介绍这些参数的设置，详细的问题在后文具体例如中展开。

  部份可能会干扰VASP运行的参数在那个地址被刻意隐去了，需要的同窗仍是请查看VASP自带

  的帮忙文档原文。

  参数列表如下：

  >SYSTEMnameofSystem

  任务的名字***

  >NWRITEverbositywrite-flag（howmuchiswritten）

  输出内容详细程度0-3缺省2

  若是是做长时刻动力学计算的话最好选0或1（首末步/每步核运动输出）

  听说也能够结合shell的tail或grep命令手动输出

  >ISTARTstartjob:

  restart选项0-3缺省0/1for无/有前次计算的WAVECAR（波函数）

  1'restartwithconstantenergycut-off'

  2'restartwithconstantbasisset'

  3'fullrestartincludingwavefunctionandchargeprediction'

  ICHARGcharge:

1-file2-atom10-constDefault:

ifISTART=02else0

  ISPINspinpolarizedcalculation（2-yes1-no）default2

  MAGMOMinitialmagmoment/atomDefaultNIONS*1

  INIWAVinitialelectrwf.:

0-lowe1-rand

  Default1onlyusedforstartjobs（ISTART=0）

  IDIPOLcalculatemonopole/dipoleandquadrupolecorrections

  1-3只计算第一/二/三晶矢方向适于slab的计算

  4全数计算尤其适于就算孤立分子

  >PRECprecession:

medium,highorlowalso:

normal,accurate）

  Default:

Medium+采纳了优化的accurate来替代high,因此一样不推荐利用

  high。

只是high能够确保'绝对收敛'，作为参考值有时也是必要的。

  一样受推荐的是normal，作为日常计算选项，可惜的是说明文档提供的信息不足。

  受PREC阻碍的参数有四类：

ENCUT;NGX,NGY,NGZ;NGXF,NGYF,NGZF;ROPT

  若是设置了PREC，这些参数就都不需要显现了

  固然直接设置相应的参数也是一样成效的，那个地址不展开了，随后详释

  >ENCUTenergycutoffineV:

defaulttakenfromPOTCAR-file

  important!

重要到几乎最好不要手工去设置

  除非文献告知你要用多少，或通过结果靠得住性的验证

  固然，为了测试一下提交的任务，也不妨先设个较小的值

  附加说明：

  当且仅当POTCAR里头没有设置ENCUT时（其实貌似没有才是常态），才受PREC设置影

  响从POTCAR里找出相应的ENMAX/ENMIN值来设置。

  PREC=LowMediumAccurateHigh

  ENCUT=ENMINENMAXENMAX130%ENMAX

  关于多个元素的POTCAR不同的ENMAX/ENMIN，都取最大值

  >NGX,NGY,NGZ:

FFTmeshforwavefunctions

  >NGFX,NGFY,NGFZ:

FFTmeshforcharges

  也是两类重要的最好不要去动的参数，PREC设置将从POTCAR中自动读取。

  PREC=High,Accurate2倍值，用来幸免wraparounderrors取得精准解

  PREC=Low,Medium,Normal3/4也已经足够精准到1meV/atom

  >LREAL:

Default=.FALSE.

  赝势的非局域部份用到的一个积分在倒格空间或实空间都能够求值。

那个选项就

  是决定是在哪个空间里求。

在倒格空间里，采纳平面波基组求解，在实空间里，采纳积

  分球求解。

  缺省是.FALSE，即不在实空间求。

但效率会低一些。

  其他选项是OorOn,AorAuto和.True.。

  On和.TRUE.的不同在于是不是利用King-Smith算法优化，Auto那么自动选择，推荐。

  >ROPT:

优化操纵每一个核周围的积分球内的格点数，LREAL=AutoorOn

  ForLREAL=On

  PREC=Low700pointsintherealspacesphere（ROPT=

  PREC=Med1000pointsintherealspacesphere（ROPT=

  PREC=High1500pointsintherealspacesphere（ROPT=

  ForLREAL=Auto

  PREC=Lowaccuracy1e-2（ROPT=

  PREC=Medaccuracy2e-3（ROPT=

  PREC=Highaccuracy2e-4（ROPT=2E-4）

  >NELM,NELMINandNELMDLnr.ofelectronicsteps

  Default

  最大电子自洽循环次数NELM=60

  最小次数NELMIN=2

  弛豫次数NELMDL=-5ifISTART=0,INIWAV=1,andIALGO=8

  -12ifISTART=0,INIWAV=1,andIALGO=48

  0else

  若是初始的波函数采取随机赋值，即ISTART=0,INIWAV=1，那么极可能开始的值比较离

  谱，那么在第一步核运动循环之前采纳NELMDL（负值）步的非自洽（保留初始的H）步计算将

  减少计算所需的时刻。

  若是NELMDL取正值，将在每次核运动以后附加指定次数的弛豫步，目前不明白能够干吗

  >EDIFF电子SC循环的收敛精度缺省：

1e-4

  注意，即便EDIFF=0，NELM步也会执行

  >EDIFFG核运动的收敛精度缺省：

EDIFF*10（总能量）

  EDIFFG<0那么在所有的力都小于EDIFFG时停止

  EDIFFG=0那么在NSW步后停止

  此参数不支持MD，仅用于Relax

  >NSW指定核运动步数缺省:

0

  NBLOCKandKBLOCKinnerblock;outerblock

  DefaultNBLOCK=1KBLOCK=NSW

  >IBRIONionicrelaxation:

-1-Fixed0-MD1-quasi-New2-CG3-Damp5-freq

  DefaultifNSW=0or1IBRION=-1elseIBRION=0

  那个参数是和ISIF;IALGO/ALGO一路决定怎么算的最重要的参数

  1-3是三种Relax的方式，受ISIF决定是不是固定核位置、晶胞大小和形状

  0是标准的ab-initioMD，不受ISIF阻碍，即不改变晶胞大小和形状

  5可能是和0差不多吧？

支持Hessian和Freq（仅Г点）的计算和部份固定的MD

  详细的要在例如中具体情形具体分析了。

  >ISIFcalculatestressandwhattorelax

  DefaultifIBRION=0（MD）0else2

  ISIF│calculate│calculate│relax│change│change

  │force│stresstensor│ions│cellshape│cellvolume

  ──┼─────┼───────┼───┼──────┼──────

  0│yes│no│yes│no│no

  1│yes│traceonly│yes│no│no

  2│yes│yes│yes│no│no

  3│yes│yes│yes│yes│yes

  4│yes│yes│yes│yes│no

  5│yes│yes│no│yes│no

  6│yes│yes│no│yes│yes

  7│yes│yes│no│no│yes

  Traceonlymeansthatonlythetotalpressure

  IWAVPRpredictionofwf.:

0-non1-charg2-wave3-comb

  DefaultifIBRION=0（MD）2

  ifIBRION=1,2（relaxation）1

  else（staticcalculation）0

  以上选项保留TMPCAR+10那么全数利用内存，不保留此文件

  IWAVPRdetermineshowwavefunctionsand/orchargedensityare

  extrapolatedfromoneionicconfigurationtothenextconfiguration.

  >ISYMsymmetry:

0-nonsym1-usesym是不是利用对称性Default1

  SYMPRECdeterminesprecisionofthepositionsinPOSCARfile.Default1e-5

  LCORRHarris-correctiontoforces.Default.TRUE.

  >POTIMtime-stepforion-motion（fs）

  Default

  IBRION=0（MD）nodefault,必需指定，MD每步步长

  IBRION=1,2,3（relaxation）最小化的'scalingconstant',尤其是IBRION=1

  >TEBEG,TEENDtemperatureduringrun（MD有效）

  Default:

TEBEG=0TEEND=TEBEG

  注意VASP的温度概念与实际温度有细微的不同，因此

  TEBEG=T×（N-1）/NT为实际温度，N为原子数

  SMASS操纵MD中的速度模拟方式

  default-3微正那么系综（总自由能不变）

  -2维持初速度不变

  -1每NBLOCK步伐整速度，来保证动能持续

  >=0Nosé算法模拟正那么系综，（不懂-,-）

  NPACOandAPACO

  NPACO:

numberofslotsforpaircorrelation（PC）function.Default256

  APACO:

maximumdistancefortheevaluationofPCfunctioninA.Default16

  简单说确实是在不超过APACO的NPACO个距离上求成对相关函PCF

  RWIGSWigner-Seitz半径DOS计算用

  >NELECT总电子数

  若是系统不是电中性的就必需设置，所带电荷作为均一的背景电子气考虑

  NUPDOWNdefault不考虑电子自旋态改变的可能

  EMIN,EMAXenergy-rangeforDOSCARfile

  >ISMEARpart.occupancies:

-5tetwithBlochl-4-tet-1-fermi0-gaus>0MP

  采纳所谓部份占有波函数，用一个函数来滑腻积分，尤其是关于金属体系可减少k点

  DefaultISMEAR=1若是在KPOINTS里利用了tetrahedra方式推荐ISMEAR=5

  SIGMAdeterminesthewidthofthesmearingineV

  DefaultSIGMA=

  >ALGOalgorithm:

Normal（Davidson）|Fast（mixed）|Very_Fast（RMM-DIIS）

  >IALGOalgorithm:

useonly48（RMM-DIIS）or38（Davidson）or8（CG）

  DefaultIALGO=38for

  算法是最重要的参数之一。

一样VASP推荐利用的是以上三种算法，一样来讲8/38是初期

  比较快收敛，在接近平稳时采纳48较快，在初期或MD时利用48可能会碰到不收敛的情形

  。

也能够利用ALGO参数来替代IALGO，设置Fast，VASP会先用38，再自动切换到48。

  各类算法只要收敛，结果应该一致。

  另一个可能有效的选项是-1。

不进行实际的计算，只对重要的步骤做计算测试，并将测

  试取得的各部份耗时输出在OUTPUT里。

  VOSKOWNuseVWNinterpolation算法，default0不用，若是利用了PW91或需要计

  算磁性质，不妨设为1用

  mixingtag&MAXMIX

  IMIX=typeofmixing

  AMIX=linearmixingparameter

  AMIN=minimalmixingparameter

  BMIX=cutoffwavevectorforKerkermixingscheme

  AMIX_MAG=linearmixingparameterformagnetization

  BMIX_MAG=cutoffwavevectorforKerkermixingschemeformag.

  WC=weightfactorforeachstepinBroydenmixingscheme

  INIMIX=typeofinitialmixinginBroydenmixingscheme

  MIXPRE=typeofpreconditioninginBroydenmixingscheme

  MAXMIX=maximumnumberstepsstoredinBroydenmixer

  Default

  forUS-PPnon-magnetic

  IMIX=4

  AMIX=

  BMIX=

  WC=1000.

  INIMIX=1

  MIXPRE=1

  MAXMIX=-45

  值得注意的是，在MD或Relax的时候，设置MAXMIX（>0,一样约3倍的电子SC步数），可能

  会大大减少核运动步数。

可是同时会增加对内存的要求。

  LWAVE,LCHARGandLVTOTcreateWAVECAR/CHGCAR/LOCPOT

  LELFcreateELFCAR

  LORBITcreatePROOUT

  输出文件的选项

  NPAR并行计算band的节点数，每一个节点计算一个band固然能够提高并行效率，减少通

  讯量，只是貌似此刻硬件的要紧限制仍是内存，而那个选项的利用可能会大幅增加内存

  的需求

  >NBANDS总能带数。

  之因此把它放在最后，是因为它关于解决内存需求的重要性。

计算需要大量的

  能带（空带），至少要1个空带（不然VASP会给出警告）。

一样NBANDS=NELECT/2+NIONS/2

  以上可取得较精准的结果，若是内存不够就只好减少NBANDS，在捐躯精度和体系大小之

  间平稳了。

  最后提示一下大多数参数的首字母代表了参数的性质

  I初始化L逻辑开关E能量N数量T温度

  便于经历:

）