关于如何准确计算IrO2态密度中的分轨道分布（复现一篇PRL）

由于二氧化铱中的铱为6配位，属于八面体，因此在准确计算八面体的分轨道dxy，dzy，dxz，dx2-y2，dz2的态密度分布时，晶体坐标轴的方向应当与IrO2中的IrO6八面体的取向一致，materials studio导入IrO2的cif文件时，其xyz应该调整为下图中的xyz取向，否则计算出来的分轨道就会由于分解（如dxycosα）得到不准确的结果。在原始晶胞基础上通过创建√2 *√2*1晶胞，可以调整为正确的xyz方向：

操作如下：

在Redefine lattice 中，把A设置为0 0 1，B设置为1 -1 0，C设置为1 1 0. 这是由于新设置的xyz的每一个坐标轴的向量方向在原始的abc中的向量分别为（0 0 1）（1 -1 0）（1 1 0）.即创建一个√2 *√2*1晶胞

新生成的√2 *√2*1晶胞样子为这样:

再通过调整其xyz方向沿着晶胞的abc，即可，如下图在lattice Parameters的Advanced中 选择A along X，B in XY plane，然后点击 Re-orient to standard. 即得到了计算中的正确的xyz方向.

通过VASP计算得到的下图的两种不同xyz设置的总态密度可以看到明显的不同，上图为未经过xyz矫正为沿着八面体方向的态密度，下图为矫正过后的态密度

可以看出，经过调整后的态密度明显毛刺少了，这是由于未经过调整的分轨道，会在xyz上分解，即dxy×cosα，dxy×sinα…..从而造成大量毛刺出现。

利用p4vasp处理态密度时候应该注意，应当选取同八面体方向的IrO6中的Ir分析其态密度，这是由于，在√2 *√2*1晶胞中存在相互垂直的IrO6八面体，其中对于一方的dz2轨道而言，其dz2等同于垂直方向的IrO6的Ir的dx2-y2. 反之也成立，即一方的dx2-y2轨道等于另一方的dz2轨道。如果选取了全部的Ir原子分析其DOS会发生错误的，错误表现为两个反键轨道合并到了成键键轨道中了，造成两个反键轨道为空轨道的假象。注意到这一点我们就可以选取正确的Ir原子得到正确的DOS图啦

最终得到的DOS与文献中（10.1103/PhysRevLett.112.117601）的对比

差不多比较完美了，可能提高KPOINT或者ENCUT的数值应该可以得到更精确的DOS了。