将优化好的CONTCAR复制为POSCAR进行下面的计算，后面的POSCAR都不会动 自洽计算相关参数： IBRION = -1 !原子不移动 ISIF = 2 !不改变元胞形状和体积 NSW=0 !原子弛豫步数为0
!................................................ ISTART = 0 ICHARG = 2 !0，从初始波函数计算电荷密度，1，从ACHRGCAR文件读入，2，用原子电荷，密度组装 NSW = 0 !不再进行原子弛豫 !................................................. ENCUT = 400 ISMEAR = 0 !ISMEAR=1或2时，是针对于金属而言的，对于绝缘体或半导体来说一般选取的是0或-5,而取-5时候，主要是能量值得差异 SIGMA = 0.05 ISPIN = 2 !没有自旋就注释掉 IBRION = -1 !原子不移动 ISIF = 2 !不改变元胞形状和体积 !................................................... EDIFF = 1E-7 EDIFFG = -0.001 LORBIT = 11 !是否输出投影波函数 POTIM = 0.4 !只影响收敛速度，不影响结果 LREAL=Auto !.FALSE ：倒空间 .TRUE 实空间,一般超过20个原子使用 .ON和Auto都是实空间 !................................................. PREC = Accurate !用来选择ENCUT，Low，取potcar中的最小值，Medium或Accurate取最大值，High取最大值1.3倍 !可能影响运算效率....................................... LPLANE = .TRUE. LSCALU = .FALSE. NPAR=2 !节点数，且有了NPAR，NCORE就不必要了，两个乘起来等于总核数 !可能影响运算效率....................................... ALGO = Fast !简化电子优化算法的开关，Normal对应38，Very_Fast对应48，Fast更稳定 ISYM = 2 !是否具有对称性，1具有，0不具有，2内存对电荷密度的处理更有效 NELM=250 !自洽中最多迭代的次数 !....................................................... LSORBIT = T ! soc GGA_COMPAT = F ! !!!!! SAXIS = 0 0 1 ! !!!!! !....................................................... LMAXMIX = 4 !.........有d电子以后的轨道需要加U LDAU=T LDAUTYPE=2 LDAUL= 2 -1 -1 !2时对d加u，3是对f加u LDAUU= 2.7 0.0 0.0 LDAUJ= 0.0 0.0 0.0 LASPH=T 常见错误： （1）自洽无法收敛：减小ENCUT （2）自洽计算时提交上去后没有反应，且不报错，也不生成文件，去掉u，减小encut （3）自洽出现 VERY BAD NEWS! internal error in subroutine SGRCON: Found some non-integer element in rotation matrix 注意对称性有没有开 （4）自洽不收敛，加上以下几行： AMIX = 0.02 BMIX = 0.0001 AMIX_MAG = 0.8 BMIX_MAG = 0.0001 注意vasp6计算不能加，否则出现奇怪的错误！！！！！！！！！！！！！ 出现Error EDDDAV: Call to ZHEGV failed. 可能是阶段能ENCUTT太小 如果达到精度，依然计算不停止，是因为ENCUT太大
No Comments
Leave a comment Cancel