VASP+LOBSTER绘制COHP与PDOS合并图：从原理到实践

1. 项目概述：从COHP与PDOS到合并图的绘制逻辑

在材料计算和电子结构分析领域，COHP（Crystal Orbital Hamilton Population，晶体轨道哈密顿布居）和PDOS（Projected Density of States，投影态密度）是两个揭示化学键本质的利器。很多刚接触这个领域的朋友，在分别画出COHP和PDOS图后，常常会问：如何把这两张图合并到一张图上，进行更直观的对比分析？这看似是一个简单的绘图技巧问题，背后实则是对电子结构分析逻辑的深度整合。

简单来说，COHP告诉你原子轨道之间成键、反键和非键相互作用的“强度”和“能量位置”，它直接量化了化学键的强弱和性质。而PDOS则告诉你特定原子或轨道的电子态在能量上的分布情况。将两者合并，相当于把“键的相互作用信息”和“轨道的电子分布信息”放在同一个能量坐标轴下进行对照。这样，你就能一眼看出：在某个能量区间（比如费米能级附近的成键区域），是哪些原子轨道在主导成键？这些成键相互作用的强度如何？这对于理解材料的稳定性、导电性、磁性等性质至关重要。

我自己在分析过渡金属化合物、二维材料以及异质结界面时，频繁使用这种合并图。它不仅能验证计算结果的合理性，更能从电子层面揭示许多实验现象背后的机理。本文将基于VASP计算输出的常见数据格式，手把手带你走通从数据处理到最终绘图的完整流程，并分享一些让图表更专业、分析更到位的实战心得。

2. 核心概念解析：COHP与PDOS为何需要“同台竞技”

在深入操作之前，我们必须厘清COHP和PDOS各自扮演的角色，以及合并它们能带来什么额外的信息维度。这是正确解读合并图的前提。

2.1 PDOS：轨道的“能量户口本”

PDOS可以理解为特定原子、特定轨道（如s, p, d, f）甚至特定磁量子数轨道（如dxy, dyz, dz2, dxz, dx2-y2）的电子态密度在能量轴上的分布。它回答的问题是：“在某个能量E处，这个轨道有多少电子态可以占据？” 它的数值通常是正值（有时为了区分自旋向上向下会取正负），单位是states/eV。

通过PDOS，我们可以：

• 识别特征峰 ：找到来自特定轨道的电子态集中的能量区域。
• 判断轨道杂化 ：比较不同轨道PDOS峰的相对位置和重叠程度，判断它们是否发生了杂化。
• 分析费米能级处的态密度 ：这对于理解金属、半导体或绝缘体的导电性非常关键。

但是，PDOS有一个局限性：它只告诉我们电子态“在哪里”和“有多少”，但没有直接告诉我们这些态之间是否相互作用，以及相互作用的性质。

2.2 COHP：化学键的“能量体检报告”

COHP则向前迈进了一大步。它通过将能带或态密度投影到成键/反键轨道上，量化了特定原子对之间轨道重叠的哈密顿矩阵元。简单类比，如果两个原子轨道重叠形成分子轨道，COHP会计算这个过程中能量降低（成键）或升高（反键）的“贡献量”。

• 负值COHP (-pCOHP) ：代表成键相互作用。数值越负（在图表中通常表现为向下的峰），表示在该能量区间内，这两个轨道之间的成键作用越强，对体系稳定的贡献越大。
• 正值COHP ：代表反键相互作用。数值越正（在图表中通常表现为向上的峰），表示在该能量区间内，轨道间的相互作用是 destabilizing（去稳定化）的。
• COHP积分 (ICOHP) ：将COHP曲线在某个能量范围内（通常是从能量最低点到费米能级）进行积分，得到一个总的数值。负的ICOHP绝对值越大，表示这对原子之间的键总体越强。

因此，COHP直接回答了：“原子A和原子B之间，在哪个能量区间成键或反键？强度如何？”

2.3 合并图的威力：建立“轨道-成键”的直接关联

将PDOS和-COHP（通常展示负值，即成键贡献）绘制在同一张图上，共享X轴（能量，单位eV，费米能级通常设为0点），就实现了信息的叠加。

合并图的核心解读逻辑 ：

1. 定位成键区域 ：首先在-COHP曲线上找到显著的负峰区域，这些就是强成键相互作用的能量区间。
2. 追溯贡献轨道 ：然后，将视线垂直向上移动到同一能量区间的PDOS部分。观察在这个区间内，哪些原子、哪些轨道的PDOS有显著的峰。
3. 建立对应关系 ：此时，PDOS峰所对应的轨道，很可能就是贡献于该成键相互作用的主要轨道。例如，如果你在-5 eV到-2 eV区间看到一个很强的-COHP负峰，同时发现在这个区间，原子A的d_z2轨道和原子B的p_z轨道PDOS峰都很显著且重叠很好，那么你就可以很有信心地得出结论：A-d_z2和B-p_z轨道在此能量区间形成了强σ键。

这种直观的对应关系，使得合并图成为阐释复杂体系成键机制不可替代的工具。它把抽象的能带或态密度信息，转化为了易于理解的化学图像。

3. 工具链与数据准备：计算后处理的关键步骤

绘制合并图不是一个单纯的绘图问题，而是一个计算后处理流程。通常，这个流程始于第一性原理计算软件（如VASP, Quantum ESPRESSO），经过专门的后处理程序（如LOBSTER, pymatgen），最后在绘图软件（如Origin, Python Matplotlib, Gnuplot）中完成。

3.1 计算与后处理软件选择

主流方案：VASP + LOBSTER 这是目前最成熟、最通用的方案。

• VASP ：负责完成结构优化和静态自洽计算，产生 WAVECAR 和 CHGCAR 等文件。计算时，需要设置足够高的精度（尤其是 ENCUT 和 KPOINTS ）和完整的投影子（ LORBIT = 11 或使用 PREC = Accurate 配合适当的设置），以确保波函数信息足够用于后续投影。
• LOBSTER ：这是一个专门用于从VASP输出中计算COHP、COOP（晶体轨道重叠布居）和投影到局域轨道态密度的强大工具。它通过将平面波基组投影到用户定义的