polt-包管
INCAR
是决定
how to
do
的文件
限于能力,只对部分最基
本的一些参数
(>,
没有这个标志的参数都是可以不
出现的
)
详细说明,在这里只是简单介绍这些参
数的设置,详细的问题在后文具体
示例中展开。
部分可能会干扰
VASP
运行的参数在这里被刻意隐
去了,需要的同学还是请
查看
VASP
自带的帮助文档原文。
参数列表如下:
>SYSTEM
name of System
任务的名字
***
>NWRITE verbosity write-flag (how much
is written)
输出内容详细程度
0-3
缺省
2
如果是做长时间动力学计算的话最好选<
/p>
0
或
1(
首末步
/
每步核运动输出
)
据说也可以结合
shell
的
tail
或
grep
命令手动输出
p>
>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
1
/
10
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
< br>的计算
4
全部计算尤其适于就算孤立分子
>PREC precession: medium, high or
low(VASP.4.5+ also: normal, accurate)
Default: Medium VASP4.5+
采用了优
化的
accurate
来替代
high
,
所以一般不推
荐使用
high
。不过
high
可以确保
'
绝对收敛
'
,作为参考值有时也
是必要的。
同样受推荐的是
norm
al
,作为日常计算选项,可惜的是说明文档提供的信
息不足。
受
PREC
影响的参数有四类:
ENCUT; NGX,NGY,NGZ; NGXF,
NGYF, NGZF; ROPT
如果设置了
PREC
p>
,这些参数就都不需要出现了
当然直接设
置相应的参数也是同样效果的,这里不展开了,随后详释
>ENCUT energy cutoff in eV : default
taken from POTCAR-file
important!
重要到几乎最好不要手工去设置
除非文献告诉你要用多少,或者经过结果可靠性的验证
当然,为了测试一下提交的任务,也不妨先设个较小的值
附加说明:
当且仅当
POTCAR
里头没有设置
ENCUT
时
(
其实貌似没有才是常态
)
,才受
PREC
设置影响从
POTCAR
里找出相应的
ENMAX/ENMI
N
值来设置。
PREC= Low
Medium Accurate High
ENCUT= ENMIN ENMAX
ENMAX 130%ENMAX
2
/
10
<
/p>
对于多个元素的
POTCAR
不同的
p>
ENMAX/ENMIN
,都取最大值
>NGX,NGY,NGZ: FFT mesh for
wavefunctions
>NGFX, NGFY, NGFZ: FFT
mesh for charges
也是两类重要的最好不要去动的参数,
PREC
设置将从
POTCAR
< br>中自动读
取。
PREC=High,Accurate
2
倍值,用来避免
wrap around
errors
得到精确解
PREC=Low,Medium,Normal
3/4
也已经足够精确到
1 meV/atom
>LREAL: Default= .FALSE.
赝势的
非局域部分用到的一个积分在倒格空间或者实空间都可以求值。这
个选项就
是决定是在哪个空间里求。在倒格空间里,采用平面波基组求解,在实空<
/p>
间里,采用积
分球求解。
缺省是
< br>.FALSE
,即不在实空间求。但效率会低一些。
其他选项是
O or On,A or Auto
和
.True.
。
On
和
.TRUE.
的差别在于是否使用
King-Smith
算法优化,
Auto
则自动选择,
推荐。
< br>
>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)
3
/
10
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, IN
IWAV=1
,那么很可能开始
的值比较离谱,那么在第一步核
运动循环之前采用
NELMDL(
负值
)
步的非自洽
(
保
留初始的
H)
步计算将减少计算所需的时间。
如果
NELMDL
取正值
,将在每次核运动之后附加指定次数的弛豫步,目前
不知道可以干嘛
>EDIFF
电子
SC
循环的收敛精度缺省:
1e-4
注意,即使<
/p>
EDIFF=0
,
NELM
步也会执行
>EDIFFG
核运动的收敛精度缺省:
EDIFF*10
(
总能量
)
EDIFFG<0
则在所有的力都小于
EDIFFG
时停止
p>
EDIFFG=0
则在
< br>NSW
步后停止
此参数不支持
MD
,仅用于
Relax
4
/
10
>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(
仅
Г
点
)
的计算以及部分固定
的<
/p>
MD
详细的要在示例中具体情况具体分析了。
>ISIF calculate stress and what to
relax
Default if IBRION=0 (MD) 0 else 2
ISIF│calculate│calculate│relax│change│c
hange
│force│stress
tensor│ions│cell shape│cell volume
──
┼
─────
┼<
/p>
───────
┼
───
┼
──────
┼
─────
─
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
5
/
10