文档介绍：vasp计算总结1PPT课件加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拷贝至同一目录下3PPT课件用VASP计算能带ISTART=1ICHARG=11NBANDS一般可以用默认,有时候可以适当增加ISMARE=-5不可用。半导体或绝缘体用0,金属用1NAME20Line-:,g++编译,,删除前七行grep-vEEIGENVAL>./,,第二列数据是以Ferimlevel为参考的本征值。+---------011------,第四个数是相应k点的权重。第一行,前三个整数无意义,第四个整数,如果是2,表示是自旋极化的计算,如果是1,表示非自旋极化的计算。第2至5行的数据含义不大明确,可以不管它。第6行的数据表示:第一个数表示体系总的价电子数目,第二个数表示的计算能带时总的k点数目,第三个数表示的是计算能带时计算了多少条能带。第9行给出的是该k点对应的本征值的序号(即第几条能带),及相应的本征值。+./#!/bin/(改为要计算的值)docat>POSCAR<<!(POSCAR文件,要计算的量改为$i)!echo"$i";mpirun-np10vasp>=`grep"TOTEN"OUTCAR|tail-1|awk'{printf"%\n",$5}'`-$imvOUTCAROUTCAR-$imvDOSCARDOSCAR-$iecho$i$E>>commentdone找最低能量的晶格参数C7PPT课件VASP中,用BerryPhase的方法计算极化值INCAR中要添加的参数:LCALCPOL=.=1E-6精度DIPOL=,任意选的,但是不要和离子重合首先,,磁性有影响时要考虑进去。MS查看对称性,找到所对应的群。/#14找到对应的点群铁电体物理,钟维烈著,P651。附录一,30个晶体点群。10个极性点群是非常有可能产生极化的,11个非极性中兴对称群是可能有极化的,11个中心对称点群是没有极化的。点电荷估计:晶格某一方向所有原子坐标*离子价态相加,与优化后相加得到的值对比。8PPT课件Spinorbitalcoupling(SO