当前位置:文档之家› VASP控制参数文件INCAR的简单介绍

VASP控制参数文件INCAR的简单介绍

限于能力,只对部分最基本的一些参数(>,没有这个标志的参数都是可以不出现的)详细说明,在这里只是简单介绍这些参数的设置,详细的问题在后文具体示例中展开。

部分可能会干扰VASP运行的参数在这里被刻意隐去了,需要的同学还是请查看VASP自带的帮助文档原文。

参数列表如下:>SYSTEM name of System任务的名字***>NWRITE verbosity write-flag (how much is written)输出内容详细程度0-3 缺省2如果是做长时间动力学计算的话最好选0或1(首末步/每步核运动输出)据说也可以结合shell的tail或grep命令手动输出>ISTART startjob:restart选项0-3 缺省0/1 for 无/有前次计算的WAVECAR(波函数)1 'restart with constant energy cut-off'2 'restart with constant basis set'3 'full restart including wave function and charge prediction'ICHARG charge: 1-file 2-atom 10-const Default:if ISTART=0 2 else 0ISPIN spin polarized calculation (2-yes 1-no) default 2MAGMOM initial mag moment / atom Default NIONS*1INIWAV initial electr wf. : 0-lowe 1-randDefault 1 only used for start jobs (ISTART=0)IDIPOL calculate monopole/dipole and quadrupole corrections1-3 只计算第一/二/三晶矢方向适于slab的计算4 全部计算尤其适于就算孤立分子>PREC precession: medium, high or low(VASP.4.5+ also: normal, accurate)Default: Medium VASP4.5+采用了优化的accurate来替代high,所以一般不推荐使用high。

不过high可以确保'绝对收敛',作为参考值有时也是必要的。

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

受PREC影响的参数有四类:ENCUT; NGX,NGY,NGZ; NGXF, NGYF, NGZF; ROPT如果设置了PREC,这些参数就都不需要出现了当然直接设置相应的参数也是同样效果的,这里不展开了,随后详释>ENCUT energy cutoff in eV : default taken from POTCAR-fileimportant! 重要到几乎最好不要手工去设置除非文献告诉你要用多少,或者经过结果可靠性的验证当然,为了测试一下提交的任务,也不妨先设个较小的值附加说明:当且仅当POTCAR里头没有设置ENCUT时(其实貌似没有才是常态),才受PREC设置影响从POTCAR里找出相应的ENMAX/ENMIN值来设置。

PREC= Low Medium Accurate HighENCUT= ENMIN ENMAX ENMAX 130%ENMAX对于多个元素的POTCAR不同的ENMAX/ENMIN,都取最大值>NGX,NGY,NGZ: FFT mesh for wavefunctions>NGFX, NGFY, NGFZ: FFT mesh for charges也是两类重要的最好不要去动的参数,PREC设置将从POTCAR中自动读取。

PREC=High,Accurate 2倍值,用来避免wrap around errors得到精确解PREC=Low,Medium,Normal 3/4 也已经足够精确到1 meV/atom>LREAL: Default= .FALSE.赝势的非局域部分用到的一个积分在倒格空间或者实空间都可以求值。

这个选项就是决定是在哪个空间里求。

在倒格空间里,采用平面波基组求解,在实空间里,采用积分球求解。

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

但效率会低一些。

其他选项是O or On,A or Auto 和.True.。

On和.TRUE.的差别在于是否使用King-Smith算法优化,Auto则自动选择,推荐。

