博文
VASP INCAR 参数设置简要说明
||
SYSTEM #系统名称
NWRITE #OUTCAR 输出详细程度 越高越详细，defalut=2
ENCUT #平面波基组能量最大值，默认从POTCAR 中选择最大的原子，对于不同的k点，有着不同的平面波基组个数。
ENAUG #
PREC #影响ENCUT、NGX、NGXF、ROPT，Low、Medium、High、Normal、Accurate，对于ENCUT：Low和Medium选择POTCAR 中的ENMIN，N、A选择ENMAX，HIgh选择1.30*ENMAX； FFT格子：High和Accurate通过包含两倍于Gmax的向量来避免wrap-around error，augmentation charges； ROPT real space projection???
ISPIN #1-自旋关闭 2-自旋开启
MAGMON #
ISTART #0|1|2|3 0-新工作，从INIWAV设定的方式初始化波函数 1-从WAVECAR读取参数，根据新的shape和ENCUT来设置平面波基组。 2-从WAVECAR读取波函数，但使用相同的基组，体积改变时会产生Pulay stress，但是体积形状驰豫的时候应选2 3-在2的基础上从TMPCAR读取前两部离子位置和波函数并继续
ICHARG #0|1|2|11|12 0-从初始波函数计算 1-从CHGCAR读取并插值到新的坐标 2-原子电荷密度叠加 11-电荷密度在非自洽计算中不改变，从CHGCAR 中读取并获得特征值和DOS 12-原子电荷密度叠加的非自洽计算，用于Harris-Foulkes 泛函
INIWAV #0|1 在ISTART=0的时候设置初始波函数 0-自由电子气最低动能的特征函数 1-随机数 Default-1
NELM #电子自洽步数上限， default=60
NELMIN #电子自洽的最少步数 default=2
NELMDL #初始非自洽步数——初始波函数离收敛值很远，需要经过多步非自洽计算之后，再开始更新电荷密度和哈密顿量。取正值为每一步离子驰豫后均进行非自洽计算，取负值意味着只在开始计算时进行非自洽。
EDIFF #电子自洽循环的终止条件，结束循环时体系能量与前一步的差值。 default=10^-4
EDIFFG #离子驰豫的终止条件，取正值为能量收敛阈值，default=10*EDIFF，取负值为原子力收敛阈值
NSW #离子驰豫的最大步数 default=0
NBLOCK&KBLOCK：离子驰豫多少步之后进行电子对相关和DOS计算并写入OUTCAR，默认不需更改
IBRION #决定离子的移动和更新方式 -1-不移动离子，但NSW外循环仍然进行 0-标准分子动力学 1-准牛顿方法RMM-DIIS 在接近最小值点处收敛快 2-CG法，稳妥的方法 3、5-
POTIM #对于IBRION=0(分子动力学）POTIM是以fs为单位的时间步长 对于IBRION=1，2，3 POTIM是一个系数
ISIF #应力张量和离子位置、晶胞形状的设置 0~7 改变晶胞体积的时候，必须提高ENCUT或者设置PREC=High
PSTRESS #外加压力
IWAVPR # 离子移动位置时候波函数核电荷密度的差值更新方式，不建议更改。0-不进行外插 1-简单差值 2-二阶插值，MD必须使用 3-波函数二阶插值，电荷密度简单差值 11,12
ISYM #0|1|2 对称分析算法选择，PAW默认为2，US-PP默认为1， 0为无对称分析
LCORR #Harris 力修正
TEBEG，TEEND
SMASS #控制MD时的速度、系宗
NPACO，APACO
POMASS，ZVAL
RWIGS
LORBIT #PROCAR文件、spd、投射波函数是否输出
NELECT #电子数目，均匀电子气
NUPDOWN #正反自旋电子数之差
EMIN，EMAX #DOS计算的范围，单位为eV
ISMEAR #占有数f的设置，-1-费米 0-高斯 1~N-N阶MP -2-从WAVECAR或INCAR读入并保持不变 -3-在SMEARINGS指定的参数之间循环 -4-四面体无布洛赫修正 -5-四面体有布洛赫修正 对于半导体、绝缘体考虑-5，有时使用0
SIGMA #smearing 宽度
LREAL #赝势非局域积分Cink的计算方法，是空间或倒空间
ROPT #积分球内点的个数，为每种原子设置一个数值
GGA #顾名思义，GGA方法的选择
VOSKOWN #交换关联势关联部分的插值方法，0为标准，1为V-W-N方法，PW91时应选1
ALGO #电子迭代计算波函数时使用的算法，Normal->IALGO=38 ， Very_Fast->IALGO=48 , Fast->初相使用38，之后使用48
IALGO #普通用户应该设置为38（Davidson blocked iteration scheme）或48（RMM-DIIS），高性能、大体系、并行计算时应该选择48.
LDIAG #是否进行子空间对角化
NSIM #同时优化的能带个数，可能能够提高速度。
MIXING Tags
WEIMIN，EBREAK，DEPER #能带对角化迭代的精细调节
LWAVE，LCHARG #决定是否写入WAVECAR、CHGCAR .TRUE.|.FALSE.
LVTOT #是否写入LOCPOT文件，局域势
LELF #是否写入ELFCAR ，electron localization function
并行化：
可以提供能带并行化和平面波系数并行化，两者推荐同时进行
NPAR #1-只对平面波系数并行化 N-N=nodes，每个节点计算一个band， else-nodes/NPAR个节点计算一个band
LPLANE #实空间数据分布， .TRUE.可以节省带宽，但是恶化负载均衡，当NGZ足够大的时候开启这个选项
松紧带方法：elastic band method
用来计算不同构型之间的势垒
POTCAR、INCAR、KPOINTS必须在文件夹里放好，00,01,02等文件夹里放置POSCAR
IMAGES #images数量，必须能够整出节点数量
SPRING #弹簧系数？ 0-image 只能沿相邻images连线切线的垂直方向移动，iamges之间的距离一阶不变 取负值比如-5的时候，进行NEB方法计算
分波电荷密度：
LPARD #必须先进行预先的VASP计算，获得已经收敛的波函数
IBAND #指定要计算的band标号
EINT #指定要计算的能态的能量范围，应设置两个参数，若只设置一个，则第二个默认为EF
NBMOD #进行分波电荷密度计算的能态个数，若已经设置IBAND，NBMOD将自动设置，若为0则是计算全部band的电荷密度 -1|-2|-3
https://m.sciencenet.cn/blog-708035-708093.html
上一篇：作为本征态和C代码的中文
下一篇：深度优先读文献
扫一扫，分享此博文
全部作者的其他最新博文
- • 深度优先读文献
- • 作为本征态和C代码的中文
- • 分析力学笔记
- • 蹲墙角的薛定谔方程
- • 数理方程与复变函数闲谈