发布时间 : 星期一 文章VASP参数设置详解更新完毕开始阅读5983f9671ed9ad51f01df2b2
? LCORR：决定是否对非完全自洽计算中的力进行Harris修正。取
值：.TRUE./.FALSE.，缺省值：.TRUE.
? TEBEG, TEEND：控制从头分子动力学计算中的起始温度和最终温度（MD有效）。
取值：实数，缺省值：TEBEG = 0 TEEND = TEBEG
注意VASP的温度定义与实际温度有细微的差别，所以TEBEG=T×(N-1)/N，T为实际要求的温度，N为原子数。
? SMASS：控制从头MD中的速度。取值：-3~0，缺省值：-3微正则系综（总自由
能不变），-2保持初速度不变，-1每NBLOCK步调整速度，来保证动能连续，>=0 Nosé算法模拟正则系综
? NPACO：成对相关函数的槽数。取值：整数，缺省值：256
? APACO：成对相关函数求值中的最大距离（?）。取值：整数，缺省值：16
简单说就是在不超过APACO的NPACO个距离上求成对相关函数PCF。
? RWIGS：给出Wigner-Seitz半径，DOS计算用。取值：实数数组，缺省值：从POTCAR
文件中读取
? NELECT：总电子数，如果系统不是电中性的就必须设置，所带电荷作为均一的背
景电子气考虑。取值：实数，缺省值：-（价电子数），由POSCAR和POTCAR文件自动决定（通常不必给出）。
? NUPDOWN上下自旋成分间的电子数之差。取值：整数，缺省值：未设置（此时
将进行完全弛豫）
? EMIN, EMAX：DOS求值的最小/最大能量。取值：实数，缺省值：EMIN = - (lowest
KS-eigenvalue -Δ)，EMAX = - (highest KS-eigenvalue -Δ)
? ISMEAR：决定每个波函数的部分占位fnk如何设置。取值：-5 | -4 | -3 | -2 | 0 | N（-5-带有Blochl修正的四面体方法, -4-不带Blochl修正的四面体方法, -3-根据INCAR文件中提供的smearing参数执行循环，-2-从WAVECAR文件中读取,
-1-Fermi-smearing, 0-Gaussian smearing, >0-method of Methfessel-Paxton orderN，缺省值：1
采用部分占位波函数，用一个函数来平滑积分，尤其是对于金属体系可减少k点。
? SIGMA：决定smearing的宽度（eV）。取值：实数，缺省值：0.2
? ALGO：指定电子最小化算法。取值：Normal(blocked Davidson block iteration
scheme)/VeryFast(RMM-DIIS)/Fast (前两个算法的混合)/All(波函数的所有带同时更新)/ Damped(damped velocity friction algorithm)，缺省值: Normal
? IALGO：指定主算法（整数选择算法）。取值：8(共轭梯度算法)/38(Davidson block
iteration scheme)/48(RMM-DIIS)，缺省值：38
算法是最重要的参数之一。一般VASP推荐使用的是以上三种算法，一般来说8/38是初期比较快收敛，在接近平衡时采用48较快，在初期或MD时使用48可能会遇到不收敛的情况。也可以使用ALGO参数来替代IALGO，设置Fast，VASP会先用38，再自动切换到48。各种算法只要收敛，结果应该一致。
另一个可能有用的选项是-1。不进行实际的计算，只对重要的步骤做计算测试，并将测试得到的各部分耗时输出到OUTPUT中。
? VOSKOWN：决定是否使用VWN插值算法。取值：0(不使用)/1(使用)，缺省值：
0
如果使用了PW91泛函或需要计算磁性质时可以设为1使用。
? Mixing-tags：
? IMIX：混合的类型，取值：整数，缺省值：4
? AMIX：线性混合参数。取值：实数，缺省值：0.8（US-PP），0.4（PAW） ? AMIN：最小混合参数。取值：实数，缺省值：0.1
? BMIX：Kerker混合方案的截断波向量。取值：实数，缺省值：1.0 ? AMIX_MAG：磁化过程的线性混合参数。取值：实数，缺省值：1.6 ? BMIX_MAG：磁化过程的Kerker混合方案的截断波向量。取值：实数，缺
省值：1.0
? WC：Broyden混合方案中每步的加权因子。取值：实数，缺省值：1000.0 ? INIMIX：Broyden混合方案中的初始混合类型。取值：整数，缺省值：1 ? MIXPRE：Broyden混合方案中的预处理类型。取值：整数，缺省值：1 ? MAXMIX：Broyden混合器中存储的最大步数。取值：整数，缺省值：-45
值得注意的是，在MD或者弛豫的时候，设置MAXMIX（>0,一般约3倍的电子SC步数）可能会大大减少核运动步数，但同时也会增加对内存的要求。
? LWAVE,LCHARG：决定是否把波函数（或电荷密度）写入外部文件WAVECAR
（或CHGCAR和CHG）中。取值：.TRUE./.FALSE.，缺省值：.TRUE.
? LVTOT：决定是否把总局域势写入外部文件LOCPOT中。取值：.TRUE./.FALSE.，
缺省值：. FALSE.
? LELF：决定是否创建ELFCAR文件。取值：.TRUE./.FALSE.，缺省值：. FALSE.
ELFCAR用于保存ELF(electron localization function)。
? LORBIT：和适当的RWIGS一起决定是否创建PROCAR或PROOUT文件。取值：
0 | 1 | 2 | 5 | 10 | 11 | 12，缺省值：0(创建DOSCAR和PROCAR文件)
? NPAR：用于控制VASP切换到带间并行模式。取值：整数，缺省值：总节点数。
每一个节点计算一个带，可以提高并行效率，减少通讯量，但可能会大幅增加内存的需求。
? NBANDS：给出计算中实际的总能带数。取值：整数，缺省值：NELECT/2+NIONS/2
（非自旋），0.6*NELECT + NMAG（自旋）
NBANDS的取值应使计算中包含相当数量的空带。因为计算需要大量的空带，至少要求1个空带（否则VASP会给出警告）。
NBANDS对于解决内存需求非常重要。一般NBANDS在NELECT/2+NIONS/2以上可以得到较为精确的结果，但如果内存不够就只能减少NBANDS，在牺牲精度和体系大小之间平衡了。
POTCAR文件
赝势文件，最重要的输入文件之一。可以理解为分子力学模拟中的力场文件，但包含的信息更多。VASP将各元素优化的INCAR里的参数也包含在该文件中，作为支持PREC的缺省选择。
通常各元素的POTCAR文件已经包含在VASP软件包中的赝势库里了，用户需要做的是确定自己具体需要哪几种赝势，然后按照POSCAR文件里的顺序，将所选择的POTCAR文件按顺序连接起来就可以了。如以下命令： cat file1 file2 file3 > POTCAR
赝势库中的赝势文件可以分为以下几类：
? 根据产生方法的不同有Ultra-soft赝势（USPP）和投影扩充波赝势（PAW）； ? 根据交换关联函数的不同有LDA（local density approximation）和GGA（generalized
gradient approximation，又可以再分为PW91和PBE）；
? 根据半芯态的不同元素A的赝势文件还可以分为A，A_sv和A_pv；
? 根据选取的不同截断能量（ENMAX）元素A的赝势文件还可以分为A，A_s和
A_h。