Va sp的输入输出文件.pptVIP

* * * * * * * * * * EDIFFG=EDIFF×10 当数值为负数时，表示以力作为收敛标准，多数情况均采用 力作为收敛标准。 ALGO=38|48 该关键词确定能量计算迭代方法 38-采用Davidson优化方法；(可靠，但速度慢) 48-采用RMM-DIIS算法；(常用，速度快) ISYM=0|1|2 该关键词确定能量和构型优化时是否使用对称性(将影响到K 点数目和计算量大小) 0-不使用对称性； 1-采用对称性； 2-用于PAW型赝势； (7) (8) NELM=整数 该关键词确定能量自洽场最大迭代轮数，缺省为60轮； NELMIN=整数 在构型优化中，计算每个构象能量时最少迭代轮数，一般为3~4，以保证能量和力的稳定性； 定义DFT泛函类型，注意要与POTCAR中的赝势类型一致。 (9) (10) ISPIN=1|2 1-非自旋极化计算(缺省) 2-自旋极化计算, 将给出体系磁矩大小(对含有过渡金属原 子体系，一般均要采用自旋极化方法)。 (11) (12) (13) ISMEAR选择： 1) 对半导体或绝缘体选取-5，如果单胞较大时，或者所选取k 点数目少时，用0； 2) 对导体，通常用0； SIGMA取值： SIGMA取值的原则是使得计算得到的TS项(OUTCAR中)， 分摊到每个原子上时小于1meV，否则得到的总能量不准确， 对导体尤其要注意该参数的选择。 以下为构型优化所用关键词： NSW=整数 构型优化的最大轮数 IBRION = -1|0|1|2 构型优化方法： -1-构型不变更； 0-分子动力学模拟； 1-采用准牛顿方法确定新的构型(当初始构型较合理时使用)； 2-采用CG方法确定构型(当初始构型离平衡位置较远时使用)。 POTIM=数值 控制构型优化步长，缺省为0.5，对动力学模拟则为时 间步长(单位为fs) (14) 输出控制关键词： LCHARG = .FALSE. (输出电荷密度？) LWAVE = .FALSE. (输出波函数？) LVTOT = .FALSE.(输出静电势，求功函时使用) 其他关键词： NPAR = 8 (CPU数目，并行计算时使用) LPLANE = .TRUE.(与并行算法有关) (15) (16) 实例： SYSTEM = Silicon bulk NPAR = 2 LPLANE = .TRUE. Elecronic minimisation ISTART = 0 LREAL= .FALSE. PREC = Medium precission: Mediun/High/Low EDIFF = 1e-4 converge criterion: default = 1e-4 EDIFFG = -0.02 converge criterion for relation loop IALGO = 48 algorithm (8-CG, 48-RMM) NELMIN = 3 the minimum number of electronic SC steps ISYM = 1 symmetry (2-PAW on, 1-US-PPs on, 0-off) ISIF = 3 Relax ions ISPIN = 1 ISMEAR = -5 (tetrahedron/gaussian/m-p) SIGMA = 0.1 OUTPUT CONTROL LCHARG = .FALSE. LWAVE = .FALSE. the i/o cost is not worth it. LVTOT = .FALSE. IONIC RELAXATION NBLOCK = 1 steps for inner block NSW = 300 number of steps for IOM IBRION = 1 -1:no update 0-MD 1-quasi-New 2-CG POTIM = 0.50 default 0.5 of IBRION=1-3 练习： 对M