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 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一起决定怎么算的最重要的参数
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
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
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 控制MD中的速度模拟方法
default -3 微正则系综(总自由能不变)
-2 保持初速度不变
-1 每NBLOCK步调整速度，来保证动能连续
>=0 Nosé算法模拟正则系综，(不懂-,-)
NPACO and APACO
NPACO : number of slots for pair correlation (PC) function. Default 256
APACO : maximum distance for the evaluation of PC function in A. Default 16 简单说就是在不超过APACO的NPACO个距离上求成对相关函PCF
RWIGS 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=5
SIGMA determines the width of the smearing in eV
Default 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
算法是最重要的参数之一。一般VASP推荐使用的是以上三种算法，一般来说8/38是初期比较快收敛，在接近平衡时采用48较快，在初期或MD时使用48可能会遇到不收敛的情况。也可以使用ALGO参数来替代IALGO，设置Fast，VASP会先用38，再自动切换到48。各种算法只要收敛，结果应该一致。
另一个可能有用的选项是-1。不进行实际的计算，只对重要的步骤做计算测试，并将测
试得到的各部分耗时输出在OUTPUT里。
VOSKOWN use VWN interpolation 算法，default 0 不用，如果使用了PW91或者需要计算磁性质，不妨设为1 用
mixing tag & MAXMIX
IMIX = type of mixing
AMIX = linear mixing parameter
AMIN = minimal mixing parameter
BMIX = cutoff wave vector for Kerker mixing scheme
AMIX_MAG = linear mixing parameter for magnetization
BMIX_MAG = cutoff wave vector for Kerker mixing scheme for mag.
WC = weight factor for each step in Broyden mixing scheme
INIMIX = type of initial mixing in Broyden mixing scheme
MIXPRE = type of preconditioning in Broyden mixing scheme
MAXMIX = maximum number steps stored in Broyden mixer
Default
for US-PP non-magnetic
IMIX = 4
AMIX = 0.8
BMIX = 1.0
WC = 1000.
INIMIX = 1
MIXPRE = 1
MAXMIX = -45
值得注意的是，在MD或者Relax的时候，设置MAXMIX(>0,一般约3倍的电子SC步数)，可能
会大大减少核运动步数。但是同时会增加对内存的要求。
LWAVE,LCHARG and LVTOT create WAVECAR/CHGCAR/LOCPOT
LELF create ELFCAR
LORBIT create PROOUT
输出文件的选项
NPAR 并行计算band的节点数，每一个节点计算一个band当然可以提高并行效率，减少通讯量，不过貌似现在硬件的主要限制还是内存，而这个选项的使用可能会大幅增加内存的需求
>NBANDS 总能带数。之所以把它放在最后，是因为它对于解决内存需求的重要性。计算需要大量的能带(空带)，至少要1个空带(否则VASP会给出警告)。一般NBANDS=NELECT/2+NIONS/2 以上可得到较精确的结果，如果内存不够就只好减少NBANDS，在牺牲精度和体系大小之间平衡了
最后提示一下大多数参数的首字母代表了参数的性质
I 初始化L 逻辑开关E 能量N 数目T 温度便于记忆:)
vasp参数详解
转载]VASP参数设置详解 已有1283 次阅读2011-9-15 14:26 |个人分类:vasp|系统分类:科研笔记 转自小木虫，略有增减 软件主要功能： 采用周期性边界条件（或超原胞模型）处理原子、分子、团簇、纳米线（或管）、薄膜、晶体、准晶和无定性材料，以及表面体系和固体 l 计算材料的结构参数（键长、键角、晶格常数、原子位置等）和构型 l 计算材料的状态方程和力学性质（体弹性模量和弹性常数） l 计算材料的电子结构（能级、电荷密度分布、能带、电子态密度和ELF） l 计算材料的光学性质 l 计算材料的磁学性质 l 计算材料的晶格动力学性质（声子谱等） l 表面体系的模拟（重构、表面态和STM模拟） l 从头分子动力学模拟 l 计算材料的激发态（GW准粒子修正） 计算主要的四个参数文件：INCAR ,POSCAR,POTCAR ,KPOINTS,下面简要介绍，详细权威的请参照手册 INCAR文件： 该文件控制VASP进行何种性质的计算，并设置了计算方法中一些重要的参数，这些参数主要包括以下几类： l 对所计算的体系进行注释：SYSTEM l 定义如何输入或构造初始的电荷密度和波函数：ISTART，ICHARG，INIWAV l 定义电子的优化 –平面波切断动能和缀加电荷时的切断值：ENCUT，ENAUG –电子部分优化的方法：ALGO，IALGO，LDIAG –电荷密度混合的方法：IMIX，AMIX，AMIN，BMIX，AMIX_MAG，BMIX_MAG，WC，INIMIX，MIXPRE，MAXMIX –自洽迭代步数和收敛标准：NELM，NELMIN，NELMDL，EDIFF l 定义离子或原子的优化 –原子位置优化的方法、移动的步长和步数：IBRION，NFREE，POTIM，NSW –分子动力学相关参数：SMASS，TEBEG，TEEND，POMASS，NBLOCK，KBLOCK，PSTRESS –离子弛豫收敛标准：EDIFFG l 定义态密度积分的方法和参数 – smearing方法和参数：ISMEAR，SIGMA –计算态密度时能量范围和点数：EMIN，EMAX，NEDOS –计算分波态密度的参数：RWIGS，LORBIT l 其它 –计算精度控制：PREC –磁性计算：ISPIN，MAGMOM，NUPDOWN –交换关联函数：GGA，VOSKOWN –计算ELF和总的局域势：LELF，LVTOT –结构优化参数：ISIF –等等。 主要参数说明如下：
VASP参数设置详解
VASP参数设置详解 计算材料2010-11-30 20:11:32 阅读197 评论0 字号：大中小订阅 转自小木虫，略有增减 软件主要功能： 采用周期性边界条件（或超原胞模型）处理原子、分子、团簇、纳米线（或管）、薄膜、晶体、准晶和无定性材料，以及表面体系和固体 l 计算材料的结构参数（键长、键角、晶格常数、原子位置等）和构型 l 计算材料的状态方程和力学性质（体弹性模量和弹性常数） l 计算材料的电子结构（能级、电荷密度分布、能带、电子态密度和ELF） l 计算材料的光学性质 l 计算材料的磁学性质 l 计算材料的晶格动力学性质（声子谱等） l 表面体系的模拟（重构、表面态和STM模拟） l 从头分子动力学模拟 l 计算材料的激发态（GW准粒子修正） 计算主要的四个参数文件：INCAR ,POSCAR,POTCAR ,KPOINTS,下面简要介绍，详细权威的请参照手册 INCAR文件： 该文件控制VASP进行何种性质的计算，并设置了计算方法中一些重要的参数，这些参数主要包括以下几类： 对所计算的体系进行注释：SYSTEM
●定义如何输入或构造初始的电荷密度和波函数：ISTART，ICHARG，INIWAV ●定义电子的优化 –平面波切断动能和缀加电荷时的切断值：ENCUT，ENAUG –电子部分优化的方法：ALGO，IALGO，LDIAG –电荷密度混合的方法：IMIX，AMIX，AMIN，BMIX，AMIX_MAG，BMIX_MAG，WC，INIMIX，MIXPRE，MAXMIX –自洽迭代步数和收敛标准：NELM，NELMIN，NELMDL，EDIFF ●定义离子或原子的优化 –原子位置优化的方法、移动的步长和步数：IBRION，NFREE，POTIM，NSW –分子动力学相关参数：SMASS，TEBEG，TEEND，POMASS，NBLOCK，KBLOCK，PSTRESS –离子弛豫收敛标准：EDIFFG ●定义态密度积分的方法和参数 –smearing方法和参数：ISMEAR，SIGMA –计算态密度时能量范围和点数：EMIN，EMAX，NEDOS –计算分波态密度的参数：RWIGS，LORBIT ●其它 –计算精度控制：PREC –磁性计算：ISPIN，MAGMOM，NUPDOWN –交换关联函数：GGA，VOSKOWN –计算ELF和总的局域势：LELF，LVTOT –结构优化参数：ISIF –等等。 主要参数说明如下： ?SYSTEM：该输入文件所要执行的任务的名字。取值：字符串，缺省值：SYSTEM ?NWRITE：输出内容详细程度。取值：0~4，缺省值：2
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 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，这些参数就都不需要出现了 当然直接设置相应的参数也是同样效果的，这里不展开了，随后详释
VASP-INCAR参数设置
V A S P-I N C A R参数设置-CAL-FENGHAI-(2020YEAR-YICAI)_JINGBIAN
1. 结构优化 (Opt) SYSTEM = opt ISTART = 0 INIWAV = 1 ICHARG = 2 ISPIN = 2 LREAL = Auto ENCUT = 400 PREC = high NSW= 600 NELM = 60 IBRION = 2 ISIF = 2 POTIM = 0.1 ALGO= Fast LVDW = .TRUE. EDIFF = 1E-5 EDIFFG = 1E-4 or -0.05 # 体系需计算TS时，全部结构优化EDIFFG均设置为-0.05 ISMEAR = 0 SIGMA = 0.2 LCHARG = .FALSE. LWAVE = .FALSE.
2. 过渡态搜索 (TS): 计算时先进行低精度计算，再进行高精度计算 SYSTEM= TS ISTART = 0 INIWAV = 1 ICHARG = 2 ISPIN = 2 LREAL = Auto ENCUT = 400 PREC = high NSW = 600 NELMIN = 6 IBRION = 3 or 1 # 过渡态计算低精度为3，高精度为1 ISIF = 2 POTIM = 0.01 ALGO = Fast LVDW = .TRUE. EDIFF = 1E-5 EDIFFG = -1 or -0.05 # 过渡态计算低精度为-1，高精度为-0.05 ISMEAR = 0 SIGMA = 0.05 LCHARG= .FALSE. LWAVE= .FALSE. IMAGES=8 # TS专属设置 SPRING=-5 # TS专属设置 LCLIMB=.TRUE. # TS专属设置
vasp计算参数设置
软件主要功能： 采用周期性边界条件（或超原胞模型）处理原子、分子、团簇、纳米线（或管）、薄膜、晶体、准晶和无定性材料，以及表面体系和固体 l 计算材料的结构参数（键长、键角、晶格常数、原子位置等）和构型 l 计算材料的状态方程和力学性质（体弹性模量和弹性常数） l 计算材料的电子结构（能级、电荷密度分布、能带、电子态密度和ELF） l 计算材料的光学性质 l 计算材料的磁学性质 l 计算材料的晶格动力学性质（声子谱等） l 表面体系的模拟（重构、表面态和STM模拟） l 从头分子动力学模拟 l 计算材料的激发态（GW准粒子修正） 计算主要的四个参数文件：INCAR ,POSCAR,POTCAR ,KPOINTS,下面简要介绍，详细权威的请参照手册 INCAR文件： 该文件控制VASP进行何种性质的计算，并设置了计算方法中一些重要的参数，这些参数主要包括以下几类： l 对所计算的体系进行注释：SYSTEM l 定义如何输入或构造初始的电荷密度和波函数：ISTART，ICHARG，INIWA V l 定义电子的优化 –平面波切断动能和缀加电荷时的切断值：ENCUT，ENAUG –电子部分优化的方法：ALGO，IALGO，LDIAG –电荷密度混合的方法：IMIX，AMIX，AMIN，BMIX，AMIX_MAG，BMIX_MAG，WC，INIMIX，MIXPRE，MAXMIX –自洽迭代步数和收敛标准：NELM，NELMIN，NELMDL，EDIFF l 定义离子或原子的优化 –原子位置优化的方法、移动的步长和步数：IBRION，NFREE，POTIM，NSW –分子动力学相关参数：SMASS，TEBEG，TEEND，POMASS，NBLOCK，KBLOCK，PSTRESS –离子弛豫收敛标准：EDIFFG l 定义态密度积分的方法和参数 –smearing方法和参数：ISMEAR，SIGMA –计算态密度时能量范围和点数：EMIN，EMAX，NEDOS –计算分波态密度的参数：RWIGS，LORBIT l 其它 –计算精度控制：PREC –磁性计算：ISPIN，MAGMOM，NUPDOWN –交换关联函数：GGA，VOSKOWN –计算ELF和总的局域势：LELF，LVTOT –结构优化参数：ISIF –等等。 主要参数说明如下： ? SYSTEM：该输入文件所要执行的任务的名字。取值：字符串，缺省值：SYSTEM
初学VASP中电子态密度计算设置参考
初学VASP中电子态密度计算基本设置参考主要分成三步：一、结构优化；二、静态自洽计算；三、非自洽计算以Al－FCC为例子 第一步结构优化 输入文件（INCAR, POTCAR, POSCAR, KPOINT） INCAR文件 System=Al ISTART=0 ISMEAR=1 SIGMA=0.2 ISPIN=2 GGA=91; VOSKOWN=1; EDIFF=0.1E-05; EDIFFG=-0.01 IBRION=2 NSW=50 ISIF=2 (OR 3) NPAR=10 POTCAR 文件直接在势库中拷贝 POSCAR文件 Al 4.05 1.0 0.0 0.0 0.0 1.0 0.0
0.0 0.0 1.0 4 Direct 0.0 0.0 0.0 0.5 0.5 0.0 0.5 0.0 0.5 0.0 0.5 0.5 KPOINT 文件 Automatic generation Mohkorst Pack 15 15 15 0.0 0.0 0.0 第二步静态自洽计算 INCAR：PREC = Medium，ISTART = 0，ICHARG = 2，ISMEAR = -5输入文件（INCAR, POTCAR, POSCAR, KPOINT） INCAR文件 System=Al ISTART=0 ISMEAR=1 SIGMA=0.2 ISPIN=2
GGA=91; VOSKOWN=1; EDIFF=0.1E-05; EDIFFG=-0.01 #IBRION=2 #NSW=50 #ISIF=2 (OR 3) NPAR=10 POTCAR 文件直接在势库中拷贝 POSCAR文件 Al 4.05 1.0 0.0 0.0 0.0 1.0 0.0 0.0 0.0 1.0 4 Selective Dynamic Direct 0.0 0.0 0.0 T T T 0.5 0.5 0.0 T T T 0.5 0.0 0.5 T T T 0.0 0.5 0.5 T T T KPOINT 文件 Automatic generation
VASP遇到小总结问题
VASP 计算的过程遇到的问题 01、第一原理计算的一些心得 （1）第一性原理其实是包括基于密度泛函的从头算和基于Hartree-Fock自洽计算的从头算，前者以电子密度作为基本变量（霍亨伯格-科洪定理），通过求解Kohn-Sham方程，迭代自洽得到体系的基态电子密度，然后求体系的基态性质；后者则通过自洽求解Hartree-Fock方程，获得体系的波函数，求基态性质； 评述：K-S方程的计算水平达到了H-F水平，同时还考虑了电子间的交换关联作用。 （2）关于DFT中密度泛函的Functional，其实是交换关联泛函 包括LDA，GGA，杂化泛函等等 一般LDA为局域密度近似，在空间某点用均匀电子气密度作为交换关联泛函的唯一变量，多数为参数化的CA-PZ方案； GGA为广义梯度近似，不仅将电子密度作为交换关联泛函的变量，也考虑了密度的梯度为变量，包括PBE,PW,RPBE等方案，BL YP泛函也属于GGA； 此外还有一些杂化泛函，B3L YP等。 （3）关于赝势 在处理计算体系中原子的电子态时，有两种方法，一种是考虑所有电子，叫做全电子法，比如WIEN2K中的FLAPW方法(线性缀加平面波)；此外还有一种方法是只考虑价电子，而把芯电子和原子核构成离子实放在一起考虑，即赝势法，一般赝势法是选取一个截断半径，截断半径以内，波函数变化较平滑，和真实的不同，截断半径以外则和真实情况相同，而且赝势法得到的能量本征值和全电子法应该相同。 赝势包括模守恒和超软，模守恒较硬，一般需要较大的截断能，超软势则可以用较小的截断能即可。另外，模守恒势的散射特性和全电子相同，因此一般红外，拉曼等光谱的计算需要用模守恒势。 赝势的测试标准应是赝势与全电子法计算结果的匹配度，而不是赝势与实验结果的匹配度，因为和实验结果的匹配可能是偶然的。 （4）关于收敛测试 （a）Ecut，也就是截断能，一般情况下，总能相对于不同Ecut做计算，当Ecut增大时总能变化不明显了即可；然而，在需要考虑体系应力时，还需对应力进行收敛测试，而且应力相对于Ecut的收敛要比总能更为苛刻，也就是某个截断能下总能已经收敛了，但应力未必收敛。 （b）K-point，即K网格，一般金属需要较大的K网格，采用超晶胞时可以选用相对较小的K网格，但实际上还是要经过测试。 （5）关于磁性 一般何时考虑自旋呢？举例子，例如BaTiO3中，Ba、Ti和O分别为+2，+4和-2价，离子全部为各个轨道满壳层的结构，就不必考虑自旋了；对于BaMnO3中，由于Mn+3价时d 轨道还有电子，但未满，因此需考虑Mn的自旋，至于Ba和O则不必考虑。其实设定自旋就是给定一个原子磁矩的初始值，只在刚开始计算时作为初始值使用，具体的可参照磁性物理。 （6）关于几何优化 包括很多种了，比如晶格常数和原子位置同时优化，只优化原子位置，只优化晶格常数，还有晶格常数和原子位置分开优化等等。
VASP经典学习教程,有用
VASP 学习教程太原理工大学量子化学课题组 2012/5/25 太原
目录 第一章Linux命令 (1) 1.1 常用命令 (1) 1.1.1 浏览目录 (1) 1.1.2 浏览文件 (1) 1.1.3 目录操作 (1) 1.1.4 文件操作 (1) 1.1.5 系统信息 (1) 第二章SSH软件使用 (2) 2.1 软件界面 (2) 2.2 SSH transfer的应用 (3) 2.2.1 文件传输 (3) 2.2.2 简单应用 (3) 第三章V ASP的四个输入文件 (3) 3.1 INCAR (3) 3.2 KPOINTS (4) 3.3 POSCAR (4) 3.4 POTCAR (5) 第四章实例 (5) 4.1 模型的构建 (5) 4.2 V ASP计算 (8) 4.2.1 参数测试 (8) 4.2.2 晶胞优化(Cu) (13) 4.2.3 Cu(100)表面的能量 (2) 4.2.4 吸附分子CO、H、CHO的结构优化 (2) 4.2.5 CO吸附于Cu100表面H位 (4) 4.2.6 H吸附于Cu100表面H位 (5) 4.2.7 CHO吸附于Cu100表面B位 (6) 4.2.8 CO和H共吸附于Cu100表面 (7) 4.2.9 过渡态计算 (8)
第一章Linux命令 1.1 常用命令 1.1.1 浏览目录 cd: 进入某个目录。如：cd /home/songluzhi/vasp/CH4 cd .. 上一层目录；cd / 根目录； ls: 显示目录下的文件。 注：输入目录名时，可只输入前3个字母，按Tab键补全。1.1.2 浏览文件 cat：显示文件内容。如：cat INCAR 如果文件较大，可用：cat INCAR | more (可以按上下键查看) 合并文件：cat A B > C (A和B的内容合并，A在前，B在后) 1.1.3 目录操作 mkdir：建立目录；rmdir：删除目录。 如：mkdir T-CH3-Rh111 1.1.4 文件操作 rm：删除文件；vi：编辑文件；cp：拷贝文件 mv：移动文件；pwd：显示当前路径。 如：rm INCAR rm a* (删除以a开头的所有文件) rm -rf abc (强制删除文件abc) tar：解压缩文件。压缩文件？？rar 1.1.5 系统信息 df：分区占用大小。如：df -h du：各级目录的大小。 top：运行的任务。 ps ax：查看详细任务。 kill：杀死任务。如：kill 12058 （杀死PID为12058的任务）注：PID为top命令的第一列数字。
初学VASP最重要的INCAR参数
初学VASP(六) 最重要的INCAR参数 初学VASP(六) 最重要的INCAR参数 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)
VASP-INCAR参数设置
1. 结构优化(Opt) SYSTEM = opt ISTART = 0 INIWAV = 1 ICHARG = 2 ISPIN = 2 LREAL = Auto ENCUT = 400 PREC = high NSW= 600 NELM = 60 IBRION = 2 ISIF = 2 POTIM = 0.1 ALGO= Fast LVDW = .TRUE. EDIFF = 1E-5 EDIFFG = 1E-4 or -0.05 # 体系需计算TS时，全部结构优化EDIFFG均设置为-0.05 ISMEAR = 0 SIGMA = 0.2 LCHARG = .FALSE. LWAVE = .FALSE.
2. 过渡态搜索(TS): 计算时先进行低精度计算，再进行高精度计算 SYSTEM= TS ISTART = 0 INIWAV = 1 ICHARG = 2 ISPIN = 2 LREAL = Auto ENCUT = 400 PREC = high NSW = 600 NELMIN = 6 IBRION = 3 or 1 # 过渡态计算低精度为3，高精度为1 ISIF = 2 POTIM = 0.01 ALGO = Fast LVDW = .TRUE. EDIFF = 1E-5 EDIFFG = -1 or -0.05 # 过渡态计算低精度为-1，高精度为-0.05 ISMEAR = 0 SIGMA = 0.05 LCHARG= .FALSE. LWAVE= .FALSE. IMAGES=8 # TS专属设置 SPRING=-5 # TS专属设置 LCLIMB=.TRUE. # TS专属设置
vasp参数介绍
内容 描述体系：SYSTEM 设置如何输入或构造初始的电荷密度和波函数：ISTART, ICHARG, INIWA V 设置电子的优化： 平面波切断动能和缀加电荷时的切断值：ENCUT, ENAUG 电子部分优化的方法：ALGO, IALGO, LDIAG 电荷密度混合的方法：IMIX, AMIX, AMIN, BMIX, AMIX_MAG, BMIX_MAG, WC, INIMIX, MIXPRE, MAXMIX 自洽迭代步数和收敛标准：NELM, NELMIN, NELMDL, EDIFF 设置原子的驰豫： 原子如何移动以及步长和步数：IBRION, NFREE, POTIM, NSW 分子动力学相关参数：SMASS, TEBEG, TEEND, POMASS,NBLOCK, KBLOCK, PSTRESS 原子驰豫收敛标准:EDIFFG 定义态密度积分的方法和参数: smearing方法和参数：ISMEAR, SIGMA 计算态密度时能量范围和点数:EMIN, EMAX, NEDOS 计算分波态密度的参数：RWIGS, LORBIT 其他： 计算精度控制:PREC 磁性计算:ISPIN, MAGMOM, NUPDOWN 交换关联函数:GGA, VOSKOWN 计算ELF和总的局域势:LELF, LVTOT 结构优化参数:ISIF 初始化 SYSTEM: 注释所计算的体系，以示说明。 NWRITE: 默认值为2，可赋予值为0| 1 | 2 |3 |4。决定OUTCAR中输出内容的详细程度ISTART: 如果计算目录中有W A VECAR文件，则默认值为1，否则为0。可赋予值为0| 1 | 2 | 3 。决定是否读入W A VECAR： 0: 开始新的计算，按INIWA V初始化波函数 1: 接着计算，通常用在测试ENCUT的收敛性以及计算结合能曲线(也就是体积和总能的关系) 2: 接着计算，通常用在希望保持基矢不变的计算中 3: 接着计算，读入上一次计算得到的电荷密度和波函数，不推荐用 ICHARG: 如果ISTART=0，则默认值为2，否则为0。可赋予值为0 | 1 | 2。决定了如何构造初始的电荷密度 0: 从初始的波函数构造 1：读入CHGCAR读入，并同原子密度进行线性插值 2: 构造原子密度 11：读入自洽的CHGCAR，并进行能带计算或态密度的非自洽计算 12: 非自洽的原子密度计算 电子自洽收敛 INIWAV: 默认值为1，可赋予值为0 | 1。只在开始新的计算(也就是ISTART=0)中有效。决定了如何初始化波函数 0: 采用”jellium 波函数”
VASP全参数设置详解
VASP参数设置详解软件主要功能： 采用周期性边界条件（或超原胞模型）处理原子、分子、团簇、纳米线（或管）、薄膜、晶体、准晶和无定性材料，以及表面体系和固体 l 计算材料的结构参数（键长、键角、晶格常数、原子位置等）和构型 l 计算材料的状态方程和力学性质（体弹性模量和弹性常数） l 计算材料的电子结构（能级、电荷密度分布、能带、电子态密度和ELF） l 计算材料的光学性质 l 计算材料的磁学性质 l 计算材料的晶格动力学性质（声子谱等） l 表面体系的模拟（重构、表面态和STM模拟） l 从头分子动力学模拟 l 计算材料的激发态（GW准粒子修正） 计算主要的四个参数文件：INCAR ,POSCAR,POTCAR ,KPOINTS,下面简要介绍，详细权威的请参照手册 INCAR文件： 该文件控制VASP进行何种性质的计算，并设置了计算方法中一些重要的参数，这些参数主要包括以下几类： 对所计算的体系进行注释：SYSTEM 定义如何输入或构造初始的电荷密度和波函数：ISTART，ICHARG，INIWAV 定义电子的优化 –平面波切断动能和缀加电荷时的切断值：ENCUT，ENAUG –电子部分优化的方法：ALGO，IALGO，LDIAG –电荷密度混合的方法：IMIX，AMIX，AMIN，BMIX，AMIX_MAG，BMIX_MAG，
WC，INIMIX，MIXPRE，MAXMIX –自洽迭代步数和收敛标准：NELM，NELMIN，NELMDL，EDIFF 定义离子或原子的优化 –原子位置优化的方法、移动的步长和步数：IBRION，NFREE，POTIM，NSW –分子动力学相关参数：SMASS，TEBEG，TEEND，POMASS，NBLOCK，KBLOCK，PSTRESS –离子弛豫收敛标准：EDIFFG 定义态密度积分的方法和参数 –smearing方法和参数：ISMEAR，SIGMA –计算态密度时能量范围和点数：EMIN，EMAX，NEDOS –计算分波态密度的参数：RWIGS，LORBIT 其它 –计算精度控制：PREC –磁性计算：ISPIN，MAGMOM，NUPDOWN –交换关联函数：GGA，VOSKOWN –计算ELF和总的局域势：LELF，LVTOT –结构优化参数：ISIF –等等。 主要参数说明如下： SYSTEM：该输入文件所要执行的任务的名字。取值：字符串，缺省值：SYSTEM NWRITE：输出内容详细程度。取值：0~4，缺省值：2 如果是做长时间动力学计算的话，最好选0或1(首末步/每步核运动输出)，短时运算用2，选3则会在出错的时候给出说明信息。
vasp在计算磁性的实例和讨论
兄弟，问3个问题 1，vasp在计算磁性的时候，oszicar中得到的磁矩和outcar中得到各原子磁矩之和不一致，在投稿的是否曾碰到有审稿人质疑，对于这个不一致你们一般是怎么解释的了？ 2，另外，磁性计算应该比较负责。你应该还使用别的程序计算过磁性，与vasp结果比较是否一致，对磁性计算采用的程序有什么推荐。 ps：由于曾使用vasp和dmol算过非周期体系磁性，结构对磁性影响非常大，因此使用这两个程序计算的磁性要一致很麻烦。还不敢确定到底是哪个程序可能不可靠。 3，如果采用vasp计算磁性，对采用的方法和设置有什么推荐。 1，OSZICAR中得到的磁矩是OUTCAR中最后一步得到的总磁矩是相等的。总磁矩和各原子的磁矩(RMT球内的磁矩)之和之差就是间隙区的磁矩。因为有间隙区存在,不一致是正常的。 2，如果算磁性，全电子的结果更精确，我的一些计算结果显示磁性原子对在最近邻的位置时，PAW与FPLAW给出的能量差不一致，在长程时符合的很好。虽然并没有改变定性结论。感觉PAW似乎不能很好地描述较强耦合。我试图在找出原因，主要使用exciting和vasp做比较。计算磁性推荐使用FP-LAPW, FP-LMTO, FPLO很吸引人(不过是商业的)，后者是O(N)算法。 3，使用vasp计算磁性，注意不同的初始磁矩是否收敛为同一个磁矩。倒没有特别要注意的地方，个人认为。 归根结底，需要一个优秀的交换关联形式出现 VASP计算是否也是像计算DOS和能带一样要进行三步（结构优化，静态自洽计算，非自洽计算），然后看最后一步的出的磁矩呢？ 一直想计算固体中某个原子的磁矩，根据OUTCAR的结果似乎不能分析，因为它里面总磁矩跟OSZICAR的值有一定的差别，据说是OUTCAR中只考虑WS半径内磁矩造成的。最近看到一个帖子说是可以用bader电荷分析方法分析原子磁矩。如法炮制之后发现给出的总磁矩与OSZICAR的结果符合的甚好，可是觉得没有根据，有谁知道这样做的依据吗，欢迎讨论！ 设置ISPIN=2计算得到的态密度成为自旋态密度。 设置ISPIN=2就可以计算磁性，铁磁和反铁磁在MAG里设置。最后得到的DOS是分up和down的。 磁性计算 (2006-12-03 21:02) 标签： - 分类：Vasp ·磁性计算
VASP参数设置详解(精)
VASP参数设置详解 计算材料 2010-11-30 20:11:32 阅读197 评论0 字号：大中小订阅 转自小木虫，略有增减 软件主要功能： 采用周期性边界条件（或超原胞模型）处理原子、分子、团簇、纳米线（或管）、薄膜、晶体、准晶和无定性材料，以及表面体系和固体 l 计算材料的结构参数（键长、键角、晶格常数、原子位置等）和构型 l 计算材料的状态方程和力学性质（体弹性模量和弹性常数） l 计算材料的电子结构（能级、电荷密度分布、能带、电子态密度和ELF） l 计算材料的光学性质 l 计算材料的磁学性质 l 计算材料的晶格动力学性质（声子谱等） l 表面体系的模拟（重构、表面态和STM模拟） l 从头分子动力学模拟 l 计算材料的激发态（GW准粒子修正） 计算主要的四个参数文件：INCAR ,POSCAR,POTCAR ,KPOINTS,下面简要介绍，详细权威的请参照手册 INCAR文件： 该文件控制VASP进行何种性质的计算，并设置了计算方法中一些重要的参数，这些参数主要包括以下几类： l对所计算的体系进行注释：SYSTEM
l定义如何输入或构造初始的电荷密度和波函数：ISTART，ICHARG，INIWAV l定义电子的优化 –平面波切断动能和缀加电荷时的切断值：ENCUT，ENAUG –电子部分优化的方法：ALGO，IALGO，LDIAG –电荷密度混合的方法：IMIX，AMIX，AMIN，BMIX，AMIX_MAG，BMIX_MAG，WC，INIMIX，MIXPRE，MAXMIX –自洽迭代步数和收敛标准：NELM，NELMIN，NELMDL，EDIFF l定义离子或原子的优化 –原子位置优化的方法、移动的步长和步数：IBRION，NFREE，POTIM，NSW –分子动力学相关参数：SMASS，TEBEG，TEEND，POMASS，NBLOCK，KBLOCK，PSTRESS –离子弛豫收敛标准：EDIFFG l定义态密度积分的方法和参数 –smearing方法和参数：ISMEAR，SIGMA –计算态密度时能量范围和点数：EMIN，EMAX，NEDOS –计算分波态密度的参数：RWIGS，LORBIT l其它 –计算精度控制：PREC –磁性计算：ISPIN，MAGMOM，NUPDOWN –交换关联函数：GGA，VOSKOWN –计算ELF和总的局域势：LELF，LVTOT –结构优化参数：ISIF –等等。 主要参数说明如下： ?SYSTEM：该输入文件所要执行的任务的名字。取值：字符串，缺省值：SYSTEM ?NWRITE：输出内容详细程度。取值：0~4，缺省值：2
vasp常见问题
最近在学vasp，这篇文章是百度文库找到的，看了不错，转载一把。另外附上vasp程序，linux中下载后无须安装即可使用。单机中可能会出现内存溢出问题，可以放机群上使用。 01、第一原理计算的一些心得 （1）第一性原理其实是包括基于密度泛函的从头算和基于Hartree-Fock自洽计算的从头算，前者以电子密度作为基本变量（霍亨伯格-科洪定理），通过求解Kohn-Sham方程，迭代自洽得到体系的基态电子密度，然后求体系的基态性质；后者则通过自洽求解Hartree-Fock方程，获得体系的波函数，求基态性质；评述：K-S方程的计算水平达到了H-F 水平，同时还考虑了电子间的交换关联作用。 （2）关于DFT中密度泛函的Functional，其实是交换关联泛函，包括LDA，GGA，杂化泛函等等一般LDA为局域密度近似，在空间某点用均匀电子气密度作为交换关联泛函的唯一变量，多数为参数化的CA-PZ方案；GGA为广义梯度近似，不仅将电子密度作为交换关联泛函的变量，也考虑了密度的梯度为变量，包括PBE,PW,RPBE等方案，BL YP泛函也属于GGA；此外还有一些杂化泛函，B3L YP等。（3）关于赝势在处理计算体系中原子的电子态时，有两种方法，一种是考虑所有电子，叫做全电子法，比如WIEN2K中的FLAPW方法(线性缀加平面波)；此外还有一种方法是只考虑价电子，而把芯电子和原子核构成离子实放在一起考虑，即赝势法。 一般赝势法是选取一个截断半径，截断半径以内，波函数变化较平滑，和真实的不同，截断半径以外则和真实情况相同，而且赝势法得到的能量本征值和全电子法应该相同。赝势包括模守恒和超软，模守恒较硬，一般需要较大的截断能，超软势则可以用较小的截断能即可。另外，模守恒势的散射特性和全电子相同，因此一般红外，拉曼等光谱的计算需要用模守恒势。 赝势的测试标准应是赝势与全电子法计算结果的匹配度，而不是赝势与实验结果的匹配度，因为和实验结果的匹配可能是偶然的。（4）关于收敛测试（a）Ecut，也就是截断能，一般情况下，总能相对于不同Ecut做计算，当Ecut增大时总能变化不明显了即可；然而，在需要考虑体系应力时，还需对应力进行收敛测试，而且应力相对于Ecut的收敛要比总能更为苛刻，也就是某个截断能下总能已经收敛了，但应力未必收敛。（b）K-point，即K网格，一般金属需要较大的K网格，采用超晶胞时可以选用相对较小的K网格，但实际上还是要经过测试。（5）关于磁性一般何时考虑自旋呢？举例子，例如BaTiO3中，Ba、Ti和O分别为+2，+4和-2价，离子全部为各个轨道满壳层的结构，就不必考虑自旋了；对于BaMnO3中，由于Mn+3价时d轨道还有电子，但未满，因此需考虑Mn的自旋，至于Ba和O则不必考虑。其实设定自旋就是给定一个原子磁矩的初始值，只在刚开始计算时作为初始值使用，具体的可参照磁性物理。（6）关于几何优化包括很多种了，比如晶格常数和原子位置同时优化，只优化原子位置，只优化晶格常数，还有晶格常数和原子位置分开优化等等。在PRL一篇文章中见到过只优化原子位置，晶格常数用实验值的例子（PRL 100, 186402 (2008)）；也见到过晶格常数先优化，之后固定晶格常数优化原子位置的情况；更多的情况则是Full geometry optimization。 一般情况下，也有不优化几何结构直接计算电子结构的，但是对于缺陷形成能的计算则往往要优化。（7）关于软件软件大致分为基于平面波的软件，如CASTEP、PWSCF 和ABINIT等等，计算量大概和体系原子数目的三次方相关；还有基于原子轨道线性组合的软件(LCAO)，比如openmx，siesta，dmol等，计算量和体系原子数目相关，一般可模拟较多原子数目的体系。 V ASP是使用赝势和平面波基组，进行从头量子力学分子动力学计算的软件包，它基于CASTEP1989版开发。V AMP/V ASP中的方法基于有限温度下的局域密度近似（用自由能作
VASP使用总结
VASP计算的理论及实践总结 一、赝势的选取 二、收敛测试 1、VASP测试截断能和K 点 2、MS测试 三、结构弛豫 四、VASP的使用流程（计算性质） 1、VASP的四个输入文件的设置 2、输出文件的查看及指令 3、计算单电能 (1) 测试截断能 (2) 测试K点 4、进行结构优化 5、计算弹性常数 6、一些常用指令
一、赝势的选取 VASP赝势库中分为：PP和PAW两种势，PP又分为SP（标准）和USPP（超软）。 交换关联函数分为：LDA（局域密度近似）和GGA（广义梯度近似）。GGA 又分为PW91和PBE。 在VASP中，其中pot ,pot-gga是属于超软势（使用较少）。Paw, paw-pbe ,和paw-gga是属于PAW。采用较多的是PAW-pbe 和PAW-gga。 此外vasp 中的赝势分为几种，包扩标准赝势(没有下标的)、还有硬（harder）赝势(_h)、软（softer）赝势(_s), 所谓的硬（难以赝化），就是指该元素原子的截断动能比较大，假想的势能与实际比较接近，计算得到的结果准确，但比较耗时，难以收敛。软（容易赝化），表示该元素原子的截断动能比较小，赝势模型比较粗糙，但相对简单，可以使计算很快收敛（比如VASP开发的超软赝势）。即硬的赝势精度高，但计算耗时。软的精度低，容易收敛，但节省计算时间。 另一种情况：如Gd_3，这是把f电子放入核内处理，对于Gd来说，f电子恰好半满。所以把f电子作为价电子处理的赝势还是蛮好的（类似还有Lu，全满）。（相对其他的4f元素来说，至于把f电子作为芯内处理，是以前对4f元素的通用做法。计算结果挺好） 常用的做法是：用两种赝势测试一下对自己所关心的问题的影响情况。在影响不大的情况下，选用不含4f电子的赝势（即后缀是3），一来减少计算量，二来避免DFT对4f电子的处理。 【1.赝势的选择： vasp的赝势文件放在目录~/vasp/potentials 下，可以看到该目录又包含五个子目录pot pot_GGA potpaw potpaw_GGA potpaw_PBE ，其中每一个子目录对应一种赝势形式。
VASP参数
SYSTERM=hafnium oxide ISTART=0 ICHARG=2 NWRITE=2 LWAVE=.False. PREC=high???? 默认值Medium在4版本，在5版本normal specify high --mannual 46 #LVTOT=.True. #IDIPOL=4 ALGO=Fast LREAL=Auto #LELF=.TRUE. LCHARG=.FALSE. ################electronic relaxation############ NELM=80 ENCUT=500????? EDIFF=1E-05 #ISPIN=2 ###################ion relaxtion################ NSW= IBRION=2 ISIF=2 EDIFFG=-0.005 #################################################### ISMEAR=0 SIGMA=0.05 #LORBIT=11 #NEDOS= #NBANDS = ISYM=0 #NELECT ###################bader charge analysis#################### #LAECHG=.Ture. #NGXF #NGYF #NGZF ###################NEB######################################## #ICHAIN= #IMAGES= #SPRING= #LCLIMB=
#POTIM= #IOPT= vasp incar 主要参数设置 SYSTEM=B N －－－－－－作业说明语句,即任务的名字 IALGO=38 －－－－－－指定电子自洽计算的算法 38 Davidson 48 RMM-DIIS 算法 NELM=200 －－－－－－自洽迭代次数，一般默认40次 NELMIN=4 －－－－－－最小迭代次数，表面或动力学计算是应增大默认为2次 #NELMDL=-12 －－－－－－开始几步不自洽处理，默认IALGO=8时为－5，IALGO=48时为－12，其他0 #ISTART=0 －－－－－－是否使用已有波函数。有WAVECAR是默认为1使用，否则为0 新的作业 #ALGO=very_fast －－－－同IALGO EDIFF=1E-05 －－－－总能计算允许误差，默认1E-04 ISMEAR=0 －－－－－－确定电子的部分占据数，半导体一般取小于0，金属大于1 LREAL=Auto －－－－－－投影操作是在实空间还是在倒空间，.FALSE.倒，.TRUE.实。AUTO实优化ISYM=0 －－－－－－是否具有对称性，默认1，具有，0 不具有。 #ICHARG=11 －－－－－－如何产生初始电荷密度，默认ISTART＝0时取2，LCAO。否则0，从初始波函 #LDIAG = .FALSE. －－－－－是否执行子空间对角化。IALGO=8时为.TRUE. NELECT －－－－－－总电子数，如果系统不是电中性的就必须设置，所带电荷作为均一的背景电子气考虑。 # ! plot charge density #LPARD = .TRUE. --------是否计算部分电荷。 #IBAND = 176 177 －－－－指定要计算的带 #KPUSE = 1 7 －－－－－要考虑的K点 #LSEPB = .TRUE. #LSEPK = .TRUE. # ! change configuration #FERWE = 159*1.0 1.0 0.0 80*0.0－－电子占据数从INCAR中读入并在计算中不变 #FERDO = 159*1.0 0.0 1.0 80*0.0－－电子自旋极化计算（按每个K点每个能带排列） #SIGMA=0.2 －－－－－－－－－－－展开宽度，默认ISMEAR=1 SIGMA=0.2 #ISPIN=2 －－－－－－－－－－－－指定电子体系的总自旋。 # ! optimize IBRION=2 －－－－－－离子是否运动，－1不运动但做NSW外循环。0动力学模拟，1准牛顿法离子弛豫－－－－－－2 CG法离子弛豫，3 采用衰减二阶运动方程离子弛豫， ISIF=2 －－－－－是否离子弛豫及是否改变原胞。 NSW=200 －－－－－离子运动的步数，默认0