《vasp计算参数设置》由会员分享，可在线阅读，更多相关《vasp计算参数设置（10页珍藏版）》请在金锄头文库上搜索。
1、软件主要功能：采用周期性边界条件（或超原胞模型）处理原子、分子、团簇、纳米线（或管） 、薄膜、 晶体、准晶和无定性材料，以及表面体系和固体 l 计算材料的结构参数（键长、键角、晶格常数、原子位置等）和构型 l 计算材料的状态方程和力学性质（体弹性模量和弹性常数） l 计算材料的电子结构（能级、电荷密度分布、能带、电子态密度和 ELF） l 计算材料的光学性质 l 计算材料的磁学性质 l 计算材料的晶格动力学性质（声子谱等） l 表面体系的模拟（重构、表面态和 STM 模拟） l 从头分子动力学模拟 l 计算材料的激发态（GW 准粒子修正） 计算主要的四个参数文件：INCAR ,POSCAR,P
2、OTCAR ,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， MAXM
3、IX 自洽迭代步数和收敛标准：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，MAGMO
4、M，NUPDOWN 交换关联函数：GGA，VOSKOWN 计算 ELF 和总的局域势：LELF，LVTOT 结构优化参数：ISIF 等等。主要参数说明如下：? SYSTEM：该输入文件所要执行的任务的名字。取值：字符串，缺省值：SYSTEM? NWRITE：输出内容详细程度。取值：04，缺省值：2 如果是做长时间动力学计算的话，最好选 0 或 1(首末步/每步核运动输出)，短时运算用 2， 选 3 则会在出错的时候给出说明信息。? ISTART：决定是否读取 WAVECAR 文件。取值：02，缺省 0/1 for 无/有 前次计算 的 WAVECAR（波函数） 0：begin from scr
5、atch，根据 INIWAV 初始化波函数 1：restart with constant energy cut-off，从 WAVECAR 读取波函数（重定义平面波集） 2：restart with constant basis set，从 WAVECAR 读取波函数（平面波集不变）? ICHARG：决定如何建立初始电荷密度。取值：02，缺省值: if ISTART=0 2 else 0 0：由初始波函数计算电荷密度 1：从 CHGCAR 文件读取电荷密度 2：使用原子电荷密度的叠加 +10 非自洽计算? ISPIN：是否进行 spin polarized calculation。取值：1，
6、2（1-no，2-yes） ，缺省值：2? MAGMOM：在 ICHARG=2 或在 CHGCAR 中未包含磁化密度（ICHARG=1）时，指 定每个原子的初始磁化时刻。取值：实数数组，缺省值: 对 ISPIN=2 NIONS*1.0，对非共 线型磁化体系 3*NIONS*1.0? INIWAV 如何设置初始波函数，只在 ISTART=0 时使用。取值：0，1（0-最低动能的 平面波，1-随机数） ，缺省值：1。? IDIPOL 控制计算单极、偶极和四极修正。取值：14。 13 只计算第一/二/三晶矢方向，适于厚板（slab）的计算 4 所有方向都计算，适于计算孤立分子? PREC 进动（pr
7、ecession） 。取值：low/medium/high/normal/ accurate/single)，缺省值: Normal（VASP.4.X） ；Medium（VASP.5.X） VASP4.5+采用了优化的 accurate 来替代 high，所以一般不推荐使用 high。不过 high 可以确 保“绝对收敛”，作为参考值有时也是必要的。同样受推荐的是 normal，适于作为日常计算 选项。 受 PREC 影响的参数有四类： ENCUT；NGX，NGY，NGZ；NGXF，NGYF，NGZF；ROPT。如果设置了 PREC，这些 参数就都不需要出现了，当然直接设置相应的参数也有同样效
8、果。具体影响效果见 p5354。? ENCUT 平面波基组的截断能量（eV） 。取值：实数，缺省值：受 PREC 设置影响，从 POTCAR 文件中找出相应的 ENMAX/ENMIN 值来设置。 PREC = Low Medium Accurate NormalENCUT = ENMIN ENMAX ENMAX ENMAXSingle High ENMAX ENMAX*1.3 对于多个元素不同的 ENMAX/ENMIN，都取最大值。 该参数非常重要，最好不要手工去设置，除非文献告诉你要用多少，或者经过结果可靠性 的验证。当然，为了测试一下提交的任务，也不妨先设个较小的值。? NGX，NGY，N
9、GZ：控制 FFT 网格在三个晶矢方向上的格点数量。 ? NGFX，NGFY，NGFZ：控制第二次更精确的 FFT 网格的格点数量。 也是两类重要的最好不要去动的参数。在未指定的情况下将根据 PREC 的设置从 POTCAR 中自动读取。 PREC=High/Accurate，基组中向量的 2 倍值，用来避免 wrap around errors，得到精确解。 PREC=Low/Medium/Normal，基组中向量的 3/4 倍值（已足够精确到 1 meV/atom） 。? LREAL: 决定投射是在实空间还是倒易空间进行。取值：.TRUE.（实空间） /.FALSE.（倒易空间） ，缺省值
10、：.FALSE. 用于求解赝势的非局域部分用到的一个积分，在倒格空间里采用平面波基组求解，在实空 间里则采用积分球求解。 其他还有两个选项：O or On，A or Auto。 On 和.TRUE.的差别在于是否使用 King-Smith 算法优化，设为 Auto 则进行自动选择，推荐 使用。? ROPT: 在 LREAL=Auto or On 时，优化控制每个核周围的积分球内的格点数。取值： 实数数组 For LREAL=OnPREC=Low, 700 points in the real space sphere (ROPT=0.67)PREC=Med, 1000 points in th
11、e real space sphere (ROPT=1.0)PREC=High, 1500 points in the real space sphere (ROPT=1.5) For LREAL=AutoPREC=Low, accuracy 10-2 (ROPT=0.01)PREC=Med, accuracy 2*10-3 (ROPT=0.002)PREC=High accuracy 2*10-4 (ROPT=2E-4)? NELM, NELMIN and NELMDL：控制电子自洽循环步数。取值：整数 NELM：电子自洽循环最大次数。缺省值：60 NELMIN：电子自洽循环最小次数。缺省值
12、：2 NELMDL：弛豫次数。缺省值：if ISTART=0, INIWAV=1, and IALGO=8，-5，if ISTART=0, INIWAV=1, and IALGO=48，-12，else 0 NELMDL 可以取负值。如果初始波函数采用随机赋值，即 ISTART=0，INIWAV=1，那么 很可能开始的值比较离谱，那么在第一步核运动循环之前采用 NELMDL（负值）步的非自 洽（保留初始的 H）步计算将减少计算所需的时间。? EDIFF：指定电子自洽循环的全局中断条件，用于控制收敛精度。取值：实数，缺省 值：10-4注意，即使 EDIFF=0，NELM 步电子自洽循环也会执行。
13、? EDIFFG：指定离子弛豫循环的中断条件，用于控制核运动的收敛精度。取值，实数， 缺省值：10*EDIFF EDIFFG0 在两个离子步的总自由能之差小于 EDIFFG 时停止 EDIFFG=0 Nos 算法 模拟正则系综? NPACO：成对相关函数的槽数。取值：整数，缺省值：256 ? APACO：成对相关函数求值中的最大距离（?） 。取值：整数，缺省值：16 简单说就是在不超过 APACO 的 NPACO 个距离上求成对相关函数 PCF。? RWIGS：给出 Wigner-Seitz 半径，DOS 计算用。取值：实数数组，缺省值：从 POTCAR 文件中读取? NELECT：总电子数，
14、如果系统不是电中性的就必须设置，所带电荷作为均一的背景 电子气考虑。取值：实数，缺省值：-（价电子数） ，由 POSCAR 和 POTCAR 文件自动决 定（通常不必给出） 。? NUPDOWN 上下自旋成分间的电子数之差。取值：整数，缺省值：未设置（此时将 进行完全弛豫）? EMIN, EMAX：DOS 求值的最小/最大能量。取值：实数，缺省值：EMIN = - (lowest KS-eigenvalue - )，EMAX = - (highest KS-eigenvalue - )? ISMEAR：决定每个波函数的部分占位 fnk 如何设置。取值：-5 | -4 | -3 | -2 | 0
15、 | N（-5- 带有 Blochl 修正的四面体方法, -4-不带 Blochl 修正的四面体方法, -3-根据 INCAR 文件中提供的 smearing 参数执行循环，-2-从 WAVECAR 文件中读取, -1-Fermi-smearing, 0-Gaussian smearing, 0-method of Methfessel-Paxton order N，缺省值：1 采用部分占位波函数，用一个函数来平滑积分，尤其是对于金属体系可减少 k 点。? SIGMA：决定 smearing 的宽度（eV） 。取值：实数，缺省值：0.2? ALGO：指定电子最小化算法。取值：Normal(blocked Davidson block iteration scheme) /VeryFast(RMM-DIIS)