>ROPT: 优化控制每个核周围的积分球内的格点数,LREAL=Auto or OnFor LREAL=OnPREC= Low 700 points in the real space sphere ( ROPT=0.67)PREC= Med 1000 points in the real space sphere ( ROPT=1.0)PREC= High 1500 points in the real space sphere ( ROPT=1.5)For LREAL=AutoPREC= Low accuracy 1e-2 ( ROPT=0.01)PREC= Med accuracy 2e-3 ( ROPT=0.002)PREC= High accuracy 2e-4 ( ROPT=2E-4)>NELM, NELMIN and NELMDL nr. of electronic stepsDefault最大电子自洽循环次数NELM = 60最小次数NELMIN = 2弛豫次数NELMDL = -5 if ISTART=0, INIWAV=1, and IALGO=8-12 if ISTART=0, INIWAV=1, and IALGO=480 else如果初始的波函数采取随机赋值,即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 指定核运动步数缺省: 0NBLOCK and KBLOCK inner block; outer blockDefault NBLOCK = 1 KBLOCK = NSW>IBRION ionic relaxation: -1-Fixed 0-MD 1-quasi-New 2-CG 3-Damp 5-freqDefault if NSW=0 or 1 IBRION=-1 else IBRION=0这个参数是和ISIF;IALGO/ALGO一起决定怎么算的最重要的参数1-3 是三种Relax的方法,受ISIF决定是否固定核位置、晶胞大小和形状0 是标准的ab-initio MD,不受ISIF影响,即不改变晶胞大小和形状5 大概是和0差不多吧?支持Hessian和Freq(仅Г点)的计算以及部分固定的MD详细的要在示例中具体情况具体分析了。

>ISIF calculate stress and what to relaxDefault if IBRION=0 (MD) 0 else 2ISIF│calculate │calculate │relax │change │change│force │stress tensor │ions │cell shape │cell volume──┼─────┼───────┼───┼──────┼──────0 │yes │no │yes │no │no1 │yes │trace only │yes │no │no2 │yes │yes │yes │no │no3 │yes │yes │yes │yes │yes4 │yes │yes │yes │yes │no5 │yes │yes │no │yes │no6 │yes │yes │no │yes │yes7 │yes │yes │no │no │yesTrace only means that only the total pressureIWAVPR prediction of wf.: 0-non 1-charg 2-wave 3-combDefault if IBRION=0 (MD) 2if IBRION=1,2 (relaxation) 1else (static calculation) 0以上选项保存TMPCAR +10 则全部使用内存,不保存此文件IWAVPR determines how wave functions and/or charge density are extrapolated from one ionic configuration to the next configuration.>ISYM symmetry: 0-nonsym 1-usesym 是否使用对称性Default 1SYMPREC determines precision of the positions in POSCAR file. Default 1e-5LCORR Harris-correction to forces. Default .TRUE.>POTIM time-step for ion-motion (fs)DefaultIBRION=0 (MD) no default,必须指定,MD每步步长IBRION=1,2,3 (relaxation) 0.5 最小化的'scaling constant',尤其是IBRION=1>TEBEG, TEEND temperature during run (MD有效)Default: TEBEG = 0 TEEND = TEBEG注意VASP的温度定义与实际温度有细微的差别,所以TEBEG=T×(N-1)/N T为实际温度,N为原子数SMASS 控制MD中的速度模拟方法default -3 微正则系综(总自由能不变)-2 保持初速度不变-1 每NBLOCK步调整速度,来保证动能连续>=0 Nosé算法模拟正则系综,(不懂-,-)NPACO and APACONPACO : number of slots for pair correlation (PC) function. Default 256APACO : maximum distance for the evaluation of PC function in A. Default 16 简单说就是在不超过APACO的NPACO个距离上求成对相关函PCFRWIGS Wigner-Seitz半径DOS计算用>NELECT 总电子数如果系统不是电中性的就必须设置,所带电荷作为均一的背景电子气考虑NUPDOWN default不考虑电子自旋态改变的可能EMIN, EMAX energy-range for DOSCAR file>ISMEAR part. occupancies: -5 tet with Blochl -4-tet -1-fermi 0-gaus >0 MP采用所谓部分占有波函数,用一个函数来平滑积分,尤其是对于金属体系可减少k点Default ISMEAR = 1 如果在KPOINTS里使用了tetrahedra方法推荐ISMEAR=5SIGMA determines the width of the smearing in eVDefault SIGMA = 0.2>ALGO algorithm: Normal (Davidson) | Fast (mixed)| Very_Fast (RMM-DIIS)>IALGO algorithm: use only 48 (RMM-DIIS) or 38(Davidson) or 8(CG)Default IALGO = 38 for VASP4.5算法是最重要的参数之一。

相关主题