文档介绍：vasp计算总结优选1加U(LMCO为例)LDAU=.=2LDAUL=-122-1LDAUU==#Thesimplified(rotationallyinvariant)approachtotheLSDA+U#specifiesthel-quantumnumberforwhichtheon-siteinteractionisadded(-1=noon-sitetermsadded,1=p,2=d,3=f,Default:LDAUL=2)#specifiestheeffectiveon-siteCoulombinteractionparameters#specifiestheeffectiveon-=1ICHARG=11ISMEAR=-5LORBIT=11用VASP计算DOS准备好KPOINTS文件,增加k点网格将上一步自洽计算得到的CHG、CHGCAR拷贝至同一目录下优选3用VASP计算能带ISTART=1ICHARG=11NBANDS一般可以用默认,有时候可以适当增加ISMARE=-5不可用。半导体或绝缘体用0,金属用1NAME20Line-:,g++编译,,删除前七行grep-vEEIGENVAL>./,,第二列数据是以Ferimlevel为参考的本征值。+./#!/bin/(改为要计算的值)docat>POSCAR<<!(POSCAR文件,要计算的量改为$i)!echo"$i";mpirun-np10vasp>=`grep"TOTEN"OUTCAR|tail-1|awk'{printf"%\n",$5}'`-$imvOUTCAROUTCAR-$imvDOSCARDOSCAR-$iecho$i$E>>commentdone找最低能量的晶格参数C优选7VASP中,用BerryPhase的方法计算极化值INCAR中要添加的参数:LCALCPOL=.=1E-6精度DIPOL=,任意选的,但是不要和离子重合首先,,磁性有影响时要考虑进去。MS查看对称性,找到所对应的群。/#14找到对应的点群铁电体物理,钟维烈著,P651。附录一,30个晶体点群。10个极性点群是非常有可能产生极化的,11个非极性中兴对称群是可能有极化的,11个中心对称点群是没有极化的。点电荷估计:晶格某一方向所有原子坐标*离子价态相加,与优化后相加得到的值对比。优选8Spinorbitalcoupling(SOC)的计算(LMCO为例)ISTART=1要读取WAVCARICHARG=1最好选取1而不用11,对结果影响是比较大的LSORBIT=.=6计算soc的时候一定要有此参数,d电子4,f电子6MAGMOM=6*000400400400-400-400-400300300300300-300-300-300-372*0SAXIS=111磁矩的方向(见vasp手册,有两种设置方法)NBANDS=408能带数是线性计算的二倍ISYM=PAT=.,作为soc计算的输入文件优选9VASP难收敛的问题于一些磁性体系、镧系和锕系元素及相关化合物的静态计算(电子迭代),经常会遇到“难收敛”的问题。1、LMAXMIXDefault:LMAXMIX=2AnadditionalflagcontrolsuptowhichlquantumnumbertheonsitePAWchargedensitiesarepasse