VASP入门到精通[33]分子轨道-partial charge

优化完体系要做一步静态计算，注意加上LWAVE=.TRUE.LCHARG = .TRUE.参数如下：

PREC = AccurateENCUT = 400EDIFF = 1e-5EDIFFG=-0.01IBRION = -1ISIF = 2NSW = 0ISMEAR = 0 LCHARG = .TRUE.LWAVE = .TRUE.NELM = 100LORBIT=11ICHARG=2ISYM = 0

计算完成后，加入如下INCAR参数：ICHARG=1LPARD=.TRUE.LSEPB= .TRUE.EINT= -17.0这里要注意一些问题，参数的意义：比如小编想看PROCAR中能量从-17.0eV以上的能级的partial charge 就可以设置成这种形式，

这样的话，从-17开始以上到费米能级的所有能级都会出来；计算完成后会多出很多CHG文件如下图：

可以看到小编的文件夹里有48个partial charge文件，是因为小编在-17eV以上有48个能级，从PROCAR里面也可以看出来

接下来看一下charge的情况。比如看第一个：PARCHG.0001.ALLK改名字，让vesta能识别，比如可以改成CHGCAR01

结合分子轨道理论，小伙伴可以自己分析性质，是什么杂化形成的轨道

其中还有别的方法看partial charge，比如看某些能级范围的，可以只看HOMO-LUMO的能级轨道图如：LPARD=.TRUE.                                 计算partial charge的开关IBAND= 20 21 22 23                  选择看哪些带，这个可以从PROCAR里确定KPUSE= 1 2 3 4                          选择K点1 2 3 4的带LSEPB=.TRUE.                             分开输出文件，就像小编截图的48个文件LSEPK=.TRUE.                            分开输出k点的带