ICOHPLIST.lobster输出文件中ICOHP的轨道角动量分解非常细,自旋极化体系还有两套轨道,有时还要处理多对原子间的相互作用,数据处理较为繁琐
例如,单原子FeN4体系中Fe与第一壳层4个N原子的ICOHP计算
Fe原子考虑了3dx2-y2, 3dxy, 3dxz. 3dyz. 3dz2, 4s总计6个轨道
4个N原子考虑了2s, 2px, 2py, 2pz总计4*4个轨道
由于自旋极化,自旋向上和自旋向下是分别的2套轨道
此时要是想比较Fe的5个d轨道中那个轨道和N原子总的成键作用最强,直接在Excel里挨个加和将会非常繁琐
在这种情况下,可以考虑使用Excel的“数据透视表”功能
具体操作如下:
文件 - 打开 - ICOHPLIST.lobster - 另存为为.xlsx格式
将Fe轨道所在列“atomMU”和ICOHP计算结果“ICOHP (at) eF”所在列并在一起,选中
插入 - 数据透视表 - 确定
在“数据透视表字段”区域勾选上“atomMU”和“ICOHP (at) eF”,并将“atomMU”拖入“行”区域,将“ICOHP (at) eF”拖入“值”区域
点击“值”区域“计数项:****”右侧的倒三角 - 值字段设置 - 求和 - 确定
即可完成对所有相同行标签项的ICOHP求和
从求和结果可以看到,Fe的5个d轨道中对成键作用贡献最大的是dxy
排序:dxy(-3.983eV) >> dz2(-0.622eV) > dxz(-0.536eV) > dyz(-0.370eV) > dx2-y2(-0.210eV)
这在与定性分子轨道分析相符的基础上,进一步作了定量化