一、Vasp 部分: 1. 在计算 band 的时候,INCAR 中要设置 LORBIT = 11 这个参数,如图 1。 图1 2. 在计算 band 结束后,把 vasprun.xml 文件拷贝到桌面留作后面使用,如图 2。 图2 二、P4vasp 部分: 1. 用 p4vasp 软件打开刚刚拷贝出来的 vasprun.xml 文件,选择菜单栏中 Electronic 选项中 DOS+bands,就可以显示出大致的态密度图,如图 3,4。 图3 图4 2. 此时,选择图形上面 Show 选项中的 Bands,就可以显示出能带图,如图 5。选择菜单 栏中 Electronic 选项中 local DOS+bands control,会出现这样一个对话框,如图 6。 图5 图6 3. 图 7 中给出了一些参数设置。在对话框中最上面输入要画的元素名称,例如 Cr,因为 这里没有开自旋,所以选择了 both,根据自己需要自己调整,在轨道选择中选择自己要 画的轨道,例如要画 Cr 的 dxy 轨道,则只需要选择它就可以,其他的取消。接下来可 以选择划线的颜色(便于观察和其他轨道的区别),在颜色选项下面是设置大小,手动 输入数值,例如 2(可默认),最后单击 add new line 进行画线。若出现画错线的情况, 也可删除刚刚画的线。图 8 给出了 Cr 的 dxy 轨道的参考图。同理可以画出其它轨道的 图线或则不同原子的轨道(在最上面更改元素的名称即可)。切记画各个轨道的顺序, 后面会用到 图7 图8 4. 画好自己需要画的各个轨道后。这里我们共画了了 Cr 原子的 5 个不同 d 轨道和 S 原子 的 p 轨道,共计 6 个。我们选择 Graph 中的 Export,会出现一个小的对话框,输入文件 的保存地址和类型,例如 E:\CrS2.dat,即把文件输出到 E 盘中,文件名称保存为 CrS2.dat 格式,如图 9。 图9 三、Vim 部分: 1. 找到刚刚保存的 CrS2.dat 文件,用 Vim 编译器打开,如图 10,切换到英文输入法, 输入:set nu ,回车即可显示出我们数据有多少行。前面说我们共画了 6 个轨道,我们 的数据共可以分为 7 个大块,每一块占据的行数是相等的。自身能带占据第一个大块(有 两列数据,接下来处理数据时要删除),后面一个轨道占一个大块(有三列数据,前两 列的内容和第一个大块的数据是一样的,最有用的是第三列数据),占据的顺序就是刚 刚在 p4vasp 中画轨道的顺序,按 shift+g(或 G)跳至底行可发现在本案例中共有 24303 行,如图 11。故每一大块占据 3472 行(具体数值视实际情况而定,在后面会经常用到 这个数值和它的倍数)。输入数字 3472 后按 shift+g 可以跳至指定的第 3472 行(第一 大块结尾),验证了第一大块的行数(包含空行)为 3472 行,如图 12。连按两次 g 可 跳至顶行,输入指令 3472dd 即可删除第一大块的内容。 图 10 图 11 图 12 2. 此时共有 5 个大块的数据,我们把后面 5 个大块中的第三列里的内容复制到第一大块的 内容后面。方法如下:跳至第 3472 行,把光标一直 3473 行(第二大块的开始处)最后 面 0 处,按 Crtl+q 进入可视块,如图 13。在可视块下输入数值 6943(3472*2-1),按 shift+g 跳到 6943 行(第二大块的尾行),按下 y 键复制选中的内容,如图 14。双按 g 跳至顶行,在第一大块的后面按 p 键粘贴刚刚复制的内容,如图 15。删去刚刚复制内 容的前两列,方法与刚刚复制内容的方法一样,即把光标移至第二块内容中第二列第一 行的后面,按按 Crtl+q 进入可视块,在可视块下输入数值 3472,按 shift+g 跳到 3472 行,按 d 键即可删除所选内容。 图 13 图 14 图 15 3. 由于我们画了 Cr 原子的 5 个 d 轨道和 S 原子的 p 轨道共 6 个轨道,如图 16 所示,共 有 8 列数据,前两列对应着能带能带图,后面分别对应着开始在 p4vasp 中各个轨道的 顺序,输入 :wq 即可退出保存。 图 16 四、Original 部分 1. 把刚才用 Vim 编译器处理后的 CrS2.dat 文件用 original 打开,选中前两列画出能带图, 格式、刻度等自己调节,再选中前两列加上第 3 列(或 4、5、6…列,分别对应于各个 轨道)右击→plot→Symbol→Buble,如图 17。 图 17 2. 右击刚刚画好的能带图,选择 Copy Format→All,再粘贴到刚刚画的泡泡图上,就可以 使得泡泡图与先前的能带图具有一样的规格,如图 18。 图 18 3. 双击泡泡图,可以更改泡泡图的大小、颜色、形状等参数,如图 19。 图 19 4. 调整好泡泡图各项参数后,单击菜单栏中 Format→Size/Speed,把最下面的两个复选框 中的对号取消,如图 20。 图 20 5. 同理,画出其他轨道的泡泡图,处理方法与前面一样。画好所有轨道的泡泡图后,单击 菜单栏中的 Graph→Merge Graph Windows→Open Dialog ..,把前两行空格里的数字改 为“1”就可以将所有图合并到一起了,最终图形如图 21(仅供参考)。 图 21