INCAR简单介绍讲课教案

I N C A R简单介绍

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 0

ISPIN spin polarized calculation (2-yes 1-no) default 2 MAGMOM initial mag moment / atom Default NIONS*1 INIWAV initial electr wf. : 0-lowe 1-rand Default 1 only used for start jobs (ISTART=0)

IDIPOL corrections

calculate monopole/dipole and quadrupole

1-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-file

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

除非文献告诉你要用多少，或者经过结果可靠性的验证 当然，为了测试一下提交的任务，也不妨先设个较小的值 附加说明：

当且仅当POTCAR里头没有设置ENCUT时(其实貌似没有才是常态)，才受PREC设置影响从POTCAR里找出相应的ENMAX/ENMIN值来设置。

PREC= Low Medium Accurate High

ENCUT= 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 On

For LREAL=On

PREC= 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=Auto

PREC= 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 steps Default

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

弛豫次数 NELMDL = -5 if ISTART=0, INIWAV=1, and IALGO=8 -12 if ISTART=0, INIWAV=1, and IALGO=48 0 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 指定核运动步数缺省: 0

NBLOCK and KBLOCK inner block; outer block Default NBLOCK = 1 KBLOCK = NSW

>IBRION ionic relaxation: -1-Fixed 0-MD 1-quasi-New 2-CG 3-Damp 5-freq

Default 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 relax Default if IBRION=0 (MD) 0 else 2

ISIF│calculate │ calculate │relax │ change │ change │ force │stress tensor │ions │ cell shape │cell volume ──┼─────┼───────┼───┼──────┼────── 0 │ yes │ no │yes │ no │no 1 │ yes│ trace only │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

Trace only means that only the total pressure

IWAVPR prediction of wf.: 0-non 1-charg 2-wave 3-comb Default if IBRION=0 (MD) 2 if IBRION=1,2 (relaxation) 1