转http://blog.chinaunix.net/uid-20676695-id-1899605.html
INCAR是决定how
to do 的文件
限于能力，只对部分最基本的一些参数(>,没有这个标志的参数都是可以不出现的)
详细说明，在这里只是简单介绍这些参数的设置，详细的问题在后文具体示例中展开。
部分可能会干扰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 calculate monopole/dipole and quadrupole
corrections
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一起决定怎么算的最重要的参数
加载中…
(2013-03-12 10:25:48)