Chemistry 版 (精华区)

发信人: zjliu (秋天的萝卜), 信区: Chemistry
标  题: 初学VASP(六) 最重要的INCAR参数
发信站: BBS 哈工大紫丁香站 (Thu Jun 17 20:50:08 2004)

发信站: 日月光华  

INCAR是决定how to do 的文件

限于能力，只对部分最基本的一些参数(>)，也是计算必须涉及的最重要的参数详细说
明，在这里只是简单介绍这些参数的设置，详细的问题在后文具体示例中展开

参数列表如下：

>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'
    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)

>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
    当然直接设置相应的参数也是同样效果的，这里不展开了，到相应的参数处再详释

>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
>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 电子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
    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
             else (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 1
>ISYM symmetry: 0-nonsym 1-usesym 是否使用对称性 Default 1

SYMPREC determines precision of the positions in POSCAR file. Default 1e-5

LCORR Harris-correction to forces. Default .TRUE.

>POTIM time-step for ion-motion (fs)
    Default
     IBRION=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 Nose mass-parameter (am)

NPACO and APACO distance and nr. of slots for P.C.

POMASS mass of ions in am

ZVAL ionic valence
ZVAL ionic valence

RWIGS Wigner-Seitz radii

NELECT total number of electrons

NUPDOWN fix spin moment to specified value

EMIN, EMAX energy-range for DOSCAR file

ISMEAR part. occupancies: -5 Bl?chl -4-tet -1-fermi 0-gaus >0 MP

SIGMA broadening in eV -4-tet -1-fermi 0-gaus

ALGO algorithm: Normal (Davidson) | Fast (mixed)| Very_Fast (RMM-DIIS)

IALGO algorithm: use only 8 (CG) or 48 (RMM-DIIS) or 38(Davidson)

VOSKOWN use Vosko, Wilk, Nusair interpolation

DIPOL center of cell for dipol

AMIX, BMIX tags for mixing


WEIMIN, EBREAK, DEPER special control tags

TIME special control tag

LWAVE,LCHARG and LVTOT create WAVECAR/CHGCAR/LOCPOT

LELF create ELFCAR

LORBIT create PROOUT

NPAR parallelization over bands

LSCALAPACK switch of scaLAPACK

LSCALU switch of LU decomposition

LASYNC overlap communcation with calculations
--
╔═══════════════════╗
║★★★★★友谊第一  比赛第二★★★★★║
╚═══════════════════╝


※ 来源:·哈工大紫丁香 http://bbs.hit.edu.cn·[FROM: 202.118.229.*]