得到分波电荷密度的详细过程

第一步：计算做一次自洽计算〔静态计算〕，得到收敛的CHGCAR、CHG、WAVECAR〔注：不计算分波的话，第一步得到的CHGCAR，就可以经过处理得到电荷密度图，需要自旋极化的情况的话，就在计算时INCAR中参加ISPIN=2.〕在自洽计算的根底上，做一次非自洽计算，读入文件：第一步得到的：CHGCAR、CHG、WAVECARPOSCAR、POTCAR、KPOINTS、INCAR但是INCAR需要修改：在INCAR中参加参数：LPARD=.TURE.!表示计算分波电荷密度，具体分以下三种情况来控制需要计算的分波情况：IBAND=20212223！表示计算这四条能带对应的分波电荷密度NBMOD=4！这里代表能带数，可以不手动设置EINT=-0.10.1！表示计算〔-0.1，0.1〕能量区间的分波电荷密度图NBMOD=-2！可以不手动设置EINT=1！表示计算1~费米能级能量区间的分波电荷密度图NBMOD=-3KPUSE=1234！计算这四个K点对应的分波电荷密度以上，1〕和3〕可以同时使用，而2〕一般单独使用。此外，还有两个参数控制输出的情况：LSEPB：指定是不是要把计算的partialchargedensity按每个带分别写到各自对应的文件PARCHG.nb.(设置为.TRUE.)中，还是把它们合并写到一个文件中(相当于把各个带对应的partialchargedensity加起来，设置为.FALSE.)。默认值为.FALSE.。

LSEPK：指定是不是把要计算的partialchargedensity按每个k分别写到各自对应的文件PARCHG.nk.(设置为.TRUE.)中，还是把它们合并写到一个文件中(相当于把各个带对应的partialchargedensity加起来，设置为.FALSE.)。默认值为.FALSE.。第二步：处理数据1〕将计算结果中的POSCAR和CHGCAR下载到本地，分别用小程序处理。POSCAR---POStoOKes11.exe〔输入版本、原子种类数和原子符号〕---得到MS可以翻开的结构文件.res文件CHGCAR〔删除原子符号那一行〕---chgcar2grd111.exe(输入版本，是否自旋极化)---可以得到.grd文件。2〕将结构文件导入MS导入电荷密度的数据，Edit—Insertform—找到文件位置，翻开—OK！看上去没什么变化，需要点击鼠标右键DisplayStyle—field—勾选Visible电荷密度，就在图中显示出来，但还需要根据自己的需要调整显示的图片：1〕右键--ColorMaps—Bands可以调整图片的细腻程度—From…to可以调整区间的大小，Legend可以在右下角显示图例。2〕切片操作，View—Toolbars—Modules显示切片操作需要的图标得到一个平行于BC面的切面：我们可以对这个切面进行操作：选定这个切面，变成黄色，点击下面的图标AlignLeft/Right：左对齐/右对齐AlignVertical:垂直对齐Alignin/out:对齐轴AlignWithScreen：对齐屏幕AlignWithViewYZPlane：对齐YZ平面AlignHorizontal：